Changeset 18447 for trunk/Mars/msim/MPhotonEvent.cc

Timestamp:

03/08/16 13:11:42 (9 years ago)

Author:

ftemme

Message:

Implemented the method AsciiWrite for MPhotonEvent

File:

• trunk/Mars/msim/MPhotonEvent.cc (modified) (1 diff)

trunk/Mars/msim/MPhotonEvent.cc

@@ -307 +307 @@
 }
 
+// Overload the AsciiWrite method to store the informations of the photons onto disc.
+// For each corsika/pedestal... event there is a simple header written. After this
+// the informations for each photon are written in one line.
+Bool_t MPhotonEvent::AsciiWrite(ostream &out) const
+{
+    // Write a simple header. Be aware there will be written one header for each
+    // corsika/pedestal... event. Each photon of the event will be written in one
+    // line.
+    out << "#";
+    out << "photonID" << ",";
+    out << "tag" << ",";
+    out << "mirrorTag" << ",";
+    out << "posX" << ",";
+    out << "posY" << ",";
+    out << "cosU" << ",";
+    out << "cosV" << ",";
+    out << "time" << ",";
+    out << "weight" << ",";
+    out << "wavelength" << ",";
+    out << "productionHeight" << ",";
+    out << "primaryID" << endl;
+
+    // Get number of photons
+    const Int_t num = GetNumPhotons();
+
+    // Loop over all photons
+    for (Int_t i=0; i<num; i++)
+    {
+        // Get i-th photon
+        const MPhotonData &ph = operator[](i);
+
+        out << i << "," ;
+        out << ph.GetTag() << "," ;
+        out << ph.GetMirrorTag() << "," ;
+        out << ph.GetPosX() << "," ;
+        out << ph.GetPosY() << "," ;
+        out << ph.GetCosU() << "," ;
+        out << ph.GetCosV() << "," ;
+        out << ph.GetTime() << "," ;
+        out << ph.GetWeight() << "," ;
+        out << ph.GetWavelength() << "," ;
+        out << ph.GetProductionHeight() << "," ;
+        out << ph.GetPrimary() << endl ;
+    }
+    return kTRUE;
+}
+
 // --------------------------------------------------------------------------
 //