change default max distance R to 0.1

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

diff --git a/EMCAL/AliEMCALRawUtils.h b/EMCAL/AliEMCALRawUtils.h
index da327fe0de6ca9f6fe5187631604d79a65984acf..e3324629c97a592d6d1f09fb1a7d379c721d540c 100644 (file)
--- a/EMCAL/AliEMCALRawUtils.h
+++ b/EMCAL/AliEMCALRawUtils.h
@@ -1,14 +1,12 @@
+// -*- mode: c++ -*-
 #ifndef ALIEMCALRAWUTILS_H
 #define ALIEMCALRAWUTILS_H
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice     */
 
 /* $Id$ */
-/* History of cvs commits:
- *
- * $Log$
- *
- */
+
+
 //_________________________________________________________________________
 //  Utility Class for handling Raw data
 //  Does all transitions from Digits to Raw and vice versa, 
@@ -22,46 +20,69 @@
 //*-- Author: Marco van Leeuwen (LBL)
 //
 #include "TObject.h" // for ROOT types
+#include <TString.h>
+
 
+class AliAnalysisManager;
+class AliCaloRawStreamV3;
+class AliAltroMapping;
 class TGraph;
-class TF1;
 class AliRawReader;
+class AliEMCALGeometry;
+class AliCaloCalibPedestal;
+class AliCaloRawAnalyzer;
+class AliCaloRawAnalyzerLMSOffline;
+class AliEMCALTriggerRawDigitMaker;
+class AliEMCALTriggerData;
+
+#include "AliCaloConstants.h"
+
 
 class AliEMCALRawUtils : public TObject {
  public:
-  AliEMCALRawUtils();
+  AliEMCALRawUtils(Algo::fitAlgorithm fitAlgo = Algo::kStandard);
   virtual ~AliEMCALRawUtils();
-
   void Digits2Raw();
-  void Raw2Digits(AliRawReader *reader);
-
-  // Signal shape parameters
-  Double_t GetRawFormatHighLowGainFactor() const {return fHighLowGainFactor ;}
-  static Int_t GetRawFormatOrder() { return fgOrder ; }   
-  static Int_t GetRawFormatTimeBins() { return fgkTimeBins ; }    
-  static Double_t GetRawFormatTimeMax() { return fgTimeMax ; }   
-  Double_t GetRawFormatTau() const { return fgTau ; }    
-  Double_t GetRawFormatTimeTrigger() const { return fgTimeTrigger ; }
-  Int_t GetRawFormatThreshold() const { return fgThreshold ; }       
-  Int_t GetRawFormatDDLPerSuperModule() const { return fgDDLPerSuperModule ; } 
-
-  // Signal shape functions
-  void FitRaw(TGraph * gSig, TF1* signalF, Double_t & amp, Double_t & time);
-  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;  
-
-  ClassDef(AliEMCALRawUtils,0)
-
- private:
-  Double_t fHighLowGainFactor ;         // high to low gain factor for the raw RO signal
-  static Int_t fgOrder ;                // order of the gamma function for the RO signal
-  static Double_t fgTau ;                // tau parameter of gamma function for the RO signal
-  static Double_t fgTimeTrigger ;       // time of the trigger for the RO signal 
+  void Raw2Digits(AliRawReader *reader, TClonesArray *digitsArr, const AliCaloCalibPedestal* pedbadmap,
+                                 TClonesArray *digitsTRG=0x0, AliEMCALTriggerData* trgData = 0x0);
+  void AddDigit(TClonesArray *digitsArr, Int_t id, Int_t lowGain, Float_t amp, Float_t time, Float_t chi2, Int_t ndf);
+  void TrimDigits(TClonesArray *digitsArr);
+  Int_t    GetNoiseThreshold()             const { return fNoiseThreshold; }
+  Int_t    GetNPedSamples()                const { return fNPedSamples; }
+  Bool_t   GetRemoveBadChannels() const {return fRemoveBadChannels;}
+  Int_t    GetFittingAlgorithm()  const {return fFittingAlgorithm; }
+  Float_t  GetTimeMax()           const {return fTimeMax ;}
+  Float_t  GetTimeMin()           const {return fTimeMin ;}
+  Bool_t   UseFALTRO()            const {return fUseFALTRO; }
+  void SetNoiseThreshold(Int_t val)                {fNoiseThreshold=val; }
+  void SetNPedSamples(Int_t val)                   {fNPedSamples=val; }
+  void SetRemoveBadChannels(Bool_t val)            {fRemoveBadChannels=val; }
+  void SetFittingAlgorithm(Int_t val) ;             
+  void SetTimeMin(Float_t t)                       {fTimeMin   = t          ;}
+  void SetTimeMax(Float_t t)                       {fTimeMax   = t          ;}
+  void SetFALTROUsage(Bool_t val)                  {fUseFALTRO=val; }
+  AliCaloRawAnalyzer *GetRawAnalyzer()  const { return fRawAnalyzer;}
+  virtual Option_t* GetOption() const { return fOption.Data(); }
+  void SetOption(const Option_t* opt) { fOption = opt; }
+  
+private:
+  AliEMCALRawUtils(const AliEMCALRawUtils& rawUtils);  //copy ctor
+  AliEMCALRawUtils& operator =(const AliEMCALRawUtils& rawUtils);
+  Int_t fNoiseThreshold;                // threshold to consider signal or noise
+  Int_t fNPedSamples;                   // number of samples to use in pedestal calculation
+  AliEMCALGeometry* fGeom;              // geometry
+  AliAltroMapping*  fMapping[4];        // only two for now
+  TString fOption;                      // option passed from Reconstructor
+  Bool_t  fRemoveBadChannels;           // select if bad channels are removed before fitting
+  Int_t   fFittingAlgorithm;            // select the fitting algorithm
+  Float_t fTimeMin;                     // minimum threshold for the time of the signal
+  Float_t fTimeMax;                     // maximum threshold for the time of the signal
+  Bool_t  fUseFALTRO;                  // use FALTRO and pass it to the digits
+  AliCaloRawAnalyzer *fRawAnalyzer;     // e.g. for sample selection for fits
+  AliEMCALTriggerRawDigitMaker* fTriggerRawDigitMaker; // trigger raw digit info
+ 
+  ClassDef(AliEMCALRawUtils,7)          // utilities for raw signal fitting
 
-  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 Int_t fgThreshold;             // threshold
-  static Int_t fgDDLPerSuperModule;     // number of DDL per SuperModule
 };
 
 #endif