#include <DQM/SiStripMonitorApproximateCluster/plugins/SiStripMonitorApproximateCluster.cc>

Inheritance diagram for SiStripMonitorApproximateCluster:

Public Member Functions
	SiStripMonitorApproximateCluster (const edm::ParameterSet &)

	~SiStripMonitorApproximateCluster () override=default

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

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	analyze (const edm::Event &, const edm::EventSetup &) override

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

Private Attributes
siStripRawPrime::monitorApproxCluster	allClusters {}

edm::EDGetTokenT< SiStripApproximateClusterCollection >	approxClustersToken_

bool	compareClusters_

std::string	folder_

MonitorElement *	h2_CompareBarycenter_ {nullptr}

MonitorElement *	h2_CompareCharge_ {nullptr}

MonitorElement *	h2_CompareEndStrip_ {nullptr}

MonitorElement *	h2_CompareFirstStrip_ {nullptr}

MonitorElement *	h2_CompareSize_ {nullptr}

MonitorElement *	h_deltaBarycenter_ {nullptr}

MonitorElement *	h_deltaCharge_ {nullptr}

MonitorElement *	h_deltaEndStrip_ {nullptr}

MonitorElement *	h_deltaFirstStrip_ {nullptr}

MonitorElement *	h_deltaSize_ {nullptr}

MonitorElement *	h_isMatched_ {nullptr}

MonitorElement *	h_nclusters_

siStripRawPrime::monitorApproxCluster	matchedClusters {}

const edmNew::DetSetVector< SiStripCluster > *	stripClusterCollection_

edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > >	stripClustersToken_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

siStripRawPrime::monitorApproxCluster	unMatchedClusters {}

Additional Inherited Members
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

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: Monitor SiStripApproximateClusters and on-demand compare properties with original SiStripClusters

Definition at line 79 of file SiStripMonitorApproximateCluster.cc.

Constructor & Destructor Documentation

◆ SiStripMonitorApproximateCluster()

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

explicit

Definition at line 126 of file SiStripMonitorApproximateCluster.cc.

References compareClusters_, deDxTools::esConsumes(), edm::ParameterSet::getParameter(), stripClusterCollection_, stripClustersToken_, and tkGeomToken_.

     : folder_(iConfig.getParameter<std::string>("folder")),
       compareClusters_(iConfig.getParameter<bool>("compareClusters")),
       // Poducer name of input StripClusterCollection
       approxClustersToken_(consumes<SiStripApproximateClusterCollection>(
           iConfig.getParameter<edm::InputTag>("ApproxClustersProducer"))) {
   tkGeomToken_ = esConsumes();
   if (compareClusters_) {
     stripClustersToken_ =
         consumes<edmNew::DetSetVector<SiStripCluster>>(iConfig.getParameter<edm::InputTag>("ClustersProducer"));
   }
   stripClusterCollection_ = nullptr;
 }

◆ ~SiStripMonitorApproximateCluster()

SiStripMonitorApproximateCluster::~SiStripMonitorApproximateCluster ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 145 of file SiStripMonitorApproximateCluster.cc.

References funct::abs(), allClusters, approxClustersToken_, SiStripCluster::barycenter(), ALCARECOSiPixelCalSingleMuonDQM_cff::clusterCollection, compareClusters_, HLT_2023v12_cff::distance, edmNew::DetSet< T >::empty(), siStripRawPrime::monitorApproxCluster::fill(), dqm::impl::MonitorElement::Fill(), edmNew::DetSetVector< T >::find(), edm::EventSetup::getData(), h2_CompareBarycenter_, h2_CompareCharge_, h2_CompareEndStrip_, h2_CompareFirstStrip_, h2_CompareSize_, h_deltaBarycenter_, h_deltaCharge_, h_deltaEndStrip_, h_deltaFirstStrip_, h_deltaSize_, h_isMatched_, h_nclusters_, edmNew::DetSet< T >::id(), iEvent, edm::HandleBase::isValid(), matchedClusters, me0TriggerPseudoDigis_cff::nStrips, AlCaHLTBitMon_ParallelJobs::p, edm::Handle< T >::product(), StripGeomDetUnit::specificTopology(), TrackingMonitor_cfi::stripCluster, stripClusterCollection_, stripClustersToken_, tkGeomToken_, and unMatchedClusters.

                                                                                                   {
   using namespace edm;
 
   const auto& tkGeom = &iSetup.getData(tkGeomToken_);
   const auto tkDets = tkGeom->dets();
 
   // get collection of DetSetVector of clusters from Event
   edm::Handle<SiStripApproximateClusterCollection> approx_cluster_detsetvector;
   iEvent.getByToken(approxClustersToken_, approx_cluster_detsetvector);
   if (!approx_cluster_detsetvector.isValid()) {
     edm::LogError("SiStripMonitorApproximateCluster")
         << "SiStripApproximate cluster collection is not valid!" << std::endl;
 
     // if approximate clusters collection not available, then early return
     return;
   }
 
   // if requested to perform the comparison
   if (compareClusters_) {
     // get collection of DetSetVector of clusters from Event
     edm::Handle<edmNew::DetSetVector<SiStripCluster>> cluster_detsetvector;
     iEvent.getByToken(stripClustersToken_, cluster_detsetvector);
     if (!cluster_detsetvector.isValid()) {
       edm::LogError("SiStripMonitorApproximateCluster")
           << "Requested to perform comparison, but regular SiStrip cluster collection is not valid!" << std::endl;
       return;
     } else {
       stripClusterCollection_ = cluster_detsetvector.product();
     }
   }
 
   int nApproxClusters{0};
   const SiStripApproximateClusterCollection* clusterCollection = approx_cluster_detsetvector.product();
 
   for (const auto& detClusters : *clusterCollection) {
     edmNew::DetSet<SiStripCluster> strip_clusters_detset;
     const auto& detid = detClusters.id();  // get the detid of the current detset
 
     // starts here comaparison with regular clusters
     if (compareClusters_) {
       edmNew::DetSetVector<SiStripCluster>::const_iterator isearch =
           stripClusterCollection_->find(detid);  // search clusters of same detid
       strip_clusters_detset = (*isearch);
     }
 
     for (const auto& cluster : detClusters) {
       nApproxClusters++;
 
       // fill the full cluster collection
       allClusters.fill(cluster);
 
       if (compareClusters_ && !strip_clusters_detset.empty()) {
         // build the converted cluster for the matching
         uint16_t nStrips{0};
         auto det = std::find_if(tkDets.begin(), tkDets.end(), [detid](auto& elem) -> bool {
           return (elem->geographicalId().rawId() == detid);
         });
         const StripTopology& p = dynamic_cast<const StripGeomDetUnit*>(*det)->specificTopology();
         nStrips = p.nstrips() - 1;
 
         const auto convertedCluster = SiStripCluster(cluster, nStrips);
 
         float distance{9999.};
         const SiStripCluster* closestCluster{nullptr};
         for (const auto& stripCluster : strip_clusters_detset) {
           // by construction the approximated cluster width has same
           // size as the original cluster
 
           if (cluster.width() != stripCluster.size()) {
             continue;
           }
 
           float deltaBarycenter = convertedCluster.barycenter() - stripCluster.barycenter();
           if (std::abs(deltaBarycenter) < distance) {
             closestCluster = &stripCluster;
             distance = std::abs(deltaBarycenter);
           }
         }
 
         // Matching criteria:
         // - if exists a closest cluster in the DetId
         // - the size coincides with the original one
         if (closestCluster) {
           // comparisong plots
           h_deltaBarycenter_->Fill(closestCluster->barycenter() - convertedCluster.barycenter());
           h_deltaSize_->Fill(closestCluster->size() - convertedCluster.size());
           h_deltaCharge_->Fill(closestCluster->charge() - convertedCluster.charge());
           h_deltaFirstStrip_->Fill(closestCluster->firstStrip() - convertedCluster.firstStrip());
           h_deltaEndStrip_->Fill(closestCluster->endStrip() - convertedCluster.endStrip());
 
           h2_CompareBarycenter_->Fill(closestCluster->barycenter(), convertedCluster.barycenter());
           h2_CompareSize_->Fill(closestCluster->size(), convertedCluster.size());
           h2_CompareCharge_->Fill(closestCluster->charge(), convertedCluster.charge());
           h2_CompareFirstStrip_->Fill(closestCluster->firstStrip(), convertedCluster.firstStrip());
           h2_CompareEndStrip_->Fill(closestCluster->endStrip(), convertedCluster.endStrip());
 
           // monitoring plots
           matchedClusters.fill(cluster);
           h_isMatched_->Fill(1);
         } else {
           // monitoring plots
           unMatchedClusters.fill(cluster);
           h_isMatched_->Fill(-1);
         }
       }  // if we're doing the comparison cluster by cluster
 
     }  // loop on clusters in a detset
   }    // loop on the detset vector
 
   h_nclusters_->Fill(nApproxClusters);
 }

◆ bookHistograms()

void SiStripMonitorApproximateCluster::bookHistograms	(	DQMStore::IBooker &	ibook,
		edm::Run const &	run,
		edm::EventSetup const &	iSetup
	)

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 257 of file SiStripMonitorApproximateCluster.cc.

References allClusters, siStripRawPrime::monitorApproxCluster::book(), dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2I(), compareClusters_, f, folder_, dqm-mbProfile::format, dqm::impl::MonitorElement::getTH1F(), h2_CompareBarycenter_, h2_CompareCharge_, h2_CompareEndStrip_, h2_CompareFirstStrip_, h2_CompareSize_, h_deltaBarycenter_, h_deltaCharge_, h_deltaEndStrip_, h_deltaFirstStrip_, h_deltaSize_, h_isMatched_, h_nclusters_, matchedClusters, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, sistrip::STRIPS_PER_APV, and unMatchedClusters.

                                                                                    {
   ibook.setCurrentFolder(folder_);
   h_nclusters_ = ibook.book1D("numberOfClusters", "total N. of clusters;N. of clusters;#clusters", 500., 0., 500000.);
   allClusters.book(ibook, folder_);
 
   //  for comparisons
   if (compareClusters_) {
     // book monitoring for matche and unmatched clusters separately
     matchedClusters.book(ibook, fmt::format("{}/MatchedClusters", folder_));
     unMatchedClusters.book(ibook, fmt::format("{}/UnmatchedClusters", folder_));
 
     // 1D histograms
     ibook.setCurrentFolder(fmt::format("{}/ClusterComparisons", folder_));
     h_deltaBarycenter_ =
         ibook.book1D("deltaBarycenter", "#Delta barycenter;#Delta barycenter;cluster pairs", 101, -50.5, 50.5);
     h_deltaSize_ = ibook.book1D("deltaSize", "#Delta size;#Delta size;cluster pairs", 201, -100.5, 100.5);
     h_deltaCharge_ = ibook.book1D("deltaCharge", "#Delta charge;#Delta charge;cluster pairs", 401, -200.5, 200.5);
 
     h_deltaFirstStrip_ =
         ibook.book1D("deltaFirstStrip", "#Delta FirstStrip; #Delta firstStrip;cluster pairs", 101, -50.5, 50.5);
     h_deltaEndStrip_ =
         ibook.book1D("deltaEndStrip", "#Delta EndStrip; #Delta endStrip; cluster pairs", 101, -50.5, 50.5);
 
     // geometric constants
     constexpr int maxNStrips = 6 * sistrip::STRIPS_PER_APV;
     constexpr float minStrip = -0.5f;
     constexpr float maxStrip = maxNStrips - 0.5f;
 
     // cluster constants
     constexpr float maxCluSize = 50;
     constexpr float maxCluCharge = 700;
 
     // 2D histograms (use TH2I for counts to limit memory allocation)
     std::string toRep = "SiStrip Cluster Barycenter";
     h2_CompareBarycenter_ = ibook.book2I("compareBarycenter",
                                          fmt::sprintf("; %s;Approx %s", toRep, toRep),
                                          maxNStrips,
                                          minStrip,
                                          maxStrip,
                                          maxNStrips,
                                          minStrip,
                                          maxStrip);
 
     toRep = "SiStrip Cluster Size";
     h2_CompareSize_ = ibook.book2I("compareSize",
                                    fmt::sprintf("; %s;Approx %s", toRep, toRep),
                                    maxCluSize,
                                    -0.5f,
                                    maxCluSize - 0.5f,
                                    maxCluSize,
                                    -0.5f,
                                    maxCluSize - 0.5f);
 
     toRep = "SiStrip Cluster Charge";
     h2_CompareCharge_ = ibook.book2I("compareCharge",
                                      fmt::sprintf("; %s;Approx %s", toRep, toRep),
                                      maxCluCharge,
                                      -0.5f,
                                      maxCluCharge - 0.5f,
                                      maxCluCharge,
                                      -0.5f,
                                      maxCluCharge - 0.5f);
 
     toRep = "SiStrip Cluster First Strip number";
     h2_CompareFirstStrip_ = ibook.book2I("compareFirstStrip",
                                          fmt::sprintf("; %s;Approx %s", toRep, toRep),
                                          maxNStrips,
                                          minStrip,
                                          maxStrip,
                                          maxNStrips,
                                          minStrip,
                                          maxStrip);
 
     toRep = "SiStrip Cluster Last Strip number";
     h2_CompareEndStrip_ = ibook.book2I("compareLastStrip",
                                        fmt::sprintf("; %s;Approx %s", toRep, toRep),
                                        maxNStrips,
                                        minStrip,
                                        maxStrip,
                                        maxNStrips,
                                        minStrip,
                                        maxStrip);
 
     h_isMatched_ = ibook.book1D("isClusterMatched", "cluster matching;is matched?;#clusters", 3, -1.5, 1.5);
     h_isMatched_->getTH1F()->GetXaxis()->SetBinLabel(1, "Not matched");
     h_isMatched_->getTH1F()->GetXaxis()->SetBinLabel(3, "Matched");
   }
 }

◆ fillDescriptions()

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

static

Definition at line 349 of file SiStripMonitorApproximateCluster.cc.

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

                                                                                                   {
   edm::ParameterSetDescription desc;
   desc.setComment("Monitor SiStripApproximateCluster collection and compare with regular SiStrip clusters");
   desc.add<bool>("compareClusters", false)->setComment("if true, will compare with regualr Strip clusters");
   desc.add<edm::InputTag>("ApproxClustersProducer", edm::InputTag("hltSiStripClusters2ApproxClusters"))
       ->setComment("approxmate clusters collection");
   desc.add<edm::InputTag>("ClustersProducer", edm::InputTag("hltSiStripClusterizerForRawPrime"))
       ->setComment("regular clusters collection");
   desc.add<std::string>("folder", "SiStripApproximateClusters")->setComment("Top Level Folder");
   descriptions.addWithDefaultLabel(desc);
 }