Changeset 19874

Timestamp:

11/18/19 11:33:03 (5 years ago)

Author:

tbretz

Message:

Just found out that this was the wrong version of the macro. The correct one (which also contains the correct function from the paper) is this one.

File:

• trunk/Mars/fact/analysis/gain/fit_spectra.C (modified) (7 diffs)

trunk/Mars/fact/analysis/gain/fit_spectra.C

@@ -23 +23 @@
 
 // --------------------------------------------------------------------------
+// After fit_spectra_bt2b.C
 
 // Fit function for a single pe spectrum
@@ -35 +36 @@
     const Double_t expo  = par[6];
 
+    const Double_t P = cross*TMath::Exp(-cross);
+
     Double_t y = 0;
     for (int N=1; N<14; N++)
@@ -41 +44 @@
         const Double_t sigN = TMath::Sqrt(N*sigma*sigma + noise*noise);
 
-        const Double_t p = TMath::Power(cross, N-1) * TMath::Power(N, -expo);
+        const Double_t p =
+            TMath::Power(N*P, N-1)/TMath::Power(TMath::Factorial(N-1), expo);
 
         y += TMath::Gaus(xx[0], muN, sigN) * p / sigN;
@@ -56 +60 @@
     Double_t expo  = f.GetParameter(6);
 
+    const Double_t P = cross*TMath::Exp(-cross);
+
     Double_t y = 0;
     for (int N=2; N<14; N++)
-        y +=  TMath::Power(cross,  N-1) * TMath::Power(N, -expo);
+        y += 
+            TMath::Power(N*P, N-1)/TMath::Power(TMath::Factorial(N-1), expo);
 
     return y / (y + 1);
@@ -71 +78 @@
     const Double_t expo  = func.GetParameter(6);
 
+    const Double_t P = cross*TMath::Exp(-cross);
+
     Double_t sum = 0;
     for (int N=1; N<14; N++)
-        sum += TMath::Power(cross, N-1) * TMath::Power(N, -expo);
+        sum +=
+            TMath::Power(N*P, N-1)/TMath::Power(TMath::Factorial(N-1), expo);
 
     const Double_t scale = hist.GetBinWidth(1);
@@ -383 +393 @@
     if (!fixednoise)
         func.SetParLimits(5, 0, 150);          // Noise
+    func.SetParLimits(6, 0, 2);            // Expo
 
     func.SetParameter(0, ampl);                         // Amplitude
@@ -394 +405 @@
         func.SetParameter(5, 0.1*maxpos);               // Noise
 
-    func.SetParameter(6, 0.4);                          // Expo
+    func.SetParameter(6, 0.95);                         // Expo
 
     func.SetRange(fitmin, fitmax);