remove setting of distance to bad channel in Omega analsys

diff --git a/PWG4/PartCorrDep/AliAnaOmegaToPi0Gamma.cxx b/PWG4/PartCorrDep/AliAnaOmegaToPi0Gamma.cxx
index 7328991812ce83833bc078bf228922b9c8e685b2..f5cbef6cc352347ac4d2433368f556fe7475796d 100644 (file)
--- a/PWG4/PartCorrDep/AliAnaOmegaToPi0Gamma.cxx
+++ b/PWG4/PartCorrDep/AliAnaOmegaToPi0Gamma.cxx
@@ -43,8 +43,8 @@ ClassImp(AliAnaOmegaToPi0Gamma)
 //______________________________________________________________________________\r
 AliAnaOmegaToPi0Gamma::AliAnaOmegaToPi0Gamma() : AliAnaPartCorrBaseClass(),\r
 fInputAODGamma(0),fInputAODPi0(0), fInputAODGammaName(""),\r
-fEventsList(0),fNVtxZBin(1), fNCentBin(1), fNRpBin(1), fNBadChDistBin(1), fNpid(3),\r
-fNmaxMixEv(4), fVtxZCut(0), fCent(0), fRp(0), fBadChDist(0),\r
+fEventsList(0),fNVtxZBin(1), fNCentBin(1), fNRpBin(1), fNBadChDistBin(3), fNpid(3),\r
+fNmaxMixEv(4), fVtxZCut(0), fCent(0), fRp(0), \r
 fPi0Mass(0.13498),fPi0MassWindow(0.015),fPi0OverOmegaPtCut(0.8),\r
 fGammaOverOmegaPtCut(0.3),\r
 fhEtalon(0),\r
@@ -68,7 +68,7 @@ fEventsList(ex.fEventsList),
 fNVtxZBin(ex.fNVtxZBin), fNCentBin(ex.fNCentBin), fNRpBin(ex.fNRpBin),\r
 fNBadChDistBin(ex.fNBadChDistBin),fNpid(ex.fNpid),\r
 fNmaxMixEv(ex.fNmaxMixEv),\r
-fVtxZCut(ex.fVtxZCut), fCent(ex.fCent), fRp(ex.fRp), fBadChDist(ex.fBadChDist),\r
+fVtxZCut(ex.fVtxZCut), fCent(ex.fCent), fRp(ex.fRp), \r
 fPi0Mass(ex.fPi0Mass),\r
 fPi0MassWindow(ex.fPi0MassWindow),\r
 fPi0OverOmegaPtCut(ex.fPi0OverOmegaPtCut),\r
@@ -110,7 +110,6 @@ AliAnaOmegaToPi0Gamma & AliAnaOmegaToPi0Gamma::operator = (const AliAnaOmegaToPi
    fVtxZCut=ex.fVtxZCut;\r
    fCent=ex.fCent;\r
    fRp=ex.fRp;\r
-   fBadChDist=ex.fBadChDist;\r
 \r
    fPi0Mass=ex.fPi0Mass;\r
    fPi0MassWindow=ex.fPi0MassWindow;\r
@@ -165,7 +164,6 @@ AliAnaOmegaToPi0Gamma::~AliAnaOmegaToPi0Gamma() {
  delete [] fVtxZCut;\r
  delete [] fCent;\r
  delete [] fRp;\r
- delete [] fBadChDist;\r
 \r
 }\r
 \r
@@ -183,9 +181,6 @@ void AliAnaOmegaToPi0Gamma::InitParameters()
   fRp=new Double_t[fNRpBin];\r
   for(int i = 0;i<fNRpBin;i++)fRp[i]=0;\r
 \r
-  fBadChDist=new Int_t [fNBadChDistBin];\r
-  for(Int_t j=0;j<fNBadChDistBin;j++)   fBadChDist[j] = j+2;\r
-\r
 }\r
 \r
 \r
@@ -249,7 +244,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      Int_t index=idim*fNpid*fNBadChDistBin+ipid*fNBadChDistBin+idist;\r
 \r
                     sprintf(key,"RealToPi0Gamma_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s Real Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s Real Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fRealOmega0[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fRealOmega0[index]->SetName(key) ;\r
@@ -257,7 +252,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fRealOmega0[index]);\r
 \r
                     sprintf(key,"MixAToPi0Gamma_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s MixA Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s MixA Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fMixAOmega0[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fMixAOmega0[index]->SetName(key) ;\r
@@ -265,7 +260,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fMixAOmega0[index]);\r
 \r
                     sprintf(key,"MixBToPi0Gamma_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s MixB Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s MixB Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fMixBOmega0[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fMixBOmega0[index]->SetName(key) ;\r
@@ -273,7 +268,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fMixBOmega0[index]);\r
 \r
                     sprintf(key,"MixCToPi0Gamma_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s MixC Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s MixC Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fMixCOmega0[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fMixCOmega0[index]->SetName(key) ;\r
@@ -281,7 +276,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fMixCOmega0[index]);\r
 \r
                     sprintf(key,"RealToPi0Gamma1_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s Real Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s Real Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fRealOmega1[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fRealOmega1[index]->SetName(key) ;\r
@@ -289,7 +284,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fRealOmega1[index]);\r
 \r
                      sprintf(key,"MixAToPi0Gamma1_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s MixA Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s MixA Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fMixAOmega1[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fMixAOmega1[index]->SetName(key) ;\r
@@ -297,7 +292,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fMixAOmega1[index]);\r
 \r
                      sprintf(key,"MixBToPi0Gamma1_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s MixB Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s MixB Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fMixBOmega1[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fMixBOmega1[index]->SetName(key) ;\r
@@ -305,7 +300,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fMixBOmega1[index]);\r
 \r
                      sprintf(key,"MixCToPi0Gamma1_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s MixC Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s MixC Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fMixCOmega1[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fMixCOmega1[index]->SetName(key) ;\r
@@ -313,7 +308,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fMixCOmega1[index]);\r
 \r
                     sprintf(key,"RealToPi0Gamma2_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s Real Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s Real Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fRealOmega2[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fRealOmega2[index]->SetName(key) ;\r
@@ -321,7 +316,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fRealOmega2[index]);\r
 \r
                      sprintf(key,"MixAToPi0Gamma2_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s MixA Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s MixA Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fMixAOmega2[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fMixAOmega2[index]->SetName(key) ;\r
@@ -329,7 +324,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fMixAOmega2[index]);\r
 \r
                      sprintf(key,"MixBToPi0Gamma2_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s MixB Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s MixB Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fMixBOmega2[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fMixBOmega2[index]->SetName(key) ;\r
@@ -337,7 +332,7 @@ TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()
                      outputContainer->Add(fMixBOmega2[index]);\r
 \r
                      sprintf(key,"MixCToPi0Gamma2_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
-                     sprintf(title, "%s MixC Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%dcm",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,fBadChDist[idist]);\r
+                     sprintf(title, "%s MixC Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
                      fhEtalon->Clone(key);\r
                      fMixCOmega2[index]=(TH2F*)fhEtalon->Clone(key) ;\r
                      fMixCOmega2[index]->SetName(key) ;\r
@@ -491,7 +486,6 @@ void AliAnaOmegaToPi0Gamma::MakeAnalysisFillHistograms()
              else continue;\r
          }\r
     }\r
-\r
     //caculate the asy and dist of the two photon from pi0 decay\r
     TLorentzVector dph1(photon1->Px(),photon1->Py(),photon1->Pz(),photon1->E());\r
     TLorentzVector dph2(photon2->Px(),photon2->Py(),photon2->Pz(),photon2->E());\r
@@ -544,9 +538,9 @@ void AliAnaOmegaToPi0Gamma::MakeAnalysisFillHistograms()
                        if(photon1->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
                           photon2->IsPIDOK(ipid,AliCaloPID::kPhoton) && \r
                           photon3->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
-                          photon1->DistToBad()>=fBadChDist[idist] &&\r
-                          photon2->DistToBad()>=fBadChDist[idist] &&\r
-                          photon3->DistToBad()>=fBadChDist[idist] ){\r
+                          photon1->DistToBad()>=idist &&\r
+                          photon2->DistToBad()>=idist &&\r
+                          photon3->DistToBad()>=idist ){\r
                           //fill the histograms\r
                            if(GetDebug() > 2) printf("Real: index  %d  pt  %2.3f  mass   %2.3f \n", index, pi0gammapt, pi0gammamass);\r
                           fRealOmega0[index]->Fill(pi0gammapt,pi0gammamass);\r
@@ -585,9 +579,9 @@ void AliAnaOmegaToPi0Gamma::MakeAnalysisFillHistograms()
                        if(photon1->IsPIDOK(ipid,AliCaloPID::kPhoton)&&\r
                           photon2->IsPIDOK(ipid,AliCaloPID::kPhoton)&&\r
                           mix1ph->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
-                          photon1->DistToBad()>=fBadChDist[idist] &&\r
-                          photon2->DistToBad()>=fBadChDist[idist] &&\r
-                          mix1ph->DistToBad()>=fBadChDist[idist] ){\r
+                          photon1->DistToBad()>=idist &&\r
+                          photon2->DistToBad()>=idist &&\r
+                          mix1ph->DistToBad()>=idist ){\r
                           if(GetDebug() > 2) printf("MixA: index  %d   pt  %2.3f  mass   %2.3f \n",index, pi0gammapt, pi0gammamass);\r
                          //fill the histograms\r
                          fMixAOmega0[index]->Fill(pi0gammapt,pi0gammamass);\r
@@ -639,9 +633,9 @@ void AliAnaOmegaToPi0Gamma::MakeAnalysisFillHistograms()
                            if(ph1->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
                               ph2->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
                                ph3->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
-                              ph1->DistToBad()>=fBadChDist[idist] &&\r
-                               ph2->DistToBad()>=fBadChDist[idist] &&\r
-                               ph3->DistToBad()>=fBadChDist[idist] ){\r
+                              ph1->DistToBad()>=idist &&\r
+                               ph2->DistToBad()>=idist &&\r
+                               ph3->DistToBad()>=idist ){\r
                                if(GetDebug() > 2) printf("MixB: index  %d   pt  %2.3f  mass   %2.3f \n", index, pi0gammapt, pi0gammamass);\r
                               //fill histograms\r
                               fMixBOmega0[index]->Fill(pi0gammapt,pi0gammamass);\r
@@ -677,9 +671,9 @@ void AliAnaOmegaToPi0Gamma::MakeAnalysisFillHistograms()
                                 if(ph1->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
                                    ph2->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
                                    ph3->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
-                                  ph1->DistToBad()>=fBadChDist[idist] &&\r
-                                   ph2->DistToBad()>=fBadChDist[idist] &&\r
-                                   ph3->DistToBad()>=fBadChDist[idist] ){\r
+                                  ph1->DistToBad()>=idist &&\r
+                                   ph2->DistToBad()>=idist &&\r
+                                   ph3->DistToBad()>=idist ){\r
                                    if(GetDebug() > 2) printf("MixC: index  %d  pt  %2.3f  mass   %2.3f \n", index, pi0gammapt, pi0gammamass);\r
                                   //fill histograms\r
                                  fMixCOmega0[index]->Fill(pi0gammapt,pi0gammamass);\r

diff --git a/PWG4/PartCorrDep/AliAnaOmegaToPi0Gamma.h b/PWG4/PartCorrDep/AliAnaOmegaToPi0Gamma.h
index 967e7cb0d85d2289124ffba72849df5a2f155dfa..13bd37ed232b66e8d9f4bf437244130f8d95ad00 100644 (file)
--- a/PWG4/PartCorrDep/AliAnaOmegaToPi0Gamma.h
+++ b/PWG4/PartCorrDep/AliAnaOmegaToPi0Gamma.h
@@ -42,7 +42,6 @@ class AliAnaOmegaToPi0Gamma : public AliAnaPartCorrBaseClass {
 
   void SetNPID(Int_t pid) {fNpid=pid;} //N pid cut 
   void SetNVtxZ(Int_t vtx){fNVtxZBin=vtx;} //N vertex Z cut
-  void SetNBadChDist(Int_t dist){fNBadChDistBin=dist;} //N dist cut
   void SetPi0MassPeakWidthCut(Double_t win){fPi0MassWindow=win;} 
 
   void SetPi0OverOmegaPtCut(Double_t cut){fPi0OverOmegaPtCut=cut;}
@@ -66,7 +65,6 @@ class AliAnaOmegaToPi0Gamma : public AliAnaPartCorrBaseClass {
   Double_t *fVtxZCut;            //![fNVtxZBin] vtertx z cut
   Double_t *fCent;               //![fNCentBin] centrality cut
   Double_t *fRp;                 //![fNRpBin] reaction plane cut
-  Int_t *fBadChDist;             //![fNBadChDistBin] bad channel dist
   
   Double_t fPi0Mass;             //!nominal pi0 mass
   Double_t fPi0MassWindow;       //!pi0 mass windows

diff --git a/PWG4/macros/AddTaskPartCorr.C b/PWG4/macros/AddTaskPartCorr.C
index 42b91ad52b5b5e7cdd8bd6d83b3815b643f51f98..77a6eba567ec867889b29c0334675a3088ea9da6 100644 (file)
--- a/PWG4/macros/AddTaskPartCorr.C
+++ b/PWG4/macros/AddTaskPartCorr.C
@@ -476,14 +476,13 @@ AliAnalysisTaskParticleCorrelation *AddTaskPartCorr(TString data, TString calori
   anaomegaToPi0Gamma->SetDebug(-1);//10 for lots of messages
   anaomegaToPi0Gamma->SetInputAODName(Form("Pi0s%s",calorimeter.Data()));
   anaomegaToPi0Gamma->SetInputAODPhotonName(Form("Photons%s",calorimeter.Data()));
-  anaomegaToPi0Gamma->SetNPID(2);
-  anaomegaToPi0Gamma->SetNVtxZ(2);
-  anaomegaToPi0Gamma->SetNBadChDist(3);
+  anaomegaToPi0Gamma->SetNPID(1);
+  anaomegaToPi0Gamma->SetNVtxZ(1);
   anaomegaToPi0Gamma->SetNEventsMixed(4);
   if(calorimeter=="PHOS")
-           anaomegaToPi0Gamma->SetPi0MassPeakWidthCut(0.008); // PHOS
+    anaomegaToPi0Gamma->SetPi0MassPeakWidthCut(0.008); // PHOS
   else if(calorimeter=="EMCAL")
-           anaomegaToPi0Gamma->SetPi0MassPeakWidthCut(0.012); // EMCAL 
+    anaomegaToPi0Gamma->SetPi0MassPeakWidthCut(0.012); // EMCAL 
   anaomegaToPi0Gamma->SetHistoPtRangeAndNBins(0, 20, 200) ;
   anaomegaToPi0Gamma->SetHistoMassRangeAndNBins(0, 1, 200) ;
   anaomegaToPi0Gamma->SetPi0OverOmegaPtCut(0.8);
@@ -508,7 +507,7 @@ AliAnalysisTaskParticleCorrelation *AddTaskPartCorr(TString data, TString calori
   maker->AddAnalysis(anaisolpi0,8);
   maker->AddAnalysis(anacorrhadronpi0,9);
   maker->AddAnalysis(anacorrhadronisopi0,10);
-//  maker->AddAnalysis(anaomegaToPi0Gamma,11);   
+  maker->AddAnalysis(anaomegaToPi0Gamma,11);   
   maker->SetAnaDebug(-1)  ;
   maker->SwitchOnHistogramsMaker()  ;
   if(data.Contains("delta")) maker->SwitchOffAODsMaker()  ;