Add some more cuts

diff --git a/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.cxx b/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.cxx
index d4d282ccab5535c61ebf7c35cc0d47d2e3bc307f..33eae5ba03d9fbbd57bf04c0b2f467c6537d513f 100644 (file)
--- a/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.cxx
+++ b/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.cxx
@@ -25,7 +25,6 @@ ClassImp(AliAnalysisTaskEMCALPi0PbPb)
 //________________________________________________________________________
 AliAnalysisTaskEMCALPi0PbPb::AliAnalysisTaskEMCALPi0PbPb(const char *name) 
   : AliAnalysisTaskSE(name),
-    fAsymMax(1),
     fCentVar("V0M"),
     fCentFrom(0),
     fCentTo(100),
@@ -35,7 +34,10 @@ AliAnalysisTaskEMCALPi0PbPb::AliAnalysisTaskEMCALPi0PbPb(const char *name)
     fClusName(),
     fDoNtuple(0),
     fDoAfterburner(0),
+    fAsymMax(1),
     fNminCells(1),
+    fMinErat(0),
+    fMinEcc(0),
     fNEvs(0),
     fGeom(0),
     fOutput(0),
@@ -393,7 +395,7 @@ void AliAnalysisTaskEMCALPi0PbPb::Terminate(Option_t *)
       fNtuple->Write();
   }
 
-  AliInfo(Form("\nAccepted %lld events", fNEvs));
+  AliInfo(Form("\n%s: Accepted %lld events", GetName(), fNEvs));
 }
 
 //________________________________________________________________________
@@ -469,7 +471,7 @@ void AliAnalysisTaskEMCALPi0PbPb::FillClusHists()
       GetSigma(clus,maxAxis,minAxis);
       if (maxAxis > 0)
         clusterEcc = TMath::Sqrt(1.0 - minAxis*minAxis/(maxAxis*maxAxis));
-
+      clus->SetChi2(clusterEcc); // store ecc in chi2
       fHClustEccentricity->Fill(clusterEcc); 
       fHClustEtaPhi->Fill(clusterVec.Eta(),clusterVec.Phi());
       fHClustEnergyPt->Fill(clusterVec.E(),clusterVec.Pt());
@@ -531,6 +533,10 @@ void AliAnalysisTaskEMCALPi0PbPb::FillPionHists()
         continue;
       if (clus1->GetNCells()<fNminCells)
         continue;
+      if (GetMaxCellEnergy(clus1)/clus1->E()<fMinErat)
+        continue;
+      if (clus1->Chi2()<fMinEcc)
+        continue;
       clus1->GetMomentum(clusterVec1,vertex);
       for (Int_t j = i+1; j<nclus; ++j) {
         AliVCluster *clus2 = static_cast<AliVCluster*>(clusters->At(j));
@@ -542,6 +548,10 @@ void AliAnalysisTaskEMCALPi0PbPb::FillPionHists()
           continue;
         if (clus2->GetNCells()<fNminCells)
           continue;
+        if (GetMaxCellEnergy(clus2)/clus2->E()<fMinErat)
+          continue;
+        if (clus2->Chi2()<fMinEcc)
+          continue;
         clus2->GetMomentum(clusterVec2,vertex);
         pionVec = clusterVec1 + clusterVec2;
         Double_t pionZgg = TMath::Abs(clusterVec1.E()-clusterVec2.E())/pionVec.E();

diff --git a/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.h b/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.h
index e6f86c25814e0846d4239b36fc9eaf91e7048f06..69cf27277ee2c234823602f74096524fe88cea94 100644 (file)
--- a/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.h
+++ b/PWG4/CaloCalib/AliAnalysisTaskEMCALPi0PbPb.h
@@ -34,6 +34,8 @@ class AliAnalysisTaskEMCALPi0PbPb : public AliAnalysisTaskSE {
   void         SetClusName(const char *name)                  { fClusName = name;           }
   void         SetDoAfterburner(Bool_t b)                     { fDoAfterburner = b;         }
   void         SetFillNtuple(Bool_t b)                        { fDoNtuple = b;              }
+  void         SetMinEcc(Double_t ecc)                        { fMinEcc = ecc;              }
+  void         SetMinErat(Double_t erat)                      { fMinErat = erat;            }
   void         SetNminCells(Int_t n)                          { fNminCells = n;             }
   void         SetUseQualFlag(Bool_t b)                       { fUseQualFlag = b;           }
   void         SetVertexRange(Double_t z1, Double_t z2)       { fVtxZMin=z1; fVtxZMax=z2;   }
@@ -48,7 +50,6 @@ class AliAnalysisTaskEMCALPi0PbPb : public AliAnalysisTaskSE {
   void         GetSigma(AliVCluster *c, Double_t &sigmaMax, Double_t &sigmaMin);
 
     // input members
-  Double_t               fAsymMax;                // energy asymmetry max (def=1)
   TString                fCentVar;                // variable for centrality determination
   Double_t               fCentFrom;               // min centrality (def=0)
   Double_t               fCentTo;                 // max centrality (def=100)
@@ -58,7 +59,10 @@ class AliAnalysisTaskEMCALPi0PbPb : public AliAnalysisTaskSE {
   TString                fClusName;               // cluster branch name (def="")
   Bool_t                 fDoNtuple;               // if true write out ntuple
   Bool_t                 fDoAfterburner;          // if true run after burner
-  Int_t                  fNminCells;              // minimum number of cells attached to cluster
+  Double_t               fAsymMax;                // maximum energy asymmetry (def=1)
+  Int_t                  fNminCells;              // minimum number of cells attached to cluster (def=1)
+  Double_t               fMinErat;                // minimum emax/ec ratio (def=0)
+  Double_t               fMinEcc;                 // minimum eccentricity (def=0)
     // derived members (ie with ! after //)
   ULong64_t              fNEvs;                   //!accepted events 
   AliEMCALGeoUtils      *fGeom;                   //!geometry utils