#include <GlobalHitsAnalyzer.h>

Inheritance diagram for GlobalHitsAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

	GlobalHitsAnalyzer (const edm::ParameterSet &)

	~GlobalHitsAnalyzer () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Member Functions
void	fillECal (const edm::Event &, const edm::EventSetup &)

void	fillG4MC (const edm::Event &)

void	fillHCal (const edm::Event &, const edm::EventSetup &)

void	fillMuon (const edm::Event &, const edm::EventSetup &)

void	fillTrk (const edm::Event &, const edm::EventSetup &)

Private Attributes
unsigned int	count

edm::InputTag	ECalEBSrc_

edm::EDGetTokenT< edm::PCaloHitContainer >	ECalEBSrc_Token_

edm::InputTag	ECalEESrc_

edm::EDGetTokenT< edm::PCaloHitContainer >	ECalEESrc_Token_

edm::InputTag	ECalESSrc_

edm::EDGetTokenT< edm::PCaloHitContainer >	ECalESSrc_Token_

std::string	fName

int	frequency

edm::InputTag	G4TrkSrc_

edm::EDGetTokenT< edm::SimTrackContainer >	G4TrkSrc_Token_

edm::InputTag	G4VtxSrc_

edm::EDGetTokenT< edm::SimVertexContainer >	G4VtxSrc_Token_

bool	getAllProvenances

edm::InputTag	HCalSrc_

edm::EDGetTokenT< edm::PCaloHitContainer >	HCalSrc_Token_

std::string	label

MonitorElement *	meCaloEcal [2]

MonitorElement *	meCaloEcalE [2]

MonitorElement *	meCaloEcalEta

MonitorElement *	meCaloEcalPhi

MonitorElement *	meCaloEcalToF [2]

MonitorElement *	meCaloHcal [2]

MonitorElement *	meCaloHcalE [2]

MonitorElement *	meCaloHcalEta

MonitorElement *	meCaloHcalPhi

MonitorElement *	meCaloHcalToF [2]

MonitorElement *	meCaloPreSh [2]

MonitorElement *	meCaloPreShE [2]

MonitorElement *	meCaloPreShEta

MonitorElement *	meCaloPreShPhi

MonitorElement *	meCaloPreShToF [2]

MonitorElement *	meGeantTrkE

MonitorElement *	meGeantTrkPt

MonitorElement *	meGeantVtxEta

MonitorElement *	meGeantVtxMulti

MonitorElement *	meGeantVtxPhi

MonitorElement *	meGeantVtxRad [2]

MonitorElement *	meGeantVtxX [2]

MonitorElement *	meGeantVtxY [2]

MonitorElement *	meGeantVtxZ [2]

MonitorElement *	meMCG4Trk [2]

MonitorElement *	meMCG4Vtx [2]

MonitorElement *	meMCRGP [2]

MonitorElement *	meMuon [2]

MonitorElement *	meMuonCscToF [2]

MonitorElement *	meMuonCscZ

MonitorElement *	meMuonDtR

MonitorElement *	meMuonDtToF [2]

MonitorElement *	meMuonEta

MonitorElement *	meMuonPhi

MonitorElement *	meMuonRpcBR

MonitorElement *	meMuonRpcBToF [2]

MonitorElement *	meMuonRpcFToF [2]

MonitorElement *	meMuonRpcFZ

MonitorElement *	meTrackerPx [2]

MonitorElement *	meTrackerPxBR

MonitorElement *	meTrackerPxBToF

MonitorElement *	meTrackerPxEta

MonitorElement *	meTrackerPxFToF

MonitorElement *	meTrackerPxFZ

MonitorElement *	meTrackerPxPhi

MonitorElement *	meTrackerSi [2]

MonitorElement *	meTrackerSiBR

MonitorElement *	meTrackerSiBToF

MonitorElement *	meTrackerSiEta

MonitorElement *	meTrackerSiFToF

MonitorElement *	meTrackerSiFZ

MonitorElement *	meTrackerSiPhi

edm::InputTag	MuonCscSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	MuonCscSrc_Token_

edm::InputTag	MuonDtSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	MuonDtSrc_Token_

edm::InputTag	MuonRpcSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	MuonRpcSrc_Token_

int	nMuonHits

int	nPxlHits

int	nRawGenPart

int	nSiHits

bool	printProvenanceInfo

edm::InputTag	PxlBrlHighSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	PxlBrlHighSrc_Token_

edm::InputTag	PxlBrlLowSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	PxlBrlLowSrc_Token_

edm::InputTag	PxlFwdHighSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	PxlFwdHighSrc_Token_

edm::InputTag	PxlFwdLowSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	PxlFwdLowSrc_Token_

edm::InputTag	SiTECHighSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	SiTECHighSrc_Token_

edm::InputTag	SiTECLowSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	SiTECLowSrc_Token_

edm::InputTag	SiTIBHighSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	SiTIBHighSrc_Token_

edm::InputTag	SiTIBLowSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	SiTIBLowSrc_Token_

edm::InputTag	SiTIDHighSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	SiTIDHighSrc_Token_

edm::InputTag	SiTIDLowSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	SiTIDLowSrc_Token_

edm::InputTag	SiTOBHighSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	SiTOBHighSrc_Token_

edm::InputTag	SiTOBLowSrc_

edm::EDGetTokenT< edm::PSimHitContainer >	SiTOBLowSrc_Token_

bool	testNumber

bool	validEB

bool	validEE

bool	validG4trkContainer

bool	validG4VtxContainer

bool	validHcal

bool	validHepMCevt

bool	validMuonCSC

bool	validMuonDt

bool	validMuonRPC

bool	validPresh

bool	validPxlBrlHigh

bool	validPxlBrlLow

bool	validPxlFwdHigh

bool	validPxlFwdLow

bool	validSiTECHigh

bool	validSiTECLow

bool	validSiTIBHigh

bool	validSiTIBLow

bool	validSiTIDHigh

bool	validSiTIDLow

bool	validSiTOBHigh

bool	validSiTOBLow

int	verbosity

int	vtxunit

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Class to fill dqm monitor elements from existing EDM file

Author: M. Strang SUNY-Buffalo

Definition at line 74 of file GlobalHitsAnalyzer.h.

Constructor & Destructor Documentation

GlobalHitsAnalyzer::GlobalHitsAnalyzer ( const edm::ParameterSet & iPSet )

explicit

Definition at line 18 of file GlobalHitsAnalyzer.cc.

References ECalEBSrc_, ECalEBSrc_Token_, ECalEESrc_, ECalEESrc_Token_, ECalESSrc_, ECalESSrc_Token_, fName, frequency, getAllProvenances, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HCalSrc_, HCalSrc_Token_, mps_fire::i, edm::InputTag::instance(), edm::InputTag::label(), meCaloEcal, meCaloEcalE, meCaloEcalEta, meCaloEcalPhi, meCaloEcalToF, meCaloHcal, meCaloHcalE, meCaloHcalEta, meCaloHcalPhi, meCaloHcalToF, meCaloPreSh, meCaloPreShE, meCaloPreShEta, meCaloPreShPhi, meCaloPreShToF, meGeantTrkE, meGeantTrkPt, meGeantVtxEta, meGeantVtxMulti, meGeantVtxPhi, meGeantVtxRad, meGeantVtxX, meGeantVtxY, meGeantVtxZ, meMCG4Trk, meMCG4Vtx, meMCRGP, meMuon, meMuonCscToF, meMuonCscZ, meMuonDtR, meMuonDtToF, meMuonEta, meMuonPhi, meMuonRpcBR, meMuonRpcBToF, meMuonRpcFToF, meMuonRpcFZ, meTrackerPx, meTrackerPxBR, meTrackerPxBToF, meTrackerPxEta, meTrackerPxFToF, meTrackerPxFZ, meTrackerPxPhi, meTrackerSi, meTrackerSiBR, meTrackerSiBToF, meTrackerSiEta, meTrackerSiFToF, meTrackerSiFZ, meTrackerSiPhi, MuonCscSrc_, MuonCscSrc_Token_, MuonDtSrc_, MuonDtSrc_Token_, MuonRpcSrc_, MuonRpcSrc_Token_, printProvenanceInfo, PxlBrlHighSrc_, PxlBrlHighSrc_Token_, PxlBrlLowSrc_, PxlBrlLowSrc_Token_, PxlFwdHighSrc_, PxlFwdHighSrc_Token_, PxlFwdLowSrc_, PxlFwdLowSrc_Token_, SiTECHighSrc_, SiTECHighSrc_Token_, SiTECLowSrc_, SiTECLowSrc_Token_, SiTIBHighSrc_, SiTIBHighSrc_Token_, SiTIBLowSrc_, SiTIBLowSrc_Token_, SiTIDHighSrc_, SiTIDHighSrc_Token_, SiTIDLowSrc_, SiTIDLowSrc_Token_, SiTOBHighSrc_, SiTOBHighSrc_Token_, SiTOBLowSrc_, SiTOBLowSrc_Token_, AlCaHLTBitMon_QueryRunRegistry::string, testNumber, validEB, validEE, validG4trkContainer, validG4VtxContainer, validHcal, validHepMCevt, validMuonCSC, validMuonDt, validMuonRPC, validPresh, validPxlBrlHigh, validPxlBrlLow, validPxlFwdHigh, validPxlFwdLow, validSiTECHigh, validSiTECLow, validSiTIBHigh, validSiTIBLow, validSiTIDHigh, validSiTIDLow, validSiTOBHigh, validSiTOBLow, verbosity, and vtxunit.

                                                                    :
   fName(""), verbosity(0), frequency(0), vtxunit(0), label(""), 
   getAllProvenances(false), printProvenanceInfo(false), testNumber(false),
   G4VtxSrc_Token_( consumes<edm::SimVertexContainer>((iPSet.getParameter<edm::InputTag>("G4VtxSrc"))) ),
   G4TrkSrc_Token_( consumes<edm::SimTrackContainer>(iPSet.getParameter<edm::InputTag>("G4TrkSrc")) ),
   count(0)
 {
   consumesMany<edm::HepMCProduct>();
   std::string MsgLoggerCat = "GlobalHitsAnalyzer_GlobalHitsAnalyzer";
 
   // get information from parameter set
   fName = iPSet.getUntrackedParameter<std::string>("Name");
   verbosity = iPSet.getUntrackedParameter<int>("Verbosity");
   frequency = iPSet.getUntrackedParameter<int>("Frequency");
   vtxunit = iPSet.getUntrackedParameter<int>("VtxUnit");
   edm::ParameterSet m_Prov =
     iPSet.getParameter<edm::ParameterSet>("ProvenanceLookup");
   getAllProvenances = 
     m_Prov.getUntrackedParameter<bool>("GetAllProvenances");
   printProvenanceInfo = 
     m_Prov.getUntrackedParameter<bool>("PrintProvenanceInfo");
   testNumber = iPSet.getUntrackedParameter<bool>("testNumber");
 
   //get Labels to use to extract information
   PxlBrlLowSrc_ = iPSet.getParameter<edm::InputTag>("PxlBrlLowSrc");
   PxlBrlHighSrc_ = iPSet.getParameter<edm::InputTag>("PxlBrlHighSrc");
   PxlFwdLowSrc_ = iPSet.getParameter<edm::InputTag>("PxlFwdLowSrc");
   PxlFwdHighSrc_ = iPSet.getParameter<edm::InputTag>("PxlFwdHighSrc");
 
   SiTIBLowSrc_ = iPSet.getParameter<edm::InputTag>("SiTIBLowSrc");
   SiTIBHighSrc_ = iPSet.getParameter<edm::InputTag>("SiTIBHighSrc");
   SiTOBLowSrc_ = iPSet.getParameter<edm::InputTag>("SiTOBLowSrc");
   SiTOBHighSrc_ = iPSet.getParameter<edm::InputTag>("SiTOBHighSrc");
   SiTIDLowSrc_ = iPSet.getParameter<edm::InputTag>("SiTIDLowSrc");
   SiTIDHighSrc_ = iPSet.getParameter<edm::InputTag>("SiTIDHighSrc");
   SiTECLowSrc_ = iPSet.getParameter<edm::InputTag>("SiTECLowSrc");
   SiTECHighSrc_ = iPSet.getParameter<edm::InputTag>("SiTECHighSrc");
 
   MuonCscSrc_ = iPSet.getParameter<edm::InputTag>("MuonCscSrc");
   MuonDtSrc_ = iPSet.getParameter<edm::InputTag>("MuonDtSrc");
   MuonRpcSrc_ = iPSet.getParameter<edm::InputTag>("MuonRpcSrc");
 
   ECalEBSrc_ = iPSet.getParameter<edm::InputTag>("ECalEBSrc");
   ECalEESrc_ = iPSet.getParameter<edm::InputTag>("ECalEESrc");
   ECalESSrc_ = iPSet.getParameter<edm::InputTag>("ECalESSrc");
 
   HCalSrc_ = iPSet.getParameter<edm::InputTag>("HCalSrc");
 
   // fix for consumes
   PxlBrlLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("PxlBrlLowSrc"));
   PxlBrlHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("PxlBrlHighSrc"));
   PxlFwdLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("PxlFwdLowSrc"));
   PxlFwdHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("PxlFwdHighSrc"));
 
   SiTIBLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTIBLowSrc"));
   SiTIBHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTIBHighSrc"));
   SiTOBLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTOBLowSrc"));
   SiTOBHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTOBHighSrc"));
   SiTIDLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTIDLowSrc"));
   SiTIDHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTIDHighSrc"));
   SiTECLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTECLowSrc"));
   SiTECHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTECHighSrc"));
 
   MuonCscSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("MuonCscSrc"));
   MuonDtSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("MuonDtSrc"));
   MuonRpcSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("MuonRpcSrc"));
 
   ECalEBSrc_Token_ = consumes<edm::PCaloHitContainer>(iPSet.getParameter<edm::InputTag>("ECalEBSrc"));
   ECalEESrc_Token_ = consumes<edm::PCaloHitContainer>(iPSet.getParameter<edm::InputTag>("ECalEESrc"));
   ECalESSrc_Token_ = consumes<edm::PCaloHitContainer>(iPSet.getParameter<edm::InputTag>("ECalESSrc"));
   HCalSrc_Token_ = consumes<edm::PCaloHitContainer>(iPSet.getParameter<edm::InputTag>("HCalSrc"));
 
   // determine whether to process subdetector or not
   validHepMCevt = iPSet.getUntrackedParameter<bool>("validHepMCevt");
   validG4VtxContainer = 
     iPSet.getUntrackedParameter<bool>("validG4VtxContainer");
   validG4trkContainer = 
     iPSet.getUntrackedParameter<bool>("validG4trkContainer");
   validPxlBrlLow = iPSet.getUntrackedParameter<bool>("validPxlBrlLow",true);
   validPxlBrlHigh = iPSet.getUntrackedParameter<bool>("validPxlBrlHigh",true);
   validPxlFwdLow = iPSet.getUntrackedParameter<bool>("validPxlFwdLow",true);
   validPxlFwdHigh = iPSet.getUntrackedParameter<bool>("validPxlFwdHigh",true);
   validSiTIBLow = iPSet.getUntrackedParameter<bool>("validSiTIBLow",true);
   validSiTIBHigh = iPSet.getUntrackedParameter<bool>("validSiTIBHigh",true);
   validSiTOBLow = iPSet.getUntrackedParameter<bool>("validSiTOBLow",true);
   validSiTOBHigh = iPSet.getUntrackedParameter<bool>("validSiTOBHigh",true);
   validSiTIDLow = iPSet.getUntrackedParameter<bool>("validSiTIDLow",true);
   validSiTIDHigh = iPSet.getUntrackedParameter<bool>("validSiTIDHigh",true);
   validSiTECLow = iPSet.getUntrackedParameter<bool>("validSiTECLow",true);
   validSiTECHigh = iPSet.getUntrackedParameter<bool>("validSiTECHigh",true);
   validMuonCSC = iPSet.getUntrackedParameter<bool>("validMuonCSC",true);
   validMuonDt = iPSet.getUntrackedParameter<bool>("validMuonDt",true);
   validMuonRPC = iPSet.getUntrackedParameter<bool>("validMuonRPC",true);
   validEB = iPSet.getUntrackedParameter<bool>("validEB",true);
   validEE = iPSet.getUntrackedParameter<bool>("validEE",true);
   validPresh = iPSet.getUntrackedParameter<bool>("validPresh",true);
   validHcal = iPSet.getUntrackedParameter<bool>("validHcal",true);  
 
   // use value of first digit to determine default output level (inclusive)
   // 0 is none, 1 is basic, 2 is fill output, 3 is gather output
   verbosity %= 10;
 
   // print out Parameter Set information being used
   if (verbosity >= 0) {
     edm::LogInfo(MsgLoggerCat) 
       << "\n===============================\n"
       << "Initialized as EDAnalyzer with parameter values:\n"
       << "    Name                  = " << fName << "\n"
       << "    Verbosity             = " << verbosity << "\n"
       << "    Frequency             = " << frequency << "\n"
       << "    VtxUnit               = " << vtxunit << "\n"
       << "    GetProv               = " << getAllProvenances << "\n"
       << "    PrintProv             = " << printProvenanceInfo << "\n"
       << "    PxlBrlLowSrc          = " << PxlBrlLowSrc_.label() 
       << ":" << PxlBrlLowSrc_.instance() << "\n"
       << "    PxlBrlHighSrc         = " << PxlBrlHighSrc_.label() 
       << ":" << PxlBrlHighSrc_.instance() << "\n"
       << "    PxlFwdLowSrc          = " << PxlFwdLowSrc_.label() 
       << ":" << PxlBrlLowSrc_.instance() << "\n"
       << "    PxlFwdHighSrc         = " << PxlFwdHighSrc_.label() 
       << ":" << PxlBrlHighSrc_.instance() << "\n"
       << "    SiTIBLowSrc           = " << SiTIBLowSrc_.label() 
       << ":" << SiTIBLowSrc_.instance() << "\n"
       << "    SiTIBHighSrc          = " << SiTIBHighSrc_.label() 
       << ":" << SiTIBHighSrc_.instance() << "\n"
       << "    SiTOBLowSrc           = " << SiTOBLowSrc_.label() 
       << ":" << SiTOBLowSrc_.instance() << "\n"
       << "    SiTOBHighSrc          = " << SiTOBHighSrc_.label() 
       << ":" << SiTOBHighSrc_.instance() << "\n"
       << "    SiTIDLowSrc           = " << SiTIDLowSrc_.label() 
       << ":" << SiTIDLowSrc_.instance() << "\n"
       << "    SiTIDHighSrc          = " << SiTIDHighSrc_.label() 
       << ":" << SiTIDHighSrc_.instance() << "\n"
       << "    SiTECLowSrc           = " << SiTECLowSrc_.label() 
       << ":" << SiTECLowSrc_.instance() << "\n"
       << "    SiTECHighSrc          = " << SiTECHighSrc_.label() 
       << ":" << SiTECHighSrc_.instance() << "\n"
       << "    MuonCscSrc            = " << MuonCscSrc_.label() 
       << ":" << MuonCscSrc_.instance() << "\n"
       << "    MuonDtSrc             = " << MuonDtSrc_.label() 
       << ":" << MuonDtSrc_.instance() << "\n"
       << "    MuonRpcSrc            = " << MuonRpcSrc_.label() 
       << ":" << MuonRpcSrc_.instance() << "\n"
       << "    ECalEBSrc             = " << ECalEBSrc_.label() 
       << ":" << ECalEBSrc_.instance() << "\n"
       << "    ECalEESrc             = " << ECalEESrc_.label() 
       << ":" << ECalEESrc_.instance() << "\n"
       << "    ECalESSrc             = " << ECalESSrc_.label() 
       << ":" << ECalESSrc_.instance() << "\n"
       << "    HCalSrc               = " << HCalSrc_.label() 
       << ":" << HCalSrc_.instance() << "\n"
       << "\n"
       << "    validHepMCevt         = "
       << ":" <<  validHepMCevt << "\n"
       << "    validG4VtxContainer   = "
       << ":" <<  validG4VtxContainer << "\n"
       << "    validG4trkContainer   = "
       << ":" <<  validG4trkContainer << "\n"
       << "    validPxlBrlLow        = "
       << ":" <<  validPxlBrlLow << "\n"
       << "    validPxlBrlHigh       = "
       << ":" <<  validPxlBrlHigh << "\n"
       << "    validPxlFwdLow        = "
       << ":" <<  validPxlFwdLow << "\n"
       << "    validPxlFwdHigh       = "
       << ":" <<  validPxlFwdHigh << "\n"
       << "    validSiTIBLow         = "
       << ":" <<  validSiTIBLow << "\n"
       << "    validSiTIBHigh        = "
       << ":" <<  validSiTIBHigh << "\n"
       << "    validSiTOBLow         = "
       << ":" <<  validSiTOBLow << "\n"
       << "    validSiTOBHigh        = "
       << ":" <<  validSiTOBHigh << "\n"
       << "    validSiTIDLow         = "
       << ":" <<  validSiTIDLow << "\n"
       << "    validSiTIDHigh        = "
       << ":" <<  validSiTIDHigh << "\n"
       << "    validSiTECLow         = "
       << ":" <<  validSiTECLow << "\n"
       << "    validSiTECHigh        = "
       << ":" <<  validSiTECHigh << "\n"
       << "    validMuonCSC          = "
       << ":" <<  validMuonCSC << "\n"
       << "    validMuonDt           = "
       << ":" <<  validMuonDt << "\n"
       << "    validMuonRPC          = "
       << ":" <<  validMuonRPC << "\n"
       << "    validEB               = "
       << ":" <<  validEB << "\n"
       << "    validEE               = "
       << ":" <<  validEE << "\n"
       << "    validPresh            = "
       << ":" <<  validPresh << "\n"
       << "    validHcal             = "
       << ":" <<  validHcal << "\n"
       << "===============================\n";
   }
 
   // initialize monitor elements
   for (Int_t i = 0; i < 2; ++i) {
     meMCRGP[i]     = nullptr;
     meMCG4Vtx[i]   = nullptr;
     meGeantVtxX[i] = nullptr;
     meGeantVtxY[i] = nullptr;
     meGeantVtxZ[i] = nullptr; 
     meMCG4Trk[i]   = nullptr;
     meCaloEcal[i]  = nullptr;
     meCaloEcalE[i] = nullptr;
     meCaloEcalToF[i] = nullptr;
     meCaloPreSh[i]   = nullptr;
     meCaloPreShE[i]  = nullptr;
     meCaloPreShToF[i]= nullptr;
     meCaloHcal[i]    = nullptr;
     meCaloHcalE[i]   = nullptr;
     meCaloHcalToF[i] = nullptr;
     meTrackerPx[i]   = nullptr;
     meTrackerSi[i]   = nullptr;
     meMuon[i]        = nullptr;
     meMuonDtToF[i]   = nullptr;
     meMuonCscToF[i]  = nullptr;
     meMuonRpcFToF[i] = nullptr;
     meMuonRpcBToF[i] = nullptr;
     meGeantVtxRad[i] = nullptr;
   }
   meGeantTrkPt   = nullptr;
   meGeantTrkE    = nullptr;
   meGeantVtxEta  = nullptr;
   meGeantVtxPhi  = nullptr;
   meGeantVtxMulti= nullptr;
   meCaloEcalPhi  = nullptr;
   meCaloEcalEta  = nullptr;
   meCaloPreShPhi = nullptr;
   meCaloPreShEta = nullptr;
   meCaloHcalPhi  = nullptr;
   meCaloHcalEta  = nullptr;
   meTrackerPxPhi = nullptr;
   meTrackerPxEta = nullptr;
   meTrackerPxBToF= nullptr;
   meTrackerPxBR  = nullptr;
   meTrackerPxFToF= nullptr;
   meTrackerPxFZ  = nullptr;
   meTrackerSiPhi = nullptr;
   meTrackerSiEta = nullptr;
   meTrackerSiBToF= nullptr;
   meTrackerSiBR  = nullptr;
   meTrackerSiFToF= nullptr;
   meTrackerSiFZ  = nullptr;
   meMuonPhi      = nullptr;
   meMuonEta      = nullptr;
   meMuonDtR      = nullptr;
   meMuonCscZ     = nullptr;
   meMuonRpcBR    = nullptr;
   meMuonRpcFZ    = nullptr;
 
 }

GlobalHitsAnalyzer::~GlobalHitsAnalyzer ( )

override

Definition at line 275 of file GlobalHitsAnalyzer.cc.

275 {}

Member Function Documentation

void GlobalHitsAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 690 of file GlobalHitsAnalyzer.cc.

References count, edm::EventID::event(), fillECal(), fillG4MC(), fillHCal(), fillMuon(), fillTrk(), frequency, getAllProvenances, edm::Event::getAllStableProvenance(), mps_fire::i, edm::EventBase::id(), nevt, printProvenanceInfo, edm::EventID::run(), AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

 {
   std::string MsgLoggerCat = "GlobalHitsAnalyzer_analyze";
 
   // keep track of number of events processed
   ++count;
 
   // get event id information
   edm::RunNumber_t nrun = iEvent.id().run();
   edm::EventNumber_t nevt = iEvent.id().event();
 
   if (verbosity > 0) {
     edm::LogInfo(MsgLoggerCat)
       << "Processing run " << nrun << ", event " << nevt
       << " (" << count << " events total)";
   } else if (verbosity == 0) {
     if (nevt%frequency == 0 || nevt == 1) {
       edm::LogInfo(MsgLoggerCat)
     << "Processing run " << nrun << ", event " << nevt
     << " (" << count << " events total)";
     }
   }
 
   // look at information available in the event
   if (getAllProvenances) {
 
     std::vector<const edm::StableProvenance*> AllProv;
     iEvent.getAllStableProvenance(AllProv);
 
     if (verbosity >= 0)
       edm::LogInfo(MsgLoggerCat)
     << "Number of Provenances = " << AllProv.size();
 
     if (printProvenanceInfo && (verbosity >= 0)) {
       TString eventout("\nProvenance info:\n");      
 
       for (unsigned int i = 0; i < AllProv.size(); ++i) {
     eventout += "\n       ******************************";
     eventout += "\n       Module       : ";
     eventout += AllProv[i]->moduleLabel();
     eventout += "\n       ProductID    : ";
     eventout += AllProv[i]->productID().id();
     eventout += "\n       ClassName    : ";
     eventout += AllProv[i]->className();
     eventout += "\n       InstanceName : ";
     eventout += AllProv[i]->productInstanceName();
     eventout += "\n       BranchName   : ";
     eventout += AllProv[i]->branchName();
       }
       eventout += "\n       ******************************\n";
       edm::LogInfo(MsgLoggerCat) << eventout << "\n";
       printProvenanceInfo = false;
     }
     getAllProvenances = false;
   }
 
   // call fill functions
   //gather G4MC information from event
   fillG4MC(iEvent);
   // gather Tracker information from event
   fillTrk(iEvent,iSetup);
   // gather muon information from event
   fillMuon(iEvent, iSetup);
   // gather Ecal information from event
   fillECal(iEvent, iSetup);
   // gather Hcal information from event
   fillHCal(iEvent, iSetup);
 
   if (verbosity > 0)
     edm::LogInfo (MsgLoggerCat)
       << "Done gathering data from event.";
 
   return;
 }

void GlobalHitsAnalyzer::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	run,
		edm::EventSetup const &	es
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 277 of file GlobalHitsAnalyzer.cc.

References DQMStore::IBooker::book1D(), mps_fire::i, meCaloEcal, meCaloEcalE, meCaloEcalEta, meCaloEcalPhi, meCaloEcalToF, meCaloHcal, meCaloHcalE, meCaloHcalEta, meCaloHcalPhi, meCaloHcalToF, meCaloPreSh, meCaloPreShE, meCaloPreShEta, meCaloPreShPhi, meCaloPreShToF, meGeantTrkE, meGeantTrkPt, meGeantVtxEta, meGeantVtxMulti, meGeantVtxPhi, meGeantVtxRad, meGeantVtxX, meGeantVtxY, meGeantVtxZ, meMCG4Trk, meMCG4Vtx, meMCRGP, meMuon, meMuonCscToF, meMuonCscZ, meMuonDtR, meMuonDtToF, meMuonEta, meMuonPhi, meMuonRpcBR, meMuonRpcBToF, meMuonRpcFToF, meMuonRpcFZ, meTrackerPx, meTrackerPxBR, meTrackerPxBToF, meTrackerPxEta, meTrackerPxFToF, meTrackerPxFZ, meTrackerPxPhi, meTrackerSi, meTrackerSiBR, meTrackerSiBToF, meTrackerSiEta, meTrackerSiFToF, meTrackerSiFZ, meTrackerSiPhi, MonitorElement::setAxisTitle(), and DQMStore::IBooker::setCurrentFolder().

                                                                                                             {
   // book histograms
   Char_t hname[200];
   Char_t htitle[200];
 
   // MCGeant
   iBooker.setCurrentFolder("GlobalHitsV/MCGeant");
   sprintf(hname,"hMCRGP1");
   sprintf(htitle,"RawGenParticles");
   meMCRGP[0] = iBooker.book1D(hname,htitle,100,0.,5000.);
   sprintf(hname,"hMCRGP2");
   meMCRGP[1] = iBooker.book1D(hname,htitle,100,0.,500.);  
   for (Int_t i = 0; i < 2; ++i) {
     meMCRGP[i]->setAxisTitle("Number of Raw Generated Particles",1);
     meMCRGP[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hMCG4Vtx1");
   sprintf(htitle,"G4 Vertices");
   meMCG4Vtx[0] = iBooker.book1D(hname,htitle,150,0.,15000.);
   sprintf(hname,"hMCG4Vtx2");
   meMCG4Vtx[1] = iBooker.book1D(hname,htitle,100,-0.5,99.5); 
   for (Int_t i = 0; i < 2; ++i) {
     meMCG4Vtx[i]->setAxisTitle("Number of Vertices",1);
     meMCG4Vtx[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hMCG4Trk1");
   sprintf(htitle,"G4 Tracks");
   meMCG4Trk[0] = iBooker.book1D(hname,htitle,150,0.,15000.);
   sprintf(hname,"hMCG4Trk2");
   meMCG4Trk[1] = iBooker.book1D(hname,htitle,150,-0.5,99.5);    
   for (Int_t i = 0; i < 2; ++i) {
     meMCG4Trk[i]->setAxisTitle("Number of Tracks",1);
     meMCG4Trk[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hGeantVtxX1");
   sprintf(htitle,"Geant vertex x/micrometer");
   meGeantVtxX[0] = iBooker.book1D(hname,htitle,100,-8000000.,8000000.);
   sprintf(hname,"hGeantVtxX2");
   meGeantVtxX[1] = iBooker.book1D(hname,htitle,100,-50.,50.); 
   for (Int_t i = 0; i < 2; ++i) {
     meGeantVtxX[i]->setAxisTitle("x of Vertex (um)",1);
     meGeantVtxX[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hGeantVtxY1");
   sprintf(htitle,"Geant vertex y/micrometer");
   meGeantVtxY[0] = iBooker.book1D(hname,htitle,100,-8000000,8000000.);
   sprintf(hname,"hGeantVtxY2");
   meGeantVtxY[1] = iBooker.book1D(hname,htitle,100,-50.,50.); 
   for (Int_t i = 0; i < 2; ++i) {
     meGeantVtxY[i]->setAxisTitle("y of Vertex (um)",1);
     meGeantVtxY[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hGeantVtxZ1");
   sprintf(htitle,"Geant vertex z/millimeter");
   meGeantVtxZ[0] = iBooker.book1D(hname,htitle,100,-11000.,11000.);
   sprintf(hname,"hGeantVtxZ2");
   meGeantVtxZ[1] = iBooker.book1D(hname,htitle,200,-500.,500.);
   //meGeantVtxZ[1] = iBooker.book1D(hname,htitle,100,-250.,250.);
   for (Int_t i = 0; i < 2; ++i) {
     meGeantVtxZ[i]->setAxisTitle("z of Vertex (mm)",1);
     meGeantVtxZ[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hGeantTrkPt");
   sprintf(htitle,"Log10 Geant track pt/GeV");
   meGeantTrkPt = iBooker.book1D(hname,htitle,80,-4.,4.);
   meGeantTrkPt->setAxisTitle("Log10 pT of Track (GeV)",1);
   meGeantTrkPt->setAxisTitle("Count",2);
 
   sprintf(hname,"hGeantTrkE");
   sprintf(htitle,"Log10 Geant track E/GeV");
   meGeantTrkE = iBooker.book1D(hname,htitle,80,-4.,4.);
   meGeantTrkE->setAxisTitle("Log10 E of Track (GeV)",1);
   meGeantTrkE->setAxisTitle("Count",2);
 
   sprintf(hname,"hGeantVtxEta");
   sprintf(htitle,"Geant vertices eta");
   meGeantVtxEta = iBooker.book1D(hname,htitle,220,-5.5,5.5);
   meGeantVtxEta->setAxisTitle("eta of SimVertex",1);
   meGeantVtxEta->setAxisTitle("Count",2);
 
   sprintf(hname,"hGeantVtxPhi");
   sprintf(htitle,"Geant vertices phi/rad");
   meGeantVtxPhi = iBooker.book1D(hname,htitle,100,-3.2,3.2);
   meGeantVtxPhi->setAxisTitle("phi of SimVertex (rad)",1);
   meGeantVtxPhi->setAxisTitle("Count",2);
 
   sprintf(hname,"hGeantVtxRad1");
   sprintf(htitle,"Geant vertices radius/cm");
   meGeantVtxRad[0] = iBooker.book1D(hname,htitle,130,0.,130.);
   sprintf(hname,"hGeantVtxRad2");
   meGeantVtxRad[1] = iBooker.book1D(hname,htitle,100,0.,1000.);
   for (Int_t i = 0; i < 2; ++i) {
     meGeantVtxRad[i]->setAxisTitle("radius of SimVertex (cm)",1);
     meGeantVtxRad[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hGeantVtxMulti");
   sprintf(htitle,"Geant vertices outgoing multiplicity");
   meGeantVtxMulti = iBooker.book1D(hname,htitle,20,0.,20);
   meGeantVtxMulti->setAxisTitle("multiplicity of particles attached to a SimVertex",1);
   meGeantVtxMulti->setAxisTitle("Count",2);
 
   // ECal
   iBooker.setCurrentFolder("GlobalHitsV/ECals");
   sprintf(hname,"hCaloEcal1");
   sprintf(htitle,"Ecal hits");
   meCaloEcal[0] = iBooker.book1D(hname,htitle,100,0.,10000.);
   sprintf(hname,"hCaloEcal2");
   meCaloEcal[1] = iBooker.book1D(hname,htitle,100,-0.5,99.5);
 
   sprintf(hname,"hCaloEcalE1");
   sprintf(htitle,"Ecal hits, energy/GeV");
   meCaloEcalE[0] = iBooker.book1D(hname,htitle,100,0.,10.);
   sprintf(hname,"hCaloEcalE2");
   meCaloEcalE[1] = iBooker.book1D(hname,htitle,100,0.,0.1);
   sprintf(hname,"hCaloEcalToF1");
   sprintf(htitle,"Ecal hits, ToF/ns");
   meCaloEcalToF[0] = iBooker.book1D(hname,htitle,100,0.,1000.);
   sprintf(hname,"hCaloEcalToF2");
   meCaloEcalToF[1] = iBooker.book1D(hname,htitle,100,0.,100.);
  
   for (Int_t i = 0; i < 2; ++i) {
     meCaloEcal[i]->setAxisTitle("Number of Hits",1);
     meCaloEcal[i]->setAxisTitle("Count",2);
     meCaloEcalE[i]->setAxisTitle("Energy of Hits (GeV)",1);
     meCaloEcalE[i]->setAxisTitle("Count",2);
     meCaloEcalToF[i]->setAxisTitle("Time of Flight of Hits (ns)",1);
     meCaloEcalToF[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hCaloEcalPhi");
   sprintf(htitle,"Ecal hits, phi/rad");
   meCaloEcalPhi = iBooker.book1D(hname,htitle,100,-3.2,3.2);
   meCaloEcalPhi->setAxisTitle("Phi of Hits (rad)",1);
   meCaloEcalPhi->setAxisTitle("Count",2);
 
   sprintf(hname,"hCaloEcalEta");
   sprintf(htitle,"Ecal hits, eta");
   meCaloEcalEta = iBooker.book1D(hname,htitle,100,-5.5,5.5);
   meCaloEcalEta->setAxisTitle("Eta of Hits",1);
   meCaloEcalEta->setAxisTitle("Count",2);
 
   sprintf(hname,"hCaloPreSh1");
   sprintf(htitle,"PreSh hits");
   meCaloPreSh[0] = iBooker.book1D(hname,htitle,100,0.,10000.);
   sprintf(hname,"hCaloPreSh2");
   meCaloPreSh[1] = iBooker.book1D(hname,htitle,100,-0.5,99.5);
 
   sprintf(hname,"hCaloPreShE1");
   sprintf(htitle,"PreSh hits, energy/GeV");
   meCaloPreShE[0] = iBooker.book1D(hname,htitle,100,0.,10.);
   sprintf(hname,"hCaloPreShE2");
   meCaloPreShE[1] = iBooker.book1D(hname,htitle,100,0.,0.1);
 
   sprintf(hname,"hCaloPreShToF1");
   sprintf(htitle,"PreSh hits, ToF/ns");
   meCaloPreShToF[0] = iBooker.book1D(hname,htitle,100,0.,1000.);
   sprintf(hname,"hCaloPreShToF2");
   meCaloPreShToF[1] = iBooker.book1D(hname,htitle,100,0.,100.);
 
   for (Int_t i = 0; i < 2; ++i) {
     meCaloPreSh[i]->setAxisTitle("Number of Hits",1);
     meCaloPreSh[i]->setAxisTitle("Count",2);
     meCaloPreShE[i]->setAxisTitle("Energy of Hits (GeV)",1);
     meCaloPreShE[i]->setAxisTitle("Count",2);
     meCaloPreShToF[i]->setAxisTitle("Time of Flight of Hits (ns)",1);
     meCaloPreShToF[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hCaloPreShPhi");
   sprintf(htitle,"PreSh hits, phi/rad");
   meCaloPreShPhi = iBooker.book1D(hname,htitle,100,-3.2,3.2);
   meCaloPreShPhi->setAxisTitle("Phi of Hits (rad)",1);
   meCaloPreShPhi->setAxisTitle("Count",2);
 
   sprintf(hname,"hCaloPreShEta");
   sprintf(htitle,"PreSh hits, eta");
   meCaloPreShEta = iBooker.book1D(hname,htitle,100,-5.5,5.5);
   meCaloPreShEta->setAxisTitle("Eta of Hits",1);
   meCaloPreShEta->setAxisTitle("Count",2);
 
   // Hcal
   iBooker.setCurrentFolder("GlobalHitsV/HCals");
   sprintf(hname,"hCaloHcal1");
   sprintf(htitle,"Hcal hits");
   meCaloHcal[0] = iBooker.book1D(hname,htitle,100,0.,10000.);
   sprintf(hname,"hCaloHcal2");
   meCaloHcal[1] = iBooker.book1D(hname,htitle,100,-0.5,99.5);
 
   sprintf(hname,"hCaloHcalE1");
   sprintf(htitle,"Hcal hits, energy/GeV");
   meCaloHcalE[0] = iBooker.book1D(hname,htitle,100,0.,10.);
   sprintf(hname,"hCaloHcalE2");
   meCaloHcalE[1] = iBooker.book1D(hname,htitle,100,0.,0.1);
 
   sprintf(hname,"hCaloHcalToF1");
   sprintf(htitle,"Hcal hits, ToF/ns");
   meCaloHcalToF[0] = iBooker.book1D(hname,htitle,100,0.,1000.);
   sprintf(hname,"hCaloHcalToF2");
   meCaloHcalToF[1] = iBooker.book1D(hname,htitle,100,0.,100.);
 
   for (Int_t i = 0; i < 2; ++i) {
     meCaloHcal[i]->setAxisTitle("Number of Hits",1);
     meCaloHcal[i]->setAxisTitle("Count",2);
     meCaloHcalE[i]->setAxisTitle("Energy of Hits (GeV)",1);
     meCaloHcalE[i]->setAxisTitle("Count",2);
     meCaloHcalToF[i]->setAxisTitle("Time of Flight of Hits (ns)",1);
     meCaloHcalToF[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hCaloHcalPhi");
   sprintf(htitle,"Hcal hits, phi/rad");
   meCaloHcalPhi = iBooker.book1D(hname,htitle,100,-3.2,3.2);
   meCaloHcalPhi->setAxisTitle("Phi of Hits (rad)",1);
   meCaloHcalPhi->setAxisTitle("Count",2);
 
   sprintf(hname,"hCaloHcalEta");
   sprintf(htitle,"Hcal hits, eta");
   meCaloHcalEta = iBooker.book1D(hname,htitle,100,-5.5,5.5);
   meCaloHcalEta->setAxisTitle("Eta of Hits",1);
   meCaloHcalEta->setAxisTitle("Count",2);
   
   // SiPixels
   iBooker.setCurrentFolder("GlobalHitsV/SiPixels");
   sprintf(hname,"hTrackerPx1");
   sprintf(htitle,"Pixel hits");
   meTrackerPx[0] = iBooker.book1D(hname,htitle,100,0.,10000.);
   sprintf(hname,"hTrackerPx2");
   meTrackerPx[1] = iBooker.book1D(hname,htitle,100,-0.5,99.5);
   for (Int_t i = 0; i < 2; ++i) {
     meTrackerPx[i]->setAxisTitle("Number of Pixel Hits",1);
     meTrackerPx[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hTrackerPxPhi");
   sprintf(htitle,"Pixel hits phi/rad");
   meTrackerPxPhi = iBooker.book1D(hname,htitle,100,-3.2,3.2);
   meTrackerPxPhi->setAxisTitle("Phi of Hits (rad)",1);
   meTrackerPxPhi->setAxisTitle("Count",2);
 
   sprintf(hname,"hTrackerPxEta");
   sprintf(htitle,"Pixel hits eta");
   meTrackerPxEta = iBooker.book1D(hname,htitle,100,-3.5,3.5);
   meTrackerPxEta->setAxisTitle("Eta of Hits",1);
   meTrackerPxEta->setAxisTitle("Count",2);
 
   sprintf(hname,"hTrackerPxBToF");
   sprintf(htitle,"Pixel barrel hits, ToF/ns");
   meTrackerPxBToF = iBooker.book1D(hname,htitle,100,0.,40.);
   meTrackerPxBToF->setAxisTitle("Time of Flight of Hits (ns)",1);
   meTrackerPxBToF->setAxisTitle("Count",2);
 
   sprintf(hname,"hTrackerPxBR");
   sprintf(htitle,"Pixel barrel hits, R/cm");
   meTrackerPxBR = iBooker.book1D(hname,htitle,100,0.,50.);
   meTrackerPxBR->setAxisTitle("R of Hits (cm)",1);
   meTrackerPxBR->setAxisTitle("Count",2);
 
   sprintf(hname,"hTrackerPxFToF");
   sprintf(htitle,"Pixel forward hits, ToF/ns");
   meTrackerPxFToF = iBooker.book1D(hname,htitle,100,0.,50.);
   meTrackerPxFToF->setAxisTitle("Time of Flight of Hits (ns)",1);
   meTrackerPxFToF->setAxisTitle("Count",2);
 
   sprintf(hname,"hTrackerPxFZ");
   sprintf(htitle,"Pixel forward hits, Z/cm");
   meTrackerPxFZ = 
     iBooker.book1D(hname,htitle,200,-100.,100.);
   meTrackerPxFZ->setAxisTitle("Z of Hits (cm)",1);
   meTrackerPxFZ->setAxisTitle("Count",2);
 
   // SiStrips
   iBooker.setCurrentFolder("GlobalHitsV/SiStrips");
   sprintf(hname,"hTrackerSi1");
   sprintf(htitle,"Silicon hits");
   meTrackerSi[0] = iBooker.book1D(hname,htitle,100,0.,10000.);
   sprintf(hname,"hTrackerSi2");
   meTrackerSi[1] = iBooker.book1D(hname,htitle,100,-0.5,99.5);
   for (Int_t i = 0; i < 2; ++i) { 
     meTrackerSi[i]->setAxisTitle("Number of Silicon Hits",1);
     meTrackerSi[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hTrackerSiPhi");
   sprintf(htitle,"Silicon hits phi/rad");
   meTrackerSiPhi = iBooker.book1D(hname,htitle,100,-3.2,3.2);
   meTrackerSiPhi->setAxisTitle("Phi of Hits (rad)",1);
   meTrackerSiPhi->setAxisTitle("Count",2);
 
   sprintf(hname,"hTrackerSiEta");
   sprintf(htitle,"Silicon hits eta");
   meTrackerSiEta = iBooker.book1D(hname,htitle,100,-3.5,3.5);
   meTrackerSiEta->setAxisTitle("Eta of Hits",1);
   meTrackerSiEta->setAxisTitle("Count",2);
 
   sprintf(hname,"hTrackerSiBToF");
   sprintf(htitle,"Silicon barrel hits, ToF/ns");
   meTrackerSiBToF = iBooker.book1D(hname,htitle,100,0.,50.);
   meTrackerSiBToF->setAxisTitle("Time of Flight of Hits (ns)",1);
   meTrackerSiBToF->setAxisTitle("Count",2);
 
   sprintf(hname,"hTrackerSiBR");
   sprintf(htitle,"Silicon barrel hits, R/cm");
   meTrackerSiBR = iBooker.book1D(hname,htitle,100,0.,200.);
   meTrackerSiBR->setAxisTitle("R of Hits (cm)",1);
   meTrackerSiBR->setAxisTitle("Count",2);
 
   sprintf(hname,"hTrackerSiFToF");
   sprintf(htitle,"Silicon forward hits, ToF/ns");
   meTrackerSiFToF = iBooker.book1D(hname,htitle,100,0.,75.);
   meTrackerSiFToF->setAxisTitle("Time of Flight of Hits (ns)",1);
   meTrackerSiFToF->setAxisTitle("Count",2);
 
   sprintf(hname,"hTrackerSiFZ");
   sprintf(htitle,"Silicon forward hits, Z/cm");
   meTrackerSiFZ = iBooker.book1D(hname,htitle,200,-300.,300.);
   meTrackerSiFZ->setAxisTitle("Z of Hits (cm)",1);
   meTrackerSiFZ->setAxisTitle("Count",2);
 
   // Muon
   iBooker.setCurrentFolder("GlobalHitsV/Muons");
   sprintf(hname,"hMuon1");
   sprintf(htitle,"Muon hits");
   meMuon[0] = iBooker.book1D(hname,htitle,100,0.,10000.);
   sprintf(hname,"hMuon2");
   meMuon[1] = iBooker.book1D(hname,htitle,100,-0.5,99.5);
   for (Int_t i = 0; i < 2; ++i) { 
     meMuon[i]->setAxisTitle("Number of Muon Hits",1);
     meMuon[i]->setAxisTitle("Count",2);
   }  
 
   sprintf(hname,"hMuonPhi");
   sprintf(htitle,"Muon hits phi/rad");
   meMuonPhi = iBooker.book1D(hname,htitle,100,-3.2,3.2);
   meMuonPhi->setAxisTitle("Phi of Hits (rad)",1);
   meMuonPhi->setAxisTitle("Count",2);
 
   sprintf(hname,"hMuonEta");
   sprintf(htitle,"Muon hits eta");
   meMuonEta = iBooker.book1D(hname,htitle,100,-3.5,3.5);
   meMuonEta->setAxisTitle("Eta of Hits",1);
   meMuonEta->setAxisTitle("Count",2);
 
   sprintf(hname,"hMuonCscToF1");
   sprintf(htitle,"Muon CSC hits, ToF/ns");
   meMuonCscToF[0] = iBooker.book1D(hname,htitle,100,0.,250.);
   sprintf(hname,"hMuonCscToF2");
   meMuonCscToF[1] = iBooker.book1D(hname,htitle,100,0.,50.);
   for (Int_t i = 0; i < 2; ++i) {   
     meMuonCscToF[i]->setAxisTitle("Time of Flight of Hits (ns)",1);
     meMuonCscToF[i]->setAxisTitle("Count",2);
   }  
 
   sprintf(hname,"hMuonCscZ");
   sprintf(htitle,"Muon CSC hits, Z/cm");
   meMuonCscZ = iBooker.book1D(hname,htitle,200,-1500.,1500.);
   meMuonCscZ->setAxisTitle("Z of Hits (cm)",1);
   meMuonCscZ->setAxisTitle("Count",2);
 
   sprintf(hname,"hMuonDtToF1");
   sprintf(htitle,"Muon DT hits, ToF/ns");
   meMuonDtToF[0] = iBooker.book1D(hname,htitle,100,0.,250.);
   sprintf(hname,"hMuonDtToF2");
   meMuonDtToF[1] = iBooker.book1D(hname,htitle,100,0.,50.);
   for (Int_t i = 0; i < 2; ++i) {   
     meMuonDtToF[i]->setAxisTitle("Time of Flight of Hits (ns)",1);
     meMuonDtToF[i]->setAxisTitle("Count",2);
   } 
 
   sprintf(hname,"hMuonDtR");
   sprintf(htitle,"Muon DT hits, R/cm");
   meMuonDtR = iBooker.book1D(hname,htitle,100,0.,1500.); 
   meMuonDtR->setAxisTitle("R of Hits (cm)",1);
   meMuonDtR->setAxisTitle("Count",2);
 
   sprintf(hname,"hMuonRpcFToF1");
   sprintf(htitle,"Muon RPC forward hits, ToF/ns");
   meMuonRpcFToF[0] = iBooker.book1D(hname,htitle,100,0.,250.);
   sprintf(hname,"hMuonRpcFToF2");
   meMuonRpcFToF[1] = iBooker.book1D(hname,htitle,100,0.,50.);
   for (Int_t i = 0; i < 2; ++i) {   
     meMuonRpcFToF[i]->setAxisTitle("Time of Flight of Hits (ns)",1);
     meMuonRpcFToF[i]->setAxisTitle("Count",2);
   }  
   sprintf(hname,"hMuonRpcFZ");
   sprintf(htitle,"Muon RPC forward hits, Z/cm");
   meMuonRpcFZ = iBooker.book1D(hname,htitle,201,-1500.,1500.);
   meMuonRpcFZ->setAxisTitle("Z of Hits (cm)",1);
   meMuonRpcFZ->setAxisTitle("Count",2);
 
   sprintf(hname,"hMuonRpcBToF1");
   sprintf(htitle,"Muon RPC barrel hits, ToF/ns");
   meMuonRpcBToF[0] = iBooker.book1D(hname,htitle,100,0.,250.);
   sprintf(hname,"hMuonRpcBToF2");
   meMuonRpcBToF[1] = iBooker.book1D(hname,htitle,100,0.,50.);
   for (Int_t i = 0; i < 2; ++i) {   
     meMuonRpcBToF[i]->setAxisTitle("Time of Flight of Hits (ns)",1);
     meMuonRpcBToF[i]->setAxisTitle("Count",2);
   }
 
   sprintf(hname,"hMuonRpcBR");
   sprintf(htitle,"Muon RPC barrel hits, R/cm");
   meMuonRpcBR = iBooker.book1D(hname,htitle,100,0.,1500.);
   meMuonRpcBR->setAxisTitle("R of Hits (cm)",1);
   meMuonRpcBR->setAxisTitle("Count",2); 
 }

void GlobalHitsAnalyzer::fillECal	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 1647 of file GlobalHitsAnalyzer.cc.

References dEcal, gamEcalExtractorBlocks_cff::detector, ECalEBSrc_Token_, ECalEESrc_Token_, ECalESSrc_Token_, PV3DBase< T, PVType, FrameType >::eta(), MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), ecaldqm::getGeometry(), CaloGeometry::getSubdetectorGeometry(), mps_fire::i, edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), LogDebug, meCaloEcal, meCaloEcalE, meCaloEcalEta, meCaloEcalPhi, meCaloEcalToF, meCaloPreSh, meCaloPreShE, meCaloPreShEta, meCaloPreShPhi, meCaloPreShToF, PV3DBase< T, PVType, FrameType >::phi(), sdEcalBrl, sdEcalFwd, sdEcalPS, AlCaHLTBitMon_QueryRunRegistry::string, subdetector, validEB, validEE, validPresh, and verbosity.

Referenced by analyze().

 {
   std::string MsgLoggerCat = "GlobalHitsAnalyzer_fillECal";
 
   TString eventout;
   if (verbosity > 0)
     eventout = "\nGathering info:";  
   
   // access the calorimeter geometry
   edm::ESHandle<CaloGeometry> theCaloGeometry;
   iSetup.get<CaloGeometryRecord>().get(theCaloGeometry);
   if (!theCaloGeometry.isValid()) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to find CaloGeometryRecord in event!";
     return;
   }
   const CaloGeometry& theCalo(*theCaloGeometry);
     
   // iterator to access containers
   edm::PCaloHitContainer::const_iterator itHit;
 
   // get  ECal information
   edm::PCaloHitContainer theECalHits;
   // extract EB container
   edm::Handle<edm::PCaloHitContainer> EBContainer;
   iEvent.getByToken(ECalEBSrc_Token_,EBContainer);
   if (!EBContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find EcalHitsEB in event!";
     validEB = false;
   }
   // extract EE container
   edm::Handle<edm::PCaloHitContainer> EEContainer;
   iEvent.getByToken(ECalEESrc_Token_,EEContainer);
   if (!EEContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find EcalHitsEE in event!";
     validEE = false;
   }
   // place both containers into new container
   if (validEB)
     theECalHits.insert(theECalHits.end(),EBContainer->begin(),
                EBContainer->end());
   if (validEE)
     theECalHits.insert(theECalHits.end(),EEContainer->begin(),
                EEContainer->end());
 
   // cycle through new container
   int i = 0, j = 0;
   for (itHit = theECalHits.begin(); itHit != theECalHits.end(); ++itHit) {
 
     ++i;
 
     // create a DetId from the detUnitId
     DetId theDetUnitId(itHit->id());
     int detector = theDetUnitId.det();
     int subdetector = theDetUnitId.subdetId();
 
     // check that expected detector is returned
     if ((detector == dEcal) && 
     ((subdetector == sdEcalBrl) ||
      (subdetector == sdEcalFwd))) {
 
       // get the Cell geometry
       auto theDet = (theCalo.getSubdetectorGeometry(theDetUnitId))->getGeometry(theDetUnitId);
 
       if (!theDet) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to get CaloCellGeometry from ECalHits for Hit " << i;
     continue;
       }
 
       ++j;
 
       // get the global position of the cell
       const GlobalPoint& globalposition = theDet->getPosition();
 
       if (meCaloEcalE[0]) meCaloEcalE[0]->Fill(itHit->energy());
       if (meCaloEcalE[1]) meCaloEcalE[1]->Fill(itHit->energy());
       if (meCaloEcalToF[0]) meCaloEcalToF[0]->Fill(itHit->time());
       if (meCaloEcalToF[1]) meCaloEcalToF[1]->Fill(itHit->time());
       if (meCaloEcalPhi) meCaloEcalPhi->Fill(globalposition.phi());
       if (meCaloEcalEta) meCaloEcalEta->Fill(globalposition.eta());
 
     } else {
       edm::LogWarning(MsgLoggerCat)
     << "ECal PCaloHit " << i 
     << " is expected to be (det,subdet) = (" 
     << dEcal << "," << sdEcalBrl
     << " || " << sdEcalFwd << "); value returned is: ("
     << detector << "," << subdetector << ")";
       continue;
     } // end detector type check
   } // end loop through ECal Hits
 
   if (verbosity > 1) {
     eventout += "\n          Number of ECal Hits collected:............. ";
     eventout += j;
   }  
 
   if (meCaloEcal[0]) meCaloEcal[0]->Fill((float)j);
   if (meCaloEcal[1]) meCaloEcal[1]->Fill((float)j); 
 
   // Get Preshower information
   // extract PreShower container
   edm::Handle<edm::PCaloHitContainer> PreShContainer;
   iEvent.getByToken(ECalESSrc_Token_,PreShContainer);
   if (!PreShContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find EcalHitsES in event!";
     validPresh = false;
   }
 
   if (validPresh) {
     // cycle through container
     int i = 0, j = 0;
     for (itHit = PreShContainer->begin(); 
      itHit != PreShContainer->end(); ++itHit) {
       
       ++i;
       
       // create a DetId from the detUnitId
       DetId theDetUnitId(itHit->id());
       int detector = theDetUnitId.det();
       int subdetector = theDetUnitId.subdetId();
       
       // check that expected detector is returned
       if ((detector == dEcal) && 
       (subdetector == sdEcalPS)) {
     
     // get the Cell geometry
     auto theDet = (theCalo.getSubdetectorGeometry(theDetUnitId))->getGeometry(theDetUnitId);
     
     if (!theDet) {
       edm::LogWarning(MsgLoggerCat)
         << "Unable to get CaloCellGeometry from PreShContainer for Hit " 
         << i;
       continue;
     }
     
     ++j;
     
     // get the global position of the cell
     const GlobalPoint& globalposition = theDet->getPosition();
     
     if (meCaloPreShE[0]) meCaloPreShE[0]->Fill(itHit->energy());
     if (meCaloPreShE[1]) meCaloPreShE[1]->Fill(itHit->energy());
     if (meCaloPreShToF[0]) meCaloPreShToF[0]->Fill(itHit->time());
     if (meCaloPreShToF[1]) meCaloPreShToF[1]->Fill(itHit->time());
     if (meCaloPreShPhi) meCaloPreShPhi->Fill(globalposition.phi());
     if (meCaloPreShEta) meCaloPreShEta->Fill(globalposition.eta());
     
       } else {
     edm::LogWarning(MsgLoggerCat)
       << "PreSh PCaloHit " << i 
       << " is expected to be (det,subdet) = (" 
       << dEcal << "," << sdEcalPS
       << "); value returned is: ("
       << detector << "," << subdetector << ")";
     continue;
       } // end detector type check
     } // end loop through PreShower Hits
     
     if (verbosity > 1) {
       eventout += "\n          Number of PreSh Hits collected:............ ";
       eventout += j;
     }  
     
     if (meCaloPreSh[0]) meCaloPreSh[0]->Fill((float)j);
     if (meCaloPreSh[1]) meCaloPreSh[1]->Fill((float)j); 
   }
   
   if (verbosity > 0)
     edm::LogInfo(MsgLoggerCat) << eventout << "\n";
   
   return;
 }

void GlobalHitsAnalyzer::fillG4MC ( const edm::Event & iEvent )

private

Definition at line 767 of file GlobalHitsAnalyzer.cc.

References edm::Provenance::branchDescription(), MonitorElement::Fill(), HcalObjRepresent::Fill(), G4TrkSrc_Token_, G4VtxSrc_Token_, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), edm::Event::getManyByType(), mps_fire::i, edm::HandleBase::isValid(), LogDebug, SiStripPI::max, meGeantTrkE, meGeantTrkPt, meGeantVtxEta, meGeantVtxMulti, meGeantVtxPhi, meGeantVtxRad, meGeantVtxX, meGeantVtxY, meGeantVtxZ, meMCG4Trk, meMCG4Vtx, meMCRGP, trackerHitRTTI::multi, nRawGenPart, edm::HandleBase::provenance(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, csvLumiCalc::unit, validG4trkContainer, validG4VtxContainer, validHepMCevt, verbosity, and vtxunit.

Referenced by analyze().

 {
 
   std::string MsgLoggerCat = "GlobalHitsAnalyzer_fillG4MC";
  
   TString eventout;
   if (verbosity > 0)
     eventout = "\nGathering info:";
 
   // get MC information
   edm::Handle<edm::HepMCProduct> HepMCEvt;
   std::vector<edm::Handle<edm::HepMCProduct> > AllHepMCEvt;
   iEvent.getManyByType(AllHepMCEvt);
 
   // loop through products and extract VtxSmearing if available. Any of them
   // should have the information needed
   for (unsigned int i = 0; i < AllHepMCEvt.size(); ++i) {
     HepMCEvt = AllHepMCEvt[i];
     if ((HepMCEvt.provenance()->branchDescription()).moduleLabel() == "generatorSmeared")
       break;
   }
 
   if (!HepMCEvt.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find HepMCProduct in event!";
     validHepMCevt = false;
   } else {
     eventout += "\n          Using HepMCProduct: ";
     eventout += (HepMCEvt.provenance()->branchDescription()).moduleLabel();
   }
   if (validHepMCevt) {
     const HepMC::GenEvent* MCEvt = HepMCEvt->GetEvent();
     nRawGenPart = MCEvt->particles_size();
     
     
     if (verbosity > 1) {
       eventout += "\n          Number of Raw Particles collected:......... ";
       eventout += nRawGenPart;
     }      
     
     if (meMCRGP[0]) meMCRGP[0]->Fill((float)nRawGenPart);
     if (meMCRGP[1]) meMCRGP[1]->Fill((float)nRawGenPart);  
   }
   
   
   // get G4Vertex information
   // convert unit stored in SimVertex to mm
   float unit = 0.;
   if (vtxunit == 0) unit = 1.;  // already in mm
   if (vtxunit == 1) unit = 10.; // stored in cm, convert to mm
 
   edm::Handle<edm::SimVertexContainer> G4VtxContainer;
   iEvent.getByToken(G4VtxSrc_Token_, G4VtxContainer);
 
   // needed here by vertex multiplicity
   edm::Handle<edm::SimTrackContainer> G4TrkContainer;
   iEvent.getByToken(G4TrkSrc_Token_, G4TrkContainer);
 
   if (!G4VtxContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find SimVertex in event!";
     validG4VtxContainer = false;
   }
   if (validG4VtxContainer) {
     int i = 0;
     edm::SimVertexContainer::const_iterator itVtx;
     for (itVtx = G4VtxContainer->begin(); itVtx != G4VtxContainer->end(); 
      ++itVtx) {
       
       ++i;
       
       const math::XYZTLorentzVector G4Vtx1(itVtx->position().x(),
                        itVtx->position().y(),
                        itVtx->position().z(),
                        itVtx->position().e());
       
       double G4Vtx[4];
       G4Vtx1.GetCoordinates(G4Vtx);
       
       if (meGeantVtxX[0]) meGeantVtxX[0]->Fill((G4Vtx[0]*unit)/micrometer);
       if (meGeantVtxX[1]) meGeantVtxX[1]->Fill((G4Vtx[0]*unit)/micrometer);
       
       if (meGeantVtxY[0]) meGeantVtxY[0]->Fill((G4Vtx[1]*unit)/micrometer);
       if (meGeantVtxY[1]) meGeantVtxY[1]->Fill((G4Vtx[1]*unit)/micrometer);
       
       if (meGeantVtxZ[0]) meGeantVtxZ[0]->Fill((G4Vtx[2]*unit)/millimeter);
       if (meGeantVtxZ[1]) meGeantVtxZ[1]->Fill((G4Vtx[2]*unit)/millimeter); 
 
       if (meGeantVtxEta) meGeantVtxEta->Fill(G4Vtx1.eta());
       if (meGeantVtxPhi) meGeantVtxPhi->Fill(G4Vtx1.phi());
       if (meGeantVtxRad[0]) meGeantVtxRad[0]->Fill(G4Vtx1.rho());
       if (meGeantVtxRad[1]) meGeantVtxRad[1]->Fill(G4Vtx1.rho());
 
       if (meGeantVtxMulti) { 
         int multi = 0;
         if ( G4TrkContainer.isValid() ) {
           edm::SimTrackContainer::const_iterator itTrk;
           for (itTrk = G4TrkContainer->begin(); itTrk != G4TrkContainer->end(); 
                ++itTrk) {
             if ( (*itTrk).vertIndex() == i ) { multi++; }
           }
         }
         meGeantVtxMulti->Fill(((double)multi+0.5));
     }
       
     }
     
     if (verbosity > 1) {
       eventout += "\n          Number of G4Vertices collected:............ ";
       eventout += i;
     }  
     
     if (meMCG4Vtx[0]) meMCG4Vtx[0]->Fill((float)i);
     if (meMCG4Vtx[1]) meMCG4Vtx[1]->Fill((float)i);  
 
   }
 
   // get G4Track information
   if (!G4TrkContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find SimTrack in event!";
     validG4trkContainer = false;
   }
   if (validG4trkContainer) {
     int i = 0;
     edm::SimTrackContainer::const_iterator itTrk;
     for (itTrk = G4TrkContainer->begin(); itTrk != G4TrkContainer->end(); 
      ++itTrk) {
       
       ++i;
       
       const math::XYZTLorentzVector G4Trk1(itTrk->momentum().x(),
                        itTrk->momentum().y(),
                        itTrk->momentum().z(),
                        itTrk->momentum().e());
       double G4Trk[4];
       G4Trk1.GetCoordinates(G4Trk);
       
       if (meGeantTrkPt) meGeantTrkPt->
               Fill(std::log10(std::max(sqrt(G4Trk[0]*G4Trk[0]+G4Trk[1]*G4Trk[1]),-9.)));
       if (meGeantTrkE) meGeantTrkE->Fill(std::log10(std::max(G4Trk[3],-9.)));
     } 
     
     if (verbosity > 1) {
       eventout += "\n          Number of G4Tracks collected:.............. ";
       eventout += i;
     }  
     
     if (meMCG4Trk[0]) meMCG4Trk[0]->Fill((float)i);
     if (meMCG4Trk[1]) meMCG4Trk[1]->Fill((float)i); 
   }
 
   if (verbosity > 0)
     edm::LogInfo(MsgLoggerCat) << eventout << "\n";
     
   return;
 }

void GlobalHitsAnalyzer::fillHCal	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 1830 of file GlobalHitsAnalyzer.cc.

References DetId::det(), gamEcalExtractorBlocks_cff::detector, dHcal, PV3DBase< T, PVType, FrameType >::eta(), MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), HcalGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), HCalSrc_Token_, mps_fire::i, edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), LogDebug, meCaloHcal, meCaloHcalE, meCaloHcalEta, meCaloHcalPhi, meCaloHcalToF, PV3DBase< T, PVType, FrameType >::phi(), HcalHitRelabeller::relabel(), sdHcalBrl, sdHcalEC, sdHcalFwd, sdHcalOut, AlCaHLTBitMon_QueryRunRegistry::string, subdetector, DetId::subdetId(), testNumber, validHcal, and verbosity.

Referenced by analyze().

 {
   std::string MsgLoggerCat = "GlobalHitsAnalyzer_fillHCal";
 
   TString eventout;
   if (verbosity > 0)
     eventout = "\nGathering info:";  
   
   // access the calorimeter geometry
   edm::ESHandle<CaloGeometry> theCaloGeometry;
   iSetup.get<CaloGeometryRecord>().get(theCaloGeometry);
   if (!theCaloGeometry.isValid()) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to find CaloGeometryRecord in event!";
     return;
   }
   const CaloGeometry& theCalo(*theCaloGeometry);
 
   edm::ESHandle<HcalDDDRecConstants> pHRNDC;
   iSetup.get<HcalRecNumberingRecord>().get( pHRNDC );
   const HcalDDDRecConstants* hcons = &(*pHRNDC);
 
   // iterator to access containers
   edm::PCaloHitContainer::const_iterator itHit;
 
   // get  HCal information
   // extract HCal container
   edm::Handle<edm::PCaloHitContainer> HCalContainer;
   iEvent.getByToken(HCalSrc_Token_,HCalContainer);
   if (!HCalContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find HCalHits in event!";
     validHcal = false;
   }
 
   if (validHcal) {
     // cycle through container
     int i = 0, j = 0;
     for (itHit = HCalContainer->begin(); 
      itHit != HCalContainer->end(); ++itHit) {
       
       ++i;
 
       // create a DetId from the detUnitId      
       DetId theDetUnitId;
       unsigned int id_ = itHit->id();
       if(testNumber) theDetUnitId = HcalHitRelabeller::relabel(id_,hcons);
       else theDetUnitId = id_;
 
       int detector = theDetUnitId.det();
       int subdetector = theDetUnitId.subdetId();
       
       // check that expected detector is returned
       if ((detector == dHcal) && 
       ((subdetector == sdHcalBrl) ||
        (subdetector == sdHcalEC) ||
        (subdetector == sdHcalOut) ||
        (subdetector == sdHcalFwd))) {
     
     // get the Cell geometry
     const HcalGeometry *theDet = dynamic_cast<const HcalGeometry*>
       (theCalo.getSubdetectorGeometry(theDetUnitId));
     
     if (!theDet) {
       edm::LogWarning(MsgLoggerCat)
         << "Unable to get HcalGeometry from HCalContainer for Hit " 
         << i;
       continue;
     }
     
     ++j;
     
     // get the global position of the cell
     const GlobalPoint& globalposition = theDet->getPosition(theDetUnitId);
     
     if (meCaloHcalE[0]) meCaloHcalE[0]->Fill(itHit->energy());
     if (meCaloHcalE[1]) meCaloHcalE[1]->Fill(itHit->energy());
     if (meCaloHcalToF[0]) meCaloHcalToF[0]->Fill(itHit->time());
     if (meCaloHcalToF[1]) meCaloHcalToF[1]->Fill(itHit->time());
     if (meCaloHcalPhi) meCaloHcalPhi->Fill(globalposition.phi());
     if (meCaloHcalEta) meCaloHcalEta->Fill(globalposition.eta());
     
       } else {
     edm::LogWarning(MsgLoggerCat)
       << "HCal PCaloHit " << i 
       << " is expected to be (det,subdet) = (" 
       << dHcal << "," << sdHcalBrl
       << " || " << sdHcalEC << " || " << sdHcalOut << " || " << sdHcalFwd
       << "); value returned is: ("
       << detector << "," << subdetector << ")";
     continue;
       } // end detector type check
     } // end loop through HCal Hits
     
     if (verbosity > 1) {
       eventout += "\n          Number of HCal Hits collected:............. ";
       eventout += j;
     }  
     
     if (meCaloHcal[0]) meCaloHcal[0]->Fill((float)j);
     if (meCaloHcal[1]) meCaloHcal[1]->Fill((float)j); 
   }
 
   if (verbosity > 0)
     edm::LogInfo(MsgLoggerCat) << eventout << "\n";
   
   return;
 }

void GlobalHitsAnalyzer::fillMuon	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 1344 of file GlobalHitsAnalyzer.cc.

References gamEcalExtractorBlocks_cff::detector, dMuon, PVValHelper::eta, MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, RPCGeometry::idToDetUnit(), CSCGeometry::idToDetUnit(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), DTGeometry::layer(), LogDebug, meMuon, meMuonCscToF, meMuonCscZ, meMuonDtR, meMuonDtToF, meMuonEta, meMuonPhi, meMuonRpcBR, meMuonRpcBToF, meMuonRpcFToF, meMuonRpcFZ, MuonCscSrc_Token_, MuonDtSrc_Token_, MuonRpcSrc_Token_, nMuonHits, perp(), phi, RPCDetId::region(), sdMuonCSC, sdMuonDT, sdMuonRPC, sdMuonRPCRgnBrl, sdMuonRPCRgnFwdn, sdMuonRPCRgnFwdp, AlCaHLTBitMon_QueryRunRegistry::string, subdetector, GeomDet::surface(), validMuonCSC, validMuonDt, validMuonRPC, verbosity, and z.

Referenced by analyze().

 {
   nMuonHits = 0;
   std::string MsgLoggerCat = "GlobalHitsAnalyzer_fillMuon";
 
   TString eventout;
   if (verbosity > 0)
     eventout = "\nGathering info:";  
 
   // iterator to access containers
   edm::PSimHitContainer::const_iterator itHit;
 
   // access the CSC Muon
   // access the CSC Muon geometry
   edm::ESHandle<CSCGeometry> theCSCGeometry;
   iSetup.get<MuonGeometryRecord>().get(theCSCGeometry);
   if (!theCSCGeometry.isValid()) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to find MuonGeometryRecord for the CSCGeometry in event!";
     return;
   }
   const CSCGeometry& theCSCMuon(*theCSCGeometry);
 
   // get Muon CSC information
   edm::Handle<edm::PSimHitContainer> MuonCSCContainer;
   iEvent.getByToken(MuonCscSrc_Token_,MuonCSCContainer);
   if (!MuonCSCContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find MuonCSCHits in event!";
     validMuonCSC = false;
   }
 
   if (validMuonCSC) {
     // cycle through container
     int i = 0, j = 0;
     for (itHit = MuonCSCContainer->begin(); itHit != MuonCSCContainer->end(); 
      ++itHit) {
       
       ++i;
       
       // create a DetId from the detUnitId
       DetId theDetUnitId(itHit->detUnitId());
       int detector = theDetUnitId.det();
       int subdetector = theDetUnitId.subdetId();
       
       // check that expected detector is returned
       if ((detector == dMuon) && 
       (subdetector == sdMuonCSC)) {
     
     // get the GeomDetUnit from the geometry using theDetUnitID
     const GeomDetUnit *theDet = theCSCMuon.idToDetUnit(theDetUnitId);
     
     if (!theDet) {
       edm::LogWarning(MsgLoggerCat)
         << "Unable to get GeomDetUnit from theCSCMuon for hit " << i;
       continue;
     }
     
     ++j;
     
     // get the Surface of the hit (knows how to go from local <-> global)
     const BoundPlane& bSurface = theDet->surface();
     
     if (meMuonCscToF[0]) meMuonCscToF[0]->Fill(itHit->tof());
     if (meMuonCscToF[1]) meMuonCscToF[1]->Fill(itHit->tof());
     if (meMuonCscZ) 
       meMuonCscZ->Fill(bSurface.toGlobal(itHit->localPosition()).z());
     if (meMuonPhi)
       meMuonPhi->Fill(bSurface.toGlobal(itHit->localPosition()).phi());
     if (meMuonEta)
       meMuonEta->Fill(bSurface.toGlobal(itHit->localPosition()).eta());
     
       } else {
     edm::LogWarning(MsgLoggerCat)
       << "MuonCsc PSimHit " << i 
       << " is expected to be (det,subdet) = (" 
       << dMuon << "," << sdMuonCSC
       << "); value returned is: ("
       << detector << "," << subdetector << ")";
     continue;
       } // end detector type check
     } // end loop through CSC Hits
     
     if (verbosity > 1) {
       eventout += "\n          Number of CSC muon Hits collected:......... ";
       eventout += j;
     }  
 
     nMuonHits += j;
   }
 
   // access the DT Muon
   // access the DT Muon geometry
   edm::ESHandle<DTGeometry> theDTGeometry;
   iSetup.get<MuonGeometryRecord>().get(theDTGeometry);
   if (!theDTGeometry.isValid()) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to find MuonGeometryRecord for the DTGeometry in event!";
     return;
   }
   const DTGeometry& theDTMuon(*theDTGeometry);
 
   // get Muon DT information
   edm::Handle<edm::PSimHitContainer> MuonDtContainer;
   iEvent.getByToken(MuonDtSrc_Token_,MuonDtContainer);
   if (!MuonDtContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find MuonDTHits in event!";
     validMuonDt = false;
   }
 
   if (validMuonDt) {
     // cycle through container
     int i = 0, j = 0;
     for (itHit = MuonDtContainer->begin(); itHit != MuonDtContainer->end(); 
      ++itHit) {
       
       ++i;
       
       // create a DetId from the detUnitId
       DetId theDetUnitId(itHit->detUnitId());
       int detector = theDetUnitId.det();
       int subdetector = theDetUnitId.subdetId();
       
       // check that expected detector is returned
       if ((detector == dMuon) && 
       (subdetector == sdMuonDT)) {
     
     // CSC uses wires and layers rather than the full detID
     // get the wireId
     DTWireId wireId(itHit->detUnitId());
     
     // get the DTLayer from the geometry using the wireID
     const DTLayer *theDet = theDTMuon.layer(wireId.layerId());
     
     if (!theDet) {
       edm::LogWarning(MsgLoggerCat)
         << "Unable to get GeomDetUnit from theDtMuon for hit " << i;
       continue;
     }
     
     ++j;
     
     // get the Surface of the hit (knows how to go from local <-> global)
     const BoundPlane& bSurface = theDet->surface();
     
     if (meMuonDtToF[0]) meMuonDtToF[0]->Fill(itHit->tof());
     if (meMuonDtToF[1]) meMuonDtToF[1]->Fill(itHit->tof());
     if (meMuonDtR) 
       meMuonDtR->Fill(bSurface.toGlobal(itHit->localPosition()).perp());
     if (meMuonPhi)
       meMuonPhi->Fill(bSurface.toGlobal(itHit->localPosition()).phi());
     if (meMuonEta)
       meMuonEta->Fill(bSurface.toGlobal(itHit->localPosition()).eta());
     
       } else {
     edm::LogWarning(MsgLoggerCat)
       << "MuonDt PSimHit " << i 
       << " is expected to be (det,subdet) = (" 
       << dMuon << "," << sdMuonDT
       << "); value returned is: ("
       << detector << "," << subdetector << ")";
     continue;
       } // end detector type check
     } // end loop through DT Hits
     
     if (verbosity > 1) {
       eventout += "\n          Number of DT muon Hits collected:.......... ";
       eventout += j;
     } 
     
     nMuonHits += j;
   }
 
   //int RPCBrl = 0, RPCFwd = 0;
   // access the RPC Muon
   // access the RPC Muon geometry
   edm::ESHandle<RPCGeometry> theRPCGeometry;
   iSetup.get<MuonGeometryRecord>().get(theRPCGeometry);
   if (!theRPCGeometry.isValid()) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to find MuonGeometryRecord for the RPCGeometry in event!";
     return;
   }
   const RPCGeometry& theRPCMuon(*theRPCGeometry);
 
   // get Muon RPC information
   edm::Handle<edm::PSimHitContainer> MuonRPCContainer;
   iEvent.getByToken(MuonRpcSrc_Token_,MuonRPCContainer);
   if (!MuonRPCContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find MuonRPCHits in event!";
     validMuonRPC = false;
   }
 
   if (validMuonRPC) {
     // cycle through container
     int i = 0, j = 0;
     int RPCBrl =0, RPCFwd = 0;
     for (itHit = MuonRPCContainer->begin(); itHit != MuonRPCContainer->end(); 
      ++itHit) {
       
       ++i;
       
       // create a DetID from the detUnitId
       DetId theDetUnitId(itHit->detUnitId());
       int detector = theDetUnitId.det();
       int subdetector = theDetUnitId.subdetId();
       
       // check that expected detector is returned
       if ((detector == dMuon) && 
       (subdetector == sdMuonRPC)) {
     
     // get an RPCDetID from the detUnitID
     RPCDetId RPCId(itHit->detUnitId());      
     
     // find the region of the RPC hit
     int region = RPCId.region();
     
     // get the GeomDetUnit from the geometry using the RPCDetId
     const GeomDetUnit *theDet = theRPCMuon.idToDetUnit(theDetUnitId);
     
     if (!theDet) {
       edm::LogWarning(MsgLoggerCat)
         << "Unable to get GeomDetUnit from theRPCMuon for hit " << i;
       continue;
     }
     
     ++j;
     
     // get the Surface of the hit (knows how to go from local <-> global)
     const BoundPlane& bSurface = theDet->surface();
     
     // gather necessary information
     if ((region == sdMuonRPCRgnFwdp) || (region == sdMuonRPCRgnFwdn)) {
       ++RPCFwd;
       
       if (meMuonRpcFToF[0]) meMuonRpcFToF[0]->Fill(itHit->tof());
       if (meMuonRpcFToF[1]) meMuonRpcFToF[1]->Fill(itHit->tof());
       if (meMuonRpcFZ) 
         meMuonRpcFZ->Fill(bSurface.toGlobal(itHit->localPosition()).z());
       if (meMuonPhi)
         meMuonPhi->Fill(bSurface.toGlobal(itHit->localPosition()).phi());
       if (meMuonEta)
         meMuonEta->Fill(bSurface.toGlobal(itHit->localPosition()).eta());
       
     } else if (region == sdMuonRPCRgnBrl) {
       ++RPCBrl;
       
       if (meMuonRpcBToF[0]) meMuonRpcBToF[0]->Fill(itHit->tof());
       if (meMuonRpcBToF[1]) meMuonRpcBToF[1]->Fill(itHit->tof());
       if (meMuonRpcBR) 
         meMuonRpcBR->Fill(bSurface.toGlobal(itHit->
                         localPosition()).perp());
       if (meMuonPhi)
         meMuonPhi->Fill(bSurface.toGlobal(itHit->localPosition()).phi());
       if (meMuonEta)
         meMuonEta->Fill(bSurface.toGlobal(itHit->localPosition()).eta());
       
     } else {
       edm::LogWarning(MsgLoggerCat)
         << "Invalid region for RPC Muon hit" << i;
       continue;
     } // end check of region
       } else {
     edm::LogWarning(MsgLoggerCat)
       << "MuonRpc PSimHit " << i 
       << " is expected to be (det,subdet) = (" 
         << dMuon << "," << sdMuonRPC
       << "); value returned is: ("
       << detector << "," << subdetector << ")";
     continue;
       } // end detector type check
     } // end loop through RPC Hits
     
     if (verbosity > 1) {
       eventout += "\n          Number of RPC muon Hits collected:......... ";
       eventout += j;
       eventout += "\n                    RPC Barrel muon Hits:............ ";
       eventout += RPCBrl;
       eventout += "\n                    RPC Forward muon Hits:........... ";
       eventout += RPCFwd;    
     }  
     
     nMuonHits += j;
   }
 
   if (meMuon[0]) meMuon[0]->Fill((float)nMuonHits);
   if (meMuon[1]) meMuon[1]->Fill((float)nMuonHits); 
   
   if (verbosity > 0)
     edm::LogInfo(MsgLoggerCat) << eventout << "\n";
   
   return;
 }

void GlobalHitsAnalyzer::fillTrk	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 931 of file GlobalHitsAnalyzer.cc.

References gamEcalExtractorBlocks_cff::detector, dTrk, PVValHelper::eta, MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, TrackerGeometry::idToDetUnit(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), LogDebug, meTrackerPx, meTrackerPxBR, meTrackerPxBToF, meTrackerPxEta, meTrackerPxFToF, meTrackerPxFZ, meTrackerPxPhi, meTrackerSi, meTrackerSiBR, meTrackerSiBToF, meTrackerSiEta, meTrackerSiFToF, meTrackerSiFZ, meTrackerSiPhi, nPxlHits, nSiHits, perp(), phi, PxlBrlHighSrc_Token_, PxlBrlLowSrc_Token_, PxlFwdHighSrc_Token_, PxlFwdLowSrc_Token_, sdPxlBrl, sdPxlFwd, sdSiTEC, sdSiTIB, sdSiTID, sdSiTOB, SiTECHighSrc_Token_, SiTECLowSrc_Token_, SiTIBHighSrc_Token_, SiTIBLowSrc_Token_, SiTIDHighSrc_Token_, SiTIDLowSrc_Token_, SiTOBHighSrc_Token_, SiTOBLowSrc_Token_, AlCaHLTBitMon_QueryRunRegistry::string, subdetector, GeomDet::surface(), validPxlBrlHigh, validPxlBrlLow, validPxlFwdHigh, validPxlFwdLow, validSiTECHigh, validSiTECLow, validSiTIBHigh, validSiTIBLow, validSiTIDHigh, validSiTIDLow, validSiTOBHigh, validSiTOBLow, verbosity, and z.

Referenced by analyze().

 {
 
   nPxlHits = 0;
   std::string MsgLoggerCat = "GlobalHitsAnalyzer_fillTrk";
 
   TString eventout;
   if (verbosity > 0)
     eventout = "\nGathering info:";  
   
   // access the tracker geometry
   edm::ESHandle<TrackerGeometry> theTrackerGeometry;
   iSetup.get<TrackerDigiGeometryRecord>().get(theTrackerGeometry);
   if (!theTrackerGeometry.isValid()) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to find TrackerDigiGeometryRecord in event!";
     return;
   }
   const TrackerGeometry& theTracker(*theTrackerGeometry);
     
   // iterator to access containers
   edm::PSimHitContainer::const_iterator itHit;
 
   // get Pixel Barrel information
   edm::PSimHitContainer thePxlBrlHits;
   // extract low container
   edm::Handle<edm::PSimHitContainer> PxlBrlLowContainer;
   iEvent.getByToken(PxlBrlLowSrc_Token_,PxlBrlLowContainer);
   if (!PxlBrlLowContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find TrackerHitsPixelBarrelLowTof in event!";
     validPxlBrlLow = false;
   }
   // extract high container
   edm::Handle<edm::PSimHitContainer> PxlBrlHighContainer;
   iEvent.getByToken(PxlBrlHighSrc_Token_,PxlBrlHighContainer);
   if (!PxlBrlHighContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find TrackerHitsPixelBarrelHighTof in event!";
     validPxlBrlHigh = false;
   }
   // place both containers into new container
   if (validPxlBrlLow) 
     thePxlBrlHits.insert(thePxlBrlHits.end(),PxlBrlLowContainer->begin(),
              PxlBrlLowContainer->end());
   if(validPxlBrlHigh)
     thePxlBrlHits.insert(thePxlBrlHits.end(),PxlBrlHighContainer->begin(),
              PxlBrlHighContainer->end());
 
   // cycle through new container
   int i = 0, j = 0;
   for (itHit = thePxlBrlHits.begin(); itHit != thePxlBrlHits.end(); ++itHit) {
     
     ++i;
     
     // create a DetId from the detUnitId
     DetId theDetUnitId(itHit->detUnitId());
     int detector = theDetUnitId.det();
     int subdetector = theDetUnitId.subdetId();
     
     // check that expected detector is returned
     if ((detector == dTrk) && (subdetector == sdPxlBrl)) {
       
       // get the GeomDetUnit from the geometry using theDetUnitID
       const GeomDetUnit *theDet = theTracker.idToDetUnit(theDetUnitId);
       
       if (!theDet) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to get GeomDetUnit from PxlBrlHits for Hit " << i;
     continue;
       }
       
       ++j;
       
       // get the Surface of the hit (knows how to go from local <-> global)
       const BoundPlane& bSurface = theDet->surface();
       
       if(meTrackerPxBToF) meTrackerPxBToF->Fill(itHit->tof());
       if(meTrackerPxBR) 
     meTrackerPxBR->Fill(bSurface.toGlobal(itHit->localPosition()).perp());
       if(meTrackerPxPhi) 
     meTrackerPxPhi->Fill(bSurface.toGlobal(itHit->localPosition()).phi());
       if(meTrackerPxEta) 
     meTrackerPxEta->Fill(bSurface.toGlobal(itHit->localPosition()).eta());
       
     } else {
       edm::LogWarning(MsgLoggerCat)
     << "PxlBrl PSimHit " << i 
     << " is expected to be (det,subdet) = (" 
     << dTrk << "," << sdPxlBrl
     << "); value returned is: ("
     << detector << "," << subdetector << ")";
       continue;
     } // end detector type check
   } // end loop through PxlBrl Hits
   
   if (verbosity > 1) {
     eventout += "\n          Number of Pixel Barrel Hits collected:..... ";
     eventout += j;
   }  
   
   nPxlHits += j;
   
   // get Pixel Forward information
   edm::PSimHitContainer thePxlFwdHits;
   // extract low container
   edm::Handle<edm::PSimHitContainer> PxlFwdLowContainer;
   iEvent.getByToken(PxlFwdLowSrc_Token_,PxlFwdLowContainer);
   if (!PxlFwdLowContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find TrackerHitsPixelEndcapLowTof in event!";
     validPxlFwdLow = false;
   }
   // extract high container
   edm::Handle<edm::PSimHitContainer> PxlFwdHighContainer;
   iEvent.getByToken(PxlFwdHighSrc_Token_,PxlFwdHighContainer);
   if (!PxlFwdHighContainer.isValid()) {
     LogDebug("GlobalHitsAnalyzer_fillTrk")
       << "Unable to find TrackerHitsPixelEndcapHighTof in event!";
     validPxlFwdHigh = false;
   }
   // place both containers into new container
   if (validPxlFwdLow)
     thePxlFwdHits.insert(thePxlFwdHits.end(),PxlFwdLowContainer->begin(),
              PxlFwdLowContainer->end());
   if (validPxlFwdHigh)
     thePxlFwdHits.insert(thePxlFwdHits.end(),PxlFwdHighContainer->begin(),
              PxlFwdHighContainer->end());
 
   // cycle through new container
   i = 0; j = 0;
   for (itHit = thePxlFwdHits.begin(); itHit != thePxlFwdHits.end(); ++itHit) {
 
     ++i;
 
     // create a DetId from the detUnitId
     DetId theDetUnitId(itHit->detUnitId());
     int detector = theDetUnitId.det();
     int subdetector = theDetUnitId.subdetId();
 
     // check that expected detector is returned
     if ((detector == dTrk) && (subdetector == sdPxlFwd)) {
 
       // get the GeomDetUnit from the geometry using theDetUnitID
       const GeomDetUnit *theDet = theTracker.idToDetUnit(theDetUnitId);
 
       if (!theDet) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to get GeomDetUnit from PxlFwdHits for Hit " << i;;
     continue;
       }
 
       ++j;
 
       // get the Surface of the hit (knows how to go from local <-> global)
       const BoundPlane& bSurface = theDet->surface();
 
       if(meTrackerPxFToF) meTrackerPxFToF->Fill(itHit->tof());
       if(meTrackerPxFZ) 
     meTrackerPxFZ->Fill(bSurface.toGlobal(itHit->localPosition()).z());
       if(meTrackerPxPhi) 
     meTrackerPxPhi->Fill(bSurface.toGlobal(itHit->localPosition()).phi());
       if(meTrackerPxEta) 
     meTrackerPxEta->Fill(bSurface.toGlobal(itHit->localPosition()).eta());
 
     } else {
       edm::LogWarning(MsgLoggerCat)
     << "PxlFwd PSimHit " << i 
     << " is expected to be (det,subdet) = (" 
     << dTrk << "," << sdPxlFwd
     << "); value returned is: ("
     << detector << "," << subdetector << ")";
       continue;
     } // end detector type check
   } // end loop through PxlFwd Hits
 
   if (verbosity > 1) {
     eventout += "\n          Number of Pixel Forward Hits collected:.... ";
     eventout += j;
   }  
 
   nPxlHits += j;
 
   if (meTrackerPx[0]) meTrackerPx[0]->Fill((float)nPxlHits);
   if (meTrackerPx[1]) meTrackerPx[1]->Fill((float)nPxlHits); 
 
   // get Silicon Barrel information
   nSiHits = 0;
   edm::PSimHitContainer theSiBrlHits;
   // extract TIB low container
   edm::Handle<edm::PSimHitContainer> SiTIBLowContainer;
   iEvent.getByToken(SiTIBLowSrc_Token_,SiTIBLowContainer);
   if (!SiTIBLowContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find TrackerHitsTIBLowTof in event!";
     validSiTIBLow = false;
   }
   // extract TIB high container
   edm::Handle<edm::PSimHitContainer> SiTIBHighContainer;
   iEvent.getByToken(SiTIBHighSrc_Token_,SiTIBHighContainer);
   if (!SiTIBHighContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find TrackerHitsTIBHighTof in event!";
     validSiTIBHigh = false;
   }
   // extract TOB low container
   edm::Handle<edm::PSimHitContainer> SiTOBLowContainer;
   iEvent.getByToken(SiTOBLowSrc_Token_,SiTOBLowContainer);
   if (!SiTOBLowContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find TrackerHitsTOBLowTof in event!";
     validSiTOBLow = false;
   }
   // extract TOB high container
   edm::Handle<edm::PSimHitContainer> SiTOBHighContainer;
   iEvent.getByToken(SiTOBHighSrc_Token_,SiTOBHighContainer);
   if (!SiTOBHighContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find TrackerHitsTOBHighTof in event!";
     validSiTOBHigh = false;
   }
   // place all containers into new container
   if (validSiTIBLow)
     theSiBrlHits.insert(theSiBrlHits.end(),SiTIBLowContainer->begin(),
             SiTIBLowContainer->end());
   if (validSiTIBHigh)
     theSiBrlHits.insert(theSiBrlHits.end(),SiTIBHighContainer->begin(),
             SiTIBHighContainer->end());
   if (validSiTOBLow)
     theSiBrlHits.insert(theSiBrlHits.end(),SiTOBLowContainer->begin(),
             SiTOBLowContainer->end());
   if (validSiTOBHigh)
     theSiBrlHits.insert(theSiBrlHits.end(),SiTOBHighContainer->begin(),
             SiTOBHighContainer->end());
 
   // cycle through new container
   i = 0; j = 0;
   for (itHit = theSiBrlHits.begin(); itHit != theSiBrlHits.end(); ++itHit) {
 
     ++i;
 
     // create a DetId from the detUnitId
     DetId theDetUnitId(itHit->detUnitId());
     int detector = theDetUnitId.det();
     int subdetector = theDetUnitId.subdetId();
 
     // check that expected detector is returned
     if ((detector == dTrk) && 
     ((subdetector == sdSiTIB) ||
      (subdetector == sdSiTOB))) {
 
       // get the GeomDetUnit from the geometry using theDetUnitID
       const GeomDetUnit *theDet = theTracker.idToDetUnit(theDetUnitId);
 
       if (!theDet) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to get GeomDetUnit from SiBrlHits for Hit " << i;
     continue;
       }
 
       ++j;
 
       // get the Surface of the hit (knows how to go from local <-> global)
       const BoundPlane& bSurface = theDet->surface();
 
       if(meTrackerSiBToF) meTrackerSiBToF->Fill(itHit->tof());
       if(meTrackerSiBR) 
     meTrackerSiBR->Fill(bSurface.toGlobal(itHit->localPosition()).perp());
       if(meTrackerSiPhi) 
     meTrackerSiPhi->Fill(bSurface.toGlobal(itHit->localPosition()).phi());
       if(meTrackerSiEta) 
     meTrackerSiEta->Fill(bSurface.toGlobal(itHit->localPosition()).eta());
 
     } else {
       edm::LogWarning(MsgLoggerCat)
     << "SiBrl PSimHit " << i 
     << " is expected to be (det,subdet) = (" 
     << dTrk << "," << sdSiTIB
     << " || " << sdSiTOB << "); value returned is: ("
     << detector << "," << subdetector << ")";
       continue;
     } // end detector type check
   } // end loop through SiBrl Hits
 
   if (verbosity > 1) {
     eventout += "\n          Number of Silicon Barrel Hits collected:... ";
     eventout += j;
   }  
 
   nSiHits += j;
 
   // get Silicon Forward information
   edm::PSimHitContainer theSiFwdHits;
   // extract TID low container
   edm::Handle<edm::PSimHitContainer> SiTIDLowContainer;
   iEvent.getByToken(SiTIDLowSrc_Token_,SiTIDLowContainer);
   if (!SiTIDLowContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find TrackerHitsTIDLowTof in event!";
     validSiTIDLow = false;
   }
   // extract TID high container
   edm::Handle<edm::PSimHitContainer> SiTIDHighContainer;
   iEvent.getByToken(SiTIDHighSrc_Token_,SiTIDHighContainer);
   if (!SiTIDHighContainer.isValid()) {
     LogDebug("GlobalHitsAnalyzer_fillTrk")
       << "Unable to find TrackerHitsTIDHighTof in event!";
     validSiTIDHigh = false;
   }
   // extract TEC low container
   edm::Handle<edm::PSimHitContainer> SiTECLowContainer;
   iEvent.getByToken(SiTECLowSrc_Token_,SiTECLowContainer);
   if (!SiTECLowContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find TrackerHitsTECLowTof in event!";
     validSiTECLow = false;
   }
   // extract TEC high container
   edm::Handle<edm::PSimHitContainer> SiTECHighContainer;
   iEvent.getByToken(SiTECHighSrc_Token_,SiTECHighContainer);
   if (!SiTECHighContainer.isValid()) {
     LogDebug(MsgLoggerCat)
       << "Unable to find TrackerHitsTECHighTof in event!";
     validSiTECHigh = false;
   }
   // place all containers into new container
   if (validSiTIDLow)
     theSiFwdHits.insert(theSiFwdHits.end(),SiTIDLowContainer->begin(),
             SiTIDLowContainer->end());
   if (validSiTIDHigh)
     theSiFwdHits.insert(theSiFwdHits.end(),SiTIDHighContainer->begin(),
             SiTIDHighContainer->end());
   if (validSiTECLow)
     theSiFwdHits.insert(theSiFwdHits.end(),SiTECLowContainer->begin(),
             SiTECLowContainer->end());
   if (validSiTECHigh)
     theSiFwdHits.insert(theSiFwdHits.end(),SiTECHighContainer->begin(),
             SiTECHighContainer->end());
 
   // cycle through container
   i = 0; j = 0;
   for (itHit = theSiFwdHits.begin(); itHit != theSiFwdHits.end(); ++itHit) {
 
     ++i;
 
     // create a DetId from the detUnitId
     DetId theDetUnitId(itHit->detUnitId());
     int detector = theDetUnitId.det();
     int subdetector = theDetUnitId.subdetId();
 
     // check that expected detector is returned 
     if ((detector == dTrk) && 
     ((subdetector == sdSiTID) ||
      (subdetector == sdSiTEC))) {
       
       // get the GeomDetUnit from the geometry using theDetUnitID
       const GeomDetUnit *theDet = theTracker.idToDetUnit(theDetUnitId);
       
       if (!theDet) {
     edm::LogWarning(MsgLoggerCat)
       << "Unable to get GeomDetUnit from SiFwdHits Hit " << i;
     return;
       }
       
       ++j;
 
       // get the Surface of the hit (knows how to go from local <-> global)
       const BoundPlane& bSurface = theDet->surface();
       
       if(meTrackerSiFToF) meTrackerSiFToF->Fill(itHit->tof());
       if(meTrackerSiFZ) 
     meTrackerSiFZ->Fill(bSurface.toGlobal(itHit->localPosition()).z());
       if(meTrackerSiPhi) 
     meTrackerSiPhi->Fill(bSurface.toGlobal(itHit->localPosition()).phi());
       if(meTrackerSiEta) 
     meTrackerSiEta->Fill(bSurface.toGlobal(itHit->localPosition()).eta());
 
     } else {
       edm::LogWarning(MsgLoggerCat)
     << "SiFwd PSimHit " << i 
     << " is expected to be (det,subdet) = (" 
     << dTrk << "," << sdSiTOB
     << " || " << sdSiTEC << "); value returned is: ("
     << detector << "," << subdetector << ")";
       continue;
     } // end check detector type
   } // end loop through SiFwd Hits
 
   if (verbosity > 1) {
     eventout += "\n          Number of Silicon Forward Hits collected:.. ";
     eventout += j;
   }  
 
   nSiHits +=j;
 
   if (meTrackerSi[0]) meTrackerSi[0]->Fill((float)nSiHits);
   if (meTrackerSi[1]) meTrackerSi[1]->Fill((float)nSiHits); 
 
   if (verbosity > 0)
     edm::LogInfo(MsgLoggerCat) << eventout << "\n";
 
   return;
 }

Member Data Documentation

unsigned int GlobalHitsAnalyzer::count

private

Definition at line 254 of file GlobalHitsAnalyzer.h.

Referenced by analyze().

edm::InputTag GlobalHitsAnalyzer::ECalEBSrc_

private

Definition at line 160 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> GlobalHitsAnalyzer::ECalEBSrc_Token_

private

Definition at line 162 of file GlobalHitsAnalyzer.h.

Referenced by fillECal(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::ECalEESrc_

private

Definition at line 161 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> GlobalHitsAnalyzer::ECalEESrc_Token_

private

Definition at line 163 of file GlobalHitsAnalyzer.h.

Referenced by fillECal(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::ECalESSrc_

private

Definition at line 171 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> GlobalHitsAnalyzer::ECalESSrc_Token_

private

Definition at line 172 of file GlobalHitsAnalyzer.h.

Referenced by fillECal(), and GlobalHitsAnalyzer().

std::string GlobalHitsAnalyzer::fName

private

Definition at line 101 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

int GlobalHitsAnalyzer::frequency

private

Definition at line 103 of file GlobalHitsAnalyzer.h.

Referenced by analyze(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::G4TrkSrc_

private

Definition at line 149 of file GlobalHitsAnalyzer.h.

edm::EDGetTokenT<edm::SimTrackContainer> GlobalHitsAnalyzer::G4TrkSrc_Token_

private

Definition at line 151 of file GlobalHitsAnalyzer.h.

Referenced by fillG4MC().

edm::InputTag GlobalHitsAnalyzer::G4VtxSrc_

private

Definition at line 148 of file GlobalHitsAnalyzer.h.

edm::EDGetTokenT<edm::SimVertexContainer> GlobalHitsAnalyzer::G4VtxSrc_Token_

private

Definition at line 150 of file GlobalHitsAnalyzer.h.

Referenced by fillG4MC().

bool GlobalHitsAnalyzer::getAllProvenances

private

Definition at line 106 of file GlobalHitsAnalyzer.h.

Referenced by analyze(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::HCalSrc_

private

Definition at line 181 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> GlobalHitsAnalyzer::HCalSrc_Token_

private

Definition at line 182 of file GlobalHitsAnalyzer.h.

Referenced by fillHCal(), and GlobalHitsAnalyzer().

std::string GlobalHitsAnalyzer::label

private

Definition at line 105 of file GlobalHitsAnalyzer.h.

MonitorElement* GlobalHitsAnalyzer::meCaloEcal[2]

private

Definition at line 155 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillECal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloEcalE[2]

private

Definition at line 156 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillECal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloEcalEta

private

Definition at line 159 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillECal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloEcalPhi

private

Definition at line 158 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillECal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloEcalToF[2]

private

Definition at line 157 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillECal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloHcal[2]

private

Definition at line 176 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillHCal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloHcalE[2]

private

Definition at line 177 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillHCal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloHcalEta

private

Definition at line 180 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillHCal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloHcalPhi

private

Definition at line 179 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillHCal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloHcalToF[2]

private

Definition at line 178 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillHCal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloPreSh[2]

private

Definition at line 166 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillECal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloPreShE[2]

private

Definition at line 167 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillECal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloPreShEta

private

Definition at line 170 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillECal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloPreShPhi

private

Definition at line 169 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillECal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meCaloPreShToF[2]

private

Definition at line 168 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillECal(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meGeantTrkE

private

Definition at line 141 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meGeantTrkPt

private

Definition at line 140 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meGeantVtxEta

private

Definition at line 142 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meGeantVtxMulti

private

Definition at line 145 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meGeantVtxPhi

private

Definition at line 143 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meGeantVtxRad[2]

private

Definition at line 144 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meGeantVtxX[2]

private

Definition at line 136 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meGeantVtxY[2]

private

Definition at line 137 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meGeantVtxZ[2]

private

Definition at line 138 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMCG4Trk[2]

private

Definition at line 139 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMCG4Vtx[2]

private

Definition at line 135 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMCRGP[2]

private

Definition at line 134 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillG4MC(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuon[2]

private

Definition at line 230 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuonCscToF[2]

private

Definition at line 241 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuonCscZ

private

Definition at line 242 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuonDtR

private

Definition at line 237 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuonDtToF[2]

private

Definition at line 236 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuonEta

private

Definition at line 232 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuonPhi

private

Definition at line 231 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuonRpcBR

private

Definition at line 249 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuonRpcBToF[2]

private

Definition at line 248 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuonRpcFToF[2]

private

Definition at line 246 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meMuonRpcFZ

private

Definition at line 247 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillMuon(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerPx[2]

private

Definition at line 187 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerPxBR

private

Definition at line 191 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerPxBToF

private

Definition at line 190 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerPxEta

private

Definition at line 189 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerPxFToF

private

Definition at line 192 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerPxFZ

private

Definition at line 193 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerPxPhi

private

Definition at line 188 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerSi[2]

private

Definition at line 205 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerSiBR

private

Definition at line 209 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerSiBToF

private

Definition at line 208 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerSiEta

private

Definition at line 207 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerSiFToF

private

Definition at line 210 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerSiFZ

private

Definition at line 211 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

MonitorElement* GlobalHitsAnalyzer::meTrackerSiPhi

private

Definition at line 206 of file GlobalHitsAnalyzer.h.

Referenced by bookHistograms(), fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::MuonCscSrc_

private

Definition at line 243 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::MuonCscSrc_Token_

private

Definition at line 244 of file GlobalHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::MuonDtSrc_

private

Definition at line 238 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::MuonDtSrc_Token_

private

Definition at line 239 of file GlobalHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::MuonRpcSrc_

private

Definition at line 250 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::MuonRpcSrc_Token_

private

Definition at line 251 of file GlobalHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalHitsAnalyzer().

int GlobalHitsAnalyzer::nMuonHits

private

Definition at line 233 of file GlobalHitsAnalyzer.h.

Referenced by fillMuon().

int GlobalHitsAnalyzer::nPxlHits

private

Definition at line 186 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk().

int GlobalHitsAnalyzer::nRawGenPart

private

Definition at line 146 of file GlobalHitsAnalyzer.h.

Referenced by fillG4MC().

int GlobalHitsAnalyzer::nSiHits

private

Definition at line 204 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk().

bool GlobalHitsAnalyzer::printProvenanceInfo

private

Definition at line 107 of file GlobalHitsAnalyzer.h.

Referenced by analyze(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::PxlBrlHighSrc_

private

Definition at line 195 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::PxlBrlHighSrc_Token_

private

Definition at line 199 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::PxlBrlLowSrc_

private

Definition at line 194 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::PxlBrlLowSrc_Token_

private

Definition at line 198 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::PxlFwdHighSrc_

private

Definition at line 197 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::PxlFwdHighSrc_Token_

private

Definition at line 201 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::PxlFwdLowSrc_

private

Definition at line 196 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::PxlFwdLowSrc_Token_

private

Definition at line 200 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::SiTECHighSrc_

private

Definition at line 219 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::SiTECHighSrc_Token_

private

Definition at line 227 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::SiTECLowSrc_

private

Definition at line 218 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::SiTECLowSrc_Token_

private

Definition at line 226 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::SiTIBHighSrc_

private

Definition at line 213 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::SiTIBHighSrc_Token_

private

Definition at line 221 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::SiTIBLowSrc_

private

Definition at line 212 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::SiTIBLowSrc_Token_

private

Definition at line 220 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::SiTIDHighSrc_

private

Definition at line 217 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::SiTIDHighSrc_Token_

private

Definition at line 225 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::SiTIDLowSrc_

private

Definition at line 216 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::SiTIDLowSrc_Token_

private

Definition at line 224 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::SiTOBHighSrc_

private

Definition at line 215 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::SiTOBHighSrc_Token_

private

Definition at line 223 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

edm::InputTag GlobalHitsAnalyzer::SiTOBLowSrc_

private

Definition at line 214 of file GlobalHitsAnalyzer.h.

Referenced by GlobalHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalHitsAnalyzer::SiTOBLowSrc_Token_

private

Definition at line 222 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::testNumber

private

Definition at line 108 of file GlobalHitsAnalyzer.h.

Referenced by fillHCal(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validEB

private

Definition at line 128 of file GlobalHitsAnalyzer.h.

Referenced by fillECal(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validEE

private

Definition at line 129 of file GlobalHitsAnalyzer.h.

Referenced by fillECal(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validG4trkContainer

private

Definition at line 112 of file GlobalHitsAnalyzer.h.

Referenced by fillG4MC(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validG4VtxContainer

private

Definition at line 111 of file GlobalHitsAnalyzer.h.

Referenced by fillG4MC(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validHcal

private

Definition at line 131 of file GlobalHitsAnalyzer.h.

Referenced by fillHCal(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validHepMCevt

private

Definition at line 110 of file GlobalHitsAnalyzer.h.

Referenced by fillG4MC(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validMuonCSC

private

Definition at line 125 of file GlobalHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validMuonDt

private

Definition at line 126 of file GlobalHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validMuonRPC

private

Definition at line 127 of file GlobalHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validPresh

private

Definition at line 130 of file GlobalHitsAnalyzer.h.

Referenced by fillECal(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validPxlBrlHigh

private

Definition at line 114 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validPxlBrlLow

private

Definition at line 113 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validPxlFwdHigh

private

Definition at line 116 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validPxlFwdLow

private

Definition at line 115 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validSiTECHigh

private

Definition at line 124 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validSiTECLow

private

Definition at line 123 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validSiTIBHigh

private

Definition at line 118 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validSiTIBLow

private

Definition at line 117 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validSiTIDHigh

private

Definition at line 122 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validSiTIDLow

private

Definition at line 121 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validSiTOBHigh

private

Definition at line 120 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

bool GlobalHitsAnalyzer::validSiTOBLow

private

Definition at line 119 of file GlobalHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalHitsAnalyzer().

int GlobalHitsAnalyzer::verbosity

private

Definition at line 102 of file GlobalHitsAnalyzer.h.

Referenced by analyze(), fillECal(), fillG4MC(), fillHCal(), fillMuon(), fillTrk(), and GlobalHitsAnalyzer().

int GlobalHitsAnalyzer::vtxunit

private

Definition at line 104 of file GlobalHitsAnalyzer.h.

Referenced by fillG4MC(), and GlobalHitsAnalyzer().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation