HarpoAnalysePrint.cxx

Go to the documentation of this file.

// 
// File HarpoAnalysePrint.cxx
//
#include "HarpoAnalysePrint.h"
#include "HarpoConfig.h"
#include "HarpoDetSet.h"
//#include "HarpoDetSet.h"
#include "HarpoEvent.h"

#include "TFile.h"

#include <cstdlib>
#include <cstring>
#include <cassert>
#include <iostream>

ClassImp(HarpoAnalysePrint)

void   HarpoAnalysePrint::print()
 {
   ULong_t nd; // number of detectors
   HarpoEventHeader *hdr = fEvt->GetHeader();

   assert(hdr != NULL);
   hdr->print();
   
   for (nd = 0; nd < gkNDetectors; nd++) {
     //     if (fEvt->isdataExist(nd)) {
       HarpoDccMap *plane = fEvt->GetDccMap(nd);
       if (plane != NULL )plane->print();
 }
 }

 void   HarpoAnalysePrint::process()
 {
   nEvents++;
   std::cout << " Processing Event " << nEvents<< std::endl;
   //fEvt->print();
    ULong_t ndet;
    int i,j;
    Long_t lEnum;
    
    lEnum = fEvt->GetHeader()->GetEvtNo();
    hevtno->Fill(lEnum + 0.5);
    for(ndet=0;ndet<gkNDetectors;ndet++) {
      if (gHDetSet->isExist(ndet)) {
        HarpoDccMap * m = fEvt->GetDccMap(ndet);
        if ( m != NULL ) {
          //      for(j=0;j<NCHAN;j++) // Channels
          for(j=0;j<NALL;j++) // Channels
            for(i=0;i<NADC;i++) //Time bins
              {
                h10[ndet]->Fill(j+0.5,m->GetData(j,i));
                h11[ndet]->Fill(i+0.5,m->GetData(j,i));
              }                         
          
        } else {
          std::cout << "No Data Map for plane "<< ndet << std::endl;
        }
      }
    }

 }

 void   HarpoAnalysePrint::Init()
 {
   hevtno = new TH1F("hevtno","Event Number",300,0.0,3000.0);
   h10[0] = new TProfile("h10_x","Plane X Dcc vs Chan",304,0.0,304.0,0.,4096.);
   h10[1] = new TProfile("h10_y","Plane Y Dcc vs Chan",304,0.0,304.0,0.,4096.);
   h11[0] = new TProfile("h11_x","Plane X Dcc vs Time",512,0.0,512.0,0.,4096.);
   h11[1] = new TProfile("h11_y","Plane Y Dcc vs Time",512,0.0,512.0,0.,4096.);
 }

void   HarpoAnalysePrint::Save(char * /* mode */)
 {

   OpenHistFile("print");

   hevtno->Write();
   h10[0]->Write();
   h10[1]->Write();
   h11[0]->Write();
   h11[1]->Write();
 }