Changeset 18569

Timestamp:

09/03/16 20:40:27 (8 years ago)

Author:

tbretz

Message:

Implemented function SetupFits to be able to read from a fits file.

Location:

trunk/Mars/mcorsika

Files:

• MCorsikaEvtHeader.cc (modified) (2 diffs)
• MCorsikaEvtHeader.h (modified) (1 diff)

trunk/Mars/mcorsika/MCorsikaEvtHeader.cc

@@ -44 +44 @@
 #include <fstream>
 
+#include "fits.h"
+
 #include <TMath.h>
 
@@ -196 +198 @@
     return fInFormat->Eof() ? kERROR : kTRUE;
 }
+
+Bool_t MCorsikaEvtHeader::SetupFits(fits &fin)
+{
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fEvtNumber",              fEvtNumber              ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fNumReuse",               fNumReuse               ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fTotalEnergy",            fTotalEnergy            ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fStartAltitude",          fStartAltitude          ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fFirstTargetNum",         fFirstTargetNum         ))  return kFALSE;
+    if (!fin.SetRefAddress( "CorsikaEvtHeader.fFirstInteractionHeight", fFirstInteractionHeight ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fMomentumX",              fMomentumX              ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fMomentumY",              fMomentumY              ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fMomentumZ",              fMomentumZ              ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fZd",                     fZd                     ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fAz",                     fAz                     ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fX",                      fX                      ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fY",                      fY                      ))  return kFALSE;
+    if (!fin.SetRefAddress("MCorsikaEvtHeader.fWeightedNumPhotons",     fWeightedNumPhotons     ))  return kFALSE;
+
+    return kTRUE;
+}

trunk/Mars/mcorsika/MCorsikaEvtHeader.h

@@ -83 +83 @@
     Int_t ReadEvtEnd(MCorsikaFormat *informat); // read in event end block
 
+    Bool_t SetupFits(fits &fin);
+
     ClassDef(MCorsikaEvtHeader, 3) // Parameter Conatiner for raw EVENT HEADER
 };