MDQ_IslandCounter
[Low Level Data Quality Checks]

Plots number of TPIs in each detector per MIDAS event with different normalizations. More...

Functions
INT	MDQ_IslandCounter_init (void)
INT	MDQ_IslandCounter (EVENT_HEADER *, void *)
INT	MDQ_IslandCounter_eor (INT)
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.
TH1F *	hDQ_TDCCheck_muSc
	Coount of number of muSc hits in a run in the TDC.
map< std::string, TH1F * >	DQ_IslandCounter_histograms_map
	Maps bank name to histogram of number of TPIs in MIDAS event.
map< std::string, TH1F * >	DQ_IslandCounter_histograms_normalised_map
	Same as DQ_IslandCounter_histograms_map, but normalized vertical axis to muon count according to TDC.
map< std::string, TH1F * >	DQ_IslandCounter_histograms_both_axes_normalised_map
	Same as DQ_IslandCounter_histograms_normalised_map, but additionally horizontal axis is normalized in the same way.
ANA_MODULE	MDQ_IslandCounter_module

---

Detailed Description

Plots number of TPIs in each detector per MIDAS event with different normalizations.

Author:

Andrew Edmonds

Creates several plots.

---

Function Documentation

INT MDQ_IslandCounter	(	EVENT_HEADER *	pheader,
		void *	pevent
	)

This method fills the histograms

Definition at line 178 of file MDQ_IslandCounter.cpp.

References DQ_IslandCounter_histograms_map, DQ_IslandCounter_histograms_normalised_map, TGlobalData::fPulseIslandToChannelMap, and TSetupData::GetDetectorName().

{
        // 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 mapIter = pulse_islands_map.begin(); mapIter != pulse_islands_map.end(); ++mapIter) 
        {
          std::string bankname = mapIter->first;
          std::string detname = gSetup->GetDetectorName(bankname);
          std::vector<TPulseIsland*> thePulses = mapIter->second;
                        
          // Make sure the histograms exist and then fill them
          if (DQ_IslandCounter_histograms_map.find(bankname) != DQ_IslandCounter_histograms_map.end()) {
            int n_pulse_islands = thePulses.size();
            
            DQ_IslandCounter_histograms_map[bankname]->Fill(n_pulse_islands);
            DQ_IslandCounter_histograms_normalised_map[bankname]->Fill(n_pulse_islands);
          }
        }
        return SUCCESS;
}

INT MDQ_IslandCounter_eor

(

INT

run_number

)

This method does any last minute things to the histograms at the end of the run

Definition at line 138 of file MDQ_IslandCounter.cpp.

References DQ_IslandCounter_histograms_both_axes_normalised_map, DQ_IslandCounter_histograms_normalised_map, TGlobalData::fPulseIslandToChannelMap, TSetupData::GetDetectorName(), and hDQ_TDCCheck_muSc.

                                          {

  // 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 mapIter = pulse_islands_map.begin(); mapIter != pulse_islands_map.end(); ++mapIter) {

    std::string bankname = mapIter->first;
    std::string detname = gSetup->GetDetectorName(bankname);
      
    // Make sure the histograms exist and then fill them
    if (DQ_IslandCounter_histograms_normalised_map.find(bankname) != DQ_IslandCounter_histograms_normalised_map.end()) {

      TH1F* normalised_histogram = DQ_IslandCounter_histograms_normalised_map[bankname];
      normalised_histogram->Scale(1./hDQ_TDCCheck_muSc->GetEntries());

      // to get the normalised x-axis we need to go through and fill a second histogram
      TH1F* both_axes_normalised_histogram = DQ_IslandCounter_histograms_both_axes_normalised_map[bankname];
      for (int iBin = 1; iBin < normalised_histogram->GetNbinsX(); ++iBin) {
        double normalised_histogram_bin_content = normalised_histogram->GetBinContent(iBin);
        double normalised_histogram_bin_center = normalised_histogram->GetBinCenter(iBin);
        both_axes_normalised_histogram->Fill(normalised_histogram_bin_center / hDQ_TDCCheck_muSc->GetEntries(), normalised_histogram_bin_content);
      }
    }
  }

  return SUCCESS;
}

INT MDQ_IslandCounter_init

(

)

This method initializes histograms.

Definition at line 86 of file MDQ_IslandCounter.cpp.

References DQ_IslandCounter_histograms_both_axes_normalised_map, DQ_IslandCounter_histograms_map, DQ_IslandCounter_histograms_normalised_map, TSetupData::fBankToDetectorMap, and TSetupData::GetDetectorName().

{
    // See if the DataQuality_LowLevel/ directory already exists
  if (!gDirectory->Cd("DataQuality_LowLevel")) {
    
    std::string dir_name("DataQuality_LowLevel/");
    gDirectory->mkdir(dir_name.c_str());
    gDirectory->Cd(dir_name.c_str());
  }

  // Create a histogram 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);

    // hDQ_IslandCounter_[DetName]_[BankName]
    std::string histname = "hDQ_IslandCounter_" + detname + "_" + bankname;
    std::string histtitle = "Distribution of the number of islands per event in " + detname;
    TH1F* hDQ_Histogram = new TH1F(histname.c_str(), histtitle.c_str(), 3500, 0, 3500);
    hDQ_Histogram->GetXaxis()->SetTitle("Number of TPulseIslands");
    hDQ_Histogram->GetYaxis()->SetTitle("Number of Events");
    DQ_IslandCounter_histograms_map[bankname] = hDQ_Histogram;

    // The normalised histogram
    histname += "_normalised";
    histtitle += " (normalised)";
    TH1F* hDQ_Histogram_Normalised = new TH1F(histname.c_str(), histtitle.c_str(), 3500,0,3500);
    hDQ_Histogram_Normalised->GetXaxis()->SetTitle("Number of TPulseIslands");
    std::string yaxislabel = hDQ_Histogram->GetYaxis()->GetTitle();
    yaxislabel += " per TDC muSc Hit";
    hDQ_Histogram_Normalised->GetYaxis()->SetTitle(yaxislabel.c_str());
    DQ_IslandCounter_histograms_normalised_map[bankname] = hDQ_Histogram_Normalised;

    histname += "_both_axes_normalised";
    histtitle += " (both axes normalised)";
    TH1F* hDQ_Histogram_Both_Axes_Normalised = new TH1F(histname.c_str(), histtitle.c_str(), 3500,0,0.0035);
    hDQ_Histogram_Both_Axes_Normalised->GetXaxis()->SetTitle("Number of TPulseIslands per muSc TDC Hit");
    yaxislabel = hDQ_Histogram->GetYaxis()->GetTitle();
    yaxislabel += " per TDC muSc Hit";
    hDQ_Histogram_Both_Axes_Normalised->GetYaxis()->SetTitle(yaxislabel.c_str());
    DQ_IslandCounter_histograms_both_axes_normalised_map[bankname] = hDQ_Histogram_Both_Axes_Normalised;
  }

  gDirectory->Cd("/MidasHists/");
  return SUCCESS;
}

---

Variable Documentation

DQ_IslandCounter_histograms_both_axes_normalised_map

Same as DQ_IslandCounter_histograms_normalised_map, but additionally horizontal axis is normalized in the same way.

Definition at line 68 of file MDQ_IslandCounter.cpp.

Referenced by MDQ_IslandCounter_eor(), and MDQ_IslandCounter_init().

DQ_IslandCounter_histograms_map

Maps bank name to histogram of number of TPIs in MIDAS event.

Definition at line 66 of file MDQ_IslandCounter.cpp.

Referenced by MDQ_IslandCounter(), and MDQ_IslandCounter_init().

DQ_IslandCounter_histograms_normalised_map

Same as DQ_IslandCounter_histograms_map, but normalized vertical axis to muon count according to TDC.

Definition at line 67 of file MDQ_IslandCounter.cpp.

Referenced by MDQ_IslandCounter(), MDQ_IslandCounter_eor(), and MDQ_IslandCounter_init().

TGlobalData* gData

Object to hold data used and produced by modules throughout alcapana stage of analysis.

Definition at line 76 of file analyzer.cpp.

TSetupData* gSetup

Hardware information about digitizers and detectors to be used during alcapana stage of analysis.

Definition at line 80 of file analyzer.cpp.

HNDLE hDB

TH1F* hDQ_TDCCheck_muSc

Coount of number of muSc hits in a run in the TDC.

Definition at line 84 of file MDQ_TDCCheck.cpp.

ANA_MODULE MDQ_IslandCounter_module

Initial value:

{
        "MDQ_IslandCounter",                    
        "Andrew Edmonds",              
        MDQ_IslandCounter,                      
        NULL,                          
        MDQ_IslandCounter_eor,                          
        MDQ_IslandCounter_init,                 
        NULL,                          
        NULL,                          
        0,                             
        NULL,                          
}

Definition at line 70 of file MDQ_IslandCounter.cpp.