Possibility of cut on distance to bad added
authorprsnko <Dmitri.Peressounko@cern.ch>
Sat, 27 Sep 2014 13:27:46 +0000 (17:27 +0400)
committermkrzewic <mikolaj.krzewicki@cern.ch>
Mon, 29 Sep 2014 08:47:54 +0000 (10:47 +0200)
PWGGA/PHOSTasks/PHOS_Tagging/AliAnalysisTaskTaggedPhotons.cxx
PWGGA/PHOSTasks/PHOS_Tagging/AliAnalysisTaskTaggedPhotons.h

index 015a44d..6526959 100644 (file)
@@ -73,6 +73,7 @@ AliAnalysisTaskTaggedPhotons::AliAnalysisTaskTaggedPhotons() :
   fZmin(0.),
   fPhimax(0.),
   fPhimin(0.),
+  fMinBCDistance(0.),
   fCentrality(0.),
   fCentBin(0), 
   fIsMB(0),
@@ -102,6 +103,7 @@ AliAnalysisTaskTaggedPhotons::AliAnalysisTaskTaggedPhotons(const char *name) :
   fZmin(60.),
   fPhimax(250.),
   fPhimin(320.),
+  fMinBCDistance(0.),
   fCentrality(0.),
   fCentBin(0),
   fIsMB(0),
@@ -137,6 +139,7 @@ AliAnalysisTaskTaggedPhotons::AliAnalysisTaskTaggedPhotons(const AliAnalysisTask
   fZmin(60.),
   fPhimax(250.),
   fPhimin(320.),
+  fMinBCDistance(0.),  
   fCentrality(0.),
   fCentBin(0),
   fIsMB(0),
@@ -727,6 +730,9 @@ void AliAnalysisTaskTaggedPhotons::UserExec(Option_t *)
   
     if(clu->GetMCEnergyFraction()>kEcrossCut) //Ecross cut, should be filled with Tender
      continue ;    
+    
+    if(clu->GetDistanceToBadChannel()<fMinBCDistance)
+      continue ;
 
     Float_t pos[3] ;
     clu->GetPosition(pos) ;
index c6b96ca..de3ca9f 100644 (file)
@@ -43,6 +43,7 @@ public:
   void SetMC(Bool_t isMC=kTRUE){fIsMC=isMC;}
   void SetFastMC(void){fIsFastMC=kTRUE;fIsMC=kTRUE; } //same as MC, but bypass event checks
   void SetPi0WeightParameters(TArrayD * ar) ;
+  void SetDustanceToBad(Float_t cut=2.5){fMinBCDistance=cut;}
 
 protected:
   void    FillMCHistos() ;
@@ -84,6 +85,7 @@ private:
   Float_t fZmin ;               //area
   Float_t fPhimax ;             //covered by
   Float_t fPhimin ;             //full calorimeter
+  Float_t fMinBCDistance;       //minimal distance to bad channel
   Double_t fWeightParamPi0[7] ; //Parameters to calculate weights
   //
   Double_t fCentrality;
@@ -93,6 +95,6 @@ private:
   Bool_t fIsFastMC; //This is fast MC, bypass event checks
   TH2I * fPHOSBadMap[6] ; 
     
-  ClassDef(AliAnalysisTaskTaggedPhotons, 2);   // a PHOS photon analysis task 
+  ClassDef(AliAnalysisTaskTaggedPhotons, 3);   // a PHOS photon analysis task 
 };
 #endif // ALIANALYSISTASKTAGGEDPHOTONS_H