including option to apply a cut in the distance to bad channel of clusters; changing...

diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
index 90c9bdf3c97ca7ced49bff73d134597fdfcf9652..3d21d0c824529b09d38602fcdc4ca71ead18bbf6 100644 (file)
--- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
+++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
@@ -83,14 +83,15 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton() :
   fClusIdFromTracks(""),
   fCpvFromTrack(kFALSE),
   fNBinsPt(200),
-  fPtBinLowEdge(-0.25),
-  fPtBinHighEdge(99.75),
+  fPtBinLowEdge(0),
+  fPtBinHighEdge(100),
   fRemMatchClus(kFALSE),
   fMinIsoClusE(0),
   fNCuts(5),
   fTrCoreRem(kFALSE),
   fClusTDiff(30e-9),
   fPileUpRejSPD(kFALSE),
+  fDistToBadChan(0),
   fESD(0),
   fAOD(0),
   fVEvent(0),
@@ -191,14 +192,15 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton(const char *name) :
   fClusIdFromTracks(""),
   fCpvFromTrack(kFALSE),
   fNBinsPt(200),
-  fPtBinLowEdge(-0.25),
-  fPtBinHighEdge(99.75),
+  fPtBinLowEdge(0.),
+  fPtBinHighEdge(100),
   fRemMatchClus(kFALSE),
   fMinIsoClusE(0),
   fNCuts(5),
   fTrCoreRem(kFALSE),
   fClusTDiff(30e-9),
   fPileUpRejSPD(kFALSE),
+  fDistToBadChan(0),
   fESD(0),
   fAOD(0),
   fVEvent(0),
@@ -730,6 +732,11 @@ void AliAnalysisTaskEMCALIsoPhoton::FillClusHists()
        printf("cluster is exotic! xxxx\n");
       continue;
     }
+    if(c->GetDistanceToBadChannel()<fDistToBadChan){
+      if(fDebug)
+       printf("cluster distance to bad channel is %1.1f (<%1.1f) xxxx\n",c->GetDistanceToBadChannel(),fDistToBadChan);
+      continue;
+    }
     Short_t id;
     Double_t Emax = GetMaxCellEnergy( c, id);
     if(fDebug)

diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h
index 4df15e4e1c4c232d9c66d10532874d733db203d7..2d304cdc842ca02881dac8579633f264e447c6ae 100644 (file)
--- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h
+++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h
@@ -75,7 +75,8 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE {
   void                  SetMinIsoClusE(Double_t emin)           { fMinIsoClusE        = emin;    }
   void                  SetTrCoreRemoval(Bool_t b)              { fTrCoreRem          = b;       }
   void                  SetClusTDiff(Double_t diff)             { fClusTDiff          = diff;    }
-  void                  SetPileUpRejSPD()                       { fPileUpRejSPD       = kTRUE;  }
+  void                  SetPileUpRejSPD()                       { fPileUpRejSPD       = kTRUE;   }
+  void                  SetDistanceToBadCh(Double_t d)          { fDistToBadChan      = d;       }
  protected:
   TObjArray             *fESDClusters;           //!pointer to EMCal clusters
   TObjArray             *fAODClusters;           //!pointer to EMCal clusters
@@ -123,6 +124,7 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE {
   Bool_t                 fTrCoreRem;             // flag to set the removal of the core in track isolation (true removes it, default)
   Double_t               fClusTDiff;             // variable to hold the time diff between the candidate cluster and the isolation clusters
   Bool_t                 fPileUpRejSPD;          // flag to set pile-up rejection via SPD (multiple vertices)
+  Double_t               fDistToBadChan;         // distance to bad channel
   
  private:
   AliESDEvent *fESD;      //! ESD object

diff --git a/PWGGA/EMCALTasks/macros/AddTaskEMCALIsoPhoton.C b/PWGGA/EMCALTasks/macros/AddTaskEMCALIsoPhoton.C
index 9c1dea0a852bb796b007491a2356853781f6ac5d..9d56a66e34a72b9fdce935b124773e062f636019 100644 (file)
--- a/PWGGA/EMCALTasks/macros/AddTaskEMCALIsoPhoton.C
+++ b/PWGGA/EMCALTasks/macros/AddTaskEMCALIsoPhoton.C
@@ -5,7 +5,8 @@ AliAnalysisTaskEMCALIsoPhoton *AddTaskEMCALIsoPhoton(
                                                     TString trigbitname = "kEMC7",
                                                     TString geoname="EMCAL_COMPLETEV1",
                                                     TString pathstrsel = "/",
-                                                    TString trackSelType = "standard"
+                                                    TString trackSelType = "standard",
+                                                    Int_t   distToBadCh = 0
                                                     )
 {
   // Get the pointer to the existing analysis manager via the static access method.
@@ -62,6 +63,7 @@ AliAnalysisTaskEMCALIsoPhoton *AddTaskEMCALIsoPhoton(
   ana->SetPeriod(period.Data());
   ana->SetGeoName(geoname.Data());  
   //ana->SetTrackFilterBit(128);
+  ana->SetDistanceToBadCh(distToBadCh);
   mgr->AddTask(ana);
   TString containername = "histEMCIsoPhoton."+trigbitname;
   TString containernameQA = "histosQA."+trigbitname;