Changeset 18179 for branches/Corsika7405Compatibility/mcorsika/MCorsikaRead.cc

Timestamp:

05/08/15 12:32:23 (9 years ago)

Author:

dbaack

Message:

 

File:

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

branches/Corsika7405Compatibility/mcorsika/MCorsikaRead.cc

@@ -290 +290 @@
                   status = fInFormat->Read(buffer, 272 * sizeof(Float_t));
                   status = fRunHeader->ReadEventHeader(buffer);
+                  //std::cout << "\x1b[31;1m EVTH Header: \x1b[0m";
+                  //fRunHeader->Print(nullptr);
                   break;
                   }
@@ -329 +331 @@
 Int_t MCorsikaRead::ReadTelescopePosition()
 {
-   if (!fInFormat->Read(&fNumTelescopes, 4)) return kERROR;
+        if (!fInFormat->Read(&fNumTelescopes, 4)) return kERROR;
 
    if (fTelescopeIdx>=fNumTelescopes)
@@ -474 +476 @@
 Int_t MCorsikaRead::Process()
 {
+        //std::cout << "Process \n\n" << std::endl;
    while (1)  // loop for multiple input files
    {
+           //std::cout << "OuterLoop \n\n" << std::endl;
+
       if (fInFormat)
       {
@@ -492 +497 @@
 
          Int_t status = kTRUE;
+
+         //std::cout << "||||ReadState: " << fReadState
+        //               << "  Block: " << fBlockType << "\n\n" << std::endl;
+
          switch (fBlockType)
             {
@@ -500 +509 @@
 
             case 1201:       // telescope position
+            //  std::cout << "TelescopPos ||||||" << std::endl;
                status = ReadTelescopePosition();
                break;
@@ -505 +515 @@
             case 1202:  // the event header
                Float_t buffer[272];
+               static int asd = 0;
+              // std::cout << "EventHeader |||||| " << ++asd << std::endl;
                if (!fInFormat->Read(buffer, 272 * sizeof(Float_t))) 
                   return kFALSE;
@@ -511 +523 @@
                   {
                   fRunHeader->ReadEventHeader(buffer);
-                  fRunHeader->Print();
+                  //fRunHeader->Print();
                   }
 
                status = fEvtHeader->ReadEvt(buffer);
+
+               //fEvtHeader->Print();
+               //std::cout << "TotReuse: " << fEvtHeader->GetTotReuse();
+
                if (fArrayIdx >= (Int_t)fEvtHeader->GetTotReuse())
                   {
@@ -527 +543 @@
 
             case 1204: // top level block for one array (only for eventio data)
+            //  std::cout << "EventIO1 |||||||" << std::endl;
                if (fArrayIdx < 0 || fArrayIdx == fBlockIdentifier)
                   { 
@@ -535 +552 @@
                   // skip this array of telescopes
                   fInFormat->Seek(fBlockLength);
-
+                   std::cout << "Skip Array!" << std::endl;
                break;
 
@@ -541 +558 @@
             case 1205: // eventio data
                {
+              // std::cout << "EventIO2 |||||||" << std::endl;
                Int_t telIdx   = fBlockIdentifier % 1000;
-               if (fBlockVersion == 0                               &&
-                     (fTelescopeIdx < 0 || fTelescopeIdx ==  telIdx)     )
-                  {
+           //    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;
@@ -560 +582 @@
                else
                   // skip this telescope
+               {
                   fInFormat->Seek(fBlockLength);
+
+                   std::cout << "Skip this telescope!" << std::endl;
+               }
                }
                break;
 
             case 1209:  // the event end
+           //   std::cout << "EventEnd ||||||" << std::endl;
                status = fEvtHeader->ReadEvtEnd(fInFormat);
                
@@ -581 +608 @@
 
             case 1210:  // the run end
+                //std::cout << "RunEnd ||||||" << std::endl;
                status = fRunHeader->ReadEvtEnd(fInFormat, kTRUE);
                fNumEvents += fRunHeader->GetNumEvents();
@@ -594 +622 @@
 
             case 1105:  // event block of raw format
+                static int NmrEventBlock = 0;
+                //std::cout << "EventBlock ||||||" << ++NmrEventBlock << std::endl;
                if (fReadState == 2 || fReadState == 10)
                   {
@@ -607 +637 @@
 
             case 1109:  // save corsika events
+                //std::cout << "SaveCorsika ||||||" << std::endl;
                fEvtHeader->InitXY();
                status = fEvent->ReadCorsikaEvt(&fRawEvemtBuffer[0], 
@@ -630 +661 @@
             default:
                // unknown block, we skip it
+               std::cout << "UnknownBlock " << fBlockType << std::endl;
                fReadState = 15;
                fInFormat->Seek(fBlockLength);
 
             }
+
+         std::cout << "Status: " << status << std::endl;
 
          if (status != kTRUE)
@@ -640 +674 @@
 
          Int_t headerStatus =  ReadNextBlockHeader();
+
+         std::cout << "Header Status: " << headerStatus << std::endl;
+
+
          if (headerStatus == kFALSE) 
             // end of file