Fix in spread function integral calculation (parameters order)

[u/mrichter/AliRoot.git] / ANALYSIS / AliTriggerAnalysis.cxx

diff --git a/ANALYSIS/AliTriggerAnalysis.cxx b/ANALYSIS/AliTriggerAnalysis.cxx
index 68bcf7f5b9ccce7cf12e0f4aa57ef2e7f6504911..218e31cf5a62873f9674609812bb8e8ea251b78a 100644 (file)
--- a/ANALYSIS/AliTriggerAnalysis.cxx
+++ b/ANALYSIS/AliTriggerAnalysis.cxx
@@ -249,6 +249,10 @@ const char* AliTriggerAnalysis::GetTriggerName(Trigger trigger)
     case kZDC : str = "ZDC"; break;
     case kZDCA : str = "ZDC A"; break;
     case kZDCC : str = "ZDC C"; break;
+    case kZNA : str = "ZN A"; break;
+    case kZNC : str = "ZN C"; break;
+    case kZNABG : str = "ZN A BG"; break;
+    case kZNCBG : str = "ZN C BG"; break;
     case kFMDA : str = "FMD A"; break;
     case kFMDC : str = "FMD C"; break;
     case kFPANY : str = "SPD GFO | V0 | ZDC | FMD"; break;
@@ -471,6 +475,38 @@ Int_t AliTriggerAnalysis::EvaluateTrigger(const AliESDEvent* aEsd, Trigger trigg
         decision = 1;
       break;
     }
+    case kZNA:
+    {
+      if (!offline)
+        AliFatal(Form("Online trigger not available for trigger %d", triggerNoFlags));
+      if (ZDCTDCTrigger(aEsd,kASide,kTRUE,kFALSE,kFALSE))
+       decision = 1;
+      break;
+    }
+    case kZNC:
+    {
+      if (!offline)
+        AliFatal(Form("Online trigger not available for trigger %d", triggerNoFlags));
+      if (ZDCTDCTrigger(aEsd,kCSide,kTRUE,kFALSE,kFALSE))
+       decision = 1;
+      break;
+    }
+    case kZNABG:
+    {
+      if (!offline)
+        AliFatal(Form("Online trigger not available for trigger %d", triggerNoFlags));
+      if (ZDCTimeBGTrigger(aEsd,kASide))
+       decision = 1;
+      break;
+    }
+    case kZNCBG:
+    {
+      if (!offline)
+        AliFatal(Form("Online trigger not available for trigger %d", triggerNoFlags));
+      if (ZDCTimeBGTrigger(aEsd,kCSide))
+       decision = 1;
+      break;
+    }
     case kFMDA:
     {
       if (!offline)
@@ -683,6 +719,30 @@ Bool_t AliTriggerAnalysis::IsOfflineTriggerFired(const AliESDEvent* aEsd, Trigge
         decision = kTRUE;
       break;
     }
+    case kZNA:
+    {
+      if (ZDCTDCTrigger(aEsd,kASide,kTRUE,kFALSE,kFALSE))
+       decision = kTRUE;
+      break;
+    }
+    case kZNC:
+    {
+      if (ZDCTDCTrigger(aEsd,kCSide,kTRUE,kFALSE,kFALSE))
+       decision = kTRUE;
+      break;
+    }
+    case kZNABG:
+    {
+      if (ZDCTimeBGTrigger(aEsd,kASide))
+       decision = kTRUE;
+      break;
+    }
+    case kZNCBG:
+    {
+      if (ZDCTimeBGTrigger(aEsd,kCSide))
+       decision = kTRUE;
+      break;
+    }
     case kFMDA:
     {
       if (FMDTrigger(aEsd, kASide))
@@ -1372,6 +1432,52 @@ Bool_t AliTriggerAnalysis::ZDCTimeTrigger(const AliESDEvent *aEsd, Bool_t fillHi
   return zdcAccept;
 }
 
+Bool_t AliTriggerAnalysis::ZDCTimeBGTrigger(const AliESDEvent *aEsd, AliceSide side) const
+{
+  // This method implements a selection
+  // based on the timing in of zdcN
+  // It can be used in order to flag background
+  // ** So far only implemented for the 2012 pA run **
+
+  if(fMC) return kFALSE;
+
+  AliESDZDC *zdcData = aEsd->GetESDZDC();
+  Bool_t zna = kFALSE;
+  Bool_t znc = kFALSE;
+
+  Float_t tdcC=999, tdcCcorr=999, tdcA=999, tdcAcorr=999;
+  for(Int_t i = 0; i < 4; ++i) {
+    if (zdcData->GetZDCTDCData(10,i) != 0) {
+      znc = kTRUE;
+      tdcC = 0.025*(zdcData->GetZDCTDCData(10,i)-zdcData->GetZDCTDCData(14,i));
+      tdcCcorr = zdcData->GetZDCTDCCorrected(10,i);
+    }
+  }
+  for(Int_t j = 0; j < 4; ++j) {
+    if (zdcData->GetZDCTDCData(12,j) != 0) {
+      zna = kTRUE;
+      tdcA = 0.025*(zdcData->GetZDCTDCData(12,j)-zdcData->GetZDCTDCData(14,j));
+      tdcAcorr = zdcData->GetZDCTDCCorrected(12,j);
+    }
+  }
+
+  const Int_t runNumber = aEsd->GetRunNumber();
+  if(runNumber<188124 || (runNumber>188374 && runNumber<194713)){ // FIXME: end of pA-run is not known
+    AliError(Form(" ZN BG time cut not implemented for run %d",runNumber));
+    return kFALSE;
+  }
+
+  Bool_t znabg = (zna && (TMath::Abs(tdcAcorr)>2.0));
+  Bool_t zncbg = (znc && (TMath::Abs(tdcCcorr)>5.0));
+
+  // Printf("Checking ZN background (time) for run %d, A = %d, time=%2.2f, C = %d, time=%2.2f",runNumber,(Int_t)zna,tdcAcorr,(Int_t)znc,tdcCcorr);
+  // Printf("Checking ZN background (time) for run %d, A-BG = %d, C-BG = %d",runNumber,(Int_t)znabg,(Int_t)zncbg);
+
+  if (side == kASide) return znabg;
+  if (side == kCSide) return zncbg;
+  return kFALSE;
+}
+
 Bool_t AliTriggerAnalysis::ZDCTrigger(const AliESDEvent* aEsd, AliceSide side) const
 {
   // Returns if ZDC triggered
@@ -1682,13 +1788,13 @@ AliTriggerAnalysis::T0Decision AliTriggerAnalysis::T0Trigger(const AliESDEvent*
   Float_t  tvdc[5] ;
   for (Int_t ii=0; ii<5; ii++)
     tvdc[ii] = esdT0->GetTVDC(ii);
-  Int_t trig=esdT0->GetT0Trig();
+  //  Int_t trig=esdT0->GetT0Trig();
   //  cout<<" T0 trig "<<trig<<endl;   
 
   if(fillHists) fHistT0->Fill(tvdc[0]);
     
   if (online) {
-    if(trig&1) return  kT0BB;
+    if(aEsd->GetHeader()->GetFiredTriggerInputs().Contains("0TVX") ) return kT0BB;
   }
   else {