Introduce the ZN background check for the pA run (INT5 trigger)
authorzconesa <zconesa@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 18 Sep 2012 14:53:11 +0000 (14:53 +0000)
committerzconesa <zconesa@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 18 Sep 2012 14:53:11 +0000 (14:53 +0000)
ANALYSIS/AliTriggerAnalysis.cxx
ANALYSIS/AliTriggerAnalysis.h
ANALYSIS/macros/BrowseAndFillPhysicsSelectionOADB.C
OADB/COMMON/PHYSICSSELECTION/data/physicsSelection.root

index 68bcf7f..efdcab3 100644 (file)
@@ -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,49 @@ 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
+
+  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){
+    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
index 5ad04d3..221f17f 100644 (file)
@@ -26,7 +26,7 @@ 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, 
+                  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,
                   kStartOfFlags = 0x0100, kOfflineFlag = 0x8000, kOneParticle = 0x10000, kOneTrack = 0x20000}; // MB1, MB2, MB3 definition from ALICE-INT-2005-025
@@ -63,6 +63,7 @@ class AliTriggerAnalysis : public TObject
     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);
index 2b8c385..73c39e3 100644 (file)
@@ -1228,7 +1228,7 @@ void BrowseAndFillPhysicsSelectionOADB(Bool_t fill = kFALSE) {
   oadbLHC12g->AddBGTriggerClass          ( AliVEvent::kCINT5,"+CINT5-C-NOPF-ALLNOTRD","C",triggerCount);
   oadbLHC12g->AddBGTriggerClass          ( AliVEvent::kCINT5,"+CINT5-E-NOPF-ALLNOTRD","E",triggerCount);  
   oadbLHC12g->SetHardwareTrigger         ( triggerCount,"V0A || V0C");
-  oadbLHC12g->SetOfflineTrigger          ( triggerCount,"(V0A || V0C) && !V0ABG && !V0CBG  && !TPCLaserWarmUp");
+  oadbLHC12g->SetOfflineTrigger          ( triggerCount,"(V0A || V0C) && !V0ABG && !V0CBG && !ZNABG && !ZNCBG && !TPCLaserWarmUp");
 
   triggerCount++;
   oadbLHC12g->AddCollisionTriggerClass   ( AliVEvent::kZED,"+C1ZED-B-NOPF-ALLNOTRD", "B",  triggerCount);
index 74df31e..9ee2811 100644 (file)
Binary files a/OADB/COMMON/PHYSICSSELECTION/data/physicsSelection.root and b/OADB/COMMON/PHYSICSSELECTION/data/physicsSelection.root differ