Changeset 1729

Timestamp:

01/23/03 17:25:53 (22 years ago)

Author:

rwagner

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• manalysis/MSigmabarCalc.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -1 +1 @@
 
                                                          -*-*- END -*-*-
+
+ 2003/01/23: Robert Wagner
+
+  * manalyis/MSigmabarCalc.cc
+    - MSigmabarCalc::Process()
+      Replaced fMcEvt->GetTheta() by fMcEvt->GetThetaTelescope()
+  
+
  2003/01/20: Oscar Blanch
   
   * mmc/MMcEvt.[hxx,cxx]
-    - Variable members: fElecCphFraction, fMuonCphFraction, 
+    - Data members: fElecCphFraction, fMuonCphFraction, 
       fOtherCphFraction have been introduced.
     - Class version updated to 3.
 
   * mmc/MMcCorsikaRunHeader.[h.cc]
-    - Varible members: fWobbleMode and fAtmosphericModel introduced.
-    - Clas version updated to 2.
+    - Data members: fWobbleMode and fAtmosphericModel introduced.
+    - Class version updated to 2.
+
+
 
  2003/01/19: Thomas Bretz

trunk/MagicSoft/Mars/manalysis/MSigmabarCalc.cc

@@ -16 +16 @@
 !
 !
-!   Author(s): Robert Wagner <magicdev@rwagner.de> 10/2002
+!   Author(s): Robert Wagner <rwagner@mppmu.mpg.de> 10/2002
 !
-!   Copyright: MAGIC Software Development, 2000-2002
+!   Copyright: MAGIC Software Development, 2000-2003
 !
 !
@@ -31 +31 @@
 //  in MSigmabarParam. For the time being, Theta is taken from a Monte     //
 //  Carlo container. This is preliminary and to be changed ASAP.           //
- //                                                                         //
+//                                                                         //
 /////////////////////////////////////////////////////////////////////////////
 #include "MSigmabarCalc.h"
@@ -98 +98 @@
         return kFALSE;
     }
+
+    fRun = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
+    if (!fRun)
+    {
+        *fLog << dbginf << "MRawRunHeader not found... aborting." << endl;
+        return kFALSE;
+    }
     
     fMcEvt = (MMcEvt*)pList->FindObject("MMcEvt");
@@ -123 +130 @@
   fSigmabarMin = TMath::Min((Double_t)fSig->GetSigmabar(), fSigmabarMin);
 
-  if (fMcEvt->GetTheta()*kRad2Deg < 90)
-    fThetaMax    = TMath::Max(fMcEvt->GetTheta()*kRad2Deg, fThetaMax);
-  fThetaMin    = TMath::Min(fMcEvt->GetTheta()*kRad2Deg, fThetaMin);
+  if (fMcEvt->GetThetaTelescope()*kRad2Deg < 90)
+    fThetaMax    = TMath::Max(fMcEvt->GetThetaTelescope()*kRad2Deg, fThetaMax);
+  fThetaMin    = TMath::Min(fMcEvt->GetThetaTelescope()*kRad2Deg, fThetaMin);
 
   return rc;
@@ -138 +145 @@
 {
    
-  fSigParam->SetParams(1, fSigmabarMin, fSigmabarMax, fThetaMin, fThetaMax);  
+  fSigParam->SetParams(1, fSigmabarMin, fSigmabarMax, fThetaMin, fThetaMax); 
+  fSigParam->SetRunNumber(fRun->GetRunNumber());
+ 
   Reset();