Changeset 1783 for trunk/MagicSoft/Mars/mhist/MHMcRate.cc

Timestamp:

02/21/03 18:09:36 (22 years ago)

Author:

moralejo

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHMcRate.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mhist/MHMcRate.cc

@@ -18 +18 @@
 !   Author(s): Thomas Bretz  12/2000 <mailto:tbretz@uni-sw.gwdg.de>
 !   Author(s): Harald Kornmayer 1/2001
+!   Author(s): Abelardo Moralejo 2/2003
+!
+!   Explanations on the rate calculation can be found in 
+!   chapter 7 of the following diploma thesis: 
+!   http://www.pd.infn.it/magic/tesi2.ps.gz (in Italian)
 !
 !   Copyright: MAGIC Software Development, 2000-2001
@@ -38 +43 @@
     fPartId=0;               // Type of particle
 
-    fEnergyMax=0.0;          // Maximum Energy in GeV
-    fEnergyMin=1000000.0;    // Minimum Energy in GeV
-
+    fEnergyMax=0.0;          // Maximum Energy (TeV)
+    fEnergyMin=1000000.0;    // Minimum Energy (TeV)
+
+    fSolidAngle = -1.;       // Solid angle within which incident directions
+                             // are distributed
     fThetaMax=0.0;           // Maximum theta angle of run
     fThetaMin=370.0;         // Minimum theta angle of run
@@ -119 +126 @@
 // --------------------------------------------------------------------------
 //
-//  Set the information about trigger due only to the background conditions
+//  Set the information about trigger due only to the Night Sky Background:
 //
 void MHMcRate::SetBackground (Float_t showers, Float_t triggers)
@@ -184 +191 @@
 
     if (fShowerRate <= 0)
-        fShowerRate = fFlux0/specidx*(epowmin-epowmax);
+      fShowerRate = fFlux0/specidx*(epowmax-epowmin);
+
+    if (fSolidAngle < 0.)
+      fSolidAngle = (fPhiMax-fPhiMin)*(cos(fThetaMin)-cos(fThetaMax));
 
     if (fPartId!=1)
-        fShowerRate *= (fPhiMax-fPhiMin)*(cos(fThetaMax)-cos(fThetaMin));
-
-    const Double_t impactdiff = fImpactMax-fImpactMin;
-
-    fShowerRate *= TMath::Pi()*(impactdiff/100.0*impactdiff/100.0);
+      fShowerRate *= fSolidAngle;
+
+    fShowerRate *= TMath::Pi()*(fImpactMax/100.0*fImpactMax/100.0 - 
+                                fImpactMin/100.0*fImpactMin/100.0);
 
     fShowerRateError = sqrt(fShowerRate);