]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGLF/totEt/macros/hadEt/Misidentified.C
Split: removed dirs now in AliPhysics
[u/mrichter/AliRoot.git] / PWGLF / totEt / macros / hadEt / Misidentified.C
diff --git a/PWGLF/totEt/macros/hadEt/Misidentified.C b/PWGLF/totEt/macros/hadEt/Misidentified.C
deleted file mode 100644 (file)
index 4d0b345..0000000
+++ /dev/null
@@ -1,80 +0,0 @@
-//Christine Nattrass, University of Tennessee at Knoxville
-//This macro is for investigating the number of particles misidentified by the PID algorithm used for the transverse energy analysis
-//Misidentified particles are plotted as a percentage of the *total* number of particles.
-//Uses the output of AliAnalysisTaskHadEt
-
-void Misidentified(char *prodname = "Enter Production Name", char *shortprodname = "LHC10d4", bool TPC = true, char *filename="Et.ESD.new.sim.LHC10d4.pp.merged.root"){
-  char *myname = "TPCITS";
-  if(TPC) myname = "TPC";
-  TFile *file =  new TFile(filename);
-  if(!file){
-    cerr<<"Error, no file found"<<endl;
-    return;
-  }
-  TH2F *all = out2->FindObject(Form("MisidentifiedPIDs%s",myname));
-  Float_t totaln = ((TH2F*)out2->FindObject(Form("dEdxAll%s",myname)))->GetEntries();
-  all->Scale(100.0/totaln);
-  gStyle->SetPalette(1);
-
-  all->GetXaxis()->SetRange(all->GetXaxis()->FindBin(0.5),all->GetXaxis()->FindBin(4.));
-  all->GetYaxis()->SetRange(all->GetYaxis()->FindBin(0.0));
-  all->SetMarkerSize(2);
-  all->GetXaxis()->SetLabelSize(0.0);
-  all->GetYaxis()->SetLabelSize(0.0);
-  all->GetXaxis()->SetTitleSize(.06);
-  all->GetYaxis()->SetTitleSize(.06);
-  all->GetYaxis()->SetTitle("PID real");
-  all->GetXaxis()->SetTitle("PID identified");
-
-
-  gStyle->SetOptTitle(0);
-  gStyle->SetOptStat(0);
-  gStyle->SetOptFit(0);
-  TCanvas *c = new TCanvas("c","c",600,400);
-  c->SetTopMargin(0.02);
-  c->SetRightMargin(0.124161);
-  c->SetBottomMargin(0.147849);
-  c->SetLeftMargin(0.129195);
-  c->SetBorderSize(0);
-  c->SetFillColor(0);
-  c->SetFillColor(0);
-  c->SetBorderMode(0);
-  c->SetFrameFillColor(0);
-  c->SetFrameBorderMode(0);
-  all->Draw("colz");
-  all->Draw("textsame");
-
-  TLatex *pi = new TLatex(0.937284,-0.85435,"#pi^{#pm}");
-  pi->Draw();
-  TLatex *p = new TLatex(1.937284,-0.85435,"p,#bar{p}");
-  p->Draw();
-  TLatex *K = new TLatex(2.937284,-0.85435,"K^{#pm}");
-  K->Draw();
-  TLatex *E = new TLatex(3.937284,-0.85435,"e^{#pm}");
-  E->Draw();
-  TLatex *pi2 = new TLatex(0.297751,0.937284,"#pi^{#pm}");
-  pi2->Draw();
-  TLatex *p2 = new TLatex(0.297751,1.937284,"p,#bar{p}");
-  p2->Draw();
-  TLatex *K2 = new TLatex(0.297751,2.937284,"K^{#pm}");
-  K2->Draw();
-  TLatex *E2 = new TLatex(0.297751,3.937284,"e^{#pm}");
-  E2->Draw();
-  TLatex *other = new TLatex(0.297751,0.937284-1.0,"#mu^{#pm}");
-  other->Draw();
-
-
-  TLatex *tex = new TLatex(-0.0932606,-1.19453,prodname);
- tex->SetTextSize(0.0537634);
- tex->Draw();
-
-
-  if(TPC){
-    c->SaveAs(Form("pics/%s/Misidentified.eps",shortprodname));
-    c->SaveAs(Form("pics/%s/Misidentified.png",shortprodname));
-  }
-  else{
-    c->SaveAs(Form("pics/%s/MisidentifiedITS.eps",shortprodname));
-    c->SaveAs(Form("pics/%s/MisidentifiedITS.png",shortprodname));
-  }
-}