Re: [BLAST_ANAWARE] commands for adding ntuples (fwd) (help!)

From: Jason Seely (seely@MIT.EDU)
Date: Sun Nov 03 2002 - 12:15:07 EST


hey guys,

i've done the initial zooming to find good proton peaks on which to cut. but i can't get
the individual cuts to work. if i say

%root -l show_protons.C 2450 -c prot

then i get tancredi's error below. if i try to send one of the elements of prot, say
plnn[12][4], then i get the
'no match' problem which i mentioned yesterday. so i'm hoping chris or chi has some idea
about what the error below means, and where i could look to fix it? i'll keep looking.

thanks,

jason

Tancredi Botto wrote:

> actually,
> it works ! Although my internet connection to a spud is surely slower than yours...!
>
> there is a
>
> Error: operator+= not defined for TCut FILE:/tmp/filey3Bj4z_cint LINE:1
> *** Interpreter error recovered ***
>
> which may be a typo somewhere, or just bad syntax
> chris, any ideas ?
>
> Should we just havev a chat room instead of a mailing list ??
> -- t
> ________________________________________________________________________________
> Tancredi Botto, phone: +1-617-253-9204 mobile: +1-978-490-4124
> research scientist MIT/Bates, 21 Manning Av Middleton MA, 01949
> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
>
> On Sat, 2 Nov 2002, Tancredi Botto wrote:
>
> >
> > Chirs,
> > please check jason's note about show_protons.C
> >
> > I believe the right sintax is
> > root -l show_protons.C -c prot 2450
> >
> > without brackets, as in -c prot[][]
> > -- t
> > ________________________________________________________________________________
> > Tancredi Botto, phone: +1-617-253-9204 mobile: +1-978-490-4124
> > research scientist MIT/Bates, 21 Manning Av Middleton MA, 01949
> > ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
> >
> > On Sat, 2 Nov 2002, Chris & Kim Crawford wrote:
> >
> > > hi all,
> > > here is an quick way to chain the 'ntp' ntuples together
> > > using the init.C script:
> > >
> > > bash> root -l ntp.C 1479 1480
> > > root[0] chain[0]->Draw("whatever")
> > > root[1] Draw16("ttl#t:ttl#b") //draw's 16 histograms automatically
> > > root[2] CutG() // after hitting enter, click on the bounds of your
> > > // graphical cut; double click on last point
> > > root[3] Draw("ttl2t:ttl3b") //it will now draw all points in black
> > > // and superimpose points passing the cut in red
> > > // (root still has problems w/ TCutG & TChain together
> > > // so this may not work until we upgrade root)
> > >
> > > for the 'eep' ntuple, use init.C instead of ntp.C.
> > >
> > > chaining is done internally using the 'TChain' class, which
> > > works as an ntuple across multiple files.
> > > TChain ntp("ntp")
> > > ntp->Add("new_datafile")
> > >
> > > --hope this helps, chris
> > >
> > >
> > >
> > > Tancredi Botto wrote:
> > > >
> > > > Sometimes you just need to string ntuples together.
> > > > Read on
> > > >
> > > > -- t
> > > > ________________________________________________________________________________
> > > > Tancredi Botto, phone: +1-617-253-9204 mobile: +1-978-490-4124
> > > > research scientist MIT/Bates, 21 Manning Av Middleton MA, 01949
> > > > ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
> > > >
> > > > ---------- Forwarded message ----------
> > > > Date: Thu, 31 Oct 2002 23:24:04 -0500
> > > > From: vitaliy ziskin <vziskin@yahoo.com>
> > > > To: Tancredi Botto <tancredi@lns.mit.edu>
> > > > Subject: commands for adding ntuples
> > > >
> > > > TFile run1("filename1")
> > > > TFile run1("filename2")
> > > > TFile run1("filename3")
> > > > TFile run1("filename4")
> > > >
> > > > TNtuple *ntp1
> > > > TNtuple *ntp2
> > > > TNtuple *ntp3
> > > > TNtuple *ntp4
> > > >
> > > > ntp1 = (TNtuple*)run1.Get("eep")
> > > > ntp2 = (TNtuple*)run2.Get("eep")
> > > > ntp3 = (TNtuple*)run3.Get("eep")
> > > > ntp4 = (TNtuple*)run4.Get("eep")
> > > >
> > > > ntp1->Draw("ncr>>htemp")
> > > > ntp2->Draw("ncr>>+htemp")
> > > > ntp13->Draw("ncr>>+htemp")
> > > > ntp4->Draw("ncr>>+htemp")
> > > >
> > > > htemp->Draw()
> > > >
> > > > // you don't need to declaire htemp
> > >
> >



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