[d2n-analysis-talk] DL.bit3 Variable (Followup From 9/2/10 Meeting)
David Flay
flay at jlab.org
Sat Sep 4 13:43:37 EDT 2010
Hi Brad,
So I looked into my code, and I had (for hit 1):
DL.bit3[Ndata.DL.bit3 - 2]
However, if Ndata.DL.bit3 < 2, I would have a negative index. So, to fix
that, I added the cut:
Ndata.DL.bit3!=1 to the rest of the code for that plot. So I have:
//hit 0
c1->cd(2);
gPad->SetLogx(0);
gPad->SetLogy(1);
cut3 = cut2+"&&(Ndata.DL.bit3!=0)&&(DL.bit3[Ndata.DL.bit3-1]>0)";
rootvar = "DL.bit3[Ndata.DL.bit3-1]";
range = ">>g0(100,0,2000)";
plot = rootvar+range;
ch->Draw(plot,cut3);
g0->SetTitle("DL.bit3[2] [hit 0] {Good Electron Cuts}");
g0->GetXaxis()->SetTitle(rootvar+" [hit 0]");
g0->GetXaxis()->CenterTitle();
Double_t hit0 = g0->GetEntries();
c1->Update();
//hit 1
c1->cd(3);
gPad->SetLogx(0);
gPad->SetLogy(1);
cut3 =
cut2+"&&(Ndata.DL.bit3!=0)&&(Ndata.DL.bit3!=1)&&(DL.bit3[Ndata.DL.bit3-2]>0)";
rootvar = "DL.bit3[Ndata.DL.bit3-2]";
range = ">>g1(100,-100,1900)";
plot = rootvar+range;
ch->Draw(plot,cut3);
g1->SetTitle("DL.bit3[1] [hit 1] {Good Electron Cuts}");
g1->GetXaxis()->SetTitle(rootvar+" [hit 1]");
g1->GetXaxis()->CenterTitle();
Double_t hit1 = g1->GetEntries();
c1->Update();
utilizing this code, I get the behavior I would have expected (that is, no
large peak at DL.bit3 <=20). See attached plot.
-------------------------------------------------
David Flay
Physics Department
Temple University
Philadelphia, PA 19122
office: Barton Hall, BA319
phone: (215) 204-1331
e-mail: flay at jlab.org
flay at temple.edu
website: http://www.jlab.org/~flay
http://quarks.temple.edu
-------------------------------------------------
-------------- next part --------------
A non-text attachment was scrubbed...
Name: DL-bit3_study-indiv-hits_20676_9_4_10.png
Type: image/png
Size: 31129 bytes
Desc: not available
Url : https://mailman.jlab.org/pipermail/d2n-analysis-talk/attachments/20100904/e72a3717/attachment-0001.png
More information about the d2n-analysis-talk
mailing list