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();
}