analyzer_example.cxx

Go to the documentation of this file.

// Simple example of how to define a class that derives from 
// TRootanaEventLoop, thereby abstracting a way a lot of the 
// default stuff associated with analyzer.cxx

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

#include "TRootanaEventLoop.hxx"
#include "TH1D.h"
#include "TV1190Data.hxx"
#include "TV792Data.hxx"

class MyTestLoop: public TRootanaEventLoop {

  int nnn ;

public:

  MyTestLoop() {
    nnn = 0;
  };

  virtual ~MyTestLoop() {};

  void BeginRun(int transition,int run,int time){
    std::cout << "Custom: begin run " << run << std::endl;
  }

  void EndRun(int transition,int run,int time){
    std::cout << "Custom end run " << run <<std::endl;
  }

  bool ProcessMidasEvent(TDataContainer& dataContainer){


    void *ptr;
    int size = dataContainer.GetMidasData().LocateBank(NULL, "FR10", &ptr);
    if (ptr){
      nnn++;
      if(nnn%100 == 0){
        std::cout << "Current run :  " << GetCurrentRunNumber() << std::endl;
        std::cout << "Trying to handle this event... " << size << " " << nnn << std::endl;
      }

      fOnlineHistDir->cd();
      char sname[256];
      sprintf(sname, "size%d", 0);
      TH1D* hsize = (TH1D*)fOnlineHistDir->Get(sname);
      if (!hsize){
        printf("Create [%s]\n", sname);
        hsize = new TH1D(sname, sname, 600, 0, 6000);
      }      
      hsize->Fill(size);
    }


    TV1190Data *v1190 = dataContainer.GetEventData<TV1190Data>("TDC0");
    if(v1190){ 

      std::cout << "TDC measurements for V1190" << std::endl;
      std::vector<TDCMeasurement>& measurements = v1190->GetMeasurements();
      for(unsigned int i = 0; i < measurements.size(); i++){
        TDCMeasurement tdcmeas = measurements[i];
        
        std::cout << "Measurement: " <<    tdcmeas.GetMeasurement() << " for tdc/chan "  <<
          tdcmeas.GetTDCNumber() << "/"<< tdcmeas.GetChannel() << std::endl;

      }   
      
    }

    TV792Data *v792 = dataContainer.GetEventData<TV792Data>("ADC0");
    if(v792 ){ 
 
      v792->Print();
      
    }
 
    return true;
  }


  // Describe some other command line argument
  void Usage(void){
    std::cout << "\t-D option: this is a fun new option " << std::endl;
  }

  // Define some other command line argument
  bool CheckOption(std::string option){
    const char* arg = option.c_str();
    if (strncmp(arg,"-D",2)==0){
      std::cout << arg+2 << std::endl;
      std::cout << "I'm happy with this flag!" << std::endl;
      return true;
    }

    return false;
  }

}; 






int main(int argc, char *argv[])
{

  MyTestLoop::CreateSingleton<MyTestLoop>();
  return MyTestLoop::Get().ExecuteLoop(argc, argv);

}