]>
Commit | Line | Data |
---|---|---|
2012850d | 1 | #ifndef ALIEMCAL_H |
2 | #define ALIEMCAL_H | |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
6 | /* $Id$ */ | |
2576b1b8 | 7 | /* History of cvs commits: |
8 | * | |
9 | * $Log$ | |
23ef18ac | 10 | * Revision 1.42 2007/02/24 20:42:35 pavlinov |
11 | * fixed error of Geant3 parameters initialisation | |
12 | * | |
91e5f344 | 13 | * Revision 1.41 2007/02/05 10:43:25 hristov |
14 | * Changes for correct initialization of Geant4 (Mihaela) | |
15 | * | |
7235aed2 | 16 | * Revision 1.40 2006/12/05 17:19:26 gustavo |
17 | * Updated AliEMCAL::Digits2Raw, reads first provisional RCU mapping files to make Raw data with new AliCaloAltroMapping and AliCaloRawStream | |
18 | * | |
2576b1b8 | 19 | * |
20 | */ | |
2012850d | 21 | //_________________________________________________________________________ |
22 | // Base Class for EMCAL | |
14ce0a6e | 23 | // holds all geant information of |
24 | // materials, etc. | |
2012850d | 25 | // |
23ef18ac | 26 | //*-- Author: Yves Schutz (SUBATECH) |
2012850d | 27 | |
2012850d | 28 | // --- ROOT system --- |
05a92d59 | 29 | |
2012850d | 30 | class TString ; |
05a92d59 | 31 | class TTask ; |
32 | class TFolder ; | |
f51151a0 | 33 | class TRandom ; |
2d5d9e60 | 34 | class TGraph; |
35 | class TF1; | |
23ef18ac | 36 | class AliEMCALGeometry; |
2012850d | 37 | |
38 | // --- AliRoot header files --- | |
2d5d9e60 | 39 | class AliRawReader; |
2012850d | 40 | #include "AliDetector.h" |
8367ce9a | 41 | #include "AliEMCALGeometry.h" |
0631a867 | 42 | #include "AliEMCALTrigger.h" |
b4215a15 | 43 | |
2012850d | 44 | class AliEMCAL : public AliDetector { |
45 | ||
46 | public: | |
b4215a15 | 47 | |
b13bbe81 | 48 | AliEMCAL(); |
2012850d | 49 | AliEMCAL(const char* name, const char* title=""); |
0a4cb131 | 50 | |
2012850d | 51 | virtual ~AliEMCAL() ; |
5dee926e | 52 | virtual void AddHit(Int_t, Int_t*, Float_t *) { |
d64c959b | 53 | Fatal("AddHit(Int_t, Int_t*, Float_t *", "not to be used: use AddHit( Int_t shunt, Int_t primary, Int_t track,Int_t id, Float_t *hits )") ; |
b13bbe81 | 54 | } |
8367ce9a | 55 | virtual AliDigitizer* CreateDigitizer(AliRunDigitizer* manager) const; |
f51151a0 | 56 | virtual void CreateMaterials() ; |
57 | virtual void Digits2Raw(); | |
b4215a15 | 58 | |
4800667c | 59 | using AliDetector::Raw2Digits; |
2d5d9e60 | 60 | virtual void Raw2Digits(AliRawReader *reader); |
b4215a15 | 61 | |
fdebddeb | 62 | virtual void FinishRun() {} |
8367ce9a | 63 | virtual AliEMCALGeometry * GetGeometry() const |
b4215a15 | 64 | {return AliEMCALGeometry::GetInstance(GetTitle(),"") ; } |
8367ce9a | 65 | virtual void Hits2SDigits(); |
91e5f344 | 66 | // virtual void Init(); |
fdebddeb | 67 | virtual Int_t IsVersion(void) const = 0 ; |
b4215a15 | 68 | |
69 | virtual AliTriggerDetector* CreateTriggerDetector() const | |
70 | { return new AliEMCALTrigger(); } | |
0a4cb131 | 71 | |
56088960 | 72 | // Raw Read Out |
73 | Double_t GetRawFormatCapa() const { return fgCapa ; } | |
74 | Double_t GetRawFormatHighCharge() const { return fHighCharge ; } | |
75 | Double_t GetRawFormatHighGain() const { return fHighGain ; } | |
76 | Double_t GetRawFormatHighLowGainFactor() const { return fHighLowGainFactor ; } | |
77 | Double_t GetRawFormatLowCharge() const { return ( fHighCharge * fHighLowGainFactor ) ; } | |
78 | Double_t GetRawFormatLowGain() const { return ( fHighGain / fHighLowGainFactor ) ; } | |
79 | Int_t GetRawFormatLowGainOffset() const { return fLowGainOffset ; } | |
80 | Int_t GetRawFormatOrder() const { return fgOrder ; } | |
14ce0a6e | 81 | Int_t GetRawFormatTimeBins() const { return fgkTimeBins ; } |
56088960 | 82 | Double_t GetRawFormatTimeMax() const { return fgTimeMax ; } |
83 | Double_t GetRawFormatTimePeak() const { return fgTimePeak ; } | |
2d5d9e60 | 84 | Double_t GetRawFormatTimeTrigger() const { return fgTimeTrigger ; } |
2d5d9e60 | 85 | Int_t GetRawFormatThreshold() const { return fgThreshold ; } |
2576b1b8 | 86 | Int_t GetRawFormatDDLPerSuperModule() const { return fgDDLPerSuperModule ; } |
56088960 | 87 | static Double_t RawResponseFunctionMax(Double_t charge, Double_t gain) ; |
03ecfe88 | 88 | Bool_t RawSampledResponse(Double_t dtime, Double_t damp, Int_t * adcH, Int_t * adcL) const ; |
56088960 | 89 | // |
8367ce9a | 90 | virtual AliLoader* MakeLoader(const char* topfoldername); |
fdebddeb | 91 | virtual const TString Version() const {return TString(" ") ; } |
1963b290 | 92 | |
2012850d | 93 | protected: |
56088960 | 94 | |
fac5662b | 95 | static Double_t RawResponseFunction(Double_t *x, Double_t *par) ; |
2d5d9e60 | 96 | void FitRaw(Bool_t lowGainFlag, TGraph * gLowGain, TGraph * gHighGain, TF1* signalF, Double_t & energy, Double_t & time) ; |
f51151a0 | 97 | |
7235aed2 | 98 | void InitConstants(); //initializes some params |
91e5f344 | 99 | void DefineMediumParameters(); // define tracking medium parameters |
0a4cb131 | 100 | |
f51151a0 | 101 | Int_t fBirkC0; // constants for Birk's Law implementation |
102 | Double_t fBirkC1; // constants for Birk's Law implementation | |
103 | Double_t fBirkC2; // constants for Birk's Law implementation | |
1963b290 | 104 | |
56088960 | 105 | static Double_t fgCapa ; // capacitor of the preamplifier for the raw RO signal |
106 | Double_t fHighCharge ; // high charge (to convert energy to charge) for the raw RO signal | |
107 | Double_t fHighGain ; // high gain for the raw RO signal | |
108 | Double_t fHighLowGainFactor ; // high to low gain factor for the raw RO signal | |
109 | Int_t fLowGainOffset ; // to separate high from low gain in the DDL | |
110 | static Int_t fgOrder ; // order of the gamma function for the RO signal | |
14ce0a6e | 111 | static const Int_t fgkTimeBins = 256 ; // number of sampling bins of the raw RO signal |
56088960 | 112 | static Double_t fgTimeMax ; // maximum sampled time of the raw RO signal |
113 | static Double_t fgTimePeak ; // peaking time of the raw RO signal | |
114 | static Double_t fgTimeTrigger ; // time of the trigger for the RO signal | |
14ce0a6e | 115 | static Int_t fgThreshold; // threshold |
23ef18ac | 116 | static Int_t fgDDLPerSuperModule; // number of DDL per SuperModule |
2576b1b8 | 117 | |
23ef18ac | 118 | AliEMCALGeometry* fGeometry; //! |
119 | ||
2576b1b8 | 120 | private: |
e939a978 | 121 | AliEMCAL(const AliEMCAL& emcal); |
122 | AliEMCAL & operator = (const AliEMCAL & /*rvalue*/); | |
123 | ||
56088960 | 124 | ClassDef(AliEMCAL,9) // Electromagnetic calorimeter (base class) |
125 | ||
126 | } ; | |
2012850d | 127 | |
128 | #endif // ALIEMCAL_H |