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::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

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
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeomToken_
unsigned int	count
edm::ESGetToken< CSCGeometry, MuonGeometryRecord >	cscGeomToken_
edm::ESGetToken< DTGeometry, MuonGeometryRecord >	dtGeomToken_
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::GetterOfProducts< edm::SortedCollection< HBHERecHit, edm::StrictWeakOrdering< HBHERecHit > > >	HBHERecHitgetter_
edm::InputTag	HCalSrc_
edm::EDGetTokenT< edm::PCaloHitContainer >	HCalSrc_Token_
edm::GetterOfProducts< edm::SortedCollection< HFRecHit, edm::StrictWeakOrdering< HFRecHit > > >	HFRecHitgetter_
std::string	hitsProducer
edm::GetterOfProducts< edm::SortedCollection< HORecHit, edm::StrictWeakOrdering< HORecHit > > >	HORecHitgetter_
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::ESGetToken< RPCGeometry, MuonGeometryRecord >	rpcGeomToken_
edm::InputTag	SiPxlSrc_
edm::EDGetTokenT< SiPixelRecHitCollection >	SiPxlSrc_Token_
edm::InputTag	SiStripSrc_
edm::EDGetTokenT< SiStripMatchedRecHit2DCollection >	SiStripSrc_Token_
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tGeomToken_
std::map< int, edm::PSimHitContainer >	theMap
TrackerHitAssociator::Config	trackerHitAssociatorConfig_
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_
int	verbosity

Additional Inherited Members
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 143 of file GlobalRecHitsAnalyzer.h.

Member Typedef Documentation

MapType

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

Definition at line 145 of file GlobalRecHitsAnalyzer.h.

Constructor & Destructor Documentation

GlobalRecHitsAnalyzer()

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

explicit

Definition at line 18 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(), HBHERecHitgetter_, HCalSrc_, HCalSrc_Token_, HFRecHitgetter_, hitsProducer, HORecHitgetter_, ProducerED_cfi::InputTag, edm::InputTag::instance(), edm::BranchDescription::isAnyAlias(), 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()),
      caloGeomToken_(esConsumes()),
      tTopoToken_(esConsumes()),
      tGeomToken_(esConsumes()),
      dtGeomToken_(esConsumes()),
      cscGeomToken_(esConsumes()),
      rpcGeomToken_(esConsumes()),
      count(0) {
  HBHERecHitgetter_ = edm::GetterOfProducts<edm::SortedCollection<HBHERecHit, edm::StrictWeakOrdering<HBHERecHit>>>(
      edm::ProcessMatch("*"), this);
  HFRecHitgetter_ = edm::GetterOfProducts<edm::SortedCollection<HFRecHit, edm::StrictWeakOrdering<HFRecHit>>>(
      edm::ProcessMatch("*"), this);
  HORecHitgetter_ = edm::GetterOfProducts<edm::SortedCollection<HORecHit, edm::StrictWeakOrdering<HORecHit>>>(
      edm::ProcessMatch("*"), this);
  callWhenNewProductsRegistered([this](edm::BranchDescription const& bd) {
    // in case of EDAliases, consume only the aliased-for original products
    if (bd.isAnyAlias())
      return;
    this->HBHERecHitgetter_(bd);
    this->HFRecHitgetter_(bd);
    this->HORecHitgetter_(bd);
  });
  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::~GlobalRecHitsAnalyzer ( )

override

Definition at line 139 of file GlobalRecHitsAnalyzer.cc.

{}

Member Function Documentation

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 311 of file GlobalRecHitsAnalyzer.cc.

References count, fillECal(), fillHCal(), fillMuon(), fillTrk(), frequency, getAllProvenances, mps_fire::i, iEvent, nevt, printProvenanceInfo, 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;
}

bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 141 of file GlobalRecHitsAnalyzer.cc.

References dqm::implementation::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::implementation::NavigatorBase::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.);
}

compute()

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 1355 of file GlobalRecHitsAnalyzer.cc.

References dqm::impl::MonitorElement::Fill(), findBestRecHit(), DTHitQualityUtils::findMuSimHit(), nano_mu_digi_cff::layer, 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;
}

fillECal()

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

private

Definition at line 377 of file GlobalRecHitsAnalyzer.cc.

References HLT_2023v12_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_2023v12_cff::endcapHits, ESHits_Token_, dqm::impl::MonitorElement::Fill(), iEvent, 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 = EcalUncalibRecHitEB.isValid();
  if (!validUncalibRecHitEB) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalUncalRecHitEB in event!";
  }

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

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

  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 = EcalUncalibRecHitEE.isValid();
  if (!validuncalibRecHitEE) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalUncalRecHitEE in event!";
  }

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

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

  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 = EcalRecHitES.isValid();
  if (!validRecHitES) {
    LogDebug(MsgLoggerCat) << "Unable to find EcalRecHitES in event!";
  }

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

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

fillHCal()

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

private

Definition at line 564 of file GlobalRecHitsAnalyzer.cc.

References caloGeomToken_, hcalRecHitTable_cff::detId, dqm::impl::MonitorElement::Fill(), edm::GetterOfProducts< T >::fillHandles(), edm::EventSetup::getHandle(), photonIsolationHIProducer_cfi::hbhe, HBHERecHitgetter_, HCalSrc_Token_, photonIsolationHIProducer_cfi::hf, HFRecHitgetter_, photonIsolationHIProducer_cfi::ho, HORecHitgetter_, iEvent, edm::HandleBase::isValid(), LogDebug, mehHcaln, mehHcalRes, edm::Handle< 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
  const auto& geometry = iSetup.getHandle(caloGeomToken_);
  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 = hcalHits.isValid();
  if (!validhcalHits) {
    LogDebug(MsgLoggerCat) << "Unable to find hcalHits in event!";
  }

  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();
      }
    }
  }

  // get HBHE information
  std::vector<edm::Handle<HBHERecHitCollection>> hbhe;
  HBHERecHitgetter_.fillHandles(iEvent, hbhe);
  bool validHBHE = hbhe[0].isValid();

  if (!validHBHE) {
    LogDebug(MsgLoggerCat) << "Unable to find any HBHERecHitCollections in event!";
  }

  else {
    std::vector<edm::Handle<HBHERecHitCollection>>::iterator ihbhe;
    int iHB = 0;
    int iHE = 0;
    for (ihbhe = hbhe.begin(); ihbhe != hbhe.end(); ++ihbhe) {
      for (HBHERecHitCollection::const_iterator jhbhe = (*ihbhe)->begin(); jhbhe != (*ihbhe)->end(); ++jhbhe) {
        HcalDetId cell(jhbhe->id());

        if (cell.subdet() == sdHcalBrl) {
          ++iHB;
          mehHcalRes[0]->Fill(jhbhe->energy() - fHBEnergySimHits[cell]);
        }

        else if (cell.subdet() == sdHcalEC) {
          ++iHE;
          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;
  HFRecHitgetter_.fillHandles(iEvent, hf);
  bool validHF = hf[0].isValid();
  if (!validHF) {
    LogDebug(MsgLoggerCat) << "Unable to find any HFRecHitCollections in event!";
  } else {
    std::vector<edm::Handle<HFRecHitCollection>>::iterator ihf;

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

        if (cell.subdet() == sdHcalFwd) {
          ++iHF;
          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;
  HORecHitgetter_.fillHandles(iEvent, ho);
  bool validHO = ho[0].isValid();
  if (!validHO) {
    LogDebug(MsgLoggerCat) << "Unable to find any HORecHitCollections in event!";
  } else {
    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;
          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;
}

fillMuon()

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

private

Definition at line 1072 of file GlobalRecHitsAnalyzer.cc.

References nano_mu_digi_cff::chamberType, compute(), cscGeomToken_, hcalRecHitTable_cff::detId, dtGeomToken_, dqm::impl::MonitorElement::Fill(), edm::EventSetup::getHandle(), mps_fire::i, CSCGeometry::idToDetUnit(), iEvent, edm::HandleBase::isValid(), nano_mu_digi_cff::layer, 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(), OfflineHarvestingSequence_cosmic::ptype, alignCSCRings::r, rpcPointValidation_cfi::recHit, nano_mu_digi_cff::roll, rpcGeomToken_, rpcPointValidation_cfi::simHit, FastTrackerRecHitCombiner_cfi::simHits, 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
  const auto& dtGeom = iSetup.getHandle(dtGeomToken_);
  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 = dtsimHits.isValid();
  if (!validdtsim) {
    LogDebug(MsgLoggerCat) << "Unable to find dtsimHits in event!";
  }

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

  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 = cf.isValid();
  if (!validXFrame) {
    LogDebug(MsgLoggerCat) << "Unable to find muo CSC crossingFrame in event!";
  } else {
    const MixCollection<PSimHit> simHits(cf.product());

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

  // get geometry
  const auto& hGeom = iSetup.getHandle(cscGeomToken_);
  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 = hRecHits.isValid();
  if (!validCSC) {
    LogDebug(MsgLoggerCat) << "Unable to find CSC RecHits in event!";
  } else {
    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;
  const auto& rpcGeom = iSetup.getHandle(rpcGeomToken_);
  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 = simHit.isValid();
  if (!validrpcsim) {
    LogDebug(MsgLoggerCat) << "Unable to find RPCSimHit in event!";
  }

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

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

fillTrk()

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

private

Definition at line 724 of file GlobalRecHitsAnalyzer.cc.

References TrackerHitAssociator::associateHit(), edmNew::DetSet< T >::begin(), pv::closest(), hcalRecHitTable_cff::detId, edmNew::DetSet< T >::end(), edmNew::DetSetVector< T >::end(), dqm::impl::MonitorElement::Fill(), edmNew::DetSetVector< T >::find(), TrackingRecHit::geographicalId(), edm::EventSetup::getData(), edm::EventSetup::getHandle(), mps_fire::i, iEvent, edm::HandleBase::isValid(), BaseTrackerRecHit::localPosition(), LogDebug, visualization-live-secondInstance_cfg::m, matched, mehSiPixeln, mehSiPixelResX, mehSiPixelResY, mehSiStripn, mehSiStripResX, mehSiStripResY, position, projectHit(), TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), TrackerTopology::pxfSide(), nano_mu_digi_cff::rawId, sdPxlBrl, sdPxlFwd, sdSiTEC, sdSiTIB, sdSiTID, sdSiTOB, SiPxlSrc_Token_, SiStripSrc_Token_, mathSSE::sqrt(), GluedGeomDet::stereoDet(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), GeomDet::surface(), TrackerTopology::tecWheel(), tGeomToken_, TrackerTopology::tibLayer(), TrackerTopology::tidWheel(), TrackerTopology::tobLayer(), PbPb_ZMuSkimMuonDPG_cff::tracker, trackerHitAssociatorConfig_, tTopoToken_, verbosity, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by analyze().

                                                                                       {
  //Retrieve tracker topology from geometry
  const TrackerTopology* const tTopo = &iSetup.getData(tTopoToken_);
  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 = rechitsmatched.isValid();
  if (!validstrip) {
    LogDebug(MsgLoggerCat) << "Unable to find stripmatchedrechits in event!";
  }

  TrackerHitAssociator associate(iEvent, trackerHitAssociatorConfig_);

  const auto& tGeomHandle = iSetup.getHandle(tGeomToken_);
  if (!tGeomHandle.isValid()) {
    edm::LogWarning(MsgLoggerCat) << "Unable to find TrackerDigiGeometry in event!";
    return;
  }
  const TrackerGeometry& tracker(*tGeomHandle);

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

    // loop over det units
    for (TrackerGeometry::DetContainer::const_iterator it = tGeomHandle->dets().begin();
         it != tGeomHandle->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;

        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 = recHitColl.isValid();
  if (!validpixel) {
    LogDebug(MsgLoggerCat) << "Unable to find SiPixelRecHitCollection in event!";
  } else {
    int nPxlBrl = 0, nPxlFwd = 0;
    //iterate over detunits
    for (TrackerGeometry::DetContainer::const_iterator it = tGeomHandle->dets().begin();
         it != tGeomHandle->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;
}

findBestRecHit()

template<typename type >

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

private

Definition at line 1325 of file GlobalRecHitsAnalyzer.cc.

References nano_mu_digi_cff::layer, rpcPointValidation_cfi::recHit, recHitDistFromWire(), and FastTrackerRecHitMaskProducer_cfi::recHits.

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

map1DRecHitsPerWire()

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

private

Definition at line 1299 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;
}

plotResolution()

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

private

Definition at line 1406 of file GlobalRecHitsAnalyzer.cc.

References dqm::impl::MonitorElement::Fill(), nano_mu_digi_cff::layer, mehCSCResRDPhi, PV3DBase< T, PVType, FrameType >::phi(), rpcPointValidation_cfi::recHit, and rpcPointValidation_cfi::simHit.

Referenced by fillMuon().

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

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

projectHit()

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

private

Definition at line 1272 of file GlobalRecHitsAnalyzer.cc.

References funct::cos(), DeadROC_duringRun::dir, pfClustersFromCombinedCaloHF_cfi::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);
}

recHitDistFromWire() [1/2]

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

private

Definition at line 1344 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.;
}

recHitDistFromWire() [2/2]

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

private

Definition at line 1350 of file GlobalRecHitsAnalyzer.cc.

References nano_mu_digi_cff::layer, and rpcPointValidation_cfi::recHit.

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

simHitDistFromWire()

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

private

Definition at line 1312 of file GlobalRecHitsAnalyzer.cc.

References nano_mu_digi_cff::layer, DTWireId::wire(), 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

caloGeomToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> GlobalRecHitsAnalyzer::caloGeomToken_

private

Definition at line 293 of file GlobalRecHitsAnalyzer.h.

Referenced by fillHCal().

count

unsigned int GlobalRecHitsAnalyzer::count

private

Definition at line 301 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze().

cscGeomToken_

edm::ESGetToken<CSCGeometry, MuonGeometryRecord> GlobalRecHitsAnalyzer::cscGeomToken_

private

Definition at line 297 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon().

dtGeomToken_

edm::ESGetToken<DTGeometry, MuonGeometryRecord> GlobalRecHitsAnalyzer::dtGeomToken_

private

Definition at line 296 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon().

EBHits_Token_

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

private

Definition at line 197 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

ECalEBSrc_

edm::InputTag GlobalRecHitsAnalyzer::ECalEBSrc_

private

Definition at line 187 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

ECalEBSrc_Token_

edm::EDGetTokenT<EBRecHitCollection> GlobalRecHitsAnalyzer::ECalEBSrc_Token_

private

Definition at line 192 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

ECalEESrc_

edm::InputTag GlobalRecHitsAnalyzer::ECalEESrc_

private

Definition at line 189 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

ECalEESrc_Token_

edm::EDGetTokenT<EERecHitCollection> GlobalRecHitsAnalyzer::ECalEESrc_Token_

private

Definition at line 193 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

ECalESSrc_

edm::InputTag GlobalRecHitsAnalyzer::ECalESSrc_

private

Definition at line 191 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

ECalESSrc_Token_

edm::EDGetTokenT<ESRecHitCollection> GlobalRecHitsAnalyzer::ECalESSrc_Token_

private

Definition at line 194 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

ECalUncalEBSrc_

edm::InputTag GlobalRecHitsAnalyzer::ECalUncalEBSrc_

private

Definition at line 188 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

ECalUncalEBSrc_Token_

edm::EDGetTokenT<EBUncalibratedRecHitCollection> GlobalRecHitsAnalyzer::ECalUncalEBSrc_Token_

private

Definition at line 195 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

ECalUncalEESrc_

edm::InputTag GlobalRecHitsAnalyzer::ECalUncalEESrc_

private

Definition at line 190 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

ECalUncalEESrc_Token_

edm::EDGetTokenT<EEUncalibratedRecHitCollection> GlobalRecHitsAnalyzer::ECalUncalEESrc_Token_

private

Definition at line 196 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

EEHits_Token_

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

private

Definition at line 198 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

ESHits_Token_

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

private

Definition at line 199 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

fName

std::string GlobalRecHitsAnalyzer::fName

private

Definition at line 169 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

frequency

int GlobalRecHitsAnalyzer::frequency

private

Definition at line 171 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze(), and GlobalRecHitsAnalyzer().

getAllProvenances

bool GlobalRecHitsAnalyzer::getAllProvenances

private

Definition at line 173 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze(), and GlobalRecHitsAnalyzer().

HBHERecHitgetter_

edm::GetterOfProducts<edm::SortedCollection<HBHERecHit, edm::StrictWeakOrdering<HBHERecHit> > > GlobalRecHitsAnalyzer::HBHERecHitgetter_

private

Definition at line 183 of file GlobalRecHitsAnalyzer.h.

Referenced by fillHCal(), and GlobalRecHitsAnalyzer().

HCalSrc_

edm::InputTag GlobalRecHitsAnalyzer::HCalSrc_

private

Definition at line 206 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

HCalSrc_Token_

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

private

Definition at line 207 of file GlobalRecHitsAnalyzer.h.

Referenced by fillHCal(), and GlobalRecHitsAnalyzer().

HFRecHitgetter_

edm::GetterOfProducts<edm::SortedCollection<HFRecHit, edm::StrictWeakOrdering<HFRecHit> > > GlobalRecHitsAnalyzer::HFRecHitgetter_

private

Definition at line 184 of file GlobalRecHitsAnalyzer.h.

Referenced by fillHCal(), and GlobalRecHitsAnalyzer().

hitsProducer

std::string GlobalRecHitsAnalyzer::hitsProducer

private

Definition at line 175 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

HORecHitgetter_

edm::GetterOfProducts<edm::SortedCollection<HORecHit, edm::StrictWeakOrdering<HORecHit> > > GlobalRecHitsAnalyzer::HORecHitgetter_

private

Definition at line 185 of file GlobalRecHitsAnalyzer.h.

Referenced by fillHCal(), and GlobalRecHitsAnalyzer().

label

std::string GlobalRecHitsAnalyzer::label

private

Definition at line 172 of file GlobalRecHitsAnalyzer.h.

matched

std::vector<PSimHit> GlobalRecHitsAnalyzer::matched

private

Definition at line 219 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk().

mehCSCn

MonitorElement* GlobalRecHitsAnalyzer::mehCSCn

private

Definition at line 238 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

mehCSCResRDPhi

MonitorElement* GlobalRecHitsAnalyzer::mehCSCResRDPhi

private

Definition at line 241 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and plotResolution().

mehDtMuonn

MonitorElement* GlobalRecHitsAnalyzer::mehDtMuonn

private

Definition at line 237 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

mehDtMuonRes

MonitorElement* GlobalRecHitsAnalyzer::mehDtMuonRes

private

Definition at line 240 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and compute().

mehEcaln

MonitorElement* GlobalRecHitsAnalyzer::mehEcaln[3]

private

Definition at line 180 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

mehEcalRes

MonitorElement* GlobalRecHitsAnalyzer::mehEcalRes[3]

private

Definition at line 181 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

mehHcaln

MonitorElement* GlobalRecHitsAnalyzer::mehHcaln[4]

private

Definition at line 203 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

mehHcalRes

MonitorElement* GlobalRecHitsAnalyzer::mehHcalRes[4]

private

Definition at line 204 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

mehRPCn

MonitorElement* GlobalRecHitsAnalyzer::mehRPCn

private

Definition at line 239 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

mehRPCResX

MonitorElement* GlobalRecHitsAnalyzer::mehRPCResX

private

Definition at line 242 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

mehSiPixeln

MonitorElement* GlobalRecHitsAnalyzer::mehSiPixeln[7]

private

Definition at line 227 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

mehSiPixelResX

MonitorElement* GlobalRecHitsAnalyzer::mehSiPixelResX[7]

private

Definition at line 228 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

mehSiPixelResY

MonitorElement* GlobalRecHitsAnalyzer::mehSiPixelResY[7]

private

Definition at line 229 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

mehSiStripn

MonitorElement* GlobalRecHitsAnalyzer::mehSiStripn[19]

private

Definition at line 212 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

mehSiStripResX

MonitorElement* GlobalRecHitsAnalyzer::mehSiStripResX[19]

private

Definition at line 213 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

mehSiStripResY

MonitorElement* GlobalRecHitsAnalyzer::mehSiStripResY[19]

private

Definition at line 214 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MuCSCHits_Token_

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

private

Definition at line 279 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

MuCSCSrc_

edm::InputTag GlobalRecHitsAnalyzer::MuCSCSrc_

private

Definition at line 277 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

MuCSCSrc_Token_

edm::EDGetTokenT<CSCRecHit2DCollection> GlobalRecHitsAnalyzer::MuCSCSrc_Token_

private

Definition at line 278 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

MuDTSimSrc_

edm::InputTag GlobalRecHitsAnalyzer::MuDTSimSrc_

private

Definition at line 245 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

MuDTSimSrc_Token_

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

private

Definition at line 247 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

MuDTSrc_

edm::InputTag GlobalRecHitsAnalyzer::MuDTSrc_

private

Definition at line 244 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

MuDTSrc_Token_

edm::EDGetTokenT<DTRecHitCollection> GlobalRecHitsAnalyzer::MuDTSrc_Token_

private

Definition at line 246 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

MuRPCSimSrc_

edm::InputTag GlobalRecHitsAnalyzer::MuRPCSimSrc_

private

Definition at line 289 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

MuRPCSimSrc_Token_

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

private

Definition at line 291 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

MuRPCSrc_

edm::InputTag GlobalRecHitsAnalyzer::MuRPCSrc_

private

Definition at line 288 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

MuRPCSrc_Token_

edm::EDGetTokenT<RPCRecHitCollection> GlobalRecHitsAnalyzer::MuRPCSrc_Token_

private

Definition at line 290 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

printProvenanceInfo

bool GlobalRecHitsAnalyzer::printProvenanceInfo

private

Definition at line 174 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze(), and GlobalRecHitsAnalyzer().

rpcGeomToken_

edm::ESGetToken<RPCGeometry, MuonGeometryRecord> GlobalRecHitsAnalyzer::rpcGeomToken_

private

Definition at line 298 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon().

SiPxlSrc_

edm::InputTag GlobalRecHitsAnalyzer::SiPxlSrc_

private

Definition at line 231 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

SiPxlSrc_Token_

edm::EDGetTokenT<SiPixelRecHitCollection> GlobalRecHitsAnalyzer::SiPxlSrc_Token_

private

Definition at line 232 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalRecHitsAnalyzer().

SiStripSrc_

edm::InputTag GlobalRecHitsAnalyzer::SiStripSrc_

private

Definition at line 216 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

SiStripSrc_Token_

edm::EDGetTokenT<SiStripMatchedRecHit2DCollection> GlobalRecHitsAnalyzer::SiStripSrc_Token_

private

Definition at line 217 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalRecHitsAnalyzer().

tGeomToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> GlobalRecHitsAnalyzer::tGeomToken_

private

Definition at line 295 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk().

theMap

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

private

Definition at line 281 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon().

trackerHitAssociatorConfig_

TrackerHitAssociator::Config GlobalRecHitsAnalyzer::trackerHitAssociatorConfig_

private

Definition at line 223 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk().

tTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> GlobalRecHitsAnalyzer::tTopoToken_

private

Definition at line 294 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk().

verbosity

int GlobalRecHitsAnalyzer::verbosity

private

Definition at line 170 of file GlobalRecHitsAnalyzer.h.

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