]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ANALYSIS/AliTriggerAnalysis.cxx
Fixes in the command to submit merging
[u/mrichter/AliRoot.git] / ANALYSIS / AliTriggerAnalysis.cxx
index 16fa060b827adbb60d2456704f717ccf24b20a5f..7e5b17c55ac29fd889e6a8703d71976274772ed0 100644 (file)
@@ -41,6 +41,8 @@
 #include <AliESDVZERO.h>
 #include <AliESDZDC.h>
 #include <AliESDFMD.h>
+#include <AliESDVertex.h>
+#include <AliESDtrackCuts.h>
 
 ClassImp(AliTriggerAnalysis)
 
@@ -52,6 +54,7 @@ AliTriggerAnalysis::AliTriggerAnalysis() :
   fV0HwAdcThr(2.5),
   fV0HwWinLow(61.5),
   fV0HwWinHigh(86.5),
+  fDoFMD(kTRUE),
   fFMDLowCut(0.2),
   fFMDHitCut(0.5),
   fHistBitsSPD(0),
@@ -59,12 +62,15 @@ AliTriggerAnalysis::AliTriggerAnalysis() :
   fHistV0A(0),       
   fHistV0C(0),
   fHistZDC(0),    
+  fHistTDCZDC(0),    
+  fHistTimeZDC(0),    
   fHistFMDA(0),    
   fHistFMDC(0),   
   fHistFMDSingle(0),
   fHistFMDSum(0),
   fTriggerClasses(0),
-  fMC(kFALSE)
+  fMC(kFALSE),
+  fEsdTrackCuts(0)
 {
   // constructor
 }
@@ -102,6 +108,16 @@ AliTriggerAnalysis::~AliTriggerAnalysis()
     delete fHistZDC;
     fHistZDC = 0;
   }
+  if (fHistTDCZDC)
+  {
+    delete fHistTDCZDC;
+    fHistTDCZDC = 0;
+  }
+  if (fHistTimeZDC)
+  {
+    delete fHistTimeZDC;
+    fHistTimeZDC = 0;
+  }
 
   if (fHistFMDA)
   {
@@ -133,6 +149,11 @@ AliTriggerAnalysis::~AliTriggerAnalysis()
     delete fTriggerClasses;
     fTriggerClasses = 0;
   }
+
+  if (fEsdTrackCuts){
+    delete fEsdTrackCuts;
+    fEsdTrackCuts =0;
+  }
 }
 
 void AliTriggerAnalysis::EnableHistograms()
@@ -145,9 +166,11 @@ void AliTriggerAnalysis::EnableHistograms()
   
   fHistBitsSPD = new TH2F("fHistBitsSPD", "SPD GFO;number of fired chips (offline);number of fired chips (hardware)", 1202, -1.5, 1200.5, 1202, -1.5, 1200.5);
   fHistFiredBitsSPD = new TH1F("fHistFiredBitsSPD", "SPD GFO Hardware;chip number;events", 1200, -0.5, 1199.5);
-  fHistV0A = new TH1F("fHistV0A", "V0A;leading time (ns);events", 200, 0, 100);
-  fHistV0C = new TH1F("fHistV0C", "V0C;leading time (ns);events", 200, 0, 100);
+  fHistV0A = new TH1F("fHistV0A", "V0A;leading time (ns);events", 400, -100, 100);
+  fHistV0C = new TH1F("fHistV0C", "V0C;leading time (ns);events", 400, -100, 100);
   fHistZDC = new TH1F("fHistZDC", "ZDC;trigger bits;events", 8, -1.5, 6.5);
+  fHistTDCZDC = new TH1F("fHistTDCZDC", "ZDC;TDC bits;events", 32, -0.5, 32-0.5);
+  fHistTimeZDC = new TH2F("fHistTimeZDC", "ZDC;TDC timing A+C vs C-A; events", 120,-30,30,120,-600,-540);
   
   // TODO check limits
   fHistFMDA = new TH1F("fHistFMDA", "FMDA;combinations above threshold;events", 102, -1.5, 100.5);
@@ -199,6 +222,12 @@ const char* AliTriggerAnalysis::GetTriggerName(Trigger trigger)
   if (trigger & kOfflineFlag)
     str += " OFFLINE";  
   
+  if (trigger & kOneParticle)
+    str += " OneParticle";  
+
+  if (trigger & kOneTrack)
+    str += " OneTrack";  
+
   return str;
 }
 
@@ -283,119 +312,120 @@ Bool_t AliTriggerAnalysis::IsOfflineTriggerFired(const AliESDEvent* aEsd, Trigge
 
   UInt_t triggerNoFlags = (UInt_t) trigger % (UInt_t) kStartOfFlags;
   
+  Bool_t decision = kFALSE;
   switch (triggerNoFlags)
   {
     case kAcceptAll:
     {
-      return kTRUE;
+      decision = kTRUE;
       break;
     }
     case kMB1:
     {
       if (SPDGFOTrigger(aEsd, 0) || V0Trigger(aEsd, kASide, kFALSE) == kV0BB || V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kMB2:
     {
       if (SPDGFOTrigger(aEsd, 0) && (V0Trigger(aEsd, kASide, kFALSE) == kV0BB || V0Trigger(aEsd, kCSide, kFALSE) == kV0BB))
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kMB3:
     {
       if (SPDGFOTrigger(aEsd, 0) && V0Trigger(aEsd, kASide, kFALSE) == kV0BB && V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kSPDGFO:
     {
       if (SPDGFOTrigger(aEsd, 0))
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kSPDGFOBits:
     {
       if (SPDGFOTrigger(aEsd, 1))
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kV0A:
     {
       if (V0Trigger(aEsd, kASide, kFALSE) == kV0BB)
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kV0C:
     {
       if (V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kV0OR:
     {
       if (V0Trigger(aEsd, kASide, kFALSE) == kV0BB || V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kV0AND:
     {
       if (V0Trigger(aEsd, kASide, kFALSE) == kV0BB && V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kV0ABG:
     {
       if (V0Trigger(aEsd, kASide, kFALSE) == kV0BG)
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kV0CBG:
     {
       if (V0Trigger(aEsd, kCSide, kFALSE) == kV0BG)
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kZDC:
     {
       if (ZDCTrigger(aEsd, kASide) || ZDCTrigger(aEsd, kCentralBarrel) || ZDCTrigger(aEsd, kCSide))
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kZDCA:
     {
       if (ZDCTrigger(aEsd, kASide))
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kZDCC:
     {
       if (ZDCTrigger(aEsd, kCSide))
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kFMDA:
     {
       if (FMDTrigger(aEsd, kASide))
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kFMDC:
     {
       if (FMDTrigger(aEsd, kCSide))
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kFPANY:
     {
       if (SPDGFOTrigger(aEsd, 0) || V0Trigger(aEsd, kASide, kFALSE) == kV0BB || V0Trigger(aEsd, kCSide, kFALSE) == kV0BB || ZDCTrigger(aEsd, kASide) || ZDCTrigger(aEsd, kCentralBarrel) || ZDCTrigger(aEsd, kCSide) || FMDTrigger(aEsd, kASide) || FMDTrigger(aEsd, kCSide))
-        return kTRUE;
+        decision = kTRUE;
       break;
     }
     case kNSD1:
     {
       if (SPDFiredChips(aEsd, 0) >= 5 || (V0Trigger(aEsd, kASide, kFALSE) == kV0BB && V0Trigger(aEsd, kCSide, kFALSE) == kV0BB))
-        return kTRUE;
+        decision = kTRUE;
        break;
     }
     case kMB1Prime:
@@ -409,7 +439,7 @@ Bool_t AliTriggerAnalysis::IsOfflineTriggerFired(const AliESDEvent* aEsd, Trigge
         count++;
       
       if (count >= 2)
-        return kTRUE;
+        decision = kTRUE;
         
       break;
     }
@@ -419,9 +449,66 @@ Bool_t AliTriggerAnalysis::IsOfflineTriggerFired(const AliESDEvent* aEsd, Trigge
     }
   }
   
-  return kFALSE;
-}
+  // hadron-level requirement
+  if (decision && (trigger & kOneParticle))
+  {
+    decision = kFALSE;
+    
+    const AliESDVertex* vertex = aEsd->GetPrimaryVertexSPD();
+    const AliMultiplicity* mult = aEsd->GetMultiplicity();
+
+    if (mult && vertex && vertex->GetNContributors() > 0 && (!vertex->IsFromVertexerZ() || vertex->GetDispersion() < 0.02) && TMath::Abs(vertex->GetZv()) < 5.5) 
+    {
+      for (Int_t i=0; i<mult->GetNumberOfTracklets(); ++i)
+      {
+        if (TMath::Abs(mult->GetEta(i)) < 1)
+        {
+          decision = kTRUE;
+          break;
+        }
+      }
+    }
+  }
 
+  // hadron level definition for TPC tracks
+
+  if (decision && (trigger & kOneTrack))
+  {
+    decision = kFALSE;
+    const AliESDVertex* vertex =0x0;
+    vertex = aEsd->GetPrimaryVertexTracks();
+    if (!vertex || vertex->GetNContributors() <= 0)
+    {
+      vertex = aEsd->GetPrimaryVertexSPD();
+    }
+    Float_t ptmin, ptmax;
+    fEsdTrackCuts->GetPtRange(ptmin,ptmax);
+    AliDebug(3, Form("ptmin = %f, ptmax = %f\n",ptmin, ptmax));
+
+    if (vertex && vertex->GetNContributors() > 0 && (!vertex->IsFromVertexerZ() || vertex->GetDispersion() < 0.02) && TMath::Abs(vertex->GetZv()) < 10.) {
+      AliDebug(3,Form("Check on the vertex passed\n"));
+      for (Int_t i=0; i<aEsd->GetNumberOfTracks(); ++i){
+       if (fEsdTrackCuts->AcceptTrack(aEsd->GetTrack(i))){
+         AliDebug(2, Form("pt of track = %f --> check passed\n",aEsd->GetTrack(i)->Pt()));
+         decision = kTRUE;
+         break;
+        }
+      }
+    }
+    else{
+      AliDebug(4,Form("Check on the vertex not passed\n"));
+      for (Int_t i=0; i<aEsd->GetNumberOfTracks(); ++i){
+        if (fEsdTrackCuts->AcceptTrack(aEsd->GetTrack(i))){
+         AliDebug(4,Form("pt of track = %f --> check would be passed if the vertex was ok\n",aEsd->GetTrack(i)->Pt()));
+         break;
+       }
+      }
+    }
+    if (!decision) AliDebug(3,("Check for kOneTrack NOT passed\n"));
+  }
+
+  return decision;
+}
 
 Bool_t AliTriggerAnalysis::IsTriggerClassFired(const AliESDEvent* aEsd, const Char_t* tclass) const 
 {
@@ -504,7 +591,9 @@ void AliTriggerAnalysis::FillHistograms(const AliESDEvent* aEsd)
   
   V0Trigger(aEsd, kASide, kFALSE, kTRUE);
   V0Trigger(aEsd, kCSide, kFALSE, kTRUE);
-  
+  ZDCTDCTrigger(aEsd,kASide,kFALSE,kFALSE,kTRUE);
+  ZDCTimeTrigger(aEsd,kTRUE);
+
   AliESDZDC* zdcData = aEsd->GetESDZDC();
   if (zdcData)
   {
@@ -531,8 +620,10 @@ void AliTriggerAnalysis::FillHistograms(const AliESDEvent* aEsd)
     AliError("AliESDZDC not available");
   }
   
-  fHistFMDA->Fill(FMDHitCombinations(aEsd, kASide, kTRUE));
-  fHistFMDC->Fill(FMDHitCombinations(aEsd, kCSide, kTRUE));
+  if (fDoFMD) {
+    fHistFMDA->Fill(FMDHitCombinations(aEsd, kASide, kTRUE));
+    fHistFMDC->Fill(FMDHitCombinations(aEsd, kCSide, kTRUE));
+  }
 }
   
 void AliTriggerAnalysis::FillTriggerClasses(const AliESDEvent* aEsd)
@@ -546,8 +637,6 @@ void AliTriggerAnalysis::FillTriggerClasses(const AliESDEvent* aEsd)
     fTriggerClasses->Add(new TObjString(aEsd->GetFiredTriggerClasses().Data()), count);
   }
   count->SetVal(count->GetVal() + 1);
-  
-  // TODO add first and last orbit number here
 }
 
 Int_t AliTriggerAnalysis::SSDClusters(const AliESDEvent* aEsd)
@@ -559,12 +648,15 @@ Int_t AliTriggerAnalysis::SSDClusters(const AliESDEvent* aEsd)
 }
 
 
-Int_t AliTriggerAnalysis::SPDFiredChips(const AliESDEvent* aEsd, Int_t origin, Bool_t fillHists)
+Int_t AliTriggerAnalysis::SPDFiredChips(const AliESDEvent* aEsd, Int_t origin, Bool_t fillHists, Int_t layer)
 {
   // returns the number of fired chips in the SPD
   //
   // origin = 0 --> aEsd->GetMultiplicity()->GetNumberOfFiredChips() (filled from clusters)
   // origin = 1 --> aEsd->GetMultiplicity()->TestFastOrFiredChips() (from hardware bits)
+  // layer  = 0 --> both layers
+  // layer  = 1 --> inner
+  // layer  = 2 --> outer
   
   const AliMultiplicity* mult = aEsd->GetMultiplicity();
   if (!mult)
@@ -573,12 +665,23 @@ Int_t AliTriggerAnalysis::SPDFiredChips(const AliESDEvent* aEsd, Int_t origin, B
     return -1;
   }
   
-  if (origin == 0)
-    return mult->GetNumberOfFiredChips(0) + mult->GetNumberOfFiredChips(1);
+  if (origin == 0){
+    if (layer == 0) 
+      return mult->GetNumberOfFiredChips(0) + mult->GetNumberOfFiredChips(1);
+
+    return mult->GetNumberOfFiredChips(layer-1); 
+  }
     
   if (origin == 1)
   {
     Int_t nChips = 0;
+    Int_t firstChip = 0;
+    Int_t lastChip  = 1200;
+    if(layer == 1)
+      lastChip  = 400;
+    if(layer == 2)
+      firstChip = 400;
+
     for (Int_t i=0; i<1200; i++)
     {
       if (mult->TestFastOrFiredChips(i) == kTRUE)
@@ -623,13 +726,14 @@ AliTriggerAnalysis::V0Decision AliTriggerAnalysis::V0Trigger(const AliESDEvent*
   // argument 'online' is used as a switch between online and offline trigger algorithms
   //
   // Based on an algorithm by Cvetan Cheshkov
-  
+
   AliESDVZERO* esdV0 = aEsd->GetVZEROData();
   if (!esdV0)
   {
     AliError("AliESDVZERO not available");
     return kV0Invalid;
   }
+  AliDebug(2,Form("In V0Trigger: %f %f",esdV0->GetV0ATime(),esdV0->GetV0CTime()));
 
   Int_t begin = -1;
   Int_t end = -1;
@@ -647,6 +751,37 @@ AliTriggerAnalysis::V0Decision AliTriggerAnalysis::V0Trigger(const AliESDEvent*
   else
     return kV0Invalid;
     
+   if (esdV0->TestBit(AliESDVZERO::kDecisionFilled)) {
+    if (online) {
+      if (esdV0->TestBit(AliESDVZERO::kOnlineBitsFilled)) {
+       for (Int_t i = begin; i < end; ++i) {
+         if (esdV0->GetBBFlag(i)) return kV0BB;
+       }
+       for (Int_t i = begin; i < end; ++i) {
+         if (esdV0->GetBGFlag(i)) return kV0BG;
+       }
+       return kV0Empty;
+      }
+      else {
+       AliWarning("V0 online trigger analysis is not yet available!");
+       return kV0BB;
+      }
+    }
+    else {
+
+      if (fillHists) {
+       if (side == kASide && fHistV0A)
+         fHistV0A->Fill(esdV0->GetV0ATime());
+       if (side == kCSide && fHistV0C)
+         fHistV0C->Fill(esdV0->GetV0CTime());
+      }
+
+      if (side == kASide) return (V0Decision)esdV0->GetV0ADecision();
+      else if (side == kCSide) return (V0Decision)esdV0->GetV0CDecision();
+      else return kV0Invalid;
+    }
+  }
+
   Float_t time = 0;
   Float_t weight = 0;
   if (fMC)
@@ -708,7 +843,7 @@ AliTriggerAnalysis::V0Decision AliTriggerAnalysis::V0Trigger(const AliESDEvent*
     else {
       for (Int_t i = begin; i < end; ++i) {
        if (esdV0->GetTime(i) > 1e-6 && esdV0->GetAdc(i) > fV0AdcThr) {
-         Float_t correctedTime = V0CorrectLeadingTime(i, esdV0->GetTime(i), esdV0->GetAdc(i));
+         Float_t correctedTime = V0CorrectLeadingTime(i, esdV0->GetTime(i), esdV0->GetAdc(i),aEsd->GetRunNumber());
          Float_t timeWeight = V0LeadingTimeWeight(esdV0->GetAdc(i));
          time += correctedTime*timeWeight;
             
@@ -753,7 +888,7 @@ AliTriggerAnalysis::V0Decision AliTriggerAnalysis::V0Trigger(const AliESDEvent*
   return kV0Empty;
 }
 
-Float_t AliTriggerAnalysis::V0CorrectLeadingTime(Int_t i, Float_t time, Float_t adc) const
+Float_t AliTriggerAnalysis::V0CorrectLeadingTime(Int_t i, Float_t time, Float_t adc, Int_t runNumber) const
 {
   // Correct for slewing and align the channels
   //
@@ -764,7 +899,8 @@ Float_t AliTriggerAnalysis::V0CorrectLeadingTime(Int_t i, Float_t time, Float_t
   // Time alignment
   Float_t timeShift[64] = {0.477957 , 0.0889999 , 0.757669 , 0.205439 , 0.239666 , -0.183705 , 0.442873 , -0.281366 , 0.260976 , 0.788995 , 0.974758 , 0.548532 , 0.495023 , 0.868472 , 0.661167 , 0.358307 , 0.221243 , 0.530179 , 1.26696 , 1.33082 , 1.27086 , 1.77133 , 1.10253 , 0.634806 , 2.14838 , 1.50212 , 1.59253 , 1.66122 , 1.16957 , 1.52056 , 1.47791 , 1.81905 , -1.94123 , -1.29124 , -2.16045 , -1.78939 , -3.11111 , -1.87178 , -1.57671 , -1.70311 , -1.81208 , -1.94475 , -2.53058 , -1.7042 , -2.08109 , -1.84416 , -0.61073 , -1.77145 , 0.16999 , -0.0585339 , 0.00401133 , 0.397726 , 0.851111 , 0.264187 , 0.59573 , -0.158263 , 0.584362 , 1.20835 , 0.927573 , 1.13895 , 0.64648 , 2.18747 , 1.68909 , 0.451194};
 
-  time -= timeShift[i];
+  if(runNumber < 106031)
+    time -= timeShift[i];
 
   // Slewing correction
   if (adc == 0) return time;
@@ -772,6 +908,7 @@ Float_t AliTriggerAnalysis::V0CorrectLeadingTime(Int_t i, Float_t time, Float_t
   Float_t p1 = 1.57345e1;
   Float_t p2 =-4.25603e-1;
 
+  if(runNumber >= 106031) adc *= (2.5/4.0);
   return (time - p1*TMath::Power(adc,p2));
 }
 
@@ -786,6 +923,103 @@ Float_t AliTriggerAnalysis::V0LeadingTimeWeight(Float_t adc) const
   return 1./(p1*p1*TMath::Power(adc,2.*(p2-1.))+p3*p3);
 }
 
+
+Bool_t AliTriggerAnalysis::ZDCTDCTrigger(const AliESDEvent* aEsd, AliceSide side, Bool_t useZN, Bool_t useZP, Bool_t fillHists) const
+{
+  // Returns if ZDC triggered, based on TDC information 
+  
+  AliESDZDC *esdZDC = aEsd->GetESDZDC();
+
+  Bool_t zdcNA = kFALSE;
+  Bool_t zdcNC = kFALSE;
+  Bool_t zdcPA = kFALSE;
+  Bool_t zdcPC = kFALSE;
+
+  if (fMC) {
+        // If it's MC, we use the energy
+    Double_t minEnergy = 0;
+    Double_t  zNCEnergy = esdZDC->GetZDCN1Energy();
+    Double_t  zPCEnergy = esdZDC->GetZDCP1Energy();
+    Double_t  zNAEnergy = esdZDC->GetZDCN2Energy();
+    Double_t  zPAEnergy = esdZDC->GetZDCP2Energy();
+    zdcNA = (zNAEnergy>minEnergy);
+    zdcNC = (zNCEnergy>minEnergy);
+    zdcPA = (zPAEnergy>minEnergy);
+    zdcPC = (zPCEnergy>minEnergy);
+
+  }
+  else {
+
+    Bool_t tdc[32] = {kFALSE};
+    for(Int_t itdc=0; itdc<32; itdc++){
+      for(Int_t i=0; i<4; i++){
+       if (esdZDC->GetZDCTDCData(itdc, i) != 0){
+         tdc[itdc] = kTRUE;
+       }
+      }
+      if(fillHists && tdc[itdc]) {
+       fHistTDCZDC->Fill(itdc);
+      }    
+    }
+    zdcNA = tdc[12];
+    zdcNC = tdc[10];
+    zdcPA = tdc[13];
+    zdcPC = tdc[11];
+  }
+
+  if (side == kASide) return ((useZP&&zdcPA) || (useZN&&zdcNA)); 
+  if (side == kCSide) return ((useZP&&zdcPC) || (useZN&&zdcNC)); 
+  return kFALSE;
+}
+
+Bool_t AliTriggerAnalysis::ZDCTimeTrigger(const AliESDEvent *aEsd, Bool_t fillHists) const
+{
+  // This method implements a selection
+  // based on the timing in both sides of zdcN
+  // It can be used in order to eliminate
+  // parasitic collisions
+  Bool_t zdcAccept = kFALSE;
+  AliESDZDC *esdZDC = aEsd->GetESDZDC();
+
+  if(fMC) {
+     UInt_t esdFlag =  esdZDC->GetESDQuality();
+
+     Bool_t znaFired=kFALSE, zpaFired=kFALSE;
+     Bool_t zem1Fired=kFALSE, zem2Fired=kFALSE;
+     Bool_t zncFired=kFALSE, zpcFired=kFALSE;
+     //
+     // **** Trigger patterns
+     if((esdFlag & 0x00000001) == 0x00000001) znaFired=kTRUE;
+     if((esdFlag & 0x00000002) == 0x00000002) zpaFired=kTRUE;
+     if((esdFlag & 0x00000004) == 0x00000004) zem1Fired=kTRUE;
+     if((esdFlag & 0x00000008) == 0x00000008) zem2Fired=kTRUE;
+     if((esdFlag & 0x00000010) == 0x00000010) zncFired=kTRUE;
+     if((esdFlag & 0x00000020) == 0x00000020) zpcFired=kTRUE;
+     zdcAccept = (znaFired | zncFired);
+  }
+  else {
+    const Float_t refSum = -568.5;
+    const Float_t refDelta = -2.1;
+    const Float_t sigmaSum = 3.25;
+    const Float_t sigmaDelta = 2.25;
+    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)); 
+       for(Int_t j = 0; j < 4; ++j) {
+         if (esdZDC->GetZDCTDCData(12,j) != 0) {
+           Float_t tdcA = 0.025*(esdZDC->GetZDCTDCData(12,j)-esdZDC->GetZDCTDCData(14,j));
+           if(fillHists) fHistTimeZDC->Fill(tdcC-tdcA,tdcC+tdcA);
+           if (((tdcC-tdcA-refDelta)*(tdcC-tdcA-refDelta)/(sigmaDelta*sigmaDelta) +
+                (tdcC+tdcA-refSum)*(tdcC+tdcA-refSum)/(sigmaSum*sigmaSum))< 1.0)
+             zdcAccept = kTRUE;
+         }
+       }
+      }
+    }
+  }
+  return zdcAccept;
+}
+
 Bool_t AliTriggerAnalysis::ZDCTrigger(const AliESDEvent* aEsd, AliceSide side) const
 {
   // Returns if ZDC triggered
@@ -823,6 +1057,9 @@ Int_t AliTriggerAnalysis::FMDHitCombinations(const AliESDEvent* aEsd, AliceSide
   //
   // Authors: FMD team, Hans Dalsgaard (code merged from FMD/AliFMDOfflineTrigger)
 
+  if (!fDoFMD)
+    return -1;
+
   // Workaround for AliESDEvent::GetFMDData is not const!
   const AliESDFMD* fmdData = (const_cast<AliESDEvent*>(aEsd))->GetFMDData();
   if (!fmdData)
@@ -900,7 +1137,7 @@ Long64_t AliTriggerAnalysis::Merge(TCollection* list)
   TObject* obj;
 
   // collections of all histograms
-  const Int_t nHists = 9;
+  const Int_t nHists = 11;
   TList collections[nHists];
 
   Int_t count = 0;
@@ -914,6 +1151,8 @@ Long64_t AliTriggerAnalysis::Merge(TCollection* list)
     collections[n++].Add(entry->fHistV0A);
     collections[n++].Add(entry->fHistV0C);
     collections[n++].Add(entry->fHistZDC);
+    collections[n++].Add(entry->fHistTDCZDC);
+    collections[n++].Add(entry->fHistTimeZDC);
     collections[n++].Add(entry->fHistFMDA);
     collections[n++].Add(entry->fHistFMDC);
     collections[n++].Add(entry->fHistFMDSingle);
@@ -949,6 +1188,8 @@ Long64_t AliTriggerAnalysis::Merge(TCollection* list)
   fHistV0A->Merge(&collections[n++]);
   fHistV0C->Merge(&collections[n++]);
   fHistZDC->Merge(&collections[n++]);
+  fHistTDCZDC->Merge(&collections[n++]);
+  fHistTimeZDC->Merge(&collections[n++]);
   fHistFMDA->Merge(&collections[n++]);
   fHistFMDC->Merge(&collections[n++]);
   fHistFMDSingle->Merge(&collections[n++]);
@@ -968,20 +1209,34 @@ void AliTriggerAnalysis::SaveHistograms() const
   if (!fHistBitsSPD)
     return;
     
-  fHistBitsSPD->Write();
-  fHistBitsSPD->ProjectionX();
-  fHistBitsSPD->ProjectionY();
-  fHistFiredBitsSPD->Write();
-  fHistV0A->Write();
-  fHistV0C->Write();
-  fHistZDC->Write();
-  fHistFMDA->Write();
-  fHistFMDC->Write();
-  fHistFMDSingle->Write();
-  fHistFMDSum->Write();
-  
-  if (fSPDGFOEfficiency)
-    fSPDGFOEfficiency->Write("fSPDGFOEfficiency");
+  if (fHistBitsSPD) {
+    fHistBitsSPD->Write();
+    fHistBitsSPD->ProjectionX();
+    fHistBitsSPD->ProjectionY();
+  }
+  else Printf("Cannot save fHistBitsSPD");
+  if (fHistFiredBitsSPD) fHistFiredBitsSPD->Write();
+  else Printf("Cannot save fHistFiredBitsSPD");
+  if (fHistV0A) fHistV0A->Write();
+  else Printf("Cannot save fHistV0A");
+  if (fHistV0C) fHistV0C->Write();
+  else Printf("Cannot save fHistV0C");
+  if (fHistZDC) fHistZDC->Write();
+  else Printf("Cannot save fHistZDC");
+  if (fHistTDCZDC) fHistTDCZDC->Write();
+  else Printf("Cannot save fHistTDCZDC");
+  if (fHistTimeZDC) fHistTimeZDC->Write();
+  else Printf("Cannot save fHistTimeZDC");
+  if (fHistFMDA) fHistFMDA->Write();
+  else Printf("Cannot save fHistFMDA");
+  if (fHistFMDC) fHistFMDC->Write();
+  else Printf("Cannot save fHistFMDC");
+  if (fHistFMDSingle) fHistFMDSingle->Write();
+  else Printf("Cannot save fHistFMDSingle");
+  if (fHistFMDSum) fHistFMDSum->Write();
+  else Printf("Cannot save fHistFMDSum");
+  if (fSPDGFOEfficiency) fSPDGFOEfficiency->Write("fSPDGFOEfficiency");
+  else Printf("Cannot save fSPDGFOEfficiency");
   
   fTriggerClasses->Write("fTriggerClasses", TObject::kSingleKey);
 }
@@ -1003,7 +1258,7 @@ void AliTriggerAnalysis::PrintTriggerClasses() const
   {
     TParameter<Long64_t>* param = static_cast<TParameter<Long64_t>*> (fTriggerClasses->GetValue(obj));
     
-    Printf(" %s: %ld triggers", obj->String().Data(), param->GetVal());
+    Printf(" %s: %ld triggers", obj->String().Data(), (Long_t)param->GetVal());
     
     TObjArray* tokens = obj->String().Tokenize(" ");
     for (Int_t i=0; i<tokens->GetEntries(); i++)
@@ -1028,7 +1283,7 @@ void AliTriggerAnalysis::PrintTriggerClasses() const
   {
     TParameter<Long64_t>* param = static_cast<TParameter<Long64_t>*> (singleTrigger.GetValue(obj));
     
-    Printf("  %s: %ld triggers", obj->String().Data(), param->GetVal());
+    Printf("  %s: %ld triggers", obj->String().Data(), (Long_t)param->GetVal());
   }
   delete iter;