Changeset 9337 for trunk/MagicSoft/Mars/mhflux

Timestamp:

02/15/09 13:47:00 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhflux/MMcSpectrumWeight.cc (modified) (11 diffs)

trunk/MagicSoft/Mars/mhflux/MMcSpectrumWeight.cc

@@ -41 +41 @@
 //  -------
 //
-//   - Corsika spectrun: dN/dE = A * E^(a)
+//   - Corsika spectrum: dN/dE = A * E^(a)
 //     with a = fOldSlope, and A = N/integral{E*de} from ELowLim to EUppLim
 //
@@ -90 +90 @@
 #include "MLogManip.h"
 
+#include "MString.h"
 #include "MParList.h"
 #include "MParameters.h"
@@ -204 +205 @@
 TString MMcSpectrumWeight::ReplaceX(TString str) const
 {
-    return str.ReplaceAll(Form("%s.fEnergy", fNameMcEvt.Data()), "(x)");
+    return str.ReplaceAll(MString::Format("%s.fEnergy", fNameMcEvt.Data()), "(x)");
 }
 
@@ -216 +217 @@
 TString MMcSpectrumWeight::GetFormulaSpecOld(const char *name) const
 {
-    return Form("pow(%s.fEnergy, %.3f)", name, fOldSlope);
+    return MString::Format("pow(%s.fEnergy, %.3f)", name, fOldSlope);
 }
 
@@ -231 +232 @@
 TString MMcSpectrumWeight::GetFormulaSpecNew(const char *name) const
 {
-    TString str = fFormula.IsNull() ? Form("pow(%s.fEnergy, %.3f)", name, fNewSlope) : fFormula.Data();
+    TString str = fFormula.IsNull() ? MString::Format("pow(%s.fEnergy, %.3f)", name, fNewSlope) : fFormula;
     if (!fFormula.IsNull())
-        str.ReplaceAll("X", Form("(%s.fEnergy)", name));
+        str.ReplaceAll("X", MString::Format("(%s.fEnergy)", name));
 
     return str;
@@ -265 +266 @@
 {
     if (GetFormulaSpecOld()==GetFormulaSpecNew())
-        return Form("%.16e", fNorm);
+        return MString::Format("%.16e", fNorm);
 
     const Double_t iold = fNormEnergy<0 ? GetSpecOldIntegral() : CalcSpecOld(fNormEnergy);
@@ -272 +273 @@
     const Double_t norm = fNorm*iold/inew;
 
-    return Form("%.16e*(%s)/(%s)", norm, GetFormulaSpecNew(name).Data(), GetFormulaSpecOld(name).Data());
+    return MString::Format("%.16e*(%s)/(%s)", norm, GetFormulaSpecNew(name).Data(), GetFormulaSpecOld(name).Data());
 }
 
@@ -388 +389 @@
     }
 
+    if (fNormEnergy<0 && fEnergyMin>0 && TMath::Abs(fEnergyMin-rh.GetELowLim())>1e-10)
+    {
+        *fLog << err;
+        *fLog << "You try to use changing minimum simulated Monte Carlo energies" << endl;
+        *fLog << "together with a normalization calculated from the integral." << endl;
+        *fLog << "This is not yet supported. Please switch to a normalization" << endl;
+        *fLog << "at a dedicated energy by specifying the energy" << endl;
+        *fLog << "    MMcSpectrumWeight.NormEnergy: 500" << endl;
+        *fLog << "in your sponde.rc." << endl;
+        return kFALSE;
+    }
+
     fOldSlope  = rh.GetSlopeSpec();
     fEnergyMin = rh.GetELowLim();
@@ -413 +426 @@
 
     fFunc = new TF1("", GetFormulaWeightsX().Data());
+    fFunc->SetName("SpectralWeighs");
     gROOT->GetListOfFunctions()->Remove(fFunc);
-    fFunc->SetName("SpectralWeighs");
 
     return kTRUE;
@@ -465 +478 @@
     const TAxis &axey = *h.GetYaxis();
 
-    // Find energy range between the minimum energy to be filles (emin)
+    // Find energy range between the minimum energy to be filled (emin)
     // and the minimum energy corresponding to the data filled into
     // this histogram (fEnergyMin)
     const Int_t first = axey.FindFixBin(emin);
-    const Int_t last  = axey.FindFixBin(fEnergyMin);
-    const Int_t max   = axey.FindFixBin(fEnergyMax);
+    const Int_t last  = axey.FindFixBin(fEnergyMin); // data range min energy
+    const Int_t max   = axey.FindFixBin(fEnergyMax); // data range max energy
 
     for (int x=1; x<=h.GetNbinsX(); x++)
@@ -537 +550 @@
     }
     *fLog << " Additional user norm.:    " << fNorm << endl;
-    *fLog << " Spectra are normalized:   " << (fNormEnergy<0?"by integral":Form("at %.1fGeV", fNormEnergy)) << endl;
+    *fLog << " Spectra are normalized:   " << (fNormEnergy<0?"by integral":MString::Format("at %.1fGeV", fNormEnergy)) << endl;
     if (hasold)
     {