Changeset 9359 for trunk/MagicSoft/Mars/msim

Timestamp:

02/20/09 00:39:31 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msim

Files:

• MSimMMCS.cc (modified) (4 diffs)
• MSimMMCS.h (modified) (2 diffs)

trunk/MagicSoft/Mars/msim/MSimMMCS.cc

@@ -67 +67 @@
 Int_t MSimMMCS::PreProcess(MParList *plist)
 {
-    if (!plist->FindCreateObj("MMcRunHeader"))
+    fMcRunHeader = (MMcRunHeader*)plist->FindCreateObj("MMcRunHeader");
+    if (!fMcRunHeader)
         return kFALSE;
 
@@ -105 +106 @@
     }
 
-
     return kTRUE;
 }
@@ -126 +126 @@
     // ----------------------------------------------------
 
-    MMcRunHeader *mrh = (MMcRunHeader*)plist->FindObject("MMcRunHeader");
-    if (!mrh)
-    {
-        *fLog << err << "MMcRunHeader not found... aborting." << endl;
-        return kFALSE;
-    }
-
     //    fNumPheFromDNSB        MMcPedestalNSBAdd   // Number of phe/ns from diffuse NSB
 
     // FIXME: Is there a way to write them as LAST entry in the file?
-    mrh->SetNumSimulatedShowers(fRunHeader->GetNumEvents());
-    mrh->SetImpactMax(fRunHeader->GetImpactMax());
-    mrh->SetCorsikaVersion(TMath::Nint(fRunHeader->GetProgramVersion()*100));
+    fMcRunHeader->SetNumSimulatedShowers(fRunHeader->GetNumEvents());
+    fMcRunHeader->SetCorsikaVersion(TMath::Nint(fRunHeader->GetProgramVersion()*100));
+
+    if (fRunHeader->GetImpactMax()>0)
+        fMcRunHeader->SetImpactMax(fRunHeader->GetImpactMax());
 
     // ----------------------------------------------------
@@ -185 +180 @@
     fMcEvt->SetPhotElfromShower(0);
 
+    if (fRunHeader->GetImpactMax()<0 &&
+        fEvtHeader->GetImpact()>fMcRunHeader->GetImpactMax())
+        fMcRunHeader->SetImpactMax(fEvtHeader->GetImpact());
+
     return kTRUE;
 }

trunk/MagicSoft/Mars/msim/MSimMMCS.h

@@ -12 +12 @@
 class MCorsikaEvtHeader;
 class MCorsikaRunHeader;
+class MMcRunHeader;
 
 class MSimMMCS : public MTask
@@ -21 +22 @@
     MCorsikaEvtHeader *fEvtHeader;
     MCorsikaRunHeader *fRunHeader;
+    MMcRunHeader      *fMcRunHeader;
 
     // MTask