#include <HSCPValidation/HSCPValidator/src/HSCPValidator.cc>

Inheritance diagram for HSCPValidator:

Public Member Functions
	HSCPValidator (const edm::ParameterSet &)

	~HSCPValidator () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	beginJob () override

void	endJob () override

bool	IncreasedTreshold (const trigger::TriggerEvent &trEv, const edm::InputTag &InputPath, double NewThreshold, double etaCut, int NObjectAboveThreshold, bool averageThreshold)

std::string	intToString (int num)

void	makeGenPlots (const edm::Event &iEvent)

void	makeHLTPlots (const edm::Event &iEvent)

void	makeRecoPlots (const edm::Event &iEvent)

void	makeSimDigiPlotsECAL (const edm::Event &iEvent)

void	makeSimDigiPlotsRPC (const edm::Event &iEvent)

void	makeSimTrackPlots (const edm::Event &iEvent)

Private Attributes
TH1F *	cluSizeDistribution_

edm::EDGetTokenT< edm::ValueMap< reco::DeDxData > >	dEdxTrackToken_

TH2F *	dedxVsp

TH2F *	digiOccupancyMapEB_

TH2F *	digiOccupancyMapEEM_

TH2F *	digiOccupancyMapEEP_

TH1F *	digisEcalNumHistEB_

TH1F *	digisEcalNumHistEE_

bool	doGenPlots_

bool	doHLTPlots_

bool	doRecoPlots_

bool	doSimDigiPlots_

bool	doSimTrackPlots_

edm::EDGetTokenT< EBDigiCollection >	EBDigiCollectionToken_

edm::EDGetTokenT< edm::PCaloHitContainer >	ebSimHitToken_

edm::EDGetTokenT< EEDigiCollection >	EEDigiCollectionToken_

edm::EDGetTokenT< edm::PCaloHitContainer >	eeSimHitToken_

TH1F *	efficiencyRPCRecHitSimDigis_

TH1F *	h_genhscp_met

TH1F *	h_genhscp_met_nohscp

TH1F *	h_genhscp_scaloret

TH1F *	h_genhscp_scaloret_nohscp

TH1F *	hltjet

TH1F *	hltmet

TH1F *	hltmu

TH1F *	particleBetaHist_

TH1F *	particleBetaInverseHist_

TH1F *	particleEtaHist_

std::vector< int >	particleIds_

std::map< int, int >	particleIdsFoundMap_

TH1F *	particleMassHist_

TH1F *	particlePhiHist_

TH1F *	particlePHist_

TH1F *	particlePtHist_

int	particleStatus_

TH1F *	particleStatusHist_

TH2F *	RecoHSCPPtVsGenPt

TH1F *	residualsRPCRecHitSimDigis_

TH1F *	rpcBXBarrel_ [6]

TH1F *	rpcBXEndCap_ [3]

edm::ESHandle< RPCGeometry >	rpcGeo

edm::EDGetTokenT< RPCRecHitCollection >	rpcRecHitsToken_

TH1F *	rpcTimeOfFlightBarrel_ [6]

TH1F *	rpcTimeOfFlightEndCap_ [3]

TH1F *	simHitsEcalDigiMatchEnergyHistEB_

TH1F *	simHitsEcalDigiMatchEnergyHistEE_

TH2F *	simHitsEcalDigiMatchEnergyVsTimeHistEB_

TH2F *	simHitsEcalDigiMatchEnergyVsTimeHistEE_

TH1F *	simHitsEcalDigiMatchIEtaHist_

TH1F *	simHitsEcalDigiMatchIPhiHist_

TH1F *	simHitsEcalDigiMatchTimeHistEB_

TH1F *	simHitsEcalDigiMatchTimeHistEE_

TH1F *	simHitsEcalEnergyHistEB_

TH1F *	simHitsEcalEnergyHistEE_

TH2F *	simHitsEcalEnergyVsTimeHistEB_

TH2F *	simHitsEcalEnergyVsTimeHistEE_

TH1F *	simHitsEcalNumHistEB_

TH1F *	simHitsEcalNumHistEE_

TH1F *	simHitsEcalTimeHistEB_

TH1F *	simHitsEcalTimeHistEE_

TH1F *	simTrackParticleBetaHist_

TH1F *	simTrackParticleEtaHist_

TH1F *	simTrackParticlePhiHist_

TH1F *	simTrackParticlePHist_

TH1F *	simTrackParticlePtHist_

edm::EDGetTokenT< edm::SimTrackContainer >	simTracksToken_

edm::EDGetTokenT< edm::SimTrackContainer >	simTrackToken_

edm::EDGetTokenT< reco::TrackCollection >	tkTracksToken_

edm::EDGetTokenT< edm::HepMCProduct >	token_

edm::EDGetTokenT< trigger::TriggerEvent >	trEvToken_

edm::EDGetTokenT< edm::TriggerResults >	triggerResultsToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 46 of file HSCPValidator.h.

Constructor & Destructor Documentation

HSCPValidator::HSCPValidator ( const edm::ParameterSet & iConfig )

explicit

Definition at line 96 of file HSCPValidator.cc.

     : doGenPlots_(iConfig.getParameter<bool>("MakeGenPlots")),
       doHLTPlots_(iConfig.getParameter<bool>("MakeHLTPlots")),
       doSimTrackPlots_(iConfig.getParameter<bool>("MakeSimTrackPlots")),
       doSimDigiPlots_(iConfig.getParameter<bool>("MakeSimDigiPlots")),
       doRecoPlots_(iConfig.getParameter<bool>("MakeRecoPlots")),
       token_(consumes<edm::HepMCProduct>(iConfig.getParameter<edm::InputTag>("generatorLabel"))),
       simTracksToken_(consumes<edm::SimTrackContainer>(edm::InputTag("g4SimHits"))),
       trEvToken_(consumes<trigger::TriggerEvent>(edm::InputTag("hltTriggerSummaryAOD"))),
       tkTracksToken_(consumes<reco::TrackCollection>(edm::InputTag("generalTracks"))),
       dEdxTrackToken_(consumes<edm::ValueMap<reco::DeDxData> >(edm::InputTag("dedxHarmonic2"))),
       rpcRecHitsToken_(consumes<RPCRecHitCollection>(edm::InputTag("rpcRecHits"))),
       triggerResultsToken_(consumes<edm::TriggerResults>(edm::InputTag("TriggerResults", "", "HLT"))),
       particleIds_(iConfig.getParameter<std::vector<int> >("particleIds")),
       particleStatus_(iConfig.getUntrackedParameter<int>("particleStatus", 1)),
       ebSimHitToken_(consumes<edm::PCaloHitContainer>(iConfig.getParameter<edm::InputTag>("EBSimHitCollection"))),
       eeSimHitToken_(consumes<edm::PCaloHitContainer>(iConfig.getParameter<edm::InputTag>("EESimHitCollection"))),
       simTrackToken_(consumes<edm::SimTrackContainer>(iConfig.getParameter<edm::InputTag>("SimTrackCollection"))),
       EBDigiCollectionToken_(consumes<EBDigiCollection>(iConfig.getParameter<edm::InputTag>("EBDigiCollection"))),
       EEDigiCollectionToken_(consumes<EEDigiCollection>(iConfig.getParameter<edm::InputTag>("EEDigiCollection"))) {
   //now do what ever initialization is needed
   // GEN
   particleEtaHist_ = fileService->make<TH1F>("particleEta", "Eta of gen particle", 100, -5, 5);
   particlePhiHist_ = fileService->make<TH1F>("particlePhi", "Phi of gen particle", 180, -3.15, 3.15);
   particlePHist_ = fileService->make<TH1F>("particleP", "Momentum of gen particle", 500, 0, 2000);
   particlePtHist_ = fileService->make<TH1F>("particlePt", "P_{T} of gen particle", 500, 0, 2000);
   particleMassHist_ = fileService->make<TH1F>("particleMass", "Mass of gen particle", 1000, 0, 2000);
   particleStatusHist_ = fileService->make<TH1F>("particleStatus", "Status of gen particle", 10, 0, 10);
   particleBetaHist_ = fileService->make<TH1F>("particleBeta", "Beta of gen particle", 100, 0, 1);
   particleBetaInverseHist_ = fileService->make<TH1F>("particleBetaInverse", "1/#beta of gen particle", 100, 0, 5);
 
   h_genhscp_met = fileService->make<TH1F>("hscp_met", "missing E_{T} hscp", 100, 0., 1500.);
   h_genhscp_met_nohscp = fileService->make<TH1F>("hscp_met_nohscp", "missing E_{T} w/o hscp", 100, 0., 1500.);
   h_genhscp_scaloret = fileService->make<TH1F>("hscp_scaloret", "scalor E_{T} sum", 100, 0., 1500.);
   h_genhscp_scaloret_nohscp =
       fileService->make<TH1F>("hscp_scaloret_nohscp", "scalor E_{T} sum w/o hscp", 100, 0., 1500.);
 
   //SIM track Info
   simTrackParticleEtaHist_ = fileService->make<TH1F>("simTrackParticleEta", "Eta of simTrackParticle", 100, -5, 5);
   simTrackParticlePhiHist_ =
       fileService->make<TH1F>("simTrackParticlePhi", "Phi of simTrackParticle", 180, -3.15, 3.15);
   simTrackParticlePHist_ = fileService->make<TH1F>("simTrackParticleP", "Momentum of simTrackParticle", 500, 0, 2000);
   simTrackParticlePtHist_ = fileService->make<TH1F>("simTrackParticlePt", "P_{T} of simTrackParticle", 500, 0, 2000);
   simTrackParticleBetaHist_ = fileService->make<TH1F>("simTrackParticleBeta", "Beta of simTrackParticle", 100, 0, 1);
   //reco track Info
 
   RecoHSCPPtVsGenPt = fileService->make<TH2F>("Recovsgenpt", "RecovsGen", 100, 0, 1000, 100, 0, 1000);
   dedxVsp = fileService->make<TH2F>("dedxvsp", "dedxvsp", 100, 0, 1000, 100, 0, 10);
   //HLT Info
   hltmet = fileService->make<TH1F>("HLT_MET", "MET", 3, -1, 2);
   hltjet = fileService->make<TH1F>("HLT_JET", "JET", 3, -1, 2);
   hltmu = fileService->make<TH1F>("HLT_Mu", "Mu", 3, -1, 2);
 
   // SIM-DIGI: ECAL
   simHitsEcalEnergyHistEB_ =
       fileService->make<TH1F>("ecalEnergyOfSimHitsEB", "HSCP SimTrack-matching SimHit energy EB [GeV]", 125, -1, 4);
   simHitsEcalEnergyHistEE_ =
       fileService->make<TH1F>("ecalEnergyOfSimHitsEE", "HSCP SimTrack-matching SimHit energy EE [GeV]", 125, -1, 4);
   simHitsEcalTimeHistEB_ =
       fileService->make<TH1F>("ecalTimingOfSimHitsEB", "HSCP SimTrack-matching SimHit time EB [ns]", 115, -15, 100);
   simHitsEcalTimeHistEE_ =
       fileService->make<TH1F>("ecalTimingOfSimHitsEE", "HSCP SimTrack-matching SimHit time EE [ns]", 115, -15, 100);
   simHitsEcalNumHistEB_ = fileService->make<TH1F>(
       "ecalNumberOfSimHitsEB", "Number of HSCP SimTrack-matching EB sim hits in event", 100, 0, 200);
   simHitsEcalNumHistEE_ = fileService->make<TH1F>(
       "ecalNumberOfSimHitsEE", "Number of HSCP SimTrack-matching EE sim hits in event", 100, 0, 200);
   simHitsEcalEnergyVsTimeHistEB_ =
       fileService->make<TH2F>("ecalEnergyVsTimeOfSimHitsEB",
                               "Energy vs. time of HSCP SimTrack-matching EB sim hits in event",
                               115,
                               -15,
                               100,
                               125,
                               -1,
                               4);
   simHitsEcalEnergyVsTimeHistEE_ =
       fileService->make<TH2F>("ecalEnergyVsTimeOfSimHitsEE",
                               "Energy vs. time of HSCP SimTrack-matching EE sim hits in event",
                               115,
                               -15,
                               100,
                               125,
                               -1,
                               4);
   simHitsEcalDigiMatchEnergyHistEB_ =
       fileService->make<TH1F>("ecalEnergyOfDigiMatSimHitsEB", "HSCP digi-matching SimHit energy EB [GeV]", 125, -1, 4);
   simHitsEcalDigiMatchEnergyHistEE_ =
       fileService->make<TH1F>("ecalEnergyOfDigiMatSimHitsEE", "HSCP digi-matching SimHit energy EE [GeV]", 125, -1, 4);
   simHitsEcalDigiMatchTimeHistEB_ =
       fileService->make<TH1F>("ecalTimingOfDigiMatSimHitsEB", "HSCP digi-matching SimHit time EB [ns]", 115, -15, 100);
   simHitsEcalDigiMatchTimeHistEE_ =
       fileService->make<TH1F>("ecalTimingOfDigiMatSimHitsEE", "HSCP digi-matching SimHit time EE [ns]", 115, -15, 100);
   simHitsEcalDigiMatchEnergyVsTimeHistEB_ = fileService->make<TH2F>(
       "ecalEnergyVsTimeOfDigiMatSimHitsEB", "HSCP digi-matching EB SimHit energy vs. time", 115, -15, 100, 125, -1, 4);
   simHitsEcalDigiMatchEnergyVsTimeHistEE_ = fileService->make<TH2F>(
       "ecalEnergyVsTimeOfDigiMatSimHitsEE", "HSCP digi-matching EE SimHit energy vs. time", 115, -15, 100, 125, -1, 4);
   simHitsEcalDigiMatchIEtaHist_ =
       fileService->make<TH1F>("ecalIEtaOfDigiMatchSimHits", "iEta of digi-matching Ecal simHits (EB)", 171, -85, 86);
   simHitsEcalDigiMatchIPhiHist_ =
       fileService->make<TH1F>("ecalIPhiOfDigiMatchSimHits", "iPhi of digi-matching Ecal simHits (EB)", 360, 1, 361);
   digisEcalNumHistEB_ =
       fileService->make<TH1F>("ecalDigisNumberEB", "Number of EB digis matching simhits in event", 200, 0, 1000);
   digisEcalNumHistEE_ =
       fileService->make<TH1F>("ecalDigisNumberEE", "Number of EE digis matching simhits in event", 200, 0, 1000);
   digiOccupancyMapEB_ = fileService->make<TH2F>(
       "ecalDigiOccupancyMapEB", "Occupancy of simhit-matching digis EB;i#phi;i#eta", 360, 1, 361, 171, -85, 86);
   digiOccupancyMapEEP_ = fileService->make<TH2F>(
       "ecalDigiOccupancyMapEEM", "Occupancy of simhit-matching digis EEM;ix;iy", 100, 1, 100, 100, 1, 100);
   digiOccupancyMapEEM_ = fileService->make<TH2F>(
       "ecalDigiOccupancyMapEEP", "Occupancy of simhit-matching digis EEP;ix;iy", 100, 1, 100, 100, 1, 100);
 
   // SIM-DIGI: RPC
   residualsRPCRecHitSimDigis_ =
       fileService->make<TH1F>("residualsRPCRecHitSimDigis", "HSCP SimHit - Clossest RPC RecHit", 100, -5, 5);
   efficiencyRPCRecHitSimDigis_ =
       fileService->make<TH1F>("efficiencyRPCRecHitSimDigis", "HSCP SimHits RecHits Efficiency", 2, -0.5, 1.5);
   cluSizeDistribution_ = fileService->make<TH1F>("RPCCluSizeDistro", "RPC HSCP CluSize Distribution", 11, -0.5, 10.5);
   rpcTimeOfFlightBarrel_[0] = fileService->make<TH1F>("RPCToFLayer1", "RPC HSCP Time Of Flight Layer 1", 50, 5, 100);
   rpcTimeOfFlightBarrel_[1] = fileService->make<TH1F>("RPCToFLayer2", "RPC HSCP Time Of Flight Layer 2", 50, 5, 100);
   rpcTimeOfFlightBarrel_[2] = fileService->make<TH1F>("RPCToFLayer3", "RPC HSCP Time Of Flight Layer 3", 50, 5, 100);
   rpcTimeOfFlightBarrel_[3] = fileService->make<TH1F>("RPCToFLayer4", "RPC HSCP Time Of Flight Layer 4", 50, 5, 100);
   rpcTimeOfFlightBarrel_[4] = fileService->make<TH1F>("RPCToFLayer5", "RPC HSCP Time Of Flight Layer 5", 50, 5, 100);
   rpcTimeOfFlightBarrel_[5] = fileService->make<TH1F>("RPCToFLayer6", "RPC HSCP Time Of Flight Layer 6", 50, 5, 100);
   rpcBXBarrel_[0] = fileService->make<TH1F>("RPCBXLayer1", "RPC HSCP BX Layer 1", 5, -0.5, 4.5);
   rpcBXBarrel_[1] = fileService->make<TH1F>("RPCBXLayer2", "RPC HSCP BX Layer 2", 5, -0.5, 4.5);
   rpcBXBarrel_[2] = fileService->make<TH1F>("RPCBXLayer3", "RPC HSCP BX Layer 3", 5, -0.5, 4.5);
   rpcBXBarrel_[3] = fileService->make<TH1F>("RPCBXLayer4", "RPC HSCP BX Layer 4", 5, -0.5, 4.5);
   rpcBXBarrel_[4] = fileService->make<TH1F>("RPCBXLayer5", "RPC HSCP BX Layer 5", 5, -0.5, 4.5);
   rpcBXBarrel_[5] = fileService->make<TH1F>("RPCBXLayer6", "RPC HSCP BX Layer 6", 5, -0.5, 4.5);
   rpcTimeOfFlightEndCap_[0] = fileService->make<TH1F>("RPCToFDisk1", "RPC HSCP Time Of Flight Disk 1", 50, 5, 100);
   rpcTimeOfFlightEndCap_[1] = fileService->make<TH1F>("RPCToFDisk2", "RPC HSCP Time Of Flight Disk 2", 50, 5, 100);
   rpcTimeOfFlightEndCap_[2] = fileService->make<TH1F>("RPCToFDisk3", "RPC HSCP Time Of Flight Disk 3", 50, 5, 100);
   rpcBXEndCap_[0] = fileService->make<TH1F>("RPCBXDisk1", "RPC HSCP BX Disk 1", 5, -0.5, 4.5);
   rpcBXEndCap_[1] = fileService->make<TH1F>("RPCBXDisk2", "RPC HSCP BX Disk 2", 5, -0.5, 4.5);
   rpcBXEndCap_[2] = fileService->make<TH1F>("RPCBXDisk3", "RPC HSCP BX Disk 3", 5, -0.5, 4.5);
 }

HSCPValidator::~HSCPValidator ( )

override

Definition at line 233 of file HSCPValidator.cc.

                               {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   //   particleEtaHist_ = fileService->make<TH1F>("particleEta","Eta of gen particle",100,-5,5);
   //   particlePhiHist_ = fileService->make<TH1F>("particlePhi","Phi of gen particle",180,-3.15,3.15);
   //   particlePHist_ = fileService->make<TH1F>("particleP","Momentum of gen particle",500,0,2000);
   //   particlePtHist_ = fileService->make<TH1F>("particlePt","P_{T} of gen particle",500,0,2000);
   //   particleMassHist_ = fileService->make<TH1F>("particleMass","Mass of gen particle",1000,0,2000);
   //   particleStatusHist_ = fileService->make<TH1F>("particleStatus","Status of gen particle",10,0,10);
   //   particleBetaHist_ = fileService->make<TH1F>("particleBeta","Beta of gen particle",100,0,1);
   //   particleBetaInverseHist_ = fileService->make<TH1F>("particleBetaInverse","1/#beta of gen particle",100,0,5);
 }

Member Function Documentation

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

overrideprivate

Definition at line 251 of file HSCPValidator.cc.

References doGenPlots_, doHLTPlots_, doRecoPlots_, doSimDigiPlots_, doSimTrackPlots_, edm::EventSetup::get(), makeGenPlots(), makeHLTPlots(), makeRecoPlots(), makeSimDigiPlotsECAL(), makeSimDigiPlotsRPC(), makeSimTrackPlots(), and rpcGeo.

                                                                                {
   using namespace edm;
   iSetup.get<MuonGeometryRecord>().get(rpcGeo);
 
   if (doGenPlots_)
     makeGenPlots(iEvent);
   if (doSimTrackPlots_)
     makeSimTrackPlots(iEvent);
   if (doSimDigiPlots_) {
     makeSimDigiPlotsECAL(iEvent);
     makeSimDigiPlotsRPC(iEvent);
   }
   if (doHLTPlots_) {
     makeHLTPlots(iEvent);
   }
   if (doRecoPlots_) {
     makeRecoPlots(iEvent);
   }
 }

void HSCPValidator::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 272 of file HSCPValidator.cc.

272 {}

void HSCPValidator::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 275 of file HSCPValidator.cc.

References gather_cfg::cout, intToString(), particleIdsFoundMap_, and AlCaHLTBitMon_QueryRunRegistry::string.

                            {
   std::string frequencies = "";
   for (std::map<int, int>::const_iterator itr = particleIdsFoundMap_.begin(); itr != particleIdsFoundMap_.end();
        ++itr) {
     frequencies += "PDG ID: ";
     frequencies += intToString(itr->first);
     frequencies += " Frequency: ";
     frequencies += intToString(itr->second);
     frequencies += "\n";
   }
   std::cout << "Found PDGIds: "
             << "\n\n"
             << frequencies << std::endl;
 }

bool HSCPValidator::IncreasedTreshold	(	const trigger::TriggerEvent &	trEv,
		const edm::InputTag &	InputPath,
		double	NewThreshold,
		double	etaCut,
		int	NObjectAboveThreshold,
		bool	averageThreshold
	)

private

Definition at line 839 of file HSCPValidator.cc.

References Average, DEFINE_FWK_MODULE, PVValHelper::eta, fftjetproducer_cfi::etaCut, trigger::TriggerEvent::filterIds(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::getObjects(), mps_fire::i, SiStripPI::max, dqmiodumpmetadata::n, DiDispStaMuonMonitor_cfi::pt, and trigger::TriggerEvent::sizeFilters().

Referenced by makeHLTPlots().

                                                              {
   unsigned int filterIndex = trEv.filterIndex(InputPath);
   //if(filterIndex<trEv.sizeFilters())printf("SELECTED INDEX =%i --> %s    XXX   %s\n",filterIndex,trEv.filterTag(filterIndex).label().c_str(), trEv.filterTag(filterIndex).process().c_str());
 
   if (filterIndex < trEv.sizeFilters()) {
     const trigger::Vids& VIDS(trEv.filterIds(filterIndex));
     const trigger::Keys& KEYS(trEv.filterKeys(filterIndex));
     const int nI(VIDS.size());
     const int nK(KEYS.size());
     assert(nI == nK);
     const int n(std::max(nI, nK));
     const trigger::TriggerObjectCollection& TOC(trEv.getObjects());
 
     if (!averageThreshold) {
       int NObjectAboveThresholdObserved = 0;
       for (int i = 0; i != n; ++i) {
         if (TOC[KEYS[i]].pt() > NewThreshold && fabs(TOC[KEYS[i]].eta()) < etaCut)
           NObjectAboveThresholdObserved++;
         //cout << "   " << i << " " << VIDS[i] << "/" << KEYS[i] << ": "<< TOC[KEYS[i]].id() << " " << TOC[KEYS[i]].pt() << " " << TOC[KEYS[i]].eta() << " " << TOC[KEYS[i]].phi() << " " << TOC[KEYS[i]].mass()<< endl;
       }
       if (NObjectAboveThresholdObserved >= NObjectAboveThreshold)
         return true;
 
     } else {
       std::vector<double> ObjPt;
 
       for (int i = 0; i != n; ++i) {
         ObjPt.push_back(TOC[KEYS[i]].pt());
         //cout << "   " << i << " " << VIDS[i] << "/" << KEYS[i] << ": "<< TOC[KEYS[i]].id() << " " << TOC[KEYS[i]].pt() << " " << TOC[KEYS[i]].eta() << " " << TOC[KEYS[i]].phi() << " " << TOC[KEYS[i]].mass()<< endl;
       }
       if ((int)(ObjPt.size()) < NObjectAboveThreshold)
         return false;
       std::sort(ObjPt.begin(), ObjPt.end());
 
       double Average = 0;
       for (int i = 0; i < NObjectAboveThreshold; i++) {
         Average += ObjPt[ObjPt.size() - 1 - i];
       }
       Average /= NObjectAboveThreshold;
       //cout << "AVERAGE = " << Average << endl;
 
       if (Average > NewThreshold)
         return true;
     }
   }
   return false;
 }

std::string HSCPValidator::intToString ( int num )

private

Definition at line 830 of file HSCPValidator.cc.

Referenced by endJob().

                                             {
   using namespace std;
   ostringstream myStream;
   myStream << num << flush;
   return (myStream.str());  //returns the string form of the stringstream object
 }

void HSCPValidator::makeGenPlots ( const edm::Event & iEvent )

private

Definition at line 291 of file HSCPValidator.cc.

References funct::abs(), spr::find(), edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), h_genhscp_met, h_genhscp_met_nohscp, h_genhscp_scaloret, h_genhscp_scaloret_nohscp, mag(), AlCaHLTBitMon_ParallelJobs::p, particleBetaHist_, particleBetaInverseHist_, particleEtaHist_, particleIds_, particleIdsFoundMap_, particleMassHist_, particlePhiHist_, particlePHist_, particlePtHist_, particleStatus_, particleStatusHist_, funct::pow(), mathSSE::sqrt(), and token_.

Referenced by analyze().

                                                        {
   using namespace edm;
 
   double missingpx = 0;
   double missingpy = 0;
   double missingpx_nohscp = 0;
   double missingpy_nohscp = 0;
   double scalorEt = 0;
   double scalorEt_nohscp = 0;
 
   Handle<HepMCProduct> evt;
   iEvent.getByToken(token_, evt);
 
   HepMC::GenEvent* myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
   for (HepMC::GenEvent::particle_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end(); ++p) {
     if ((*p)->status() != particleStatus_)
       continue;
     //calculate MET(neutrino as MET)
     if (abs((*p)->pdg_id()) != 12 && abs((*p)->pdg_id()) != 14 &&
         abs((*p)->pdg_id()) != 16) {  //for non-neutrino particles.
       missingpx -= (*p)->momentum().px();
       missingpy -= (*p)->momentum().py();
       scalorEt += (*p)->momentum().perp();
     }
 
     // Check if the particleId is in our R-hadron list
     std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(), particleIds_.end(), (*p)->pdg_id());
     if (partIdItr == particleIds_.end()) {
       //calculate MET(neutrino+ HSCP as MET)
       if (abs((*p)->pdg_id()) != 12 && abs((*p)->pdg_id()) != 14 &&
           abs((*p)->pdg_id()) != 16) {  //for non-neutrino particles.
         missingpx_nohscp -= (*p)->momentum().px();
         missingpy_nohscp -= (*p)->momentum().py();
         scalorEt_nohscp += (*p)->momentum().perp();
       }
     } else {
       particleStatusHist_->Fill((*p)->status());
 
       std::pair<std::map<int, int>::iterator, bool> pair =
           particleIdsFoundMap_.insert(std::make_pair<int, int>((*p)->pdg_id(), 1));
       if (!pair.second) {
         ++(pair.first->second);
       }
 
       double mag = sqrt(pow((*p)->momentum().px(), 2) + pow((*p)->momentum().py(), 2) + pow((*p)->momentum().pz(), 2));
       particleEtaHist_->Fill((*p)->momentum().eta());
       particlePhiHist_->Fill((*p)->momentum().phi());
       particlePHist_->Fill(mag);
       particlePtHist_->Fill((*p)->momentum().perp());
       particleMassHist_->Fill((*p)->generated_mass());
       float particleP = mag;
       float particleM = (*p)->generated_mass();
       particleBetaHist_->Fill(particleP / sqrt(particleP * particleP + particleM * particleM));
       particleBetaInverseHist_->Fill(sqrt(particleP * particleP + particleM * particleM) / particleP);
     }
   }
 
   h_genhscp_met->Fill(sqrt(missingpx * missingpx + missingpy * missingpy));
   h_genhscp_met_nohscp->Fill(sqrt(missingpx_nohscp * missingpx_nohscp + missingpy_nohscp * missingpy_nohscp));
   h_genhscp_scaloret->Fill(scalorEt);
   h_genhscp_scaloret_nohscp->Fill(scalorEt_nohscp);
 
   delete myGenEvent;
 }

void HSCPValidator::makeHLTPlots ( const edm::Event & iEvent )

private

Definition at line 384 of file HSCPValidator.cc.

References edm::TriggerResultsByName::accept(), gather_cfg::cout, edm::Event::getByToken(), hltjet, hltmet, hltmu, mps_fire::i, IncreasedTreshold(), HLT_2018_cff::InputTag, edm::TriggerResultsByName::isValid(), edm::HandleBase::isValid(), edm::TriggerResultsByName::size(), trEvToken_, edm::TriggerResultsByName::triggerIndex(), edm::TriggerResultsByName::triggerName(), triggerResults, edm::Event::triggerResultsByName(), and triggerResultsToken_.

Referenced by analyze().

                                                        {
   using namespace edm;
   //get HLT infos
 
   edm::Handle<edm::TriggerResults> triggerResults;
   iEvent.getByToken(triggerResultsToken_, triggerResults);
 
   edm::TriggerResultsByName tr(nullptr, nullptr);
   if (triggerResults.isValid()) {
     tr = iEvent.triggerResultsByName(*triggerResults);
   }
 
   if (!tr.isValid()) {
     std::cout << "Trigger Results not available" << std::endl;
   }
 
   edm::Handle<trigger::TriggerEvent> trEvHandle;
   iEvent.getByToken(trEvToken_, trEvHandle);
   trigger::TriggerEvent trEv = *trEvHandle;
 
   unsigned int TrIndex_Unknown = tr.size();
 
   // HLT TRIGGER BASED ON 1 MUON!
   if (TrIndex_Unknown != tr.triggerIndex("HLT_Mu40_v1")) {
     if (tr.accept(tr.triggerIndex("HLT_Mu40_v1")))
       hltmu->Fill(1);
     else {
       hltmu->Fill(0);
     }
   } else {
     if (TrIndex_Unknown != tr.triggerIndex("HLT_Mu30_v1")) {
       if (IncreasedTreshold(trEv, InputTag("hltSingleMu30L3Filtered30", "", "HLT"), 40, 2.1, 1, false))
         hltmu->Fill(1);
       else
         hltmu->Fill(0);
     } else {
       printf("BUG with HLT_Mu\n");
       std::cout << "trigger names are : ";
       for (unsigned int i = 0; i < tr.size(); i++) {
         std::cout << " " << tr.triggerName(i);
       }
       std::cout << std::endl;
     }
   }
 
   // HLT TRIGGER BASED ON MET!
   if (TrIndex_Unknown != tr.triggerIndex("HLT_PFMHT150_v3")) {
     if (tr.accept(tr.triggerIndex("HLT_PFMHT150_v3")))
       hltmet->Fill(1);
     else
       hltmet->Fill(0);
   } else {
     if (TrIndex_Unknown != tr.triggerIndex("HLT_PFMHT150_v2")) {
       if (tr.accept(tr.triggerIndex("HLT_PFMHT150_v2")))
         hltmet->Fill(1);
       else
         hltmet->Fill(0);
     } else {
       printf("BUG with HLT_MET\n");
     }
   }
 
   // HLT TRIGGER BASED ON 1 JET!
   if (TrIndex_Unknown != tr.triggerIndex("HLT_Jet370_v1")) {
     if (tr.accept(tr.triggerIndex("HLT_Jet370_v1")))
       hltjet->Fill(1);
     else
       hltjet->Fill(0);
   } else {
     if (TrIndex_Unknown != tr.triggerIndex("HLT_Jet100U")) {
       if (IncreasedTreshold(trEv, InputTag("hlt1jet100U", "", "HLT"), 140, 5., 1, false))
         hltjet->Fill(1);
       else
         hltjet->Fill(0);
     } else {
       if (TrIndex_Unknown != tr.triggerIndex("HLT_Jet70U")) {
         if (IncreasedTreshold(trEv, InputTag("hlt1jet70U", "", "HLT"), 140, 5., 1, false))
           hltjet->Fill(1);
         else
           hltjet->Fill(0);
       } else {
         if (TrIndex_Unknown != tr.triggerIndex("HLT_Jet50U")) {
           if (IncreasedTreshold(trEv, InputTag("hlt1jet50U", "", "HLT"), 140, 2.5, 1, false))
             hltjet->Fill(1);
           else
             hltjet->Fill(0);
         } else {
           printf("BUG with HLT_Jet\n");
         }
       }
     }
   }
 }

void HSCPValidator::makeRecoPlots ( const edm::Event & iEvent )

private

Definition at line 676 of file HSCPValidator.cc.

References dEdxTrackToken_, dedxVsp, HLT_2018_cff::distance, spr::find(), edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), mps_fire::i, HLT_2018_cff::minR, AlCaHLTBitMon_ParallelJobs::p, particleIds_, particleStatus_, funct::pow(), edm::Handle< T >::product(), RecoHSCPPtVsGenPt, mathSSE::sqrt(), tkTracksToken_, and token_.

Referenced by analyze().

                                                         {
   using namespace edm;
   using namespace reco;
 
   Handle<HepMCProduct> evt;
   iEvent.getByToken(token_, evt);
 
   Handle<TrackCollection> tkTracks;
   iEvent.getByToken(tkTracksToken_, tkTracks);
   const reco::TrackCollection tkTC = *(tkTracks.product());
 
   Handle<ValueMap<DeDxData> > dEdxTrackHandle;
   iEvent.getByToken(dEdxTrackToken_, dEdxTrackHandle);
   const ValueMap<DeDxData> dEdxTrack = *dEdxTrackHandle.product();
 
   for (size_t i = 0; i < tkTracks->size(); i++) {
     reco::TrackRef trkRef = reco::TrackRef(tkTracks, i);
 
     if (trkRef->pt() < 5 || trkRef->normalizedChi2() > 10)
       continue;
 
     double minR = 999;
     double hscpgenPt = -1;
 
     HepMC::GenEvent* myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
     for (HepMC::GenEvent::particle_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end(); ++p) {
       if ((*p)->status() != particleStatus_)
         continue;
       // Check if the particleId is in our R-hadron list
       std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(), particleIds_.end(), (*p)->pdg_id());
       if (partIdItr != particleIds_.end()) {
         //calculate DeltaR
         double distance =
             pow((*p)->momentum().eta() - trkRef->eta(), 2) + pow((*p)->momentum().phi() - trkRef->phi(), 2);
         distance = sqrt(distance);
         if (distance < minR) {
           minR = distance;
           hscpgenPt = (*p)->momentum().perp();
         }
       }
     }
     RecoHSCPPtVsGenPt->Fill(trkRef->pt(), hscpgenPt);
 
     delete myGenEvent;
     double dedx = dEdxTrack[trkRef].dEdx();
     dedxVsp->Fill(trkRef->p(), dedx);
   }
 }

void HSCPValidator::makeSimDigiPlotsECAL ( const edm::Event & iEvent )

private

Definition at line 479 of file HSCPValidator.cc.

References EcalMGPASample::adc(), edm::DataFrameContainer::begin(), gather_cfg::cout, digiOccupancyMapEB_, digiOccupancyMapEEM_, digiOccupancyMapEEP_, digisEcalNumHistEB_, digisEcalNumHistEE_, EBDigiCollectionToken_, ebSimHitToken_, EEDigiCollectionToken_, SimL1EmulatorDM_cff::eeDigis, hgcalTriggerNtuples_cfi::eeSimHits, eeSimHitToken_, edm::DataFrameContainer::end(), spr::find(), edm::Event::getByToken(), mps_fire::i, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, edm::HandleBase::isValid(), particleIds_, edm::Handle< T >::product(), EcalDataFrame::sample(), simHitsEcalDigiMatchEnergyHistEB_, simHitsEcalDigiMatchEnergyHistEE_, simHitsEcalDigiMatchEnergyVsTimeHistEB_, simHitsEcalDigiMatchEnergyVsTimeHistEE_, simHitsEcalDigiMatchIEtaHist_, simHitsEcalDigiMatchIPhiHist_, simHitsEcalDigiMatchTimeHistEB_, simHitsEcalDigiMatchTimeHistEE_, simHitsEcalEnergyHistEB_, simHitsEcalEnergyHistEE_, simHitsEcalEnergyVsTimeHistEB_, simHitsEcalEnergyVsTimeHistEE_, simHitsEcalNumHistEB_, simHitsEcalNumHistEE_, simHitsEcalTimeHistEB_, simHitsEcalTimeHistEE_, muonSimHitMatcherPSet::simTrack, TrackCandidateProducer_cfi::simTracks, simTrackToken_, EcalDataFrame::size(), and ecaldqm::zside().

Referenced by analyze().

                                                                {
   using namespace edm;
   // EB SimHits
   Handle<PCaloHitContainer> ebSimHits;
   iEvent.getByToken(ebSimHitToken_, ebSimHits);
   if (!ebSimHits.isValid()) {
     std::cout << "Cannot get EBSimHits from event!" << std::endl;
     return;
   }
   // EE SimHits
   Handle<PCaloHitContainer> eeSimHits;
   iEvent.getByToken(eeSimHitToken_, eeSimHits);
   if (!eeSimHits.isValid()) {
     std::cout << "Cannot get EESimHits from event!" << std::endl;
     return;
   }
   // SimTracks
   Handle<SimTrackContainer> simTracks;
   iEvent.getByToken(simTrackToken_, simTracks);
   if (!simTracks.isValid()) {
     std::cout << "Cannot get SimTracks from event!" << std::endl;
     return;
   }
   // EB Digis
   Handle<EBDigiCollection> ebDigis;
   iEvent.getByToken(EBDigiCollectionToken_, ebDigis);
   if (!ebDigis.isValid()) {
     std::cout << "Cannot get EBDigis from event!" << std::endl;
     return;
   }
   // EE Digis
   Handle<EEDigiCollection> eeDigis;
   iEvent.getByToken(EEDigiCollectionToken_, eeDigis);
   if (!eeDigis.isValid()) {
     std::cout << "Cannot get EEDigis from event!" << std::endl;
     return;
   }
 
   // EB first
   // 1) Look at SimTracks, getting only the HSCP tracks
   // 2) Match to PCaloHits
   // 3) Match to digis
   int numMatchedSimHitsEventEB = 0;
   int numMatchedDigisEventEB = 0;
   const PCaloHitContainer* phitsEB = nullptr;
   phitsEB = ebSimHits.product();
   for (SimTrackContainer::const_iterator simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack) {
     // Check if the particleId is in our list
     std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(), particleIds_.end(), simTrack->type());
     if (partIdItr == particleIds_.end())
       continue;
 
     PCaloHitContainer mySimHitsEB;
     std::vector<EBDataFrame> myDigisEB;
 
     //int particleId = simTrack->type();
     int trackId = simTrack->trackId();
     PCaloHitContainer::const_iterator simHitItr = phitsEB->begin();
     while (simHitItr != phitsEB->end()) {
       if (simHitItr->geantTrackId() == trackId)
         mySimHitsEB.push_back(*simHitItr);
       ++simHitItr;
     }
     if (mySimHitsEB.empty()) {
       std::cout << "Could not find matching EB PCaloHits for SimTrack id: " << trackId << ".  Skipping this SimTrack"
                 << std::endl;
       continue;
     }
 
     // Loop over matching PCaloHits
     for (simHitItr = mySimHitsEB.begin(); simHitItr != mySimHitsEB.end(); ++simHitItr) {
       simHitsEcalEnergyHistEB_->Fill(simHitItr->energy());
       simHitsEcalTimeHistEB_->Fill(simHitItr->time());
       simHitsEcalEnergyVsTimeHistEB_->Fill(simHitItr->time(), simHitItr->energy());
       EBDetId simHitId = EBDetId(simHitItr->id());
       std::cout << "SimHit DetId found: " << simHitId << " for PDGid: " << simTrack->type() << std::endl;
       //std::cout << "SimHit hashedIndex: " << simHitId.hashedIndex() << std::endl;
       std::cout << "SimHit energy: " << simHitItr->energy() << " time: " << simHitItr->time() << std::endl;
       ++numMatchedSimHitsEventEB;
 
       EBDigiCollection::const_iterator digiItr = ebDigis->begin();
       while (digiItr != ebDigis->end() && (digiItr->id() != simHitId))
         ++digiItr;
       if (digiItr == ebDigis->end()) {
         // Commented out for debugging ease, Aug 3 2009
         std::cout << "Could not find simHit detId: " << simHitId << "in EBDigiCollection!" << std::endl;
         continue;
       }
       std::vector<EBDataFrame>::const_iterator myDigiItr = myDigisEB.begin();
       while (myDigiItr != myDigisEB.end() && (digiItr->id() != myDigiItr->id()))
         ++myDigiItr;
       if (myDigiItr != myDigisEB.end())
         continue;  // if this digi is already in the list, skip it
 
       ++numMatchedDigisEventEB;
       EBDataFrame df = *digiItr;
       myDigisEB.push_back(df);
       std::cout << "SAMPLE ADCs: "
                 << "\t";
       for (int i = 0; i < 10; ++i)
         std::cout << i << "\t";
       std::cout << std::endl << "\t\t";
       for (int i = 0; i < df.size(); ++i) {
         std::cout << df.sample(i).adc() << "\t";
       }
       std::cout << std::endl << std::endl;
 
       simHitsEcalDigiMatchEnergyHistEB_->Fill(simHitItr->energy());
       simHitsEcalDigiMatchTimeHistEB_->Fill(simHitItr->time());
       simHitsEcalDigiMatchEnergyVsTimeHistEB_->Fill(simHitItr->time(), simHitItr->energy());
       simHitsEcalDigiMatchIEtaHist_->Fill(((EBDetId)digiItr->id()).ieta());
       simHitsEcalDigiMatchIPhiHist_->Fill(((EBDetId)digiItr->id()).iphi());
       digiOccupancyMapEB_->Fill(((EBDetId)digiItr->id()).iphi(), ((EBDetId)digiItr->id()).ieta());
     }
   }
   simHitsEcalNumHistEB_->Fill(numMatchedSimHitsEventEB);
   digisEcalNumHistEB_->Fill(numMatchedDigisEventEB);
 
   // EE next
   int numMatchedSimHitsEventEE = 0;
   int numMatchedDigisEventEE = 0;
   const PCaloHitContainer* phitsEE = nullptr;
   phitsEE = eeSimHits.product();
   for (SimTrackContainer::const_iterator simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack) {
     // Check if the particleId is in our list
     std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(), particleIds_.end(), simTrack->type());
     if (partIdItr == particleIds_.end())
       continue;
 
     PCaloHitContainer mySimHitsEE;
     std::vector<EEDataFrame> myDigisEE;
 
     //int particleId = simTrack->type();
     int trackId = simTrack->trackId();
     PCaloHitContainer::const_iterator simHitItr = phitsEE->begin();
     while (simHitItr != phitsEE->end()) {
       if (simHitItr->geantTrackId() == trackId)
         mySimHitsEE.push_back(*simHitItr);
       ++simHitItr;
     }
     if (mySimHitsEE.empty()) {
       std::cout << "Could not find matching EE PCaloHits for SimTrack id: " << trackId << ".  Skipping this SimTrack"
                 << std::endl;
       continue;
     }
 
     // Loop over matching PCaloHits
     for (simHitItr = mySimHitsEE.begin(); simHitItr != mySimHitsEE.end(); ++simHitItr) {
       simHitsEcalEnergyHistEE_->Fill(simHitItr->energy());
       simHitsEcalTimeHistEE_->Fill(simHitItr->time());
       simHitsEcalEnergyVsTimeHistEE_->Fill(simHitItr->time(), simHitItr->energy());
       EEDetId simHitId = EEDetId(simHitItr->id());
       std::cout << "SimHit DetId found: " << simHitId << " for PDGid: " << simTrack->type() << std::endl;
       //std::cout << "SimHit hashedIndex: " << simHitId.hashedIndex() << std::endl;
       std::cout << "SimHit energy: " << simHitItr->energy() << " time: " << simHitItr->time() << std::endl;
       ++numMatchedSimHitsEventEE;
 
       EEDigiCollection::const_iterator digiItr = eeDigis->begin();
       while (digiItr != eeDigis->end() && (digiItr->id() != simHitId))
         ++digiItr;
       if (digiItr == eeDigis->end()) {
         // Commented out for debugging ease, Aug 3 2009
         std::cout << "Could not find simHit detId: " << simHitId << "in EEDigiCollection!" << std::endl;
         continue;
       }
       std::vector<EEDataFrame>::const_iterator myDigiItr = myDigisEE.begin();
       while (myDigiItr != myDigisEE.end() && (digiItr->id() != myDigiItr->id()))
         ++myDigiItr;
       if (myDigiItr != myDigisEE.end())
         continue;  // if this digi is already in the list, skip it
 
       ++numMatchedDigisEventEE;
       EEDataFrame df = *digiItr;
       myDigisEE.push_back(df);
       std::cout << "SAMPLE ADCs: "
                 << "\t";
       for (int i = 0; i < 10; ++i)
         std::cout << i << "\t";
       std::cout << std::endl << "\t\t";
       for (int i = 0; i < df.size(); ++i) {
         std::cout << df.sample(i).adc() << "\t";
       }
       std::cout << std::endl << std::endl;
 
       simHitsEcalDigiMatchEnergyHistEE_->Fill(simHitItr->energy());
       simHitsEcalDigiMatchTimeHistEE_->Fill(simHitItr->time());
       simHitsEcalDigiMatchEnergyVsTimeHistEE_->Fill(simHitItr->time(), simHitItr->energy());
       if (((EEDetId)digiItr->id()).zside() > 0)
         digiOccupancyMapEEP_->Fill(((EEDetId)digiItr->id()).ix(), ((EEDetId)digiItr->id()).iy());
       else if (((EEDetId)digiItr->id()).zside() < 0)
         digiOccupancyMapEEM_->Fill(((EEDetId)digiItr->id()).ix(), ((EEDetId)digiItr->id()).iy());
     }
   }
   simHitsEcalNumHistEE_->Fill(numMatchedSimHitsEventEE);
   digisEcalNumHistEE_->Fill(numMatchedDigisEventEE);
 }

void HSCPValidator::makeSimDigiPlotsRPC ( const edm::Event & iEvent )

private

Definition at line 726 of file HSCPValidator.cc.

References l1GtPatternGenerator_cfi::bx, cluSizeDistribution_, gather_cfg::cout, efficiencyRPCRecHitSimDigis_, PV3DBase< T, PVType, FrameType >::eta(), HcalObjRepresent::Fill(), spr::find(), edm::Event::getByToken(), edm::Event::getManyByType(), mps_fire::i, createfilelist::int, RPCDetId::layer(), PV3DBase< T, PVType, FrameType >::mag(), DetId::Muon, RPCGeomServ::name(), particleIds_, PV3DBase< T, PVType, FrameType >::phi(), rpcPointValidation_cfi::recHit, RPCDetId::region(), residualsRPCRecHitSimDigis_, RPCGeometry::roll(), MuonSubdetId::RPC, rpcBXBarrel_, rpcBXEndCap_, rpcGeo, JetIDParams_cfi::rpcRecHits, rpcRecHitsToken_, rpcTimeOfFlightBarrel_, rpcTimeOfFlightEndCap_, RPCDetId::station(), GeomDet::surface(), and PV3DBase< T, PVType, FrameType >::x().

Referenced by analyze().

                                                               {
   using namespace edm;
 
   //std::cout << " Getting the SimHits " <<std::endl;
   std::vector<Handle<edm::PSimHitContainer> > theSimHitContainers;
   iEvent.getManyByType(theSimHitContainers);
   //std::cout << " The Number of sim Hits is  " << theSimHitContainers.size() <<std::endl;
 
   Handle<RPCRecHitCollection> rpcRecHits;
   iEvent.getByToken(rpcRecHitsToken_, rpcRecHits);
 
   //SimTrack Stuff
   std::vector<PSimHit> theSimHits;
 
   for (int i = 0; i < int(theSimHitContainers.size()); i++) {
     theSimHits.insert(theSimHits.end(), theSimHitContainers.at(i)->begin(), theSimHitContainers.at(i)->end());
   }
 
   for (std::vector<PSimHit>::const_iterator iHit = theSimHits.begin(); iHit != theSimHits.end(); iHit++) {
     std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(), particleIds_.end(), (*iHit).particleType());
     if (partIdItr == particleIds_.end())
       continue;
 
     DetId theDetUnitId((*iHit).detUnitId());
 
     DetId simdetid = DetId((*iHit).detUnitId());
 
     if (simdetid.det() == DetId::Muon && simdetid.subdetId() == MuonSubdetId::RPC) {  //Only RPCs
 
       RPCDetId rollId(theDetUnitId);
       RPCGeomServ rpcsrv(rollId);
 
       //std::cout << " Reading the Roll"<<std::endl;
       const RPCRoll* rollasociated = rpcGeo->roll(rollId);
 
       //std::cout << " Getting the Surface"<<std::endl;
       const BoundPlane& RPCSurface = rollasociated->surface();
 
       GlobalPoint SimHitInGlobal = RPCSurface.toGlobal((*iHit).localPosition());
 
       std::cout << "\t\t We have an RPC Sim Hit! in t=" << (*iHit).timeOfFlight() << "ns " << rpcsrv.name()
                 << " Global postition=" << SimHitInGlobal << std::endl;
 
       int layer = 0;
 
       if (rollId.station() == 1 && rollId.layer() == 1)
         layer = 1;
       else if (rollId.station() == 1 && rollId.layer() == 2)
         layer = 2;
       else if (rollId.station() == 2 && rollId.layer() == 1)
         layer = 3;
       else if (rollId.station() == 2 && rollId.layer() == 2)
         layer = 4;
       else if (rollId.station() == 3)
         layer = 5;
       else if (rollId.station() == 4)
         layer = 6;
 
       if (rollId.region() == 0) {
         rpcTimeOfFlightBarrel_[layer - 1]->Fill((*iHit).timeOfFlight());
       } else {
         rpcTimeOfFlightEndCap_[rollId.station() - 1]->Fill((*iHit).timeOfFlight());
       }
 
       std::cout << "\t\t r=" << SimHitInGlobal.mag() << " phi=" << SimHitInGlobal.phi()
                 << " eta=" << SimHitInGlobal.eta() << std::endl;
 
       int cluSize = 0;
       int bx = 100;
       float minres = 3000.;
 
       std::cout << "\t \t \t \t Getting RecHits in Roll Asociated" << std::endl;
 
       typedef std::pair<RPCRecHitCollection::const_iterator, RPCRecHitCollection::const_iterator> rangeRecHits;
       rangeRecHits recHitCollection = rpcRecHits->get(rollId);
       RPCRecHitCollection::const_iterator recHit;
 
       efficiencyRPCRecHitSimDigis_->Fill(0);
 
       for (recHit = recHitCollection.first; recHit != recHitCollection.second; recHit++) {
         LocalPoint recHitPos = recHit->localPosition();
         float res = (*iHit).localPosition().x() - recHitPos.x();
         if (fabs(res) < fabs(minres)) {
           minres = res;
           cluSize = recHit->clusterSize();
           bx = recHit->BunchX();
           std::cout << "\t New Min Res " << res << "cm." << std::endl;
         }
       }
 
       if (minres < 3000.) {
         residualsRPCRecHitSimDigis_->Fill(minres);
         efficiencyRPCRecHitSimDigis_->Fill(1);
         cluSizeDistribution_->Fill(cluSize);
         if (rollId.region() == 0)
           rpcBXBarrel_[layer - 1]->Fill(bx);
         else
           rpcBXEndCap_[rollId.station() - 1]->Fill(bx);
       }
     }
   }
 }

void HSCPValidator::makeSimTrackPlots ( const edm::Event & iEvent )

private

Definition at line 357 of file HSCPValidator.cc.

References spr::find(), edm::Event::getByToken(), particleIds_, edm::Handle< T >::product(), muonSimHitMatcherPSet::simTrack, simTrackParticleBetaHist_, simTrackParticleEtaHist_, simTrackParticlePhiHist_, simTrackParticlePHist_, simTrackParticlePtHist_, TrackCandidateProducer_cfi::simTracks, and simTracksToken_.

Referenced by analyze().

                                                             {
   using namespace edm;
   //get sim track infos
   Handle<edm::SimTrackContainer> simTracksHandle;
   iEvent.getByToken(simTracksToken_, simTracksHandle);
   const SimTrackContainer simTracks = *(simTracksHandle.product());
 
   SimTrackContainer::const_iterator simTrack;
 
   for (simTrack = simTracks.begin(); simTrack != simTracks.end(); ++simTrack) {
     // Check if the particleId is in our list
     std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(), particleIds_.end(), simTrack->type());
     if (partIdItr == particleIds_.end())
       continue;
 
     simTrackParticleEtaHist_->Fill((*simTrack).momentum().eta());
     simTrackParticlePhiHist_->Fill((*simTrack).momentum().phi());
     simTrackParticlePHist_->Fill((*simTrack).momentum().P());
 
     simTrackParticlePtHist_->Fill((*simTrack).momentum().pt());
 
     simTrackParticleBetaHist_->Fill((*simTrack).momentum().P() / (*simTrack).momentum().e());
 
     // std::cout<<"Particle:"<<simTrack->type()<<" Charge:"<<simTrack->charge()<<std::endl;
   }
 }

Member Data Documentation

TH1F* HSCPValidator::cluSizeDistribution_

private

Definition at line 138 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().

edm::EDGetTokenT<edm::ValueMap<reco::DeDxData> > HSCPValidator::dEdxTrackToken_

private

Definition at line 80 of file HSCPValidator.h.

Referenced by makeRecoPlots().

TH2F* HSCPValidator::dedxVsp

private

Definition at line 149 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeRecoPlots().

TH2F* HSCPValidator::digiOccupancyMapEB_

private

Definition at line 132 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH2F* HSCPValidator::digiOccupancyMapEEM_

private

Definition at line 134 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH2F* HSCPValidator::digiOccupancyMapEEP_

private

Definition at line 133 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::digisEcalNumHistEB_

private

Definition at line 130 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::digisEcalNumHistEE_

private

Definition at line 131 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

bool HSCPValidator::doGenPlots_

private

Definition at line 69 of file HSCPValidator.h.

Referenced by analyze().

bool HSCPValidator::doHLTPlots_

private

Definition at line 70 of file HSCPValidator.h.

Referenced by analyze().

bool HSCPValidator::doRecoPlots_

private

Definition at line 73 of file HSCPValidator.h.

Referenced by analyze().

bool HSCPValidator::doSimDigiPlots_

private

Definition at line 72 of file HSCPValidator.h.

Referenced by analyze().

bool HSCPValidator::doSimTrackPlots_

private

Definition at line 71 of file HSCPValidator.h.

Referenced by analyze().

edm::EDGetTokenT<EBDigiCollection> HSCPValidator::EBDigiCollectionToken_

private

Definition at line 110 of file HSCPValidator.h.

Referenced by makeSimDigiPlotsECAL().

edm::EDGetTokenT<edm::PCaloHitContainer> HSCPValidator::ebSimHitToken_

private

Definition at line 107 of file HSCPValidator.h.

Referenced by makeSimDigiPlotsECAL().

edm::EDGetTokenT<EEDigiCollection> HSCPValidator::EEDigiCollectionToken_

private

Definition at line 111 of file HSCPValidator.h.

Referenced by makeSimDigiPlotsECAL().

edm::EDGetTokenT<edm::PCaloHitContainer> HSCPValidator::eeSimHitToken_

private

Definition at line 108 of file HSCPValidator.h.

Referenced by makeSimDigiPlotsECAL().

TH1F* HSCPValidator::efficiencyRPCRecHitSimDigis_

private

Definition at line 137 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().

TH1F* HSCPValidator::h_genhscp_met

private

Definition at line 94 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

TH1F* HSCPValidator::h_genhscp_met_nohscp

private

Definition at line 95 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

TH1F* HSCPValidator::h_genhscp_scaloret

private

Definition at line 96 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

TH1F* HSCPValidator::h_genhscp_scaloret_nohscp

private

Definition at line 97 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

TH1F* HSCPValidator::hltjet

private

Definition at line 145 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeHLTPlots().

TH1F* HSCPValidator::hltmet

private

Definition at line 144 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeHLTPlots().

TH1F* HSCPValidator::hltmu

private

Definition at line 146 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeHLTPlots().

TH1F* HSCPValidator::particleBetaHist_

private

Definition at line 92 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

TH1F* HSCPValidator::particleBetaInverseHist_

private

Definition at line 93 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

TH1F* HSCPValidator::particleEtaHist_

private

Definition at line 86 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

std::vector<int> HSCPValidator::particleIds_

private

Definition at line 83 of file HSCPValidator.h.

Referenced by makeGenPlots(), makeRecoPlots(), makeSimDigiPlotsECAL(), makeSimDigiPlotsRPC(), and makeSimTrackPlots().

std::map<int, int> HSCPValidator::particleIdsFoundMap_

private

Definition at line 85 of file HSCPValidator.h.

Referenced by endJob(), and makeGenPlots().

TH1F* HSCPValidator::particleMassHist_

private

Definition at line 90 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

TH1F* HSCPValidator::particlePhiHist_

private

Definition at line 87 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

TH1F* HSCPValidator::particlePHist_

private

Definition at line 88 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

TH1F* HSCPValidator::particlePtHist_

private

Definition at line 89 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

int HSCPValidator::particleStatus_

private

Definition at line 84 of file HSCPValidator.h.

Referenced by makeGenPlots(), and makeRecoPlots().

TH1F* HSCPValidator::particleStatusHist_

private

Definition at line 91 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeGenPlots().

TH2F* HSCPValidator::RecoHSCPPtVsGenPt

private

Definition at line 148 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeRecoPlots().

TH1F* HSCPValidator::residualsRPCRecHitSimDigis_

private

Definition at line 136 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().

TH1F* HSCPValidator::rpcBXBarrel_[6]

private

Definition at line 140 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().

TH1F* HSCPValidator::rpcBXEndCap_[3]

private

Definition at line 142 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().

edm::ESHandle<RPCGeometry> HSCPValidator::rpcGeo

private

Definition at line 112 of file HSCPValidator.h.

Referenced by analyze(), and makeSimDigiPlotsRPC().

edm::EDGetTokenT<RPCRecHitCollection> HSCPValidator::rpcRecHitsToken_

private

Definition at line 81 of file HSCPValidator.h.

Referenced by makeSimDigiPlotsRPC().

TH1F* HSCPValidator::rpcTimeOfFlightBarrel_[6]

private

Definition at line 139 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().

TH1F* HSCPValidator::rpcTimeOfFlightEndCap_[3]

private

Definition at line 141 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().

TH1F* HSCPValidator::simHitsEcalDigiMatchEnergyHistEB_

private

Definition at line 122 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalDigiMatchEnergyHistEE_

private

Definition at line 125 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH2F* HSCPValidator::simHitsEcalDigiMatchEnergyVsTimeHistEB_

private

Definition at line 124 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH2F* HSCPValidator::simHitsEcalDigiMatchEnergyVsTimeHistEE_

private

Definition at line 127 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalDigiMatchIEtaHist_

private

Definition at line 128 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalDigiMatchIPhiHist_

private

Definition at line 129 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalDigiMatchTimeHistEB_

private

Definition at line 123 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalDigiMatchTimeHistEE_

private

Definition at line 126 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalEnergyHistEB_

private

Definition at line 114 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalEnergyHistEE_

private

Definition at line 118 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH2F* HSCPValidator::simHitsEcalEnergyVsTimeHistEB_

private

Definition at line 117 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH2F* HSCPValidator::simHitsEcalEnergyVsTimeHistEE_

private

Definition at line 121 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalNumHistEB_

private

Definition at line 116 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalNumHistEE_

private

Definition at line 120 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalTimeHistEB_

private

Definition at line 115 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simHitsEcalTimeHistEE_

private

Definition at line 119 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().

TH1F* HSCPValidator::simTrackParticleBetaHist_

private

Definition at line 104 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimTrackPlots().

TH1F* HSCPValidator::simTrackParticleEtaHist_

private

Definition at line 100 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimTrackPlots().

TH1F* HSCPValidator::simTrackParticlePhiHist_

private

Definition at line 101 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimTrackPlots().

TH1F* HSCPValidator::simTrackParticlePHist_

private

Definition at line 102 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimTrackPlots().

TH1F* HSCPValidator::simTrackParticlePtHist_

private

Definition at line 103 of file HSCPValidator.h.

Referenced by HSCPValidator(), and makeSimTrackPlots().

edm::EDGetTokenT<edm::SimTrackContainer> HSCPValidator::simTracksToken_

private

Definition at line 77 of file HSCPValidator.h.

Referenced by makeSimTrackPlots().

edm::EDGetTokenT<edm::SimTrackContainer> HSCPValidator::simTrackToken_

private

Definition at line 109 of file HSCPValidator.h.

Referenced by makeSimDigiPlotsECAL().

edm::EDGetTokenT<reco::TrackCollection> HSCPValidator::tkTracksToken_

private

Definition at line 79 of file HSCPValidator.h.

Referenced by makeRecoPlots().

edm::EDGetTokenT<edm::HepMCProduct> HSCPValidator::token_

private

Definition at line 76 of file HSCPValidator.h.

Referenced by makeGenPlots(), and makeRecoPlots().

edm::EDGetTokenT<trigger::TriggerEvent> HSCPValidator::trEvToken_

private

Definition at line 78 of file HSCPValidator.h.

Referenced by makeHLTPlots().

edm::EDGetTokenT<edm::TriggerResults> HSCPValidator::triggerResultsToken_

private

Definition at line 82 of file HSCPValidator.h.

Referenced by makeHLTPlots().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation