]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG0/AliPWG0Helper.cxx
Added forward declaration
[u/mrichter/AliRoot.git] / PWG0 / AliPWG0Helper.cxx
index 13e01024122930435e726d6d90e8a0387dca93e0..5466eec538dace8ceeb9d4bd0dc6e1a6a90bcab3 100644 (file)
@@ -15,7 +15,7 @@
 ClassImp(AliPWG0Helper)
 
 //____________________________________________________________________
-Bool_t AliPWG0Helper::IsEventTriggered(AliESD* aEsd)
+Bool_t AliPWG0Helper::IsEventTriggered(AliESD* aEsd, Trigger trigger)
 {
   // check if the event was triggered
   //
@@ -28,8 +28,21 @@ Bool_t AliPWG0Helper::IsEventTriggered(AliESD* aEsd)
 
   ULong64_t triggerMask = aEsd->GetTriggerMask();
 
-  if (triggerMask&32 && ((triggerMask&1) || (triggerMask&2)))
-    return kTRUE;
+  switch (trigger)
+  {
+    case kMB1:
+    {
+      if (triggerMask&32 || ((triggerMask&1) || (triggerMask&2)))
+        return kTRUE;
+      break;
+    }
+    case kMB2:
+    {
+      if (triggerMask&32 && ((triggerMask&1) || (triggerMask&2)))
+        return kTRUE;
+      break;
+    }
+  }
 
   return kFALSE;
 }
@@ -56,6 +69,8 @@ Bool_t AliPWG0Helper::IsVertexReconstructed(AliESD* aEsd)
   if (vtx_res[2]==0 || vtx_res[2]>0.1)
     return kFALSE;
 
+  // check Ncontributors, if <0 it means error *gna*
+
   return kTRUE;
 }
 
@@ -176,7 +191,7 @@ void AliPWG0Helper::CreateDividedProjections(TH3* hist, TH3* hist2, const char*
   TH1* division = dynamic_cast<TH1*> (proj->Clone(Form("%s_div_%s", proj->GetName(), proj2->GetName())));
   //printf("doing axis: %s, x axis has %d %d bins, min %f %f max %f %f\n", axis, division->GetNbinsX(), proj2->GetNbinsX(), division->GetXaxis()->GetBinLowEdge(1), proj2->GetXaxis()->GetBinLowEdge(1), division->GetXaxis()->GetBinUpEdge(division->GetNbinsX()), proj2->GetXaxis()->GetBinUpEdge(proj2->GetNbinsX()));
   //printf("doing axis: %s, y axis has %d %d bins, min %f %f max %f %f\n", axis, division->GetNbinsY(), proj2->GetNbinsY(), division->GetYaxis()->GetBinLowEdge(1), proj2->GetYaxis()->GetBinLowEdge(1), division->GetYaxis()->GetBinUpEdge(division->GetNbinsY()), proj2->GetYaxis()->GetBinUpEdge(proj2->GetNbinsY()));
-  division->Divide(proj2);
+  division->Divide(proj, proj2, 1, 1, "B");
   division->SetTitle(Form("%s divided %s", proj->GetTitle(), proj2->GetTitle()));
 
   if (putErrors)