MPileUpPerBlock.cpp File Reference

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

Go to the source code of this file.

Functions
INT	MPileUpPerBlock_init (void)
INT	MPileUpPerBlock (EVENT_HEADER *, void *)

Variables
HNDLE	hDB
TGlobalData *	gData
TSetupData *	gSetup
static TH1I *	hPileUpPerBlock
static const double	pileup_tw = 10000.
ANA_MODULE	MPileUpPerBlock_module

Function Documentation

INT MPileUpPerBlock	(	EVENT_HEADER *	pheader,
		void *	pevent
	)

Definition at line 67 of file MPileUpPerBlock.cpp.

References TSetupData::fBankToDetectorMap, TGlobalData::fPulseIslandToChannelMap, hPileUpPerBlock, pileup_tw, and SUCCESS.

{

  unsigned int nPileUp = 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;
  std::vector<TPulseIsland*>::iterator curr_musc_pulse, next_musc_pulse;

  for (iBank = tpis.begin(); iBank != tpis.end(); iBank++) {
    if (dets.at(iBank->first) == "muSC") {
      musc_pulses = &tpis.at(iBank->first);
      curr_musc_pulse = musc_pulses->begin();
      next_musc_pulse = curr_musc_pulse + 1;
    }
  }

  double cTime = 0;
  for (; curr_musc_pulse != musc_pulses->end();) {
    if ((*next_musc_pulse)->GetPulseTime() - (*curr_musc_pulse)->GetPulseTime() < pileup_tw)
      nPileUp++;
    curr_musc_pulse++;
    next_musc_pulse++;
  }

  hPileUpPerBlock->Fill(nPileUp);

  return SUCCESS;
}

INT MPileUpPerBlock_init

(

)

This method initializes histograms.

Definition at line 59 of file MPileUpPerBlock.cpp.

References hPileUpPerBlock, and SUCCESS.

{
  
  hPileUpPerBlock = new TH1I("hPileUpPerBlock","Number of Pile Up; Percent; Counts",100,0.,1);
  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* hPileUpPerBlock

static

Definition at line 37 of file MPileUpPerBlock.cpp.

Referenced by MPileUpPerBlock(), and MPileUpPerBlock_init().

ANA_MODULE MPileUpPerBlock_module

Initial value:

=
{
  "MPileUpPerBlock",                    
  "John R Quirk",              
  MPileUpPerBlock,                      
  NULL,                          
  NULL,                          
  MPileUpPerBlock_init,                 
  NULL,                          
  NULL,                          
  0,                             
  NULL,                          
}

Definition at line 43 of file MPileUpPerBlock.cpp.

const double pileup_tw = 10000.

static

Definition at line 40 of file MPileUpPerBlock.cpp.

Referenced by MPileUpPerBlock().