default geant cuts for EMCAL's material: cutele=cutgam=100KeV

[u/mrichter/AliRoot.git] / EMCAL / AliEMCAL.h

diff --git a/EMCAL/AliEMCAL.h b/EMCAL/AliEMCAL.h
index 61c81b3facfb42f398f836da05dda2438d3a85df..82e937732d63ce36c69e33ae9ff8df5031d0c299 100644 (file)
--- a/EMCAL/AliEMCAL.h
+++ b/EMCAL/AliEMCAL.h
@@ -10,21 +10,17 @@
 //                  
 //*-- Author: Yves Schutz (SUBATECH)
 
-#include <stdlib.h>
-
 // --- ROOT system ---
 
 class TString ;
 class TTask ;
 class TFolder ;
+class TRandom ; 
 
 // --- AliRoot header files ---
 
 #include "AliDetector.h"
-//class AliDetector;
-class AliEMCALGeometry ; 
-//class AliEMCALQAChecker ;
-
+#include "AliEMCALGeometry.h" 
 class AliEMCAL : public AliDetector {
 
  public:
@@ -37,40 +33,64 @@ class AliEMCAL : public AliDetector {
     Fatal("cpy ctor", "not implemented") ;  
   }
   virtual ~AliEMCAL() ; 
-  virtual void   AddHit(Int_t, Int_t*, Float_t *) {
-    // do not use this definition but the one below
-    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   AddHit(Int_t, Int_t*, Float_t *) const{
+    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  AddHit( Int_t shunt, Int_t primary, Int_t track, 
-                        Int_t id, Float_t *hits ) = 0 ;
+  virtual void  Copy(AliEMCAL & emcal) ; 
+  virtual AliDigitizer* CreateDigitizer(AliRunDigitizer* manager) const;
   virtual void  CreateMaterials() ;   
-  virtual void  FinishRun() {WriteQA();}                  
-  virtual AliEMCALGeometry * GetGeometry() const ;   
-  virtual Int_t   IsVersion(void) const = 0 ; 
-  //AliEMCALQAChecker * QAChecker() const {return fQATask;}  
-  virtual void  SetTreeAddress() ;
-  virtual TTree * TreeQA() const {return fTreeQA; }                
-  virtual const TString Version() const {return TString(" ") ; }  
-  virtual void WriteQA() ; 
-  AliEMCAL & operator = (const AliEMCAL & rvalue)  {
-    // assignement operator requested by coding convention
-    // but not needed
-    Fatal("operator =", "not implemented") ;  
-    return *this ; 
-  }
- 
+  //  virtual void  
+  virtual void  Digits2Raw();
+  virtual void  FinishRun() {}                  
+  virtual AliEMCALGeometry * GetGeometry() const 
+  {return AliEMCALGeometry::GetInstance(GetTitle(),"") ;  }   
+  virtual void    Hits2SDigits();
+  virtual Int_t   IsVersion(void) const = 0 ;   
+  // 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 fkTimeBins ; }    
+  Double_t GetRawFormatTimeMax() const { return fgTimeMax ; }   
+  Double_t GetRawFormatTimePeak() const { return fgTimePeak ; }    
+  Double_t GetRawFormatTimeTrigger() const { return fgTimeTrigger ; }   
+  static Double_t RawResponseFunctionMax(Double_t charge, Double_t gain) ;
+  //  
   virtual AliLoader* MakeLoader(const char* topfoldername);
-  
+  virtual void SetTreeAddress() ;              
+  virtual const TString Version() const {return TString(" ") ; }   
+  AliEMCAL & operator = (const AliEMCAL & /*rvalue*/)  {
+    Fatal("operator =", "not implemented") ;  return *this ; }
+
 protected:
+  
+  friend class AliEMCALGetter;
+  static Double_t RawResponseFunction(Double_t *x, Double_t *par) ; 
+  Bool_t   RawSampledResponse(const Double_t dtime, const Double_t damp, Int_t * adcH, Int_t * adcL) const ; 
+
 
-  //AliEMCALQAChecker * fQATask ; //! PHOS checkers container
-  TTree * fTreeQA ;            // the QA tree that contains the alarms
-  AliEMCALGeometry * fGeom ;   // the geometry object
+  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
 
-  ClassDef(AliEMCAL,4) // 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 fkTimeBins = 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 
 
-} ;
+  ClassDef(AliEMCAL,9) // Electromagnetic calorimeter (base class)
+    
+    } ;
 
 #endif // ALIEMCAL_H