Allow to reject clusters with less than n cells
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 14 Apr 2010 09:29:31 +0000 (09:29 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 14 Apr 2010 09:29:31 +0000 (09:29 +0000)
PWG4/PartCorrDep/AliAnaPhoton.cxx
PWG4/PartCorrDep/AliAnaPhoton.h

index ddb102a..30d46fe 100755 (executable)
@@ -48,7 +48,7 @@ ClassImp(AliAnaPhoton)
     AliAnaPartCorrBaseClass(), fCalorimeter(""), 
     fMinDist(0.),fMinDist2(0.),fMinDist3(0.),fRejectTrackMatch(0),
        fCheckConversion(kFALSE),fAddConvertedPairsToAOD(kFALSE), fMassCut(0),
-       fTimeCutMin(-1), fTimeCutMax(9999999),
+    fTimeCutMin(-1), fTimeCutMax(9999999), fNCellsCut(0),
        fhPtPhoton(0),fhPhiPhoton(0),fhEtaPhoton(0),
     //MC
     fhDeltaE(0), fhDeltaPt(0),fhRatioE(0), fhRatioPt(0),fh2E(0),fh2Pt(0),
@@ -74,7 +74,7 @@ AliAnaPhoton::AliAnaPhoton(const AliAnaPhoton & g) :
   fMinDist(g.fMinDist),fMinDist2(g.fMinDist2), fMinDist3(g.fMinDist3),
   fRejectTrackMatch(g.fRejectTrackMatch),
   fCheckConversion(g.fCheckConversion),fAddConvertedPairsToAOD(g.fAddConvertedPairsToAOD),
-  fMassCut(g.fMassCut), fTimeCutMin(g.fTimeCutMin), fTimeCutMax(g.fTimeCutMax),        
+  fMassCut(g.fMassCut), fTimeCutMin(g.fTimeCutMin), fTimeCutMax(g.fTimeCutMax),        fNCellsCut(g.fNCellsCut),
   fhPtPhoton(g.fhPtPhoton),fhPhiPhoton(g.fhPhiPhoton),fhEtaPhoton(g.fhEtaPhoton), 
   //MC
   fhDeltaE(g.fhDeltaE), fhDeltaPt(g.fhDeltaPt), 
@@ -109,7 +109,8 @@ AliAnaPhoton & AliAnaPhoton::operator = (const AliAnaPhoton & g)
        
   fTimeCutMin = g.fTimeCutMin; 
   fTimeCutMax = g.fTimeCutMax; 
-
+  fNCellsCut  = g.fNCellsCut;
+       
   fRejectTrackMatch       = g.fRejectTrackMatch;
   fCheckConversion        = g.fCheckConversion;
   fAddConvertedPairsToAOD = g.fAddConvertedPairsToAOD;
@@ -429,7 +430,8 @@ void AliAnaPhoton::InitParameters()
        
   fTimeCutMin  = -1;
   fTimeCutMax  = 9999999;
-
+  fNCellsCut = 0;
+       
   fRejectTrackMatch       = kTRUE ;
   fCheckConversion        = kFALSE;
   fAddConvertedPairsToAOD = kFALSE;
@@ -483,8 +485,11 @@ void  AliAnaPhoton::MakeAnalysisFillAOD()
     //If too small or big pt, skip it
     if(mom.Pt() < GetMinPt() || mom.Pt() > GetMaxPt() ) continue ; 
     Double_t tof = calo->GetTOF()*1e9;
+         
        if(tof < fTimeCutMin || tof > fTimeCutMax) continue;
-
+         
+       if(calo->GetNCells() <= fNCellsCut) continue;
+         
        //printf("AliAnaPhoton::Current Event %d; Current File Name : %s, E %f, pT %f, Ecl %f\n",GetReader()->GetEventNumber(),(GetReader()->GetCurrentFileName()).Data(), mom.E(), mom.Pt(),calo->E());
 
     //Check acceptance selection
@@ -855,7 +860,7 @@ void AliAnaPhoton::Print(const Option_t * opt) const
   printf("Add conversion pair to AOD           = %d\n",fAddConvertedPairsToAOD);
   printf("Conversion pair mass cut             = %f\n",fMassCut);
   printf("Time Cut: %3.1f < TOF  < %3.1f\n", fTimeCutMin, fTimeCutMax);
-
+  printf("Number of cells in cluster is        > %d \n", fNCellsCut);
   printf("    \n") ;
        
 } 
index 6eb7c51..91280ac 100755 (executable)
@@ -76,6 +76,10 @@ public:
   Double_t GetTimeCutMin() const {return fTimeCutMin;}
   Double_t GetTimeCutMax() const {return fTimeCutMax;} 
        
+  void SetNCellCut(Int_t n) {fNCellsCut = n;}
+  Double_t GetNCellCut() const {return fNCellsCut;}
+       
+       
   private:
  
   TString fCalorimeter ; // Calorimeter where the gamma is searched;
@@ -88,7 +92,8 @@ public:
   Float_t fMassCut;                // Mass cut for the conversion pairs selection
   Double_t fTimeCutMin  ;    // Remove clusters/cells with time smaller than this value, in ns
   Double_t fTimeCutMax  ;    // Remove clusters/cells with time larger than this value, in ns
-
+  Int_t fNCellsCut ;     // Accept for the analysis clusters with more than fNCellsCut cells
+       
   //Histograms  
   TH1F * fhPtPhoton   ; //! Number of identified photon vs transerse momentum 
   TH2F * fhPhiPhoton  ; //! Azimuthal angle of identified  photon vs transerse momentum 
@@ -134,7 +139,7 @@ public:
   TH2F * fhPhiUnknown;  //! Phi of identified  Unknown gamma
   TH2F * fhEtaUnknown;  //! eta of identified  Unknown gamma
 
-   ClassDef(AliAnaPhoton,7)
+   ClassDef(AliAnaPhoton,8)
 
 } ;