* See cxx source for full Copyright notice */
/* $Id$ */
-
+/* History of cvs commits:
+ *
+ * $Log$
+ *
+ */
//_________________________________________________________________________
// Base Class for EMCAL
+// holds all geant information of
+// materials, etc.
//
//*-- Author: Yves Schutz (SUBATECH)
class TTask ;
class TFolder ;
class TRandom ;
+class TGraph;
+class TF1;
// --- AliRoot header files ---
-
+class AliRawReader;
#include "AliDetector.h"
#include "AliEMCALGeometry.h"
+#include "AliEMCALTrigger.h"
+
class AliEMCAL : public AliDetector {
public:
-
+
AliEMCAL();
AliEMCAL(const char* name, const char* title="");
- AliEMCAL(const AliEMCAL& emcal) : AliDetector(emcal) {
- // cpy ctor: no implementation yet
- // requested by the Coding Convention
- Fatal("cpy ctor", "not implemented") ;
- }
+
virtual ~AliEMCAL() ;
- virtual void AddHit(Int_t, Int_t*, Float_t *) const{
+ virtual void AddHit(Int_t, Int_t*, Float_t *) {
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 )") ;
}
- virtual void Copy(AliEMCAL & emcal) ;
virtual AliDigitizer* CreateDigitizer(AliRunDigitizer* manager) const;
virtual void CreateMaterials() ;
+ // virtual void
virtual void Digits2Raw();
+
+ using AliDetector::Raw2Digits;
+ virtual void Raw2Digits(AliRawReader *reader);
+
virtual void FinishRun() {}
virtual AliEMCALGeometry * GetGeometry() const
- {return AliEMCALGeometry::GetInstance(GetTitle(),"") ; }
+ {return AliEMCALGeometry::GetInstance(GetTitle(),"") ; }
virtual void Hits2SDigits();
virtual Int_t IsVersion(void) const = 0 ;
- Int_t GetRawFormatHighGainFactor() const { return fHighGainFactor ; }
- Int_t GetRawFormatHighGainOffset() const { return fHighGainOffset ; }
- Int_t GetRawFormatTimeBins() const { return fkTimeBins ; }
- Double_t GetRawFormatTimeMax() const { return fTimeMax ; }
- Double_t GetRawFormatTimePeak() const { return fTimePeak ; }
- Double_t GetRawFormatTimeRes() const { return fTimeRes ; }
+
+ virtual AliTriggerDetector* CreateTriggerDetector() const
+ { return new AliEMCALTrigger(); }
+
+ // Raw Read Out
+ Double_t GetRawFormatCapa() const { return fgCapa ; }
+ Double_t GetRawFormatHighCharge() const { return fHighCharge ; }
+ Double_t GetRawFormatHighGain() const { return fHighGain ; }
+ Double_t GetRawFormatHighLowGainFactor() const { return fHighLowGainFactor ; }
+ Double_t GetRawFormatLowCharge() const { return ( fHighCharge * fHighLowGainFactor ) ; }
+ Double_t GetRawFormatLowGain() const { return ( fHighGain / fHighLowGainFactor ) ; }
+ Int_t GetRawFormatLowGainOffset() const { return fLowGainOffset ; }
+ Int_t GetRawFormatOrder() const { return fgOrder ; }
+ Int_t GetRawFormatTimeBins() const { return fgkTimeBins ; }
+ Double_t GetRawFormatTimeMax() const { return fgTimeMax ; }
+ Double_t GetRawFormatTimePeak() const { return fgTimePeak ; }
+ Double_t GetRawFormatTimeTrigger() const { return fgTimeTrigger ; }
+ Int_t GetRawFormatThreshold() const { return fgThreshold ; }
+ Int_t GetRawFormatDDLPerSuperModule() const { return fgDDLPerSuperModule ; }
+ static Double_t RawResponseFunctionMax(Double_t charge, Double_t gain) ;
+ Bool_t RawSampledResponse(Double_t dtime, Double_t damp, Int_t * adcH, Int_t * adcL) const ;
+ //
virtual AliLoader* MakeLoader(const char* topfoldername);
- static Double_t RawResponseFunction(Double_t *x, Double_t *par) ;
- virtual void SetTreeAddress() ;
virtual const TString Version() const {return TString(" ") ; }
- AliEMCAL & operator = (const AliEMCAL & /*rvalue*/) {
- Fatal("operator =", "not implemented") ; return *this ; }
-
+
protected:
+
+ static Double_t RawResponseFunction(Double_t *x, Double_t *par) ;
+ void FitRaw(Bool_t lowGainFlag, TGraph * gLowGain, TGraph * gHighGain, TF1* signalF, Double_t & energy, Double_t & time) ;
- Bool_t RawSampledResponse(const Float_t dtime, const Int_t damp, Int_t * adcH, Int_t * adcL) const ;
+ void Init(void); //initializes some params
Int_t fBirkC0; // constants for Birk's Law implementation
Double_t fBirkC1; // constants for Birk's Law implementation
Double_t fBirkC2; // constants for Birk's Law implementation
- Int_t fHighGainFactor ; // High gain attenuation factor of the raw RO signal
- Int_t fHighGainOffset ; // offset added to the module id to distinguish high and low gain data
- static const Int_t fkTimeBins = 256 ; // number of sampling bins of the raw RO signal
- Double_t fTimeMax ; // maximum sampled time of the raw RO signal
- Double_t fTimePeak ; // peaking time of the raw RO signal
- Double_t fTimeRes ; // decay rime width of the raw RO signal
- ClassDef(AliEMCAL,8) // Electromagnetic calorimeter (base class)
+ static Double_t fgCapa ; // capacitor of the preamplifier for the raw RO signal
+ Double_t fHighCharge ; // high charge (to convert energy to charge) for the raw RO signal
+ Double_t fHighGain ; // high gain for the raw RO signal
+ Double_t fHighLowGainFactor ; // high to low gain factor for the raw RO signal
+ Int_t fLowGainOffset ; // to separate high from low gain in the DDL
+ static Int_t fgOrder ; // order of the gamma function for the RO signal
+ static const Int_t fgkTimeBins = 256 ; // number of sampling bins of the raw RO signal
+ static Double_t fgTimeMax ; // maximum sampled time of the raw RO signal
+ static Double_t fgTimePeak ; // peaking time of the raw RO signal
+ static Double_t fgTimeTrigger ; // time of the trigger for the RO signal
+ static Int_t fgThreshold; // threshold
+ static Int_t fgDDLPerSuperModule; // number of DDL per SuperModule
+
+private:
+ AliEMCAL(const AliEMCAL& emcal);
+ AliEMCAL & operator = (const AliEMCAL & /*rvalue*/);
-} ;
+ ClassDef(AliEMCAL,9) // Electromagnetic calorimeter (base class)
+
+ } ;
#endif // ALIEMCAL_H
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff