MStoppedMuonsPerBlock.cpp File Reference

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

Go to the source code of this file.

Functions
INT	MStoppedMuonsPerBlock_init (void)
INT	MStoppedMuonsPerBlock (EVENT_HEADER *, void *)

Variables
HNDLE	hDB
TGlobalData *	gData
TSetupData *	gSetup
static TH1I *	hNumberOfStoppedMuonsPerBlock
static const double	tw_muveto = 100.
static const double	tw_sc = 100.
static const double	tw_si = 100.
ANA_MODULE	MStoppedMuonsPerBlock_module

Function Documentation

INT MStoppedMuonsPerBlock	(	EVENT_HEADER *	pheader,
		void *	pevent
	)

Definition at line 69 of file MStoppedMuonsPerBlock.cpp.

References TSetupData::fBankToDetectorMap, TGlobalData::fPulseIslandToChannelMap, hNumberOfStoppedMuonsPerBlock, SUCCESS, tw_muveto, tw_sc, and tw_si.

{

  unsigned int nStopped = 0;

  // 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::vector<TPulseIsland*> >::iterator iBank;
  std::vector<TPulseIsland*> *musc_pulses = NULL, *muveto_pulses = NULL, *scr_pulses = NULL, *scl_pulses = NULL;
  std::vector< std::vector<TPulseIsland*>* > sir_pulses, sil_pulses;
  std::vector<TPulseIsland*>::iterator c_musc_pulse, c_muveto_pulse, c_scr_pulse, c_scl_pulse;
  std::vector< std::vector<TPulseIsland*>::iterator > c_sir_pulse, c_sil_pulse;

  std::string cDet, cBank;
  int nSiL = 0, nSiR = 0;
  for (iBank = tpis.begin(); iBank != tpis.end(); iBank++) {
    cBank = iBank->first;
    cDet = dets[cBank];
    if (cDet == "muSc") {
      musc_pulses = &tpis.at(cBank);
      c_musc_pulse = musc_pulses->begin();
    } else if (cDet == "muVeto"){
      muveto_pulses = &tpis.at(cBank);
      c_muveto_pulse = muveto_pulses->begin();
    } else if (cDet == "ScR") {
      scr_pulses = &tpis.at(cBank);
      c_scr_pulse = scr_pulses->begin();
    } else if (cDet == "ScL") {
      scl_pulses = &tpis.at(cBank);
      c_scl_pulse = scl_pulses->begin();
    } else if (cDet.substr(0,2) == "Si" && cDet.substr(cDet.size()-5,4) == "fast") {
      if(cDet[2] == 'R') {
        sir_pulses.push_back(&tpis.at(cBank));
        c_sir_pulse.push_back(std::vector<TPulseIsland*>::iterator());
        nSiR++;
      } else {
        sil_pulses.push_back(&tpis.at(cBank));
        c_sil_pulse.push_back(std::vector<TPulseIsland*>::iterator());
        nSiL++;
      }
    }
  }


  double cTime = 0;
  for (; c_musc_pulse != musc_pulses->end(); c_musc_pulse++) {
    cTime = (*c_musc_pulse)->GetPulseTime();

    while ((*c_muveto_pulse)->GetPulseTime() < cTime)
      c_muveto_pulse++;
    while ((*c_scl_pulse)->GetPulseTime() < cTime)
      c_scl_pulse++;
    while ((*c_scr_pulse)->GetPulseTime() < cTime)
      c_scr_pulse++;
    for(int iSiR = 0; iSiR < nSiR; iSiR++)
      while ((*c_sir_pulse[iSiR])->GetPulseTime() < cTime)
        c_sir_pulse[iSiR]++;
    for (int iSiL = 0; iSiL < nSiL; iSiL++)
      while ((*c_sil_pulse[iSiL])->GetPulseTime() < cTime)
        c_sil_pulse[iSiL]++;

    if ((*c_muveto_pulse)->GetPulseTime() - cTime < tw_muveto)
      continue;
    if ((*c_scr_pulse)->GetPulseTime() - cTime < tw_sc)
      continue;
    if ((*c_scl_pulse)->GetPulseTime() - cTime < tw_sc)
      continue;
    for (int iSiR = 0; iSiR < nSiR; iSiR++)
      if ((*c_sir_pulse[iSiR])->GetPulseTime() - cTime < tw_si)
        continue;
    for (int iSiL = 0; iSiL < nSiL; iSiL++)
      if ((*c_sil_pulse[iSiL])->GetPulseTime() - cTime < tw_si)
        continue;

    nStopped++;
  }

  hNumberOfStoppedMuonsPerBlock->Fill(nStopped);

  return SUCCESS;
}

INT MStoppedMuonsPerBlock_init

(

)

This method initializes histograms.

Definition at line 61 of file MStoppedMuonsPerBlock.cpp.

References hNumberOfStoppedMuonsPerBlock, and SUCCESS.

{
  
  hNumberOfStoppedMuonsPerBlock = new TH1I("hNumberOfStoppedMuons","Number of (Possibly) Stopped Muons",2000,0.,2000.);
  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.

TH1I* hNumberOfStoppedMuonsPerBlock

static

Definition at line 37 of file MStoppedMuonsPerBlock.cpp.

Referenced by MStoppedMuonsPerBlock(), and MStoppedMuonsPerBlock_init().

ANA_MODULE MStoppedMuonsPerBlock_module

Initial value:

=
{
  "MStoppedMuonsPerBlock",                    
  "John R Quirk",              
  MStoppedMuonsPerBlock,                      
  NULL,                          
  NULL,                          
  MStoppedMuonsPerBlock_init,                 
  NULL,                          
  NULL,                          
  0,                             
  NULL,                          
}

Definition at line 45 of file MStoppedMuonsPerBlock.cpp.

const double tw_muveto = 100.

static

Definition at line 40 of file MStoppedMuonsPerBlock.cpp.

Referenced by MStoppedMuonsPerBlock().

const double tw_sc = 100.

static

Definition at line 41 of file MStoppedMuonsPerBlock.cpp.

Referenced by MStoppedMuonsPerBlock().

const double tw_si = 100.

static

Definition at line 42 of file MStoppedMuonsPerBlock.cpp.

Referenced by MStoppedMuonsPerBlock().