Changeset 1729 for trunk/MagicSoft/Mars/manalysis/MSigmabarCalc.cc

Timestamp:

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

Author:

rwagner

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MSigmabarCalc.cc (modified) (5 diffs)

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();