From ac31ee0d2db07015eca41781eddeea9aa00b2a78 Mon Sep 17 00:00:00 2001 From: dmuhlhei Date: Thu, 11 Feb 2016 10:36:14 +0100 Subject: [PATCH] adding histo to CaloPhotonCuts --- PWGGA/GammaConv/AliCaloPhotonCuts.cxx | 10 ++++++++++ PWGGA/GammaConv/AliCaloPhotonCuts.h | 3 ++- 2 files changed, 12 insertions(+), 1 deletion(-) diff --git a/PWGGA/GammaConv/AliCaloPhotonCuts.cxx b/PWGGA/GammaConv/AliCaloPhotonCuts.cxx index 6b9daba6197..2ed9bbf6773 100644 --- a/PWGGA/GammaConv/AliCaloPhotonCuts.cxx +++ b/PWGGA/GammaConv/AliCaloPhotonCuts.cxx @@ -193,6 +193,7 @@ AliCaloPhotonCuts::AliCaloPhotonCuts(Bool_t isJetJet, const char *name,const cha fHistClusterIncludedCellsAfterQA(NULL), fHistClusterEnergyFracCellsBeforeQA(NULL), fHistClusterEnergyFracCellsAfterQA(NULL), + fHistClusterIncludedCellsTimingAfterQA(NULL), fHistClusterRBeforeQA(NULL), fHistClusterRAfterQA(NULL), fHistClusterdEtadPhiBeforeQA(NULL), @@ -329,6 +330,7 @@ AliCaloPhotonCuts::AliCaloPhotonCuts(const AliCaloPhotonCuts &ref) : fHistClusterIncludedCellsAfterQA(NULL), fHistClusterEnergyFracCellsBeforeQA(NULL), fHistClusterEnergyFracCellsAfterQA(NULL), + fHistClusterIncludedCellsTimingAfterQA(NULL), fHistClusterRBeforeQA(NULL), fHistClusterRAfterQA(NULL), fHistClusterdEtadPhiBeforeQA(NULL), @@ -607,6 +609,9 @@ void AliCaloPhotonCuts::InitCutHistograms(TString name){ fHistClusterEnergyFracCellsAfterQA = new TH1F(Form("ClusterEnergyFracCells_afterClusterQA %s",GetCutNumber().Data()),"ClusterEnergyFracCells_afterClusterQA",nMaxCellsEMCAL,0,nMaxCellsEMCAL); fHistClusterEnergyFracCellsAfterQA->Sumw2(); fHistExtQA->Add(fHistClusterEnergyFracCellsAfterQA); + fHistClusterIncludedCellsTimingAfterQA = new TH1F(Form("ClusterIncludedCellsTiming_afterClusterQA %s",GetCutNumber().Data()),"ClusterIncludedCellsTiming_afterClusterQA; t (ns)",200,-500,500); + fHistClusterIncludedCellsTimingAfterQA->Sumw2(); + fHistExtQA->Add(fHistClusterIncludedCellsTimingAfterQA); } else if( fClusterType == 2 ){ //PHOS Int_t nMaxCellsPHOS = fNMaxPHOSModules*56*64; @@ -634,6 +639,9 @@ void AliCaloPhotonCuts::InitCutHistograms(TString name){ fHistClusterEnergyFracCellsAfterQA = new TH1F(Form("ClusterEnergyFracCells_afterClusterQA %s",GetCutNumber().Data()),"ClusterEnergyFracCells_afterClusterQA",nMaxCellsPHOS,0,nMaxCellsPHOS); fHistClusterEnergyFracCellsAfterQA->Sumw2(); fHistExtQA->Add(fHistClusterEnergyFracCellsAfterQA); + fHistClusterIncludedCellsTimingAfterQA = new TH1F(Form("ClusterIncludedCellsTiming_afterClusterQA %s",GetCutNumber().Data()),"ClusterIncludedCellsTiming_afterClusterQA; t (ns)",200,-500,500); + fHistClusterIncludedCellsTimingAfterQA->Sumw2(); + fHistExtQA->Add(fHistClusterIncludedCellsTimingAfterQA); } else{AliError(Form("fExtendedMatchAndQA (%i) not (yet) defined for cluster type (%i)",fExtendedMatchAndQA,fClusterType));} } @@ -1142,6 +1150,8 @@ Bool_t AliCaloPhotonCuts::ClusterQualityCuts(AliVCluster* cluster, AliVEvent *ev for(Int_t iCell=0;iCellFill(cluster->GetCellAbsId(iCell)); if(cluster->E()>0.) fHistClusterEnergyFracCellsAfterQA->Fill(cluster->GetCellAbsId(iCell),cells->GetCellAmplitude(cluster->GetCellAbsId(iCell))/cluster->E()); + if(isMC==0) fHistClusterIncludedCellsTimingAfterQA->Fill(cells->GetCellTime(cluster->GetCellAbsId(iCell))*1E9); + else fHistClusterIncludedCellsTimingAfterQA->Fill(cells->GetCellTime(cluster->GetCellAbsId(iCell))*1E8); } } diff --git a/PWGGA/GammaConv/AliCaloPhotonCuts.h b/PWGGA/GammaConv/AliCaloPhotonCuts.h index 293736efe36..0e9d623d124 100644 --- a/PWGGA/GammaConv/AliCaloPhotonCuts.h +++ b/PWGGA/GammaConv/AliCaloPhotonCuts.h @@ -348,6 +348,7 @@ class AliCaloPhotonCuts : public AliAnalysisCuts { TH1F* fHistClusterIncludedCellsAfterQA; // CellIDs in Cluster of accepted ones TH1F* fHistClusterEnergyFracCellsBeforeQA; // Energy fraction of CellIDs in Cluster TH1F* fHistClusterEnergyFracCellsAfterQA; // Energy fraction of CellIDs in Cluster of accepted ones + TH1F* fHistClusterIncludedCellsTimingAfterQA; // Timing of CellIDs in Cluster of accepted ones //Track matching histograms TH1F* fHistClusterRBeforeQA; // cluster position in R=SQRT(x^2+y^2) (before QA) @@ -375,7 +376,7 @@ class AliCaloPhotonCuts : public AliAnalysisCuts { private: - ClassDef(AliCaloPhotonCuts,20) + ClassDef(AliCaloPhotonCuts,21) }; #endif