HarpoPmm2AnalyseEvProp.cxx

Go to the documentation of this file.

// 
// File HarpoPmm2AnalyseEvProp.cxx
//
#include "HarpoPmm2AnalyseEvProp.h"
#include "HarpoConfig.h"
#include "Pmm2Event.h"

#include "TFile.h"

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

const int MaxTime = 16777216; // 2^24 Time Stamp overflow value

ClassImp(HarpoPmm2AnalyseEvProp)

void   HarpoPmm2AnalyseEvProp::print()
 {
   //TODO
 }

 void   HarpoPmm2AnalyseEvProp::process()
 {
   if( verbose ) std::cout << " Processing Event " << nEvents<< std::endl;
   int itdiff; // Difference between times stamps in events N and N -1
   //fEvt->print();

   Pmm2Event *anaEvt; //pointer to PMM2 part
   if (gHDetSet->isExist(PMM2)) 
    {
      anaEvt = static_cast<Pmm2Event *>(fEvt->GetDetEvent(PMM2));
      anaEvt->print();
    }
   else
     {
       std::cout << "No PMM2 data for this event" << std::endl;
       return;
     }

   int imes;
   int esize = anaEvt->GetHeader()->eventSize;
   int ievnum = anaEvt->GetHeader()->eventNumb;

   helen->Fill(esize+0.5);  // Number of hits
   Pmm2MesVect * pm = anaEvt->GetMesurements();
   //int esize = pm->Size();  // Also give number of Digi's

   if( verbose ) pm->print();

   for(imes=0;imes<esize;imes++) {
     // i cann't use operator[] on pointer
     Pmm2MesVect &m = *pm;
   
     // *** This to methos to get channels must be the same *********

     // access via local reference m 
     //int ch = m[imes].getChNum(); // channel

     // we can use at in plase of brackets
     int ch = pm->at(imes).getChNum(); // channel

     heprof->Fill(ch + 0.5); // Channel profile
     
     int iadc = m[imes].getCharge(); // Charge in adc counts
     int ig = m[imes].getGain(); // Adc Gain bit
     hadc[ig]->Fill(iadc + 0.5); // All adc's values
     
     int ipps = m[imes].getPPS(); // PPS counter
     int its  = m[imes].getTimeStamp(); // Time stamp not corrected
     //std::cout << "pps " << ipps << " ts " << its << std::endl;
     hpps->Fill(ipps + 0.5, its + 0.5);

     htvse->Fill(ievnum + 0.5, its + 0.5);
     
     if ( (imes == 0) && (itsold != -1) )
       {
         itdiff = its - itsold;
         if ( itdiff < 0 ) {
           itdiff = MaxTime - itdiff;
         }
         //       std::cout << "e " << nEvents << " diff " << itdiff << std::endl;
         htdiff->Fill(itdiff + 0.5);
         htdiff2->Fill(itdiff + 0.5);
       }
     
     itsold = its;
     
     int iramp = m[imes].getTRamp(); // TDC Ramp bit
     int iftime  = m[imes].getFTime(); // Fine Time
     hramp->Fill(iramp + 0.5, iftime + 0.5);
     
   } // Digi Loop
   nEvents++;
 }

 void   HarpoPmm2AnalyseEvProp::Init()
 {
   nEvents = 0;
   itsold = -1;
   verbose = gHConfig->GetVerbose();

   helen = new TH1F("hELen","Hist per Event",20,0.0,19.0);
   heprof = new TH1F("hEProf","Hits per Chan",20,0.0,19.0);
   hadc[0] = new TH1F("hAdcg0","All Adc counts Gain 0",256,0.0,1024.0);
   hadc[1] = new TH1F("hAdcg1","All Adc counts Gain 1",256,0.0,1024.0);
   hpps = new TH2F("hPps","Time Stamp vs PPS",20,0.0,19.0,256,0.0,0.0);
   htvse = new TH2F("hTvsE","Time Stamp vs Ev Num",512,0.0,0.0,512,0.0,0.0);
   hramp = new TH2F("hRamp","Fine Time vs TDC Ramp",4,0.0,3.0,256,0.0,0.0); 
   htdiff = new TH1F("hTDiff","Time Diff",1024,0.0,1000000.0);
   htdiff2 = new TH1F("hTDiff2","Time Diff low",1024,0.0,15000.0);
}

void   HarpoPmm2AnalyseEvProp::Save(char * /* mode */)
 {
   TString * hstFile = gHConfig->GetHistFile();
   if ( hstFile == NULL ) {
     std::cout << gHConfig->GetProgramName() << " " <<
       "No Hist File name gived, use default" << std::endl;
     hstFile = new TString("pmm2evprop.root");
   }

   TFile *hf = new TFile(hstFile->Data(),"RECREATE");
   helen->Write();
   heprof->Write();
   hadc[0]->Write();
   hadc[1]->Write();
   hpps->Write();
   htvse->Write();
   hramp->Write();
   htdiff->Write();
   htdiff2->Write();
   hf->Close();
 }