I meant show_protons.C works with the -c prot cut work, that can (has to)
befine tuned to cut the proton peak in each TDC L-R or R-L plot.
My suggestion is: root -l show_protons.C -c padc_cc
you get very clean proton peaks. Zoom in around the peak and decide on a
cut. There is quite some ways to get confused since we have 256 L<>R
possible coincs. Let's start with some clean ones (elastic protons are not
everywhare). Do some tests as well to e time efficient
BTW: there are sample plots of variouos show_protons.C cuts on my desk.
After you define the cuts (in cuts.C) and wish to study the L/R
asymmetry you should
root -l show_protons.C -c padc -c prot
i.e. without a Cerenkov cut (for now)
of course when we'll calibrate against (tdc) offsets this will be just easier.
For this calibration all data since thursday evening (strobe in) is fine.
Peter, you want pitch in with a quick instruction on how toi use your
script ?
-- ________________________________________________________________________________ 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:
> > 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