Change for a realistic time distribution of the raw signal
[u/mrichter/AliRoot.git] / PHOS / AliPHOS.h
1 #ifndef ALIPHOS_H
2 #define ALIPHOS_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice     */
5 /* $Id$ */
6
7 //_________________________________________________________________________
8 //  Base Class for PHOS     
9 //                  
10 //*-- Author: Laurent Aphecetche & Yves Schutz (SUBATECH)
11
12
13 // --- ROOT system ---
14 class TString ; 
15 class TTask ;
16 class TFolder ;
17 class TTree ; 
18 class TRandom ; 
19
20 // --- AliRoot header files ---
21 #include "AliDetector.h" 
22 #include "AliPHOSGeometry.h" 
23 class AliPHOSQAChecker ;
24
25 class AliPHOS : public AliDetector {
26
27 public:
28
29   AliPHOS() ;
30   AliPHOS(const char* name, const char* title="") ;  
31   AliPHOS(AliPHOS & phos) : AliDetector(phos) {
32     Copy(*this) ; 
33   }
34   virtual ~AliPHOS() ; 
35   virtual void   AddHit(Int_t, Int_t*, Float_t *) {
36     // do not use this definition but the one below
37     Fatal("AddHit(Int_t, Int_t*, Float_t *)", "do not use") ;
38     
39   }
40   virtual void   AddHit( Int_t shunt, Int_t primary, Int_t track, 
41                          Int_t id, Float_t *hits ) = 0 ;   
42   virtual void Copy(AliPHOS & phos) ; 
43   virtual AliDigitizer* CreateDigitizer(AliRunDigitizer* manager) const;
44   virtual void  CreateMaterials() ;            
45   virtual void  Digits2Raw();
46   virtual void  FinishRun() {WriteQA();}
47   virtual AliPHOSGeometry * GetGeometry() const 
48   {return AliPHOSGeometry::GetInstance(GetTitle(),"") ;  }
49   virtual void    Hits2SDigits();
50   virtual Int_t   IsVersion(void) const = 0 ;  
51   // Raw Read Out
52   Double_t GetRawFormatCapa() const { return fgCapa ; }   
53   Double_t GetRawFormatHighCharge() const { return fHighCharge ; }  
54   Double_t GetRawFormatHighGain() const { return fHighGain ; }  
55   Double_t GetRawFormatHighLowGainFactor() const { return fHighLowGainFactor ; }  
56   Double_t GetRawFormatLowCharge() const { return ( fHighCharge *  fHighLowGainFactor ) ; }  
57   Double_t GetRawFormatLowGain() const { return ( fHighGain / fHighLowGainFactor ) ; }  
58   Int_t GetRawFormatLowGainOffset() const { return fLowGainOffset ; }  
59   Int_t GetRawFormatOrder() const { return fgOrder ; }   
60   Int_t GetRawFormatTimeBins() const { return fkTimeBins ; }    
61   Double_t GetRawFormatTimeMax() const { return fgTimeMax ; }   
62   Double_t GetRawFormatTimePeak() const { return fgTimePeak ; }    
63   Double_t GetRawFormatTimeTrigger() const { return fgTimeTrigger ; }   
64   static Double_t RawResponseFunction(Double_t *x, Double_t *par) ; 
65   static Double_t RawResponseFunctionMax(Double_t charge, Double_t gain) ;
66   //
67   virtual AliLoader* MakeLoader(const char* topfoldername);
68   AliPHOSQAChecker * QAChecker() {return fQATask;}  
69   virtual void    SetTreeAddress();   
70   virtual TTree * TreeQA() const {return fTreeQA; } 
71   virtual const TString Version() const {return TString(" ") ; } 
72   virtual void WriteQA() ; 
73   AliPHOS & operator = (const AliPHOS & /*rvalue*/)  {
74     Fatal("operator =", "not implemented") ; return *this ; }
75
76
77 protected:
78
79     Bool_t   RawSampledResponse(const Double_t dtime, const Double_t damp, Int_t * adcH, Int_t * adcL) const ; 
80
81
82   AliPHOSQAChecker * fQATask ;          //! PHOS checkers container
83   TTree * fTreeQA ;                     // the QA tree that contains the alarms
84   static Double_t fgCapa ;              // capacitor of the preamplifier for the raw RO signal
85   Double_t fHighCharge ;                // high charge (to convert energy to charge) for the raw RO signal
86   Double_t fHighGain ;                  // high gain for the raw RO signal
87   Double_t fHighLowGainFactor ;         // high to low gain factor for the raw RO signal
88   Int_t    fLowGainOffset ;             // to separate high from low gain in the DDL
89   static Int_t fgOrder ;                // order of the gamma function for the RO signal
90   static const Int_t fkTimeBins = 256 ; // number of sampling bins of the raw RO signal  
91   static Double_t fgTimeMax ;           // maximum sampled time of the raw RO signal                             
92   static Double_t fgTimePeak ;          // peaking time of the raw RO signal                                    
93   static Double_t fgTimeTrigger ;       // time of the trigger for the RO signal 
94                                         
95   ClassDef(AliPHOS,4) // Photon Spectrometer Detector (base class)
96 } ;
97
98 #endif // ALIPHOS_H