#include <Phase2OTValidateCluster.cc>

Inheritance diagram for Phase2OTValidateCluster:

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_.

     : config_(iConfig),
       simtrackminpt_(config_.getParameter<double>("SimTrackMinPt")),
       simOTLinksToken_(consumes<edm::DetSetVector<PixelDigiSimLink>>(
           config_.getParameter<edm::InputTag>("OuterTrackerDigiSimLinkSource"))),
       simTracksToken_(consumes<edm::SimTrackContainer>(config_.getParameter<edm::InputTag>("simtracks"))),
       clustersToken_(
           consumes<Phase2TrackerCluster1DCollectionNew>(config_.getParameter<edm::InputTag>("ClusterSource"))),
       pSimHitSrc_(config_.getParameter<std::vector<edm::InputTag>>("PSimHitSource")),
       geomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>()),
       topoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>()) {
   edm::LogInfo("Phase2OTValidateCluster") << ">>> Construct Phase2OTValidateCluster ";
   for (const auto& itag : pSimHitSrc_)
     simHitTokens_.push_back(consumes<edm::PSimHitContainer>(itag));
 }

◆ ~Phase2OTValidateCluster()

Phase2OTValidateCluster::~Phase2OTValidateCluster ( )

override

Definition at line 118 of file Phase2OTValidateCluster.cc.

                                                   {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   edm::LogInfo("Phase2OTValidateCluster") << ">>> Destroy Phase2OTValidateCluster ";
 }

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_.

                                                                                          {
   // Getting simHits
   std::vector<edm::Handle<edm::PSimHitContainer>> simHits;
   for (const auto& itoken : simHitTokens_) {
     const auto& simHitHandle = iEvent.getHandle(itoken);
     if (!simHitHandle.isValid())
       continue;
     simHits.emplace_back(simHitHandle);
   }
   // Get the SimTracks
   const auto& simTracksRaw = iEvent.getHandle(simTracksToken_);
   // Rearrange the simTracks for ease of use <simTrackID, simTrack>
   SimTracksMap simTracks;
   for (edm::SimTrackContainer::const_iterator simTrackIt(simTracksRaw->begin()); simTrackIt != simTracksRaw->end();
        ++simTrackIt) {
     if (simTrackIt->momentum().pt() > simtrackminpt_) {
       simTracks.emplace(simTrackIt->trackId(), *simTrackIt);
     }
   }
   fillOTHistos(iEvent, simHits, simTracks);
 }

◆ 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_.

                                                                           {
   std::string top_folder = config_.getParameter<std::string>("TopFolderName");
   edm::LogInfo("Phase2OTValidateCluster") << " Booking Histograms in: " << top_folder;
 
   edm::ESWatcher<TrackerDigiGeometryRecord> theTkDigiGeomWatcher;
   if (theTkDigiGeomWatcher.check(iSetup)) {
     for (auto const& det_u : tkGeom_->detUnits()) {
       //Always check TrackerNumberingBuilder before changing this part
       if ((det_u->subDetector() == GeomDetEnumerators::SubDetector::P2PXB ||
            det_u->subDetector() == GeomDetEnumerators::SubDetector::P2PXEC))
         continue;  //continue if Pixel
       uint32_t detId_raw = det_u->geographicalId().rawId();
       bookLayerHistos(ibooker, detId_raw, top_folder);
     }
   }
 }

◆ bookLayerHistos()

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

private

Definition at line 277 of file Phase2OTValidateCluster.cc.

Referenced by bookHistograms().

                                                                                                                 {
   std::string folderName = phase2tkutil::getOTHistoId(det_id, tTopo_);
   if (folderName.empty()) {
     edm::LogWarning("Phase2OTValidateCluster") << ">>>> Invalid histo_id ";
     return;
   }
 
   if (layerMEs_.find(folderName) == layerMEs_.end()) {
     ibooker.cd();
     edm::LogInfo("Phase2TrackerValidateDigi") << " Booking Histograms in: " << subdir + '/' + folderName;
     ClusterMEs local_mes;
     if (tkGeom_->getDetectorType(det_id) == TrackerGeometry::ModuleType::Ph2PSP) {
       ibooker.setCurrentFolder(subdir + '/' + folderName);
 
       local_mes.deltaX_P =
           phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_X_Pixel"), ibooker);
 
       local_mes.deltaY_P =
           phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_Y_Pixel"), ibooker);
 
       // Puting primary digis in a subfolder
       ibooker.setCurrentFolder(subdir + '/' + folderName + "/PrimarySimHits");
 
       local_mes.deltaX_P_primary =
           phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_X_Pixel_Primary"), ibooker);
 
       local_mes.deltaY_P_primary =
           phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_Y_Pixel_Primary"), ibooker);
     }
     ibooker.setCurrentFolder(subdir + '/' + folderName);
 
     local_mes.deltaX_S =
         phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_X_Strip"), ibooker);
 
     local_mes.deltaY_S =
         phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_Y_Strip"), ibooker);
 
     // Puting primary digis in a subfolder
     ibooker.setCurrentFolder(subdir + '/' + folderName + "/PrimarySimHits");
 
     local_mes.deltaX_S_primary =
         phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_X_Strip_Primary"), ibooker);
 
     local_mes.deltaY_S_primary =
         phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Delta_Y_Strip_Primary"), ibooker);
 
     layerMEs_.emplace(folderName, local_mes);
   }
 }

◆ 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_.

                                                                                          {
   tkGeom_ = &iSetup.getData(geomToken_);
   tTopo_ = &iSetup.getData(topoToken_);
 }

◆ 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.

                                                                                          {
   edm::ParameterSetDescription desc;
   //for macro-pixel sensors
   std::string mptag = "macro-pixel sensor";
   std::string striptag = "strip sensor";
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Delta_X_Pixel");
     psd0.add<std::string>("title", "#Delta X " + mptag + ";Cluster resolution X coordinate [#mum]");
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmax", 250);
     psd0.add<double>("xmin", -250);
     psd0.add<int>("NxBins", 100);
     desc.add<edm::ParameterSetDescription>("Delta_X_Pixel", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Delta_Y_Pixel");
     psd0.add<std::string>("title", "#Delta Y " + mptag + ";Cluster resolution Y coordinate [#mum]");
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmin", -1500);
     psd0.add<double>("xmax", 1500);
     psd0.add<int>("NxBins", 100);
     desc.add<edm::ParameterSetDescription>("Delta_Y_Pixel", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Delta_X_Pixel_Primary");
     psd0.add<std::string>("title", "#Delta X " + mptag + ";cluster resolution X coordinate [#mum]");
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmin", -250);
     psd0.add<double>("xmax", 250);
     psd0.add<int>("NxBins", 100);
     desc.add<edm::ParameterSetDescription>("Delta_X_Pixel_Primary", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Delta_Y_Pixel_Primary");
     psd0.add<std::string>("title", "#Delta Y " + mptag + ";cluster resolution Y coordinate [#mum]");
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmin", -1500);
     psd0.add<double>("xmax", 1500);
     psd0.add<int>("NxBins", 100);
     desc.add<edm::ParameterSetDescription>("Delta_Y_Pixel_Primary", psd0);
   }
 
   //strip sensors
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Delta_X_Strip");
     psd0.add<std::string>("title", "#Delta X " + striptag + ";Cluster resolution X coordinate [#mum]");
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmin", -250);
     psd0.add<double>("xmax", 250);
     psd0.add<int>("NxBins", 100);
     desc.add<edm::ParameterSetDescription>("Delta_X_Strip", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Delta_Y_Strip");
     psd0.add<std::string>("title", "#Delta Y " + striptag + ";Cluster resolution Y coordinate [cm]");
     psd0.add<double>("xmin", -5.0);
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmax", 5.0);
     psd0.add<int>("NxBins", 100);
     desc.add<edm::ParameterSetDescription>("Delta_Y_Strip", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Delta_X_Strip_Primary");
     psd0.add<std::string>("title", "#Delta X " + striptag + ";Cluster resolution X coordinate [#mum]");
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmin", -250);
     psd0.add<double>("xmax", 250);
     psd0.add<int>("NxBins", 100);
     desc.add<edm::ParameterSetDescription>("Delta_X_Strip_Primary", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Delta_Y_Strip_Primary");
     psd0.add<std::string>("title", "#Delta Y " + striptag + ";Cluster resolution Y coordinate [cm]");
     psd0.add<double>("xmin", -5.0);
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmax", 5.0);
     psd0.add<int>("NxBins", 100);
     desc.add<edm::ParameterSetDescription>("Delta_Y_Strip_Primary", psd0);
   }
   desc.add<std::string>("TopFolderName", "TrackerPhase2OTClusterV");
   desc.add<edm::InputTag>("ClusterSource", edm::InputTag("siPhase2Clusters"));
   desc.add<edm::InputTag>("OuterTrackerDigiSimLinkSource", edm::InputTag("simSiPixelDigis", "Tracker"));
   desc.add<edm::InputTag>("simtracks", edm::InputTag("g4SimHits"));
   desc.add<double>("SimTrackMinPt", 0.0);
   desc.add<std::vector<edm::InputTag>>("PSimHitSource",
                                        {
                                            edm::InputTag("g4SimHits:TrackerHitsTIBLowTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsTIBHighTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsTIDLowTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsTIDHighTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsTOBLowTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsTOBHighTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsTECLowTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsTECHighTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsPixelBarrelLowTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsPixelBarrelHighTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsPixelEndcapLowTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsPixelEndcapHighTof"),
                                        });
   descriptions.add("Phase2OTValidateCluster", desc);
 }

◆ 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, hcalRecHitTable_cff::detId, 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().

                                                                                             {
   // Getting the clusters
   const auto& clusterHandle = iEvent.getHandle(clustersToken_);
   // Getting PixelDigiSimLinks
   const auto& pixelSimLinksHandle = iEvent.getHandle(simOTLinksToken_);
 
   // Number of clusters
   std::map<std::string, unsigned int> nPrimarySimHits[3];
   std::map<std::string, unsigned int> nOtherSimHits[3];
   for (const auto& DSVItr : *clusterHandle) {
     // Getting the id of detector unit
     uint32_t rawid = DSVItr.detId();
     DetId detId(rawid);
     const GeomDetUnit* geomDetUnit(tkGeom_->idToDetUnit(detId));
     if (!geomDetUnit)
       continue;
     TrackerGeometry::ModuleType mType = tkGeom_->getDetectorType(detId);
 
     std::string folderkey = phase2tkutil::getOTHistoId(detId, tTopo_);
     for (const auto& clusterItr : DSVItr) {
       MeasurementPoint mpCluster(clusterItr.center(), clusterItr.column() + 0.5);
       Local3DPoint localPosCluster = geomDetUnit->topology().localPosition(mpCluster);
 
       // Get simTracks from the cluster
       std::vector<unsigned int> clusterSimTrackIds;
       for (unsigned int i(0); i < clusterItr.size(); ++i) {
         unsigned int channel(Phase2TrackerDigi::pixelToChannel(clusterItr.firstRow() + i, clusterItr.column()));
         std::vector<unsigned int> simTrackIds(getSimTrackId(pixelSimLinksHandle, detId, channel));
         for (auto it : simTrackIds) {
           bool add = true;
           for (unsigned int j = 0; j < clusterSimTrackIds.size(); ++j) {
             // only save simtrackids that are not present yet
             if (it == clusterSimTrackIds.at(j))
               add = false;
           }
           if (add)
             clusterSimTrackIds.push_back(it);
         }
       }
       std::sort(clusterSimTrackIds.begin(), clusterSimTrackIds.end());
       const PSimHit* closestSimHit = nullptr;
       float mind = 1e4;
       // Get the SimHit
       for (const auto& psimhitCont : simHits) {
         for (const auto& simhitIt : *psimhitCont) {
           if (rawid == simhitIt.detUnitId()) {
             auto it = std::lower_bound(clusterSimTrackIds.begin(), clusterSimTrackIds.end(), simhitIt.trackId());
             if (it != clusterSimTrackIds.end() && *it == simhitIt.trackId()) {
               float dx = simhitIt.localPosition().x() - localPosCluster.x();
               float dy = simhitIt.localPosition().y() - localPosCluster.y();
               float dist = dx * dx + dy * dy;
               if (!closestSimHit || dist < mind) {
                 mind = dist;
                 closestSimHit = &simhitIt;
               }
             }
           }
         }  //end loop over PSimhitcontainers
       }    //end loop over simHits
 
       if (!closestSimHit)
         continue;
       // only look at simhits from highpT tracks
       auto simTrackIt(simTracks.find(closestSimHit->trackId()));
       if (simTrackIt == simTracks.end())
         continue;
 
       Local3DPoint localPosSimHit(closestSimHit->localPosition());
       const float deltaX = localPosCluster.x() - localPosSimHit.x();
       const float deltaY = localPosCluster.y() - localPosSimHit.y();
 
       auto layerMEit = layerMEs_.find(folderkey);
       if (layerMEit == layerMEs_.end())
         continue;
 
       ClusterMEs& local_mes = layerMEit->second;
       if (mType == TrackerGeometry::ModuleType::Ph2PSP) {
         local_mes.deltaX_P->Fill(phase2tkutil::cmtomicron * deltaX);
         local_mes.deltaY_P->Fill(phase2tkutil::cmtomicron * deltaY);
       } else if (mType == TrackerGeometry::ModuleType::Ph2PSS || mType == TrackerGeometry::ModuleType::Ph2SS) {
         local_mes.deltaX_S->Fill(phase2tkutil::cmtomicron * deltaX);
         local_mes.deltaY_S->Fill(deltaY);
       }
       // Primary particles only
       if (phase2tkutil::isPrimary(simTrackIt->second, closestSimHit)) {
         if (mType == TrackerGeometry::ModuleType::Ph2PSP) {
           local_mes.deltaX_P_primary->Fill(phase2tkutil::cmtomicron * deltaX);
           local_mes.deltaY_P_primary->Fill(phase2tkutil::cmtomicron * deltaY);
         } else if (mType == TrackerGeometry::ModuleType::Ph2PSS || mType == TrackerGeometry::ModuleType::Ph2SS) {
           local_mes.deltaX_S_primary->Fill(phase2tkutil::cmtomicron * deltaX);
           local_mes.deltaY_S_primary->Fill(deltaY);
         }
       }
     }
   }
 }

◆ 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.

References hcalRecHitTable_cff::detId.

Referenced by fillOTHistos().

                                                                                                                  {
   std::vector<unsigned int> retvec;
   edm::DetSetVector<PixelDigiSimLink>::const_iterator DSViter(pixelSimLinks->find(detId));
   if (DSViter == pixelSimLinks->end())
     return retvec;
   for (edm::DetSet<PixelDigiSimLink>::const_iterator it = DSViter->data.begin(); it != DSViter->data.end(); ++it) {
     if (channel == it->channel()) {
       retvec.push_back(it->SimTrackId());
     }
   }
   return retvec;
 }

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().

Classes

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ SimHitsMap

◆ SimTracksMap

Constructor & Destructor Documentation

◆ Phase2OTValidateCluster()

◆ ~Phase2OTValidateCluster()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ bookLayerHistos()

◆ dqmBeginRun()

◆ fillDescriptions()

◆ fillOTHistos()

◆ getSimTrackId()

Member Data Documentation

◆ clustersToken_

◆ config_

◆ geomToken_

◆ layerMEs_

◆ pSimHitSrc_

◆ simHitTokens_

◆ simITLinksToken_

◆ simOTLinksToken_

◆ simtrackminpt_

◆ simTracksToken_

◆ tkGeom_

◆ topoToken_

◆ tTopo_