use all histo

diff --git a/PWGJE/EMCALJetTasks/AliHadCorrTask.cxx b/PWGJE/EMCALJetTasks/AliHadCorrTask.cxx
index db809ff863a8ba4f02a9ddf66495b2f4c3f32753..667958b2a2221e274a7af7ba7a02b8a11fca8ce2 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliHadCorrTask.cxx
+++ b/PWGJE/EMCALJetTasks/AliHadCorrTask.cxx
@@ -38,6 +38,7 @@ AliHadCorrTask::AliHadCorrTask() :
   fDoExact(kFALSE),
   fEsdMode(kTRUE),
   fOutClusters(0),
+  fHistMatchEtaPhiAll(0),
   fHistNclusvsCent(0),
   fHistNclusMatchvsCent(0),
   fHistEbefore(0),
@@ -87,6 +88,7 @@ AliHadCorrTask::AliHadCorrTask(const char *name, Bool_t histo) :
   fDoExact(kFALSE),
   fEsdMode(kTRUE),
   fOutClusters(0),
+  fHistMatchEtaPhiAll(0),
   fHistNclusvsCent(0),
   fHistNclusMatchvsCent(0),
   fHistEbefore(0),
@@ -380,6 +382,9 @@ void AliHadCorrTask::UserCreateOutputObjects()
 
   const Int_t nCentChBins = fNcentBins * 2;
 
+  fHistMatchEtaPhiAll = new TH2F("fHistMatchEtaPhiAll", "fHistMatchEtaPhiAll", fNbins, -0.1, 0.1, fNbins, -0.1, 0.1);
+  fOutput->Add(fHistMatchEtaPhiAll);
+
   for(Int_t icent=0; icent<nCentChBins; ++icent) {
     for(Int_t ipt=0; ipt<9; ++ipt) {
       for(Int_t ieta=0; ieta<2; ++ieta) {
@@ -597,6 +602,7 @@ void AliHadCorrTask::DoMatchedTracksLoop(AliEmcalParticle *emccluster, Double_t
       Int_t etabin = 0;
       if(track->Eta() > 0) etabin=1;
       fHistMatchEtaPhi[centbinch][mombin][etabin]->Fill(etadiff, phidiff);
+      fHistMatchEtaPhiAll->Fill(etadiff, phidiff);
     }
     
     Double_t etaCut   = 0.0;
@@ -745,6 +751,7 @@ Double_t AliHadCorrTask::ApplyHadCorrOneTrack(AliEmcalParticle *emccluster, Doub
       etabin = 1;
            
     fHistMatchEtaPhi[centbinch][mombin][etabin]->Fill(dEtaMin, dPhiMin);
+    fHistMatchEtaPhiAll->Fill(dEtaMin, dPhiMin);
     
     if (mom > 0) {
       fHistMatchEvsP[fCentBin]->Fill(energyclus, energyclus / mom);
@@ -866,7 +873,6 @@ Double_t AliHadCorrTask::ApplyHadCorrAllTracks(AliEmcalParticle *emccluster, Dou
          if (track) {
            Int_t centbinchm = fCentBin;
            if (track->Charge()<0) centbinchm += fNcentBins;
-           
            fHistEsubPchRat[centbinchm]->Fill(totalTrkP, Esub / totalTrkP);
            fHistEsubPch[centbinchm]->Fill(totalTrkP, Esub);
          }

diff --git a/PWGJE/EMCALJetTasks/AliHadCorrTask.h b/PWGJE/EMCALJetTasks/AliHadCorrTask.h
index 4149108709990bae5ea920aa236ed5349d868eed..c48cac73b24da715141c3e71ce834990fe21d495 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliHadCorrTask.h
+++ b/PWGJE/EMCALJetTasks/AliHadCorrTask.h
@@ -56,6 +56,7 @@ class AliHadCorrTask : public AliAnalysisTaskEmcal {
 
   // QA plots
   TH2                   *fHistMatchEtaPhi[8][9][2];  //!deta vs. dphi of matched cluster-track pairs
+  TH2                   *fHistMatchEtaPhiAll;        //!deta vs. dphi of matched cluster-track pairs
   TH2                   *fHistMatchEvsP[4];          //!cluster energy vs. track momentum of matched pairs
   TH2                   *fHistNMatchEnergy[4];       //!n matches vs. cluster energy
   TH2                   *fHistNCellsEnergy[4][4];    //!n cells vs. cluster energy
@@ -70,8 +71,6 @@ class AliHadCorrTask : public AliAnalysisTaskEmcal {
   TH1                   *fHistEsubPch[8];            //!Esub vs. total momentum of matched tracks (only 1 match)
   TH2                   *fHistEsubPchRat[8];         //!Esub/momentum of matched tracks vs. total momentum of matched tracks (only 1 match)
   TH2                   *fHistEsubPchRatAll[8];      //!Esub/momentum of matched tracks vs. total momentum of matched tracks (all number of matches)
-
-  // Embedding QA plots
   TH2                   *fHistEmbTrackMatchesOversub[4];    //!Over-subtracted energy / cluster energy with embedded track matches (non-embedded matches < 5%)
   TH2                   *fHistNonEmbTrackMatchesOversub[4]; //!Over-subtracted energy / cluster energy with non-embedded track matches (embedded matches < 5%)
   TH2                   *fHistOversubMCClusters[4];         //!Over-subtracted energy / cluster energy (cluster MC energy fraction > 95%)
@@ -82,6 +81,6 @@ class AliHadCorrTask : public AliAnalysisTaskEmcal {
   AliHadCorrTask(const AliHadCorrTask&);            // not implemented
   AliHadCorrTask &operator=(const AliHadCorrTask&); // not implemented
 
-  ClassDef(AliHadCorrTask, 13) // Hadronic correction task
+  ClassDef(AliHadCorrTask, 14) // Hadronic correction task
 };
 #endif