Changeset bc58cf5 in git for modules/MomentumSmearing.cc

Timestamp:

Aug 25, 2016, 4:33:41 PM (8 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

79d4123, b3bd4d2

Parents:

2264876

Message:

put LogNormal in MomentumSmearing

File:

• modules/MomentumSmearing.cc (modified) (2 diffs)

modules/MomentumSmearing.cc

@@ -109 +109 @@
 
     // apply smearing formula
-    pt = gRandom->Gaus(pt, fFormula->Eval(pt, eta, phi, e) * pt);
+    //pt = gRandom->Gaus(pt, fFormula->Eval(pt, eta, phi, e) * pt);
+    pt = LogNormal(pt, fFormula->Eval(pt, eta, phi, e) * pt );
     
-    if(pt <= 0.0) continue;
+    //if(pt <= 0.0) continue;
 
     mother = candidate;
@@ -124 +125 @@
   }
 }
+//----------------------------------------------------------------
+
+Double_t MomentumSmearing::LogNormal(Double_t mean, Double_t sigma)
+{
+  Double_t a, b;
+
+  if(mean > 0.0)
+  {
+    b = TMath::Sqrt(TMath::Log((1.0 + (sigma*sigma)/(mean*mean))));
+    a = TMath::Log(mean) - 0.5*b*b;
+
+    return TMath::Exp(a + b*gRandom->Gaus(0.0, 1.0));
+  }
+  else
+  {
+    return 0.0;
+  }
+}
+
 
 //------------------------------------------------------------------------------