]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGGA/CaloTrackCorrelations/macros/QA/DrawAnaCaloTrackQA.C
Merge branch 'feature-movesplit'
[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / macros / QA / DrawAnaCaloTrackQA.C
index 70a09b631d378b2eccd8430a5f03c538993c8d78..61d541d0596b277271ee52e5dd8c51ef80aae369 100644 (file)
@@ -22,7 +22,7 @@ Int_t color[]={kBlack,kRed,kOrange+1,kYellow+1,kGreen+2,kBlue,kCyan+1,kViolet,kM
 
 //_______________________________________________________________________
 void DrawAnaCaloTrackQA(TString listName = "Pi0IM_GammaTrackCorr_EMCAL_default",
-                        TString fileName = "AnalysisResultsQATrain.root",
+                        TString fileName = "AnalysisResults.root",
                         Bool_t export = kFALSE)
 {
 
@@ -96,10 +96,10 @@ void CaloQA()
   ccalo->cd(2);
   //gPad->SetLogy();
   
-  TH1F* hRaw  = (TH1F*) GetHisto("AnaPhoton_hCut_0_Open");
-  TH1F* hCorr = (TH1F*) GetHisto("AnaPhoton_hCut_4_NCells");
-  TH1F* hTM   = (TH1F*) GetHisto("AnaPhoton_hCut_7_Matching");
-  TH1F* hShSh = (TH1F*) GetHisto("AnaPhoton_hCut_9_PID");
+  TH1F* hRaw  = (TH1F*) GetHisto("AnaPhoton_hPt_Cut_0_Open");
+  TH1F* hCorr = (TH1F*) GetHisto("AnaPhoton_hPt_Cut_4_NCells");
+  TH1F* hTM   = (TH1F*) GetHisto("AnaPhoton_hPt_Cut_7_Matching");
+  TH1F* hShSh = (TH1F*) GetHisto("AnaPhoton_hPt_Cut_9_PID");
   
   hRaw->Sumw2();
   
@@ -578,9 +578,9 @@ void CorrelationQA()
   Float_t assocBins[] = {0.5,2.,5.,10.,20.};
   Int_t nAssocBins = 4;
   
-  TH1F * hLeading = (TH1F*) GetHisto("AnaPhotonHadronCorr_hPtLeading");
-  Int_t minClusterEBin = hLeading->FindBin(minClusterE);
-  Float_t nTrig = hLeading->Integral(minClusterE,100000);
+  TH1F * hTrigger = (TH1F*) GetHisto("AnaPhotonHadronCorr_hPtTrigger");
+  Int_t minClusterEBin = hTrigger->FindBin(minClusterE);
+  Float_t nTrig = hTrigger->Integral(minClusterE,100000);
   
   //Azimuthal correlation
   cCorrelation->cd(1);
@@ -694,9 +694,13 @@ void MCQA()
   if(!h2ClusterPho) return;
   
   
-  TH1F* hPrimPho = (TH1F*) GetHisto("QA_hGenMCAccE_Photon");
-  TH1F* hPrimPi0 = (TH1F*) GetHisto("QA_hGenMCAccE_Pi0");
-  TH1F* hPrimEta = (TH1F*) GetHisto("QA_hGenMCAccE_Eta");
+//  TH1F* hPrimPho = (TH1F*) GetHisto("QA_hGenMCAccE_Photon");
+//  TH1F* hPrimPi0 = (TH1F*) GetHisto("QA_hGenMCAccE_Pi0");
+//  TH1F* hPrimEta = (TH1F*) GetHisto("QA_hGenMCAccE_Eta");
+
+  TH1F* hPrimPho = (TH1F*) GetHisto("AnaPhoton_hPtPrim_MCPhoton");
+  TH1F* hPrimPi0 = (TH1F*) GetHisto("AnaPi0_hPrimPi0Pt");
+  TH1F* hPrimEta = (TH1F*) GetHisto("AnaPi0_hPrimEtaPt");
   
   TCanvas * cmc = new TCanvas(Form("MCHisto_%s",histoTag.Data()),"",1000,1000);
   cmc->Divide(2,2);
@@ -713,7 +717,8 @@ void MCQA()
   hClusterPho->SetMarkerColor(1);
   hClusterPho->SetMarkerStyle(20);
   hClusterPho->SetAxisRange(0.,50.,"X");
-  hClusterPho->SetXTitle("E_{rec,gen} (GeV)");
+  //hClusterPho->SetXTitle("E_{rec,gen} (GeV)");
+  hClusterPho->SetXTitle("E_{rec}, p_{T,gen} (GeV)");
   hClusterPho->Draw("");
 
   hClusterPi0->Sumw2();
@@ -838,10 +843,10 @@ void MCQA()
   TH2F* h2PrimPhoEta = (TH2F*) GetHisto("AnaPhoton_hYPrim_MCPhoton");
   TH2F* h2PrimPi0Eta = (TH2F*) GetHisto("AnaPi0_hPrimPi0Rapidity");
   TH2F* h2PrimEtaEta = (TH2F*) GetHisto("AnaPi0_hPrimEtaRapidity");
-  
-  hPrimPhoEta->Sumw2();
-  hPrimEtaEta->Sumw2();
-  hPrimPi0Eta->Sumw2();
+
+  h2PrimPhoEta->Sumw2();
+  h2PrimEtaEta->Sumw2();
+  h2PrimPi0Eta->Sumw2();
   
   Int_t binMin = hPrimPho->FindBin(3);