Added ability to enable/disable PHOS modules from analysis.

diff --git a/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0Flow.cxx b/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0Flow.cxx
index 9b624d2813c23bc15bd9a699321ddf2ca07f0f8f..eb0641f005a55581f16cd987dd9d3768e66e58fd 100644 (file)
--- a/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0Flow.cxx
+++ b/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0Flow.cxx
@@ -83,6 +83,7 @@ AliAnalysisTaskPi0Flow::AliAnalysisTaskPi0Flow(const char *name, Period period)
   fPeriod(period),
   fMaxAbsVertexZ(10.),
   fManualV0EPCalc(false),
+  fModuleEnabled({true}),
   fOutputContainer(0x0),
   fNonLinCorr(0),
   fEvent(0x0),
@@ -479,6 +480,13 @@ void AliAnalysisTaskPi0Flow::SetCentralityBinning(const TArrayD& edges, const TA
   fCentNMixed = nMixed;
 }
 
+//_____________________________________________________________________________
+void AliAnalysisTaskPi0Flow::SetEnablePHOSModule(int module, Bool_t enable)
+{
+  if( module < 1 || 5 < module )
+    AliFatal("PHOS Module must be between 1 and 5");
+  fModuleEnabled[module-1] = enable;
+}
 
 
 //________________________________________________________________________
@@ -599,6 +607,8 @@ void AliAnalysisTaskPi0Flow::SelectPhotonClusters()
     Int_t mod  = relId[0] ;
     Int_t cellX = relId[2];
     Int_t cellZ = relId[3] ;
+    if ( ! fModuleEnabled[mod-1] )
+      continue;
     if ( !IsGoodChannel("PHOS",mod,cellX,cellZ) )
       continue ; // reject if not.
 

diff --git a/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0Flow.h b/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0Flow.h
index 93873921408500d3a3e62b5c286b52ef40fd1abb..13c95efce692c97df8d4f19b628a428ae93aa79d 100644 (file)
--- a/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0Flow.h
+++ b/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0Flow.h
@@ -47,6 +47,7 @@ public:
     void SetEventMixingRPBinning(UInt_t nBins) { fNEMRPBins = nBins; }
     void SetMaxAbsVertexZ(Float_t z) { fMaxAbsVertexZ = z; }
     void SetManualV0EPCalc(Bool_t manCalc = true) {fManualV0EPCalc = manCalc;}
+    void SetEnablePHOSModule(int module, Bool_t enable = true);
     
     void SetPHOSBadMap(Int_t mod,TH2I * badMapHist);
     //Where to read AODB object with EP calibration if not default
@@ -158,6 +159,7 @@ private:
     Period fPeriod;
     Float_t fMaxAbsVertexZ; // in cm
     Bool_t fManualV0EPCalc;
+    Bool_t fModuleEnabled[5];
 
 
     TList * fOutputContainer;        //final histogram container