STEER/AliAODHandler.h

   1 #ifndef ALIAODHANDLER_H
   2 #define ALIAODHANDLER_H
   3 /* Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 //-------------------------------------------------------------------------
   9 //     Implementation of the Event Handler Interface for AOD
  10 //     Author: Andreas Morsch, CERN
  11 //-------------------------------------------------------------------------
  12
  13 #include "AliVEventHandler.h"
  14
  15 class AliAODEvent;
  16 class TFile;
  17 class TTree;
  18 class TObjArray;
  19 class AliMCEventHandler;
  20 class AliAODMCHeader;
  21 class AliAODExtension;
  22 class AliGenEventHeader;
  23
  24
  25
  26 class AliAODHandler : public AliVEventHandler {
  27
  28  public:
  29     AliAODHandler();
  30     AliAODHandler(const char* name, const char* title);
  31     virtual ~AliAODHandler();
  32     virtual void         SetOutputFileName(const char* fname);
  33     virtual const char*  GetOutputFileName();
  34     // Extra outputs as a string separated by commas
  35     virtual const char*  GetExtraOutputs() const;
  36     virtual Bool_t       Init(Option_t* option);
  37     virtual Bool_t       Init(TTree* /*tree*/, Option_t* /*option*/)  {return kTRUE;}
  38     virtual Bool_t       GetEntry() {return kTRUE;}
  39
  40     virtual Bool_t       BeginEvent(Long64_t /*entry*/) {fFillAOD=kFALSE; return kTRUE;}
  41     virtual Bool_t       Notify() { return AliVEventHandler::Notify(); };
  42     virtual Bool_t       Notify(const char * /* path */) {return kTRUE;}
  43     virtual Bool_t       FinishEvent();
  44     virtual Bool_t       Terminate();
  45     virtual Bool_t       TerminateIO();
  46     //
  47     virtual void         SetCreateNonStandardAOD()               {fIsStandard = kFALSE;}
  48     virtual void         SetFillAOD(Bool_t b)                    {fFillAOD = b;}
  49     virtual void         SetFillAODforRun(Bool_t b)              {fFillAODRun = b;}
  50     virtual void         SetNeedsHeaderReplication()             {fNeedsHeaderReplication             = kTRUE;}
  51     virtual void         SetNeedsTracksBranchReplication()       {fNeedsTracksBranchReplication       = kTRUE;}
  52     virtual void         SetNeedsVerticesBranchReplication()     {fNeedsVerticesBranchReplication     = kTRUE;}
  53     virtual void         SetNeedsV0sBranchReplication()          {fNeedsV0sBranchReplication          = kTRUE;}
  54     virtual void         SetNeedsCascadesBranchReplication()     {fNeedsCascadesBranchReplication     = kTRUE;}
  55     virtual void         SetNeedsTrackletsBranchReplication()    {fNeedsTrackletsBranchReplication    = kTRUE;}
  56     virtual void         SetNeedsPMDClustersBranchReplication()  {fNeedsPMDClustersBranchReplication  = kTRUE;}
  57     virtual void         SetNeedsJetsBranchReplication()         {fNeedsJetsBranchReplication         = kTRUE;}
  58     virtual void         SetNeedsFMDClustersBranchReplication()  {fNeedsFMDClustersBranchReplication  = kTRUE;}
  59     virtual void         SetNeedsCaloClustersBranchReplication() {fNeedsCaloClustersBranchReplication = kTRUE;}
  60     virtual void         SetNeedsMCParticlesBranchReplication()  {fNeedsMCParticlesBranchReplication  = kTRUE;}
  61     virtual void         SetNeedsDimuonsBranchReplication()      {fNeedsDimuonsBranchReplication      = kTRUE;}
  62     virtual void         SetAODIsReplicated() {fAODIsReplicated = kTRUE;}
  63     //
  64     AliAODEvent*         GetAOD()  {return fAODEvent;}
  65     virtual TTree*       GetTree() const {return fTreeA;}
  66     TObjArray*           GetExtensions() const {return fExtensions;}
  67     AliAODExtension*     GetExtension(const char *filename) const;
  68     TObjArray*           GetFilters() const {return fFilters;}
  69     AliAODExtension*     GetFilteredAOD(const char *filename) const;
  70     void                 CreateTree(Int_t flag);
  71     void                 FillTree();
  72     void                 AddAODtoTreeUserInfo();
  73     void                 AddBranch(const char* cname, void* addobj, const char *fname="");
  74     AliAODExtension*     AddExtension(const char *filename, const char *title="");
  75     AliAODExtension*     AddFilteredAOD(const char *filename, const char *filtername);
  76     Bool_t               IsStandard()                         const {return fIsStandard;}
  77     Bool_t               GetFillAOD()                         const {return fFillAOD;}
  78     Bool_t               NeedsHeaderReplication()             const {return  fNeedsHeaderReplication;}
  79     Bool_t               NeedsTracksBranchReplication()       const {return  fNeedsTracksBranchReplication;}
  80     Bool_t               NeedsVerticesBranchReplication()     const {return  fNeedsVerticesBranchReplication;}
  81     Bool_t               NeedsV0sBranchReplication()          const {return  fNeedsV0sBranchReplication;}
  82     Bool_t               NeedsCascadesBranchReplication()     const {return  fNeedsCascadesBranchReplication;}
  83     Bool_t               NeedsTrackletsBranchReplication()    const {return  fNeedsTrackletsBranchReplication;}
  84     Bool_t               NeedsPMDClustersBranchReplication()  const {return  fNeedsPMDClustersBranchReplication;}
  85     Bool_t               NeedsJetsBranchReplication()         const {return  fNeedsJetsBranchReplication;}
  86     Bool_t               NeedsFMDClustersBranchReplication()  const {return  fNeedsFMDClustersBranchReplication;}
  87     Bool_t               NeedsCaloClustersBranchReplication() const {return  fNeedsCaloClustersBranchReplication;}
  88     Bool_t               NeedsMCParticlesBranchReplication()  const {return  fNeedsMCParticlesBranchReplication;}
  89     Bool_t               NeedsDimuonsBranchReplication()      const {return  fNeedsDimuonsBranchReplication;}
  90     Bool_t               AODIsReplicated()                    const {return  fAODIsReplicated;}
  91     //
  92     void                 SetInputTree(TTree* /*tree*/) {;}
  93     void                 SetMCEventHandler(AliMCEventHandler* mcH) {fMCEventH = mcH;} // For internal use
  94     void StoreMCParticles(); // Store MC particles, only to be called from AliAnalyisTaskMCParticleFilter
  95
  96  private:
  97     void SetMCHeaderInfo(AliAODMCHeader *mcHeader,AliGenEventHeader *genHeader); // Utility function t catch different types of eventheaders
  98     AliAODHandler(const AliAODHandler&);             // Not implemented
  99     AliAODHandler& operator=(const AliAODHandler&);  // Not implemented
 100  private:
 101     Bool_t                   fIsStandard;                         // Flag for standard aod creation
 102     Bool_t                   fFillAOD;                            // Flag for filling of the AOD tree at the end (all or nothing evt by evt)
 103     Bool_t                   fFillAODRun;                         // Flag for filling of the AOD tree at the end (run)
 104     Bool_t                   fNeedsHeaderReplication;             // Flag for header replication
 105     Bool_t                   fNeedsTracksBranchReplication;       // Flag for tracks replication
 106     Bool_t                   fNeedsVerticesBranchReplication;     // Flag for vertices replication
 107     Bool_t                   fNeedsV0sBranchReplication;          // Flag for V0s replication
 108     Bool_t                   fNeedsCascadesBranchReplication;     // Flag for Cascade replication
 109     Bool_t                   fNeedsTrackletsBranchReplication;    // Flag for Tracklets replication
 110     Bool_t                   fNeedsPMDClustersBranchReplication;  // Flag for PMDClusters replication
 111     Bool_t                   fNeedsJetsBranchReplication;         // Flag for Jets replication
 112     Bool_t                   fNeedsFMDClustersBranchReplication;  // Flag for FMDClusters replication
 113     Bool_t                   fNeedsCaloClustersBranchReplication; // Flag for CaloClusters replication
 114     Bool_t                   fNeedsMCParticlesBranchReplication;  // Flag for MCParticles replication
 115     Bool_t                   fNeedsDimuonsBranchReplication;      // Flag for Dimuons replication
 116     Bool_t                   fAODIsReplicated;                    // Flag true if replication as been executed
 117     AliAODEvent             *fAODEvent;               //! Pointer to the AOD event
 118     AliMCEventHandler       *fMCEventH;               //! Pointer to mc event handler needed not to depend on the manager
 119     TTree                   *fTreeA;                  //! tree for AOD persistency
 120     TFile                   *fFileA;                  //! Output file
 121     TString                  fFileName;               //  Output file name
 122     TObjArray               *fExtensions;             //  List of extensions
 123     TObjArray               *fFilters;                //  List of filtered AOD's
 124     ClassDef(AliAODHandler, 6)
 125 };
 126
 127 //-------------------------------------------------------------------------
 128 //     Support class for AOD extensions. This is created by the user analysis
 129 //     that requires a separate file for some AOD branches. The name of the
 130 //     AliAODExtension object is the file name where the AOD branches will be
 131 //     stored.
 132 //     Author: Andrei Gheata, CERN
 133 //-------------------------------------------------------------------------
 134
 135 class AliAODExtension : public TNamed {
 136
 137 public:
 138
 139 enum EAliAODExtensionFlags {
 140    kFilteredAOD      = BIT(14)
 141 };
 142
 143     AliAODExtension() : TNamed(), fAODEvent(0), fTreeE(0), fFileE(0), fNtotal(0), fNpassed(0), fSelected(kFALSE) {}
 144     AliAODExtension(const char* name, const char* title, Bool_t isfilter=kFALSE);
 145     virtual ~AliAODExtension();
 146     void                 AddBranch(const char* cname, void* addobj);
 147     Bool_t               FinishEvent();
 148     Int_t                GetNtotal() const         {return fNtotal;}
 149     Int_t                GetNpassed() const        {return fNpassed;}
 150     const char*          GetOutputFileName() const {return TNamed::GetName();}
 151     AliAODEvent*         GetAOD() const            {return fAODEvent;}
 152     TTree*               GetTree() const           {return fTreeE;}
 153     Bool_t               Init(Option_t *option);
 154     Bool_t               IsFilteredAOD() const     {return TObject::TestBit(kFilteredAOD);}
 155     Bool_t               IsEventSelected() const   {return fSelected;}
 156     void                 SelectEvent(Bool_t flag=kTRUE)  {fSelected = flag;}
 157     void                 SetEvent(AliAODEvent *event);
 158     void                 SetOutputFileName(const char* fname) {TNamed::SetName(fname);}
 159     Bool_t               TerminateIO();
 160  private:
 161     AliAODExtension(const AliAODExtension&);             // Not implemented
 162     AliAODExtension& operator=(const AliAODExtension&);  // Not implemented
 163
 164  private:
 165     AliAODEvent             *fAODEvent;               //! Pointer to the AOD event
 166     TTree                   *fTreeE;                  //! tree for AOD persistency
 167     TFile                   *fFileE;                  //! Output file
 168     Int_t                    fNtotal;                 //! Number of processed events
 169     Int_t                    fNpassed;                //! Number of events that passed the filter
 170     Bool_t                   fSelected;               //! Select current event for filtered AOD's. Made false at event start.
 171     ClassDef(AliAODExtension, 1)                      // Support for extra AOD branches in a separate AOD file
 172 };
 173 #endif