MPulseTimeSeparation.cpp

Go to the documentation of this file.

/********************************************************************\

Name:         MPulseTimeSeparation
Created by:   Joe Grange

Contents:     A module to fill a histogram of the pulse times from each channel

\********************************************************************/

/* Standard includes */
#include <stdio.h>
#include <stdlib.h>
#include <string>
#include <map>
#include <utility>
#include <sstream>

/* MIDAS includes */
#include "midas.h"

/* ROOT includes */
#include <TH1.h>

/* AlCap includes */
#include "TOctalFADCIsland.h"
#include "TOctalFADCBankReader.h"
#include "TGlobalData.h"
#include "TSetupData.h"
#include "TPulseIsland.h"

using std::string;
using std::map;
using std::vector;
using std::pair;

/*-- Module declaration --------------------------------------------*/
INT  MPulseTimeSeparation_init(void);
INT  MPulseTimeSeparation(EVENT_HEADER*, void*);
vector<string> GetAllBankNames();
double GetClockTickForChannel(string bank_name);

extern HNDLE hDB;
extern TGlobalData* gData;
extern TSetupData* gSetup;

static vector<TOctalFADCBankReader*> fadc_bank_readers;
static std::map<std::string, TH1*>  time_separation_histogram_map;
static std::map<std::string, TH1*>  time_separation_wrt_first_pulse_histogram_map;

ANA_MODULE MPulseTimeSeparation_module =
{
        "MPulseTimeSeparation",                    /* module name           */
        "Joe Grange",              /* author                */
        MPulseTimeSeparation,                      /* event routine         */
        NULL,                          /* BOR routine           */
        NULL,                          /* EOR routine           */
        MPulseTimeSeparation_init,                 /* init routine          */
        NULL,                          /* exit routine          */
        NULL,                          /* parameter structure   */
        0,                             /* structure size        */
        NULL,                          /* initial parameters    */
};

INT MPulseTimeSeparation_init()
{
  // This histogram has the pulse times on the X-axis and the number of pulses on the Y-axis
  // One histogram is created for each detector

  std::map<std::string, std::string> bank_to_detector_map = gSetup->fBankToDetectorMap;
  for(std::map<std::string, std::string>::iterator mapIter = bank_to_detector_map.begin(); 
      mapIter != bank_to_detector_map.end(); mapIter++) { 

    std::string detname = gSetup->GetDetectorName(mapIter->first);
    std::string histname = "h" + detname + "_PulseSeparation";
    std::string histtitle = "Plot of the time difference between consecutive pulses for the " + detname + " detector";
    TH1I* hPulseTimeDiff = new TH1I(histname.c_str(),histtitle.c_str(),200,0,200);
    hPulseTimeDiff->GetXaxis()->SetTitle("Time Difference Between Consecutive Pulses [ms]");
    hPulseTimeDiff->GetYaxis()->SetTitle("Number of pulse pairs");
    hPulseTimeDiff->SetBit(TH1::kCanRebin);

    time_separation_histogram_map[mapIter->first] = hPulseTimeDiff;

    histname = "h" + detname + "_PulseSeparationWRTFirstPulse";
    histtitle = "Plot of the time difference wrt the first pulse for the " + detname + " detector";
    hPulseTimeDiff = new TH1I(histname.c_str(), histtitle.c_str(),200,0,200);
    hPulseTimeDiff->GetXaxis()->SetTitle("Time Difference Between wrt the First Pulse [ms]");
    hPulseTimeDiff->GetYaxis()->SetTitle("Number of pulse pairs");
    hPulseTimeDiff->SetBit(TH1::kCanRebin);

    time_separation_wrt_first_pulse_histogram_map[mapIter->first] = hPulseTimeDiff;
  }

  return SUCCESS;
}

INT MPulseTimeSeparation(EVENT_HEADER *pheader, void *pevent)
{
        // Get the event number
        int midas_event_number = pheader->serial_number;

        // Some typedefs
        typedef map<string, vector<TPulseIsland*> > TStringPulseIslandMap;
        typedef pair<string, vector<TPulseIsland*> > TStringPulseIslandPair;
        typedef map<string, vector<TPulseIsland*> >::iterator map_iterator;

        // Fetch a reference to the gData structure that stores a map
        // of (bank_name, vector<TPulseIsland*>) pairs
        TStringPulseIslandMap& pulse_islands_map =
                gData->fPulseIslandToChannelMap;

        // Loop over the map and get each bankname, vector pair
        for (map_iterator theMapIter = pulse_islands_map.begin(); theMapIter != pulse_islands_map.end(); theMapIter++) 
        {
          std::string bankname = theMapIter->first;
          std::vector<TPulseIsland*> thePulses = theMapIter->second;
                        
          // Loop over the TPulseIslands and plot the histogram
          for (std::vector<TPulseIsland*>::iterator thePulseIter = thePulses.begin(); thePulseIter != thePulses.end(); thePulseIter++) {
            
            //      if ((*thePulseIter)->GetPulseLength() == 4) {
              if (thePulseIter!=thePulses.begin()) 
                time_separation_histogram_map[bankname]->Fill( ((*thePulseIter)->GetTimeStamp() - (*(thePulseIter-1))->GetTimeStamp()) * (*thePulseIter)->GetClockTickInNs() * 1e-6); // ns->ms

              int initial_time_stamp = (*thePulses.begin())->GetTimeStamp();

              double pulse_frequency = 65; // Hz
              double pulse_period = (1 / pulse_frequency) * 1e9; // ns
         
              if (initial_time_stamp * (*thePulseIter)->GetClockTickInNs() > pulse_period) {
                initial_time_stamp += pulse_period / (*thePulseIter)->GetClockTickInNs();
              }

              time_separation_wrt_first_pulse_histogram_map[bankname]->Fill( ((*thePulseIter)->GetTimeStamp() - initial_time_stamp)  * (*thePulseIter)->GetClockTickInNs() * 1e-6); // ns->ms             
              //            }
          }
        }
        return SUCCESS;
}