Index: trunk/MagicSoft/Mars/mhist/MHMcRate.cc
===================================================================
--- trunk/MagicSoft/Mars/mhist/MHMcRate.cc	(revision 1108)
+++ trunk/MagicSoft/Mars/mhist/MHMcRate.cc	(revision 1164)
@@ -129,4 +129,24 @@
 // --------------------------------------------------------------------------
 //
+//  set the parameters to compute the incident rate 
+//
+void MHMcRate:: SetFlux(Float_t flux0, Float_t specindx)
+{
+    fFlux0=flux0;
+    fSpecIndex=specindx;
+
+}
+
+// --------------------------------------------------------------------------
+//
+//  set the incident rate 
+//
+void MHMcRate:: SetIncidentRate(Float_t showerrate)
+{
+    fShowerRate=showerrate;
+}
+
+// --------------------------------------------------------------------------
+//
 //  update the limits for energy, theta, phi and impact parameter
 //
@@ -163,5 +183,6 @@
     const Double_t epowmin = pow(fEnergyMin, specidx);
 
-    fShowerRate = fFlux0/specidx*(epowmax-epowmin);
+    if (fShowerRate <= 0)
+	fShowerRate = fFlux0/specidx*(epowmin-epowmax);
 
     if (fPartId!=1)