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

Timestamp:

06/18/15 15:22:04 (9 years ago)

Author:

dbaack

Message:

Fixed bug that creates empty or duplicated events for CSCAT=1 in Corsika

File:

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

branches/Corsika7405Compatibility/mcorsika/MCorsikaRead.cc

@@ -406 +406 @@
       return kTRUE;
 
-   if (fBlockType == 1204 && fReadState != 2)
+   if (fBlockType == 1204 && (fReadState != 2 && fReadState != 15) )
       // next is a new set of telescope arrays, we store the previous ones
       // but not if this is the first one (fReadState != 2)
@@ -478 +478 @@
 Int_t MCorsikaRead::Process()
 {
-        fEvent->Clear();
-        fEvent->Reset();
+        //fEvent->Clear();
+
+        //fEvent->Resize(0);
+
+
+        std::cout << "\n\n MCorsikaRead Process " << std::endl;
         fEvent->Resize(0);
-
-
-        std::cout << "\n\n MCorsikaRead Process " << std::endl;
    while (1)  // loop for multiple input files
    {
@@ -530 +531 @@
                status = fEvtHeader->ReadEvt(buffer);
 
+               std::cout << "----------ArrayIndex:: " << fArrayIdx << std::endl;
 
                if (fArrayIdx >= (Int_t)fEvtHeader->GetTotReuse())
@@ -550 +552 @@
                else
                   // skip this array of telescopes
+                   {
                   fInFormat->Seek(fBlockLength);
                    std::cout << "Skip Array!" << std::endl;
+                   }
                break;
 
@@ -558 +562 @@
                {
                Int_t telIdx   = fBlockIdentifier % 1000;
+
+               std::cout << "----------TelescopeIndex:: " << fTelescopeIdx << std::endl;
 
                if (fBlockVersion == 0  && (fTelescopeIdx < 0 || fTelescopeIdx ==  telIdx) )
@@ -575 +581 @@
 
 
-                  bool tmpWeight = false;
+                  /*bool tmpWeight = false;
                   std::cout << "ReadEventIO Telescope data with " << fEvent->GetNumPhotons() << std::endl;
                   for(int i=0; i<fEvent->GetNumPhotons(); i++)
@@ -585 +591 @@
                   }
                   if(tmpWeight)
-                          std::cout << "CorsikaRead Photon weight != 1" << std::endl;
+                          std::cout << "CorsikaRead Photon weight != 1" << std::endl;*/
 
                   }
@@ -629 +635 @@
             case 1105:  // event block of raw format
                 static int NmrEventBlock = 0;
-                //std::cout << "EventBlock ||||||" << ++NmrEventBlock << std::endl;
+                std::cout << "EventBlock ||||||" << ++NmrEventBlock << std::endl;
                if (fReadState == 2 || fReadState == 10)
                   {
@@ -677 +683 @@
             return status;
 
+
          Int_t headerStatus =  ReadNextBlockHeader();
+
+         std::cout << "fBlockType: " << fBlockType << "  fReadState: " << fReadState
+                         << " HeaderStatus: " << headerStatus << std::endl;
+
 
          if (headerStatus == kFALSE)