2 Double_t x[18] = { .1, .2, .3, .4, .5, .6, .7, .8, .9,
3 1., 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8 };
4 Double_t
errx[18] = { 0 };
66 TFile
f(
"ana_bhabha660a_dqa_mcPat_zy_old.root" );
67 for ( Int_t j = 5; j < 15; j++ )
70 "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_"
71 "pxy0_em<0.15+%d*0.1&&NGch>=1",
73 Bhabha->Draw(
"kal_pxy0_em",
cut );
74 pxy_ep1[j] = htemp->GetEntries();
76 "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_"
77 "pxy0_em<0.15+%d*0.1&&NGch>=2",
79 Bhabha->Draw(
"kal_pxy0_em",
cut );
80 pxy_ep2[j] = htemp->GetEntries();
82 "kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_"
83 "pxy0_ep<0.15+%d*0.1&&NGch>=1",
85 Bhabha->Draw(
"kal_pxy0_ep",
cut );
86 pxy_em1[j] = htemp->GetEntries();
88 "kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_"
89 "pxy0_ep<0.15+%d*0.1&&NGch>=2",
91 Bhabha->Draw(
"kal_pxy0_ep",
cut );
92 pxy_em2[j] = htemp->GetEntries();
95 for ( Int_t i = 5; i < 15; i++ )
104 cout <<
"end eff1" << endl;
108 TFile
f(
"ana_bhabha660a_dqa_dataPat_zy_old.root" );
109 for ( Int_t j = 5; j < 15; j++ )
112 "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_"
113 "pxy0_em<0.15+%d*0.1&&NGch>=1",
115 Bhabha->Draw(
"kal_pxy0_em",
cut );
116 pxy_ep3[j] = htemp->GetEntries();
118 "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_"
119 "pxy0_em<0.15+%d*0.1&&NGch>=2",
121 Bhabha->Draw(
"kal_pxy0_em",
cut );
122 pxy_ep4[j] = htemp->GetEntries();
125 "kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_"
126 "pxy0_ep<0.15+%d*0.1&&NGch>=1",
128 Bhabha->Draw(
"kal_pxy0_ep",
cut );
129 pxy_em3[j] = htemp->GetEntries();
131 "kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_"
132 "pxy0_ep<0.15+%d*0.1&&NGch>=2",
134 Bhabha->Draw(
"kal_pxy0_ep",
cut );
135 pxy_em4[j] = htemp->GetEntries();
138 for ( Int_t i = 5; i < 15; i++ )
146 cout <<
"end eff2" << endl;
151 TFile
f(
"ana_bhabha660a_dqa_mcPat_zy_new.root" );
152 for ( Int_t j = 5; j < 15; j++ )
154 cout <<
"j=" << j << endl;
156 "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_"
157 "pxy0_em<0.15+%d*0.1&&NGch>=1",
159 Bhabha->Draw(
"kal_pxy0_em",
cut );
160 pxy_ep5[j] = htemp->GetEntries();
162 "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_"
163 "pxy0_em<0.15+%d*0.1&&NGch>=2",
165 Bhabha->Draw(
"kal_pxy0_em",
cut );
166 pxy_ep6[j] = htemp->GetEntries();
169 "kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_"
170 "pxy0_ep<0.15+%d*0.1&&NGch>=1",
172 Bhabha->Draw(
"kal_pxy0_ep",
cut );
173 pxy_em5[j] = htemp->GetEntries();
175 "kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_"
176 "pxy0_ep<0.15+%d*0.1&&NGch>=2",
178 Bhabha->Draw(
"kal_pxy0_ep",
cut );
179 pxy_em6[j] = htemp->GetEntries();
182 for ( Int_t i = 5; i < 15; i++ )
191 cout <<
"end eff3" << endl;
197 TFile
f(
"ana_bhabha660a_dqa_dataPat_zy_new.root" );
198 for ( Int_t j = 5; j < 15; j++ )
201 "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_"
202 "pxy0_em<0.15+%d*0.1&&NGch>=1",
204 Bhabha->Draw(
"kal_pxy0_em",
cut );
205 pxy_ep7[j] = htemp->GetEntries();
207 "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_"
208 "pxy0_em<0.15+%d*0.1&&NGch>=2",
210 Bhabha->Draw(
"kal_pxy0_em",
cut );
211 pxy_ep8[j] = htemp->GetEntries();
214 "kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_"
215 "pxy0_ep<0.15+%d*0.1&&NGch>=1",
217 Bhabha->Draw(
"kal_pxy0_ep",
cut );
218 pxy_em7[j] = htemp->GetEntries();
220 "kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_"
221 "pxy0_ep<0.15+%d*0.1&&NGch>=2",
223 Bhabha->Draw(
"kal_pxy0_ep",
cut );
224 pxy_em8[j] = htemp->GetEntries();
227 for ( Int_t i = 5; i < 15; i++ )
236 cout <<
"end eff4" << endl;
243 TFile
f(
"./ana_bhabha660a_dqa_mc_Pat_newKal.root" );
244 for ( Int_t j = 5; j < 15; j++ )
247 "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_"
248 "pxy0_em<0.15+%d*0.1&&NGch>=1",
250 Bhabha->Draw(
"kal_pxy0_em",
cut );
251 pxy_ep9[j] = htemp->GetEntries();
253 "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_"
254 "pxy0_em<0.15+%d*0.1&&NGch>=2",
256 Bhabha->Draw(
"kal_pxy0_em",
cut );
260 "kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_"
261 "pxy0_ep<0.15+%d*0.1&&NGch>=1",
263 Bhabha->Draw(
"kal_pxy0_ep",
cut );
264 pxy_em9[j] = htemp->GetEntries();
266 "kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_"
267 "pxy0_ep<0.15+%d*0.1&&NGch>=2",
269 Bhabha->Draw(
"kal_pxy0_ep",
cut );
273 for ( Int_t i = 5; i < 15; i++ )
282 cout <<
"end eff5" << endl;
328 TCanvas* c1 =
new TCanvas(
"c1",
"bhabha_eff", 100, 10, 600, 400 );
329 TMultiGraph* mg =
new TMultiGraph();
333 mg->SetTitle(
"bhabha_eff" );
335 gr1->SetLineColor( 1 );
336 gr2->SetLineColor( 2 );
337 gr3->SetLineColor( 1 );
338 gr4->SetLineColor( 2 );
339 gr5->SetLineColor( 1 );
340 gr6->SetLineColor( 1 );
341 gr7->SetLineColor( 2 );
342 gr8->SetLineColor( 2 );
343 gr9->SetLineColor( 5 );
344 gr10->SetLineColor( 5 );
348 gr1->SetMarkerColor( 1 );
349 gr2->SetMarkerColor( 2 );
350 gr3->SetMarkerColor( 1 );
351 gr4->SetMarkerColor( 1 );
352 gr5->SetMarkerColor( 1 );
353 gr6->SetMarkerColor( 1 );
354 gr7->SetMarkerColor( 2 );
355 gr8->SetMarkerColor( 2 );
356 gr9->SetMarkerColor( 5 );
357 gr10->SetMarkerColor( 5 );
361 gr1->SetMarkerStyle( 8 );
362 gr2->SetMarkerStyle( 8 );
363 gr3->SetMarkerStyle( 20 );
364 gr4->SetMarkerStyle( 20 );
365 gr5->SetMarkerStyle( 8 );
366 gr6->SetMarkerStyle( 8 );
367 gr7->SetMarkerStyle( 20 );
368 gr8->SetMarkerStyle( 20 );
369 gr9->SetMarkerStyle( 28 );
370 gr10->SetMarkerStyle( 28 );
388 mg->GetYaxis()->SetRangeUser( 0.8, 1.02 );
389 mg->GetXaxis()->SetRangeUser( 0.5, 1.9 );
390 mg->GetXaxis()->SetTitle(
"Pt/GeV" );
391 mg->GetYaxis()->SetTitle(
"Eff" );
393 TLegend*
legend =
new TLegend( 0.7, 0.65, 0.85, 0.75 );
418 legend->AddEntry( gr3,
"e+ Data oldPat",
"p" );
427 legend->AddEntry( gr7,
"e+ Data newPat",
"p" );
sprintf(cut, "kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_" "pxy0_em<0.15+%d*0.1&&NGch>=2", j, j)
TFile f("ana_bhabha660a_dqa_mcPat_zy_old.root")
*********Class see also m_nmax DOUBLE PRECISION m_MasPhot DOUBLE PRECISION m_phsu DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_r2 DOUBLE PRECISION m_WtMass INTEGER m_nmax INTEGER m_Nevgen INTEGER m_IsFSR INTEGER m_MarTot *COMMON c_KarFin $ !Output file $ !Event serial number $ !alpha QED at Thomson limit $ !minimum energy at CMS for remooval $ !infrared cut