The corrected ZDC timing cuts as members of the class.
authorcvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 19 Apr 2011 12:28:26 +0000 (12:28 +0000)
committercvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 19 Apr 2011 12:28:26 +0000 (12:28 +0000)
ANALYSIS/AliTriggerAnalysis.cxx
ANALYSIS/AliTriggerAnalysis.h

index 05145d7bf0fad24c6eab2e46d79e0dc1a3c5108c..ebfdf507e8ac36b127319fd44576c009e6a7b8fd 100644 (file)
@@ -58,6 +58,10 @@ AliTriggerAnalysis::AliTriggerAnalysis() :
   fZDCCutRefDelta(-2.1),
   fZDCCutSigmaSum(3.25),
   fZDCCutSigmaDelta(2.25),
+  fZDCCutRefSumCorr(-65.5),
+  fZDCCutRefDeltaCorr(-2.1),
+  fZDCCutSigmaSumCorr(6.0),
+  fZDCCutSigmaDeltaCorr(1.2),
   fDoFMD(kTRUE),
   fFMDLowCut(0.2),
   fFMDHitCut(0.5),
@@ -1191,12 +1195,6 @@ Bool_t AliTriggerAnalysis::ZDCTimeTrigger(const AliESDEvent *aEsd, Bool_t fillHi
      zdcAccept = (znaFired | zncFired);
   }
   else {
-
-    // Cuts for the new corrected TDC values
-    const Float_t refSumCorr = -65.5;
-    const Float_t sigmaSumCorr = 6.0;
-    const Float_t sigmaDeltaCorr = 1.2;
-
     for(Int_t i = 0; i < 4; ++i) {
       if (esdZDC->GetZDCTDCData(10,i) != 0) {
        Float_t tdcC = 0.025*(esdZDC->GetZDCTDCData(10,i)-esdZDC->GetZDCTDCData(14,i)); 
@@ -1211,8 +1209,8 @@ Bool_t AliTriggerAnalysis::ZDCTimeTrigger(const AliESDEvent *aEsd, Bool_t fillHi
              fHistTimeCorrZDC->Fill(tdcCcorr-tdcAcorr,tdcCcorr+tdcAcorr);
            }
            if (esdZDC->TestBit(AliESDZDC::kCorrectedTDCFilled)) {
-             if (((tdcCcorr-tdcAcorr-fZDCCutRefDelta)*(tdcCcorr-tdcAcorr-fZDCCutRefDelta)/(sigmaDeltaCorr*sigmaDeltaCorr) +
-                  (tdcCcorr+tdcAcorr-refSumCorr)*(tdcCcorr+tdcAcorr-refSumCorr)/(sigmaSumCorr*sigmaSumCorr))< 1.0)
+             if (((tdcCcorr-tdcAcorr-fZDCCutRefDeltaCorr)*(tdcCcorr-tdcAcorr-fZDCCutRefDeltaCorr)/(fZDCCutSigmaDeltaCorr*fZDCCutSigmaDeltaCorr) +
+                  (tdcCcorr+tdcAcorr-fZDCCutRefSumCorr)*(tdcCcorr+tdcAcorr-fZDCCutRefSumCorr)/(fZDCCutSigmaSumCorr*fZDCCutSigmaSumCorr))< 1.0)
                zdcAccept = kTRUE;
            }
            else {
index b0930c7a7ad9588fd86156735bb6aa68c747e465..7dd19212d101201027176794a76358078ecee481 100644 (file)
@@ -76,6 +76,7 @@ class AliTriggerAnalysis : public TObject
     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; }
 
     Int_t GetSPDGFOThreshhold() const { return fSPDGFOThreshold; }
     Float_t GetV0TimeOffset() const { return fV0TimeOffset; }
@@ -118,6 +119,11 @@ 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
+
     Bool_t  fDoFMD;                 // If false, skips the FMD (physics selection runs much faster)
     Float_t fFMDLowCut;                    // 
     Float_t fFMDHitCut;                    // 
@@ -140,7 +146,7 @@ class AliTriggerAnalysis : public TObject
     Bool_t fMC;              // flag if MC is analyzed
     AliESDtrackCuts* fEsdTrackCuts;  //Track Cuts to select ESD tracks
 
-    ClassDef(AliTriggerAnalysis, 14)
+    ClassDef(AliTriggerAnalysis, 15)
     
   private:
     AliTriggerAnalysis(const AliTriggerAnalysis&);