1 {
2 Double_t x[10] = { .1, .2, .3, .4, .5, .6, .7, .8, .9, 1. };
3 Double_t
errx[10] = { 0 };
4 Double_t effpt[10];
5 Double_t errpt[10];
7 Char_t filename[100];
8 Char_t name[100];
9
10 for ( Int_t i = 0; i < 10; i++ )
11 {
13 strcpy( filename,
file );
14 strcat( filename, name );
15
17 runinfo->Draw( "runinfo_ngood", "runinfo_ngood>0.5" );
18 effpt[i] = htemp->GetEntries() /
eventNo;
19 }
20
21 for ( Int_t i = 0; i < 10; i++ )
22 { errpt[i] = sqrt( effpt[i] * ( 1. - effpt[i] ) /
eventNo ); }
23
24 TGraph* gr1 =
new TGraphErrors( 10, x, effpt,
errx, errpt );
25
26 TCanvas* c1 = new TCanvas( "c1", "pioneff", 100, 10, 600, 400 );
27 TMultiGraph* mg = new TMultiGraph();
28
29 c1->SetGrid();
30 mg->SetTitle( "pioneff" );
31 gr1->SetLineColor( 2 );
32
33 mg->Add( gr1 );
34 mg->Draw( "ALP" );
35 mg->GetYaxis()->SetRangeUser( 0, 1 );
36 mg->GetXaxis()->SetTitle( "Pt/GeV" );
37 mg->GetYaxis()->SetTitle( "Eff" );
38
39 TLegend*
legend =
new TLegend( 0.7, 0.65, 0.85, 0.75 );
40 legend->AddEntry( gr1,
"pion_Pt",
"l" );
42}
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")