src/common/MPulseLengths.cpp File Reference

#include <stdio.h>
#include <stdlib.h>
#include <string>
#include <map>
#include <utility>
#include <sstream>
#include "midas.h"
#include <TH2.h>
#include "TOctalFADCIsland.h"
#include "TOctalFADCBankReader.h"
#include "TGlobalData.h"
#include "TSetupData.h"
#include "TPulseIsland.h"

Go to the source code of this file.

Functions
INT	MPulseLengths_init (void)
INT	MPulseLengths (EVENT_HEADER *, void *)
double	GetClockTickForChannel (string bank_name)
Variables
HNDLE	hDB
TGlobalData *	gData
	Object to hold data used and produced by modules throughout alcapana stage of analysis.
TSetupData *	gSetup
	Hardware information about digitizers and detectors to be used during alcapana stage of analysis.
static vector < TOctalFADCBankReader * >	fadc_bank_readers
static TH2I *	average_length_histogram
map< std::string, TH1I * >	length_histograms_map
ANA_MODULE	MPulseLengths_module

---

Function Documentation

double GetClockTickForChannel

(

string

bank_name

)

INT MPulseLengths	(	EVENT_HEADER *	pheader,
		void *	pevent
	)

This method processes one MIDAS block, producing a vector of TOctalFADCIsland objects from the raw Octal FADC data.

Definition at line 99 of file MPulseLengths.cpp.

References average_length_histogram, TGlobalData::fPulseIslandToChannelMap, and length_histograms_map.

{
        // 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;
          
          if (thePulses.size() != 0) {
            // Loop over the TPulseIslands and plot the histogram
            double total_length_of_pulses = 0;
            for (std::vector<TPulseIsland*>::iterator thePulseIter = thePulses.begin(); thePulseIter != thePulses.end(); thePulseIter++) {
              length_histograms_map[bankname]->Fill((*thePulseIter)->GetPulseLength());
              total_length_of_pulses += (*thePulseIter)->GetPulseLength();
            }
            // Fill the histogram
            average_length_histogram->Fill(bankname.c_str(), midas_event_number, total_length_of_pulses / thePulses.size());
          }
        }
        return SUCCESS;
}

INT MPulseLengths_init

(

)

This method initializes histograms.

Definition at line 65 of file MPulseLengths.cpp.

References average_length_histogram, TSetupData::fBankToDetectorMap, TSetupData::GetDetectorName(), and length_histograms_map.

{
  // This histogram has the pulse lengths 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 bankname = mapIter->first;
    std::string detname = gSetup->GetDetectorName(bankname);
    std::string histname = "h" + detname + "_Lengths";
    std::string histtitle = "Plot of the pulse lengths for the " + detname + " detector";
    TH1I* hDetLengths = new TH1I(histname.c_str(), histtitle.c_str(), 100,0,100);
    hDetLengths->GetXaxis()->SetTitle("Pulse Lengths [N Samples]");
    hDetLengths->GetYaxis()->SetTitle("Arbitrary Unit");
    hDetLengths->SetBit(TH1::kCanRebin);

    length_histograms_map[bankname] = hDetLengths;
  }

  std::string histname = "hAvgPulseLengthsPerChannel";
  std::string histtitle = "Plot of the average pulse lengths per event for the each channel";
  average_length_histogram = new TH2I(histname.c_str(), histtitle.c_str(), 1,0,1, 5000,0,5000);
  average_length_histogram->GetXaxis()->SetTitle("Bank Name");
  average_length_histogram->GetYaxis()->SetTitle("MIDAS Event Number");
  average_length_histogram->SetBit(TH1::kCanRebin);

  return SUCCESS;
}

---

Variable Documentation

TH2I* average_length_histogram [static]

Definition at line 46 of file MPulseLengths.cpp.

Referenced by MPulseLengths(), and MPulseLengths_init().

vector<TOctalFADCBankReader*> fadc_bank_readers [static]

Definition at line 45 of file MPulseLengths.cpp.

HNDLE hDB

map<std::string, TH1I*> length_histograms_map

Definition at line 47 of file MPulseLengths.cpp.

Referenced by MPulseLengths(), and MPulseLengths_init().

ANA_MODULE MPulseLengths_module

Initial value:

{
        "MPulseLengths",                    
        "Andrew Edmonds",              
        MPulseLengths,                      
        NULL,                          
        NULL,                          
        MPulseLengths_init,                 
        NULL,                          
        NULL,                          
        0,                             
        NULL,                          
}

Definition at line 49 of file MPulseLengths.cpp.