ANALYSIS/AliTriggerAnalysis.h

   1 /* $Id: AliTriggerAnalysis.h 35782 2009-10-22 11:54:31Z jgrosseo $ */
   2
   3 #ifndef ALITRIGGERANALYSIS_H
   4 #define ALITRIGGERANALYSIS_H
   5
   6 #include <TObject.h>
   7
   8 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   9  * See cxx source for full Copyright notice                               */
  10
  11 //-------------------------------------------------------------------------
  12 //                      Implementation of   Class AliTriggerAnalysis
  13 //   This class provides function to check if events have been triggered based on the data in the ESD
  14 //   The trigger bits, trigger class inputs and only the data (offline trigger) can be used
  15 //   Origin: Jan Fiete Grosse-Oetringhaus, CERN
  16 //-------------------------------------------------------------------------
  17
  18 class AliESDEvent;
  19 class AliESDtrackCuts;
  20 class TH1F;
  21 class TH2F;
  22 class TCollection;
  23 class TMap;
  24
  25 class AliTriggerAnalysis : public TObject
  26 {
  27   public:
  28     enum Trigger { kAcceptAll = 1, kMB1 = 2, kMB2, kMB3, kSPDGFO, kSPDGFOBits, kV0A, kV0C, kV0OR, kV0AND,
  29                    kV0ABG, kV0CBG, kZDC, kZDCA, kZDCC, kZNA, kZNC, kZNABG, kZNCBG, kFMDA, kFMDC, kFPANY, kNSD1, kMB1Prime,
  30                    kSPDGFOL0, kSPDGFOL1, kZDCTDCA, kZDCTDCC, kZDCTime, kCTPV0A, kCTPV0C, kTPCLaserWarmUp, kSPDClsVsTrkBG,
  31                    kCentral,kSemiCentral, kT0, kT0BG, kT0Pileup, kEMCAL,kTPCHVdip,
  32                    kTRDHCO, kTRDHJT, kTRDHSE, kTRDHQU, kTRDHEE,
  33                    kStartOfFlags = 0x0100, kOfflineFlag = 0x8000, kOneParticle = 0x10000, kOneTrack = 0x20000}; // MB1, MB2, MB3 definition from ALICE-INT-2005-025
  34     enum AliceSide { kASide = 1, kCSide, kCentralBarrel };
  35     enum V0Decision { kV0Invalid = -1, kV0Empty = 0, kV0BB, kV0BG, kV0Fake };
  36     enum T0Decision { kT0Invalid = -1, kT0Empty = 0, kT0BB, kT0DecBG, kT0DecPileup };
  37
  38     AliTriggerAnalysis();
  39     virtual ~AliTriggerAnalysis();
  40
  41     void EnableHistograms();
  42     void SetAnalyzeMC(Bool_t flag = kTRUE) { fMC = flag; }
  43
  44     Bool_t IsTriggerFired(const AliESDEvent* aEsd, Trigger trigger);
  45     Int_t EvaluateTrigger(const AliESDEvent* aEsd, Trigger trigger);
  46
  47     // using trigger bits in ESD
  48     Bool_t IsTriggerBitFired(const AliESDEvent* aEsd, Trigger trigger) const;
  49     Bool_t IsTriggerBitFired(ULong64_t triggerMask, Trigger trigger) const;
  50     Bool_t IsTriggerBitFired(const AliESDEvent* aEsd, ULong64_t tclass) const;
  51
  52     // using ESD data from detectors
  53     Bool_t IsOfflineTriggerFired(const AliESDEvent* aEsd, Trigger trigger);
  54
  55     // using trigger classes in ESD
  56     Bool_t IsTriggerClassFired(const AliESDEvent* aEsd, const Char_t* tclass) const;
  57
  58     // some "raw" trigger functions
  59     Int_t SPDFiredChips(const AliESDEvent* aEsd, Int_t origin, Bool_t fillHists = kFALSE, Int_t layer = 0);
  60     Bool_t SPDGFOTrigger(const AliESDEvent* aEsd, Int_t origin);
  61     Bool_t IsSPDClusterVsTrackletBG(const AliESDEvent* esd, Bool_t fillHists = kFALSE);
  62     V0Decision V0Trigger(const AliESDEvent* aEsd, AliceSide side, Bool_t online, Bool_t fillHists = kFALSE);
  63     T0Decision T0Trigger(const AliESDEvent* aEsd, Bool_t online, Bool_t fillHists = kFALSE);
  64     Bool_t ZDCTrigger   (const AliESDEvent* aEsd, AliceSide side) const;
  65     Bool_t ZDCTDCTrigger(const AliESDEvent* aEsd, AliceSide side, Bool_t useZN=kTRUE, Bool_t useZP=kFALSE, Bool_t fillHists=kFALSE) const;
  66     Bool_t ZDCTimeTrigger(const AliESDEvent *aEsd, Bool_t fillHists=kFALSE) const;
  67     Bool_t ZDCTimeBGTrigger(const AliESDEvent *aEsd, AliceSide side) const;
  68     Bool_t FMDTrigger(const AliESDEvent* aEsd, AliceSide side);
  69     Int_t SSDClusters(const AliESDEvent* aEsd);
  70     Bool_t EMCALCellsTrigger(const AliESDEvent *aEsd);
  71     static const char* GetTriggerName(Trigger trigger);
  72
  73     Bool_t IsLaserWarmUpTPCEvent(const AliESDEvent* esd);
  74         Bool_t IsHVdipTPCEvent(const AliESDEvent* esd);
  75     Bool_t TRDTrigger(const AliESDEvent* esd, Trigger trigger);
  76
  77     void FillHistograms(const AliESDEvent* aEsd);
  78     void FillTriggerClasses(const AliESDEvent* aEsd);
  79
  80     void SetSPDGFOThreshhold(Int_t t) { fSPDGFOThreshold = t; }
  81     void SetSPDGFOEfficiency(TH1F* hist) { fSPDGFOEfficiency = hist; }
  82     void SetSPDClustersVsTrackletsParameters(Float_t a, Float_t b) { fASPDCvsTCut = a; fBSPDCvsTCut =b;}
  83     void SetV0TimeOffset(Float_t offset) { fV0TimeOffset = offset; }
  84     void SetV0AdcThr(Float_t thr) { fV0AdcThr = thr; }
  85     void SetV0HwPars(Float_t thr, Float_t winLow, Float_t winHigh) { fV0HwAdcThr = thr; fV0HwWinLow = winLow; fV0HwWinHigh = winHigh; }
  86     void SetFMDThreshold(Float_t low, Float_t hit) { fFMDLowCut = low; fFMDHitCut = hit; }
  87     void SetDoFMD(Bool_t flag = kTRUE) {fDoFMD = flag;}
  88     void SetZDCCutParams(Float_t refSum, Float_t refDelta, Float_t sigmaSum, Float_t sigmaDelta) { fZDCCutRefSum = refSum; fZDCCutRefDelta = refDelta; fZDCCutSigmaSum = sigmaSum; fZDCCutSigmaDelta = sigmaDelta; }
  89     void SetCorrZDCCutParams(Float_t refSum, Float_t refDelta, Float_t sigmaSum, Float_t sigmaDelta) { fZDCCutRefSumCorr = refSum; fZDCCutRefDeltaCorr = refDelta; fZDCCutSigmaSumCorr = sigmaSum; fZDCCutSigmaDeltaCorr = sigmaDelta; }
  90     void SetZNCorrCutParams(Float_t znaTimeCorrMin, Float_t znaTimeCorrMax, Float_t zncTimeCorrMin, Float_t zncTimeCorrMax)
  91     { fZDCCutZNATimeCorrMin = znaTimeCorrMin; fZDCCutZNATimeCorrMax = znaTimeCorrMax;
  92       fZDCCutZNCTimeCorrMin = zncTimeCorrMin; fZDCCutZNCTimeCorrMax = zncTimeCorrMax; }
  93
  94     void SetTRDTriggerParameters(Float_t ptHSE, UChar_t pidHSE, Float_t ptHQU, UChar_t pidHQU, Float_t ptHEE, UChar_t pidHEE, UChar_t minSectorHEE, UChar_t maxSectorHEE, Float_t ptHJT, UChar_t nHJT) {
  95       fTRDptHSE = ptHSE; fTRDpidHSE = pidHSE;
  96       fTRDptHQU = ptHQU; fTRDpidHQU = pidHQU;
  97       fTRDptHEE = ptHEE; fTRDpidHEE = pidHEE;
  98       fTRDminSectorHEE = minSectorHEE; fTRDmaxSectorHEE = maxSectorHEE;
  99       fTRDptHJT = ptHJT; fTRDnHJT = nHJT;
 100     }
 101     void SetUseTRDMatchCut(Bool_t flag){ fUseTRDMatchCut=flag; }
 102     Bool_t GetUseTRDMatchCut(){ return fUseTRDMatchCut; }
 103
 104     Int_t GetSPDGFOThreshhold() const { return fSPDGFOThreshold; }
 105     Float_t GetV0TimeOffset() const { return fV0TimeOffset; }
 106     Float_t GetV0AdcThr()     const { return fV0AdcThr; }
 107     Float_t GetFMDLowThreshold() const { return fFMDLowCut; }
 108     Float_t GetFMDHitThreshold() const { return fFMDHitCut; }
 109     TMap * GetTriggerClasses() const { return fTriggerClasses;}
 110
 111
 112     virtual Long64_t Merge(TCollection* list);
 113     void SaveHistograms() const;
 114
 115     void PrintTriggerClasses() const;
 116     void SetESDTrackCuts(AliESDtrackCuts* cuts) { fEsdTrackCuts = cuts;}
 117     AliESDtrackCuts* GetESDTrackCuts() const  {return fEsdTrackCuts;}
 118
 119     void SetTPCOnly(Bool_t bTPCOnly) {fTPCOnly = bTPCOnly;}
 120     Bool_t GetTPCOnly() const {return fTPCOnly;}
 121
 122   protected:
 123     Bool_t IsL0InputFired(const AliESDEvent* aEsd, UInt_t input) const;
 124     Bool_t IsL1InputFired(const AliESDEvent* aEsd, UInt_t input) const;
 125     Bool_t IsL2InputFired(const AliESDEvent* aEsd, UInt_t input) const;
 126     Bool_t IsInputFired(const AliESDEvent* aEsd, Char_t level, UInt_t input) const;
 127
 128     Float_t V0CorrectLeadingTime(Int_t i, Float_t time, Float_t adc, Int_t runNumber) const;
 129     Float_t V0LeadingTimeWeight(Float_t adc) const;
 130
 131     Int_t FMDHitCombinations(const AliESDEvent* aEsd, AliceSide side, Bool_t fillHists = kFALSE);
 132
 133     Int_t fSPDGFOThreshold;         // number of chips to accept a SPD GF0 trigger
 134     TH1F* fSPDGFOEfficiency;         // SPD FASTOR efficiency - is applied in SPDFiredChips. Histogram contains efficiency as function of chip number (bin 1..400: first layer; 401..1200: second layer)
 135
 136     Float_t fV0TimeOffset;          // time offset applied to the times read from the V0 (in ns)
 137     Float_t fV0AdcThr;              // thresholds applied on V0 ADC data
 138     Float_t fV0HwAdcThr;            // online V0 trigger - thresholds applied on ADC data
 139     Float_t fV0HwWinLow;            // online V0 trigger - lower edge of time window
 140     Float_t fV0HwWinHigh;           // online V0 trigger - upper edge of time window
 141
 142     Float_t fZDCCutRefSum;          // ZDC time cut configuration
 143     Float_t fZDCCutRefDelta;        // ZDC time cut configuration
 144     Float_t fZDCCutSigmaSum;        // ZDC time cut configuration
 145     Float_t fZDCCutSigmaDelta;      // ZDC time cut configuration
 146
 147     Float_t fZDCCutRefSumCorr;      // Corrected ZDC time cut configuration
 148     Float_t fZDCCutRefDeltaCorr;    // Corrected ZDC time cut configuration
 149     Float_t fZDCCutSigmaSumCorr;    // Corrected ZDC time cut configuration
 150     Float_t fZDCCutSigmaDeltaCorr;  // Corrected ZDC time cut configuration
 151
 152     Float_t fZDCCutZNATimeCorrMin;  // Corrected ZNA time cut configuration
 153     Float_t fZDCCutZNATimeCorrMax;  // Corrected ZNA time cut configuration
 154     Float_t fZDCCutZNCTimeCorrMin;  // Corrected ZNA time cut configuration
 155     Float_t fZDCCutZNCTimeCorrMax;  // Corrected ZNA time cut configuration
 156
 157     Float_t fASPDCvsTCut; // constant for the linear cut in SPD clusters vs tracklets
 158     Float_t fBSPDCvsTCut; // slope for the linear cut in SPD  clusters vs tracklets
 159
 160     // Variables for the TRD triggers
 161     Float_t fTRDptHSE;          // pt threshold for HSE trigger
 162     UChar_t fTRDpidHSE;         // PID threshold for HSE trigger
 163     Float_t fTRDptHQU;          // pt threshold for HQU trigger
 164     UChar_t fTRDpidHQU;         // PID threshold for HQU trigger
 165     Float_t fTRDptHEE;          // pt threshold for HEE trigger
 166     UChar_t fTRDpidHEE;         // PID threshold for HEE trigger
 167     UChar_t fTRDminSectorHEE;   // min sector for HEE trigger
 168     UChar_t fTRDmaxSectorHEE;   // max sector for HEE trigger
 169     Float_t fTRDptHJT;          // pt threshold for HJT trigger
 170     UChar_t fTRDnHJT;           // no of track threshold for HJT trigger
 171     Bool_t  fUseTRDMatchCut;    // flag to decide whether to use or not the matching cut
 172
 173     Bool_t  fDoFMD;                 // If false, skips the FMD (physics selection runs much faster)
 174     Float_t fFMDLowCut;             //
 175     Float_t fFMDHitCut;             //
 176
 177     TH2F* fHistBitsSPD;        // offline trigger bits (calculated from clusters) vs hardware trigger bits
 178     TH1F* fHistFiredBitsSPD;   // fired hardware bits
 179     TH2F* fHistSPDClsVsTrk;    // histogram of clusters vs tracklet BG cut
 180     TH1F* fHistV0A;            // histograms that histogram the criterion the cut is applied on: bb triggers
 181     TH1F* fHistV0C;            // histograms that histogram the criterion the cut is applied on: bb triggers
 182     TH1F* fHistZDC;            //histograms that histogram the criterion the cut is applied on: fired bits (6 bins)
 183     TH1F* fHistTDCZDC;         // histograms that histogram the criterion the cut is applied on: TDC bits (32 bins)
 184     TH2F* fHistTimeZDC;        // histograms that histogram the criterion the cut is applied on: ZDC TDC timing
 185     TH2F* fHistTimeCorrZDC;    // histograms that histogram the criterion the cut is applied on: ZDC Corrected TDC timing
 186     TH1F* fHistFMDA;           // histograms that histogram the criterion the cut is applied on: number of hit combination above threshold
 187     TH1F* fHistFMDC;           // histograms that histogram the criterion the cut is applied on: number of hit combination above threshold
 188     TH1F* fHistFMDSingle;      // histograms that histogram the criterion the cut is applied on: single mult value (more than one entry per event)
 189     TH1F* fHistFMDSum;         // histograms that histogram the criterion the cut is applied on: summed mult value (more than one entry per event)
 190     TH1F* fHistT0;             // histograms that histogram the criterion the cut is applied on: bb triggers
 191     TMap* fTriggerClasses;    // counts the active trigger classes (uses the full string)
 192
 193     Bool_t fMC;              // flag if MC is analyzed
 194     AliESDtrackCuts* fEsdTrackCuts;  //Track Cuts to select ESD tracks
 195
 196     Bool_t fTPCOnly;         // flag to set whether TPC only tracks have to be used for the offline trigger
 197
 198     ClassDef(AliTriggerAnalysis, 20)
 199
 200   private:
 201     AliTriggerAnalysis(const AliTriggerAnalysis&);
 202     AliTriggerAnalysis& operator=(const AliTriggerAnalysis&);
 203 };
 204
 205 #endif