CMS 3D CMS Logo
List of all members | Classes | Public Types | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
Phase2OTValidateCluster Class Reference

#include <Phase2OTValidateCluster.cc>

Inheritance diagram for Phase2OTValidateCluster:
DQMEDAnalyzer edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >

Classes

struct  ClusterMEs
 

Public Types

typedef std::map< unsigned int, std::vector< PSimHit > > SimHitsMap
 
typedef std::map< unsigned int, SimTrack > SimTracksMap
 
- 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 &iEvent, const edm::EventSetup &iSetup) override
 
void bookHistograms (DQMStore::IBooker &ibooker, edm::Run const &iRun, edm::EventSetup const &iSetup) override
 
void dqmBeginRun (const edm::Run &iRun, const edm::EventSetup &iSetup) override
 
 Phase2OTValidateCluster (const edm::ParameterSet &)
 
 ~Phase2OTValidateCluster () 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
 
 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
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- 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 &)
 

Private Member Functions

void bookLayerHistos (DQMStore::IBooker &ibooker, uint32_t det_it, const std::string &subdir)
 
void fillOTHistos (const edm::Event &iEvent, const std::vector< edm::Handle< edm::PSimHitContainer >> &simHits, const std::map< unsigned int, SimTrack > &simTracks)
 
std::vector< unsigned int > getSimTrackId (const edm::Handle< edm::DetSetVector< PixelDigiSimLink >> &pixelSimLinks, const DetId &detId, unsigned int channel)
 

Private Attributes

edm::EDGetTokenT< Phase2TrackerCluster1DCollectionNew > clustersToken_
 
edm::ParameterSet config_
 
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
 
std::map< std::string, ClusterMEs > layerMEs_
 
std::vector< edm::InputTag > pSimHitSrc_
 
std::vector< edm::EDGetTokenT< edm::PSimHitContainer > > simHitTokens_
 
edm::EDGetTokenT< edm::DetSetVector< PixelDigiSimLink > > simITLinksToken_
 
edm::EDGetTokenT< edm::DetSetVector< PixelDigiSimLink > > simOTLinksToken_
 
double simtrackminpt_
 
edm::EDGetTokenT< edm::SimTrackContainer > simTracksToken_
 
const TrackerGeometry * tkGeom_ = nullptr
 
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
 
const TrackerTopology * tTopo_ = nullptr
 

Additional Inherited Members

- Protected Member Functions inherited from DQMEDAnalyzer
uint64_t meId () const
 
- Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken > lumiToken_
 
edm::EDPutTokenT< DQMToken > runToken_
 
unsigned int streamId_
 

Detailed Description

Description: Validation plots tracker clusters.

Definition at line 50 of file Phase2OTValidateCluster.cc.

Member Typedef Documentation

◆ SimHitsMap

typedef std::map<unsigned int, std::vector<PSimHit> > Phase2OTValidateCluster::SimHitsMap

Definition at line 52 of file Phase2OTValidateCluster.cc.

◆ SimTracksMap

typedef std::map<unsigned int, SimTrack> Phase2OTValidateCluster::SimTracksMap

Definition at line 53 of file Phase2OTValidateCluster.cc.

Constructor & Destructor Documentation

◆ Phase2OTValidateCluster()

Phase2OTValidateCluster::Phase2OTValidateCluster ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 102 of file Phase2OTValidateCluster.cc.

References pSimHitSrc_, and simHitTokens_.

103  : config_(iConfig),
104  simtrackminpt_(config_.getParameter<double>("SimTrackMinPt")),
105  simOTLinksToken_(consumes<edm::DetSetVector<PixelDigiSimLink>>(
106  config_.getParameter<edm::InputTag>("OuterTrackerDigiSimLinkSource"))),
107  simTracksToken_(consumes<edm::SimTrackContainer>(config_.getParameter<edm::InputTag>("simtracks"))),
108  clustersToken_(
109  consumes<Phase2TrackerCluster1DCollectionNew>(config_.getParameter<edm::InputTag>("ClusterSource"))),
110  pSimHitSrc_(config_.getParameter<std::vector<edm::InputTag>>("PSimHitSource")),
111  geomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>()),
112  topoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>()) {
113  edm::LogInfo("Phase2OTValidateCluster") << ">>> Construct Phase2OTValidateCluster ";
114  for (const auto& itag : pSimHitSrc_)
115  simHitTokens_.push_back(consumes<edm::PSimHitContainer>(itag));
116 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
Phase2OTValidateCluster::simTracksToken_
edm::EDGetTokenT< edm::SimTrackContainer > simTracksToken_
Definition: Phase2OTValidateCluster.cc:89
Phase2OTValidateCluster::pSimHitSrc_
std::vector< edm::InputTag > pSimHitSrc_
Definition: Phase2OTValidateCluster.cc:91
Phase2OTValidateCluster::simOTLinksToken_
edm::EDGetTokenT< edm::DetSetVector< PixelDigiSimLink > > simOTLinksToken_
Definition: Phase2OTValidateCluster.cc:87
Phase2OTValidateCluster::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Definition: Phase2OTValidateCluster.cc:93
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
Phase2OTValidateCluster::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: Phase2OTValidateCluster.cc:92
Phase2OTValidateCluster::simHitTokens_
std::vector< edm::EDGetTokenT< edm::PSimHitContainer > > simHitTokens_
Definition: Phase2OTValidateCluster.cc:86
Phase2OTValidateCluster::clustersToken_
edm::EDGetTokenT< Phase2TrackerCluster1DCollectionNew > clustersToken_
Definition: Phase2OTValidateCluster.cc:90

◆ ~Phase2OTValidateCluster()

Phase2OTValidateCluster::~Phase2OTValidateCluster ( )
override

Definition at line 118 of file Phase2OTValidateCluster.cc.

118  {
119  // do anything here that needs to be done at desctruction time
120  // (e.g. close files, deallocate resources etc.)
121  edm::LogInfo("Phase2OTValidateCluster") << ">>> Destroy Phase2OTValidateCluster ";
122 }
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125

Member Function Documentation

◆ analyze()

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

Reimplemented from DQMEDAnalyzer.

Definition at line 133 of file Phase2OTValidateCluster.cc.

References fillOTHistos(), iEvent, FastTrackerRecHitCombiner_cfi::simHits, simHitTokens_, simtrackminpt_, TrackCandidateProducer_cfi::simTracks, and simTracksToken_.

133  {
134  // Getting simHits
135  std::vector<edm::Handle<edm::PSimHitContainer>> simHits;
136  for (const auto& itoken : simHitTokens_) {
137  const auto& simHitHandle = iEvent.getHandle(itoken);
138  if (!simHitHandle.isValid())
139  continue;
140  simHits.emplace_back(simHitHandle);
141  }
142  // Get the SimTracks
143  const auto& simTracksRaw = iEvent.getHandle(simTracksToken_);
144  // Rearrange the simTracks for ease of use <simTrackID, simTrack>
145  SimTracksMap simTracks;
146  for (edm::SimTrackContainer::const_iterator simTrackIt(simTracksRaw->begin()); simTrackIt != simTracksRaw->end();
147  ++simTrackIt) {
148  if (simTrackIt->momentum().pt() > simtrackminpt_) {
149  simTracks.emplace(simTrackIt->trackId(), *simTrackIt);
150  }
151  }
152  fillOTHistos(iEvent, simHits, simTracks);
153 }
Phase2OTValidateCluster::SimTracksMap
std::map< unsigned int, SimTrack > SimTracksMap
Definition: Phase2OTValidateCluster.cc:53
Phase2OTValidateCluster::simTracksToken_
edm::EDGetTokenT< edm::SimTrackContainer > simTracksToken_
Definition: Phase2OTValidateCluster.cc:89
Phase2OTValidateCluster::fillOTHistos
void fillOTHistos(const edm::Event &iEvent, const std::vector< edm::Handle< edm::PSimHitContainer >> &simHits, const std::map< unsigned int, SimTrack > &simTracks)
Definition: Phase2OTValidateCluster.cc:155
iEvent
int iEvent
Definition: GenABIO.cc:224
TrackCandidateProducer_cfi.simTracks
simTracks
Definition: TrackCandidateProducer_cfi.py:15
Phase2OTValidateCluster::simHitTokens_
std::vector< edm::EDGetTokenT< edm::PSimHitContainer > > simHitTokens_
Definition: Phase2OTValidateCluster.cc:86
FastTrackerRecHitCombiner_cfi.simHits
simHits
Definition: FastTrackerRecHitCombiner_cfi.py:5

◆ bookHistograms()

void Phase2OTValidateCluster::bookHistograms ( DQMStore::IBooker &  ibooker,
edm::Run const &  iRun,
edm::EventSetup const &  iSetup 
)
overridevirtual

Implements DQMEDAnalyzer.

Definition at line 257 of file Phase2OTValidateCluster.cc.

References bookLayerHistos(), edm::ESWatcher< T >::check(), config_, TrackerGeometry::detUnits(), edm::ParameterSet::getParameter(), GeomDetEnumerators::P2PXB, GeomDetEnumerators::P2PXEC, AlCaHLTBitMon_QueryRunRegistry::string, and tkGeom_.

259  {
260  std::string top_folder = config_.getParameter<std::string>("TopFolderName");
261  edm::LogInfo("Phase2OTValidateCluster") << " Booking Histograms in: " << top_folder;
262 
263  edm::ESWatcher<TrackerDigiGeometryRecord> theTkDigiGeomWatcher;
264  if (theTkDigiGeomWatcher.check(iSetup)) {
265  for (auto const& det_u : tkGeom_->detUnits()) {
266  //Always check TrackerNumberingBuilder before changing this part
267  if ((det_u->subDetector() == GeomDetEnumerators::SubDetector::P2PXB ||
268  det_u->subDetector() == GeomDetEnumerators::SubDetector::P2PXEC))
269  continue; //continue if Pixel
270  uint32_t detId_raw = det_u->geographicalId().rawId();
271  bookLayerHistos(ibooker, detId_raw, top_folder);
272  }
273  }
274 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
TrackerGeometry::detUnits
const DetContainer & detUnits() const override
Returm a vector of all GeomDet.
Definition: TrackerGeometry.h:61
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
Phase2OTValidateCluster::tkGeom_
const TrackerGeometry * tkGeom_
Definition: Phase2OTValidateCluster.cc:94
Phase2OTValidateCluster::bookLayerHistos
void bookLayerHistos(DQMStore::IBooker &ibooker, uint32_t det_it, const std::string &subdir)
Definition: Phase2OTValidateCluster.cc:277
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
edm::ESWatcher::check
bool check(const edm::EventSetup &iSetup)
Definition: ESWatcher.h:57

◆ bookLayerHistos()

void Phase2OTValidateCluster::bookLayerHistos ( DQMStore::IBooker &  ibooker,
uint32_t  det_it,
const std::string &  subdir 
)
private

Definition at line 277 of file Phase2OTValidateCluster.cc.

References phase2tkutil::book1DFromPSet(), dqm::implementation::NavigatorBase::cd(), config_, Phase2OTValidateCluster::ClusterMEs::deltaX_P, Phase2OTValidateCluster::ClusterMEs::deltaX_P_primary, Phase2OTValidateCluster::ClusterMEs::deltaX_S, Phase2OTValidateCluster::ClusterMEs::deltaX_S_primary, Phase2OTValidateCluster::ClusterMEs::deltaY_P, Phase2OTValidateCluster::ClusterMEs::deltaY_P_primary, Phase2OTValidateCluster::ClusterMEs::deltaY_S, Phase2OTValidateCluster::ClusterMEs::deltaY_S_primary, ctppsRandomDQMSource_cfi::folderName, TrackerGeometry::getDetectorType(), phase2tkutil::getOTHistoId(), edm::ParameterSet::getParameter(), layerMEs_, TrackerGeometry::Ph2PSP, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, mkLumiAveragedPlots::subdir, tkGeom_, and tTopo_.

Referenced by bookHistograms().

277  {
278  std::string folderName = phase2tkutil::getOTHistoId(det_id, tTopo_);
279  if (folderName.empty()) {
280  edm::LogWarning("Phase2OTValidateCluster") << ">>>> Invalid histo_id ";
281  return;
282  }
283 
284  if (layerMEs_.find(folderName) == layerMEs_.end()) {
285  ibooker.cd();
286  edm::LogInfo("Phase2TrackerValidateDigi") << " Booking Histograms in: " << subdir + '/' + folderName;
287  ClusterMEs local_mes;
288  if (tkGeom_->getDetectorType(det_id) == TrackerGeometry::ModuleType::Ph2PSP) {
289  ibooker.setCurrentFolder(subdir + '/' + folderName);
290 
291  local_mes.deltaX_P =
292  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_X_Pixel"), ibooker);
293 
294  local_mes.deltaY_P =
295  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_Y_Pixel"), ibooker);
296 
297  // Puting primary digis in a subfolder
298  ibooker.setCurrentFolder(subdir + '/' + folderName + "/PrimarySimHits");
299 
300  local_mes.deltaX_P_primary =
301  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_X_Pixel_Primary"), ibooker);
302 
303  local_mes.deltaY_P_primary =
304  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_Y_Pixel_Primary"), ibooker);
305  }
306  ibooker.setCurrentFolder(subdir + '/' + folderName);
307 
308  local_mes.deltaX_S =
309  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_X_Strip"), ibooker);
310 
311  local_mes.deltaY_S =
312  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_Y_Strip"), ibooker);
313 
314  // Puting primary digis in a subfolder
315  ibooker.setCurrentFolder(subdir + '/' + folderName + "/PrimarySimHits");
316 
317  local_mes.deltaX_S_primary =
318  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_X_Strip_Primary"), ibooker);
319 
320  local_mes.deltaY_S_primary =
321  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_Y_Strip_Primary"), ibooker);
322 
323  layerMEs_.emplace(folderName, local_mes);
324  }
325 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:36
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
Phase2OTValidateCluster::tkGeom_
const TrackerGeometry * tkGeom_
Definition: Phase2OTValidateCluster.cc:94
Phase2OTValidateCluster::tTopo_
const TrackerTopology * tTopo_
Definition: Phase2OTValidateCluster.cc:95
phase2tkutil::getOTHistoId
std::string getOTHistoId(uint32_t det_id, const TrackerTopology *tTopo)
Definition: TrackerPhase2DQMUtil.cc:25
TrackerGeometry::getDetectorType
ModuleType getDetectorType(DetId) const
Definition: TrackerGeometry.cc:247
phase2tkutil::book1DFromPSet
MonitorElement * book1DFromPSet(const edm::ParameterSet &hpars, DQMStore::IBooker &ibooker)
Definition: TrackerPhase2DQMUtil.cc:50
ctppsRandomDQMSource_cfi.folderName
folderName
Definition: ctppsRandomDQMSource_cfi.py:7
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
Phase2OTValidateCluster::layerMEs_
std::map< std::string, ClusterMEs > layerMEs_
Definition: Phase2OTValidateCluster.cc:82
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122

◆ dqmBeginRun()

void Phase2OTValidateCluster::dqmBeginRun ( const edm::Run &  iRun,
const edm::EventSetup &  iSetup 
)
overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 126 of file Phase2OTValidateCluster.cc.

References geomToken_, edm::EventSetup::getData(), tkGeom_, topoToken_, and tTopo_.

126  {
127  tkGeom_ = &iSetup.getData(geomToken_);
128  tTopo_ = &iSetup.getData(topoToken_);
129 }
edm::EventSetup::getData
T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)
Definition: EventSetup.h:119
Phase2OTValidateCluster::tkGeom_
const TrackerGeometry * tkGeom_
Definition: Phase2OTValidateCluster.cc:94
Phase2OTValidateCluster::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Definition: Phase2OTValidateCluster.cc:93
Phase2OTValidateCluster::tTopo_
const TrackerTopology * tTopo_
Definition: Phase2OTValidateCluster.cc:95
Phase2OTValidateCluster::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: Phase2OTValidateCluster.cc:92

◆ fillDescriptions()

void Phase2OTValidateCluster::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 341 of file Phase2OTValidateCluster.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

341  {
342  edm::ParameterSetDescription desc;
343  //for macro-pixel sensors
344  std::string mptag = "macro-pixel sensor";
345  std::string striptag = "strip sensor";
346  {
347  edm::ParameterSetDescription psd0;
348  psd0.add<std::string>("name", "Delta_X_Pixel");
349  psd0.add<std::string>("title", "#Delta X " + mptag + ";Cluster resolution X coordinate [#mum]");
350  psd0.add<bool>("switch", true);
351  psd0.add<double>("xmax", 250);
352  psd0.add<double>("xmin", -250);
353  psd0.add<int>("NxBins", 100);
354  desc.add<edm::ParameterSetDescription>("Delta_X_Pixel", psd0);
355  }
356  {
357  edm::ParameterSetDescription psd0;
358  psd0.add<std::string>("name", "Delta_Y_Pixel");
359  psd0.add<std::string>("title", "#Delta Y " + mptag + ";Cluster resolution Y coordinate [#mum]");
360  psd0.add<bool>("switch", true);
361  psd0.add<double>("xmin", -1500);
362  psd0.add<double>("xmax", 1500);
363  psd0.add<int>("NxBins", 100);
364  desc.add<edm::ParameterSetDescription>("Delta_Y_Pixel", psd0);
365  }
366  {
367  edm::ParameterSetDescription psd0;
368  psd0.add<std::string>("name", "Delta_X_Pixel_Primary");
369  psd0.add<std::string>("title", "#Delta X " + mptag + ";cluster resolution X coordinate [#mum]");
370  psd0.add<bool>("switch", true);
371  psd0.add<double>("xmin", -250);
372  psd0.add<double>("xmax", 250);
373  psd0.add<int>("NxBins", 100);
374  desc.add<edm::ParameterSetDescription>("Delta_X_Pixel_Primary", psd0);
375  }
376  {
377  edm::ParameterSetDescription psd0;
378  psd0.add<std::string>("name", "Delta_Y_Pixel_Primary");
379  psd0.add<std::string>("title", "#Delta Y " + mptag + ";cluster resolution Y coordinate [#mum]");
380  psd0.add<bool>("switch", true);
381  psd0.add<double>("xmin", -1500);
382  psd0.add<double>("xmax", 1500);
383  psd0.add<int>("NxBins", 100);
384  desc.add<edm::ParameterSetDescription>("Delta_Y_Pixel_Primary", psd0);
385  }
386 
387  //strip sensors
388  {
389  edm::ParameterSetDescription psd0;
390  psd0.add<std::string>("name", "Delta_X_Strip");
391  psd0.add<std::string>("title", "#Delta X " + striptag + ";Cluster resolution X coordinate [#mum]");
392  psd0.add<bool>("switch", true);
393  psd0.add<double>("xmin", -250);
394  psd0.add<double>("xmax", 250);
395  psd0.add<int>("NxBins", 100);
396  desc.add<edm::ParameterSetDescription>("Delta_X_Strip", psd0);
397  }
398  {
399  edm::ParameterSetDescription psd0;
400  psd0.add<std::string>("name", "Delta_Y_Strip");
401  psd0.add<std::string>("title", "#Delta Y " + striptag + ";Cluster resolution Y coordinate [cm]");
402  psd0.add<double>("xmin", -5.0);
403  psd0.add<bool>("switch", true);
404  psd0.add<double>("xmax", 5.0);
405  psd0.add<int>("NxBins", 100);
406  desc.add<edm::ParameterSetDescription>("Delta_Y_Strip", psd0);
407  }
408  {
409  edm::ParameterSetDescription psd0;
410  psd0.add<std::string>("name", "Delta_X_Strip_Primary");
411  psd0.add<std::string>("title", "#Delta X " + striptag + ";Cluster resolution X coordinate [#mum]");
412  psd0.add<bool>("switch", true);
413  psd0.add<double>("xmin", -250);
414  psd0.add<double>("xmax", 250);
415  psd0.add<int>("NxBins", 100);
416  desc.add<edm::ParameterSetDescription>("Delta_X_Strip_Primary", psd0);
417  }
418  {
419  edm::ParameterSetDescription psd0;
420  psd0.add<std::string>("name", "Delta_Y_Strip_Primary");
421  psd0.add<std::string>("title", "#Delta Y " + striptag + ";Cluster resolution Y coordinate [cm]");
422  psd0.add<double>("xmin", -5.0);
423  psd0.add<bool>("switch", true);
424  psd0.add<double>("xmax", 5.0);
425  psd0.add<int>("NxBins", 100);
426  desc.add<edm::ParameterSetDescription>("Delta_Y_Strip_Primary", psd0);
427  }
428  desc.add<std::string>("TopFolderName", "TrackerPhase2OTClusterV");
429  desc.add<edm::InputTag>("ClusterSource", edm::InputTag("siPhase2Clusters"));
430  desc.add<edm::InputTag>("OuterTrackerDigiSimLinkSource", edm::InputTag("simSiPixelDigis", "Tracker"));
431  desc.add<edm::InputTag>("simtracks", edm::InputTag("g4SimHits"));
432  desc.add<double>("SimTrackMinPt", 0.0);
433  desc.add<std::vector<edm::InputTag>>("PSimHitSource",
434  {
435  edm::InputTag("g4SimHits:TrackerHitsTIBLowTof"),
436  edm::InputTag("g4SimHits:TrackerHitsTIBHighTof"),
437  edm::InputTag("g4SimHits:TrackerHitsTIDLowTof"),
438  edm::InputTag("g4SimHits:TrackerHitsTIDHighTof"),
439  edm::InputTag("g4SimHits:TrackerHitsTOBLowTof"),
440  edm::InputTag("g4SimHits:TrackerHitsTOBHighTof"),
441  edm::InputTag("g4SimHits:TrackerHitsTECLowTof"),
442  edm::InputTag("g4SimHits:TrackerHitsTECHighTof"),
443  edm::InputTag("g4SimHits:TrackerHitsPixelBarrelLowTof"),
444  edm::InputTag("g4SimHits:TrackerHitsPixelBarrelHighTof"),
445  edm::InputTag("g4SimHits:TrackerHitsPixelEndcapLowTof"),
446  edm::InputTag("g4SimHits:TrackerHitsPixelEndcapHighTof"),
447  });
448  descriptions.add("Phase2OTValidateCluster", desc);
449 }
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57

◆ fillOTHistos()

void Phase2OTValidateCluster::fillOTHistos ( const edm::Event &  iEvent,
const std::vector< edm::Handle< edm::PSimHitContainer >> &  simHits,
const std::map< unsigned int, SimTrack > &  simTracks 
)
private

Definition at line 155 of file Phase2OTValidateCluster.cc.

References PVValHelper::add(), clustersToken_, phase2tkutil::cmtomicron, Phase2OTValidateCluster::ClusterMEs::deltaX_P, Phase2OTValidateCluster::ClusterMEs::deltaX_P_primary, Phase2OTValidateCluster::ClusterMEs::deltaX_S, Phase2OTValidateCluster::ClusterMEs::deltaX_S_primary, Phase2OTValidateCluster::ClusterMEs::deltaY_P, Phase2OTValidateCluster::ClusterMEs::deltaY_P_primary, Phase2OTValidateCluster::ClusterMEs::deltaY_S, Phase2OTValidateCluster::ClusterMEs::deltaY_S_primary, PVValHelper::dx, PVValHelper::dy, vertexPlots::e4, dqm::impl::MonitorElement::Fill(), TrackerGeometry::getDetectorType(), phase2tkutil::getOTHistoId(), getSimTrackId(), mps_fire::i, TrackerGeometry::idToDetUnit(), iEvent, phase2tkutil::isPrimary(), dqmiolumiharvest::j, layerMEs_, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, TrackerGeometry::Ph2PSP, TrackerGeometry::Ph2PSS, TrackerGeometry::Ph2SS, Phase2TrackerDigi::pixelToChannel(), FastTrackerRecHitCombiner_cfi::simHits, simOTLinksToken_, TrackCandidateProducer_cfi::simTracks, jetUpdater_cfi::sort, AlCaHLTBitMon_QueryRunRegistry::string, tkGeom_, tTopo_, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by analyze().

157  {
158  // Getting the clusters
159  const auto& clusterHandle = iEvent.getHandle(clustersToken_);
160  // Getting PixelDigiSimLinks
161  const auto& pixelSimLinksHandle = iEvent.getHandle(simOTLinksToken_);
162 
163  // Number of clusters
164  std::map<std::string, unsigned int> nPrimarySimHits[3];
165  std::map<std::string, unsigned int> nOtherSimHits[3];
166  for (const auto& DSVItr : *clusterHandle) {
167  // Getting the id of detector unit
168  uint32_t rawid = DSVItr.detId();
169  DetId detId(rawid);
170  const GeomDetUnit* geomDetUnit(tkGeom_->idToDetUnit(detId));
171  if (!geomDetUnit)
172  continue;
173  TrackerGeometry::ModuleType mType = tkGeom_->getDetectorType(detId);
174 
175  std::string folderkey = phase2tkutil::getOTHistoId(detId, tTopo_);
176  for (const auto& clusterItr : DSVItr) {
177  MeasurementPoint mpCluster(clusterItr.center(), clusterItr.column() + 0.5);
178  Local3DPoint localPosCluster = geomDetUnit->topology().localPosition(mpCluster);
179 
180  // Get simTracks from the cluster
181  std::vector<unsigned int> clusterSimTrackIds;
182  for (unsigned int i(0); i < clusterItr.size(); ++i) {
183  unsigned int channel(Phase2TrackerDigi::pixelToChannel(clusterItr.firstRow() + i, clusterItr.column()));
184  std::vector<unsigned int> simTrackIds(getSimTrackId(pixelSimLinksHandle, detId, channel));
185  for (auto it : simTrackIds) {
186  bool add = true;
187  for (unsigned int j = 0; j < clusterSimTrackIds.size(); ++j) {
188  // only save simtrackids that are not present yet
189  if (it == clusterSimTrackIds.at(j))
190  add = false;
191  }
192  if (add)
193  clusterSimTrackIds.push_back(it);
194  }
195  }
196  std::sort(clusterSimTrackIds.begin(), clusterSimTrackIds.end());
197  const PSimHit* closestSimHit = nullptr;
198  float mind = 1e4;
199  // Get the SimHit
200  for (const auto& psimhitCont : simHits) {
201  for (const auto& simhitIt : *psimhitCont) {
202  if (rawid == simhitIt.detUnitId()) {
203  auto it = std::lower_bound(clusterSimTrackIds.begin(), clusterSimTrackIds.end(), simhitIt.trackId());
204  if (it != clusterSimTrackIds.end() && *it == simhitIt.trackId()) {
205  float dx = simhitIt.localPosition().x() - localPosCluster.x();
206  float dy = simhitIt.localPosition().y() - localPosCluster.y();
207  float dist = dx * dx + dy * dy;
208  if (!closestSimHit || dist < mind) {
209  mind = dist;
210  closestSimHit = &simhitIt;
211  }
212  }
213  }
214  } //end loop over PSimhitcontainers
215  } //end loop over simHits
216 
217  if (!closestSimHit)
218  continue;
219  // only look at simhits from highpT tracks
220  auto simTrackIt(simTracks.find(closestSimHit->trackId()));
221  if (simTrackIt == simTracks.end())
222  continue;
223 
224  Local3DPoint localPosSimHit(closestSimHit->localPosition());
225  const float deltaX = localPosCluster.x() - localPosSimHit.x();
226  const float deltaY = localPosCluster.y() - localPosSimHit.y();
227 
228  auto layerMEit = layerMEs_.find(folderkey);
229  if (layerMEit == layerMEs_.end())
230  continue;
231 
232  ClusterMEs& local_mes = layerMEit->second;
233  if (mType == TrackerGeometry::ModuleType::Ph2PSP) {
234  local_mes.deltaX_P->Fill(phase2tkutil::cmtomicron * deltaX);
235  local_mes.deltaY_P->Fill(phase2tkutil::cmtomicron * deltaY);
236  } else if (mType == TrackerGeometry::ModuleType::Ph2PSS || mType == TrackerGeometry::ModuleType::Ph2SS) {
237  local_mes.deltaX_S->Fill(phase2tkutil::cmtomicron * deltaX);
238  local_mes.deltaY_S->Fill(deltaY);
239  }
240  // Primary particles only
241  if (phase2tkutil::isPrimary(simTrackIt->second, closestSimHit)) {
242  if (mType == TrackerGeometry::ModuleType::Ph2PSP) {
243  local_mes.deltaX_P_primary->Fill(phase2tkutil::cmtomicron * deltaX);
244  local_mes.deltaY_P_primary->Fill(phase2tkutil::cmtomicron * deltaY);
245  } else if (mType == TrackerGeometry::ModuleType::Ph2PSS || mType == TrackerGeometry::ModuleType::Ph2SS) {
246  local_mes.deltaX_S_primary->Fill(phase2tkutil::cmtomicron * deltaX);
247  local_mes.deltaY_S_primary->Fill(deltaY);
248  }
249  }
250  }
251  }
252 }
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:30
mps_fire.i
i
Definition: mps_fire.py:429
TrackerGeometry::idToDetUnit
const TrackerGeomDet * idToDetUnit(DetId) const override
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: TrackerGeometry.cc:183
phase2tkutil::cmtomicron
static constexpr float cmtomicron
Definition: TrackerPhase2ValidationUtil.h:9
phase2tkutil::isPrimary
bool isPrimary(const SimTrack &simTrk, const PSimHit *simHit)
Definition: TrackerPhase2ValidationUtil.cc:2
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
Phase2OTValidateCluster::tkGeom_
const TrackerGeometry * tkGeom_
Definition: Phase2OTValidateCluster.cc:94
Phase2OTValidateCluster::simOTLinksToken_
edm::EDGetTokenT< edm::DetSetVector< PixelDigiSimLink > > simOTLinksToken_
Definition: Phase2OTValidateCluster.cc:87
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
pfDeepBoostedJetPreprocessParams_cfi.lower_bound
lower_bound
Definition: pfDeepBoostedJetPreprocessParams_cfi.py:15
iEvent
int iEvent
Definition: GenABIO.cc:224
Phase2OTValidateCluster::tTopo_
const TrackerTopology * tTopo_
Definition: Phase2OTValidateCluster.cc:95
phase2tkutil::getOTHistoId
std::string getOTHistoId(uint32_t det_id, const TrackerTopology *tTopo)
Definition: TrackerPhase2DQMUtil.cc:25
TrackCandidateProducer_cfi.simTracks
simTracks
Definition: TrackCandidateProducer_cfi.py:15
TrackerGeometry::getDetectorType
ModuleType getDetectorType(DetId) const
Definition: TrackerGeometry.cc:247
Phase2TrackerDigi::pixelToChannel
static PackedDigiType pixelToChannel(unsigned int row, unsigned int col)
Definition: Phase2TrackerDigi.h:43
DetId
Definition: DetId.h:17
Phase2OTValidateCluster::layerMEs_
std::map< std::string, ClusterMEs > layerMEs_
Definition: Phase2OTValidateCluster.cc:82
PVValHelper::add
void add(std::map< std::string, TH1 *> &h, TH1 *hist)
Definition: PVValidationHelpers.cc:12
Phase2OTValidateCluster::clustersToken_
edm::EDGetTokenT< Phase2TrackerCluster1DCollectionNew > clustersToken_
Definition: Phase2OTValidateCluster.cc:90
Phase2OTValidateCluster::getSimTrackId
std::vector< unsigned int > getSimTrackId(const edm::Handle< edm::DetSetVector< PixelDigiSimLink >> &pixelSimLinks, const DetId &detId, unsigned int channel)
Definition: Phase2OTValidateCluster.cc:327
FastTrackerRecHitCombiner_cfi.simHits
simHits
Definition: FastTrackerRecHitCombiner_cfi.py:5

◆ getSimTrackId()

std::vector< unsigned int > Phase2OTValidateCluster::getSimTrackId ( const edm::Handle< edm::DetSetVector< PixelDigiSimLink >> &  pixelSimLinks,
const DetId &  detId,
unsigned int  channel 
)
private

Definition at line 327 of file Phase2OTValidateCluster.cc.

Referenced by fillOTHistos().

328  {
329  std::vector<unsigned int> retvec;
330  edm::DetSetVector<PixelDigiSimLink>::const_iterator DSViter(pixelSimLinks->find(detId));
331  if (DSViter == pixelSimLinks->end())
332  return retvec;
333  for (edm::DetSet<PixelDigiSimLink>::const_iterator it = DSViter->data.begin(); it != DSViter->data.end(); ++it) {
334  if (channel == it->channel()) {
335  retvec.push_back(it->SimTrackId());
336  }
337  }
338  return retvec;
339 }
edm::DetSetVector::find
iterator find(det_id_type id)
Definition: DetSetVector.h:264
edm::DetSetVector::end
iterator end()
Return the off-the-end iterator.
Definition: DetSetVector.h:325
edm::DetSet::const_iterator
collection_type::const_iterator const_iterator
Definition: DetSet.h:31
edm::DetSetVector::const_iterator
collection_type::const_iterator const_iterator
Definition: DetSetVector.h:102

Member Data Documentation

◆ clustersToken_

edm::EDGetTokenT<Phase2TrackerCluster1DCollectionNew> Phase2OTValidateCluster::clustersToken_
private

Definition at line 90 of file Phase2OTValidateCluster.cc.

Referenced by fillOTHistos().

◆ config_

edm::ParameterSet Phase2OTValidateCluster::config_
private

Definition at line 84 of file Phase2OTValidateCluster.cc.

Referenced by bookHistograms(), and bookLayerHistos().

◆ geomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> Phase2OTValidateCluster::geomToken_
private

Definition at line 92 of file Phase2OTValidateCluster.cc.

Referenced by dqmBeginRun().

◆ layerMEs_

std::map<std::string, ClusterMEs> Phase2OTValidateCluster::layerMEs_
private

Definition at line 82 of file Phase2OTValidateCluster.cc.

Referenced by bookLayerHistos(), and fillOTHistos().

◆ pSimHitSrc_

std::vector<edm::InputTag> Phase2OTValidateCluster::pSimHitSrc_
private

Definition at line 91 of file Phase2OTValidateCluster.cc.

Referenced by Phase2OTValidateCluster().

◆ simHitTokens_

std::vector<edm::EDGetTokenT<edm::PSimHitContainer> > Phase2OTValidateCluster::simHitTokens_
private

Definition at line 86 of file Phase2OTValidateCluster.cc.

Referenced by analyze(), and Phase2OTValidateCluster().

◆ simITLinksToken_

edm::EDGetTokenT<edm::DetSetVector<PixelDigiSimLink> > Phase2OTValidateCluster::simITLinksToken_
private

Definition at line 88 of file Phase2OTValidateCluster.cc.

◆ simOTLinksToken_

edm::EDGetTokenT<edm::DetSetVector<PixelDigiSimLink> > Phase2OTValidateCluster::simOTLinksToken_
private

Definition at line 87 of file Phase2OTValidateCluster.cc.

Referenced by fillOTHistos().

◆ simtrackminpt_

double Phase2OTValidateCluster::simtrackminpt_
private

Definition at line 85 of file Phase2OTValidateCluster.cc.

Referenced by analyze().

◆ simTracksToken_

edm::EDGetTokenT<edm::SimTrackContainer> Phase2OTValidateCluster::simTracksToken_
private

Definition at line 89 of file Phase2OTValidateCluster.cc.

Referenced by analyze().

◆ tkGeom_

const TrackerGeometry* Phase2OTValidateCluster::tkGeom_ = nullptr
private

Definition at line 94 of file Phase2OTValidateCluster.cc.

Referenced by bookHistograms(), bookLayerHistos(), dqmBeginRun(), and fillOTHistos().

◆ topoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> Phase2OTValidateCluster::topoToken_
private

Definition at line 93 of file Phase2OTValidateCluster.cc.

Referenced by dqmBeginRun().

◆ tTopo_

const TrackerTopology* Phase2OTValidateCluster::tTopo_ = nullptr
private

Definition at line 95 of file Phase2OTValidateCluster.cc.

Referenced by bookLayerHistos(), dqmBeginRun(), and fillOTHistos().

Generated for CMSSW Reference Manual by doxygen 1.8.14