Add class to access simulation parameters, AliEMCALSimParam, to be set in configurati...
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALDigitizer.h
1 #ifndef ALIEMCALDIGITIZER_H
2 #define ALIEMCALDIGITIZER_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 // Class that performs digitization of Summable digits from simulated data
10 //  
11 // In addition it performs mixing of summable digits from different events.
12 //
13 // For each event two branches are created in TreeD:
14 //   "EMCAL" - list of digits
15 //   "AliEMCALDigitizer" - AliEMCALDigitizer with all parameters used in digitization
16 //
17 // Note, that one cset title for new digits branch, and repeat digitization with
18 // another set of parameters.
19 //
20 //*-- Author: Sahal Yacoob (LBL)
21 // based on : AliPHOSDigit
22 // July     2003 Yves Schutz : NewIO 
23 // November 2003 Aleksei Pavlinov : Shish-Kebab geometry 
24 //_________________________________________________________________________ 
25
26
27 // --- ROOT system ---
28 class TClonesArray ; 
29 class TBrowser;
30
31 // --- Standard library ---
32
33 // --- AliRoot header files ---
34 #include "AliDigitizer.h"
35 #include "AliConfig.h"
36 class AliEMCALCalibData ;
37 class AliEMCALSDigitizer ;
38 class AliRunDigitizer ;
39
40 class AliEMCALDigitizer: public AliDigitizer {
41
42 public:
43   AliEMCALDigitizer() ;          // ctor
44   AliEMCALDigitizer(TString alirunFileNameFile, TString eventFolderName = AliConfig::GetDefaultEventFolderName()) ;  
45   AliEMCALDigitizer(const AliEMCALDigitizer & dtizer) ;
46   AliEMCALDigitizer(AliRunDigitizer * manager) ;
47   virtual ~AliEMCALDigitizer() ;       
48
49   void    Digitize(Int_t event);          // Make Digits from SDigits stored in fSDigits
50   void    Exec(Option_t *option);               // Supervising method
51
52   Float_t GetDigitThreshold() const { return fDigitThreshold;}
53   //Float_t GetPedestal()       const { return fPedestal; }
54   Float_t GetPinNoise()       const { return fPinNoise;}
55   //Float_t GetSlope()          const { return fSlope; }
56   Float_t GetTimeResolution() const { return fTimeResolution ; }
57   Float_t GetECAchannel()     const { return fADCchannelEC ; }
58   Float_t GetECApedestal()    const { return fADCpedestalEC ; }
59   void    SetEventRange(Int_t first=0, Int_t last=-1) {fFirstEvent=first; fLastEvent=last; }
60   void    SetDigitThreshold(Float_t EMCThreshold)  {fDigitThreshold = EMCThreshold;}
61   void    SetPinNoise(Float_t PinNoise )         {fPinNoise = PinNoise;}
62
63   //General
64   Int_t GetDigitsInRun()  const { return fDigitsInRun; } 
65   void  MixWith(TString alirunFileName, 
66                 TString eventFolderName = AliConfig::GetDefaultEventFolderName()) ; // Add another one file to mix
67   void  Print(Option_t* option="") const ;
68   void  Print1(Option_t * option); // *MENU*
69  
70   AliEMCALDigitizer & operator = (const AliEMCALDigitizer & /*rvalue*/)  {
71     // assignement operator requested by coding convention but not needed
72    Fatal("operator =", "not implemented") ;  
73    return *this ; 
74   }
75
76   virtual void Browse(TBrowser* b);
77
78 private:
79
80   Bool_t  Init();                   
81   void    InitParameters() ; 
82   void    PrintDigits(Option_t * option) ;
83   void    Unload() ; 
84   void    WriteDigits() ;         // Writes Digits the current event
85   Float_t TimeOfNoise(void) ;     // Calculate time signal generated by noise
86
87   //Calculate the time of crossing of the threshold by front edge
88   //Float_t FrontEdgeTime(TClonesArray * ticks) ;
89         
90   Int_t   DigitizeEnergy(Float_t energy, Int_t AbsId) ;
91
92 private:
93   
94   Bool_t  fDefaultInit;           //! Says if the task was created by defaut ctor (only parameters are initialized)
95   Int_t   fDigitsInRun ;          //! Total number of digits in one run
96   Bool_t  fInit ;                 //! To avoid overwriting existing files
97
98   Int_t   fInput ;                // Number of files to merge
99   TString * fInputFileNames ;     //[fInput] List of file names to merge 
100   TString * fEventNames ;         //[fInput] List of event names to merge
101
102   Float_t fDigitThreshold  ;      // Threshold for storing digits in EMC
103   Int_t   fMeanPhotonElectron ;   // number of photon electrons per GeV deposited energy 
104   //Float_t fPedestal ;           // Calibration parameters //Not used, remove?
105   //Float_t fSlope ;              // read from SDigitizer   //Not used, remove?
106   Float_t fPinNoise ;             // Electronics noise in EMC
107   Float_t fTimeResolution ;       // Time resolution of FEE electronics
108   //Float_t fTimeThreshold ;        // Threshold to start timing for given crystall //Not used, remove?
109   //Float_t fTimeSignalLength ;     // Length of the timing signal //Not used, remove?
110   Float_t fADCchannelEC ;         // width of one ADC channel in EC section (GeV)
111   Float_t fADCpedestalEC ;        // pedestal for one ADC channel
112   Int_t   fNADCEC ;               // number of channels in EC section ADC
113
114   TString fEventFolderName;         // skowron: name of EFN to read data from in stand alone mode
115   Int_t   fFirstEvent;        // first event to process
116   Int_t   fLastEvent;         // last  event to process
117
118   AliEMCALCalibData * fCalibData; //Calibration data pointer
119
120   ClassDef(AliEMCALDigitizer,8)  // description 
121 };
122
123
124 #endif // AliEMCALDIGITIZER_H