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::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::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 141 of file GlobalRecHitsAnalyzer.h.

Member Typedef Documentation

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

Definition at line 143 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_FULL_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()),
      caloGeomToken_(esConsumes()),
      tTopoToken_(esConsumes()),
      tGeomToken_(esConsumes()),
      dtGeomToken_(esConsumes()),
      cscGeomToken_(esConsumes()),
      rpcGeomToken_(esConsumes()),
      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 126 of file GlobalRecHitsAnalyzer.cc.

{}

Member Function Documentation

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 298 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 128 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.);
}

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

References dqm::impl::MonitorElement::Fill(), findBestRecHit(), DTHitQualityUtils::findMuSimHit(), DTGeometry::layer(), phase1PixelTopology::layer, mehDtMuonRes, recHitDistFromWire(), HLT_FULL_cff::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 364 of file GlobalRecHitsAnalyzer.cc.

References HLT_FULL_cff::barrelHits, edm::SortedCollection< T, SORT >::begin(), EBHits_Token_, ECalEBSrc_Token_, ECalEESrc_Token_, ECalESSrc_Token_, ECalUncalEBSrc_Token_, ECalUncalEESrc_Token_, EEHits_Token_, edm::SortedCollection< T, SORT >::end(), HLT_FULL_cff::endcapHits, ESHits_Token_, dqm::impl::MonitorElement::Fill(), edm::SortedCollection< T, SORT >::find(), edm::Event::getByToken(), edm::HandleBase::isValid(), LogDebug, mehEcaln, mehEcalRes, ootPhotons_cff::preshowerHits, edm::Handle< T >::product(), DetId::rawId(), 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;
}

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

private

Definition at line 551 of file GlobalRecHitsAnalyzer.cc.

References caloGeomToken_, dqm::impl::MonitorElement::Fill(), relativeConstraints::geometry, edm::Event::getByToken(), edm::EventSetup::getHandle(), edm::Event::getManyByType(), HCalSrc_Token_, hcalSimParameters_cfi::ho, ihf, edm::HandleBase::isValid(), LogDebug, mehHcaln, mehHcalRes, edm::Handle< T >::product(), sdHcalBrl, sdHcalEC, sdHcalFwd, sdHcalOut, 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;
  iEvent.getManyByType(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;
  iEvent.getManyByType(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;
  iEvent.getManyByType(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;
}

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

private

Definition at line 1060 of file GlobalRecHitsAnalyzer.cc.

References CSCLayer::chamber(), CSCChamberSpecs::chamberType(), compute(), CSCDetId, cscGeomToken_, dtGeomToken_, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), edm::EventSetup::getHandle(), mps_fire::i, CSCGeometry::idToDetUnit(), RPCRoll::isForward(), edm::HandleBase::isValid(), phase1PixelTopology::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(), alignCSCRings::r, RPCDetId, rpcGeomToken_, trackerHits::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
  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;
}

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

private

Definition at line 712 of file GlobalRecHitsAnalyzer.cc.

References TrackerHitAssociator::associateHit(), edmNew::DetSet< T >::begin(), edmNew::DetSet< T >::end(), dqm::impl::MonitorElement::Fill(), TrackingRecHit::geographicalId(), edm::Event::getByToken(), edm::EventSetup::getData(), edm::EventSetup::getHandle(), mps_fire::i, TrackerGeometry::idToDet(), 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(), 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(), 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;
}

template<typename type >

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

private

Definition at line 1313 of file GlobalRecHitsAnalyzer.cc.

References 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 1287 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 1394 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 1260 of file GlobalRecHitsAnalyzer.cc.

References funct::cos(), DeadROC_duringRun::dir, PSimHit::localDirection(), PSimHit::localPosition(), pileupReCalc_HLTpaths::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 1332 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 1338 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 1300 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

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

private

Definition at line 287 of file GlobalRecHitsAnalyzer.h.

Referenced by fillHCal().

unsigned int GlobalRecHitsAnalyzer::count

private

Definition at line 295 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 291 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon().

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

private

Definition at line 290 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon().

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

private

Definition at line 191 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::ECalEBSrc_

private

Definition at line 181 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<EBRecHitCollection> GlobalRecHitsAnalyzer::ECalEBSrc_Token_

private

Definition at line 186 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::ECalEESrc_

private

Definition at line 183 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<EERecHitCollection> GlobalRecHitsAnalyzer::ECalEESrc_Token_

private

Definition at line 187 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::ECalESSrc_

private

Definition at line 185 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<ESRecHitCollection> GlobalRecHitsAnalyzer::ECalESSrc_Token_

private

Definition at line 188 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::ECalUncalEBSrc_

private

Definition at line 182 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<EBUncalibratedRecHitCollection> GlobalRecHitsAnalyzer::ECalUncalEBSrc_Token_

private

Definition at line 189 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::ECalUncalEESrc_

private

Definition at line 184 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<EEUncalibratedRecHitCollection> GlobalRecHitsAnalyzer::ECalUncalEESrc_Token_

private

Definition at line 190 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

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

private

Definition at line 192 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

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

private

Definition at line 193 of file GlobalRecHitsAnalyzer.h.

Referenced by fillECal(), and GlobalRecHitsAnalyzer().

std::string GlobalRecHitsAnalyzer::fName

private

Definition at line 167 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

int GlobalRecHitsAnalyzer::frequency

private

Definition at line 169 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze(), and GlobalRecHitsAnalyzer().

bool GlobalRecHitsAnalyzer::getAllProvenances

private

Definition at line 171 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::HCalSrc_

private

Definition at line 200 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

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

private

Definition at line 201 of file GlobalRecHitsAnalyzer.h.

Referenced by fillHCal(), and GlobalRecHitsAnalyzer().

std::string GlobalRecHitsAnalyzer::hitsProducer

private

Definition at line 173 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

std::string GlobalRecHitsAnalyzer::label

private

Definition at line 170 of file GlobalRecHitsAnalyzer.h.

std::vector<PSimHit> GlobalRecHitsAnalyzer::matched

private

Definition at line 213 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehCSCn

private

Definition at line 232 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalRecHitsAnalyzer::mehCSCResRDPhi

private

Definition at line 235 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and plotResolution().

MonitorElement* GlobalRecHitsAnalyzer::mehDtMuonn

private

Definition at line 231 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalRecHitsAnalyzer::mehDtMuonRes

private

Definition at line 234 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and compute().

MonitorElement* GlobalRecHitsAnalyzer::mehEcaln[3]

private

Definition at line 178 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalRecHitsAnalyzer::mehEcalRes[3]

private

Definition at line 179 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillECal().

MonitorElement* GlobalRecHitsAnalyzer::mehHcaln[4]

private

Definition at line 197 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

MonitorElement* GlobalRecHitsAnalyzer::mehHcalRes[4]

private

Definition at line 198 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillHCal().

MonitorElement* GlobalRecHitsAnalyzer::mehRPCn

private

Definition at line 233 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalRecHitsAnalyzer::mehRPCResX

private

Definition at line 236 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillMuon().

MonitorElement* GlobalRecHitsAnalyzer::mehSiPixeln[7]

private

Definition at line 221 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehSiPixelResX[7]

private

Definition at line 222 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehSiPixelResY[7]

private

Definition at line 223 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehSiStripn[19]

private

Definition at line 206 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehSiStripResX[19]

private

Definition at line 207 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

MonitorElement* GlobalRecHitsAnalyzer::mehSiStripResY[19]

private

Definition at line 208 of file GlobalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and fillTrk().

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

private

Definition at line 273 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::MuCSCSrc_

private

Definition at line 271 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<CSCRecHit2DCollection> GlobalRecHitsAnalyzer::MuCSCSrc_Token_

private

Definition at line 272 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::MuDTSimSrc_

private

Definition at line 239 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

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

private

Definition at line 241 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::MuDTSrc_

private

Definition at line 238 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<DTRecHitCollection> GlobalRecHitsAnalyzer::MuDTSrc_Token_

private

Definition at line 240 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::MuRPCSimSrc_

private

Definition at line 283 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

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

private

Definition at line 285 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::MuRPCSrc_

private

Definition at line 282 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<RPCRecHitCollection> GlobalRecHitsAnalyzer::MuRPCSrc_Token_

private

Definition at line 284 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon(), and GlobalRecHitsAnalyzer().

bool GlobalRecHitsAnalyzer::printProvenanceInfo

private

Definition at line 172 of file GlobalRecHitsAnalyzer.h.

Referenced by analyze(), and GlobalRecHitsAnalyzer().

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

private

Definition at line 292 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon().

edm::InputTag GlobalRecHitsAnalyzer::SiPxlSrc_

private

Definition at line 225 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<SiPixelRecHitCollection> GlobalRecHitsAnalyzer::SiPxlSrc_Token_

private

Definition at line 226 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalRecHitsAnalyzer().

edm::InputTag GlobalRecHitsAnalyzer::SiStripSrc_

private

Definition at line 210 of file GlobalRecHitsAnalyzer.h.

Referenced by GlobalRecHitsAnalyzer().

edm::EDGetTokenT<SiStripMatchedRecHit2DCollection> GlobalRecHitsAnalyzer::SiStripSrc_Token_

private

Definition at line 211 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk(), and GlobalRecHitsAnalyzer().

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

private

Definition at line 289 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk().

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

private

Definition at line 275 of file GlobalRecHitsAnalyzer.h.

Referenced by fillMuon().

TrackerHitAssociator::Config GlobalRecHitsAnalyzer::trackerHitAssociatorConfig_

private

Definition at line 217 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk().

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

private

Definition at line 288 of file GlobalRecHitsAnalyzer.h.

Referenced by fillTrk().

int GlobalRecHitsAnalyzer::verbosity

private

Definition at line 168 of file GlobalRecHitsAnalyzer.h.

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