GlobalRecHitsAnalyzer Class Reference

#include <GlobalRecHitsAnalyzer.h>

Inheritance diagram for GlobalRecHitsAnalyzer:

Public Types
typedef std::map< uint32_t, float, std::less< uint32_t > >	MapType
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
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

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	GlobalRecHitsAnalyzer (const edm::ParameterSet &)
	~GlobalRecHitsAnalyzer () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () 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
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)

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
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)

Private Member Functions
template<typename type >
int	compute (const DTGeometry *dtGeom, const std::map< DTWireId, std::vector< PSimHit >> &simHitsPerWire, const std::map< DTWireId, std::vector< type >> &recHitsPerWire, int step)
void	fillECal (const edm::Event &, const edm::EventSetup &)
void	fillHCal (const edm::Event &, const edm::EventSetup &)
void	fillMuon (const edm::Event &, const edm::EventSetup &)
void	fillTrk (const edm::Event &, const edm::EventSetup &)
template<typename type >
const type *	findBestRecHit (const DTLayer *layer, DTWireId wireId, const std::vector< type > &recHits, const float simHitDist)
std::map< DTWireId, std::vector< DTRecHit1DPair > >	map1DRecHitsPerWire (const DTRecHitCollection *dt1DRecHitPairs)
void	plotResolution (const PSimHit &simHit, const CSCRecHit2D &recHit, const CSCLayer *layer, int chamberType)
std::pair< LocalPoint, LocalVector >	projectHit (const PSimHit &hit, const StripGeomDetUnit *stripDet, const BoundPlane &plane)
float	recHitDistFromWire (const DTRecHit1DPair &hitPair, const DTLayer *layer)
float	recHitDistFromWire (const DTRecHit1D &recHit, const DTLayer *layer)
float	simHitDistFromWire (const DTLayer *layer, DTWireId wireId, const PSimHit &hit)

Private Attributes
unsigned int	count
edm::EDGetTokenT< CrossingFrame< PCaloHit > >	EBHits_Token_
edm::InputTag	ECalEBSrc_
edm::EDGetTokenT< EBRecHitCollection >	ECalEBSrc_Token_
edm::InputTag	ECalEESrc_
edm::EDGetTokenT< EERecHitCollection >	ECalEESrc_Token_
edm::InputTag	ECalESSrc_
edm::EDGetTokenT< ESRecHitCollection >	ECalESSrc_Token_
edm::InputTag	ECalUncalEBSrc_
edm::EDGetTokenT< EBUncalibratedRecHitCollection >	ECalUncalEBSrc_Token_
edm::InputTag	ECalUncalEESrc_
edm::EDGetTokenT< EEUncalibratedRecHitCollection >	ECalUncalEESrc_Token_
edm::EDGetTokenT< CrossingFrame< PCaloHit > >	EEHits_Token_
edm::EDGetTokenT< CrossingFrame< PCaloHit > >	ESHits_Token_
std::string	fName
int	frequency
bool	getAllProvenances
edm::InputTag	HCalSrc_
edm::EDGetTokenT< edm::PCaloHitContainer >	HCalSrc_Token_
std::string	hitsProducer
std::string	label
std::vector< PSimHit >	matched
MonitorElement *	mehCSCn
MonitorElement *	mehCSCResRDPhi
MonitorElement *	mehDtMuonn
MonitorElement *	mehDtMuonRes
MonitorElement *	mehEcaln [3]
MonitorElement *	mehEcalRes [3]
MonitorElement *	mehHcaln [4]
MonitorElement *	mehHcalRes [4]
MonitorElement *	mehRPCn
MonitorElement *	mehRPCResX
MonitorElement *	mehSiPixeln [7]
MonitorElement *	mehSiPixelResX [7]
MonitorElement *	mehSiPixelResY [7]
MonitorElement *	mehSiStripn [19]
MonitorElement *	mehSiStripResX [19]
MonitorElement *	mehSiStripResY [19]
edm::EDGetTokenT< CrossingFrame< PSimHit > >	MuCSCHits_Token_
edm::InputTag	MuCSCSrc_
edm::EDGetTokenT< CSCRecHit2DCollection >	MuCSCSrc_Token_
edm::InputTag	MuDTSimSrc_
edm::EDGetTokenT< edm::PSimHitContainer >	MuDTSimSrc_Token_
edm::InputTag	MuDTSrc_
edm::EDGetTokenT< DTRecHitCollection >	MuDTSrc_Token_
edm::InputTag	MuRPCSimSrc_
edm::EDGetTokenT< edm::PSimHitContainer >	MuRPCSimSrc_Token_
edm::InputTag	MuRPCSrc_
edm::EDGetTokenT< RPCRecHitCollection >	MuRPCSrc_Token_
bool	printProvenanceInfo
edm::InputTag	SiPxlSrc_
edm::EDGetTokenT< SiPixelRecHitCollection >	SiPxlSrc_Token_
edm::InputTag	SiStripSrc_
edm::EDGetTokenT< SiStripMatchedRecHit2DCollection >	SiStripSrc_Token_
std::map< int, edm::PSimHitContainer >	theMap
TrackerHitAssociator::Config	trackerHitAssociatorConfig_
int	verbosity

Additional Inherited Members
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

Definition at line 140 of file GlobalRecHitsAnalyzer.h.

Member Typedef Documentation

typedef std::map<uint32_t, float, std::less<uint32_t> > GlobalRecHitsAnalyzer::MapType

Definition at line 142 of file GlobalRecHitsAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 16 of file GlobalRecHitsAnalyzer.cc.

References EBHits_Token_, ECalEBSrc_, ECalEBSrc_Token_, ECalEESrc_, ECalEESrc_Token_, ECalESSrc_, ECalESSrc_Token_, ECalUncalEBSrc_, ECalUncalEBSrc_Token_, ECalUncalEESrc_, ECalUncalEESrc_Token_, EEHits_Token_, ESHits_Token_, fName, frequency, getAllProvenances, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HCalSrc_, HCalSrc_Token_, hitsProducer, HLT_2018_cff::InputTag, edm::InputTag::instance(), edm::InputTag::label(), MuCSCHits_Token_, MuCSCSrc_, MuCSCSrc_Token_, MuDTSimSrc_, MuDTSimSrc_Token_, MuDTSrc_, MuDTSrc_Token_, MuRPCSimSrc_, MuRPCSimSrc_Token_, MuRPCSrc_, MuRPCSrc_Token_, printProvenanceInfo, SiPxlSrc_, SiPxlSrc_Token_, SiStripSrc_, SiStripSrc_Token_, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

    : fName(""),
      verbosity(0),
      frequency(0),
      label(""),
      getAllProvenances(false),
      printProvenanceInfo(false),
      trackerHitAssociatorConfig_(iPSet, consumesCollector()),
      count(0) {
  consumesMany<edm::SortedCollection<HBHERecHit, edm::StrictWeakOrdering<HBHERecHit>>>();
  consumesMany<edm::SortedCollection<HFRecHit, edm::StrictWeakOrdering<HFRecHit>>>();
  consumesMany<edm::SortedCollection<HORecHit, edm::StrictWeakOrdering<HORecHit>>>();
  std::string MsgLoggerCat = "GlobalRecHitsAnalyzer_GlobalRecHitsAnalyzer";

  // get information from parameter set
  fName = iPSet.getUntrackedParameter<std::string>("Name");
  verbosity = iPSet.getUntrackedParameter<int>("Verbosity");
  frequency = iPSet.getUntrackedParameter<int>("Frequency");
  edm::ParameterSet m_Prov = iPSet.getParameter<edm::ParameterSet>("ProvenanceLookup");
  getAllProvenances = m_Prov.getUntrackedParameter<bool>("GetAllProvenances");
  printProvenanceInfo = m_Prov.getUntrackedParameter<bool>("PrintProvenanceInfo");
  hitsProducer = iPSet.getParameter<std::string>("hitsProducer");

  //get Labels to use to extract information
  ECalEBSrc_ = iPSet.getParameter<edm::InputTag>("ECalEBSrc");
  ECalUncalEBSrc_ = iPSet.getParameter<edm::InputTag>("ECalUncalEBSrc");
  ECalEESrc_ = iPSet.getParameter<edm::InputTag>("ECalEESrc");
  ECalUncalEESrc_ = iPSet.getParameter<edm::InputTag>("ECalUncalEESrc");
  ECalESSrc_ = iPSet.getParameter<edm::InputTag>("ECalESSrc");
  HCalSrc_ = iPSet.getParameter<edm::InputTag>("HCalSrc");
  SiStripSrc_ = iPSet.getParameter<edm::InputTag>("SiStripSrc");
  SiPxlSrc_ = iPSet.getParameter<edm::InputTag>("SiPxlSrc");
  MuDTSrc_ = iPSet.getParameter<edm::InputTag>("MuDTSrc");
  MuDTSimSrc_ = iPSet.getParameter<edm::InputTag>("MuDTSimSrc");
  MuCSCSrc_ = iPSet.getParameter<edm::InputTag>("MuCSCSrc");
  MuRPCSrc_ = iPSet.getParameter<edm::InputTag>("MuRPCSrc");
  MuRPCSimSrc_ = iPSet.getParameter<edm::InputTag>("MuRPCSimSrc");

  // fix for consumes
  ECalUncalEBSrc_Token_ = consumes<EBUncalibratedRecHitCollection>(iPSet.getParameter<edm::InputTag>("ECalUncalEBSrc"));
  ECalUncalEESrc_Token_ = consumes<EEUncalibratedRecHitCollection>(iPSet.getParameter<edm::InputTag>("ECalUncalEESrc"));
  ECalEBSrc_Token_ = consumes<EBRecHitCollection>(iPSet.getParameter<edm::InputTag>("ECalEBSrc"));
  ECalEESrc_Token_ = consumes<EERecHitCollection>(iPSet.getParameter<edm::InputTag>("ECalEESrc"));
  ECalESSrc_Token_ = consumes<ESRecHitCollection>(iPSet.getParameter<edm::InputTag>("ECalESSrc"));
  HCalSrc_Token_ = consumes<edm::PCaloHitContainer>(iPSet.getParameter<edm::InputTag>("HCalSrc"));
  SiStripSrc_Token_ = consumes<SiStripMatchedRecHit2DCollection>(iPSet.getParameter<edm::InputTag>("SiStripSrc"));
  SiPxlSrc_Token_ = consumes<SiPixelRecHitCollection>(iPSet.getParameter<edm::InputTag>("SiPxlSrc"));

  MuDTSrc_Token_ = consumes<DTRecHitCollection>(iPSet.getParameter<edm::InputTag>("MuDTSrc"));
  MuDTSimSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("MuDTSimSrc"));

  MuCSCSrc_Token_ = consumes<CSCRecHit2DCollection>(iPSet.getParameter<edm::InputTag>("MuCSCSrc"));
  MuCSCHits_Token_ = consumes<CrossingFrame<PSimHit>>(
      edm::InputTag(std::string("mix"), iPSet.getParameter<std::string>("hitsProducer") + std::string("MuonCSCHits")));

  MuRPCSrc_Token_ = consumes<RPCRecHitCollection>(iPSet.getParameter<edm::InputTag>("MuRPCSrc"));
  MuRPCSimSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("MuRPCSimSrc"));

  EBHits_Token_ = consumes<CrossingFrame<PCaloHit>>(
      edm::InputTag(std::string("mix"), iPSet.getParameter<std::string>("hitsProducer") + std::string("EcalHitsEB")));
  EEHits_Token_ = consumes<CrossingFrame<PCaloHit>>(
      edm::InputTag(std::string("mix"), iPSet.getParameter<std::string>("hitsProducer") + std::string("EcalHitsEE")));
  ESHits_Token_ = consumes<CrossingFrame<PCaloHit>>(
      edm::InputTag(std::string("mix"), iPSet.getParameter<std::string>("hitsProducer") + std::string("EcalHitsES")));

  // 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;

  // create persistent object
  // produces<PGlobalRecHit>(label);

  // print out Parameter Set information being used
  if (verbosity >= 0) {
    edm::LogInfo(MsgLoggerCat) << "\n===============================\n"
                               << "Initialized as EDProducer with parameter values:\n"
                               << "    Name           = " << fName << "\n"
                               << "    Verbosity      = " << verbosity << "\n"
                               << "    Frequency      = " << frequency << "\n"
                               << "    GetProv        = " << getAllProvenances << "\n"
                               << "    PrintProv      = " << printProvenanceInfo << "\n"
                               << "    ECalEBSrc      = " << ECalEBSrc_.label() << ":" << ECalEBSrc_.instance() << "\n"
                               << "    ECalUncalEBSrc = " << ECalUncalEBSrc_.label() << ":"
                               << ECalUncalEBSrc_.instance() << "\n"
                               << "    ECalEESrc      = " << ECalEESrc_.label() << ":" << ECalUncalEESrc_.instance()
                               << "\n"
                               << "    ECalUncalEESrc = " << ECalUncalEESrc_.label() << ":" << ECalEESrc_.instance()
                               << "\n"
                               << "    ECalESSrc      = " << ECalESSrc_.label() << ":" << ECalESSrc_.instance() << "\n"
                               << "    HCalSrc        = " << HCalSrc_.label() << ":" << HCalSrc_.instance() << "\n"
                               << "    SiStripSrc     = " << SiStripSrc_.label() << ":" << SiStripSrc_.instance()
                               << "\n"
                               << "    SiPixelSrc     = " << SiPxlSrc_.label() << ":" << SiPxlSrc_.instance() << "\n"
                               << "    MuDTSrc        = " << MuDTSrc_.label() << ":" << MuDTSrc_.instance() << "\n"
                               << "    MuDTSimSrc     = " << MuDTSimSrc_.label() << ":" << MuDTSimSrc_.instance()
                               << "\n"
                               << "    MuCSCSrc       = " << MuCSCSrc_.label() << ":" << MuCSCSrc_.instance() << "\n"
                               << "    MuRPCSrc       = " << MuRPCSrc_.label() << ":" << MuRPCSrc_.instance() << "\n"
                               << "    MuRPCSimSrc    = " << MuRPCSimSrc_.label() << ":" << MuRPCSimSrc_.instance()
                               << "\n"
                               << "===============================\n";
  }
}

GlobalRecHitsAnalyzer::~GlobalRecHitsAnalyzer ( )

override

Definition at line 120 of file GlobalRecHitsAnalyzer.cc.

{}

Member Function Documentation

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

override

Definition at line 292 of file GlobalRecHitsAnalyzer.cc.

References count, edm::EventID::event(), fillECal(), 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 = "GlobalRecHitsAnalyzer_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 Ecal information from event
  fillECal(iEvent, iSetup);
  // gather Hcal information from event
  fillHCal(iEvent, iSetup);
  // gather Track information from event
  fillTrk(iEvent, iSetup);
  // gather Muon information from event
  fillMuon(iEvent, iSetup);

  if (verbosity > 0)
    edm::LogInfo(MsgLoggerCat) << "Done gathering data from event.";

  return;
}

void GlobalRecHitsAnalyzer::bookHistograms ( DQMStore::IBooker &  iBooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 122 of file GlobalRecHitsAnalyzer.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), mps_fire::i, dqmiolumiharvest::j, mehCSCn, mehCSCResRDPhi, mehDtMuonn, mehDtMuonRes, mehEcaln, mehEcalRes, mehHcaln, mehHcalRes, mehRPCn, mehRPCResX, mehSiPixeln, mehSiPixelResX, mehSiPixelResY, mehSiStripn, mehSiStripResX, mehSiStripResY, dqm::impl::MonitorElement::setAxisTitle(), and dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder().

                                                                                                        {
  // Si Strip
  string SiStripString[19] = {"TECW1",
                              "TECW2",
                              "TECW3",
                              "TECW4",
                              "TECW5",
                              "TECW6",
                              "TECW7",
                              "TECW8",
                              "TIBL1",
                              "TIBL2",
                              "TIBL3",
                              "TIBL4",
                              "TIDW1",
                              "TIDW2",
                              "TIDW3",
                              "TOBL1",
                              "TOBL2",
                              "TOBL3",
                              "TOBL4"};
  for (int i = 0; i < 19; ++i) {
    mehSiStripn[i] = nullptr;
    mehSiStripResX[i] = nullptr;
    mehSiStripResY[i] = nullptr;
  }
  string hcharname, hchartitle;
  iBooker.setCurrentFolder("GlobalRecHitsV/SiStrips");
  for (int amend = 0; amend < 19; ++amend) {
    hcharname = "hSiStripn_" + SiStripString[amend];
    hchartitle = SiStripString[amend] + "  rechits";
    mehSiStripn[amend] = iBooker.book1D(hcharname, hchartitle, 200, 0., 200.);
    mehSiStripn[amend]->setAxisTitle("Number of hits in " + SiStripString[amend], 1);
    mehSiStripn[amend]->setAxisTitle("Count", 2);
    hcharname = "hSiStripResX_" + SiStripString[amend];
    hchartitle = SiStripString[amend] + " rechit x resolution";
    mehSiStripResX[amend] = iBooker.book1D(hcharname, hchartitle, 200, -0.02, .02);
    mehSiStripResX[amend]->setAxisTitle("X-resolution in " + SiStripString[amend], 1);
    mehSiStripResX[amend]->setAxisTitle("Count", 2);
    hcharname = "hSiStripResY_" + SiStripString[amend];
    hchartitle = SiStripString[amend] + " rechit y resolution";
    mehSiStripResY[amend] = iBooker.book1D(hcharname, hchartitle, 200, -0.02, .02);
    mehSiStripResY[amend]->setAxisTitle("Y-resolution in " + SiStripString[amend], 1);
    mehSiStripResY[amend]->setAxisTitle("Count", 2);
  }

  //HCal
  //string hcharname, hchartitle;
  string HCalString[4] = {"HB", "HE", "HF", "HO"};
  float HCalnUpper[4] = {3000., 3000., 3000., 3000.};
  float HCalnLower[4] = {0., 0., 0., 0.};
  for (int j = 0; j < 4; ++j) {
    mehHcaln[j] = nullptr;
    mehHcalRes[j] = nullptr;
  }

  iBooker.setCurrentFolder("GlobalRecHitsV/HCals");
  for (int amend = 0; amend < 4; ++amend) {
    hcharname = "hHcaln_" + HCalString[amend];
    hchartitle = HCalString[amend] + "  rechits";
    mehHcaln[amend] = iBooker.book1D(hcharname, hchartitle, 1000, HCalnLower[amend], HCalnUpper[amend]);
    mehHcaln[amend]->setAxisTitle("Number of RecHits", 1);
    mehHcaln[amend]->setAxisTitle("Count", 2);
    hcharname = "hHcalRes_" + HCalString[amend];
    hchartitle = HCalString[amend] + "  rechit resolution";
    mehHcalRes[amend] = iBooker.book1D(hcharname, hchartitle, 25, -2., 2.);
    mehHcalRes[amend]->setAxisTitle("RecHit E - SimHit E", 1);
    mehHcalRes[amend]->setAxisTitle("Count", 2);
  }

  //Ecal
  string ECalString[3] = {"EB", "EE", "ES"};
  int ECalnBins[3] = {1000, 3000, 150};
  float ECalnUpper[3] = {20000., 62000., 3000.};
  float ECalnLower[3] = {0., 0., 0.};
  int ECalResBins[3] = {200, 200, 200};
  float ECalResUpper[3] = {1., 0.3, .0002};
  float ECalResLower[3] = {-1., -0.3, -.0002};
  for (int i = 0; i < 3; ++i) {
    mehEcaln[i] = nullptr;
    mehEcalRes[i] = nullptr;
  }
  iBooker.setCurrentFolder("GlobalRecHitsV/ECals");

  for (int amend = 0; amend < 3; ++amend) {
    hcharname = "hEcaln_" + ECalString[amend];
    hchartitle = ECalString[amend] + "  rechits";
    mehEcaln[amend] = iBooker.book1D(hcharname, hchartitle, ECalnBins[amend], ECalnLower[amend], ECalnUpper[amend]);
    mehEcaln[amend]->setAxisTitle("Number of RecHits", 1);
    mehEcaln[amend]->setAxisTitle("Count", 2);
    hcharname = "hEcalRes_" + ECalString[amend];
    hchartitle = ECalString[amend] + "  rechit resolution";
    mehEcalRes[amend] =
        iBooker.book1D(hcharname, hchartitle, ECalResBins[amend], ECalResLower[amend], ECalResUpper[amend]);
    mehEcalRes[amend]->setAxisTitle("RecHit E - SimHit E", 1);
    mehEcalRes[amend]->setAxisTitle("Count", 2);
  }

  //Si Pixels
  string SiPixelString[7] = {"BRL1", "BRL2", "BRL3", "FWD1n", "FWD1p", "FWD2n", "FWD2p"};
  for (int j = 0; j < 7; ++j) {
    mehSiPixeln[j] = nullptr;
    mehSiPixelResX[j] = nullptr;
    mehSiPixelResY[j] = nullptr;
  }

  iBooker.setCurrentFolder("GlobalRecHitsV/SiPixels");
  for (int amend = 0; amend < 7; ++amend) {
    hcharname = "hSiPixeln_" + SiPixelString[amend];
    hchartitle = SiPixelString[amend] + " rechits";
    mehSiPixeln[amend] = iBooker.book1D(hcharname, hchartitle, 200, 0., 200.);
    mehSiPixeln[amend]->setAxisTitle("Number of hits in " + SiPixelString[amend], 1);
    mehSiPixeln[amend]->setAxisTitle("Count", 2);
    hcharname = "hSiPixelResX_" + SiPixelString[amend];
    hchartitle = SiPixelString[amend] + " rechit x resolution";
    mehSiPixelResX[amend] = iBooker.book1D(hcharname, hchartitle, 200, -0.02, .02);
    mehSiPixelResX[amend]->setAxisTitle("X-resolution in " + SiPixelString[amend], 1);
    mehSiPixelResX[amend]->setAxisTitle("Count", 2);
    hcharname = "hSiPixelResY_" + SiPixelString[amend];
    hchartitle = SiPixelString[amend] + " rechit y resolution";

    mehSiPixelResY[amend] = iBooker.book1D(hcharname, hchartitle, 200, -0.02, .02);
    mehSiPixelResY[amend]->setAxisTitle("Y-resolution in " + SiPixelString[amend], 1);
    mehSiPixelResY[amend]->setAxisTitle("Count", 2);
  }

  //Muons
  iBooker.setCurrentFolder("GlobalRecHitsV/Muons");

  mehDtMuonn = nullptr;
  mehCSCn = nullptr;
  mehRPCn = nullptr;

  string n_List[3] = {"hDtMuonn", "hCSCn", "hRPCn"};
  string hist_string[3] = {"Dt", "CSC", "RPC"};

  for (int amend = 0; amend < 3; ++amend) {
    hchartitle = hist_string[amend] + " rechits";
    if (amend == 0) {
      mehDtMuonn = iBooker.book1D(n_List[amend], hchartitle, 50, 0., 500.);
      mehDtMuonn->setAxisTitle("Number of Rechits", 1);
      mehDtMuonn->setAxisTitle("Count", 2);
    }
    if (amend == 1) {
      mehCSCn = iBooker.book1D(n_List[amend], hchartitle, 50, 0., 500.);
      mehCSCn->setAxisTitle("Number of Rechits", 1);
      mehCSCn->setAxisTitle("Count", 2);
    }
    if (amend == 2) {
      mehRPCn = iBooker.book1D(n_List[amend], hchartitle, 50, 0., 500.);
      mehRPCn->setAxisTitle("Number of Rechits", 1);
      mehRPCn->setAxisTitle("Count", 2);
    }
  }

  mehDtMuonRes = nullptr;
  mehCSCResRDPhi = nullptr;
  mehRPCResX = nullptr;

  hcharname = "hDtMuonRes";
  hchartitle = "DT wire distance resolution";
  mehDtMuonRes = iBooker.book1D(hcharname, hchartitle, 200, -0.2, 0.2);
  hcharname = "CSCResRDPhi";
  hchartitle = "CSC perp*dphi resolution";
  mehCSCResRDPhi = iBooker.book1D(hcharname, hchartitle, 200, -0.2, 0.2);
  hcharname = "hRPCResX";
  hchartitle = "RPC rechits x resolution";
  mehRPCResX = iBooker.book1D(hcharname, hchartitle, 50, -5., 5.);
}

template<typename type >

int GlobalRecHitsAnalyzer::compute ( const DTGeometry *  dtGeom, const std::map< DTWireId, std::vector< PSimHit >> &  simHitsPerWire, const std::map< DTWireId, std::vector< type >> &  recHitsPerWire, int  step  )

private

Definition at line 1479 of file GlobalRecHitsAnalyzer.cc.

References dqm::impl::MonitorElement::Fill(), findBestRecHit(), DTHitQualityUtils::findMuSimHit(), DTGeometry::layer(), genParticles_cff::map, mehDtMuonRes, recHitDistFromWire(), FastTrackerRecHitMaskProducer_cfi::recHits, and simHitDistFromWire().

Referenced by fillMuon().

                                             {
  std::map<DTWireId, std::vector<PSimHit>> simHitsPerWire = _simHitsPerWire;
  std::map<DTWireId, std::vector<type>> recHitsPerWire = _recHitsPerWire;
  int nDt = 0;
  // Loop over cells with a muon SimHit
  for (std::map<DTWireId, std::vector<PSimHit>>::const_iterator wireAndSHits = simHitsPerWire.begin();
       wireAndSHits != simHitsPerWire.end();
       wireAndSHits++) {
    DTWireId wireId = (*wireAndSHits).first;
    std::vector<PSimHit> simHitsInCell = (*wireAndSHits).second;

    // Get the layer
    const DTLayer* layer = dtGeom->layer(wireId);

    // Look for a mu hit in the cell
    const PSimHit* muSimHit = DTHitQualityUtils::findMuSimHit(simHitsInCell);
    if (muSimHit == nullptr) {
      continue;  // Skip this cell
    }

    // Find the distance of the simhit from the wire
    float simHitWireDist = simHitDistFromWire(layer, wireId, *muSimHit);
    // Skip simhits out of the cell
    if (simHitWireDist > 2.1) {
      continue;  // Skip this cell
    }

    // Look for RecHits in the same cell
    if (recHitsPerWire.find(wireId) == recHitsPerWire.end()) {
      continue;  // No RecHit found in this cell
    } else {
      std::vector<type> recHits = recHitsPerWire[wireId];

      // Find the best RecHit
      const type* theBestRecHit = findBestRecHit(layer, wireId, recHits, simHitWireDist);

      float recHitWireDist = recHitDistFromWire(*theBestRecHit, layer);

      ++nDt;

      mehDtMuonRes->Fill(recHitWireDist - simHitWireDist);

    }  // find rechits
  }    // loop over simhits

  return nDt;
}

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

private

Definition at line 358 of file GlobalRecHitsAnalyzer.cc.

References HLT_2018_cff::barrelHits, edm::SortedCollection< T, SORT >::begin(), EBHits_Token_, CollectionTags_cfi::EBRecHit, CollectionTags_cfi::EBUncalibRecHit, ECalEBSrc_Token_, ECalEESrc_Token_, ECalESSrc_Token_, ECalUncalEBSrc_Token_, ECalUncalEESrc_Token_, EEHits_Token_, CollectionTags_cfi::EERecHit, CollectionTags_cfi::EEUncalibRecHit, edm::SortedCollection< T, SORT >::end(), HLT_2018_cff::endcapHits, ESHits_Token_, dqm::impl::MonitorElement::Fill(), edm::SortedCollection< T, SORT >::find(), edm::Event::getByToken(), edm::HandleBase::isValid(), LogDebug, mehEcaln, mehEcalRes, gedPhotons_cfi::preshowerHits, edm::Handle< T >::product(), DetId::rawId(), rpcPointValidation_cfi::recHit, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

Referenced by analyze().

                                                                                        {
  std::string MsgLoggerCat = "GlobalRecHitsAnalyzer_fillECal";

  TString eventout;
  if (verbosity > 0)
    eventout = "\nGathering info:";

  // extract crossing frame from event
  edm::Handle<CrossingFrame<PCaloHit>> crossingFrame;

  //extract EB information
  edm::Handle<EBUncalibratedRecHitCollection> EcalUncalibRecHitEB;
  iEvent.getByToken(ECalUncalEBSrc_Token_, EcalUncalibRecHitEB);
  bool validUncalibRecHitEB = true;
  if (!EcalUncalibRecHitEB.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalUncalRecHitEB in event!";
    validUncalibRecHitEB = false;
  }

  edm::Handle<EBRecHitCollection> EcalRecHitEB;
  iEvent.getByToken(ECalEBSrc_Token_, EcalRecHitEB);
  bool validRecHitEB = true;
  if (!EcalRecHitEB.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalRecHitEB in event!";
    validRecHitEB = false;
  }

  // loop over simhits
  iEvent.getByToken(EBHits_Token_, crossingFrame);
  bool validXFrame = true;
  if (!crossingFrame.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find cal barrel crossingFrame in event!";
    validXFrame = false;
  }

  MapType ebSimMap;
  if (validXFrame) {
    const MixCollection<PCaloHit> barrelHits(crossingFrame.product());
    // keep track of sum of simhit energy in each crystal
    for (auto const& iHit : barrelHits) {
      EBDetId ebid = EBDetId(iHit.id());

      uint32_t crystid = ebid.rawId();
      ebSimMap[crystid] += iHit.energy();
    }
  }

  int nEBRecHits = 0;
  // loop over RecHits
  if (validUncalibRecHitEB && validRecHitEB) {
    const EBUncalibratedRecHitCollection* EBUncalibRecHit = EcalUncalibRecHitEB.product();
    const EBRecHitCollection* EBRecHit = EcalRecHitEB.product();

    for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit = EBUncalibRecHit->begin();
         uncalibRecHit != EBUncalibRecHit->end();
         ++uncalibRecHit) {
      EBDetId EBid = EBDetId(uncalibRecHit->id());

      EcalRecHitCollection::const_iterator myRecHit = EBRecHit->find(EBid);

      if (myRecHit != EBRecHit->end()) {
        ++nEBRecHits;
        mehEcalRes[1]->Fill(myRecHit->energy() - ebSimMap[EBid.rawId()]);
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of EBRecHits collected:............ ";
      eventout += nEBRecHits;
    }
    mehEcaln[1]->Fill((float)nEBRecHits);
  }

  //extract EE information
  edm::Handle<EEUncalibratedRecHitCollection> EcalUncalibRecHitEE;
  iEvent.getByToken(ECalUncalEESrc_Token_, EcalUncalibRecHitEE);
  bool validuncalibRecHitEE = true;
  if (!EcalUncalibRecHitEE.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalUncalRecHitEE in event!";
    validuncalibRecHitEE = false;
  }

  edm::Handle<EERecHitCollection> EcalRecHitEE;
  iEvent.getByToken(ECalEESrc_Token_, EcalRecHitEE);
  bool validRecHitEE = true;
  if (!EcalRecHitEE.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalRecHitEE in event!";
    validRecHitEE = false;
  }

  // loop over simhits
  iEvent.getByToken(EEHits_Token_, crossingFrame);
  validXFrame = true;
  if (!crossingFrame.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find cal endcap crossingFrame in event!";
    validXFrame = false;
  }

  MapType eeSimMap;
  if (validXFrame) {
    const MixCollection<PCaloHit> endcapHits(crossingFrame.product());
    // keep track of sum of simhit energy in each crystal
    for (auto const& iHit : endcapHits) {
      EEDetId eeid = EEDetId(iHit.id());

      uint32_t crystid = eeid.rawId();
      eeSimMap[crystid] += iHit.energy();
    }
  }

  int nEERecHits = 0;
  if (validuncalibRecHitEE && validRecHitEE) {
    // loop over RecHits
    const EEUncalibratedRecHitCollection* EEUncalibRecHit = EcalUncalibRecHitEE.product();
    const EERecHitCollection* EERecHit = EcalRecHitEE.product();

    for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit = EEUncalibRecHit->begin();
         uncalibRecHit != EEUncalibRecHit->end();
         ++uncalibRecHit) {
      EEDetId EEid = EEDetId(uncalibRecHit->id());

      EcalRecHitCollection::const_iterator myRecHit = EERecHit->find(EEid);

      if (myRecHit != EERecHit->end()) {
        ++nEERecHits;
        mehEcalRes[0]->Fill(myRecHit->energy() - eeSimMap[EEid.rawId()]);
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of EERecHits collected:............ ";
      eventout += nEERecHits;
    }
    mehEcaln[0]->Fill((float)nEERecHits);
  }

  //extract ES information
  edm::Handle<ESRecHitCollection> EcalRecHitES;
  iEvent.getByToken(ECalESSrc_Token_, EcalRecHitES);
  bool validRecHitES = true;
  if (!EcalRecHitES.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalRecHitES in event!";
    validRecHitES = false;
  }

  // loop over simhits
  iEvent.getByToken(ESHits_Token_, crossingFrame);
  validXFrame = true;
  if (!crossingFrame.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find cal preshower crossingFrame in event!";
    validXFrame = false;
  }

  MapType esSimMap;
  if (validXFrame) {
    const MixCollection<PCaloHit> preshowerHits(crossingFrame.product());
    // keep track of sum of simhit energy in each crystal
    for (auto const& iHit : preshowerHits) {
      ESDetId esid = ESDetId(iHit.id());

      uint32_t crystid = esid.rawId();
      esSimMap[crystid] += iHit.energy();
    }
  }

  int nESRecHits = 0;
  if (validRecHitES) {
    // loop over RecHits
    const ESRecHitCollection* ESRecHit = EcalRecHitES.product();
    for (EcalRecHitCollection::const_iterator recHit = ESRecHit->begin(); recHit != ESRecHit->end(); ++recHit) {
      ESDetId ESid = ESDetId(recHit->id());

      ++nESRecHits;
      mehEcalRes[2]->Fill(recHit->energy() - esSimMap[ESid.rawId()]);
    }

    if (verbosity > 1) {
      eventout += "\n          Number of ESRecHits collected:............ ";
      eventout += nESRecHits;
    }
    mehEcaln[2]->Fill(float(nESRecHits));
  }

  if (verbosity > 0)
    edm::LogInfo(MsgLoggerCat) << eventout << "\n";

  return;
}

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

private

Definition at line 553 of file GlobalRecHitsAnalyzer.cc.

References dqm::impl::MonitorElement::Fill(), geometry, edm::EventSetup::get(), edm::Event::getByToken(), CaloSubdetectorGeometry::getGeometry(), edm::Event::getManyByType(), HcalGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), photonIsolationHIProducer_cfi::hbhe, DetId::Hcal, HCalSrc_Token_, photonIsolationHIProducer_cfi::hf, photonIsolationHIProducer_cfi::ho, ihf, edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), LogDebug, mehHcaln, mehHcalRes, PV3DBase< T, PVType, FrameType >::phi(), PI, edm::Handle< T >::product(), edm::ESHandle< T >::product(), sdHcalBrl, sdHcalEC, sdHcalFwd, sdHcalOut, HLTBitAnalyser_cfi::simhits, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

Referenced by analyze().

                                                                                        {
  std::string MsgLoggerCat = "GlobalRecHitsAnalyzer_fillHCal";

  TString eventout;
  if (verbosity > 0)
    eventout = "\nGathering info:";

  // get geometry
  edm::ESHandle<CaloGeometry> geometry;
  iSetup.get<CaloGeometryRecord>().get(geometry);
  if (!geometry.isValid()) {
    edm::LogWarning(MsgLoggerCat) << "Unable to find CaloGeometry in event!";
    return;
  }

  // extract simhit info
  edm::Handle<edm::PCaloHitContainer> hcalHits;
  iEvent.getByToken(HCalSrc_Token_, hcalHits);
  bool validhcalHits = true;
  if (!hcalHits.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find hcalHits in event!";
    validhcalHits = false;
  }

  std::map<HcalDetId, float> fHBEnergySimHits;
  std::map<HcalDetId, float> fHEEnergySimHits;
  std::map<HcalDetId, float> fHOEnergySimHits;
  std::map<HcalDetId, float> fHFEnergySimHits;
  if (validhcalHits) {
    const edm::PCaloHitContainer* simhitResult = hcalHits.product();

    for (std::vector<PCaloHit>::const_iterator simhits = simhitResult->begin(); simhits != simhitResult->end();
         ++simhits) {
      HcalDetId detId(simhits->id());

      if (detId.subdet() == sdHcalBrl) {
        fHBEnergySimHits[detId] += simhits->energy();
      }
      if (detId.subdet() == sdHcalEC) {
        fHEEnergySimHits[detId] += simhits->energy();
      }
      if (detId.subdet() == sdHcalOut) {
        fHOEnergySimHits[detId] += simhits->energy();
      }
      if (detId.subdet() == sdHcalFwd) {
        fHFEnergySimHits[detId] += simhits->energy();
      }
    }
  }

  // max values to be used (HO is found in HB)
  Double_t maxHBEnergy = 0.;
  Double_t maxHEEnergy = 0.;
  Double_t maxHFEnergy = 0.;

  Double_t maxHBPhi = -1000.;
  Double_t maxHEPhi = -1000.;
  Double_t maxHOPhi = -1000.;
  Double_t maxHFPhi = -1000.;

  Double_t PI = 3.141592653589;

  // get HBHE information
  std::vector<edm::Handle<HBHERecHitCollection>> hbhe;
  iEvent.getManyByType(hbhe);
  bool validHBHE = true;
  if (!hbhe[0].isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find any HBHERecHitCollections in event!";
    validHBHE = false;
  }

  if (validHBHE) {
    std::vector<edm::Handle<HBHERecHitCollection>>::iterator ihbhe;
    const CaloGeometry* geo = geometry.product();

    int iHB = 0;
    int iHE = 0;
    for (ihbhe = hbhe.begin(); ihbhe != hbhe.end(); ++ihbhe) {
      // find max values
      for (HBHERecHitCollection::const_iterator jhbhe = (*ihbhe)->begin(); jhbhe != (*ihbhe)->end(); ++jhbhe) {
        HcalDetId cell(jhbhe->id());

        if (cell.subdet() == sdHcalBrl) {
          const HcalGeometry* cellGeometry =
              dynamic_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal, cell.subdet()));
          double fPhi = cellGeometry->getPosition(cell).phi();
          if ((jhbhe->energy()) > maxHBEnergy) {
            maxHBEnergy = jhbhe->energy();
            maxHBPhi = fPhi;
            maxHOPhi = maxHBPhi;
          }
        }

        if (cell.subdet() == sdHcalEC) {
          const HcalGeometry* cellGeometry =
              dynamic_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal, cell.subdet()));
          double fPhi = cellGeometry->getPosition(cell).phi();
          if ((jhbhe->energy()) > maxHEEnergy) {
            maxHEEnergy = jhbhe->energy();
            maxHEPhi = fPhi;
          }
        }
      }  // end find max values

      for (HBHERecHitCollection::const_iterator jhbhe = (*ihbhe)->begin(); jhbhe != (*ihbhe)->end(); ++jhbhe) {
        HcalDetId cell(jhbhe->id());

        if (cell.subdet() == sdHcalBrl) {
          ++iHB;

          const HcalGeometry* cellGeometry =
              dynamic_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal, cell.subdet()));
          double fPhi = cellGeometry->getPosition(cell).phi();

          float deltaphi = maxHBPhi - fPhi;
          if (fPhi > maxHBPhi) {
            deltaphi = fPhi - maxHBPhi;
          }
          if (deltaphi > PI) {
            deltaphi = 2.0 * PI - deltaphi;
          }

          mehHcalRes[0]->Fill(jhbhe->energy() - fHBEnergySimHits[cell]);
        }

        if (cell.subdet() == sdHcalEC) {
          ++iHE;

          const HcalGeometry* cellGeometry =
              dynamic_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal, cell.subdet()));
          double fPhi = cellGeometry->getPosition(cell).phi();

          float deltaphi = maxHEPhi - fPhi;
          if (fPhi > maxHEPhi) {
            deltaphi = fPhi - maxHEPhi;
          }
          if (deltaphi > PI) {
            deltaphi = 2.0 * PI - deltaphi;
          }
          mehHcalRes[1]->Fill(jhbhe->energy() - fHEEnergySimHits[cell]);
        }
      }
    }  // end loop through collection

    if (verbosity > 1) {
      eventout += "\n          Number of HBRecHits collected:............ ";
      eventout += iHB;
    }

    if (verbosity > 1) {
      eventout += "\n          Number of HERecHits collected:............ ";
      eventout += iHE;
    }
    mehHcaln[0]->Fill((float)iHB);
    mehHcaln[1]->Fill((float)iHE);
  }

  // get HF information
  std::vector<edm::Handle<HFRecHitCollection>> hf;
  iEvent.getManyByType(hf);
  bool validHF = true;
  if (!hf[0].isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find any HFRecHitCollections in event!";
    validHF = false;
  }
  if (validHF) {
    std::vector<edm::Handle<HFRecHitCollection>>::iterator ihf;

    int iHF = 0;
    for (ihf = hf.begin(); ihf != hf.end(); ++ihf) {
      // find max values
      for (HFRecHitCollection::const_iterator jhf = (*ihf)->begin(); jhf != (*ihf)->end(); ++jhf) {
        HcalDetId cell(jhf->id());

        if (cell.subdet() == sdHcalFwd) {
          auto cellGeometry = geometry->getSubdetectorGeometry(cell)->getGeometry(cell);
          double fPhi = cellGeometry->getPosition().phi();
          if ((jhf->energy()) > maxHFEnergy) {
            maxHFEnergy = jhf->energy();
            maxHFPhi = fPhi;
          }
        }
      }  // end find max values

      for (HFRecHitCollection::const_iterator jhf = (*ihf)->begin(); jhf != (*ihf)->end(); ++jhf) {
        HcalDetId cell(jhf->id());

        if (cell.subdet() == sdHcalFwd) {
          ++iHF;

          auto cellGeometry = geometry->getSubdetectorGeometry(cell)->getGeometry(cell);
          double fPhi = cellGeometry->getPosition().phi();

          float deltaphi = maxHBPhi - fPhi;
          if (fPhi > maxHFPhi) {
            deltaphi = fPhi - maxHFPhi;
          }
          if (deltaphi > PI) {
            deltaphi = 2.0 * PI - deltaphi;
          }

          mehHcalRes[2]->Fill(jhf->energy() - fHFEnergySimHits[cell]);
        }
      }
    }  // end loop through collection

    if (verbosity > 1) {
      eventout += "\n          Number of HFDigis collected:.............. ";
      eventout += iHF;
    }
    mehHcaln[2]->Fill((float)iHF);
  }

  // get HO information
  std::vector<edm::Handle<HORecHitCollection>> ho;
  iEvent.getManyByType(ho);
  bool validHO = true;
  if (!ho[0].isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find any HORecHitCollections in event!";
    validHO = false;
  }

  if (validHO) {
    std::vector<edm::Handle<HORecHitCollection>>::iterator iho;

    int iHO = 0;
    for (iho = ho.begin(); iho != ho.end(); ++iho) {
      for (HORecHitCollection::const_iterator jho = (*iho)->begin(); jho != (*iho)->end(); ++jho) {
        HcalDetId cell(jho->id());

        if (cell.subdet() == sdHcalOut) {
          ++iHO;

          auto cellGeometry = geometry->getSubdetectorGeometry(cell)->getGeometry(cell);
          double fPhi = cellGeometry->getPosition().phi();

          float deltaphi = maxHOPhi - fPhi;
          if (fPhi > maxHOPhi) {
            deltaphi = fPhi - maxHOPhi;
          }
          if (deltaphi > PI) {
            deltaphi = 2.0 * PI - deltaphi;
          }
          mehHcalRes[3]->Fill(jho->energy() - fHOEnergySimHits[cell]);
        }
      }
    }  // end loop through collection

    if (verbosity > 1) {
      eventout += "\n          Number of HODigis collected:.............. ";
      eventout += iHO;
    }
    mehHcaln[3]->Fill((float)iHO);
  }

  if (verbosity > 0)
    edm::LogInfo(MsgLoggerCat) << eventout << "\n";

  return;
}

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

private

Definition at line 1183 of file GlobalRecHitsAnalyzer.cc.

References CSCLayer::chamber(), CSCChamberSpecs::chamberType(), compute(), dqm::impl::MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, CSCGeometry::idToDetUnit(), RPCRoll::isForward(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), LogDebug, map1DRecHitsPerWire(), DTHitQualityUtils::mapSimHitsPerWire(), mehCSCn, mehDtMuonn, mehRPCn, mehRPCResX, MuCSCHits_Token_, MuCSCSrc_Token_, MuDTSimSrc_Token_, MuDTSrc_Token_, MuRPCSimSrc_Token_, MuRPCSrc_Token_, plotResolution(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), OfflineHarvestingSequence_cosmic::ptype, alignCSCRings::r, rpcPointValidation_cfi::recHit, RPCGeometry::roll(), rpcPointValidation_cfi::simHit, FastTrackerRecHitCombiner_cfi::simHits, CSCChamber::specs(), AlCaHLTBitMon_QueryRunRegistry::string, theMap, verbosity, and PV3DBase< T, PVType, FrameType >::x().

Referenced by analyze().

                                                                                        {
  std::string MsgLoggerCat = "GlobalRecHitsAnalyzer_fillMuon";

  TString eventout;
  if (verbosity > 0)
    eventout = "\nGathering info:";

  // get DT information
  edm::ESHandle<DTGeometry> dtGeom;
  iSetup.get<MuonGeometryRecord>().get(dtGeom);
  if (!dtGeom.isValid()) {
    edm::LogWarning(MsgLoggerCat) << "Unable to find DTMuonGeometryRecord in event!";
    return;
  }

  edm::Handle<edm::PSimHitContainer> dtsimHits;
  iEvent.getByToken(MuDTSimSrc_Token_, dtsimHits);
  bool validdtsim = true;
  if (!dtsimHits.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find dtsimHits in event!";
    validdtsim = false;
  }

  edm::Handle<DTRecHitCollection> dtRecHits;
  iEvent.getByToken(MuDTSrc_Token_, dtRecHits);
  bool validdtrec = true;
  if (!dtRecHits.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find dtRecHits in event!";
    validdtrec = false;
  }

  if (validdtsim && validdtrec) {
    std::map<DTWireId, edm::PSimHitContainer> simHitsPerWire =
        DTHitQualityUtils::mapSimHitsPerWire(*(dtsimHits.product()));

    std::map<DTWireId, std::vector<DTRecHit1DPair>> recHitsPerWire = map1DRecHitsPerWire(dtRecHits.product());

    int nDt = compute(dtGeom.product(), simHitsPerWire, recHitsPerWire, 1);

    if (verbosity > 1) {
      eventout += "\n          Number of DtMuonRecHits collected:........ ";
      eventout += nDt;
    }
    mehDtMuonn->Fill(float(nDt));
  }

  // get CSC Strip information
  // get map of sim hits
  theMap.clear();
  edm::Handle<CrossingFrame<PSimHit>> cf;

  iEvent.getByToken(MuCSCHits_Token_, cf);
  bool validXFrame = true;
  if (!cf.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find muo CSC crossingFrame in event!";
    validXFrame = false;
  }
  if (validXFrame) {
    const MixCollection<PSimHit> simHits(cf.product());

    // arrange the hits by detUnit
    for (auto const& iHit : simHits) {
      theMap[iHit.detUnitId()].push_back(iHit);
    }
  }

  // get geometry
  edm::ESHandle<CSCGeometry> hGeom;
  iSetup.get<MuonGeometryRecord>().get(hGeom);
  if (!hGeom.isValid()) {
    edm::LogWarning(MsgLoggerCat) << "Unable to find CSCMuonGeometryRecord in event!";
    return;
  }
  const CSCGeometry* theCSCGeometry = &*hGeom;

  // get rechits
  edm::Handle<CSCRecHit2DCollection> hRecHits;
  iEvent.getByToken(MuCSCSrc_Token_, hRecHits);
  bool validCSC = true;
  if (!hRecHits.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find CSC RecHits in event!";
    validCSC = false;
  }

  if (validCSC) {
    const CSCRecHit2DCollection* cscRecHits = hRecHits.product();

    int nCSC = 0;
    for (CSCRecHit2DCollection::const_iterator recHitItr = cscRecHits->begin(); recHitItr != cscRecHits->end();
         ++recHitItr) {
      int detId = (*recHitItr).cscDetId().rawId();

      edm::PSimHitContainer simHits;
      std::map<int, edm::PSimHitContainer>::const_iterator mapItr = theMap.find(detId);
      if (mapItr != theMap.end()) {
        simHits = mapItr->second;
      }

      if (simHits.size() == 1) {
        ++nCSC;

        const GeomDetUnit* detUnit = theCSCGeometry->idToDetUnit(CSCDetId(detId));
        const CSCLayer* layer = dynamic_cast<const CSCLayer*>(detUnit);

        int chamberType = layer->chamber()->specs()->chamberType();
        plotResolution(simHits[0], *recHitItr, layer, chamberType);
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of CSCRecHits collected:........... ";
      eventout += nCSC;
    }
    mehCSCn->Fill((float)nCSC);
  }

  // get RPC information
  std::map<double, int> mapsim, maprec;
  std::map<int, double> nmapsim, nmaprec;

  edm::ESHandle<RPCGeometry> rpcGeom;
  iSetup.get<MuonGeometryRecord>().get(rpcGeom);
  if (!rpcGeom.isValid()) {
    edm::LogWarning(MsgLoggerCat) << "Unable to find RPCMuonGeometryRecord in event!";
    return;
  }

  edm::Handle<edm::PSimHitContainer> simHit;
  iEvent.getByToken(MuRPCSimSrc_Token_, simHit);
  bool validrpcsim = true;
  if (!simHit.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find RPCSimHit in event!";
    validrpcsim = false;
  }

  edm::Handle<RPCRecHitCollection> recHit;
  iEvent.getByToken(MuRPCSrc_Token_, recHit);
  bool validrpc = true;
  if (!simHit.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find RPCRecHit in event!";
    validrpc = false;
  }

  if (validrpc) {
    int nRPC = 0;
    RPCRecHitCollection::const_iterator recIt;
    int nrec = 0;
    for (recIt = recHit->begin(); recIt != recHit->end(); ++recIt) {
      RPCDetId Rid = (RPCDetId)(*recIt).rpcId();
      const RPCRoll* roll = dynamic_cast<const RPCRoll*>(rpcGeom->roll(Rid));
      if (roll->isForward()) {
        if (verbosity > 1) {
          eventout += "\n          Number of RPCRecHits collected:........... ";
          eventout += nRPC;
        }

        if (verbosity > 0)
          edm::LogInfo(MsgLoggerCat) << eventout << "\n";
        return;
      }
      nrec = nrec + 1;
      LocalPoint rhitlocal = (*recIt).localPosition();
      double rhitlocalx = rhitlocal.x();
      maprec[rhitlocalx] = nrec;
    }

    int i = 0;
    for (std::map<double, int>::iterator iter = maprec.begin(); iter != maprec.end(); ++iter) {
      i = i + 1;
      nmaprec[i] = (*iter).first;
    }

    int nsim = 0;
    if (validrpcsim) {
      edm::PSimHitContainer::const_iterator simIt;
      for (simIt = simHit->begin(); simIt != simHit->end(); simIt++) {
        int ptype = (*simIt).particleType();
        if (ptype == 13 || ptype == -13) {
          nsim = nsim + 1;
          LocalPoint shitlocal = (*simIt).localPosition();
          double shitlocalx = shitlocal.x();
          mapsim[shitlocalx] = nsim;
        }
      }

      i = 0;
      for (std::map<double, int>::iterator iter = mapsim.begin(); iter != mapsim.end(); ++iter) {
        i = i + 1;
        nmapsim[i] = (*iter).first;
      }
    }

    if (nsim == nrec) {
      for (int r = 0; r < nsim; r++) {
        ++nRPC;
        mehRPCResX->Fill(nmaprec[r + 1] - nmapsim[r + 1]);
      }
    }

    if (verbosity > 1) {
      eventout += "\n          Number of RPCRecHits collected:........... ";
      eventout += nRPC;
    }
    mehRPCn->Fill((float)nRPC);
  }

  if (verbosity > 0)
    edm::LogInfo(MsgLoggerCat) << eventout << "\n";

  return;
}

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

private

Definition at line 822 of file GlobalRecHitsAnalyzer.cc.

References TrackerHitAssociator::associateHit(), edmNew::DetSet< T >::begin(), TrackerGeometry::dets(), edmNew::DetSet< T >::end(), edmNew::DetSetVector< T >::end(), dqm::impl::MonitorElement::Fill(), edmNew::DetSetVector< T >::find(), TrackingRecHit::geographicalId(), relativeConstraints::geom, edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, TrackerGeometry::idToDet(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), BaseTrackerRecHit::localPosition(), LogDebug, visualization-live-secondInstance_cfg::m, matched, mehSiPixeln, mehSiPixelResX, mehSiPixelResY, mehSiStripn, mehSiStripResX, mehSiStripResY, position, edm::ESHandle< T >::product(), projectHit(), TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), TrackerTopology::pxfSide(), sdPxlBrl, sdPxlFwd, sdSiTEC, sdSiTIB, sdSiTID, sdSiTOB, SiPxlSrc_Token_, SiStripSrc_Token_, mathSSE::sqrt(), GluedGeomDet::stereoDet(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), GeomDet::surface(), TrackerTopology::tecWheel(), TrackerTopology::tibLayer(), TrackerTopology::tidWheel(), TrackerTopology::tobLayer(), PbPb_ZMuSkimMuonDPG_cff::tracker, trackerHitAssociatorConfig_, verbosity, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by analyze().

                                                                                       {
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle;
  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
  const TrackerTopology* const tTopo = tTopoHandle.product();

  std::string MsgLoggerCat = "GlobalRecHitsAnalyzer_fillTrk";

  TString eventout;
  if (verbosity > 0)
    eventout = "\nGathering info:";

  // get strip information
  edm::Handle<SiStripMatchedRecHit2DCollection> rechitsmatched;
  iEvent.getByToken(SiStripSrc_Token_, rechitsmatched);
  bool validstrip = true;
  if (!rechitsmatched.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find stripmatchedrechits in event!";
    validstrip = false;
  }

  TrackerHitAssociator associate(iEvent, trackerHitAssociatorConfig_);

  edm::ESHandle<TrackerGeometry> pDD;
  iSetup.get<TrackerDigiGeometryRecord>().get(pDD);
  if (!pDD.isValid()) {
    edm::LogWarning(MsgLoggerCat) << "Unable to find TrackerDigiGeometry in event!";
    return;
  }
  const TrackerGeometry& tracker(*pDD);

  if (validstrip) {
    int nStripBrl = 0, nStripFwd = 0;

    // loop over det units
    for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); ++it) {
      uint32_t myid = ((*it)->geographicalId()).rawId();
      DetId detid = ((*it)->geographicalId());

      //loop over rechits-matched in the same subdetector
      SiStripMatchedRecHit2DCollection::const_iterator rechitmatchedMatch = rechitsmatched->find(detid);

      if (rechitmatchedMatch != rechitsmatched->end()) {
        SiStripMatchedRecHit2DCollection::DetSet rechitmatchedRange = *rechitmatchedMatch;
        SiStripMatchedRecHit2DCollection::DetSet::const_iterator rechitmatchedRangeIteratorBegin =
            rechitmatchedRange.begin();
        SiStripMatchedRecHit2DCollection::DetSet::const_iterator rechitmatchedRangeIteratorEnd =
            rechitmatchedRange.end();
        SiStripMatchedRecHit2DCollection::DetSet::const_iterator itermatched = rechitmatchedRangeIteratorBegin;

        for (itermatched = rechitmatchedRangeIteratorBegin; itermatched != rechitmatchedRangeIteratorEnd;
             ++itermatched) {
          SiStripMatchedRecHit2D const rechit = *itermatched;
          LocalPoint position = rechit.localPosition();

          float mindist = 999999.;
          float distx = 999999.;
          float disty = 999999.;
          float dist = 999999.;
          std::pair<LocalPoint, LocalVector> closestPair;
          matched.clear();

          float rechitmatchedx = position.x();
          float rechitmatchedy = position.y();

          matched = associate.associateHit(rechit);

          if (!matched.empty()) {
            //project simhit;
            const GluedGeomDet* gluedDet = (const GluedGeomDet*)tracker.idToDet(rechit.geographicalId());
            const StripGeomDetUnit* partnerstripdet = (StripGeomDetUnit*)gluedDet->stereoDet();
            std::pair<LocalPoint, LocalVector> hitPair;

            for (std::vector<PSimHit>::const_iterator m = matched.begin(); m != matched.end(); m++) {
              //project simhit;
              hitPair = projectHit((*m), partnerstripdet, gluedDet->surface());
              distx = fabs(rechitmatchedx - hitPair.first.x());
              disty = fabs(rechitmatchedy - hitPair.first.y());
              dist = sqrt(distx * distx + disty * disty);

              if (dist < mindist) {
                mindist = dist;
                closestPair = hitPair;
              }
            }

            // get TIB
            if (detid.subdetId() == sdSiTIB) {
              ++nStripBrl;

              if (tTopo->tibLayer(myid) == 1) {
                mehSiStripResX[8]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[8]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tibLayer(myid) == 2) {
                mehSiStripResX[9]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[9]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tibLayer(myid) == 3) {
                mehSiStripResX[10]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[10]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tibLayer(myid) == 4) {
                mehSiStripResX[11]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[11]->Fill(rechitmatchedy - closestPair.first.y());
              }
            }

            // get TOB
            if (detid.subdetId() == sdSiTOB) {
              ++nStripBrl;

              if (tTopo->tobLayer(myid) == 1) {
                mehSiStripResX[15]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[15]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tobLayer(myid) == 2) {
                mehSiStripResX[16]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[16]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tobLayer(myid) == 3) {
                mehSiStripResX[17]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[17]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tobLayer(myid) == 4) {
                mehSiStripResX[18]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[18]->Fill(rechitmatchedy - closestPair.first.y());
              }
            }

            // get TID
            if (detid.subdetId() == sdSiTID) {
              ++nStripFwd;

              if (tTopo->tidWheel(myid) == 1) {
                mehSiStripResX[12]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[12]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tidWheel(myid) == 2) {
                mehSiStripResX[13]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[13]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tidWheel(myid) == 3) {
                mehSiStripResX[14]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[14]->Fill(rechitmatchedy - closestPair.first.y());
              }
            }

            // get TEC
            if (detid.subdetId() == sdSiTEC) {
              ++nStripFwd;

              if (tTopo->tecWheel(myid) == 1) {
                mehSiStripResX[0]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[0]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tecWheel(myid) == 2) {
                mehSiStripResX[1]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[1]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tecWheel(myid) == 3) {
                mehSiStripResX[2]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[2]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tecWheel(myid) == 4) {
                mehSiStripResX[3]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[3]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tecWheel(myid) == 5) {
                mehSiStripResX[4]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[4]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tecWheel(myid) == 6) {
                mehSiStripResX[5]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[5]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tecWheel(myid) == 7) {
                mehSiStripResX[6]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[6]->Fill(rechitmatchedy - closestPair.first.y());
              }
              if (tTopo->tecWheel(myid) == 8) {
                mehSiStripResX[7]->Fill(rechitmatchedx - closestPair.first.x());
                mehSiStripResY[7]->Fill(rechitmatchedy - closestPair.first.y());
              }
            }

          }  // end if matched empty
        }
      }
    }  // end loop over det units

    if (verbosity > 1) {
      eventout += "\n          Number of BrlStripRecHits collected:...... ";
      eventout += nStripBrl;
    }

    for (int i = 8; i < 12; ++i) {
      mehSiStripn[i]->Fill((float)nStripBrl);
    }
    for (int i = 16; i < 19; ++i) {
      mehSiStripn[i]->Fill((float)nStripBrl);
    }

    if (verbosity > 1) {
      eventout += "\n          Number of FrwdStripRecHits collected:..... ";
      eventout += nStripFwd;
    }
    for (int i = 0; i < 8; ++i) {
      mehSiStripn[i]->Fill((float)nStripFwd);
    }
    for (int i = 12; i < 16; ++i) {
      mehSiStripn[i]->Fill((float)nStripFwd);
    }
  }

  // get pixel information
  //Get RecHits
  edm::Handle<SiPixelRecHitCollection> recHitColl;
  iEvent.getByToken(SiPxlSrc_Token_, recHitColl);
  bool validpixel = true;
  if (!recHitColl.isValid()) {
    LogDebug(MsgLoggerCat) << "Unable to find SiPixelRecHitCollection in event!";
    validpixel = false;
  }

  //Get event setup
  edm::ESHandle<TrackerGeometry> geom;
  iSetup.get<TrackerDigiGeometryRecord>().get(geom);
  if (!geom.isValid()) {
    edm::LogWarning(MsgLoggerCat) << "Unable to find TrackerDigiGeometry in event!";
    return;
  }

  if (validpixel) {
    int nPxlBrl = 0, nPxlFwd = 0;
    //iterate over detunits
    for (TrackerGeometry::DetContainer::const_iterator it = geom->dets().begin(); it != geom->dets().end(); ++it) {
      uint32_t myid = ((*it)->geographicalId()).rawId();
      DetId detId = ((*it)->geographicalId());
      int subid = detId.subdetId();

      if (!((subid == sdPxlBrl) || (subid == sdPxlFwd)))
        continue;

      SiPixelRecHitCollection::const_iterator pixeldet = recHitColl->find(detId);
      if (pixeldet == recHitColl->end())
        continue;
      SiPixelRecHitCollection::DetSet pixelrechitRange = *pixeldet;
      SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.begin();
      SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.end();
      SiPixelRecHitCollection::DetSet::const_iterator pixeliter = pixelrechitRangeIteratorBegin;

      std::vector<PSimHit> matched;

      //----Loop over rechits for this detId
      for (; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter) {
        matched.clear();
        matched = associate.associateHit(*pixeliter);

        if (!matched.empty()) {
          float closest = 9999.9;
          LocalPoint lp = pixeliter->localPosition();
          float rechit_x = lp.x();
          float rechit_y = lp.y();

          float sim_x = 0.;
          float sim_y = 0.;

          //loop over sim hits and fill closet
          for (std::vector<PSimHit>::const_iterator m = matched.begin(); m != matched.end(); ++m) {
            float sim_x1 = (*m).entryPoint().x();
            float sim_x2 = (*m).exitPoint().x();
            float sim_xpos = 0.5 * (sim_x1 + sim_x2);

            float sim_y1 = (*m).entryPoint().y();
            float sim_y2 = (*m).exitPoint().y();
            float sim_ypos = 0.5 * (sim_y1 + sim_y2);

            float x_res = fabs(sim_xpos - rechit_x);
            float y_res = fabs(sim_ypos - rechit_y);

            float dist = sqrt(x_res * x_res + y_res * y_res);

            if (dist < closest) {
              closest = dist;
              sim_x = sim_xpos;
              sim_y = sim_ypos;
            }
          }  // end sim hit loop

          // get Barrel pixels ***************Pixel STuff******************
          if (subid == sdPxlBrl) {
            ++nPxlBrl;

            if (tTopo->pxbLayer(myid) == 1) {
              mehSiPixelResX[0]->Fill(rechit_x - sim_x);
              mehSiPixelResY[0]->Fill(rechit_y - sim_y);
            }
            if (tTopo->pxbLayer(myid) == 2) {
              mehSiPixelResX[1]->Fill(rechit_x - sim_x);
              mehSiPixelResY[1]->Fill(rechit_y - sim_y);
            }
            if (tTopo->pxbLayer(myid) == 3) {
              mehSiPixelResX[2]->Fill(rechit_x - sim_x);
              mehSiPixelResY[2]->Fill(rechit_y - sim_y);
            }
          }

          // get Forward pixels
          if (subid == sdPxlFwd) {
            ++nPxlFwd;

            if (tTopo->pxfDisk(myid) == 1) {
              if (tTopo->pxfSide(myid) == 1) {
                mehSiPixelResX[3]->Fill(rechit_x - sim_x);
                mehSiPixelResY[3]->Fill(rechit_y - sim_y);
              }
              if (tTopo->pxfSide(myid) == 2) {
                mehSiPixelResX[4]->Fill(rechit_x - sim_x);
                mehSiPixelResY[4]->Fill(rechit_y - sim_y);
              }
            }
            if (tTopo->pxfDisk(myid) == 2) {
              if (tTopo->pxfSide(myid) == 1) {
                mehSiPixelResX[5]->Fill(rechit_x - sim_x);
                mehSiPixelResY[5]->Fill(rechit_y - sim_y);
              }
              if (tTopo->pxfSide(myid) == 2) {
                mehSiPixelResX[6]->Fill(rechit_x - sim_x);
                mehSiPixelResY[6]->Fill(rechit_y - sim_y);
              }
            }
          }
        }  // end matched emtpy
      }    // <-----end rechit loop
    }      // <------ end detunit loop

    if (verbosity > 1) {
      eventout += "\n          Number of BrlPixelRecHits collected:...... ";
      eventout += nPxlBrl;
    }
    for (int i = 0; i < 3; ++i) {
      mehSiPixeln[i]->Fill((float)nPxlBrl);
    }

    if (verbosity > 1) {
      eventout += "\n          Number of FrwdPixelRecHits collected:..... ";
      eventout += nPxlFwd;
    }

    for (int i = 3; i < 7; ++i) {
      mehSiPixeln[i]->Fill((float)nPxlFwd);
    }
  }

  if (verbosity > 0)
    edm::LogInfo(MsgLoggerCat) << eventout << "\n";

  return;
}

template<typename type >

const type * GlobalRecHitsAnalyzer::findBestRecHit ( const DTLayer *  layer, DTWireId  wireId, const std::vector< type > &  recHits, const float  simHitDist  )

private

Definition at line 1449 of file GlobalRecHitsAnalyzer.cc.

References rpcPointValidation_cfi::recHit, and recHitDistFromWire().

Referenced by compute().

                                                                          {
  float res = 99999;
  const type* theBestRecHit = nullptr;
  // Loop over RecHits within the cell
  for (typename std::vector<type>::const_iterator recHit = recHits.begin(); recHit != recHits.end(); recHit++) {
    float distTmp = recHitDistFromWire(*recHit, layer);
    if (fabs(distTmp - simHitDist) < res) {
      res = fabs(distTmp - simHitDist);
      theBestRecHit = &(*recHit);
    }
  }  // End of loop over RecHits within the cell

  return theBestRecHit;
}

std::map< DTWireId, std::vector< DTRecHit1DPair > > GlobalRecHitsAnalyzer::map1DRecHitsPerWire ( const DTRecHitCollection *  dt1DRecHitPairs )

private

Definition at line 1423 of file GlobalRecHitsAnalyzer.cc.

References runTheMatrix::ret.

Referenced by fillMuon().

                                               {
  std::map<DTWireId, std::vector<DTRecHit1DPair>> ret;

  for (DTRecHitCollection::const_iterator rechit = dt1DRecHitPairs->begin(); rechit != dt1DRecHitPairs->end();
       rechit++) {
    ret[(*rechit).wireId()].push_back(*rechit);
  }

  return ret;
}

void GlobalRecHitsAnalyzer::plotResolution ( const PSimHit &  simHit, const CSCRecHit2D &  recHit, const CSCLayer *  layer, int  chamberType  )

private

Definition at line 1530 of file GlobalRecHitsAnalyzer.cc.

References dqm::impl::MonitorElement::Fill(), PSimHit::localPosition(), CSCRecHit2D::localPosition(), mehCSCResRDPhi, PV3DBase< T, PVType, FrameType >::phi(), and GeomDet::toGlobal().

Referenced by fillMuon().

                                                            {
  GlobalPoint simHitPos = layer->toGlobal(simHit.localPosition());
  GlobalPoint recHitPos = layer->toGlobal(recHit.localPosition());

  mehCSCResRDPhi->Fill(recHitPos.phi() - simHitPos.phi());
}

std::pair< LocalPoint, LocalVector > GlobalRecHitsAnalyzer::projectHit ( const PSimHit &  hit, const StripGeomDetUnit *  stripDet, const BoundPlane &  plane  )

private

Definition at line 1396 of file GlobalRecHitsAnalyzer.cc.

References funct::cos(), DeadROC_duringRun::dir, PSimHit::localDirection(), PSimHit::localPosition(), Scenarios_cff::scale, funct::sin(), StripGeomDetUnit::specificTopology(), StripTopology::strip(), StripTopology::stripAngle(), GeomDet::surface(), Surface::toGlobal(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by fillTrk().

                                                                                              {
  const StripTopology& topol = stripDet->specificTopology();
  GlobalPoint globalpos = stripDet->surface().toGlobal(hit.localPosition());
  LocalPoint localHit = plane.toLocal(globalpos);
  //track direction
  LocalVector locdir = hit.localDirection();
  //rotate track in new frame

  GlobalVector globaldir = stripDet->surface().toGlobal(locdir);
  LocalVector dir = plane.toLocal(globaldir);
  float scale = -localHit.z() / dir.z();

  LocalPoint projectedPos = localHit + scale * dir;

  float selfAngle = topol.stripAngle(topol.strip(hit.localPosition()));

  // vector along strip in hit frame
  LocalVector stripDir(sin(selfAngle), cos(selfAngle), 0);

  LocalVector localStripDir(plane.toLocal(stripDet->surface().toGlobal(stripDir)));

  return std::pair<LocalPoint, LocalVector>(projectedPos, localStripDir);
}

float GlobalRecHitsAnalyzer::recHitDistFromWire ( const DTRecHit1DPair &  hitPair, const DTLayer *  layer  )

private

Definition at line 1468 of file GlobalRecHitsAnalyzer.cc.

References DTEnums::Left, DTRecHit1DPair::localPosition(), DTEnums::Right, and PV3DBase< T, PVType, FrameType >::x().

Referenced by compute(), and findBestRecHit().

                                                                                                   {
  // Compute the rechit distance from wire
  return fabs(hitPair.localPosition(DTEnums::Left).x() - hitPair.localPosition(DTEnums::Right).x()) / 2.;
}

float GlobalRecHitsAnalyzer::recHitDistFromWire ( const DTRecHit1D &  recHit, const DTLayer *  layer  )

private

Definition at line 1474 of file GlobalRecHitsAnalyzer.cc.

References DTRecHit1D::localPosition(), DTLayer::specificTopology(), DTWireId::wire(), DTRecHit1D::wireId(), DTTopology::wirePosition(), and PV3DBase< T, PVType, FrameType >::x().

                                                                                              {
  return fabs(recHit.localPosition().x() - layer->specificTopology().wirePosition(recHit.wireId().wire()));
}

float GlobalRecHitsAnalyzer::simHitDistFromWire ( const DTLayer *  layer, DTWireId  wireId, const PSimHit &  hit  )

private

Definition at line 1436 of file GlobalRecHitsAnalyzer.cc.

References PSimHit::entryPoint(), PSimHit::exitPoint(), DTLayer::specificTopology(), DTWireId::wire(), DTTopology::wirePosition(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by compute().

                                                                                                         {
  float xwire = layer->specificTopology().wirePosition(wireId.wire());
  LocalPoint entryP = hit.entryPoint();
  LocalPoint exitP = hit.exitPoint();
  float xEntry = entryP.x() - xwire;
  float xExit = exitP.x() - xwire;

  //FIXME: check...
  return fabs(xEntry - (entryP.z() * (xExit - xEntry)) / (exitP.z() - entryP.z()));
}

Member Data Documentation

unsigned int GlobalRecHitsAnalyzer::count

private

Definition at line 287 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<CrossingFrame<PCaloHit> > GlobalRecHitsAnalyzer::EBHits_Token_

private

Definition at line 190 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::ECalEBSrc_

private

Definition at line 180 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<EBRecHitCollection> GlobalRecHitsAnalyzer::ECalEBSrc_Token_

private

Definition at line 185 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::ECalEESrc_

private

Definition at line 182 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<EERecHitCollection> GlobalRecHitsAnalyzer::ECalEESrc_Token_

private

Definition at line 186 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::ECalESSrc_

private

Definition at line 184 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<ESRecHitCollection> GlobalRecHitsAnalyzer::ECalESSrc_Token_

private

Definition at line 187 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::ECalUncalEBSrc_

private

Definition at line 181 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<EBUncalibratedRecHitCollection> GlobalRecHitsAnalyzer::ECalUncalEBSrc_Token_

private

Definition at line 188 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::ECalUncalEESrc_

private

Definition at line 183 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<EEUncalibratedRecHitCollection> GlobalRecHitsAnalyzer::ECalUncalEESrc_Token_

private

Definition at line 189 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::EDGetTokenT<CrossingFrame<PCaloHit> > GlobalRecHitsAnalyzer::EEHits_Token_

private

Definition at line 191 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::EDGetTokenT<CrossingFrame<PCaloHit> > GlobalRecHitsAnalyzer::ESHits_Token_

private

Definition at line 192 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

std::string GlobalRecHitsAnalyzer::fName

private

Definition at line 166 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

int GlobalRecHitsAnalyzer::frequency

private

Definition at line 168 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze(), and GlobalRecHitsAnalyzer().

bool GlobalRecHitsAnalyzer::getAllProvenances

private

Definition at line 170 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::HCalSrc_

private

Definition at line 199 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> GlobalRecHitsAnalyzer::HCalSrc_Token_

private

Definition at line 200 of file GlobalRecHitsAnalyzer.h.

Referenced by fillHCal(), and GlobalRecHitsAnalyzer().

std::string GlobalRecHitsAnalyzer::hitsProducer

private

Definition at line 172 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

std::string GlobalRecHitsAnalyzer::label

private

Definition at line 169 of file GlobalRecHitsAnalyzer.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::_sort_list(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::foundIn(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::fullFilename(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputEventContent(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), Vispa.Plugins.ConfigEditor.ToolDataAccessor.ToolDataAccessor::properties(), Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::properties(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::readConnections(), Vispa.Plugins.ConfigEditor.ToolDataAccessor.ToolDataAccessor::updateProcess(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::usedBy().

std::vector<PSimHit> GlobalRecHitsAnalyzer::matched

private

Definition at line 212 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehCSCn

private

Definition at line 231 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalRecHitsAnalyzer::mehCSCResRDPhi

private

Definition at line 234 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and plotResolution().

MonitorElement* GlobalRecHitsAnalyzer::mehDtMuonn

private

Definition at line 230 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalRecHitsAnalyzer::mehDtMuonRes

private

Definition at line 233 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and compute().

MonitorElement* GlobalRecHitsAnalyzer::mehEcaln[3]

private

Definition at line 177 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalRecHitsAnalyzer::mehEcalRes[3]

private

Definition at line 178 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalRecHitsAnalyzer::mehHcaln[4]

private

Definition at line 196 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

MonitorElement* GlobalRecHitsAnalyzer::mehHcalRes[4]

private

Definition at line 197 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

MonitorElement* GlobalRecHitsAnalyzer::mehRPCn

private

Definition at line 232 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalRecHitsAnalyzer::mehRPCResX

private

Definition at line 235 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalRecHitsAnalyzer::mehSiPixeln[7]

private

Definition at line 220 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehSiPixelResX[7]

private

Definition at line 221 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehSiPixelResY[7]

private

Definition at line 222 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehSiStripn[19]

private

Definition at line 205 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehSiStripResX[19]

private

Definition at line 206 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehSiStripResY[19]

private

Definition at line 207 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

edm::EDGetTokenT<CrossingFrame<PSimHit> > GlobalRecHitsAnalyzer::MuCSCHits_Token_

private

Definition at line 272 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::MuCSCSrc_

private

Definition at line 270 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<CSCRecHit2DCollection> GlobalRecHitsAnalyzer::MuCSCSrc_Token_

private

Definition at line 271 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::MuDTSimSrc_

private

Definition at line 238 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalRecHitsAnalyzer::MuDTSimSrc_Token_

private

Definition at line 240 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::MuDTSrc_

private

Definition at line 237 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<DTRecHitCollection> GlobalRecHitsAnalyzer::MuDTSrc_Token_

private

Definition at line 239 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::MuRPCSimSrc_

private

Definition at line 282 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<edm::PSimHitContainer> GlobalRecHitsAnalyzer::MuRPCSimSrc_Token_

private

Definition at line 284 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::MuRPCSrc_

private

Definition at line 281 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<RPCRecHitCollection> GlobalRecHitsAnalyzer::MuRPCSrc_Token_

private

Definition at line 283 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

bool GlobalRecHitsAnalyzer::printProvenanceInfo

private

Definition at line 171 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::SiPxlSrc_

private

Definition at line 224 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<SiPixelRecHitCollection> GlobalRecHitsAnalyzer::SiPxlSrc_Token_

private

Definition at line 225 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::SiStripSrc_

private

Definition at line 209 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<SiStripMatchedRecHit2DCollection> GlobalRecHitsAnalyzer::SiStripSrc_Token_

private

Definition at line 210 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalRecHitsAnalyzer().

std::map<int, edm::PSimHitContainer> GlobalRecHitsAnalyzer::theMap

private

Definition at line 274 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon().

TrackerHitAssociator::Config GlobalRecHitsAnalyzer::trackerHitAssociatorConfig_

private

Definition at line 216 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk().

int GlobalRecHitsAnalyzer::verbosity

private

Definition at line 167 of file GlobalRecHitsAnalyzer.h.

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