Changeset 1783 for trunk/MagicSoft/Mars/mhist

Timestamp:

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

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

• MHMcRate.cc (modified) (4 diffs)
• MHMcRate.h (modified) (3 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);

trunk/MagicSoft/Mars/mhist/MHMcRate.h

@@ -12 +12 @@
     UShort_t fPartId;           // Type of particle
 
-    Float_t fEnergyMax;         // Maximum Energy in GeV
-    Float_t fEnergyMin;         // Minimum Energy in GeV
+    Float_t fEnergyMax;         // Maximum Energy in TeV
+    Float_t fEnergyMin;         // Minimum Energy in TeV
 
     Float_t fThetaMax;          // Maximum theta angle of run
@@ -20 +20 @@
     Float_t fPhiMin;            // Minimum phi angle of run
 
-    Float_t fImpactMax;         // Maximum impact parameter
-    Float_t fImpactMin;         // Minimum impact parameter
+    Float_t fSolidAngle;        // Solid angle within which incident directions
+                                // are distributed (sr)
+
+    Float_t fImpactMax;         // Maximum impact parameter (cm)
+    Float_t fImpactMin;         // Minimum impact parameter (cm)
 
     Float_t fBackTrig;          // Number of triggers from background
@@ -50 +53 @@
     void SetIncidentRate(Float_t showerrate);
 
+    void SetImpactMax(Float_t Impact) {fImpactMax=Impact;}
+    void SetImpactMin(Float_t Impact) {fImpactMin=Impact;}
+
+    void SetThetaMax(Float_t Theta) {fThetaMax=Theta;}
+    void SetThetaMin(Float_t Theta) {fThetaMin=Theta;}
+    void SetPhiMax(Float_t Phi) {fPhiMax=Phi;}
+    void SetPhiMin(Float_t Phi) {fPhiMin=Phi;}
+
+    void SetSolidAngle(Float_t Solid) {fSolidAngle=Solid;}
+    void SetEnergyMax(Float_t Energy) {fEnergyMax=Energy;}
+    void SetEnergyMin(Float_t Energy) {fEnergyMin=Energy;}
+
     void UpdateBoundaries(Float_t energy, Float_t theta, Float_t phi, Float_t impact);