Merge branch 'master', remote branch 'origin' into TPCdev
[u/mrichter/AliRoot.git] / ANALYSIS / AliTriggerAnalysis.h
index 7b8b191..85dc2d2 100644 (file)
@@ -26,16 +26,20 @@ class AliTriggerAnalysis : public TObject
 {
   public:
     enum Trigger { kAcceptAll = 1, kMB1 = 2, kMB2, kMB3, kSPDGFO, kSPDGFOBits, kV0A, kV0C, kV0OR, kV0AND, 
-      kV0ABG, kV0CBG, kZDC, kZDCA, kZDCC, kFMDA, kFMDC, kFPANY, kNSD1, kMB1Prime, 
-      kSPDGFOL0, kSPDGFOL1, kZDCTDCA, kZDCTDCC, kZDCTime, kCTPV0A, kCTPV0C, kTPCLaserWarmUp, 
-      kStartOfFlags = 0x0100, kOfflineFlag = 0x8000, kOneParticle = 0x10000, kOneTrack = 0x20000}; // MB1, MB2, MB3 definition from ALICE-INT-2005-025
+                  kV0ABG, kV0CBG, kZDC, kZDCA, kZDCC, kZNA, kZNC, kZNABG, kZNCBG, kFMDA, kFMDC, kFPANY, kNSD1, kMB1Prime, 
+                  kSPDGFOL0, kSPDGFOL1, kZDCTDCA, kZDCTDCC, kZDCTime, kCTPV0A, kCTPV0C, kTPCLaserWarmUp, kSPDClsVsTrkBG,
+                  kCentral,kSemiCentral, kT0, kT0BG, kT0Pileup, kEMCAL,kTPCHVdip,
+                  kTRDHCO, kTRDHJT, kTRDHSE, kTRDHQU, kTRDHEE,
+                  kIncompleteEvent,
+                  kStartOfFlags = 0x0100, kOfflineFlag = 0x8000, kOneParticle = 0x10000, kOneTrack = 0x20000}; // MB1, MB2, MB3 definition from ALICE-INT-2005-025
     enum AliceSide { kASide = 1, kCSide, kCentralBarrel };
     enum V0Decision { kV0Invalid = -1, kV0Empty = 0, kV0BB, kV0BG, kV0Fake };
+    enum T0Decision { kT0Invalid = -1, kT0Empty = 0, kT0BB, kT0DecBG, kT0DecPileup };
     
     AliTriggerAnalysis();
     virtual ~AliTriggerAnalysis();
     
-    void EnableHistograms();
+    void EnableHistograms(Bool_t isLowFlux = kFALSE);
     void SetAnalyzeMC(Bool_t flag = kTRUE) { fMC = flag; }
     
     Bool_t IsTriggerFired(const AliESDEvent* aEsd, Trigger trigger);
@@ -55,27 +59,48 @@ class AliTriggerAnalysis : public TObject
     // some "raw" trigger functions
     Int_t SPDFiredChips(const AliESDEvent* aEsd, Int_t origin, Bool_t fillHists = kFALSE, Int_t layer = 0);
     Bool_t SPDGFOTrigger(const AliESDEvent* aEsd, Int_t origin);
+    Bool_t IsSPDClusterVsTrackletBG(const AliESDEvent* esd, Bool_t fillHists = kFALSE);
     V0Decision V0Trigger(const AliESDEvent* aEsd, AliceSide side, Bool_t online, Bool_t fillHists = kFALSE);
+    T0Decision T0Trigger(const AliESDEvent* aEsd, Bool_t online, Bool_t fillHists = kFALSE);
     Bool_t ZDCTrigger   (const AliESDEvent* aEsd, AliceSide side) const;
     Bool_t ZDCTDCTrigger(const AliESDEvent* aEsd, AliceSide side, Bool_t useZN=kTRUE, Bool_t useZP=kFALSE, Bool_t fillHists=kFALSE) const;
     Bool_t ZDCTimeTrigger(const AliESDEvent *aEsd, Bool_t fillHists=kFALSE) const;
+    Bool_t ZDCTimeBGTrigger(const AliESDEvent *aEsd, AliceSide side) const;
     Bool_t FMDTrigger(const AliESDEvent* aEsd, AliceSide side);
     Int_t SSDClusters(const AliESDEvent* aEsd);
+    Bool_t EMCALCellsTrigger(const AliESDEvent *aEsd);
     static const char* GetTriggerName(Trigger trigger);
     
     Bool_t IsLaserWarmUpTPCEvent(const AliESDEvent* esd);
+    Bool_t IsHVdipTPCEvent(const AliESDEvent* esd);
+    Bool_t TRDTrigger(const AliESDEvent* esd, Trigger trigger);
+
+    Bool_t IsIncompleteEvent(const AliESDEvent* esd);
     
     void FillHistograms(const AliESDEvent* aEsd);
     void FillTriggerClasses(const AliESDEvent* aEsd);
     
     void SetSPDGFOThreshhold(Int_t t) { fSPDGFOThreshold = t; }
     void SetSPDGFOEfficiency(TH1F* hist) { fSPDGFOEfficiency = hist; }
+    void SetSPDClustersVsTrackletsParameters(Float_t a, Float_t b) { fASPDCvsTCut = a; fBSPDCvsTCut =b;}
     void SetV0TimeOffset(Float_t offset) { fV0TimeOffset = offset; }
     void SetV0AdcThr(Float_t thr) { fV0AdcThr = thr; }
     void SetV0HwPars(Float_t thr, Float_t winLow, Float_t winHigh) { fV0HwAdcThr = thr; fV0HwWinLow = winLow; fV0HwWinHigh = winHigh; }
     void SetFMDThreshold(Float_t low, Float_t hit) { fFMDLowCut = low; fFMDHitCut = hit; }
     void SetDoFMD(Bool_t flag = kTRUE) {fDoFMD = flag;}
     void SetZDCCutParams(Float_t refSum, Float_t refDelta, Float_t sigmaSum, Float_t sigmaDelta) { fZDCCutRefSum = refSum; fZDCCutRefDelta = refDelta; fZDCCutSigmaSum = sigmaSum; fZDCCutSigmaDelta = sigmaDelta; }
+    void SetCorrZDCCutParams(Float_t refSum, Float_t refDelta, Float_t sigmaSum, Float_t sigmaDelta) { fZDCCutRefSumCorr = refSum; fZDCCutRefDeltaCorr = refDelta; fZDCCutSigmaSumCorr = sigmaSum; fZDCCutSigmaDeltaCorr = sigmaDelta; }
+    void SetZNCorrCutParams(Float_t znaTimeCorrMin, Float_t znaTimeCorrMax, Float_t zncTimeCorrMin, Float_t zncTimeCorrMax)
+    { fZDCCutZNATimeCorrMin = znaTimeCorrMin; fZDCCutZNATimeCorrMax = znaTimeCorrMax; 
+      fZDCCutZNCTimeCorrMin = zncTimeCorrMin; fZDCCutZNCTimeCorrMax = zncTimeCorrMax; }
+
+    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) {
+      fTRDptHSE = ptHSE; fTRDpidHSE = pidHSE;
+      fTRDptHQU = ptHQU; fTRDpidHQU = pidHQU;
+      fTRDptHEE = ptHEE; fTRDpidHEE = pidHEE;
+      fTRDminSectorHEE = minSectorHEE; fTRDmaxSectorHEE = maxSectorHEE;
+      fTRDptHJT = ptHJT; fTRDnHJT = nHJT;
+    }
 
     Int_t GetSPDGFOThreshhold() const { return fSPDGFOThreshold; }
     Float_t GetV0TimeOffset() const { return fV0TimeOffset; }
@@ -92,7 +117,8 @@ class AliTriggerAnalysis : public TObject
     void SetESDTrackCuts(AliESDtrackCuts* cuts) { fEsdTrackCuts = cuts;}
     AliESDtrackCuts* GetESDTrackCuts() const  {return fEsdTrackCuts;}
 
-    static const UInt_t GetActiveBit(UInt_t mask) ;
+    void SetTPCOnly(Bool_t bTPCOnly) {fTPCOnly = bTPCOnly;}
+    Bool_t GetTPCOnly() const {return fTPCOnly;}
 
   protected:
     Bool_t IsL0InputFired(const AliESDEvent* aEsd, UInt_t input) const;
@@ -119,12 +145,38 @@ class AliTriggerAnalysis : public TObject
     Float_t fZDCCutSigmaSum;        // ZDC time cut configuration
     Float_t fZDCCutSigmaDelta;      // ZDC time cut configuration
 
+    Float_t fZDCCutRefSumCorr;      // Corrected ZDC time cut configuration
+    Float_t fZDCCutRefDeltaCorr;    // Corrected ZDC time cut configuration
+    Float_t fZDCCutSigmaSumCorr;    // Corrected ZDC time cut configuration
+    Float_t fZDCCutSigmaDeltaCorr;  // Corrected ZDC time cut configuration
+
+    Float_t fZDCCutZNATimeCorrMin;  // Corrected ZNA time cut configuration
+    Float_t fZDCCutZNATimeCorrMax;  // Corrected ZNA time cut configuration
+    Float_t fZDCCutZNCTimeCorrMin;  // Corrected ZNA time cut configuration
+    Float_t fZDCCutZNCTimeCorrMax;  // Corrected ZNA time cut configuration
+
+    Float_t fASPDCvsTCut; // constant for the linear cut in SPD clusters vs tracklets
+    Float_t fBSPDCvsTCut; // slope for the linear cut in SPD  clusters vs tracklets
+
+    // Variables for the TRD triggers
+    Float_t fTRDptHSE;         // pt threshold for HSE trigger
+    UChar_t fTRDpidHSE;                // PID threshold for HSE trigger
+    Float_t fTRDptHQU;         // pt threshold for HQU trigger
+    UChar_t fTRDpidHQU;                // PID threshold for HQU trigger
+    Float_t fTRDptHEE;         // pt threshold for HEE trigger
+    UChar_t fTRDpidHEE;                // PID threshold for HEE trigger
+    UChar_t fTRDminSectorHEE;  // min sector for HEE trigger
+    UChar_t fTRDmaxSectorHEE;  // max sector for HEE trigger
+    Float_t fTRDptHJT;         // pt threshold for HJT trigger
+    UChar_t fTRDnHJT;          // no of track threshold for HJT trigger
+
     Bool_t  fDoFMD;                 // If false, skips the FMD (physics selection runs much faster)
     Float_t fFMDLowCut;                    // 
     Float_t fFMDHitCut;                    // 
     
     TH2F* fHistBitsSPD;        // offline trigger bits (calculated from clusters) vs hardware trigger bits
     TH1F* fHistFiredBitsSPD;   // fired hardware bits
+    TH2F* fHistSPDClsVsTrk;    // histogram of clusters vs tracklet BG cut
     TH1F* fHistV0A;            // histograms that histogram the criterion the cut is applied on: bb triggers
     TH1F* fHistV0C;            // histograms that histogram the criterion the cut is applied on: bb triggers
     TH1F* fHistZDC;            //histograms that histogram the criterion the cut is applied on: fired bits (6 bins)
@@ -135,13 +187,15 @@ class AliTriggerAnalysis : public TObject
     TH1F* fHistFMDC;           // histograms that histogram the criterion the cut is applied on: number of hit combination above threshold
     TH1F* fHistFMDSingle;      // histograms that histogram the criterion the cut is applied on: single mult value (more than one entry per event)
     TH1F* fHistFMDSum;         // histograms that histogram the criterion the cut is applied on: summed mult value (more than one entry per event)
-    
+    TH1F* fHistT0;             // histograms that histogram the criterion the cut is applied on: bb triggers
     TMap* fTriggerClasses;    // counts the active trigger classes (uses the full string)
     
     Bool_t fMC;              // flag if MC is analyzed
     AliESDtrackCuts* fEsdTrackCuts;  //Track Cuts to select ESD tracks
 
-    ClassDef(AliTriggerAnalysis, 13)
+    Bool_t fTPCOnly;         // flag to set whether TPC only tracks have to be used for the offline trigger 
+
+    ClassDef(AliTriggerAnalysis, 22)
     
   private:
     AliTriggerAnalysis(const AliTriggerAnalysis&);