MDQ_muScTDiff
[Low Level Data Quality Checks]

Correlations of TPI timestamps with timestamps in muSc. More...

Functions
INT	MDQ_muScTDiff_init (void)
INT	MDQ_muScTDiff (EVENT_HEADER *, void *)
INT	MDQ_muScTDiff_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.
map< std::string, TH1F * >	DQ_muScTDiff_histograms_map
	Map of bank name to time stamp correlation with muSc TPI time stamps. The muSc histogram exists but is empty.
map< std::string, TH1F * >	DQ_muScTDiff_histograms_normalised_map
	Same as DQ_muScTDiff_histograms_map, but normalized to number of muon hits according to TDC.
map< std::string, TH1F * >	DQ_muScTDiff_histograms_with_time_shift_map
	Same as DQ_muScTDiff_histograms_map, except each TPI's time stamp has had the expected timing correlation peak subtracted off. If correct, the peak should be at zero in these plots.
float	axis_limit = 50000
TH1F *	hDQ_TDCCheck_muSc
	Coount of number of muSc hits in a run in the TDC.
ANA_MODULE	MDQ_muScTDiff_module

---

Detailed Description

Correlations of TPI timestamps with timestamps in muSc.

Author:

Andrew Edmonds

Creates a number of histograms. The muSc histogram exists but remains unfilled.

---

Function Documentation

INT MDQ_muScTDiff	(	EVENT_HEADER *	pheader,
		void *	pevent
	)

This method fills the histograms

Definition at line 172 of file MDQ_muScTDiff.cpp.

References axis_limit, DQ_muScTDiff_histograms_map, DQ_muScTDiff_histograms_normalised_map, DQ_muScTDiff_histograms_with_time_shift_map, TGlobalData::fPulseIslandToChannelMap, TSetupData::GetBankName(), TSetupData::GetDetectorName(), TSetupData::GetTimeShift(), and TSetupData::Instance().

{
        // 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;

        // Get the muSc pulses
        std::string muSc_bankname = gSetup->GetBankName("muSc");
        std::vector<TPulseIsland*> theMuScPulses = pulse_islands_map[muSc_bankname];

        // 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);
            
            // Don't bother comparing muSc to muSc
            if (detname == "muSc")
              continue;

            std::vector<TPulseIsland*> thePulses = mapIter->second;
                        
            // Make sure the histograms exist (put here so that it find() only called once per detector)
            if (DQ_muScTDiff_histograms_map.find(bankname) != DQ_muScTDiff_histograms_map.end()) {
              
              // Get the histograms before looping through all the pulses
              TH1F* hDQ_muScTDiff = DQ_muScTDiff_histograms_map[bankname];
              TH1F* hDQ_muScTDiff_Norm = DQ_muScTDiff_histograms_normalised_map[bankname];
              TH1F* hDQ_muScTDiff_TimeShift = DQ_muScTDiff_histograms_with_time_shift_map[bankname];
              double time_shift = TSetupData::Instance()->GetTimeShift(bankname);

              // Loop over the muSc pulses
              for (std::vector<TPulseIsland*>::iterator muScPulseIter = theMuScPulses.begin(); muScPulseIter != theMuScPulses.end(); ++muScPulseIter) {

                // Loop over the TPulseIslands and plot the histogram
                for (std::vector<TPulseIsland*>::iterator pulseIter = thePulses.begin(); pulseIter != thePulses.end(); ++pulseIter) {

                  // Get the pulse times for the muSc pulse and the detector pulse
                  double muSc_time = (*muScPulseIter)->GetPulseTime();
                  double det_time = (*pulseIter)->GetPulseTime();
                  double tdiff = muSc_time - det_time;

                  // The pulses should be time-ordered so if the tdiff goes outside of the axis range, then we can just skip to the next muSc pulse
                  if (std::fabs(tdiff) > 2*axis_limit)
                    break;

                  // Fill the histogram
                  hDQ_muScTDiff->Fill(tdiff);
                  hDQ_muScTDiff_Norm->Fill(tdiff);
                  hDQ_muScTDiff_TimeShift->Fill(tdiff+time_shift);
                }
              }
            }
          }
        return SUCCESS;
}

INT MDQ_muScTDiff_eor

(

INT

run_number

)

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

Definition at line 142 of file MDQ_muScTDiff.cpp.

References DQ_muScTDiff_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);
      
    if (DQ_muScTDiff_histograms_normalised_map.find(bankname) != DQ_muScTDiff_histograms_normalised_map.end()) {

      // Normalise to the muSc hits
      DQ_muScTDiff_histograms_normalised_map[bankname]->Scale(1./hDQ_TDCCheck_muSc->GetEntries());
    }
  }

  return SUCCESS;
}

INT MDQ_muScTDiff_init

(

)

This method initializes histograms.

Definition at line 89 of file MDQ_muScTDiff.cpp.

References axis_limit, DQ_muScTDiff_histograms_map, DQ_muScTDiff_histograms_normalised_map, DQ_muScTDiff_histograms_with_time_shift_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_muScTDiff_[DetName]_[BankName]
    std::string histname = "hDQ_muScTDiff_" + detname + "_" + bankname;
    std::string histtitle = "Time differences between muSc and " + detname;
    TH1F* hDQ_Histogram = new TH1F(histname.c_str(), histtitle.c_str(), 20000, -axis_limit, axis_limit);
    std::string axislabel = "Time Difference (muSc - " + detname + ") [ns]";
    hDQ_Histogram->GetXaxis()->SetTitle(axislabel.c_str());
    hDQ_Histogram->GetYaxis()->SetTitle("Number of TPulseIslands");
    DQ_muScTDiff_histograms_map[bankname] = hDQ_Histogram;

    // The normalised histogram
    std::string normhistname = histname + "_normalised";
    std::string normhisttitle = histtitle + " (normalised)";
    TH1F* hDQ_Histogram_Normalised = new TH1F(normhistname.c_str(), normhisttitle.c_str(), 20000, -axis_limit, axis_limit);
    hDQ_Histogram_Normalised->GetXaxis()->SetTitle(axislabel.c_str());
    std::string yaxislabel = hDQ_Histogram->GetYaxis()->GetTitle();
    yaxislabel += " per TDC muSc Hit";
    hDQ_Histogram_Normalised->GetYaxis()->SetTitle(yaxislabel.c_str());
    DQ_muScTDiff_histograms_normalised_map[bankname] = hDQ_Histogram_Normalised;

    // The histogram with the time shifts
    std::string timeshift_histname = histname + "_with-time-shift";
    std::string timeshift_histtitle = histtitle + " (with time shift)";
    TH1F* hDQ_Histogram_TimeShift = new TH1F(timeshift_histname.c_str(), timeshift_histtitle.c_str(), 20000, -axis_limit, axis_limit);
    hDQ_Histogram_TimeShift->GetXaxis()->SetTitle(axislabel.c_str());
    yaxislabel = hDQ_Histogram->GetYaxis()->GetTitle();
    hDQ_Histogram_Normalised->GetYaxis()->SetTitle(yaxislabel.c_str());
    DQ_muScTDiff_histograms_with_time_shift_map[bankname] = hDQ_Histogram_TimeShift;
  }

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

---

Variable Documentation

float axis_limit = 50000

Definition at line 69 of file MDQ_muScTDiff.cpp.

Referenced by MDQ_muScTDiff(), and MDQ_muScTDiff_init().

DQ_muScTDiff_histograms_map

Map of bank name to time stamp correlation with muSc TPI time stamps. The muSc histogram exists but is empty.

Definition at line 66 of file MDQ_muScTDiff.cpp.

Referenced by MDQ_muScTDiff(), and MDQ_muScTDiff_init().

DQ_muScTDiff_histograms_normalised_map

Same as DQ_muScTDiff_histograms_map, but normalized to number of muon hits according to TDC.

Definition at line 67 of file MDQ_muScTDiff.cpp.

Referenced by MDQ_muScTDiff(), MDQ_muScTDiff_eor(), and MDQ_muScTDiff_init().

DQ_muScTDiff_histograms_with_time_shift_map

Same as DQ_muScTDiff_histograms_map, except each TPI's time stamp has had the expected timing correlation peak subtracted off. If correct, the peak should be at zero in these plots.

Definition at line 68 of file MDQ_muScTDiff.cpp.

Referenced by MDQ_muScTDiff(), and MDQ_muScTDiff_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_muScTDiff_module

Initial value:

{
        "MDQ_muScTDiff",                    
        "Andrew Edmonds",              
        MDQ_muScTDiff,                      
        NULL,                          
        MDQ_muScTDiff_eor,                          
        MDQ_muScTDiff_init,                 
        NULL,                          
        NULL,                          
        0,                             
        NULL,                          
}

Definition at line 73 of file MDQ_muScTDiff.cpp.