Validation/SingleAlg/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");
}