//_______________________________________________________________________
void DrawAnaCaloTrackQA(TString listName = "Pi0IM_GammaTrackCorr_EMCAL_default",
- TString fileName = "AnalysisResultsQATrain.root",
+ TString fileName = "AnalysisResults.root",
Bool_t export = kFALSE)
{
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();
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);
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);
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();
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);