add protection against truncated events + coverity - Rachid

[u/mrichter/AliRoot.git] / EMCAL / AliEMCALSDigitizer.h

diff --git a/EMCAL/AliEMCALSDigitizer.h b/EMCAL/AliEMCALSDigitizer.h
index 4598365f439a70831fdc9d326855aee47274d8e1..9c88d82078e26cbc52f176f733d8b116202f16ae 100644 (file)
--- a/EMCAL/AliEMCALSDigitizer.h
+++ b/EMCAL/AliEMCALSDigitizer.h
@@ -1,72 +1,77 @@
-#ifndef ALIEMCALSDigitizer_H
-#define ALIEMCALSDigitizer_H
+#ifndef ALIEMCALSDIGITIZER_H
+#define ALIEMCALSDIGITIZER_H
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
 /* $Id$ */
 //_________________________________________________________________________
-//  Task Class for making SDigits in EMCAL      
-//                  
+// This is a TTask that makes SDigits out of Hits
+// A Summable Digits is the sum of all hits originating 
+// from one in one tower of the EMCAL 
+// A threshold for assignment of the primary to SDigit is applied 
+//
+// SDigits need to hold the energy sum of the hits, but AliEMCALDigit
+// can (should) only store amplitude.  Therefore, the SDigit energy is
+// "digitized" before being stored and must be "calibrated" back to an
+// energy before SDigits are summed to form true Digits
+//
+//
 //*-- Author: Sahal Yacoob (LBL)
 // based on : AliPHOSSDigitizer
 //_________________________________________________________________________
-//
-// Modif: 
-//  August 2002 Yves Schutz: clone PHOS as closely as possible and intoduction
-//                           of new  IO (à la PHOS)
  
 // --- ROOT system ---
 #include "TTask.h"
-#include "TString.h"
+class TBrowser;
+
 // --- Standard library ---
 
 // --- AliRoot header files ---
+#include "AliConfig.h"
 
 class AliEMCALSDigitizer: public TTask {
 
 public:
   AliEMCALSDigitizer() ;          // ctor
-  AliEMCALSDigitizer(const char* headerFile,const char* hdigitsTitle = "Default", const Bool_t toSplit = kFALSE) ; 
-  virtual ~AliEMCALSDigitizer() ; // dtor
+  AliEMCALSDigitizer(const char * alirunFileName, const char * eventFolderName = AliConfig::GetDefaultEventFolderName()) ; 
+  AliEMCALSDigitizer(const AliEMCALSDigitizer & sd) ;
+  virtual ~AliEMCALSDigitizer(); // dtor
+
+  Float_t       Digitize(Float_t energy)const; //convert energy in GeV to int amplitude
+  Float_t       Calibrate(Float_t amp)const;  //opposite of Digitize()
 
-  Float_t  Calibrate(Int_t amp)const {return (amp - fA)/fB ; }
-  Int_t    Digitize(Float_t Energy)const { return (Int_t ) ( fA + Energy*fB); }
   virtual void  Exec(Option_t *option); 
-  const char *  GetSDigitsBranch()const{return GetName();}  
-  const Int_t   GetSDigitsInRun() const {return fSDigitsInRun ;}  
-  const Float_t GetPedestalParameter()const {return fA;}
-  const Float_t GetCalibrationParameter()const{return fB;}
-  virtual void  Print(Option_t* option) const ;
-  void SetSDigitsBranch(const char * title ) ;
-  void SetPedestalParameter(Float_t A){fA = A ;}
-  void SetSlopeParameter(Float_t B){fB = B ;}
-  void UseHitsFrom(const char * filename) ;      
+  Int_t         GetSDigitsInRun() const {return fSDigitsInRun ;}  
+  virtual void  Print(Option_t *option="") const;
+  void          Print1(Option_t *option="all");  // *MENU*
+  void          SetEventFolderName(TString name) { fEventFolderName = name ; }
+  void          SetEventRange(Int_t first=0, Int_t last=-1) {fFirstEvent=first; fLastEvent=last; }
+
   Bool_t operator == (const AliEMCALSDigitizer & sd) const ;
-  const Int_t Segment2TowerID(Int_t SegmentID){
-    return Layer2TowerID(SegmentID,kFALSE);
-}
+  const AliEMCALSDigitizer & operator = (const AliEMCALSDigitizer & /*sd*/) {return *this ;}
+
+  virtual void Browse(TBrowser* b);
 
 private:
   void     Init() ;
   void     InitParameters() ; 
   void     PrintSDigits(Option_t * option) ;
-  const Int_t Layer2TowerID(Int_t,Bool_t) ;
+  void     Unload() const ;
 
 private:
   Float_t fA ;                     // Pedestal parameter
   Float_t fB ;                     // Slope Digitizition parameters
-  Float_t fPhotonElectronFactor ;  // number of photon electrons per GeV
-                                   // should be calculated independently for each layer as : 
-                                   // LightYield*LightCollectionEfficiency*LightAttenuation*APDPhotoElectronEfficiency*APDGain
-  Float_t fTowerPrimThreshold ;    // To store primary in Tower if Elos > threshold
-  Float_t fPreShowerPrimThreshold ;// To store primary if Pre Shower Elos > threshold
-   Bool_t  fDefaultInit;           //! Says if the task was created by defaut ctor (only parameters are initialized)
+  Float_t fECPrimThreshold ;       // To store primary if EC Shower Elos > threshold
+  Bool_t  fDefaultInit;            //! Says if the task was created by defaut ctor (only parameters are initialized)
+  TString fEventFolderName;        // event folder name
+  Bool_t  fInit ;                  //! tells if initialisation went OK, will revent exec if not
   Int_t   fSDigitsInRun ;          //! Total number of sdigits in one run
-  TFile * fSplitFile ;             //! file in which SDigits will eventually be stored
-  Bool_t  fToSplit ;               //! Says that sigits should be written into splip file
-
-  ClassDef(AliEMCALSDigitizer,3)  // description 
-
+  Int_t   fFirstEvent;             // first event to process
+  Int_t   fLastEvent;              // last  event to process
+  Float_t fSampling;               // See AliEMCALGeometry
+  TClonesArray* fHits;             //-> Temporal array with hits
+       
+  ClassDef(AliEMCALSDigitizer,7)  // description 
 };
 
-#endif // AliEMCALSDigitizer_H
+#endif // AliEMCALSDIGITIZER_H