Changeset 7522 for trunk/MagicSoft/Mars

Timestamp:

02/23/06 12:54:27 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• NEWS (modified) (1 diff)
• mjobs/MJSpectrum.cc (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -18 +18 @@
 
                                                  -*-*- END OF LINE -*-*-
+ 2006/02/23 Thomas Bretz
+
+   * mjobs/MJSpectrum.cc:
+     - fixed plot for Eý.dN/dE. The spectrum was multiplied with two
+       times the bin-width instead of the absolute energy
+
+
+
  2006/02/22 Daniela Dorner
 

trunk/MagicSoft/Mars/NEWS

@@ -7 +7 @@
    - ganymed: now stops in requested source isn't found in catalog file
 
+   - ganymed: fixed plot for Eý.dN/dE. The spectrum was multiplied with two
+     times the bin-width instead of the absolute energy
 
 

trunk/MagicSoft/Mars/mjobs/MJSpectrum.cc

@@ -700 +700 @@
     for (int i=0; i<spectrum.GetNbinsX(); i++)
     {
-        spectrum.SetBinContent(i+1, spectrum.GetBinContent(i+1)*spectrum.GetBinWidth(i+1)*spectrum.GetBinWidth(i+1));
-        spectrum.SetBinError(  i+1, spectrum.GetBinError(i+1)  *spectrum.GetBinWidth(i+1)*spectrum.GetBinWidth(i+1));
+        const Double_t e = TMath::Sqrt(spectrum.GetBinLowEdge(i+1)*spectrum.GetBinLowEdge(i+2))*1e-3;
+
+        spectrum.SetBinContent(i+1, spectrum.GetBinContent(i+1)*e*e);
+        spectrum.SetBinError(  i+1, spectrum.GetBinError(i+1)  *e*e);
     }