/* 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 $ */
///////////////////////////////////////////////
// //
#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 */
};