offsetadd.cxx

Go to the documentation of this file.

// g++ offsetadd.cxx `root-config --cflags --libs` -I${GRSISYS}/include -L${GRSISYS}/libraries -lMidasFormat -lXMLParser
// -ooffsetadd
 
#include "TMidasFile.h"
#include "TMidasEvent.h"
 
#include <iostream>
 
// void WriteEventToFile(TMidasFile*,std::shared_ptr<TMidasEvent>,Option_t);
 
void ProcessEvent(const std::shared_ptr<TMidasEvent>&, TMidasFile*);
 
int main(int argc, char** argv)
{
 
   if(argc != 3) {
      printf("Usage: ./offsetadd <input.mid> <output.mid>\n");
      return 1;
   }
 
   auto*                        infile  = new TMidasFile;
   auto*                        outfile = new TMidasFile;
   std::shared_ptr<TMidasEvent> event   = std::make_shared<TMidasEvent>();
 
   if(argv[1] == argv[2]) {
      printf("ERROR: Cannot overwrite midas file %s\n", argv[1]);
   }
 
   infile->Open(argv[1]);
   outfile->OutOpen(argv[2]);
 
   while(infile->Read(event) > 0) {
      switch(event->GetEventId()) {
      case 0x8000:
         printf("start of run\n");
         outfile->Write(event);
         event->Print();
         break;
      case 0x8001:
         printf("end of run\n");
         event->Print();
         outfile->Write(event);
         break;
      default: ProcessEvent(event, outfile); break;
      };
   }
   infile->Close();
   outfile->Close();
   return 0;
}
 
// void WriteEventToFile(TMidasFile *outfile,std::shared_ptr<TMidasEvent> event) {
//   outfile->Write(event);
//}
 
void ProcessEvent(const std::shared_ptr<TMidasEvent>& event, TMidasFile* outfile)
{
   if(event->GetEventId() != 1) {
      outfile->Write(event, "q");
      return;
   }
   event->SetBankList();
   // int size;
   // int data[1024];
 
   void* ptr      = nullptr;
   int   banksize = event->LocateBank(nullptr, "GRF1", &ptr);
 
   uint32_t type  = 0xffffffff;
   int      value = 0xffffffff;
 
   int dettype = 0;
   int chanadd = 0;
 
   int timelow  = 0;
   int timehigh = 0;
 
   int64_t time = 0;
 
   for(int x = 0; x < banksize; x++) {
      value = *(reinterpret_cast<int*>(ptr) + x);
      type  = value & 0xf0000000;
 
      switch(type) {
      case 0x80000000:
         dettype = value & 0x0000000f;
         chanadd = (value & 0x0003fff0) >> 4;
         break;
      case 0xa0000000: timelow = value & 0x0fffffff; break;
      case 0xb0000000: timehigh = value & 0x00003fff; break;
      };
   }
 
   // printf("chanadd = 0x%08x
   // event->Print("a");
 
   // printf("dettype  = 0x%08x\n",dettype);
   // printf("chanadd  = 0x%08x\n",chanadd);
 
   // printf("timelow  = 0x%08x\n",timelow);
   // printf("timehigh = 0x%08x\n",timehigh);
 
   if((dettype == 1) || (dettype == 5)) {   // 1 for GRIFFIN, 5 for PACES
                                            // do nothing.
   } else {
      outfile->Write(event, "q");
      return;
   }
 
   if(((chanadd & 0x0000ff00) == 0x00000000) || ((chanadd & 0x0000ff00) == 0x00001000) ||
      ((chanadd & 0x0000ff00) == 0x00001100) || ((chanadd & 0x0000ff00) == 0x00001200)) {
      // do nothing.
   } else {
      outfile->Write(event, "q");
      return;
   }
 
   time = timehigh;
   time = time << 28;
   time |= timelow & 0x0fffffff;
 
   // printf("time = 0x%016x\n",time);
   // std::cout<<"time    = "<<std::hex<<time<<std::endl;
 
   // Here's where we change the values of the time stamps!!!!
   switch(chanadd & 0x0000ff00) {
   case 0x00000000:   // if the first GRIF-16
                      //         time -= 10919355323; // run 2369 correction
      time -= 87;     // run 2394 correction
      break;
   case 0x00000100:   // if the second GRIF-16
      break;
   case 0x00001000:   // if the third GRIF-16
                      //         time -= 10919355323; // run 2369 correction
      break;
   case 0x00001100:   // if the fourth GRIF-16
                      //         time -= 10919355239; // run 2369 correction
      break;
   case 0x00001200:   // if the fifth GRIF-16
                      //         time += 7;
      time -= 87;     // run 2394 correction
      break;
   };
   if(time < 0) {
      time += 0x3ffffffffff;
   } else if(time > 0x3ffffffffff) {
      time -= 0x3ffffffffff;
   }
 
   // moving these inside the next switch, to account for doubly printed words.
   // (hey, it happens.)
   // -JKS, 14 January 2015
   // timelow = time&0x0fffffff;
   // timehigh = (time&0x3fff0000000) >> 28;
 
   //   printf(DRED);
   //   event->Print("a");
   //   printf(RESET_COLOR);
 
   std::shared_ptr<TMidasEvent> copyevent = std::make_shared<TMidasEvent>(*event);
   copyevent->SetBankList();
 
   banksize = copyevent->LocateBank(nullptr, "GRF1", &ptr);
   for(int x = 0; x < banksize; x++) {
      value = *(reinterpret_cast<int*>(ptr) + x);
      type  = value & 0xf0000000;
 
      switch(type) {
      case 0xa0000000:
         timelow = time & 0x0fffffff;
         timelow += 0xa0000000;
         *(reinterpret_cast<int*>(ptr) + x) = timelow;
         break;
      case 0xb0000000: {
         timehigh     = (time & 0x3fff0000000) >> 28;
         int tempdead = value & 0xffffc000;
         timehigh += tempdead;
         *(reinterpret_cast<int*>(ptr) + x) = timehigh;
         break;
      }
      };
 
      // printf( "0x%08x ",*((int*)ptr+x));
      // if(x!=0 && (x%7)==0)
      //   printf("\n");
   }
   // printf("===================\n");
 
   outfile->Write(copyevent, "q");
 
   //   printf(DBLUE);
   //   copyevent->Print("a");
   //   printf(RESET_COLOR);
}