Re: [BLAST_ANAWARE] chain and charge (fwd)

From: Chris Crawford (chris2@lns.mit.edu)
Date: Thu May 27 2004 - 09:37:51 EDT


hi tavi,
  'flrpm' should be identical to 'flr' and 'lr', except for the number
of events, of course. you can verify that with cuts on 'at least one
wc track' and 'two tracks, one + and one -'.
  with 'init.C' et al, you can have multiple datasets.

   root flr.C 1 2 3 4 5 -d 6 7

where 1-5 are data, and 6-7 are empty target runs. then in your macro,

  TString h("h")
  for (int nd=0;nd<2;nsets++) chain[nd]->Draw("twl>>"+h+nd, "pwl<1");
  TH1F h=h1-h2;

or of the likes.

try 'root init.C --help' for more details
--chris
 

Octavian F Filoti wrote:

>Hi Chris,
>
>why should not ask you also?! :-)
>Please, read the following message I sent to Chi.
>Thanks
>
>---------- Forwarded message ----------
>Date: Thu, 27 May 2004 00:20:09 -0400 (EDT)
>From: Octavian F Filoti <ofiloti@cisunix.unh.edu>
>To: Chi Zhang <zhangchi@MIT.EDU>
>Subject: chain and charge
>
>
>Hi Chi,
>
>I have been working for a while with 'flrpm', since it was the only one
>that accepted "chain[0,1,..]->' and TMatrix charge[0].
>Mistakenly I have discovered I have got "not good" results with it so far!
>(-:
>I am going to use 'lr' for a while, even though I have to do the charge
>part "manually", but my question is, how can I use 'chain' for the empty
>runs (to subtract them)?! I guess it should be a way of running them
>separately and then 'include' that file and those histograms, bla, bla,
>bla! How you (concretely) do it?! :-)
>
>P.S. I will look more carefull into the difference between 'lr' and
>'flrpm' and then 'bip' it to others! :-)
>
>-----------------------
>Octavian F Filoti
>Nuclear Physics Group
>Univ. of New Hampshire
>9 Library Way
>Durham, NH 03824
>phone: (603)862-1691
>FAX: (603)862-2998
>email: ofiloti@unh.edu
>
>



This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:31 EST