new histos for embedding (Salvatore)

diff --git a/PWGJE/EMCALJetTasks/AliHadCorrTask.cxx b/PWGJE/EMCALJetTasks/AliHadCorrTask.cxx
index 8202cd0e86368081c0319b40e410affcd34ef464..891bdfb7da2f680b4c921b30989a0cb7e65efe68 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliHadCorrTask.cxx
+++ b/PWGJE/EMCALJetTasks/AliHadCorrTask.cxx
@@ -58,8 +58,9 @@ AliHadCorrTask::AliHadCorrTask() :
 
       fHistEmbTrackMatchesOversub[i] = 0;
       fHistNonEmbTrackMatchesOversub[i] = 0;
+      fHistOversubMCClusters[i] = 0;
+      fHistOversubNonMCClusters[i] = 0;
       fHistOversub[i] = 0;
-      fHistOversubOverClusE[i] = 0;
 
       for(Int_t j=0; j<4; j++)
        fHistNCellsEnergy[i][j] = 0;
@@ -106,8 +107,9 @@ AliHadCorrTask::AliHadCorrTask(const char *name, Bool_t histo) :
 
       fHistEmbTrackMatchesOversub[i] = 0;
       fHistNonEmbTrackMatchesOversub[i] = 0;
+      fHistOversubMCClusters[i] = 0;
+      fHistOversubNonMCClusters[i] = 0;
       fHistOversub[i] = 0;
-      fHistOversubOverClusE[i] = 0;
 
       for(Int_t j=0; j<4; j++)
        fHistNCellsEnergy[i][j] = 0;
@@ -438,28 +440,34 @@ void AliHadCorrTask::UserCreateOutputObjects()
   if (fIsEmbedded) {
     for(Int_t icent=0; icent<fNcentBins; ++icent) {
       name = Form("fHistEmbTrackMatchesOversub_%d",icent);
-      fHistEmbTrackMatchesOversub[icent] = new TH2F(name, name, fNbins, fMinBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
+      fHistEmbTrackMatchesOversub[icent] = new TH2F(name, name, fNbins, fMinBinPt, fMaxBinPt, fNbins, 0, 1.2);
       fHistEmbTrackMatchesOversub[icent]->GetXaxis()->SetTitle("E_{clus}^{raw} (GeV)");
-      fHistEmbTrackMatchesOversub[icent]->GetYaxis()->SetTitle("E_{oversub} (GeV)");
+      fHistEmbTrackMatchesOversub[icent]->GetYaxis()->SetTitle("E_{oversub} / E_{clus}^{raw}");
       fOutput->Add(fHistEmbTrackMatchesOversub[icent]);
 
       name = Form("fHistNonEmbTrackMatchesOversub_%d",icent);
-      fHistNonEmbTrackMatchesOversub[icent] = new TH2F(name, name, fNbins, fMinBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
+      fHistNonEmbTrackMatchesOversub[icent] = new TH2F(name, name, fNbins, fMinBinPt, fMaxBinPt, fNbins, 0, 1.2);
       fHistNonEmbTrackMatchesOversub[icent]->GetXaxis()->SetTitle("E_{clus}^{raw} (GeV)");
-      fHistNonEmbTrackMatchesOversub[icent]->GetYaxis()->SetTitle("E_{oversub} (GeV)");
+      fHistNonEmbTrackMatchesOversub[icent]->GetYaxis()->SetTitle("E_{oversub} / E_{clus}^{raw}");
       fOutput->Add(fHistNonEmbTrackMatchesOversub[icent]);
 
+      name = Form("fHistOversubMCClusters_%d",icent);
+      fHistOversubMCClusters[icent] = new TH2F(name, name, fNbins, fMinBinPt, fMaxBinPt, fNbins, 0, 1.2);
+      fHistOversubMCClusters[icent]->GetXaxis()->SetTitle("E_{clus}^{raw} (GeV)");
+      fHistOversubMCClusters[icent]->GetYaxis()->SetTitle("E_{oversub} / E_{clus}^{raw}");
+      fOutput->Add(fHistOversubMCClusters[icent]);
+
+      name = Form("fHistOversubNonMCClusters_%d",icent);
+      fHistOversubNonMCClusters[icent] = new TH2F(name, name, fNbins, fMinBinPt, fMaxBinPt, fNbins, 0, 1.2);
+      fHistOversubNonMCClusters[icent]->GetXaxis()->SetTitle("E_{clus}^{raw} (GeV)");
+      fHistOversubNonMCClusters[icent]->GetYaxis()->SetTitle("E_{oversub} / E_{clus}^{raw}");
+      fOutput->Add(fHistOversubNonMCClusters[icent]);
+
       name = Form("fHistOversub_%d",icent);
-      fHistOversub[icent] = new TH2F(name, name, fNbins, fMinBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
+      fHistOversub[icent] = new TH2F(name, name, fNbins, fMinBinPt, fMaxBinPt, fNbins, 0, 1.2);
       fHistOversub[icent]->GetXaxis()->SetTitle("E_{clus}^{raw} (GeV)");
-      fHistOversub[icent]->GetYaxis()->SetTitle("E_{oversub} (GeV)");
+      fHistOversub[icent]->GetYaxis()->SetTitle("E_{oversub} / E_{clus}^{raw}");
       fOutput->Add(fHistOversub[icent]);
-
-      name = Form("fHistOversubOverClusE_%d",icent);
-      fHistOversubOverClusE[icent] = new TH2F(name, name, fNbins, fMinBinPt, fMaxBinPt, fNbins, 0, 1.2);
-      fHistOversubOverClusE[icent]->GetXaxis()->SetTitle("E_{clus}^{raw} (GeV)");
-      fHistOversubOverClusE[icent]->GetYaxis()->SetTitle("E_{oversub} / E_{clus}^{raw}");
-      fOutput->Add(fHistOversubOverClusE[icent]);
     }
   }
 
@@ -878,13 +886,17 @@ Double_t AliHadCorrTask::ApplyHadCorrAllTracks(AliEmcalParticle *emccluster, Dou
       }
 
       if (fIsEmbedded) {
-       fHistOversub[fCentBin]->Fill(energyclus, overSub);
-       fHistOversubOverClusE[fCentBin]->Fill(energyclus, overSub / energyclus);
+       fHistOversub[fCentBin]->Fill(energyclus, overSub / energyclus);
+
+       if (cluster->GetMCEnergyFraction() > 0.95)
+         fHistOversubMCClusters[fCentBin]->Fill(energyclus, overSub / energyclus);
+       else if (cluster->GetMCEnergyFraction() < 0.05)
+         fHistOversubNonMCClusters[fCentBin]->Fill(energyclus, overSub / energyclus);
 
        if (trkPMCfrac < 0.05)
-         fHistNonEmbTrackMatchesOversub[fCentBin]->Fill(energyclus, overSub);
+         fHistNonEmbTrackMatchesOversub[fCentBin]->Fill(energyclus, overSub / energyclus);
        else if (trkPMCfrac > 0.95)
-         fHistEmbTrackMatchesOversub[fCentBin]->Fill(energyclus, overSub);
+         fHistEmbTrackMatchesOversub[fCentBin]->Fill(energyclus, overSub / energyclus);
       }
     }
   }

diff --git a/PWGJE/EMCALJetTasks/AliHadCorrTask.h b/PWGJE/EMCALJetTasks/AliHadCorrTask.h
index 0e99f071c4559d9b67834f9a37f7c83f19dbf192..b2b1d5c52061f273b5acba2081d45450de0a6795 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliHadCorrTask.h
+++ b/PWGJE/EMCALJetTasks/AliHadCorrTask.h
@@ -72,10 +72,11 @@ class AliHadCorrTask : public AliAnalysisTaskEmcal {
   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 with embedded track matches (non-embedded matches < 5%)
-  TH2                   *fHistNonEmbTrackMatchesOversub[4]; //!Over-subtracted energy with non-embedded track matches (embedded matches < 5%)
-  TH2                   *fHistOversub[4];                   //!Over-subtracted energy
-  TH2                   *fHistOversubOverClusE[4];          //!Over-subtracted energy / cluster energy
+  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%)
+  TH2                   *fHistOversubNonMCClusters[4];      //!Over-subtracted energy / cluster energy (cluster MC energy fraction < 5%)
+  TH2                   *fHistOversub[4];                   //!Over-subtracted energy / cluster energy
 
  private:
   AliHadCorrTask(const AliHadCorrTask&);            // not implemented