PHOS/AliPHOSLoader.h

   1 #ifndef ALIPHOSLOADER_H
   2 #define ALIPHOSLOADER_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 //_________________________________________________________________________
   9 //  A singleton that returns various objects
  10 //  Should be used on the analysis stage to avoid confusing between different
  11 //  branches of reconstruction tree: e.g. reading RecPoints and TS made from
  12 //  another set of RecPoints.
  13 //
  14 //  The objects are retrived from folders.
  15 //*-- Author: Yves Schutz (SUBATECH) & Dmitri Peressounko (RRC KI & SUBATECH)
  16 //
  17
  18
  19 // --- ROOT system ---
  20 #include "TClonesArray.h"
  21 class TString ;
  22 class TParticle ;
  23 class TTask ;
  24
  25 // --- Standard library ---
  26
  27 // --- AliRoot header files ---
  28
  29 #include "AliRun.h"
  30 #include "AliLoader.h"
  31 #include "AliRunLoader.h"
  32 #include "AliPHOSClusterizer.h"
  33 #include "AliPHOSTrackSegmentMaker.h"
  34 #include "AliPHOSPID.h"
  35 class AliPHOS ;
  36 class AliPHOSGeometry ;
  37 class AliPHOSHit ;
  38 class AliPHOSDigit ;
  39 class AliPHOSEmcRecPoint ;
  40 class AliPHOSCpvRecPoint ;
  41 class AliPHOSTrackSegment ;
  42 class AliPHOSRecParticle ;
  43 class AliPHOSCalibrationDB ;
  44 class AliPHOSSDigitizer ;
  45 class AliPHOSDigitizer ;
  46
  47
  48
  49 //
  50
  51 class AliPHOSLoader : public AliLoader {
  52
  53  public:
  54
  55   AliPHOSLoader();
  56   AliPHOSLoader(const AliPHOSLoader & obj) : AliLoader(obj) {}
  57   AliPHOSLoader(const Char_t *detname,const Char_t *eventfoldername);
  58
  59   virtual ~AliPHOSLoader() ;
  60
  61   // assignement operator requested by coding convention, but not needed
  62   AliPHOSLoader & operator = (const AliPHOSLoader & ) {return *this;}
  63
  64   Int_t   GetEvent();//extends the method on PHOS RecPart posting
  65   Int_t   SetEvent();//extends the method on PHOS RecPart posting
  66
  67   Bool_t  BranchExists(const TString& recName);
  68   Int_t   LoadHits(Option_t* opt=""); //reads  from disk and sends them to folder; array as well as tree
  69   Int_t   LoadSDigits(Option_t* opt="");
  70   Int_t   LoadDigits(Option_t* opt=""); //reads Digits from disk and sends them to folder; array as well as tree
  71   Int_t   LoadRecPoints(Option_t* opt=""); //reads RecPoints from disk and sends them to folder; array as well as tree
  72   Int_t   LoadTracks(Option_t* opt="");  //reads Tracks from disk and sends them to folder; array as well as tree
  73   Int_t   LoadRecParticles(Option_t* opt="");
  74
  75   Int_t   PostHits();  //Posts the
  76   Int_t   PostSDigits();
  77   Int_t   PostDigits();
  78   Int_t   PostRecPoints();
  79   Int_t   PostTracks();
  80   Int_t   PostRecParticles();
  81
  82   void    CleanFolders();//cleans all the stuff loaded by this detector + calls AliLoader::Clean
  83
  84   void    CleanHits();
  85   void    CleanSDigits();
  86   void    CleanDigits();
  87   void    CleanRecPoints();
  88   void    CleanTracks();
  89   void    CleanRecParticles();
  90
  91 //up to now it is only here -> no definition about global/incremental tracking/PID
  92
  93 //   Int_t   WriteRecParticles(Option_t* opt="");//writes the reconstructed particles
  94 //   Int_t   WritePID(Option_t* opt="");//writes the task for PID to file
  95 //   Bool_t  PostPID  (AliPHOSPID * pid) const {return kTRUE;}
  96 //  Bool_t  PostQA   (void) const ; //it was empty anyway
  97
  98 /*******************************************************************/
  99 /*******************************************************************/
 100 /*******************************************************************/
 101
 102   TObject** HitsRef(){return GetDetectorDataRef(Hits());}
 103   TObject** SDigitsRef(){return GetDetectorDataRef(SDigits());}
 104   TObject** DigitsRef(){return GetDetectorDataRef(Digits());}
 105   TObject** EmcRecPointsRef(){return GetDetectorDataRef(EmcRecPoints());}
 106   TObject** CpvRecPointsRef(){return GetDetectorDataRef(CpvRecPoints());}
 107   TObject** TracksRef(){return GetDetectorDataRef(TrackSegments());}
 108   TObject** RecParticlesRef(){return GetDetectorDataRef(RecParticles());}
 109   TObject** AlarmsRef(){return GetDetectorDataRef(Alarms());}
 110   void   Track(Int_t itrack) ;
 111
 112   static AliPHOSGeometry* GetPHOSGeometry();
 113   static AliPHOSLoader* GetPHOSLoader(const  char* eventfoldername);
 114
 115   //Method to be used when digitizing under AliRunDigitizer, who opens all files etc.
 116   Int_t  EventNumber()       { return (Int_t) GetRunLoader()->GetEventNumber();}
 117   Int_t  MaxEvent()          { return (Int_t) GetRunLoader()->TreeE()->GetEntries();}
 118
 119   const AliPHOS *         PHOS();
 120   const AliPHOSGeometry  *PHOSGeometry() ;
 121   // Alarms
 122   // TFolder * Alarms() const { return (TFolder*)(ReturnO("Alarms", 0)); }
 123   TObjArray *  Alarms();
 124
 125   /*********************************************/
 126   /************    TClonesArrays     ***********/
 127   /*********************************************/
 128   /****   H i t s  ****/
 129   TClonesArray*  Hits(void);
 130   const AliPHOSHit*    Hit(Int_t index);
 131   void MakeHitsArray();
 132   /****   S D i g i t s  ****/
 133   TClonesArray*  SDigits();
 134   const AliPHOSDigit*  SDigit(Int_t index);
 135   void MakeSDigitsArray();
 136   /****  D i g i t s  ****/
 137   TClonesArray*   Digits();
 138   const AliPHOSDigit *  Digit(Int_t index);
 139   void MakeDigitsArray();
 140   /****  R e c P o i n t s  ****/
 141   TObjArray * EmcRecPoints();
 142   TObjArray * CpvRecPoints();
 143   const AliPHOSEmcRecPoint * EmcRecPoint(Int_t index) ;
 144   const AliPHOSCpvRecPoint * CpvRecPoint(Int_t index) ;
 145   void MakeRecPointsArray();
 146   /****   T r a c k S e g m e n t s ****/
 147   TClonesArray * TrackSegments();
 148   const AliPHOSTrackSegment * TrackSegment(Int_t index);
 149   void MakeTrackSegmentsArray();
 150   /****  R e c P a r t ic l e s   ****/
 151   TClonesArray * RecParticles() ;
 152   const AliPHOSRecParticle * RecParticle(Int_t index);
 153   void MakeRecParticlesArray();
 154
 155   /*********************************************/
 156   /************    T A S K S      **************/
 157   /*********************************************/
 158   //
 159   //  AliPHOSSDigitizer*  PHOSSDigitizer(TString name = AliConfig::fgkDefaultEventFolderName);
 160   //AliPHOSDigitizer*   PHOSDigitizer()  { return  dynamic_cast<AliPHOSDigitizer*>(Digitizer()) ;}
 161
 162   AliPHOSClusterizer* Clusterizer ()  {return dynamic_cast<AliPHOSClusterizer*>(Reconstructioner()) ;}
 163   Int_t PostClusterizer(TTask* clust){return PostReconstructioner(clust);}
 164   Int_t LoadClusterizer(Option_t * opt="") {return LoadReconstructioner(opt);}
 165   Int_t WriteClusterizer(Option_t * opt="") {return WriteReconstructioner(opt);}
 166
 167   AliPHOSPID * PID (){return dynamic_cast<AliPHOSPID*>(PIDTask()) ;}
 168   Int_t PostPID(TTask* pid){return PostPIDTask(pid);}
 169   Int_t LoadPID(Option_t * opt="") {return LoadPIDTask(opt);}
 170   Int_t WritePID(Option_t * opt="") {return WritePIDTask(opt);}
 171
 172
 173   AliPHOSTrackSegmentMaker * TrackSegmentMaker ()  { return dynamic_cast<AliPHOSTrackSegmentMaker *>(Tracker()) ;}
 174   Int_t PostTrackSegmentMaker(TTask* segmaker){return PostTracker(segmaker);}
 175   Int_t LoadTrackSegmentMaker(Option_t * opt="") {return LoadTracker(opt);}
 176   Int_t WriteTrackSegmentMaker(Option_t * opt="") {return WriteTracker(opt);}
 177
 178
 179   void   SetDebug(Int_t level) {fDebug = level;} // Set debug level
 180   void   SetBranchTitle(const TString& btitle);
 181
 182   AliPHOSCalibrationDB * CalibrationDB(){return  fcdb; }
 183   void ReadCalibrationDB(const char * name, const char * filename);
 184
 185  public:
 186
 187   static const TString fgkHitsName;//Name for TClonesArray with hits from one event
 188   static const TString fgkSDigitsName;//Name for TClonesArray
 189   static const TString fgkDigitsName;//Name for TClonesArray
 190   static const TString fgkEmcRecPointsName;//Name for TClonesArray
 191   static const TString fgkCpvRecPointsName;//Name for TClonesArray
 192   static const TString fgkTracksName;//Name for TClonesArray
 193   static const TString fgkRecParticlesName;//Name for TClonesArray
 194
 195   static const TString fgkEmcRecPointsBranchName;//Name for branch
 196   static const TString fgkCpvRecPointsBranchName;//Name for branch
 197   static const TString fgkTrackSegmentsBranchName;//Name for branch
 198   static const TString fgkRecParticlesBranchName;//Name for branch
 199
 200 protected:
 201   TString fBranchTitle;            //Title of the branch
 202   AliPHOSCalibrationDB * fcdb ;       //!
 203
 204 private:
 205
 206   Int_t ReadHits();
 207   Int_t ReadDigits();
 208   Int_t ReadSDigits();
 209   Int_t ReadRecPoints();
 210   Int_t ReadTracks();
 211   Int_t ReadRecParticles();
 212
 213   void  ReadTreeQA() ;
 214   Int_t  fDebug ;             // Debug level
 215
 216   ClassDef(AliPHOSLoader,3)  // Algorithm class that provides methods to retrieve objects from a list knowing the index
 217
 218 };
 219
 220 /******************************************************************************/
 221 /****************    I N L I N E S     ****************************************/
 222 /******************************************************************************/
 223
 224 inline TClonesArray* AliPHOSLoader::Hits()
 225 {
 226  return (TClonesArray*)GetDetectorData(fgkHitsName);
 227 }
 228 /******************************************************************************/
 229
 230 inline const AliPHOSHit* AliPHOSLoader::Hit(Int_t index)
 231 {
 232   const TClonesArray* tcarr = Hits();
 233   if (tcarr)
 234     return (const AliPHOSHit*) tcarr->At(index);
 235   return 0x0;
 236 }
 237 /******************************************************************************/
 238
 239 inline TClonesArray* AliPHOSLoader::SDigits()
 240 {
 241    return dynamic_cast<TClonesArray*>(GetDetectorData(fgkSDigitsName));
 242 }
 243 /******************************************************************************/
 244
 245 inline const AliPHOSDigit*  AliPHOSLoader::SDigit(Int_t index)
 246 {
 247   const TClonesArray* tcarr = SDigits();
 248   if (tcarr)
 249     return (const AliPHOSDigit*) tcarr->At(index);
 250   return 0x0;
 251 }
 252 /******************************************************************************/
 253
 254 inline TClonesArray* AliPHOSLoader::Digits()
 255 {
 256  return dynamic_cast<TClonesArray*>(GetDetectorData(fgkDigitsName));
 257 }
 258 /******************************************************************************/
 259
 260 inline const AliPHOSDigit*  AliPHOSLoader::Digit(Int_t index)
 261 {
 262   const TClonesArray* tcarr = Digits();
 263   if (tcarr)
 264     return (const AliPHOSDigit*) tcarr->At(index);
 265   return 0x0;
 266 }
 267 /******************************************************************************/
 268
 269 inline TObjArray * AliPHOSLoader::EmcRecPoints()
 270 {
 271  return dynamic_cast<TObjArray*>(GetDetectorData(fgkEmcRecPointsName));
 272 }
 273 /******************************************************************************/
 274
 275 inline const AliPHOSEmcRecPoint * AliPHOSLoader::EmcRecPoint(Int_t index)
 276 {
 277   TObjArray* tcarr = EmcRecPoints();
 278   if (tcarr)
 279     return (const AliPHOSEmcRecPoint*) tcarr->At(index);
 280   return 0x0;
 281 }
 282 /******************************************************************************/
 283
 284 inline TObjArray * AliPHOSLoader::CpvRecPoints()
 285 {
 286  return dynamic_cast<TObjArray*>(GetDetectorData(fgkCpvRecPointsName));
 287 }
 288 /******************************************************************************/
 289
 290 inline const AliPHOSCpvRecPoint * AliPHOSLoader::CpvRecPoint(Int_t index)
 291 {
 292   TObjArray* tcarr = CpvRecPoints();
 293   if (tcarr)
 294     return (const AliPHOSCpvRecPoint*) tcarr->At(index);
 295   return 0x0;
 296 }
 297 /******************************************************************************/
 298
 299 inline TClonesArray * AliPHOSLoader::TrackSegments()
 300 {
 301  return dynamic_cast<TClonesArray*>(GetDetectorData(fgkTracksName));
 302 }
 303 /******************************************************************************/
 304
 305 inline const AliPHOSTrackSegment * AliPHOSLoader::TrackSegment(Int_t index)
 306 {
 307   const TClonesArray* tcarr = TrackSegments();
 308   if (tcarr)
 309     return (const AliPHOSTrackSegment*) tcarr->At(index);
 310   return 0x0;
 311 }
 312 /******************************************************************************/
 313
 314 inline TClonesArray * AliPHOSLoader::RecParticles()
 315 {
 316  return dynamic_cast<TClonesArray*>(GetDetectorData(fgkRecParticlesName));
 317 }
 318 /******************************************************************************/
 319
 320 inline const AliPHOSRecParticle* AliPHOSLoader::RecParticle(Int_t index)
 321 {
 322   TClonesArray* tcarr = RecParticles();
 323   if (tcarr)
 324     return (const AliPHOSRecParticle*) tcarr->At(index);
 325   return 0x0;
 326 }
 327 /******************************************************************************/
 328 inline TObjArray *  AliPHOSLoader::Alarms()
 329 { return (TObjArray*)(GetQAFolder()->FindObject(fDetectorName));}
 330
 331 #endif // AliPHOSLOADER_H