Changeset 18569 for trunk/Mars/mcorsika/MCorsikaEvtHeader.cc

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.

File:

• trunk/Mars/mcorsika/MCorsikaEvtHeader.cc (modified) (2 diffs)

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;
+}