analysis/ScEventBuilder.cxx (r4864/r595)

Go to the documentation of this file.

#include <iostream>
#include <stdio.h>

#include "nxyter/Data.h"


#include "LmdFile.h"


#include "TBranch.h"
#include "TClonesArray.h"
#include "TFile.h"
#include "TTree.h"


#include "ScDigi.h"
#include "ScEvent.h"
#include "ScEventBuilder.h"

using namespace std;


ScEventBuilder::ScEventBuilder() {
  fTlimit = 100.;
}


ScEventBuilder::~ScEventBuilder() {
}


Int_t ScEventBuilder::Start() {


  // Event counter
  Int_t nMbsEvents = 0;


  // Open input LMD file
  mbs::LmdFile inFile;
  if ( ! inFile.OpenRead(fInputFile.Data()) ) {
    cout << "--E-- EventBuilder: Could not open input file " << fInputFile << endl;
    return -1;
  }


  // Open output ROOT file
  TFile* outFile = new TFile(fOutputFile, "RECREATE");
  if ( outFile->IsZombie() ) {
    cout << "--E-- EventBuilder: Could not open output file " << fOutputFile << endl;
    return -1;
  }




  // Create output tree and branches
  ScEvent* event = new ScEvent();
  TTree* outTree = new TTree("cbmdata", "data tree");
  outTree->Branch("event", &event);



  // Create pointers to MBS event, subevent and data
  mbs::EventHeader*    mbsEvt    = NULL;
  mbs::SubeventHeader* mbsSubEvt = NULL;
  nxyter::Data*         data      = NULL;


  // Global counters
  Long64_t nTotEpochs = 0;
  Long64_t nTotHits   = 0;
  Long64_t nTotAux    = 0;
  Long64_t nTotEvents = 0;


  // Variables for transient storage
  Int_t    nData            = 0;
  Int_t    currentEpoch     = 0;
  Long64_t currentEpochTime = 0;
  Long64_t hitTime          = 0;



  // MBS event loop
  while ( (mbsEvt = inFile.ReadEvent()) ) {
    nMbsEvents++;
    mbsSubEvt         = NULL;
    Int_t nMbsSubEvts = 0;
    Int_t nHits       = 0;
    Int_t nAux        = 0;
    Int_t nEpochs     = 0;
    Int_t nEvents     = 0;


    // MBS subevent loop
    while ( (mbsSubEvt = mbsEvt->NextSubEvent(mbsSubEvt)) ) {
      nMbsSubEvts++;


      // Get  messages in subevent
      data = (nxyter::Data*) mbsSubEvt->RawData();
      nData = mbsSubEvt->RawDataSize() / 6;
      while (nData--) {


        // Treat epoch marker
        if ( data->isEpochMsg() ) {
          nEpochs++;
          currentEpoch     = data->getEpoch();
          currentEpochTime = Long64_t(currentEpoch) << 14;
        }

        // Treat hit message
        else if ( data->isHitMsg() || data->isAuxMsg() ) {

          // Calculate hit time from epoch and time stamp
          if ( data->isHitMsg() ) {
            hitTime = currentEpochTime + Long64_t(data->getNxTs());
            nHits++;
          }
          else {
            hitTime = currentEpochTime + Long64_t(data->getAuxTs());
            nAux++;
          }

          // Add hit if event is empty or hit time is within limits
          if ( event->IsEmpty() || (Double_t(hitTime) - event->GetTime() < fTlimit) ) {
            event->AddData(data, hitTime, currentEpoch);
          }

          // Else create new event and add hit to that one
          else {
            outTree->Fill();
            event->Clear();
            nEvents++;
            event->AddData(data, hitTime, currentEpoch);
          }

        }              // Message type
        data++;        // Next message

      }                // Loop over messages

    }                  // MBS subevent loop


    cout << "MBS Event " << nMbsEvents << " Subevents " << nMbsSubEvts << " Epochs "
         << nEpochs << " Hits " << nHits << " Aux " << nAux << " Events "
         << nEvents << endl;

    nTotEpochs += nEpochs;
    nTotHits   += nHits;
    nTotAux    += nAux;
    nTotEvents += nEvents;


  }            // MBS event loop




  inFile.Close();
  outFile->Write();
  outFile->Close();

  cout << endl << endl;
  cout << "===============================================================" << endl;
  cout << "Total MBS events   " << nMbsEvents << endl;
  cout << "Total epochs       " << nTotEpochs << endl;
  cout << "Total hit messages " << nTotHits << endl;
  cout << "Total aux messages " << nTotAux << endl;
  cout << "Total events       " << nTotEvents << endl;
  cout << "===============================================================" << endl;






  return 0;

}



ClassImp(ScEventBuilder)