author agheata Mon, 23 Jul 2012 12:49:42 +0000 (12:49 +0000) committer agheata Mon, 23 Jul 2012 12:49:42 +0000 (12:49 +0000)
orbits = 275952095
nb*orbits = 12969748465 = 3 * 2^32 + 84846577

index 10b42da..ca0ba42 100644 (file)
@@ -593,7 +593,7 @@ Bool_t AliTriggerRunScalers::CalculateMu(Double_t &mu, Double_t &errmu, ULong64_
{

if (nB!=0 && orbits!=0)  {
-      Double_t pB = (Double_t)countsB/(nB*orbits); // probability for B trigger
+      Double_t pB = (Double_t)countsB/((ULong64_t)nB*orbits); // probability for B trigger
if (!bkgCorr || nAC==0 ) {
mu = -log(1-pB)/triggerEff;
errmu = TMath::Sqrt(pB/((1-pB)*nB*orbits) + mu*mu*errorEff*errorEff/(triggerEff*triggerEff)); //
@@ -601,7 +601,7 @@ Bool_t AliTriggerRunScalers::CalculateMu(Double_t &mu, Double_t &errmu, ULong64_
}
else
{
-         Double_t pAC = (Double_t)countsAC/(nAC*orbits); // probability for AC trigger (background)
+         Double_t pAC = (Double_t)countsAC/((ULong64_t)nAC*orbits); // probability for AC trigger (background)
mu = ( log(1.-pAC) - log(1.-pB) )/triggerEff;
// error
errmu =  TMath::Sqrt(pB/((1.-pB)*nB*orbits) + pAC/((1.-pAC)*nAC*orbits) /*- 2*TMath::Sqrt(pB*pAC/(nB*nAC*(1.-pB)*(1.-pAC)))/orbits*/ + mu*mu*errorEff*errorEff/(triggerEff*triggerEff)); // assume no correlation between B and AC rates, hence no cov term in error
@@ -623,7 +623,7 @@ Bool_t AliTriggerRunScalers::CalculateMu(Double_t &mu, Double_t &errmu, ULong64_
}
else
{
-         Double_t pAC = (Double_t)countsAC/(nAC*beamB/nB); // probability for AC trigger (background)
+         Double_t pAC = (Double_t)countsAC/((ULong64_t)nAC*beamB/(Double_t)nB); // probability for AC trigger (background)
mu = ( log(1-pAC) - log(1-pB) )/triggerEff;
// error
errmu =  TMath::Sqrt(pB/((1-pB)*beamB) + pAC/((1-pAC)*nAC*beamB/nB) + mu*mu*errorEff*errorEff/(triggerEff*triggerEff));