MCommonOnlineDisplayPlots.cpp File Reference

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

Go to the source code of this file.

Functions
INT	MCommonOnlineDisplayPlots_init (void)
INT	MCommonOnlineDisplayPlots_bor (INT)
INT	MCommonOnlineDisplayPlots (EVENT_HEADER *, void *)

Variables
HNDLE	hDB
TGlobalData *	gData
TSetupData *	gSetup
map< std::string, TH1I * >	height_histograms_map
map< std::string, TH1I * >	time_histograms_map
map< std::string, TH2D * >	shape_histograms_map
map< std::string, TH1I * >	latest_pulse_histograms_map
static TH1I *	hPulseRawCount
ANA_MODULE	MCommonOnlineDisplayPlots_module

Function Documentation

INT MCommonOnlineDisplayPlots	(	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 158 of file MCommonOnlineDisplayPlots.cpp.

References TGlobalData::fPulseIslandToChannelMap, TSetupData::GetDetectorName(), hPulseRawCount, and SUCCESS.

{
        // 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::string detname = gSetup->GetDetectorName(bankname);
          std::vector<TPulseIsland*> thePulses = theMapIter->second;
                        
          // Loop over the TPulseIslands and plot the histogram
          for (std::vector<TPulseIsland*>::iterator pulseIter = thePulses.begin(); pulseIter != thePulses.end(); pulseIter++) {

            // Make sure the histograms exist and then fill them
            // Also check that this pulse didn't underflow (i.e. has a sample value at any point of 4096)
            bool underflow = false;
            if (shape_histograms_map.find(bankname) != shape_histograms_map.end()) {
              TH2* shape_histogram = shape_histograms_map[bankname];
              TH1* latest_pulse_histogram = latest_pulse_histograms_map[bankname];
              
              latest_pulse_histogram->Reset();  
              
              std::vector<int> theSamples = (*pulseIter)->GetSamples();
              for (std::vector<int>::iterator sampleIter = theSamples.begin(); sampleIter != theSamples.end(); sampleIter++) {
                int sample_number = sampleIter - theSamples.begin();
                int sample_value = *sampleIter;

                if (sample_value == 4096) {
                  underflow = true;
                  break;
                }
                shape_histogram->Fill(sample_number, sample_value);
                latest_pulse_histogram->SetBinContent(sample_number, sample_value);
              }
            }

            if (underflow == false) {
              if (height_histograms_map.find(bankname) != height_histograms_map.end())
                height_histograms_map[bankname]->Fill((*pulseIter)->GetPulseHeight());

              if (time_histograms_map.find(bankname) != time_histograms_map.end())
                time_histograms_map[bankname]->Fill((*pulseIter)->GetPulseTime());
            }
          }

          hPulseRawCount->Fill(bankname.c_str(), thePulses.size());
        }
        return SUCCESS;
}

INT MCommonOnlineDisplayPlots_bor

(

INT

run_number)

Definition at line 138 of file MCommonOnlineDisplayPlots.cpp.

References TSetupData::fBankToDetectorMap, and hPulseRawCount.

                                                  {

  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;
    height_histograms_map[bankname]->Reset();
    time_histograms_map[bankname]->Reset();
    shape_histograms_map[bankname]->Reset();
    latest_pulse_histograms_map[bankname]->Reset();

  }

  hPulseRawCount->Reset();
}

INT MCommonOnlineDisplayPlots_init

(

)

This method initializes histograms.

Definition at line 67 of file MCommonOnlineDisplayPlots.cpp.

References TSetupData::fBankToDetectorMap, TSetupData::GetDetectorName(), hPulseRawCount, TSetupData::IsBostonCAEN(), TSetupData::IsFADC(), TSetupData::IsHoustonCAEN(), and SUCCESS.

{
  // The following histograms are created for each channel:
  // hPulseHeights: ADC value (x-axis) vs number of pulses (y-axis)
  // hPulseTimes: time stamp (x-axis) vs number of pulses (y-axis)
  // hPulseRawCount: number of pulses (y-axis) - channels on x-axis?
  // hPulseShapes: sample number (x-axis) vs ADC value (y-axis) vs pulse (z-axis)

  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);

    int n_digitizer_bits = 0;
    if (TSetupData::IsFADC(bankname))
      n_digitizer_bits = 12;
    else if (TSetupData::IsHoustonCAEN(bankname))
      n_digitizer_bits = 14; //?
    else if (TSetupData::IsBostonCAEN(bankname))
      n_digitizer_bits = 12;
    
    long max_adc_value = std::pow(2, n_digitizer_bits);

    // hPulseHeights
    std::string histname = "h" + bankname + "_Heights";
    std::string histtitle = "Plot of the pulse heights in the " + detname + " channels";
    TH1I* hPulseHeights = new TH1I(histname.c_str(), histtitle.c_str(), max_adc_value,0,max_adc_value);
    hPulseHeights->GetXaxis()->SetTitle("Pulse Height [ADC value]");
    hPulseHeights->GetYaxis()->SetTitle("Number of Pulses");
    height_histograms_map[bankname] = hPulseHeights;

    // hPulseTimes
    histname = "h" + bankname + "_Times";
    histtitle = "Plot of the pulse times in the " + detname + " channels";
    TH1I* hPulseTimes = new TH1I(histname.c_str(), histtitle.c_str(), 2000, 0, 2e8);
    hPulseTimes->GetXaxis()->SetTitle("Time");
    hPulseTimes->GetYaxis()->SetTitle("Number of Pulses");
    time_histograms_map[bankname] = hPulseTimes;

    // hPulseShapes
    histname = "h" + bankname + "_Shapes";
    histtitle = "Plot of the pulse shapes in the " + detname + " channels";
    TH2D* hPulseShapes = new TH2D(histname.c_str(), histtitle.c_str(), 400,-0.5,399.5,max_adc_value+1,0,max_adc_value+1);      
    hPulseShapes->GetXaxis()->SetTitle("Time Stamp");
    hPulseShapes->GetYaxis()->SetTitle("ADC Value");
    shape_histograms_map[bankname] = hPulseShapes;

    //hLatestPulse
    histname = "h" + bankname + "_LatestPulse";
    histtitle = "Plot of the latest pulse in the " + detname + " channels";
    TH1I* hLatestPulse = new TH1I(histname.c_str(), histtitle.c_str(), 64,-0.5,63.5);
    hLatestPulse->GetXaxis()->SetTitle("Time Stamp");
    hLatestPulse->GetYaxis()->SetTitle("ADC Value");
    latest_pulse_histograms_map[bankname] = hLatestPulse;

  }

  // hPulseRawCount
  std::string histname = "hPulseRawCount";
  std::string histtitle = "Plot of the raw counts in each channels";
  hPulseRawCount = new TH1I(histname.c_str(), histtitle.c_str(), 1,0,1);
  hPulseRawCount->GetXaxis()->SetTitle("Channel");
  hPulseRawCount->GetYaxis()->SetTitle("Number of Pulses");
  hPulseRawCount->SetBit(TH1::kCanRebin);

  return SUCCESS;
}

Variable Documentation

TGlobalData* gData

Definition at line 54 of file analyzer.cpp.

TSetupData* gSetup

Definition at line 55 of file analyzer.cpp.

HNDLE hDB

Definition at line 11 of file mucap_compress.cpp.

map<std::string, TH1I*> height_histograms_map

Definition at line 45 of file MCommonOnlineDisplayPlots.cpp.

TH1I* hPulseRawCount

static

Definition at line 49 of file MCommonOnlineDisplayPlots.cpp.

Referenced by MCommonOnlineDisplayPlots(), MCommonOnlineDisplayPlots_bor(), and MCommonOnlineDisplayPlots_init().

map<std::string, TH1I*> latest_pulse_histograms_map

Definition at line 48 of file MCommonOnlineDisplayPlots.cpp.

ANA_MODULE MCommonOnlineDisplayPlots_module

Initial value:

=
{
        "MCommonOnlineDisplayPlots",                    
        "Andrew Edmonds",              
        MCommonOnlineDisplayPlots,                      
        MCommonOnlineDisplayPlots_bor,                          
        NULL,                          
        MCommonOnlineDisplayPlots_init,                 
        NULL,                          
        NULL,                          
        0,                             
        NULL,                          
}

Definition at line 51 of file MCommonOnlineDisplayPlots.cpp.

map<std::string, TH2D*> shape_histograms_map

Definition at line 47 of file MCommonOnlineDisplayPlots.cpp.

map<std::string, TH1I*> time_histograms_map

Definition at line 46 of file MCommonOnlineDisplayPlots.cpp.