OpenCL tracker bug fix

[u/mrichter/AliRoot.git] / TOF / AliTOFT0maker.h

diff --git a/TOF/AliTOFT0maker.h b/TOF/AliTOFT0maker.h
index 46c537d70c93a2005095ac21bcbe230d8eb108db..1530ba0ec9ecd72df388b675c95a7d134fd8fe54 100644 (file)
--- a/TOF/AliTOFT0maker.h
+++ b/TOF/AliTOFT0maker.h
@@ -3,7 +3,7 @@
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
 
-/* $Id:  $ */
+/* $Id: AliTOFT0maker.h,v 1.8 2010/01/19 16:32:20 noferini Exp $ */
 
 ///////////////////////////////////////////////
 //                                          //
@@ -14,43 +14,83 @@
 
 
 #include "TObject.h"
-#include "TString.h"
-#include "AliESDEvent.h"
-#include "AliStack.h"
 
-class AliTOFcalibHisto;
+class TH1F;
+class AliESDEvent;
+
+class AliESDpid;
 class AliTOFT0v1;
+class AliTOFcalib;
 
 class AliTOFT0maker : public TObject {
 public:
   
-  AliTOFT0maker() ;
+  AliTOFT0maker() ; // default constructor
+  AliTOFT0maker(AliESDpid *externalPID, AliTOFcalib *tofCalib=NULL); // overloaded constructor
   virtual ~AliTOFT0maker() ; // dtor
-  AliTOFT0maker(const AliTOFT0maker & t);
-  AliTOFT0maker & operator=(const AliTOFT0maker & t);
- 
-  void SetESDdata(Bool_t val=kTRUE){fESDswitch=val;};
-
-  // return (...[0]=event time -- ...[1]=sigma event time in ps -- mean event time for each fill) if you can subtruct the event time; return NULL if there is no event time
-  Double_t *RemakePID(AliESDEvent *esd,Double_t t0time=0.,Double_t t0sigma=1000.); // t0time and t0sigma in ps
+  
+  // return (fCalculated[0]=event time -- fCalculated[1]=sigma event time in ps -- fCalculated[2]=mean event time for each fill -- fCalculated[3]=number of tracks at the TOF level) if you can subtruct the event time; return NULL if there is no event time
+  Double_t *ComputeT0TOF(AliESDEvent *esd,Double_t t0time=0.,Double_t t0sigma=1000.); // t0time and t0sigma in ps
+  void ApplyT0TOF(AliESDEvent *esd);
+  Float_t GetExpectedSigma(Float_t mom, Float_t tof, Float_t mass);
+  Double_t  *GetT0p(Float_t p);
+  
+  void      SetTimeResolution(Double_t timeresolution){fTimeResolution=timeresolution;};// TOF timeresolution in [ps]
+  Double_t  GetTimeResolution() const {return fTimeResolution;} // Get TOF Time Resolution
+  void SetT0FillWidth(Float_t width){if(width > 50) fT0width = width; else fT0width=150;}; // in ps
+  
+  void LoadChannelMap(const char *filename="$ALICE_ROOT/TOF/enableMap.104892.root"); //load the enable channel map
+  void ApplyMask(AliESDEvent * const esd); // Apply the channel mask
+  
+  void SetNoTOFT0(Bool_t status=kTRUE){fNoTOFT0=status;}; // disable the TOF T0 info
+  void SetMaskOffChannel(Bool_t status=kTRUE){fKmask=status;}; // switch for the map off channel
+  
+  Float_t  TuneForMC(AliESDEvent *esd); // set enabled channeld, add a smeared t0, add a TOF smearing, return true(smeared) T0 event
+  
+  void SetT0spread(Float_t t0spread){fT0spreadExt=t0spread;}; // get T0spread
+  Float_t GetT0spread() const {return fT0spreadExt;} // get T0spread
 
-  void      SetTimeResolution(Double_t timeresolution){fTimeResolution=timeresolution;};// TOF timeresolution in [s] e.g. for 120 ps -> 1.2e-10
-  Double_t  GetTimeResolution() const {return fTimeResolution;}
+  void SetT0fill(Float_t t0fill){fT0fillExt=t0fill;};
   
+  void WriteInESD(AliESDEvent *esd);
+
+  void SetTOFT0algorithm(Int_t flag=1){fTOFT0algorithm=flag;};
+  Int_t GetTOFT0algorithm() const {return fTOFT0algorithm;};
+
  private:
-  void TakeTimeRawCorrection(AliESDEvent * const esd);
-  void RemakeTOFpid(AliESDEvent *esd,Float_t timezero);
-  Double_t GetT0Fill(Int_t nrun) const ;
+  
+  AliTOFT0maker(const AliTOFT0maker &);
+  AliTOFT0maker & operator=(const AliTOFT0maker &);
+  
+  void SetTOFResponse();
 
-  AliTOFcalibHisto *fCalib; // TOF calibration object pointer
+  AliTOFT0v1 *fT0TOF; // T0-TOF
+  AliESDpid *fPIDesd; // PID esd
+  Bool_t fExternalPIDFlag; // external PID flag
+  AliTOFcalib *fTOFcalib; // TOF calibration
 
-  Bool_t fESDswitch; // if you want take the ESD time instead of the raw + time slewing correction
+  Bool_t fNoTOFT0;   // switch to avoid T0-TOF is used
 
-  Double_t fTimeResolution;  // global time resolution used to calculate T0
+  Int_t fNmomBins;
 
+  Double_t fCalculated[10]; // contains the parameters with the event time
+
+  Double_t fT0cur[2]; // current T0 and T0 sigma
+ 
+  Double_t fTimeResolution;  // global time resolution used to calculate T0
+  
   Float_t fT0sigma; // T0 resolution
   
-  ClassDef(AliTOFT0maker,1);  // Calculate the time zero using TOF detector */
+  TH1F *fHmapChannel; // histo with the channel map
+  Bool_t fKmask; // switch if you want apply a channel filter
+  
+  Float_t fT0width; // T0 FILL width
+
+  Float_t fT0spreadExt;
+  Float_t fT0fillExt; // t0spread if set 
+  Int_t  fTOFT0algorithm; // select algorithm (1 default, 2 speed up)
+
+  ClassDef(AliTOFT0maker,3);  // Calculate the time zero using TOF detector */
   
 };