EMCAL/AliEMCAL.h

   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 /* History of cvs commits:
   8  *
   9  * $Log$
  10  * Revision 1.41  2007/02/05 10:43:25  hristov
  11  * Changes for correct initialization of Geant4 (Mihaela)
  12  *
  13  * Revision 1.40  2006/12/05 17:19:26  gustavo
  14  * Updated AliEMCAL::Digits2Raw, reads first provisional RCU mapping files to make Raw data with new AliCaloAltroMapping and AliCaloRawStream
  15  *
  16  *
  17  */
  18 //_________________________________________________________________________
  19 //  Base Class for EMCAL
  20 //  holds all geant information of
  21 //  materials, etc.
  22 //
  23 //*-- Author: Yves Schutz (SUBATECH)
  24
  25 // --- ROOT system ---
  26
  27 class TString ;
  28 class TTask ;
  29 class TFolder ;
  30 class TRandom ;
  31 class TGraph;
  32 class TF1;
  33
  34 // --- AliRoot header files ---
  35 class AliRawReader;
  36 #include "AliDetector.h"
  37 #include "AliEMCALGeometry.h"
  38 #include "AliEMCALTrigger.h"
  39
  40 class AliEMCAL : public AliDetector {
  41
  42  public:
  43
  44   AliEMCAL();
  45   AliEMCAL(const char* name, const char* title="");
  46
  47   virtual ~AliEMCAL() ;
  48   virtual void   AddHit(Int_t, Int_t*, Float_t *) {
  49     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 )") ;
  50   }
  51   virtual AliDigitizer* CreateDigitizer(AliRunDigitizer* manager) const;
  52   virtual void  CreateMaterials() ;
  53   virtual void  Digits2Raw();
  54
  55   using AliDetector::Raw2Digits;
  56   virtual void  Raw2Digits(AliRawReader *reader);
  57
  58   virtual void  FinishRun() {}
  59   virtual AliEMCALGeometry * GetGeometry() const
  60     {return AliEMCALGeometry::GetInstance(GetTitle(),"") ;  }
  61   virtual void    Hits2SDigits();
  62   //  virtual void    Init();
  63   virtual Int_t   IsVersion(void) const = 0 ;
  64
  65   virtual AliTriggerDetector* CreateTriggerDetector() const
  66     { return new AliEMCALTrigger(); }
  67
  68   // Raw Read Out
  69   Double_t GetRawFormatCapa() const { return fgCapa ; }
  70   Double_t GetRawFormatHighCharge() const { return fHighCharge ; }
  71   Double_t GetRawFormatHighGain() const { return fHighGain ; }
  72   Double_t GetRawFormatHighLowGainFactor() const { return fHighLowGainFactor ; }
  73   Double_t GetRawFormatLowCharge() const { return ( fHighCharge *  fHighLowGainFactor ) ; }
  74   Double_t GetRawFormatLowGain() const { return ( fHighGain / fHighLowGainFactor ) ; }
  75   Int_t GetRawFormatLowGainOffset() const { return fLowGainOffset ; }
  76   Int_t GetRawFormatOrder() const { return fgOrder ; }
  77   Int_t GetRawFormatTimeBins() const { return fgkTimeBins ; }
  78   Double_t GetRawFormatTimeMax() const { return fgTimeMax ; }
  79   Double_t GetRawFormatTimePeak() const { return fgTimePeak ; }
  80   Double_t GetRawFormatTimeTrigger() const { return fgTimeTrigger ; }
  81   Int_t GetRawFormatThreshold() const { return fgThreshold ; }
  82   Int_t GetRawFormatDDLPerSuperModule() const { return fgDDLPerSuperModule ; }
  83   static Double_t RawResponseFunctionMax(Double_t charge, Double_t gain) ;
  84   Bool_t   RawSampledResponse(Double_t dtime, Double_t damp, Int_t * adcH, Int_t * adcL) const ;
  85   //
  86   virtual AliLoader* MakeLoader(const char* topfoldername);
  87   virtual const TString Version() const {return TString(" ") ; }
  88
  89 protected:
  90
  91   static Double_t RawResponseFunction(Double_t *x, Double_t *par) ;
  92   void FitRaw(Bool_t lowGainFlag, TGraph * gLowGain, TGraph * gHighGain, TF1* signalF, Double_t & energy, Double_t & time) ;
  93
  94   void InitConstants();  //initializes some params
  95   void DefineMediumParameters();  // define tracking medium parameters
  96
  97   Int_t fBirkC0;    // constants for Birk's Law implementation
  98   Double_t fBirkC1; // constants for Birk's Law implementation
  99   Double_t fBirkC2; // constants for Birk's Law implementation
 100
 101   static Double_t fgCapa ;              // capacitor of the preamplifier for the raw RO signal
 102   Double_t fHighCharge ;                // high charge (to convert energy to charge) for the raw RO signal
 103   Double_t fHighGain ;                  // high gain for the raw RO signal
 104   Double_t fHighLowGainFactor ;         // high to low gain factor for the raw RO signal
 105   Int_t    fLowGainOffset ;             // to separate high from low gain in the DDL
 106   static Int_t fgOrder ;                // order of the gamma function for the RO signal
 107   static const Int_t fgkTimeBins = 256 ; // number of sampling bins of the raw RO signal
 108   static Double_t fgTimeMax ;           // maximum sampled time of the raw RO signal
 109   static Double_t fgTimePeak ;          // peaking time of the raw RO signal
 110   static Double_t fgTimeTrigger ;       // time of the trigger for the RO signal
 111   static Int_t fgThreshold;             // threshold
 112   static Int_t fgDDLPerSuperModule;        // number of DDL per SuperModule
 113
 114 private:
 115   AliEMCAL(const AliEMCAL& emcal);
 116   AliEMCAL & operator = (const AliEMCAL & /*rvalue*/);
 117
 118   ClassDef(AliEMCAL,9) // Electromagnetic calorimeter (base class)
 119
 120     } ;
 121
 122 #endif // ALIEMCAL_H