]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ANALYSIS/AliTriggerAnalysis.cxx
warnings fixed
[u/mrichter/AliRoot.git] / ANALYSIS / AliTriggerAnalysis.cxx
index 6d098438807b534ea43a021914eca2ee676cc077..21742b8665ac4b5d92e34d28667897878b4e1a93 100644 (file)
@@ -29,6 +29,7 @@
 #include <TIterator.h>
 #include "TParameter.h"
 #include <TMap.h>
+#include <TRandom.h>
 
 #include <AliTriggerAnalysis.h>
 
 #include <AliESDVZERO.h>
 #include <AliESDZDC.h>
 #include <AliESDFMD.h>
+#include <AliESDVertex.h>
 
 ClassImp(AliTriggerAnalysis)
 
 AliTriggerAnalysis::AliTriggerAnalysis() :
   fSPDGFOThreshold(2),
+  fSPDGFOEfficiency(0),
   fV0TimeOffset(0),
+  fV0AdcThr(0),
+  fV0HwAdcThr(2.5),
+  fV0HwWinLow(61.5),
+  fV0HwWinHigh(86.5),
   fFMDLowCut(0.2),
   fFMDHitCut(0.5),
   fHistBitsSPD(0),
@@ -57,7 +64,8 @@ AliTriggerAnalysis::AliTriggerAnalysis() :
   fHistFMDC(0),   
   fHistFMDSingle(0),
   fHistFMDSum(0),
-  fTriggerClasses(0)
+  fTriggerClasses(0),
+  fMC(kFALSE)
 {
   // constructor
 }
@@ -122,6 +130,7 @@ AliTriggerAnalysis::~AliTriggerAnalysis()
 
   if (fTriggerClasses)
   {
+    fTriggerClasses->DeleteAll();
     delete fTriggerClasses;
     fTriggerClasses = 0;
   }
@@ -137,8 +146,8 @@ 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);
   
   // TODO check limits
@@ -173,6 +182,8 @@ const char* AliTriggerAnalysis::GetTriggerName(Trigger trigger)
     case kSPDGFOBits : str = "SPD GFO Bits"; break;
     case kV0A : str = "V0 A BB"; break;
     case kV0C : str = "V0 C BB"; break;
+    case kV0OR : str = "V0 OR BB"; break;
+    case kV0AND : str = "V0 AND BB"; break;
     case kV0ABG : str = "V0 A BG"; break;
     case kV0CBG : str = "V0 C BG"; break;
     case kZDC : str = "ZDC"; break;
@@ -181,12 +192,17 @@ const char* AliTriggerAnalysis::GetTriggerName(Trigger trigger)
     case kFMDA : str = "FMD A"; break;
     case kFMDC : str = "FMD C"; break;
     case kFPANY : str = "SPD GFO | V0 | ZDC | FMD"; break;
+    case kNSD1 : str = "NSD1"; break;
+    case kMB1Prime: str = "MB1prime"; break;
     default: str = ""; break;
   }
    
   if (trigger & kOfflineFlag)
     str += " OFFLINE";  
   
+  if (trigger & kOneParticle)
+    str += " OneParticle";  
+
   return str;
 }
 
@@ -271,101 +287,135 @@ 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) == kV0BB || V0Trigger(aEsd, kCSide) == kV0BB)
-        return kTRUE;
+      if (SPDGFOTrigger(aEsd, 0) || V0Trigger(aEsd, kASide, kFALSE) == kV0BB || V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
+        decision = kTRUE;
       break;
     }
     case kMB2:
     {
-      if (SPDGFOTrigger(aEsd, 0) && (V0Trigger(aEsd, kASide) == kV0BB || V0Trigger(aEsd, kCSide) == kV0BB))
-        return kTRUE;
+      if (SPDGFOTrigger(aEsd, 0) && (V0Trigger(aEsd, kASide, kFALSE) == kV0BB || V0Trigger(aEsd, kCSide, kFALSE) == kV0BB))
+        decision = kTRUE;
       break;
     }
     case kMB3:
     {
-      if (SPDGFOTrigger(aEsd, 0) && V0Trigger(aEsd, kASide) == kV0BB && V0Trigger(aEsd, kCSide) == kV0BB)
-        return kTRUE;
+      if (SPDGFOTrigger(aEsd, 0) && V0Trigger(aEsd, kASide, kFALSE) == kV0BB && V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
+        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) == kV0BB)
-        return kTRUE;
+      if (V0Trigger(aEsd, kASide, kFALSE) == kV0BB)
+        decision = kTRUE;
       break;
     }
     case kV0C:
     {
-      if (V0Trigger(aEsd, kCSide) == kV0BB)
-        return kTRUE;
+      if (V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
+        decision = kTRUE;
+      break;
+    }
+    case kV0OR:
+    {
+      if (V0Trigger(aEsd, kASide, kFALSE) == kV0BB || V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
+        decision = kTRUE;
+      break;
+    }
+    case kV0AND:
+    {
+      if (V0Trigger(aEsd, kASide, kFALSE) == kV0BB && V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
+        decision = kTRUE;
       break;
     }
     case kV0ABG:
     {
-      if (V0Trigger(aEsd, kASide) == kV0BG)
-        return kTRUE;
+      if (V0Trigger(aEsd, kASide, kFALSE) == kV0BG)
+        decision = kTRUE;
       break;
     }
     case kV0CBG:
     {
-      if (V0Trigger(aEsd, kCSide) == kV0BG)
-        return kTRUE;
+      if (V0Trigger(aEsd, kCSide, kFALSE) == kV0BG)
+        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) == kV0BB || V0Trigger(aEsd, kCSide) == kV0BB || ZDCTrigger(aEsd, kASide) || ZDCTrigger(aEsd, kCentralBarrel) || ZDCTrigger(aEsd, kCSide) || FMDTrigger(aEsd, kASide) || FMDTrigger(aEsd, kCSide))
-        return kTRUE;
+      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))
+        decision = kTRUE;
+      break;
+    }
+    case kNSD1:
+    {
+      if (SPDFiredChips(aEsd, 0) >= 5 || (V0Trigger(aEsd, kASide, kFALSE) == kV0BB && V0Trigger(aEsd, kCSide, kFALSE) == kV0BB))
+        decision = kTRUE;
+       break;
+    }
+    case kMB1Prime:
+    {
+      Int_t count = 0;
+      if (SPDGFOTrigger(aEsd, 0))
+        count++;
+      if (V0Trigger(aEsd, kASide, kFALSE) == kV0BB)
+        count++;
+      if (V0Trigger(aEsd, kCSide, kFALSE) == kV0BB)
+        count++;
+      
+      if (count >= 2)
+        decision = kTRUE;
+        
       break;
     }
     default:
@@ -374,9 +424,29 @@ 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;
+        }
+      }
+    }
+  }
 
+  return decision;
+}
 
 Bool_t AliTriggerAnalysis::IsTriggerClassFired(const AliESDEvent* aEsd, const Char_t* tclass) const 
 {
@@ -457,8 +527,8 @@ void AliTriggerAnalysis::FillHistograms(const AliESDEvent* aEsd)
   
   fHistBitsSPD->Fill(SPDFiredChips(aEsd, 0), SPDFiredChips(aEsd, 1, kTRUE));
   
-  V0Trigger(aEsd, kASide, kTRUE);
-  V0Trigger(aEsd, kCSide, kTRUE);
+  V0Trigger(aEsd, kASide, kFALSE, kTRUE);
+  V0Trigger(aEsd, kCSide, kFALSE, kTRUE);
   
   AliESDZDC* zdcData = aEsd->GetESDZDC();
   if (zdcData)
@@ -505,6 +575,15 @@ void AliTriggerAnalysis::FillTriggerClasses(const AliESDEvent* aEsd)
   // TODO add first and last orbit number here
 }
 
+Int_t AliTriggerAnalysis::SSDClusters(const AliESDEvent* aEsd)
+{
+  // returns the number of clusters in the SSD
+  const AliMultiplicity* mult = aEsd->GetMultiplicity();
+  Int_t clusters = mult->GetNumberOfITSClusters(4)+mult->GetNumberOfITSClusters(5);
+  return clusters;
+}
+
+
 Int_t AliTriggerAnalysis::SPDFiredChips(const AliESDEvent* aEsd, Int_t origin, Bool_t fillHists)
 {
   // returns the number of fired chips in the SPD
@@ -526,12 +605,21 @@ Int_t AliTriggerAnalysis::SPDFiredChips(const AliESDEvent* aEsd, Int_t origin, B
   {
     Int_t nChips = 0;
     for (Int_t i=0; i<1200; i++)
+    {
       if (mult->TestFastOrFiredChips(i) == kTRUE)
       {
+        // efficiency simulation (if enabled)
+        if (fSPDGFOEfficiency)
+        {
+          if (gRandom->Uniform() > fSPDGFOEfficiency->GetBinContent(i+1))
+            continue;
+        }
+        
         nChips++;
         if (fillHists)
           fHistFiredBitsSPD->Fill(i);
       }
+    }
     return nChips;
   }
   
@@ -549,11 +637,17 @@ Bool_t AliTriggerAnalysis::SPDGFOTrigger(const AliESDEvent* aEsd, Int_t origin)
   return kFALSE;
 }
 
-AliTriggerAnalysis::V0Decision AliTriggerAnalysis::V0Trigger(const AliESDEvent* aEsd, AliceSide side, Bool_t fillHists)
+AliTriggerAnalysis::V0Decision AliTriggerAnalysis::V0Trigger(const AliESDEvent* aEsd, AliceSide side, Bool_t online, Bool_t fillHists)
 {
   // Returns the V0 trigger decision in V0A | V0C
   //
-  // Based on algorithm by Cvetan Cheshkov
+  // Returns kV0Fake if the calculated average time is in a window where neither BB nor BG is expected. 
+  // The rate of such triggers can be used to estimate the background. Note that the rate has to be 
+  // rescaled with the size of the windows (numerical values see below in the code)
+  //
+  // 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)
@@ -561,7 +655,27 @@ AliTriggerAnalysis::V0Decision AliTriggerAnalysis::V0Trigger(const AliESDEvent*
     AliError("AliESDVZERO not available");
     return kV0Invalid;
   }
-  
+
+  if (esdV0->TestBit(AliESDVZERO::kDecisionFilled)) {
+    if (online) {
+      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;
+    }
+  }
+
   Int_t begin = -1;
   Int_t end = -1;
   
@@ -577,16 +691,75 @@ AliTriggerAnalysis::V0Decision AliTriggerAnalysis::V0Trigger(const AliESDEvent*
   }
   else
     return kV0Invalid;
-  
+    
   Float_t time = 0;
   Float_t weight = 0;
-  for (Int_t i = begin; i < end; ++i) {
-    if (esdV0->GetTime(i) > 1e-6) {
-      Float_t correctedTime = V0CorrectLeadingTime(i, esdV0->GetTime(i), esdV0->GetAdc(i));
-      Float_t timeWeight = V0LeadingTimeWeight(esdV0->GetAdc(i));
-      time += correctedTime*timeWeight;
-      
-      weight += timeWeight;
+  if (fMC)
+  {
+    Int_t runRange;
+    if (aEsd->GetRunNumber() <= 104803) runRange = 0;
+    else if (aEsd->GetRunNumber() <= 104876) runRange = 1;
+    else runRange = 2;
+
+    Float_t factors[3][64] = {
+      // runs: 104792-104803
+      {4.6,5.9,6.3,6.0,4.7,5.9,4.9,5.4,4.8,4.1,4.9,4.6,4.5,5.5,5.1,5.8,4.3,4.0,4.0,3.3,3.1,2.9,3.0,5.6,3.3,4.9,3.9,5.3,4.1,4.4,3.9,5.5,5.7,9.5,5.1,5.3,6.6,7.1,8.9,4.4,4.1,5.9,9.0,4.5,4.1,6.0,4.7,7.1,4.2,4.7,3.9,6.3,5.9,4.8,4.7,4.5,4.7,5.4,5.8,5.0,5.1,5.9,5.3,3.6},
+      // runs: 104841-104876
+      {4.6,4.8,4.9,4.8,4.3,4.9,4.4,4.5,4.6,5.0,4.7,4.6,4.7,4.6,4.6,5.5,4.7,4.5,4.7,5.0,6.5,7.6,5.3,4.9,5.5,4.8,4.6,4.9,4.5,4.5,4.6,4.9,5.7,9.8,4.9,5.2,7.1,7.1,8.1,4.4,4.0,6.0,8.3,4.6,4.2,5.6,4.6,6.4,4.4,4.7,4.5,6.5,6.0,4.7,4.5,4.4,4.8,5.5,5.9,5.3,5.0,5.7,5.1,3.6},
+      // runs: 104890-92
+      {4.7,5.2,4.8,5.0,4.4,5.0,4.4,4.6,4.6,4.5,4.4,4.6,4.5,4.6,4.8,5.5,4.8,4.5,4.4,4.3,5.4,7.7,5.6,5.0,5.4,4.3,4.5,4.8,4.5,4.5,4.6,5.3,5.7,9.6,4.9,5.4,6.1,7.2,8.6,4.4,4.0,5.4,8.8,4.4,4.2,5.8,4.7,6.7,4.3,4.7,4.0,6.1,6.0,4.9,4.8,4.6,4.7,5.2,5.7,5.0,5.0,5.8,5.3,3.6}
+    };
+    Float_t dA = 77.4 - 11.0;
+    Float_t dC = 77.4 - 2.9;
+    // Time misalignment
+    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};
+    Float_t dA2 = 2.8, dC2 = 3.3;
+
+    if (online) {
+      for (Int_t i = begin; i < end; ++i) {
+        Float_t tempAdc = esdV0->GetAdc(i)/factors[runRange][i];
+        Float_t tempTime = (i >= 32) ? esdV0->GetTime(i)+dA+timeShift[i]+dA2 : esdV0->GetTime(i)+dC+timeShift[i]+dC2;
+        if (esdV0->GetTime(i) >= 1e-6 &&
+            tempTime > fV0HwWinLow && tempTime < fV0HwWinHigh &&
+            tempAdc > fV0HwAdcThr)
+          return kV0BB;
+      }
+      return kV0Empty;
+    }
+    else {
+      for (Int_t i = begin; i < end; ++i) {
+        Float_t tempAdc = esdV0->GetAdc(i)/factors[runRange][i];
+        Float_t tempTime = (i >= 32) ? esdV0->GetTime(i)+dA : esdV0->GetTime(i)+dC;
+        Float_t tempRawTime = (i >= 32) ? esdV0->GetTime(i)+dA+timeShift[i]+dA2 : esdV0->GetTime(i)+dC+timeShift[i]+dC2;
+        if (esdV0->GetTime(i) >= 1e-6 &&
+            tempRawTime < 125.0 &&
+            tempAdc > fV0AdcThr) {
+          weight += 1.0;
+          time += tempTime;
+        }
+      }
+    }
+  }
+  else {
+    if (online) {
+      for (Int_t i = begin; i < end; ++i) {
+        if (esdV0->GetTime(i) >= 1e-6 &&
+            esdV0->GetTime(i) > fV0HwWinLow && esdV0->GetTime(i) < fV0HwWinHigh &&
+            esdV0->GetAdc(i) > fV0HwAdcThr)
+          return kV0BB;
+      }
+      return kV0Empty;
+    }
+    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),aEsd->GetRunNumber());
+         Float_t timeWeight = V0LeadingTimeWeight(esdV0->GetAdc(i));
+         time += correctedTime*timeWeight;
+            
+         weight += timeWeight;
+       }
+      }
     }
   }
 
@@ -604,24 +777,28 @@ AliTriggerAnalysis::V0Decision AliTriggerAnalysis::V0Trigger(const AliESDEvent*
   
   if (side == kASide)
   {
-    if (time > 75 && time < 83)
+    if (time > 68 && time < 100)
       return kV0BB;
     if (time > 54 && time < 57.5) 
       return kV0BG;
+    if (time > 57.5 && time < 68)
+      return kV0Fake;
   }
   
   if (side == kCSide)
   {
-    if (time > 75.5 && time < 82)
+    if (time > 75.5 && time < 100)
       return kV0BB;
     if (time > 69.5 && time < 73)
       return kV0BG; 
+    if (time > 55 && time < 69.5)
+      return kV0Fake;
   }
   
   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
   //
@@ -632,7 +809,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;
@@ -640,6 +818,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));
 }
 
@@ -791,12 +970,12 @@ Long64_t AliTriggerAnalysis::Merge(TCollection* list)
 
     // merge fTriggerClasses
     TIterator* iter2 = entry->fTriggerClasses->MakeIterator();
-    TObjString* obj = 0;
-    while ((obj = dynamic_cast<TObjString*> (iter2->Next())))
+    TObjString* obj2 = 0;
+    while ((obj2 = dynamic_cast<TObjString*> (iter2->Next())))
     {
-      TParameter<Long64_t>* param2 = dynamic_cast<TParameter<Long64_t>*> (entry->fTriggerClasses->GetValue(obj));
+      TParameter<Long64_t>* param2 = dynamic_cast<TParameter<Long64_t>*> (entry->fTriggerClasses->GetValue(obj2));
       
-      TParameter<Long64_t>* param1 = dynamic_cast<TParameter<Long64_t>*> (fTriggerClasses->GetValue(obj));
+      TParameter<Long64_t>* param1 = dynamic_cast<TParameter<Long64_t>*> (fTriggerClasses->GetValue(obj2));
       if (param1)
       {
         param1->SetVal(param1->GetVal() + param2->GetVal());
@@ -804,7 +983,7 @@ Long64_t AliTriggerAnalysis::Merge(TCollection* list)
       else
       {
         param1 = dynamic_cast<TParameter<Long64_t>*> (param2->Clone());
-        fTriggerClasses->Add(new TObjString(obj->String()), param1);
+        fTriggerClasses->Add(new TObjString(obj2->String()), param1);
       }
     }
     
@@ -848,6 +1027,9 @@ void AliTriggerAnalysis::SaveHistograms() const
   fHistFMDSingle->Write();
   fHistFMDSum->Write();
   
+  if (fSPDGFOEfficiency)
+    fSPDGFOEfficiency->Write("fSPDGFOEfficiency");
+  
   fTriggerClasses->Write("fTriggerClasses", TObject::kSingleKey);
 }
 
@@ -857,12 +1039,45 @@ void AliTriggerAnalysis::PrintTriggerClasses() const
   
   Printf("Trigger Classes:");
   
+  Printf("Event count for trigger combinations:");
+  
+  TMap singleTrigger;
+  singleTrigger.SetOwner();
+  
   TIterator* iter = fTriggerClasses->MakeIterator();
   TObjString* obj = 0;
   while ((obj = dynamic_cast<TObjString*> (iter->Next())))
   {
-    TParameter<Long64_t>* param = dynamic_cast<TParameter<Long64_t>*> (fTriggerClasses->GetValue(obj));
+    TParameter<Long64_t>* param = static_cast<TParameter<Long64_t>*> (fTriggerClasses->GetValue(obj));
+    
+    Printf(" %s: %lld triggers", obj->String().Data(), param->GetVal());
+    
+    TObjArray* tokens = obj->String().Tokenize(" ");
+    for (Int_t i=0; i<tokens->GetEntries(); i++)
+    {
+      TParameter<Long64_t>* count = dynamic_cast<TParameter<Long64_t>*> (singleTrigger.GetValue(((TObjString*) tokens->At(i))->String().Data()));
+      if (!count)
+      {
+        count = new TParameter<Long64_t>(((TObjString*) tokens->At(i))->String().Data(), 0);
+        singleTrigger.Add(new TObjString(((TObjString*) tokens->At(i))->String().Data()), count);
+      }
+      count->SetVal(count->GetVal() + param->GetVal());
+    }
+    
+    delete tokens;
+  }
+  delete iter;
+  
+  Printf("Event count for single trigger:");
+  
+  iter = singleTrigger.MakeIterator();
+  while ((obj = dynamic_cast<TObjString*> (iter->Next())))
+  {
+    TParameter<Long64_t>* param = static_cast<TParameter<Long64_t>*> (singleTrigger.GetValue(obj));
     
-    Printf("%s: %ld triggers", obj->String().Data(), param->GetVal());
+    Printf("  %s: %lld triggers", obj->String().Data(), param->GetVal());
   }
+  delete iter;
+  
+  singleTrigger.DeleteAll();
 }