]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG4/PartCorrBase/AliCaloTrackReader.h
protection agains null clusters
[u/mrichter/AliRoot.git] / PWG4 / PartCorrBase / AliCaloTrackReader.h
index 372f6171fd6d4e91a6af38f6119f4f54b8ce3115..5aac1f933e97a9662a4a61c3fcb7eefc16710978 100755 (executable)
 // Mother class of : AliCaloTrackESDReader: Fills ESD data in 3 TObjArrays (PHOS, EMCAL, CTS)
 //                 : AliCaloTrackMCReader: Fills Kinematics data in 3 TObjArrays (PHOS, EMCAL, CTS)
 //                 : AliCaloTrackAODReader: Fills AOD data in 3 TObjArrays (PHOS, EMCAL, CTS) 
-//                          
+//  
+// This part is commented: Mixing analysis can be done, input AOD with events
+// is opened in the AliCaloTrackReader::Init()
+
 // -- Author: Gustavo Conesa (INFN-LNF)
 
 // --- ROOT system ---
 #include "TObject.h" 
 class TObjArray ; 
-class TLorentzVector ;
 #include "TString.h"
-#include "TObjArray.h"
-class TArrayF;  
 class TTree ;
-#include "TH2I.h"
 
 //--- ANALYSIS system ---
 class AliStack ; 
@@ -31,22 +30,26 @@ class AliGenEventHeader ;
 class AliVEvent;
 class AliAODEvent;
 class AliMCEvent;
-class AliFiducialCut;
+class AliMixedEvent;
+#include "AliVCaloCells.h"
+#include "AliFiducialCut.h"
 class AliAODMCHeader;
-#include "AliPHOSGeoUtils.h"
-#include "AliEMCALGeoUtils.h"
+#include "AliCalorimeterUtils.h"
+class AliESDtrackCuts;
+class AliESDCentrality;
 
 class AliCaloTrackReader : public TObject {
 
- public: 
-  
+public: 
   AliCaloTrackReader() ; // ctor
+  virtual ~AliCaloTrackReader() ;//virtual dtor
+private:
   AliCaloTrackReader(const AliCaloTrackReader & g) ; // cpy ctor
   AliCaloTrackReader & operator = (const AliCaloTrackReader & g) ;//cpy assignment
-  virtual ~AliCaloTrackReader() ;//virtual dtor
 
+public:
   enum inputDataType {kESD, kAOD, kMC};
-  
+    
   //Select generated events, depending on comparison of pT hard and jets.
   virtual Bool_t ComparePtHardAndJetPt() ;
   virtual Bool_t IsPtHardAndJetPtComparisonSet() const {return  fComparePtHardAndJetPt ;}
@@ -75,39 +78,45 @@ class AliCaloTrackReader : public TObject {
   virtual void SetCTSPtMin(Float_t  pt)   { fCTSPtMin = pt ; }
   
   //Input setters and getters
-
+  Bool_t IsEMCALCluster(AliVCluster *clus) const;
+  Bool_t IsPHOSCluster (AliVCluster *clus)  const;
+  void SwitchOnOldAODs()   {fOldAOD = kTRUE  ; }
+  void SwitchOffOldAODs()  {fOldAOD = kFALSE ; }
+  
   Bool_t IsCTSSwitchedOn()  const { return fFillCTS ; }
-  void SwitchOnCTS()    {fFillCTS = kTRUE ; }
+  void SwitchOnCTS()    {fFillCTS = kTRUE  ; }
   void SwitchOffCTS()   {fFillCTS = kFALSE ; }
 
   Bool_t IsEMCALSwitchedOn() const { return fFillEMCAL ; }
-  void SwitchOnEMCAL()  {fFillEMCAL = kTRUE ; }
+  void SwitchOnEMCAL()  {fFillEMCAL = kTRUE  ; }
   void SwitchOffEMCAL() {fFillEMCAL = kFALSE ; }
 
   Bool_t IsPHOSSwitchedOn()  const { return fFillPHOS ; }
-  void SwitchOnPHOS()   {fFillPHOS = kTRUE ; }
+  void SwitchOnPHOS()   {fFillPHOS = kTRUE  ; }
   void SwitchOffPHOS()  {fFillPHOS = kFALSE ; }
 
   Bool_t IsEMCALCellsSwitchedOn() const { return fFillEMCALCells ; }
-  void SwitchOnEMCALCells()  {fFillEMCALCells = kTRUE ; }
+  void SwitchOnEMCALCells()  {fFillEMCALCells = kTRUE  ; }
   void SwitchOffEMCALCells() {fFillEMCALCells = kFALSE ; }
 
   Bool_t IsPHOSCellsSwitchedOn()  const { return fFillPHOSCells ; }
-  void SwitchOnPHOSCells()   {fFillPHOSCells = kTRUE ; }
+  void SwitchOnPHOSCells()   {fFillPHOSCells = kTRUE  ; }
   void SwitchOffPHOSCells()  {fFillPHOSCells = kFALSE ; }
 
   virtual Bool_t FillInputEvent(const Int_t iEntry, const char *currentFileName)  ;
-  virtual void FillInputCTS()   {;}
-  virtual void FillInputEMCAL() {;}
-  virtual void FillInputPHOS()  {;}
-  virtual void FillInputEMCALCells() {;}
-  virtual void FillInputPHOSCells()  {;}
+  virtual void FillInputCTS() ;
+  virtual void FillInputEMCAL() ;
+  virtual void FillInputPHOS() ;
+  virtual void FillInputEMCALCells() ;
+  virtual void FillInputPHOSCells() ;
+  
+  virtual TList * GetAODBranchList() const { return fAODBranchList ; }
 
-  virtual TObjArray* GetAODCTS()   const {return fAODCTS ;}
+  virtual TObjArray* GetAODCTS()   const {return fAODCTS   ;}
   virtual TObjArray* GetAODEMCAL() const {return fAODEMCAL ;}
-  virtual TObjArray* GetAODPHOS()  const {return fAODPHOS ;}
-  virtual TNamed* GetEMCALCells()  const {return fEMCALCells ;}
-  virtual TNamed* GetPHOSCells()   const {return fPHOSCells ;}
+  virtual TObjArray* GetAODPHOS()  const {return fAODPHOS  ;}
+  virtual AliVCaloCells* GetEMCALCells()  const { return fEMCALCells ;}
+  virtual AliVCaloCells* GetPHOSCells()   const { return fPHOSCells  ;}
 
   //Get MC  informatio
   //Kinematics and galice.root available 
@@ -118,49 +127,63 @@ class AliCaloTrackReader : public TObject {
   virtual TClonesArray*   GetAODMCParticles(Int_t input = 0) const ;
   virtual AliAODMCHeader* GetAODMCHeader(Int_t input = 0)    const ;
        
-  virtual AliVEvent*   GetInputEvent()  const {return fInputEvent;}
-  virtual AliAODEvent* GetOutputEvent() const {return fOutputEvent;}
-  virtual AliMCEvent*  GetMC()          const {return fMC;}
-  virtual void         GetVertex(Double_t *) const {;}
-  virtual void         GetSecondInputAODVertex(Double_t *) const {;}
-  virtual Double_t     GetBField() const { return 0.;}
+  virtual AliVEvent*     GetInputEvent()  const {return fInputEvent;}
+  virtual AliAODEvent*   GetOutputEvent() const {return fOutputEvent;}
+  virtual AliMCEvent*    GetMC()          const {return fMC;}
+  virtual AliMixedEvent* GetMixedEvent()  const {return fMixedEvent;}
+  virtual Int_t          GetNMixedEvent() const {return fNMixedEvent ; } 
+
+  virtual Double_t       GetBField() const { return 0.;}
        
   virtual void Init();
-       
-  virtual void SetInputEvent(AliVEvent* const input)  {fInputEvent  = input;}
+
+//  virtual void SetInputEvent(AliVEvent* const input)  {fInputEvent  = input;}
+  virtual void SetInputEvent(AliVEvent* const input) ;
   virtual void SetOutputEvent(AliAODEvent* const aod) {fOutputEvent = aod;}
   virtual void SetMC(AliMCEvent* const mc)            {fMC  = mc;}
 
   virtual void ResetLists();
 
-  virtual AliFiducialCut * GetFiducialCut() const {return  fFiducialCut ;}
+  virtual AliFiducialCut * GetFiducialCut() {if(!fFiducialCut) fFiducialCut = new AliFiducialCut(); 
+         return  fFiducialCut ;}
   virtual void SetFiducialCut(AliFiducialCut * const fc) { fFiducialCut = fc ;}
+  virtual Bool_t IsFiducialCutOn()       {return fCheckFidCut ; }
+  virtual void SwitchOnFiducialCut()     { fCheckFidCut = kTRUE;  fFiducialCut = new AliFiducialCut();}
+  virtual void SwitchOffFiducialCut()    { fCheckFidCut = kFALSE;}
        
   virtual void SetInputOutputMCEvent(AliVEvent* /*esd*/, AliAODEvent* /*aod*/, AliMCEvent* /*mc*/) {;}
        
   //Methods for mixing with external input file (AOD)
-  virtual TTree* GetSecondInputAODTree() const {return  fSecondInputAODTree ; } 
+  //virtual TTree* GetSecondInputAODTree() const {return  fSecondInputAODTree ; } 
   //virtual void SetSecondInputAODTree(TTree * tree) {fSecondInputAODTree = tree ;
   //                                                                                             fSecondInputAODEvent->ReadFromTree(tree);}//Connect tree and AOD event.
                                        
-  virtual AliAODEvent* GetSecondInputAODEvent() const { return fSecondInputAODEvent ; } 
+  //virtual AliAODEvent* GetSecondInputAODEvent() const { return fSecondInputAODEvent ; } 
        
-  TString GetSecondInputFileName() const    {return fSecondInputFileName ; }
-  void SetSecondInputFileName(TString name) { fSecondInputFileName = name ; }
+  //TString GetSecondInputFileName() const    {return fSecondInputFileName ; }
+  //void SetSecondInputFileName(TString name) { fSecondInputFileName = name ; }
 
-  Int_t GetSecondInputFirstEvent() const    {return fSecondInputFirstEvent ; }
-  void SetSecondInputFirstEvent(Int_t iEvent0) { fSecondInputFirstEvent = iEvent0 ; }  
+  //Int_t GetSecondInputFirstEvent() const    {return fSecondInputFirstEvent ; }
+  //void SetSecondInputFirstEvent(Int_t iEvent0) { fSecondInputFirstEvent = iEvent0 ; }        
        
-  Int_t GetAODCTSNormalInputEntries()   {if(!fSecondInputAODTree) { fAODCTSNormalInputEntries   = fAODCTS->GetEntriesFast()  ;}
-                                                                                return fAODCTSNormalInputEntries ; }
-  Int_t GetAODEMCALNormalInputEntries() {if(!fSecondInputAODTree) { fAODEMCALNormalInputEntries = fAODEMCAL->GetEntriesFast();}
-                                                                                return fAODEMCALNormalInputEntries ; }
-  Int_t GetAODPHOSNormalInputEntries()  {if(!fSecondInputAODTree) { fAODPHOSNormalInputEntries  = fAODPHOS->GetEntriesFast() ;}
-                                                                                return fAODPHOSNormalInputEntries ; }
+//  Int_t GetAODCTSNormalInputEntries()   {if(!fSecondInputAODTree) { fAODCTSNormalInputEntries   = fAODCTS->GetEntriesFast()  ;}
+//                                                                              return fAODCTSNormalInputEntries ; }
+//  Int_t GetAODEMCALNormalInputEntries() {if(!fSecondInputAODTree) { fAODEMCALNormalInputEntries = fAODEMCAL->GetEntriesFast();}
+//                                                                              return fAODEMCALNormalInputEntries ; }
+//  Int_t GetAODPHOSNormalInputEntries()  {if(!fSecondInputAODTree) { fAODPHOSNormalInputEntries  = fAODPHOS->GetEntriesFast() ;}
+//                                                                              return fAODPHOSNormalInputEntries ; }
        
-  ULong_t GetTrackStatus() const    {return fTrackStatus ; }
+  // Track selection
+  ULong_t GetTrackStatus() const   {return fTrackStatus ; }
   void SetTrackStatus(ULong_t bit) { fTrackStatus = bit ; }            
        
+  AliESDtrackCuts* GetTrackCuts()          const  { return fESDtrackCuts    ; }
+  void    SetTrackCuts(AliESDtrackCuts * cuts)    { fESDtrackCuts = cuts    ; }                  
+  Int_t   GetTrackMultiplicity()           const  { return fTrackMult       ; }
+  Float_t GetTrackMultiplicityEtaCut()     const  { return fTrackMultEtaCut ; }
+  void    SetTrackMultiplicityEtaCut(Float_t eta) { fTrackMultEtaCut = eta  ; }                
+  
+  //MC switchs
   void SwitchOnStack()              { fReadStack          = kTRUE  ; }
   void SwitchOffStack()             { fReadStack          = kFALSE ; }
   void SwitchOnAODMCParticles()     { fReadAODMCParticles = kTRUE  ; }
@@ -168,111 +191,132 @@ class AliCaloTrackReader : public TObject {
   Bool_t ReadStack()          const { return fReadStack            ; }
   Bool_t ReadAODMCParticles() const { return fReadAODMCParticles   ; }
        
-  void SwitchOnCleanStdAOD()  {fCleanOutputStdAOD = kTRUE;}
-  void SwitchOffCleanStdAOD() {fCleanOutputStdAOD = kFALSE;}
-       
   void SetDeltaAODFileName(TString name ) {fDeltaAODFileName = name ; }
   TString GetDeltaAODFileName() const {return fDeltaAODFileName ; }
 
   void SetFiredTriggerClassName(TString name ) {fFiredTriggerClassName = name ; }
   TString GetFiredTriggerClassName() const {return fFiredTriggerClassName ; }
   virtual TString GetFiredTriggerClasses() {return "";}
-       
-  //Calorimeters Geometry Methods
-  void SetEMCALGeometryName(TString name)   { fEMCALGeoName = name ; }
-  TString EMCALGeometryName() const { return fEMCALGeoName ; }
-  void InitEMCALGeometry() ; 
-  AliEMCALGeoUtils * GetEMCALGeometry() const {return fEMCALGeo;}
-  Bool_t IsEMCALGeoMatrixSet() {return fEMCALGeoMatrixSet; }
-
-  void SetPHOSGeometryName(TString name)   { fPHOSGeoName = name ; }
-  TString PHOSGeometryName() const { return fPHOSGeoName ; }
-  void InitPHOSGeometry() ; 
-  AliPHOSGeoUtils * GetPHOSGeometry() const {return fPHOSGeo;} 
-  Bool_t IsPHOSGeoMatrixSet()  {return fPHOSGeoMatrixSet ; }
 
   void AnalyzeOnlyLED()     {fAnaLED = kTRUE;}
   void AnalyzeOnlyPhysics() {fAnaLED = kFALSE;}
        
-  // Bad channels
-       
-  Bool_t IsBadChannelsRemovalSwitchedOn()  const { return fRemoveBadChannels ; }
-  void SwitchOnBadChannelsRemoval()    {fRemoveBadChannels = kTRUE ; }
-  void SwitchOffBadChannelsRemoval()   {fRemoveBadChannels = kFALSE ; }
-       
-  Int_t GetEMCALChannelStatus(Int_t iSM , Int_t iCol, Int_t iRow) const { return (Int_t) ((TH2I*)fEMCALBadChannelMap[iSM])->GetBinContent(iCol,iRow);}
-  Int_t GetPHOSChannelStatus (Int_t imod, Int_t iCol, Int_t iRow) const { return (Int_t) ((TH2I*)fPHOSBadChannelMap[imod])->GetBinContent(iCol,iRow);}
-       
-  void SetEMCALChannelStatus(Int_t iSM , Int_t iCol, Int_t iRow, Double_t c = 1) { ((TH2I*)fEMCALBadChannelMap[iSM])->SetBinContent(iCol,iRow,c);}
-  void SetPHOSChannelStatus (Int_t imod, Int_t iCol, Int_t iRow, Double_t c = 1) { ((TH2I*)fPHOSBadChannelMap[imod])->SetBinContent(iCol,iRow,c);}
-       
-  TH2I * GetEMCALChannelStatusMap(Int_t iSM) const {return (TH2I*)fEMCALBadChannelMap[iSM];}
-  TH2I * GetPHOSChannelStatusMap(Int_t imod) const {return (TH2I*)fPHOSBadChannelMap[imod];}
+  TString  GetTaskName() const {return fTaskName;}
+  void SetTaskName(TString name) {fTaskName = name;}
 
-  void SetEMCALChannelStatusMap(TObjArray map) {fEMCALBadChannelMap = map;}
-  void SetPHOSChannelStatusMap (TObjArray map) {fPHOSBadChannelMap  = map;}
-       
-  Bool_t ClusterContainsBadChannel(TString calorimeter,UShort_t* cellList, Int_t nCells);
-       
+  AliCalorimeterUtils * GetCaloUtils() const {return fCaloUtils ; }
+  void SetCaloUtils(AliCalorimeterUtils * caloutils) { fCaloUtils = caloutils ; }
+  
+  virtual void GetVertex(Double_t v[3]) const ;
+  virtual Double_t* GetVertex(const Int_t evtIndex) const {return fVertex[evtIndex];}
+  virtual void GetVertex(Double_t vertex[3], const Int_t evtIndex) const ;
+  virtual void FillVertexArray();
+ // virtual void       GetSecondInputAODVertex(Double_t *) const {;}
+  
+  void SwitchOnWriteDeltaAOD()  {fWriteOutputDeltaAOD = kTRUE ;  }
+  void SwitchOffWriteDeltaAOD() {fWriteOutputDeltaAOD = kFALSE ; }
+  Bool_t WriteDeltaAODToFile() const {return fWriteOutputDeltaAOD ; } 
+  
+  virtual void FillInputVZERO(){;}
+  Int_t GetV0Signal(Int_t i) const { return fV0ADC[i];}
+  Int_t GetV0Multiplicity(Int_t i)   const { return fV0Mul[i];}
+
+  void SwitchOnCaloFilterPatch()  {fCaloFilterPatch = kTRUE ; fFillCTS = kFALSE; }
+  void SwitchOffCaloFilterPatch() {fCaloFilterPatch = kFALSE ; }
+  Bool_t IsCaloFilterPatchOn()    {if(fDataType == kAOD) { return fCaloFilterPatch ; } 
+                                   else                  { return kFALSE ; } }
+  
+  virtual AliESDCentrality* GetCentrality() const {return 0x0;}
+
+  
+  //MC reader methods:
+  
+  virtual void AddNeutralParticlesArray(TArrayI & /*array*/)  { ; }  
+  virtual void AddChargedParticlesArray(TArrayI & /*array*/)  { ; } 
+  virtual void AddStatusArray(TArrayI & /*array*/)            { ; }
+  
+  virtual void SwitchOnPi0Decay()         { ; } 
+  virtual void SwitchOffPi0Decay()        { ; } 
+  virtual void SwitchOnStatusSelection()  { ; }
+  virtual void SwitchOffStatusSelection() { ; }
+  virtual void SwitchOnOverlapCheck()     { ; }
+  virtual void SwitchOffOverlapCheck()    { ; }
+  
+  virtual void SetEMCALOverlapAngle(Float_t /*angle*/)  { ; }
+  virtual void SetPHOSOverlapAngle(Float_t /*angle*/)   { ; }
+  
+
+  
  protected:
-  Int_t                   fEventNumber; // Event number
-  TString          fCurrentFileName; // Current file name under analysis
-  Int_t            fDataType ;   // Select MC:Kinematics, Data:ESD/AOD, MCData:Both
-  Int_t            fDebug;       // Debugging level
-  AliFiducialCut * fFiducialCut; //! Acceptance cuts
-       
+  Int_t                   fEventNumber;    // Event number
+  TString          fCurrentFileName;// Current file name under analysis
+  Int_t            fDataType ;      // Select MC:Kinematics, Data:ESD/AOD, MCData:Both
+  Int_t            fDebug;          // Debugging level
+  AliFiducialCut * fFiducialCut;    //! Acceptance cuts
+  Bool_t           fCheckFidCut ;   // Do analysis for clusters in defined region         
+
   Bool_t           fComparePtHardAndJetPt;  // In MonteCarlo, jet events, reject fake events with wrong jet energy.
   Float_t          fPtHardAndJetPtFactor;   // Factor between ptHard and jet pT to reject/accept event.
 
-  Float_t        fCTSPtMin;      // pT Threshold on charged particles 
-  Float_t        fEMCALPtMin;    // pT Threshold on emcal clusters
-  Float_t        fPHOSPtMin;     // pT Threshold on phos clusters
-
-  TObjArray *    fAODCTS ;        //! temporal referenced array with tracks
-  TObjArray *    fAODEMCAL ;      //! temporal referenced array with EMCAL CaloClusters
-  TObjArray *    fAODPHOS ;       //! temporal referenced array with PHOS CaloClusters
-  TNamed *       fEMCALCells ;    //! temporal array with EMCAL CaloCells, ESD or AOD
-  TNamed *       fPHOSCells ;     //! temporal array with PHOS CaloCells, ESD or AOD
-
-  AliVEvent   *  fInputEvent;     //! pointer to esd or aod input
-  AliAODEvent *  fOutputEvent;    //! pointer to aod output
-  AliMCEvent  *  fMC;             //! Monte Carlo Event Handler  
-
-  Bool_t         fFillCTS;        // use data from CTS
-  Bool_t         fFillEMCAL;      // use data from EMCAL
-  Bool_t         fFillPHOS;       // use data from PHOS
-  Bool_t         fFillEMCALCells; // use data from EMCAL
-  Bool_t         fFillPHOSCells;  // use data from PHOS
-
-  TTree *        fSecondInputAODTree;    // Tree with second input AOD, for mixing analysis.   
-  AliAODEvent*   fSecondInputAODEvent;   //! pointer to second input AOD event.
-  TString        fSecondInputFileName;   // File with AOD data to mix with normal stream of data.
-  Int_t          fSecondInputFirstEvent; // First event to be considered in the mixing.
+  Float_t          fCTSPtMin;       // pT Threshold on charged particles 
+  Float_t          fEMCALPtMin;     // pT Threshold on emcal clusters
+  Float_t          fPHOSPtMin;      // pT Threshold on phos clusters
+
+  TList          * fAODBranchList ; //! List with AOD branches created and needed in analysis  
+  TObjArray      * fAODCTS ;        //! temporal referenced array with tracks
+  TObjArray      * fAODEMCAL ;      //! temporal referenced array with EMCAL CaloClusters
+  TObjArray      * fAODPHOS ;       //! temporal referenced array with PHOS CaloClusters
+  AliVCaloCells  * fEMCALCells ;    //! temporal array with EMCAL CaloCells, ESD or AOD
+  AliVCaloCells  * fPHOSCells ;     //! temporal array with PHOS CaloCells, ESD or AOD
+
+  AliVEvent      * fInputEvent;     //! pointer to esd or aod input
+  AliAODEvent    * fOutputEvent;    //! pointer to aod output
+  AliMCEvent     * fMC;             //! Monte Carlo Event Handler  
+
+  Bool_t           fFillCTS;        // use data from CTS
+  Bool_t           fFillEMCAL;      // use data from EMCAL
+  Bool_t           fFillPHOS;       // use data from PHOS
+  Bool_t           fFillEMCALCells; // use data from EMCAL
+  Bool_t           fFillPHOSCells;  // use data from PHOS
+
+//  TTree *        fSecondInputAODTree;    // Tree with second input AOD, for mixing analysis. 
+//  AliAODEvent*   fSecondInputAODEvent;   //! pointer to second input AOD event.
+//  TString        fSecondInputFileName;   // File with AOD data to mix with normal stream of data.
+//  Int_t          fSecondInputFirstEvent; // First event to be considered in the mixing.
+//     
+//  Int_t          fAODCTSNormalInputEntries;   // Number of entries in CTS   in case of standard input, larger with mixing.
+//  Int_t          fAODEMCALNormalInputEntries; // Number of entries in EMCAL in case of standard input, larger with mixing.
+//  Int_t          fAODPHOSNormalInputEntries;  // Number of entries in PHOS  in case of standard input, larger with mixing.
        
-  Int_t          fAODCTSNormalInputEntries;   // Number of entries in CTS   in case of standard input, larger with mixing.
-  Int_t          fAODEMCALNormalInputEntries; // Number of entries in EMCAL in case of standard input, larger with mixing.
-  Int_t          fAODPHOSNormalInputEntries;  // Number of entries in PHOS  in case of standard input, larger with mixing.
+  ULong_t          fTrackStatus        ; // Track selection bit, select tracks refitted in TPC, ITS ...
+  AliESDtrackCuts *fESDtrackCuts       ; // Track cut  
+  Int_t            fTrackMult          ; // Track multiplicity
+  Float_t          fTrackMultEtaCut    ; // Track multiplicity eta cut
+  Bool_t           fReadStack          ; // Access kine information from stack
+  Bool_t                fReadAODMCParticles ; // Access kine information from filtered AOD MC particles
        
-  ULong_t        fTrackStatus        ; // Track selection bit, select tracks refitted in TPC, ITS ...
-  Bool_t         fReadStack          ; // Access kine information from stack
-  Bool_t            fReadAODMCParticles ; // Access kine information from filtered AOD MC particles
+  TString          fDeltaAODFileName   ; // Delta AOD file name
+  TString          fFiredTriggerClassName; // Name of trigger event type used to do the analysis
+
+  Bool_t           fAnaLED;             // Analyze LED data only.
+
+  TString          fTaskName;           // Name of task that executes the analysis
        
-  Bool_t            fCleanOutputStdAOD;   // clean the written standard tracks and caloclusters in output AOD
-  TString        fDeltaAODFileName ;   // Delta AOD file name
-  TString        fFiredTriggerClassName  ;  // Name of trigger event type used to do the analysis
-
-  TString        fEMCALGeoName;       // Name of geometry to use for EMCAL.
-  TString        fPHOSGeoName;        // Name of geometry to use for PHOS.     
-  AliEMCALGeoUtils * fEMCALGeo ;      //! EMCAL geometry pointer
-  AliPHOSGeoUtils  * fPHOSGeo  ;      //! PHOS  geometry pointer  
-  Bool_t         fEMCALGeoMatrixSet;  // Check if the transformation matrix is set for EMCAL
-  Bool_t         fPHOSGeoMatrixSet ;  // Check if the transformation matrix is set for PHOS
-  Bool_t         fAnaLED;             // Analyze LED data only.
-  Bool_t         fRemoveBadChannels;  // Check the channel status provided and remove clusters with bad channels
-  TObjArray      fEMCALBadChannelMap; //! Array of histograms with map of bad channels, EMCAL
-  TObjArray      fPHOSBadChannelMap;  //! Array of histograms with map of bad channels, PHOS
-
-  ClassDef(AliCaloTrackReader,12)
+  AliCalorimeterUtils * fCaloUtils ;    //  Pointer to CalorimeterUtils
+
+  AliMixedEvent  * fMixedEvent  ;       //! mixed event object. This class is not the owner
+  Int_t            fNMixedEvent ;       // number of events in mixed event buffer
+  Double_t      ** fVertex      ;       //! vertex array 3 dim for each mixed event buffer
+  
+  Bool_t           fWriteOutputDeltaAOD;// Write the created delta AOD objects into file  
+       Bool_t           fOldAOD;             // Old AODs, before revision 4.20
+  
+  Int_t            fV0ADC[2]    ;       // Integrated V0 signal
+  Int_t            fV0Mul[2]    ;       // Integrated V0 Multiplicity
+
+  Bool_t           fCaloFilterPatch;    // CaloFilter patch
+  
+  ClassDef(AliCaloTrackReader,21)
 } ;