Changeset 8030 for trunk/MagicSoft

Timestamp:

10/09/06 10:28:06 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• NEWS (modified) (1 diff)
• mhflux/MMcSpectrumWeight.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -32 +32 @@
      - fixed a bug which caused the spectra to be normalized at 1GeV
        by default instead of the integral
+     - allow the MC slope to be -1. With the current code this should
+       not be a problem anymore.
 
    * mbase/MMath.cc, mbase/MAGIC.cc:

trunk/MagicSoft/Mars/NEWS

@@ -162 +162 @@
      was not using the integral anymore but a normalization energy of
      1GeV. It is now the integral again.
+
+   - sponde: Fixed a problem in the weighting which did not allow
+     a MC slope of -1.
 
 

trunk/MagicSoft/Mars/mhflux/MMcSpectrumWeight.cc

@@ -107 +107 @@
     fNameMcEvt   = "MMcEvt";
 
-    fNewSlope    = -1;
-    fOldSlope    = -1;
+    fNewSlope    = -9;
+    fOldSlope    = -9;
 
     fEnergyMin   = -1;
@@ -379 +379 @@
     }
 
-    //
-    // Sanity checks to be sure that we won't divide by zero later on
-    //
-    if (rh.GetSlopeSpec()==-1)
-    {
-        *fLog << err << "The MC Slope of the power law must be different of -1... exit" << endl;
-        return kFALSE;
-    }
-
     fOldSlope  = rh.GetSlopeSpec();
     fEnergyMin = rh.GetELowLim();
     fEnergyMax = rh.GetEUppLim();
 
-    if (fNewSlope==-1)
-    {
-        *fLog << inf << "The new slope of the power law is -1... no weighting applied." << endl;
+    if (fNewSlope==-9)
+    {
+        *fLog << inf << "The new slope of the power law is undefined (-9)... no weighting applied." << endl;
         fNewSlope = fOldSlope;
     }
 
     TString form(GetFormulaWeightsX());
-
     if (fFunc)
         delete fFunc;