Changeset 18180 for branches/Corsika7405Compatibility/mcorsika

Timestamp:

05/12/15 19:33:34 (10 years ago)

Author:

dbaack

Message:

Fixed bug with non cleared "PhotonEvent"

File:

• branches/Corsika7405Compatibility/mcorsika/MCorsikaRead.cc (modified) (13 diffs)

branches/Corsika7405Compatibility/mcorsika/MCorsikaRead.cc

@@ -54 +54 @@
 #include "MCorsikaEvtHeader.h"
 
+#include "../msim/MPhotonData.h"
+
 #include "MPhotonEvent.h"
 
@@ -476 +478 @@
 Int_t MCorsikaRead::Process()
 {
-        //std::cout << "Process \n\n" << std::endl;
+        fEvent->Clear();
+        fEvent->Reset();
+        fEvent->Resize(0);
+
+
+        std::cout << "\n\n MCorsikaRead Process " << std::endl;
    while (1)  // loop for multiple input files
    {
-           //std::cout << "OuterLoop \n\n" << std::endl;
 
       if (fInFormat)
@@ -498 +504 @@
          Int_t status = kTRUE;
 
-         //std::cout << "||||ReadState: " << fReadState
-        //               << "  Block: " << fBlockType << "\n\n" << std::endl;
-
          switch (fBlockType)
             {
@@ -516 +519 @@
                Float_t buffer[272];
                static int asd = 0;
-              // std::cout << "EventHeader |||||| " << ++asd << std::endl;
                if (!fInFormat->Read(buffer, 272 * sizeof(Float_t))) 
                   return kFALSE;
@@ -528 +530 @@
                status = fEvtHeader->ReadEvt(buffer);
 
-               //fEvtHeader->Print();
-               //std::cout << "TotReuse: " << fEvtHeader->GetTotReuse();
 
                if (fArrayIdx >= (Int_t)fEvtHeader->GetTotReuse())
@@ -543 +543 @@
 
             case 1204: // top level block for one array (only for eventio data)
-            //  std::cout << "EventIO1 |||||||" << std::endl;
-               if (fArrayIdx < 0 || fArrayIdx == fBlockIdentifier)
+              if (fArrayIdx < 0 || fArrayIdx == fBlockIdentifier)
                   { 
                   fReadState = 4;
@@ -558 +557 @@
             case 1205: // eventio data
                {
-              // std::cout << "EventIO2 |||||||" << std::endl;
                Int_t telIdx   = fBlockIdentifier % 1000;
-           //    std::cout << "TelIDX: " << telIdx << " fTelescopeIdx:" << fTelescopeIdx
-            //             << " fBlockVersion" << fBlockVersion << std::endl;
 
                if (fBlockVersion == 0  && (fTelescopeIdx < 0 || fTelescopeIdx ==  telIdx) )
                {
                   status = fEvent->ReadEventIoEvt(fInFormat);
-
-                  //std::cout << " ArrayIdx:" << fBlockIdentifier / 1000.0 << std::endl;
 
                   Int_t arrayIdx = fBlockIdentifier / 1000;
@@ -579 +573 @@
                   fEvent->SimWavelength(fRunHeader->GetWavelengthMin(), 
                                           fRunHeader->GetWavelengthMax());
+
+
+                  bool tmpWeight = false;
+                  std::cout << "ReadEventIO Telescope data with " << fEvent->GetNumPhotons() << std::endl;
+                  for(int i=0; i<fEvent->GetNumPhotons(); i++)
+                  {
+
+                          const MPhotonData &ph = (*fEvent)[i];
+                          if(ph.GetWeight() != 1)
+                                  tmpWeight = true;
+                  }
+                  if(tmpWeight)
+                          std::cout << "CorsikaRead Photon weight != 1" << std::endl;
+
                   }
                else
@@ -591 +599 @@
 
             case 1209:  // the event end
-           //   std::cout << "EventEnd ||||||" << std::endl;
                status = fEvtHeader->ReadEvtEnd(fInFormat);
                
@@ -608 +615 @@
 
             case 1210:  // the run end
-                //std::cout << "RunEnd ||||||" << std::endl;
-               status = fRunHeader->ReadEvtEnd(fInFormat, kTRUE);
+                status = fRunHeader->ReadEvtEnd(fInFormat, kTRUE);
                fNumEvents += fRunHeader->GetNumEvents();
                //fRunHeader->SetReadyToSave();
@@ -637 +643 @@
 
             case 1109:  // save corsika events
-                //std::cout << "SaveCorsika ||||||" << std::endl;
+                std::cout << "Block1109" << std::endl;
                fEvtHeader->InitXY();
                status = fEvent->ReadCorsikaEvt(&fRawEvemtBuffer[0], 
@@ -667 +673 @@
             }
 
-         std::cout << "Status: " << status << std::endl;
-
          if (status != kTRUE)
             // there was an error while data were read
@@ -674 +678 @@
 
          Int_t headerStatus =  ReadNextBlockHeader();
-
-         std::cout << "Header Status: " << headerStatus << std::endl;
-
 
          if (headerStatus == kFALSE)