Single_trkeff_p.cxx

Go to the documentation of this file.

void Single_trkeff_p(char *filename="Single.root"){ 
Double_t  x[10]={.1,.2,.3,.4,.5,.6,.7,.8,.9,1.};
Double_t  errx[10]={0};
 
Int_t  kalp[10];  
Int_t  kalpp[10];  
Int_t  mcp[10];
 
Double_t  effp[10];
Double_t  effpp[10];
 
Double_t  errp[10];
Double_t   errpp[10];
 
Double_t eventNo=5000.;
Char_t  cut[100];
 
 
TFile f(filename);
 
for (Int_t j=0; j<10; j++){
  sprintf(cut,"kal_pxy>=0.05+%d*0.1&&kal_pxy<0.15+%d*0.1",j,j);
  pid_kal->Draw("kal_pxy",cut);
  kalp[j]=htemp->GetEntries(); 
  sprintf(cut,"kal_n>=0.05+%d*0.1&&kal_n<0.15+%d*0.1",j,j);
  pid_kal->Draw("kal_n",cut);
  kalpp[j]=htemp->GetEntries(); 
  sprintf(cut,"mcpart_pxy>=0.05+%d*0.1&&mcpart_pxy<0.15+%d*0.1",j,j );
  mcpart->Draw("mcpart_pxy",cut); 
  mcp[j]=htemp->GetEntries();
  }
 
 
 
for (Int_t i=0;i<10;i++){
  effp[i]=kalp[i]*.1/(mcp[i]*.1);
  errp[i]=sqrt(effp[i]*(1.-effp[i])/mcp[i]);
 
  effpp[i]=kalpp[i]*.1/(mcp[i]*.1);
  errpp[i]=sqrt(effpp[i]*(1.-effpp[i])/mcp[i]);
}
 
 
TGraph *gr1 = new TGraphErrors(10,x,effp,errx,errp);
TGraph *gr2 = new TGraphErrors(10,x,effpp,errx,errpp);
 
 
TCanvas *c1 = new TCanvas ("c1","Singleeff",100,10,600,400);
TMultiGraph *mg = new TMultiGraph();
 
    
c1->SetGrid();
 
mg->SetTitle("pioneff");
 
gr1->SetLineColor(1);
gr2->SetLineColor(2);
 
mg->Add(gr1);
mg->Add(gr2);
 
mg->Draw("ALP");
mg->GetYaxis()->SetRangeUser(0,1);
mg->GetXaxis()->SetTitle("Pt/GeV");
mg->GetYaxis()->SetTitle("Eff");
 
 
TLegend *legend=new TLegend(0.7,0.65,0.85,0.75);
legend->AddEntry(gr1,"Single_P","l");
legend->AddEntry(gr2,"Single_P'","l");
 
 
legend->Draw();
}