#include <Phase2ITValidateCluster.cc>

Inheritance diagram for Phase2ITValidateCluster:

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

	Phase2ITValidateCluster (const edm::ParameterSet &)

	~Phase2ITValidateCluster () 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	fillITHistos (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 < edmNew::DetSetVector < SiPixelCluster > >	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_

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

bookLayer

Description: Validation plots tracker clusters.

Definition at line 49 of file Phase2ITValidateCluster.cc.

Member Typedef Documentation

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

Definition at line 51 of file Phase2ITValidateCluster.cc.

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

Definition at line 52 of file Phase2ITValidateCluster.cc.

Constructor & Destructor Documentation

Phase2ITValidateCluster::Phase2ITValidateCluster ( const edm::ParameterSet & iConfig )

explicit

Definition at line 98 of file Phase2ITValidateCluster.cc.

References pSimHitSrc_, and simHitTokens_.

     : config_(iConfig),
       simtrackminpt_(config_.getParameter<double>("SimTrackMinPt")),
       simITLinksToken_(consumes<edm::DetSetVector<PixelDigiSimLink>>(
           config_.getParameter<edm::InputTag>("InnerTrackerDigiSimLinkSource"))),
       simTracksToken_(consumes<edm::SimTrackContainer>(config_.getParameter<edm::InputTag>("simtracks"))),
       clustersToken_(
           consumes<edmNew::DetSetVector<SiPixelCluster>>(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("Phase2ITValidateCluster") << ">>> Construct Phase2ITValidateCluster ";
   for (const auto& itag : pSimHitSrc_)
     simHitTokens_.push_back(consumes<edm::PSimHitContainer>(itag));
 }

Phase2ITValidateCluster::~Phase2ITValidateCluster ( )

override

Definition at line 114 of file Phase2ITValidateCluster.cc.

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

Member Function Documentation

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 129 of file Phase2ITValidateCluster.cc.

References fillITHistos(), edm::Event::getHandle(), trackerHits::simHits, simHitTokens_, simtrackminpt_, 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);
     }
   }
   fillITHistos(iEvent, simHits, simTracks);
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 239 of file Phase2ITValidateCluster.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("Phase2ITValidateCluster") << " 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 not Pixel
       uint32_t detId_raw = det_u->geographicalId().rawId();
       bookLayerHistos(ibooker, detId_raw, top_folder);
     }
   }
 }

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

private

Definition at line 259 of file Phase2ITValidateCluster.cc.

References phase2tkutil::book1DFromPSet(), dqm::implementation::NavigatorBase::cd(), config_, Phase2ITValidateCluster::ClusterMEs::deltaX_P, Phase2ITValidateCluster::ClusterMEs::deltaX_P_primary, Phase2ITValidateCluster::ClusterMEs::deltaY_P, Phase2ITValidateCluster::ClusterMEs::deltaY_P_primary, phase2tkutil::getITHistoId(), edm::ParameterSet::getParameter(), layerMEs_, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and tTopo_.

Referenced by bookHistograms().

                                                                                                                 {
   std::string folderName = phase2tkutil::getITHistoId(det_id, tTopo_);
   if (folderName.empty()) {
     edm::LogWarning("Phase2ITValidateCluster") << ">>>> Invalid histo_id ";
     return;
   }
 
   if (layerMEs_.find(folderName) == layerMEs_.end()) {
     ibooker.cd();
     ibooker.setCurrentFolder(subdir + '/' + folderName);
     edm::LogInfo("Phase2TrackerValidateDigi") << " Booking Histograms in: " << subdir + '/' + folderName;
     ClusterMEs local_mes;
 
     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);
     layerMEs_.emplace(folderName, local_mes);
   }
 }

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 122 of file Phase2ITValidateCluster.cc.

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

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

void Phase2ITValidateCluster::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 304 of file Phase2ITValidateCluster.cc.

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

                                                                                          {
   edm::ParameterSetDescription desc;
   //for macro-pixel sensors
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Delta_X_Pixel");
     psd0.add<std::string>("title", "#Delta X;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 ;Cluster resolution Y coordinate [#mum]");
     psd0.add<double>("xmin", -250.);
     psd0.add<double>("xmax", 250.);
     psd0.add<bool>("switch", true);
     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 ;cluster resolution X coordinate [#mum]");
     psd0.add<double>("xmin", -250.);
     psd0.add<double>("xmax", 250.);
     psd0.add<bool>("switch", true);
     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 ;cluster resolution Y coordinate [#mum]");
     psd0.add<double>("xmin", -250.);
     psd0.add<double>("xmax", 250.);
     psd0.add<bool>("switch", true);
     psd0.add<int>("NxBins", 100);
     desc.add<edm::ParameterSetDescription>("Delta_Y_Pixel_Primary", psd0);
   }
 
   desc.add<std::string>("TopFolderName", "TrackerPhase2ITClusterV");
   desc.add<edm::InputTag>("ClusterSource", edm::InputTag("siPixelClusters"));
   desc.add<edm::InputTag>("InnerTrackerDigiSimLinkSource", edm::InputTag("simSiPixelDigis", "Pixel"));
   desc.add<edm::InputTag>("simtracks", edm::InputTag("g4SimHits"));
   desc.add<double>("SimTrackMinPt", 0.0);
   desc.add<std::vector<edm::InputTag>>("PSimHitSource",
                                        {
                                            edm::InputTag("g4SimHits:TrackerHitsPixelBarrelLowTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsPixelBarrelHighTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsPixelEndcapLowTof"),
                                            edm::InputTag("g4SimHits:TrackerHitsPixelEndcapHighTof"),
                                        });
   descriptions.add("Phase2ITValidateCluster", desc);
 }

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

private

Definition at line 151 of file Phase2ITValidateCluster.cc.

References PVValHelper::add(), clustersToken_, phase2tkutil::cmtomicron, Phase2ITValidateCluster::ClusterMEs::deltaX_P, Phase2ITValidateCluster::ClusterMEs::deltaX_P_primary, Phase2ITValidateCluster::ClusterMEs::deltaY_P, Phase2ITValidateCluster::ClusterMEs::deltaY_P_primary, PVValHelper::dx, PVValHelper::dy, dqm::impl::MonitorElement::Fill(), edm::Event::getHandle(), phase2tkutil::getITHistoId(), getSimTrackId(), TrackerGeometry::idToDetUnit(), phase2tkutil::isPrimary(), dqmiolumiharvest::j, layerMEs_, cuda_std::lower_bound(), PixelChannelIdentifier::pixelToChannel(), trackerHits::simHits, simITLinksToken_, AlCaHLTBitMon_QueryRunRegistry::string, tkGeom_, and tTopo_.

Referenced by analyze().

                                                                                             {
   // Getting the clusters
   const auto& clusterHandle = iEvent.getHandle(clustersToken_);
   // Getting PixelDigiSimLinks
   const auto& pixelSimLinksHandle = iEvent.getHandle(simITLinksToken_);
 
   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;
 
     std::string folderkey = phase2tkutil::getITHistoId(detId, tTopo_);
     for (const auto& clusterItr : DSVItr) {
       MeasurementPoint mpCluster(clusterItr.x(), clusterItr.y());
       Local3DPoint localPosCluster = geomDetUnit->topology().localPosition(mpCluster);
       // Get simTracks from the cluster
       std::vector<unsigned int> clusterSimTrackIds;
       for (int irow = clusterItr.minPixelRow(); irow <= clusterItr.maxPixelRow(); ++irow) {
         for (int icol = clusterItr.minPixelCol(); icol <= clusterItr.maxPixelCol(); ++icol) {
           uint32_t channel = PixelChannelIdentifier::pixelToChannel(irow, icol);
           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 = 10000.;
       // 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 double deltaX = phase2tkutil::cmtomicron * (localPosCluster.x() - localPosSimHit.x());
       const double deltaY = phase2tkutil::cmtomicron * (localPosCluster.y() - localPosSimHit.y());
 
       auto layerMEIt = layerMEs_.find(folderkey);
       if (layerMEIt == layerMEs_.end())
         continue;
 
       ClusterMEs& local_mes = layerMEIt->second;
       local_mes.deltaX_P->Fill(deltaX);
       local_mes.deltaY_P->Fill(deltaY);
       // Primary particles only
       if (phase2tkutil::isPrimary(simTrackIt->second, closestSimHit)) {
         local_mes.deltaX_P_primary->Fill(deltaX);
         local_mes.deltaY_P_primary->Fill(deltaY);
       }
     }
   }
 }

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

private

Definition at line 290 of file Phase2ITValidateCluster.cc.

Referenced by fillITHistos().

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