1 TCut
nhitC(
"(nhitx+nhity)>9&&nhitx>4&&nhity>4");
2 TCut
taveC(
"tave>0.0&&tave<550.0");
12 TCut
vetoCt(
"y1<130.0&&abs(x1)<100.0&&z1>50&&z1<770");
14 TCut
cosCt(
"cosnumi>0.85");
22 TCut
veto1(
"abs(x1)<110.0&&abs(y1)<140&&z1>30&&z1<790");
23 TCut
veto2(
"abs(x2)<110.0&&abs(y2)<140&&z2>30&&z2<790");
32 TCut
tCl(
"tave>216&&tave<230");
33 TCut
tCt(
"tave>216.1&&tave<228.3");
39 if (gROOT->IsBatch()) gDirectory->cd(
"Cana");
56 cana->Scan(
"evt:run:subrun:tave",
tightC);
60 cana->Scan(
"evt:run:subrun:tave",
tightC&&
tCt);
64 cana->Scan(
"evt:run:subrun:tave",
looseC&&
tCt);
72 TCanvas*
c1 =
new TCanvas(
"tave",
"tave");
73 TH1F*
h1 =
new TH1F(
"taveh",
"NuMI neutrino candidates;t-t_{0} (usec);Time Windows per 2 usec",350,-100,600);
75 c1->Print(
"tave.pdf");
76 c1->Print(
"tave.png");
79 TCanvas*
cs =
new TCanvas(
"taveS",
"taveS");
80 TH1F* hs =
new TH1F(
"taves",
"NuMI neutrino candidates - zoomed;t-t_{0} (usec);Time Windows per .5 usec",100,195,245);
82 cs->Print(
"tave-peak.png");
84 gStyle->SetOptStat(1111);
85 TCanvas*
c2 =
new TCanvas(
"tave-zoom",
"tave-zoom");
87 c2->Print(
"tave-zoom.pdf");
88 c2->Print(
"tave-zoom.png");
90 TH1F*
h2 =
new TH1F(
"tavefineh",
";t-t_{0} (usec)",5500,0,550);
91 cana->Draw(
"tave >> tavefineh",
tightC&&
tCl);
93 double tc = h2->GetMean();
94 double rms = h2->GetRMS();
95 double n = h2->Integral();
97 double sig = rms/
sqrt(n-1);
104 double ts = tc - (41./9.) - 3*sig;
105 double te = tc + (49./9.) + 3*sig;
107 std::cout <<
"Time window between " << ts <<
" and " << te
108 <<
" at 3 sigma (" << te-ts <<
" wide)" <<
std::endl;
116 TCanvas*
c2 =
new TCanvas(
"phib",
"phib");
117 TH1F*
h2 =
new TH1F(
"phib",
";#phi_{B}",40,-200,200);
119 c2->Print(
"phib.pdf");
120 c2->Print(
"phib.png");
122 TCanvas*
c1 =
new TCanvas(
"costhetab",
"cosnumi");
123 TH1F*
h1 =
new TH1F(
"costhetab",
";cos#theta_{B}",50,0.0,1.0);
125 c1->Print(
"costhetab.pdf");
126 c1->Print(
"costhetab.png");
133 gStyle->SetOptStat(0);
135 TCanvas*
c1 =
new TCanvas(
"x1pos",
"x1pos");
136 TH1F*
h1 =
new TH1F(
"x1pos",
"x1pos",40,-200,200);
138 c1->Print(
"x1pos.pdf");
139 c1->Print(
"x1pos.png");
141 TCanvas*
c2 =
new TCanvas(
"y1pos",
"y1pos");
142 TH1F*
h2 =
new TH1F(
"y1pos",
"y1pos",50,-250,250);
144 c2->Print(
"y1pos.pdf");
145 c2->Print(
"y1pos.png");
147 TCanvas*
c3 =
new TCanvas(
"z1pos",
"z1pos");
148 TH1F*
h3 =
new TH1F(
"z1pos",
"z1pos",90,-50,850);
150 c3->Print(
"z1pos.pdf");
151 c3->Print(
"z1pos.png");
155 TCanvas*
c4 =
new TCanvas(
"xz",
"xz");
156 TH2F* h4 =
new TH2F(
"x1z1",
";z1 (cm);x1 (cm)",
159 cana->Draw(
"x1:z1>>x1z1",
nhitC&&
taveC,
"colz");
160 TBox* b4 =
new TBox(50,-100,770,100);
167 TCanvas* c5 =
new TCanvas(
"yz",
"yz");
168 TH2F*
h5 =
new TH2F(
"y1z1",
";z1 (cm);y1 (cm)",
171 cana->Draw(
"y1:z1>>y1z1",
nhitC&&
taveC,
"colz");
172 TBox* b5 =
new TBox(50,-300,770,130);
179 TCanvas* c6 =
new TCanvas(
"yx",
"yx");
180 TH2F* h6 =
new TH2F(
"y1x1",
";x1 (cm);y1 (cm)",
183 cana->Draw(
"y1:x1>>y1x1",
nhitC&&
taveC,
"colz");
184 TBox* b6 =
new TBox(-100,-300,100,130);
197 gStyle->SetOptStat(0);
198 TCanvas*
c1 =
new TCanvas(
"angles1",
"angles1",600,700);
199 TH1F*
h1 =
new TH1F(
"zenith",
";cos(#theta_{Z});Ntracks",100,0.0,1.0);
200 h1->SetMarkerStyle(20);
202 cana->Draw(
"coscosmic >> zenith",
nhitC&&
taveC,
"e");
204 double sum1 = h1->Integral(1,80);
205 double sum2 =
pow(0.8,3)/3.0;
207 sprintf(buff,
"(%f/%f)*x**2",sum1,100*sum2);
208 TF1*
f2 =
new TF1(
"f2",buff,0,1);
211 c1->Print(
"angles1.pdf");
212 c1->Print(
"angles1.png");
214 TCanvas*
c2 =
new TCanvas(
"angles2",
"angles2",600,700);
215 TH1F*
h2 =
new TH1F(
"phi",
";#phi (deg.);Ntracks",72,-180,180);
216 h2->SetMarkerStyle(20);
217 cana->Draw(
"phicosmic >> phi",
nhitC&&
taveC,
"e");
219 c2->Print(
"phi.pdf");
220 c2->Print(
"phi.png");
227 gStyle->SetOptStat(1111);
229 TCanvas*
c1 =
new TCanvas(
"nhit1",
"nhit1",400,800);
230 TH1F*
h1 =
new TH1F(
"Nhits",
"Nhits;Nhits;Ntracks",100,0.0,100.0);
232 cana->Draw(
"nhit >> Nhits");
234 c1->Print(
"nhit1.pdf");
235 c1->Print(
"nhit1.png");
237 TCanvas*
c2 =
new TCanvas(
"nhit2",
"nhit2",400,800);
241 TH1F*
h2 =
new TH1F(
"NhitsX",
"X view;Nhits;Ntracks",50,0.0,50.0);
243 cana->Draw(
"nhitx >> NhitsX");
246 TH1F*
h3 =
new TH1F(
"NhitsY",
"Y view;Nhits;Ntracks",50,0.0,50.0);
248 cana->Draw(
"nhity >> NhitsY");
250 c2->Print(
"nhit2.pdf");
251 c2->Print(
"nhit2.png");
258 TCanvas*
c1 =
new TCanvas(
"run-summary",
"run-summary",800,800);
261 fNhitTotalVsEvent->Draw();
262 TProfile*
p = fNhitTotalVsEvent->ProfileX();
263 p->SetMarkerStyle(20);
264 p->SetMarkerColor(2);
269 fNtimeWinVsEvent->Draw();
270 p = fNtimeWinVsEvent->ProfileX();
271 p->SetMarkerStyle(20);
272 p->SetMarkerColor(2);
TCut veto1("abs(x1)<110.0&&abs(y1)<140&&z1>30&&z1<790")
TCut vetoCt("y1<130.0&&abs(x1)<100.0&&z1>50&&z1<770")
TCut tCl("tave>216&&tave<230")
TCut vetoCl(veto1||veto2)
TCut veto2("abs(x2)<110.0&&abs(y2)<140&&z2>30&&z2<790")
TCut nhitC("(nhitx+nhity)>9&&nhitx>4&&nhity>4")
TCut tightC_noTime(nhitC &&vetoCt &&phiCt &&cosCt)
TCut looseC(nhitC &&taveC &&vetoCl &&cosCl)
TCut cosCt("cosnumi>0.85")
TCut tCt("tave>216.1&&tave<228.3")
TCut taveC("tave>0.0&&tave<550.0")
TCut tightC(nhitC &&taveC &&vetoCt &&phiCt &&cosCt)