Option for patching of the omega dalitz decay.

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

diff --git a/EMCAL/AliEMCALRawUtils.h b/EMCAL/AliEMCALRawUtils.h
index d554aeb22412083fc8cf387c1b3171d4de2cc208..6d456764bc65d66841f705b8b0b26743aa37ba58 100644 (file)
--- a/EMCAL/AliEMCALRawUtils.h
+++ b/EMCAL/AliEMCALRawUtils.h
@@ -40,9 +40,11 @@ class AliEMCALRawUtils : public TObject {
   AliEMCALRawUtils& operator =(const AliEMCALRawUtils& rawUtils);
 
   void Digits2Raw();
-  void Raw2Digits(AliRawReader *reader, TClonesArray *digitsArr, const AliCaloCalibPedestal* pedbadmap);
+  void Raw2Digits(AliRawReader *reader, TClonesArray *digitsArr, const AliCaloCalibPedestal* pedbadmap,
+                                 TClonesArray *digitsTRG=0x0);
 
   void AddDigit(TClonesArray *digitsArr, Int_t id, Int_t lowGain, Int_t amp, Float_t time);
+  void AddDigit(TClonesArray *digitsArr, Int_t id, Int_t timeSamples[], Int_t nSamples);
 
   // Signal shape parameters
   Double_t GetRawFormatHighLowGainFactor() const { return fHighLowGainFactor ;}
@@ -54,8 +56,9 @@ class AliEMCALRawUtils : public TObject {
   Double_t GetPedestalValue()  const {return fgPedestalValue;}
   Double_t GetFEENoise()       const {return fgFEENoise;}
 
-  Bool_t GetRemoveBadChannels()   const {return fRemoveBadChannels;}
-  Int_t  GetFittingAlgorithm()    const {return fFittingAlgorithm; }
+  Bool_t GetRemoveBadChannels() const {return fRemoveBadChannels;}
+  Int_t  GetFittingAlgorithm()  const {return fFittingAlgorithm; }
+  Bool_t UseFALTRO()            const {return fUseFALTRO; }
        
   void SetRawFormatHighLowGainFactor(Double_t val) {fHighLowGainFactor=val;}
   void SetRawFormatOrder(Int_t val)                {fOrder=val; }   
@@ -63,8 +66,9 @@ class AliEMCALRawUtils : public TObject {
   void SetNoiseThreshold(Int_t val)                {fNoiseThreshold=val; }
   void SetNPedSamples(Int_t val)                   {fNPedSamples=val; }
   void SetRemoveBadChannels(Bool_t val)            {fRemoveBadChannels=val; }
-  // not enough to set this variable to switch between algorithms, so comment it out for now..  void SetFittingAlgorithm(Int_t val)              {fFittingAlgorithm=val; }
-
+  void SetFittingAlgorithm(Int_t val) ;             
+  void SetFALTROUsage(Bool_t val)                  {fUseFALTRO=val; }
+       
   // set methods for fast fit simulation
   void SetFEENoise(Double_t val)                   {fgFEENoise = val;}
   void SetRawFormatTimeBins(Int_t val)             {fgTimeBins = val;}
@@ -87,7 +91,8 @@ class AliEMCALRawUtils : public TObject {
   void FitRaw(const Int_t firstTimeBin, const Int_t lastTimeBin, Float_t & amp, Float_t & time) const ;
   void FitParabola(const TGraph *gSig, Float_t & amp) const ; 
   static Double_t RawResponseFunction(Double_t *x, Double_t *par); 
-  Bool_t   RawSampledResponse(Double_t dtime, Double_t damp, Int_t * adcH, Int_t * adcL) const;  
+  static Double_t RawResponseFunctionLog(Double_t *x, Double_t *par); 
+  Bool_t   RawSampledResponse(Double_t dtime, Double_t damp, Int_t * adcH, Int_t * adcL, const Int_t keyErr=0) const;  
 
 
  private:
@@ -115,9 +120,11 @@ class AliEMCALRawUtils : public TObject {
 
   Bool_t fRemoveBadChannels;            // select if bad channels are removed before fitting
   Int_t  fFittingAlgorithm;             // select the fitting algorithm
+  Bool_t fUseFALTRO;                                   // use FALTRO and pass it to the digits
+       
   AliCaloRawAnalyzer *fRawAnalyzer;     // e.g. for sample selection for fits
 
-  ClassDef(AliEMCALRawUtils,5)          // utilities for raw signal fitting
+  ClassDef(AliEMCALRawUtils,6)          // utilities for raw signal fitting
 };
 
 #endif