Validation/SingleAlg/SingleAlg-00-01-00/share/eff.cxx

Go to the documentation of this file.

{
gStyle->SetFrameFillColor(0);
gStyle->SetCanvasColor(0);
gStyle->SetFillColor(0);
 
gStyle->SetPadTickX(1);
gStyle->SetPadTickY(1);
 
gStyle->SetTitleXOffset(1.);
gStyle->SetTitleYOffset(1.);
 
gStyle->SetTitleXSize(0.06);
gStyle->SetTitleYSize(0.06);
 
gStyle->SetLabelSize(0.055,"x");
gStyle->SetLabelSize(0.055,"y");
gStyle->SetMarkerSize(1.8);
gStyle->SetLineWidth(2.5);
// gStyle->SetPadBorderMode(0);
 
// gStyle->SetPadBorderMode(0);
 
 
gStyle->SetPadLeftMargin(0.15);
gStyle->SetPadTopMargin(0.1);
gStyle->SetPadBottomMargin(0.15);
 
//Double_t  x[10]={.05,.15,.25,.35,.45,.55,.65,.75,.85,.95};
//Double_t  errx[20]={0.025,0.025,0.025,0.025,0.025,0.025,0.025,0.025,0.05,0.05,0.05,0.05,0.05,0.05,0.025,0.025,0.025,0.025,0.025,0.025};
Double_t  errx[20]={0};
//Double_t errx[10]={0.05,0.05,0.05,0.05,0.05,0.05,0.05,0.05,0.05,0.05};
 Double_t  x[20]={0.025,0.075,0.125,0.175,0.225,0.275,0.325,0.375,0.425,0.475,0.525,0.575,0.625,0.675,0.725,0.775,0.825,0.875,0.925,0.975};
 
Double_t eventNo=5000.;
Char_t  cut[200];
 
Int_t  kal_single[20];  
Int_t  mc_single[20];
Double_t  eff_single[20];
Double_t   err_single[20];
 
 
Int_t  kal_single2[20];  
Int_t  mc_single2[20];
Double_t  eff_single2[20];
Double_t   err_single2[20];
 
Int_t  kal_single3[20];  
Int_t  mc_single3[20];
Double_t  eff_single3[20];
Double_t   err_single3[20];
 
 
//////////////////////////////  1  /////////////////////////////////////////////////////////
//TFile f("../test.root");
TFile f("all10.root");
//TFile f("/ihepbatch/besd22/jixb/yuanwl/single/MC/p651p.root");
for (Int_t j=1; j<18; j++){
  sprintf(cut,"(runinfo_ngood==0||runinfo_ngood==1)&&kal_n>=0.+%d*0.05&&kal_n<0.05+%d*0.05&&costheta_n>-0.7&&costheta_n<0.7",j,j);
  runinfo->Draw("kal_n",cut);
  mc_single[j]=htemp->GetEntries();
 
  sprintf(cut,"runinfo_ngood==1&&kal_n>=0.+%d*0.05&&kal_n<0.05+%d*0.05&&costheta_n>-0.7&&costheta_n<0.7",j,j);
  runinfo->Draw("kal_n",cut);
  kal_single[j]=htemp->GetEntries();
 
}
 
 
 
 
//////////////////////////////  2  /////////////////////////////////////////////////////////
TFile f("all09.root");
//TFile f("../single_proton_gold0.3mm.root");
for (Int_t j=1; j<18; j++){
sprintf(cut,"(runinfo_ngood==0||runinfo_ngood==1)&&kal_n>=0.+%d*0.05&&kal_n<0.05+%d*0.05&&costheta_n>-0.7&&costheta_n<0.7",j,j);
  runinfo->Draw("kal_n",cut);
  mc_single2[j]=htemp->GetEntries();
  sprintf(cut,"runinfo_ngood==1&&kal_n>=0.+%d*0.05&&kal_n<0.05+%d*0.05&&costheta_n>-0.7&&costheta_n<0.7",j,j);
  runinfo->Draw("kal_n",cut);
  kal_single2[j]=htemp->GetEntries();
}
 
 
 
 
//////////////////////////////  3  /////////////////////////////////////////////////////////
TFile f("all09.root");
for (Int_t j=1; j<18; j++){
 
sprintf(cut,"(runinfo_ngood==0||runinfo_ngood==1)&&kal_n>=0.+%d*0.05&&kal_n<0.05+%d*0.05&&costheta_n>-0.7&&costheta_n<0.7",j,j);
  runinfo->Draw("kal_n",cut);
  mc_single3[j]=htemp->GetEntries();
  sprintf(cut,"runinfo_ngood==1&&kal_n>=0.+%d*0.05&&kal_n<0.05+%d*0.05&&costheta_n>-0.7&&costheta_n<0.7",j,j);
  runinfo->Draw("kal_n",cut);
  kal_single3[j]=htemp->GetEntries();
}
 
 
 
for (Int_t i=1;i<18;i++){
  cout<<mc_single[i]<<endl;
  eff_single[i]=kal_single[i]*.1/(mc_single[i]*.1);
  err_single[i]=sqrt(eff_single[i]*(1.-eff_single[i])/mc_single[i]);
 
  
  eff_single2[i]=kal_single2[i]*.1/(mc_single2[i]*.1);
  err_single2[i]=sqrt(eff_single2[i]*(1.-eff_single2[i])/mc_single2[i]);
 
  
  eff_single3[i]=kal_single3[i]*.1/(mc_single3[i]*.1);
  err_single3[i]=sqrt(eff_single3[i]*(1.-eff_single3[i])/mc_single3[i]);
 
}
 cout<<"............"<<endl;
 
 
TGraph *gr1 = new TGraphErrors(18,x,eff_single,errx,err_single);
TGraph *gr2 = new TGraphErrors(18,x,eff_single2,errx,err_single2);
TGraph *gr3 = new TGraphErrors(18,x,eff_single3,errx,err_single3);
 
 
 
TCanvas *c1 = new TCanvas ("c1","p",100,10,600,400);
TMultiGraph *mg = new TMultiGraph();
 
    
//c1->SetGrid();
 
//mg->SetTitle("p");
 
 
gr1->SetLineColor(2);
gr2->SetLineColor(3);
gr3->SetLineColor(4);
 
 
 
gr1->SetMarkerColor(2);
gr2->SetMarkerColor(3);
gr3->SetMarkerColor(4);
 
 
gr1->SetMarkerStyle(20);
gr2->SetMarkerStyle(22);
gr3->SetMarkerStyle(23);
 
 
// gr1->SetMarkerSize(1);
// gr2->SetMarkerSize(1.2);
// gr3->SetMarkerSize(1);
// gr4->SetMarkerSize(1.1);
 
 
mg->Add(gr1);
mg->Add(gr2);
//mg->Add(gr3);
 
 
 
mg->Draw("AP");
mg->GetYaxis()->SetRangeUser(0.8,1.02);
mg->GetXaxis()->SetRangeUser(0.05,1.05);
mg->GetXaxis()->SetTitle("Transvers Momentum/GeV");
mg->GetYaxis()->SetTitle("Tracking Efficiency");
 
 
TPaveText *text = new TPaveText(0.98,0.62,0.78,0.67,"tr");
 
text->AddText("Boss 6.5.1");
text->SetTextColor(2);
text->SetBorderSize(0);
 
//text->Draw();
 
 
 
TLegend *legend=new TLegend(0.55,0.45,0.75,0.65);
 
legend->AddEntry(gr1,"1T tsf+pat","p");
legend->AddEntry(gr2,"0.9T tsf+pat","p");
//legend->AddEntry(gr3,"proton_Paraffin0.93","p");
 
 
legend->SetBorderSize(0);
legend->Draw();
//c1->Print("p_eff_651.eps");
 
 
 
 
}