DQA/DQASelBhabha/DQASelBhabha-00-01-00/share/eff.cxx

Go to the documentation of this file.

{ 
Double_t  x[18]={.1,.2,.3,.4,.5,.6,.7,.8,.9,1.,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8};
Double_t  errx[18]={0};
 
Int_t  pxy_ep1[18]; 
Int_t  pxy_ep2[18]; 
Int_t  pxy_ep3[18]; 
Int_t  pxy_ep4[18];
Int_t  pxy_ep5[18]; 
Int_t  pxy_ep6[18];
Int_t  pxy_ep7[18]; 
Int_t  pxy_ep8[18]; 
Int_t  pxy_ep9[18]; 
Int_t  pxy_ep10[18];
Int_t  pxy_ep11[18]; 
Int_t  pxy_ep12[18];
 
Int_t  pxy_em1[18]; 
Int_t  pxy_em2[18]; 
Int_t  pxy_em3[18]; 
Int_t  pxy_em4[18]; 
Int_t  pxy_em5[18]; 
Int_t  pxy_em6[18]; 
Int_t  pxy_em7[18]; 
Int_t  pxy_em8[18]; 
Int_t  pxy_em9[18]; 
Int_t  pxy_em10[18]; 
Int_t  pxy_em11[18]; 
Int_t  pxy_em12[18];
 
Double_t  eff_ep[18]; 
Double_t  eff_em[18]; 
Double_t  err_ep[18]; 
Double_t  err_em[18];
 
Double_t  eff_ep2[18]; 
Double_t  eff_em2[18]; 
Double_t  err_ep2[18]; 
Double_t  err_em2[18];
 
Double_t  eff_ep3[18]; 
Double_t  eff_em3[18]; 
Double_t  err_ep3[18]; 
Double_t  err_em3[18];
 
Double_t  eff_ep4[18]; 
Double_t  eff_em4[18]; 
Double_t  err_ep4[18]; 
Double_t  err_em4[18];
 
Double_t  eff_ep5[18]; 
Double_t  eff_em5[18]; 
Double_t  err_ep5[18]; 
Double_t  err_em5[18];
 
Double_t  eff_ep6[18]; 
Double_t  eff_em6[18]; 
Double_t  err_ep6[18]; 
Double_t  err_em6[18];
 
 
 
Char_t  cut[200];
 
 
//TFile f("./ana_bhabha660_dqa_mc.root");
//TFile f("ana_bhabha660a_dqa_dataPat.root");
TFile f("ana_bhabha660a_dqa_mcPat_zy_old.root");
for (Int_t j=5; j<15; j++){
 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>=1",j,j);
  Bhabha->Draw("kal_pxy0_em",cut);
  pxy_ep1[j]=htemp->GetEntries();
  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);
  Bhabha->Draw("kal_pxy0_em",cut);
  pxy_ep2[j]=htemp->GetEntries();
  sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=1",j,j);
  Bhabha->Draw("kal_pxy0_ep",cut);
  pxy_em1[j]=htemp->GetEntries();
  sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=2",j,j);
  Bhabha->Draw("kal_pxy0_ep",cut);
  pxy_em2[j]=htemp->GetEntries();
}
 
 
for (Int_t i=5;i<15;i++){
  eff_ep[i]=pxy_ep2[i]*.1/(pxy_ep1[i]*.1);
  err_ep[i]=sqrt(eff_ep[i]*(1.-eff_ep[i])/pxy_ep1[i]);
 
  eff_em[i]=pxy_em2[i]*.1/(pxy_em1[i]*.1);
  err_em[i]=sqrt(eff_em[i]*(1.-eff_em[i])/pxy_em1[i]);
}
 
cout<<"end eff1"<<endl;
//////////////////////////////////////////////////////////////////////////////////
//TFile f("./ana_bhabha660a_dqa_mc_newPat.root");
//TFile f("ana_bhabha660a_dqa_dataTsf.root");
TFile f("ana_bhabha660a_dqa_dataPat_zy_old.root");
for (Int_t j=5; j<15; j++){
  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>=1",j,j);
  Bhabha->Draw("kal_pxy0_em",cut);
  pxy_ep3[j]=htemp->GetEntries();
  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);
  Bhabha->Draw("kal_pxy0_em",cut);
  pxy_ep4[j]=htemp->GetEntries();
 
  sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=1",j,j);
  Bhabha->Draw("kal_pxy0_ep",cut);
  pxy_em3[j]=htemp->GetEntries();
  sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=2",j,j);
  Bhabha->Draw("kal_pxy0_ep",cut);
  pxy_em4[j]=htemp->GetEntries();
 
}
 
 
for (Int_t i=5;i<15;i++){
  eff_ep2[i]=pxy_ep4[i]*.1/(pxy_ep3[i]*.1);
  err_ep2[i]=sqrt(eff_ep2[i]*(1.-eff_ep2[i])/pxy_ep3[i]);
 
  eff_em2[i]=pxy_em4[i]*.1/(pxy_em3[i]*.1);
  err_em2[i]=sqrt(eff_em2[i]*(1.-eff_em2[i])/pxy_em3[i]);
}
cout<<"end eff2"<<endl;
///////////////////////////////////////////////////////////////////////////////////////
//TFile f("./ana_bhabha660a_dqa_mc_Pat.root");
//TFile f("ana_bhabha660a_dqa_mc_Pat.root");
//TFile f("ana_bhabha660a_dqa_mcTsf_zy.root");
TFile f("ana_bhabha660a_dqa_mcPat_zy_new.root");
for (Int_t j=5; j<15; j++){
    cout<<"j="<<j<<endl;
  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>=1",j,j);
  Bhabha->Draw("kal_pxy0_em",cut);
  pxy_ep5[j]=htemp->GetEntries();
  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);
  Bhabha->Draw("kal_pxy0_em",cut);
  pxy_ep6[j]=htemp->GetEntries();
 
  sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=1",j,j);
  Bhabha->Draw("kal_pxy0_ep",cut);
  pxy_em5[j]=htemp->GetEntries();
  sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=2",j,j);
  Bhabha->Draw("kal_pxy0_ep",cut);
  pxy_em6[j]=htemp->GetEntries();
 
}
 
for (Int_t i=5;i<15;i++){
  eff_ep3[i]=pxy_ep6[i]*.1/(pxy_ep5[i]*.1);
  err_ep3[i]=sqrt(eff_ep3[i]*(1.-eff_ep3[i])/pxy_ep5[i]);
 
  eff_em3[i]=pxy_em6[i]*.1/(pxy_em5[i]*.1);
  err_em3[i]=sqrt(eff_em3[i]*(1.-eff_em3[i])/pxy_em5[i]);
}
 
cout<<"end eff3"<<endl;
/////////////////////////////////////////////////////////////////////////////////////
 
 
//TFile f("./ana_bhabha660a_dqa_mc_newPat_newKal.root");
//TFile f("ana_bhabha660a_dqa_mcTsf.root");
//TFile f("ana_bhabha660a_dqa_dataTsf_zy.root");
TFile f("ana_bhabha660a_dqa_dataPat_zy_new.root");
for (Int_t j=5; j<15; j++){
  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>=1",j,j);
  Bhabha->Draw("kal_pxy0_em",cut);
  pxy_ep7[j]=htemp->GetEntries();
  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);
  Bhabha->Draw("kal_pxy0_em",cut);
  pxy_ep8[j]=htemp->GetEntries();
 
  sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=1",j,j);
  Bhabha->Draw("kal_pxy0_ep",cut);
  pxy_em7[j]=htemp->GetEntries();
  sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=2",j,j);
  Bhabha->Draw("kal_pxy0_ep",cut);
  pxy_em8[j]=htemp->GetEntries();
 
}
 
for (Int_t i=5;i<15;i++){
  eff_ep4[i]=pxy_ep8[i]*.1/(pxy_ep7[i]*.1);
  err_ep4[i]=sqrt(eff_ep4[i]*(1.-eff_ep4[i])/pxy_ep7[i]);
 
  eff_em4[i]=pxy_em8[i]*.1/(pxy_em7[i]*.1);
  err_em4[i]=sqrt(eff_em4[i]*(1.-eff_em4[i])/pxy_em7[i]);
}
 
cout<<"end eff4"<<endl;
// /////////////////////////////////////////////////////////////////////////////////////
 
cout<<pxy_em2[6]<<"   "<<pxy_em2[7]<<endl;
cout<<pxy_em8[6]<<"   "<<pxy_em8[7]<<endl;
 
 
//TFile f("./ana_data653_mc_eff_nocut.root");
TFile f("./ana_bhabha660a_dqa_mc_Pat_newKal.root");
for (Int_t j=5; j<15; j++){
  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>=1",j,j);
  Bhabha->Draw("kal_pxy0_em",cut);
  pxy_ep9[j]=htemp->GetEntries();
  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);
  Bhabha->Draw("kal_pxy0_em",cut);
  pxy_ep10[j]=htemp->GetEntries();
 
  sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=1",j,j);
  Bhabha->Draw("kal_pxy0_ep",cut);
  pxy_em9[j]=htemp->GetEntries();
  sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=2",j,j);
  Bhabha->Draw("kal_pxy0_ep",cut);
  pxy_em10[j]=htemp->GetEntries();
 
}
 
for (Int_t i=5;i<15;i++){
  eff_ep5[i]=pxy_ep10[i]*.1/(pxy_ep9[i]*.1);
  err_ep5[i]=sqrt(eff_ep5[i]*(1.-eff_ep5[i])/pxy_ep9[i]);
 
  eff_em5[i]=pxy_em10[i]*.1/(pxy_em9[i]*.1);
  err_em5[i]=sqrt(eff_em5[i]*(1.-eff_em5[i])/pxy_em9[i]);
}
 
cout<<"end eff5"<<endl;
// /////////////////////////////////////////////////////////////////////////////////////
 
 
 
// //TFile f("./ana_data653_mc_eff_nomix.root");
// TFile f("./ana_data653_mc_eff_mdctest.root");
// for (Int_t j=5; j<15; j++){
//   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>=1",j,j);
//   Bhabha->Draw("kal_pxy0_em",cut);
//   pxy_ep11[j]=htemp->GetEntries();
//   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);
//   Bhabha->Draw("kal_pxy0_em",cut);
//   pxy_ep12[j]=htemp->GetEntries();
 
//   sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=1",j,j);
//   Bhabha->Draw("kal_pxy0_ep",cut);
//   pxy_em11[j]=htemp->GetEntries();
//   sprintf(cut,"kal_costheta0_ep>-0.93&&kal_costheta0_ep<0.93&&kal_pxy0_ep>=0.05+%d*0.1&&kal_pxy0_ep<0.15+%d*0.1&&NGch>=2",j,j);
//   Bhabha->Draw("kal_pxy0_ep",cut);
//   pxy_em12[j]=htemp->GetEntries();
 
// }
 
// for (Int_t i=5;i<15;i++){
//   eff_ep6[i]=pxy_ep12[i]*.1/(pxy_ep11[i]*.1);
//   err_ep6[i]=sqrt(eff_ep6[i]*(1.-eff_ep6[i])/pxy_ep11[i]);
 
//   eff_em6[i]=pxy_em12[i]*.1/(pxy_em11[i]*.1);
//   err_em6[i]=sqrt(eff_em6[i]*(1.-eff_em6[i])/pxy_em11[i]);
// }
 
// cout<<"end eff6"<<endl;
/////////////////////////////////////////////////////////////////////////////////////
 
 
 
 
 
 
 
 
 
 
 
 
TGraph *gr1 = new TGraphErrors(15,x,eff_ep,errx,err_ep);
TGraph *gr2 = new TGraphErrors(15,x,eff_em,errx,err_em);
TGraph *gr3 = new TGraphErrors(15,x,eff_ep2,errx,err_ep2);
TGraph *gr4 = new TGraphErrors(15,x,eff_em2,errx,err_em2);
TGraph *gr5 = new TGraphErrors(15,x,eff_ep3,errx,err_ep3);
TGraph *gr6 = new TGraphErrors(15,x,eff_em3,errx,err_em3);
TGraph *gr7 = new TGraphErrors(15,x,eff_ep4,errx,err_ep4);
TGraph *gr8 = new TGraphErrors(15,x,eff_em4,errx,err_em4);
TGraph *gr9 = new TGraphErrors(15,x,eff_ep5,errx,err_ep5);
TGraph *gr10 = new TGraphErrors(15,x,eff_em5,errx,err_em5);
// TGraph *gr11 = new TGraphErrors(15,x,eff_ep6,errx,err_ep6);
// TGraph *gr12 = new TGraphErrors(15,x,eff_em6,errx,err_em6);
 
 
 
 
 
TCanvas *c1 = new TCanvas ("c1","bhabha_eff",100,10,600,400);
TMultiGraph *mg = new TMultiGraph();
 
 
//c1->SetGrid();
 
mg->SetTitle("bhabha_eff");
 
gr1->SetLineColor(1);
gr2->SetLineColor(2);
gr3->SetLineColor(1);
gr4->SetLineColor(2);
gr5->SetLineColor(1);
gr6->SetLineColor(1);
gr7->SetLineColor(2);
gr8->SetLineColor(2);
gr9->SetLineColor(5);
gr10->SetLineColor(5);
// gr11->SetLineColor(2);
// gr12->SetLineColor(2);
 
gr1->SetMarkerColor(1);
gr2->SetMarkerColor(2);
gr3->SetMarkerColor(1);
gr4->SetMarkerColor(1);
gr5->SetMarkerColor(1);
gr6->SetMarkerColor(1);
gr7->SetMarkerColor(2);
gr8->SetMarkerColor(2);
gr9->SetMarkerColor(5);
gr10->SetMarkerColor(5);
// gr11->SetMarkerColor(2);
// gr12->SetMarkerColor(2);
 
gr1->SetMarkerStyle(8);
gr2->SetMarkerStyle(8);
gr3->SetMarkerStyle(20);
gr4->SetMarkerStyle(20);
gr5->SetMarkerStyle(8);
gr6->SetMarkerStyle(8);
gr7->SetMarkerStyle(20);
gr8->SetMarkerStyle(20);
gr9->SetMarkerStyle(28);
gr10->SetMarkerStyle(28);
// gr11->SetMarkerStyle(28);
// gr12->SetMarkerStyle(28);
 
 
//mg->Add(gr1);
//mg->Add(gr2);
mg->Add(gr3);
//mg->Add(gr4);
//mg->Add(gr5);
//mg->Add(gr6);
mg->Add(gr7);
//mg->Add(gr8);
//mg->Add(gr9);
//mg->Add(gr10);
//mg->Add(gr11);
//mg->Add(gr12);
 
 
mg->Draw("AP");
mg->GetYaxis()->SetRangeUser(0.8,1.02);
mg->GetXaxis()->SetRangeUser(0.5,1.9);
mg->GetXaxis()->SetTitle("Pt/GeV");
mg->GetYaxis()->SetTitle("Eff");
 
 
TLegend *legend=new TLegend(0.7,0.65,0.85,0.75);
//legend->AddEntry(gr1,"e+ 660","p");
//legend->AddEntry(gr2,"e- 660","p");
//legend->AddEntry(gr3,"e+ 660a newPat","p");
//legend->AddEntry(gr4,"e- 660a newPat","p"); 
//legend->AddEntry(gr5,"e+ 660a Pat","p");
//legend->AddEntry(gr6,"e- 660a Pat","p");
//legend->AddEntry(gr7,"e+ 660a newPat","p");
//legend->AddEntry(gr8,"e- 660a newPat","p"); 
//legend->AddEntry(gr9,"e+ 660a Pat_newKal","p");
//legend->AddEntry(gr10,"e- 660a Pat_newKal","p"); 
//legend->AddEntry(gr11,"e+ 653_MC_NoMix","p");
//legend->AddEntry(gr12,"e- 653_MC_NoMix","p");
 
//legend->AddEntry(gr1,"e+ Data PAT+TSF","p");
//legend->AddEntry(gr2,"e- Data PAT+TSF","p");
//legend->AddEntry(gr3,"e+ Data TSF","p");
//legend->AddEntry(gr4,"e- Data TSF","p");
//legend->AddEntry(gr5,"e+ MC PAT+TSF","p");
//legend->AddEntry(gr6,"e- MC PAT+TSF","p");
//legend->AddEntry(gr7,"e+ MC TSF","p");
//legend->AddEntry(gr8,"e- MC TSF","p");
 
//legend->AddEntry(gr1,"e+ MC oldPat","p");
//legend->AddEntry(gr2,"e- MC oldPat","p");
legend->AddEntry(gr3,"e+ Data oldPat","p");
//legend->AddEntry(gr4,"e- Data oldPat","p"); 
//legend->AddEntry(gr5,"e+ MC Tsf","p");
//legend->AddEntry(gr6,"e- MC Tsf","p");
//legend->AddEntry(gr7,"e+ Data Tsf","p");
//legend->AddEntry(gr8,"e- Data Tsf","p"); 
 
//legend->AddEntry(gr5,"e+ MC newPat","p");
//legend->AddEntry(gr6,"e- MC newPat","p");
legend->AddEntry(gr7,"e+ Data newPat","p");
//legend->AddEntry(gr8,"e- Data newPat","p"); 
 
 
legend->Draw();
 
//c1->Print("all_comp.eps");
 
}