Index: trunk/MagicSoft/Mars/manalysis/MSigmabarCalc.cc
===================================================================
--- trunk/MagicSoft/Mars/manalysis/MSigmabarCalc.cc	(revision 1715)
+++ trunk/MagicSoft/Mars/manalysis/MSigmabarCalc.cc	(revision 1729)
@@ -16,7 +16,7 @@
 !
 !
-!   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,5 +31,5 @@
 //  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,4 +98,11 @@
         return kFALSE;
     }
+
+    fRun = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
+    if (!fRun)
+    {
+        *fLog << dbginf << "MRawRunHeader not found... aborting." << endl;
+        return kFALSE;
+    }
     
     fMcEvt = (MMcEvt*)pList->FindObject("MMcEvt");
@@ -123,7 +130,7 @@
   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,5 +145,7 @@
 {
    
-  fSigParam->SetParams(1, fSigmabarMin, fSigmabarMax, fThetaMin, fThetaMax);  
+  fSigParam->SetParams(1, fSigmabarMin, fSigmabarMax, fThetaMin, fThetaMax); 
+  fSigParam->SetRunNumber(fRun->GetRunNumber());
+ 
   Reset();
   
