]>
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$ */ | |
7 | ||
8 | //_________________________________________________________________________ | |
9 | // Base Class for EMCAL | |
10 | // | |
11 | //*-- Author: Yves Schutz (SUBATECH) | |
12 | ||
2012850d | 13 | // --- ROOT system --- |
05a92d59 | 14 | |
2012850d | 15 | class TString ; |
05a92d59 | 16 | class TTask ; |
17 | class TFolder ; | |
f51151a0 | 18 | class TRandom ; |
2d5d9e60 | 19 | class TGraph; |
20 | class TF1; | |
2012850d | 21 | |
22 | // --- AliRoot header files --- | |
2d5d9e60 | 23 | class AliRawReader; |
2012850d | 24 | #include "AliDetector.h" |
8367ce9a | 25 | #include "AliEMCALGeometry.h" |
2012850d | 26 | class AliEMCAL : public AliDetector { |
27 | ||
28 | public: | |
29 | ||
b13bbe81 | 30 | AliEMCAL(); |
2012850d | 31 | AliEMCAL(const char* name, const char* title=""); |
39200c71 | 32 | AliEMCAL(const AliEMCAL& emcal) : AliDetector(emcal) { |
2012850d | 33 | // cpy ctor: no implementation yet |
34 | // requested by the Coding Convention | |
9859bfc0 | 35 | Fatal("cpy ctor", "not implemented") ; |
2012850d | 36 | } |
37 | virtual ~AliEMCAL() ; | |
5dee926e | 38 | virtual void AddHit(Int_t, Int_t*, Float_t *) { |
d64c959b | 39 | 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 | 40 | } |
5dee926e | 41 | virtual void Copy(TObject & emcal) const |
42 | { Copy(dynamic_cast<AliEMCAL&>(emcal)); } | |
43 | virtual void Copy(AliEMCAL & emcal) const; | |
8367ce9a | 44 | virtual AliDigitizer* CreateDigitizer(AliRunDigitizer* manager) const; |
f51151a0 | 45 | virtual void CreateMaterials() ; |
905263da | 46 | // virtual void |
f51151a0 | 47 | virtual void Digits2Raw(); |
4800667c | 48 | |
49 | using AliDetector::Raw2Digits; | |
2d5d9e60 | 50 | virtual void Raw2Digits(AliRawReader *reader); |
4800667c | 51 | |
fdebddeb | 52 | virtual void FinishRun() {} |
8367ce9a | 53 | virtual AliEMCALGeometry * GetGeometry() const |
54 | {return AliEMCALGeometry::GetInstance(GetTitle(),"") ; } | |
55 | virtual void Hits2SDigits(); | |
fdebddeb | 56 | virtual Int_t IsVersion(void) const = 0 ; |
56088960 | 57 | // Raw Read Out |
58 | Double_t GetRawFormatCapa() const { return fgCapa ; } | |
59 | Double_t GetRawFormatHighCharge() const { return fHighCharge ; } | |
60 | Double_t GetRawFormatHighGain() const { return fHighGain ; } | |
61 | Double_t GetRawFormatHighLowGainFactor() const { return fHighLowGainFactor ; } | |
62 | Double_t GetRawFormatLowCharge() const { return ( fHighCharge * fHighLowGainFactor ) ; } | |
63 | Double_t GetRawFormatLowGain() const { return ( fHighGain / fHighLowGainFactor ) ; } | |
64 | Int_t GetRawFormatLowGainOffset() const { return fLowGainOffset ; } | |
65 | Int_t GetRawFormatOrder() const { return fgOrder ; } | |
046ae904 | 66 | Int_t GetRawFormatTimeBins() const { return fkTimeBins ; } |
56088960 | 67 | Double_t GetRawFormatTimeMax() const { return fgTimeMax ; } |
68 | Double_t GetRawFormatTimePeak() const { return fgTimePeak ; } | |
2d5d9e60 | 69 | Double_t GetRawFormatTimeTrigger() const { return fgTimeTrigger ; } |
70 | Int_t GetRawFormatDDLOffset() const { return fgDDLOffset ; } | |
71 | Int_t GetRawFormatThreshold() const { return fgThreshold ; } | |
72 | Int_t GetRawFormatChannelsPerDDL() const { return fgChannelsPerDDL ; } | |
56088960 | 73 | static Double_t RawResponseFunctionMax(Double_t charge, Double_t gain) ; |
74 | // | |
8367ce9a | 75 | virtual AliLoader* MakeLoader(const char* topfoldername); |
fdebddeb | 76 | virtual const TString Version() const {return TString(" ") ; } |
9c0a4862 | 77 | AliEMCAL & operator = (const AliEMCAL & /*rvalue*/) { |
d64c959b | 78 | Fatal("operator =", "not implemented") ; return *this ; } |
1963b290 | 79 | |
2012850d | 80 | protected: |
56088960 | 81 | |
fac5662b | 82 | static Double_t RawResponseFunction(Double_t *x, Double_t *par) ; |
56088960 | 83 | Bool_t RawSampledResponse(const Double_t dtime, const Double_t damp, Int_t * adcH, Int_t * adcL) const ; |
2d5d9e60 | 84 | void FitRaw(Bool_t lowGainFlag, TGraph * gLowGain, TGraph * gHighGain, TF1* signalF, Double_t & energy, Double_t & time) ; |
f51151a0 | 85 | |
f51151a0 | 86 | Int_t fBirkC0; // constants for Birk's Law implementation |
87 | Double_t fBirkC1; // constants for Birk's Law implementation | |
88 | Double_t fBirkC2; // constants for Birk's Law implementation | |
1963b290 | 89 | |
56088960 | 90 | static Double_t fgCapa ; // capacitor of the preamplifier for the raw RO signal |
91 | Double_t fHighCharge ; // high charge (to convert energy to charge) for the raw RO signal | |
92 | Double_t fHighGain ; // high gain for the raw RO signal | |
93 | Double_t fHighLowGainFactor ; // high to low gain factor for the raw RO signal | |
94 | Int_t fLowGainOffset ; // to separate high from low gain in the DDL | |
95 | static Int_t fgOrder ; // order of the gamma function for the RO signal | |
96 | static const Int_t fkTimeBins = 256 ; // number of sampling bins of the raw RO signal | |
97 | static Double_t fgTimeMax ; // maximum sampled time of the raw RO signal | |
98 | static Double_t fgTimePeak ; // peaking time of the raw RO signal | |
99 | static Double_t fgTimeTrigger ; // time of the trigger for the RO signal | |
2d5d9e60 | 100 | static Int_t fgDDLOffset; |
101 | static Int_t fgThreshold; | |
102 | static Int_t fgChannelsPerDDL; | |
2012850d | 103 | |
56088960 | 104 | ClassDef(AliEMCAL,9) // Electromagnetic calorimeter (base class) |
105 | ||
106 | } ; | |
2012850d | 107 | |
108 | #endif // ALIEMCAL_H |