MHitTime.cpp File Reference

#include <stdio.h>
#include <stdlib.h>
#include <string>
#include <sstream>
#include <TH1I.h>
#include "TGlobalData.h"
#include "TSetupData.h"

Go to the source code of this file.

Functions
INT	MHitTime_init (void)
INT	MHitTime (EVENT_HEADER *, void *)

Variables
HNDLE	hDB
TGlobalData *	gData
TSetupData *	gSetup
static TH1I *	hTime_MuSC
static std::map< std::string, TH1I * >	hEventTimes
static std::map< std::string, TH1I * >	hBlockTimes
static const double	BLOCKLENGTH = 10000000.
static const double	CUTWINDOW = 10000.
ANA_MODULE	MHitTime_module

Function Documentation

INT MHitTime	(	EVENT_HEADER *	pheader,
		void *	pevent
	)

Definition at line 86 of file MHitTime.cpp.

References banks, CUTWINDOW, TSetupData::fBankToDetectorMap, TGlobalData::fPulseIslandToChannelMap, hTime_MuSC, and SUCCESS.

{

  // We assume the existance of all necessary TPI vectors
  std::map< std::string, std::vector<TPulseIsland*> >& tpis = gData->fPulseIslandToChannelMap;
  std::map< std::string, std::string>& dets = gSetup->fBankToDetectorMap;
  std::map< std::string, std::string>& banks = gSetup->fDetectorToBankMap;

  std::vector<TPulseIsland*> muons = tpis.at(banks.at("muSC"));

  std::map< std::string, std::vector<TPulseIsland*> >::iterator iBank;
  std::map< std::string, std::string>::iterator iDet;
  std::vector<TPulseIsland*>::iterator cMuon;
  std::map< std::string, std::vector<TPulseIsland*>::iterator > pulses;
  for (iBank = tpis.begin(); iBank != tpis.end(); iBank++)
    if(iBank->first != banks.at("muSC"))
      pulses[iBank->first] = iBank->second.begin();
      

  double t_mu;
  double t[2];
  std::string cDet, cBank;
  for (cMuon = muons.begin(); cMuon != muons.end(); cMuon++) {
    t_mu = (*cMuon)->GetPulseTime();
    hTime_MuSC->Fill(t_mu);
    for (iDet = dets.begin(); iDet != dets.end(); iDet++) {
      cBank = iDet->first;
      cDet = iDet->second;
      std::vector<TPulseIsland*>::iterator& cPulse = pulses.at(cBank);
      if (cPulse == tpis.at(cBank).end())
        continue;
      t[0] = (*cPulse)->GetPulseTime();
      t[1] = t[0] - t_mu;
      hBlockTimes.at(cDet)->Fill(t[0]);
      while (t[1] < 0.) {
        cPulse++;
        t[0] = (*cPulse)->GetPulseTime();
        t[1] = t[0] - t_mu;
        hBlockTimes.at(cDet)->Fill(t[0]);
      }
      while (t[1] < CUTWINDOW) {
        hEventTimes.at(cDet)->Fill(t[1]);
        cPulse++;
        t[0] = (*cPulse)->GetPulseTime();
        t[1] = t[0] - t_mu;
        hBlockTimes.at(cDet)->Fill(t[0]);
      }
    }
  }

  // Are there any pulses left?
  for (iDet = dets.begin(); iDet != dets.end(); iDet++) {
    cBank = iDet->first;
    cDet = iDet->second;
    std::vector<TPulseIsland*>::iterator cPulse = pulses.at(cBank);
    while(cPulse != tpis.at(cBank).end()) {
      t[0] = (*cPulse)->GetPulseTime();
      hBlockTimes[cDet]->Fill(t[0]);
      cPulse++;
    }
  }

  return SUCCESS;
}

INT MHitTime_init

(

void

)

Definition at line 60 of file MHitTime.cpp.

References BLOCKLENGTH, CUTWINDOW, TSetupData::fBankToDetectorMap, TGlobalData::fPulseIslandToChannelMap, hTime_MuSC, and SUCCESS.

{

  std::map< std::string, std::vector<TPulseIsland*> >& tpis = gData->fPulseIslandToChannelMap;
  std::map< std::string, std::vector<TPulseIsland*> >::iterator iBank;
  std::map< std::string, std::string> dets = gSetup->fBankToDetectorMap;
  
  std::string cBank, cDet;
  hBlockTimes.clear();
  hEventTimes.clear();
  hTime_MuSC = new TH1I("hTime_MuSC", "Hit Times of MuSC;Block Hit Time (ns);Counts",
                                100, 0., BLOCKLENGTH);
  for (iBank = tpis.begin(); iBank != tpis.end(); iBank++) {
    cBank = iBank->first;
    cDet = dets[cBank];
    if (cDet == "muSC")
      continue;
    hBlockTimes[cDet] = new TH1I(("hBlockTime_" + cDet).c_str(), ("Block Hit Times of " + cDet + ";Time (ns); Counts").c_str(),
                                 100, 0., BLOCKLENGTH);
    hEventTimes[cDet] = new TH1I(("hEventTime_" + cDet).c_str(), ("Event Hit Times of " + cDet +";Time (ns); Counts").c_str(),
                                 100, 0., CUTWINDOW);
  }
  return SUCCESS;
}

Variable Documentation

const double BLOCKLENGTH = 10000000.

static

Definition at line 42 of file MHitTime.cpp.

Referenced by MHitTime_init().

const double CUTWINDOW = 10000.

static

Definition at line 43 of file MHitTime.cpp.

Referenced by MHitTime(), and MHitTime_init().

TGlobalData* gData

Definition at line 54 of file analyzer.cpp.

TSetupData* gSetup

Definition at line 55 of file analyzer.cpp.

std::map< std::string, TH1I* > hBlockTimes

static

Definition at line 40 of file MHitTime.cpp.

HNDLE hDB

Definition at line 11 of file mucap_compress.cpp.

std::map< std::string, TH1I* > hEventTimes

static

Definition at line 39 of file MHitTime.cpp.

TH1I* hTime_MuSC

static

Definition at line 38 of file MHitTime.cpp.

Referenced by MHitTime(), and MHitTime_init().

ANA_MODULE MHitTime_module

Initial value:

=
{
  "MHitTime",                    
  "John R Quirk",              
  MHitTime,                      
  NULL,                          
  NULL,                          
  MHitTime_init,                 
  NULL,                          
  NULL,                          
  0,                             
  NULL,                          
}

Definition at line 45 of file MHitTime.cpp.