including more mc-checks

diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
index 6377bc8867f68ecf8e4614c09b311c4a9d80ce24..fa2376f907d7b9103d157cb6e0d04a85906b3cc6 100644 (file)
--- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
+++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
@@ -71,7 +71,7 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton() :
   fClusEtMcPt(0),
   fClusMcDetaDphi(0),
   fNClusPerPho(0),
-  fMCDirPhotonPtEtaNoClus(0),
+  fMCDirPhotonPtEtaPhiNoClus(0),
   fHnOutput(0)
 {
   // Default constructor.
@@ -116,7 +116,7 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton(const char *name) :
   fClusEtMcPt(0),
   fClusMcDetaDphi(0),
   fNClusPerPho(0),
-  fMCDirPhotonPtEtaNoClus(0),
+  fMCDirPhotonPtEtaPhiNoClus(0),
   fHnOutput(0)
 {
   // Constructor
@@ -178,8 +178,8 @@ void AliAnalysisTaskEMCALIsoPhoton::UserCreateOutputObjects()
   fNClusPerPho = new TH2F("hNClusPerPho","Number of clusters per prompt photon;p_{T}^{MC};N_{clus}",500,0,100,11,-0.5,10.5);
   fOutputList->Add(fNClusPerPho);
 
-  fMCDirPhotonPtEtaNoClus = new TH2F("hMCDirPhotonPtEtaNoClus","#eta vs. #phi prompt photons with no reco clusters;#phi;#eta",154,-0.77,0.77,130,1.38,3.20);
-  fOutputList->Add(fMCDirPhotonPtEtaNoClus);
+  fMCDirPhotonPtEtaPhiNoClus = new TH3F("hMCDirPhotonPhiEtaNoClus","p_{T}, #eta and  #phi of prompt photons with no reco clusters;p_{T};#eta;#phi",1000,0,100,154,-0.77,0.77,130,1.38,3.20);
+  fOutputList->Add(fMCDirPhotonPtEtaPhiNoClus);
 
   Int_t nEt=1000, nM02=400, nCeIso=1000, nTrIso=1000,  nAllIso=1000,  nCeIsoNoUE=1000,  nAllIsoNoUE=1000, nTrClDphi=200, nTrClDeta=100, nClEta=140, nClPhi=128, nTime=60, nMult=100, nPhoMcPt=101;
   Int_t bins[] = {nEt, nM02, nCeIso, nTrIso, nAllIso, nCeIsoNoUE, nAllIsoNoUE, nTrClDphi, nTrClDeta,nClEta,nClPhi,nTime,nMult,nPhoMcPt};
@@ -573,7 +573,7 @@ void AliAnalysisTaskEMCALIsoPhoton ::FillMcHists()
     if((imom==6 || imom==7) && pdgMom==22) {
       fMCDirPhotonPtEtaPhi->Fill(mcp->Pt(),mcp->Eta(),mcp->Phi());
       if(fNClusForDirPho==0)
-       fMCDirPhotonPtEtaNoClus->Fill(mcp->Eta(),mcp->Phi());
+       fMCDirPhotonPtEtaPhiNoClus->Fill(mcp->Pt(),mcp->Eta(),mcp->Phi());
       if(fDebug){
        printf("Found \"photonic\" parton at position %d, with pt=%1.1f, eta=%1.1f and phi=%1.1f, and status=%d,\n",imom,mcmom->Pt(), mcmom->Eta(), mcmom->Phi(), mcmom->GetStatusCode());
        printf("with a final photon at position %d, with pt=%1.1f, eta=%1.1f and phi=%1.1f, and status=%d\n",iTrack,mcp->Pt(), mcp->Eta(), mcp->Phi(),mcp->GetStatusCode());

diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h
index 07ac709e26c34e25217b0e0bc85acce245c9998b..63ee8aff8cd82e674f5c1fd6df89a511f1f0a383 100644 (file)
--- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h
+++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h
@@ -83,19 +83,19 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE {
 
   TList       *fOutputList; //! Output list
   //histograms for events with 1+ track pt>1
-  TH1F        *fEvtSel;                  //!evt selection counter: 0=all trg, 1=pv cut 
-  TH1F        *fNClusEt10;               //!number of clusters w/ Et>10 in the event
-  TH1F        *fPVtxZ;                   //!primary vertex Z before cut
-  TH1F        *fTrMultDist;              //!track multiplicity distribution
-  TH3F        *fMCDirPhotonPtEtaPhi;     //!direct produced photon pt
-  TH1F        *fDecayPhotonPtMC;         //!decay photon pt
-  TH2F        *fCellAbsIdVsAmpl;         //!cell abs id vs cell amplitude (energy)
-  TH2F        *fNClusHighClusE;          //!total number of clusters vs. highest clus energy in the event
-  TH2F        *fClusEtMcPt;              //!cluster et x mc-pt
-  TH2F        *fClusMcDetaDphi;          //!delta-eta x delta-phi(reco-mc)
-  TH2F        *fNClusPerPho;             //!delta-eta x delta-phi(reco-mc)
-  TH2F        *fMCDirPhotonPtEtaNoClus;  //!eta x phi for prompt photons that didn't produce clusters
-  THnSparse   *fHnOutput;                //!Output matrix with 7 dimensions
+  TH1F        *fEvtSel;                    //!evt selection counter: 0=all trg, 1=pv cut 
+  TH1F        *fNClusEt10;                 //!number of clusters w/ Et>10 in the event
+  TH1F        *fPVtxZ;                     //!primary vertex Z before cut
+  TH1F        *fTrMultDist;                //!track multiplicity distribution
+  TH3F        *fMCDirPhotonPtEtaPhi;       //!direct produced photon pt
+  TH1F        *fDecayPhotonPtMC;           //!decay photon pt
+  TH2F        *fCellAbsIdVsAmpl;           //!cell abs id vs cell amplitude (energy)
+  TH2F        *fNClusHighClusE;            //!total number of clusters vs. highest clus energy in the event
+  TH2F        *fClusEtMcPt;                //!cluster et x mc-pt
+  TH2F        *fClusMcDetaDphi;            //!delta-eta x delta-phi(reco-mc)
+  TH2F        *fNClusPerPho;               //!delta-eta x delta-phi(reco-mc)
+  TH3F        *fMCDirPhotonPtEtaPhiNoClus; //!pt x eta x phi for prompt photons that didn't produce clusters
+  THnSparse   *fHnOutput;                  //!Output matrix with 7 dimensions
 
   AliAnalysisTaskEMCALIsoPhoton(const AliAnalysisTaskEMCALIsoPhoton&); // not implemented
   AliAnalysisTaskEMCALIsoPhoton& operator=(const AliAnalysisTaskEMCALIsoPhoton&); // not implemented