#include <SiStripQualityHotStripIdentifier.h>

Inheritance diagram for SiStripQualityHotStripIdentifier:

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

	~SiStripQualityHotStripIdentifier () override

Public Member Functions inherited from ConditionDBWriter< SiStripBadStrip >
	ConditionDBWriter (const edm::ParameterSet &iConfig)

	~ConditionDBWriter () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	algoAnalyze (const edm::Event &, const edm::EventSetup &) override

void	algoBeginJob (const edm::EventSetup &) override

void	algoBeginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) override

void	algoBeginRun (const edm::Run &, const edm::EventSetup &) override

void	algoEndJob () override

void	bookHistos ()

void	fillHisto (uint32_t detid, float value)

std::unique_ptr< SiStripBadStrip >	getNewObject () override

void	resetHistos ()

Private Attributes
edm::InputTag	Cluster_src_

SiStrip::QualityHistosMap	ClusterPositionHistoMap

const edm::ParameterSet	conf_

std::string	dataLabel_

edm::FileInPath	fp_

unsigned short	MaxClusterWidth_

unsigned short	MinClusterWidth_

SiStripDetInfoFileReader *	reader

const SiStripQuality *	stripQuality_ = nullptr

edm::ESGetToken< SiStripQuality, SiStripQualityRcd >	stripQualityToken_

edm::ESWatcher< SiStripQualityRcd >	stripQualityWatcher_

edm::InputTag	Track_src_

bool	tracksCollection_in_EventTree

const TrackerTopology *	tTopo = nullptr

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from ConditionDBWriter< SiStripBadStrip >
void	setDoStore (const bool doStore)
	When set to false the payload will not be written to the db. More...

void	storeOnDbNow ()

cond::Time_t	timeOfLastIOV ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 24 of file SiStripQualityHotStripIdentifier.h.

Constructor & Destructor Documentation

◆ SiStripQualityHotStripIdentifier()

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

explicit

Definition at line 19 of file SiStripQualityHotStripIdentifier.cc.

     : ConditionDBWriter<SiStripBadStrip>(iConfig),
       dataLabel_(iConfig.getUntrackedParameter<std::string>("dataLabel", "")),
       conf_(iConfig),
       fp_(iConfig.getUntrackedParameter<edm::FileInPath>(
           "file", edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat"))),
       Cluster_src_(iConfig.getParameter<edm::InputTag>("Cluster_src")),
       Track_src_(iConfig.getUntrackedParameter<edm::InputTag>("Track_src")),
       tracksCollection_in_EventTree(iConfig.getUntrackedParameter<bool>("RemoveTrackClusters", false)),
       tTopoToken_(esConsumes<edm::Transition::BeginRun>()),
       stripQualityToken_(esConsumes<edm::Transition::BeginRun>()) {
   reader = new SiStripDetInfoFileReader(fp_.fullPath());
  
   edm::ParameterSet pset = iConfig.getUntrackedParameter<edm::ParameterSet>("ClusterSelection", edm::ParameterSet());
   MinClusterWidth_ = pset.getUntrackedParameter<uint32_t>("minWidth", 1);
   MaxClusterWidth_ = pset.getUntrackedParameter<uint32_t>("maxWidth", 1000);
  
   bookHistos();
 }

References bookHistos(), fp_, edm::FileInPath::fullPath(), edm::ParameterSet::getUntrackedParameter(), MaxClusterWidth_, MinClusterWidth_, muonDTDigis_cfi::pset, and reader.

◆ ~SiStripQualityHotStripIdentifier()

SiStripQualityHotStripIdentifier::~SiStripQualityHotStripIdentifier ( )

override

Definition at line 39 of file SiStripQualityHotStripIdentifier.cc.

39 {}

Member Function Documentation

◆ algoAnalyze()

void SiStripQualityHotStripIdentifier::algoAnalyze	(	const edm::Event &	e,
		const edm::EventSetup &	eSetup
	)

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 143 of file SiStripQualityHotStripIdentifier.cc.

                                                                                                  {
   edm::Handle<edm::DetSetVector<SiStripCluster> > dsv_SiStripCluster;
   e.getByLabel(Cluster_src_, dsv_SiStripCluster);
  
   edm::Handle<reco::TrackCollection> trackCollection;
   if (tracksCollection_in_EventTree) {
     e.getByLabel(Track_src_, trackCollection);
     if (!trackCollection.isValid()) {
       edm::LogError("SiStripQualityHotStripIdentifier")
           << " [SiStripQualityHotStripIdentifier::algoAnalyze] missing trackCollection with label " << Track_src_
           << std::endl;
     }
   }
  
   std::set<const void*> vPSiStripCluster;
   //Perform track study
   if (tracksCollection_in_EventTree) {
     const reco::TrackCollection tC = *(trackCollection.product());
     int i = 0;
     for (reco::TrackCollection::const_iterator track = tC.begin(); track != tC.end(); track++) {
       LogTrace("SiStripQualityHotStripIdentifier")
           << "Track number " << i + 1 << "\n\tmomentum: " << track->momentum() << "\n\tPT: " << track->pt()
           << "\n\tvertex: " << track->vertex() << "\n\timpact parameter: " << track->d0()
           << "\n\tcharge: " << track->charge() << "\n\tnormalizedChi2: " << track->normalizedChi2()
           << "\n\tFrom EXTRA : "
           << "\n\t\touter PT " << track->outerPt() << std::endl;
  
       //Loop on rechits
       for (auto const& recHit : track->recHits()) {
         if (!recHit->isValid()) {
           LogTrace("SiStripQualityHotStripIdentifier") << "\t\t Invalid Hit " << std::endl;
           continue;
         }
  
         const SiStripRecHit2D* singleHit = dynamic_cast<const SiStripRecHit2D*>(recHit);
         const SiStripMatchedRecHit2D* matchedHit = dynamic_cast<const SiStripMatchedRecHit2D*>(recHit);
         const ProjectedSiStripRecHit2D* projectedHit = dynamic_cast<const ProjectedSiStripRecHit2D*>(recHit);
  
         if (matchedHit) {
           vPSiStripCluster.insert((void*)&(matchedHit->monoCluster()));
           vPSiStripCluster.insert((void*)&(matchedHit->stereoCluster()));
         } else if (projectedHit) {
           vPSiStripCluster.insert((void*)&*(projectedHit->originalHit().cluster()));
         } else if (singleHit) {
           vPSiStripCluster.insert((void*)&*(singleHit->cluster()));
         } else {
           LogTrace("SiStripQualityHotStripIdentifier") << "NULL hit" << std::endl;
         }
       }
     }
   }
  
   std::stringstream ss;
   //Loop on Det Clusters
   edm::DetSetVector<SiStripCluster>::const_iterator DSViter = dsv_SiStripCluster->begin();
   for (; DSViter != dsv_SiStripCluster->end(); DSViter++) {
     edm::DetSet<SiStripCluster>::const_iterator ClusIter = DSViter->data.begin();
     edm::DetSet<SiStripCluster>::const_iterator ClusIterEnd = DSViter->data.end();
     for (; ClusIter != ClusIterEnd; ++ClusIter) {
       if (MinClusterWidth_ <= ClusIter->amplitudes().size() && ClusIter->amplitudes().size() <= MaxClusterWidth_) {
         if (std::find(vPSiStripCluster.begin(), vPSiStripCluster.end(), (void*)&*ClusIter) == vPSiStripCluster.end()) {
           if (edm::isDebugEnabled())
             ss << " adding cluster to histo for detid " << DSViter->id << " with barycenter " << ClusIter->barycenter()
                << std::endl;
           fillHisto(DSViter->id, ClusIter->barycenter());
         }
       }
     }
   }
   LogTrace("SiStripQualityHotStripIdentifier") << ss.str();
 }

References edm::DetSetVector< T >::begin(), SiStripRecHit2D::cluster(), Cluster_src_, MillePedeFileConverter_cfg::e, edm::DetSetVector< T >::end(), fillHisto(), spr::find(), mps_fire::i, edm::isDebugEnabled(), LogTrace, MaxClusterWidth_, SiStripMatchedRecHit2D::monoCluster(), ProjectedSiStripRecHit2D::originalHit(), rpcPointValidation_cfi::recHit, findQualityFiles::size, contentValuesCheck::ss, SiStripMatchedRecHit2D::stereoCluster(), HLT_FULL_cff::track, Track_src_, duplicaterechits_cfi::trackCollection, and tracksCollection_in_EventTree.

◆ algoBeginJob()

void SiStripQualityHotStripIdentifier::algoBeginJob ( const edm::EventSetup & )

inlineoverrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 31 of file SiStripQualityHotStripIdentifier.h.

31 {}

◆ algoBeginLuminosityBlock()

void SiStripQualityHotStripIdentifier::algoBeginLuminosityBlock	(	const edm::LuminosityBlock &	,
		const edm::EventSetup &
	)

inlineoverrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 35 of file SiStripQualityHotStripIdentifier.h.

35 { resetHistos(); }

References resetHistos().

◆ algoBeginRun()

void SiStripQualityHotStripIdentifier::algoBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 91 of file SiStripQualityHotStripIdentifier.cc.

                                                                                                   {
   tTopo = &iSetup.getData(tTopoToken_);
  
   resetHistos();
  
   if (stripQualityWatcher_.check(iSetup)) {
     stripQuality_ = &iSetup.getData(stripQualityToken_);
   }
 }

References edm::ESWatcher< T >::check(), edm::EventSetup::getData(), resetHistos(), stripQuality_, stripQualityToken_, stripQualityWatcher_, tTopo, and tTopoToken_.

◆ algoEndJob()

void SiStripQualityHotStripIdentifier::algoEndJob ( )

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 101 of file SiStripQualityHotStripIdentifier.cc.

                                                   {
   //Clear map
   ClusterPositionHistoMap.clear();
 }

References ClusterPositionHistoMap.

◆ bookHistos()

void SiStripQualityHotStripIdentifier::bookHistos ( )

private

Definition at line 115 of file SiStripQualityHotStripIdentifier.cc.

                                                   {
   edm::LogInfo("SiStripQualityHotStripIdentifier") << " [SiStripQualityHotStripIdentifier::bookHistos] " << std::endl;
   char hname[1024];
   std::map<uint32_t, SiStripDetInfoFileReader::DetInfo>::const_iterator it = reader->getAllData().begin();
   std::map<uint32_t, SiStripDetInfoFileReader::DetInfo>::const_iterator iEnd = reader->getAllData().end();
   for (; it != iEnd; ++it) {
     sprintf(hname, "h_%d", it->first);
     SiStrip::QualityHistosMap::iterator ref = ClusterPositionHistoMap.find(it->first);
     if (ref == ClusterPositionHistoMap.end()) {
       ClusterPositionHistoMap[it->first] =
           std::make_shared<TH1F>(hname, hname, it->second.nApvs * 128, -0.5, it->second.nApvs * 128 - 0.5);
     } else
       edm::LogError("SiStripQualityHotStripIdentifier")
           << " [SiStripQualityHotStripIdentifier::bookHistos] DetId " << it->first
           << " already found in map. Ignoring new data" << std::endl;
   }
 }

References ClusterPositionHistoMap, SiStripDetInfoFileReader::getAllData(), and reader.

Referenced by SiStripQualityHotStripIdentifier().

◆ fillHisto()

void SiStripQualityHotStripIdentifier::fillHisto	(	uint32_t	detid,
		float	value
	)

private

Definition at line 133 of file SiStripQualityHotStripIdentifier.cc.

                                                                             {
   SiStrip::QualityHistosMap::iterator ref = ClusterPositionHistoMap.find(detid);
   if (ref != ClusterPositionHistoMap.end())
     ref->second->Fill(value);
   else
     edm::LogError("SiStripQualityHotStripIdentifier")
         << " [SiStripQualityHotStripIdentifier::fillHisto] Histogram not found in the list for DetId " << detid
         << " Ignoring data value " << value << std::endl;
 }

References ClusterPositionHistoMap.

Referenced by algoAnalyze().

◆ getNewObject()

std::unique_ptr< SiStripBadStrip > SiStripQualityHotStripIdentifier::getNewObject ( )

overrideprivatevirtual

Implements ConditionDBWriter< SiStripBadStrip >.

Definition at line 41 of file SiStripQualityHotStripIdentifier.cc.

                                                                               {
   auto obj = std::make_unique<SiStripBadStrip>();
  
   edm::ParameterSet parameters = conf_.getParameter<edm::ParameterSet>("AlgoParameters");
   std::string AlgoName = parameters.getParameter<std::string>("AlgoName");
   if (AlgoName == "SiStripHotStripAlgorithmFromClusterOccupancy") {
     edm::LogInfo("SiStripQualityHotStripIdentifier")
         << " [SiStripQualityHotStripIdentifier::getNewObject] call to SiStripHotStripAlgorithmFromClusterOccupancy"
         << std::endl;
  
     SiStripHotStripAlgorithmFromClusterOccupancy theIdentifier(conf_, tTopo);
     theIdentifier.setProbabilityThreshold(parameters.getUntrackedParameter<double>("ProbabilityThreshold", 1.E-7));
     theIdentifier.setMinNumEntries(parameters.getUntrackedParameter<uint32_t>("MinNumEntries", 100));
     theIdentifier.setMinNumEntriesPerStrip(parameters.getUntrackedParameter<uint32_t>("MinNumEntriesPerStrip", 5));
  
     SiStripQuality* qobj = new SiStripQuality();
     theIdentifier.extractBadStrips(qobj, ClusterPositionHistoMap, stripQuality_);
  
     edm::LogInfo("SiStripQualityHotStripIdentifier")
         << " [SiStripQualityHotStripIdentifier::getNewObject] copy SiStripObject in SiStripBadStrip" << std::endl;
  
     std::stringstream ss;
  
     SiStripBadStrip::RegistryIterator rIter = qobj->getRegistryVectorBegin();
     SiStripBadStrip::RegistryIterator rIterEnd = qobj->getRegistryVectorEnd();
     for (; rIter != rIterEnd; ++rIter) {
       SiStripBadStrip::Range range(qobj->getDataVectorBegin() + rIter->ibegin,
                                    qobj->getDataVectorBegin() + rIter->iend);
       if (!obj->put(rIter->detid, range))
         edm::LogError("SiStripQualityHotStripIdentifier")
             << "[SiStripQualityHotStripIdentifier::getNewObject] detid already exists" << std::endl;
     }
     edm::LogInfo("SiStripQualityHotStripIdentifier")
         << " [SiStripQualityHotStripIdentifier::getNewObject] " << ss.str() << std::endl;
  
   } else {
     edm::LogError("SiStripQualityHotStripIdentifier")
         << " [SiStripQualityHotStripIdentifier::getNewObject] call for a unknow HotStrip identification algoritm"
         << std::endl;
  
     std::vector<uint32_t> a;
     SiStripBadStrip::Range range(a.begin(), a.end());
     if (!obj->put(0xFFFFFFFF, range))
       edm::LogError("SiStripQualityHotStripIdentifier")
           << "[SiStripQualityHotStripIdentifier::getNewObject] detid already exists" << std::endl;
   }
  
   return obj;
 }

References a, HLTSiStripMonitoring_cff::AlgoName, ClusterPositionHistoMap, conf_, SiStripHotStripAlgorithmFromClusterOccupancy::extractBadStrips(), SiStripBadStrip::getDataVectorBegin(), edm::ParameterSet::getParameter(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), getGTfromDQMFile::obj, FastTimerService_cff::range, SiStripHotStripAlgorithmFromClusterOccupancy::setMinNumEntries(), SiStripHotStripAlgorithmFromClusterOccupancy::setMinNumEntriesPerStrip(), SiStripHotStripAlgorithmFromClusterOccupancy::setProbabilityThreshold(), AlCaHarvesting_cff::SiStripQuality, contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, stripQuality_, and tTopo.

◆ resetHistos()

void SiStripQualityHotStripIdentifier::resetHistos ( )

private

Definition at line 106 of file SiStripQualityHotStripIdentifier.cc.

                                                    {
   edm::LogInfo("SiStripQualityHotStripIdentifier") << " [SiStripQualityHotStripIdentifier::resetHistos] " << std::endl;
   SiStrip::QualityHistosMap::iterator it = ClusterPositionHistoMap.begin();
   SiStrip::QualityHistosMap::iterator iEnd = ClusterPositionHistoMap.end();
   for (; it != iEnd; ++it) {
     it->second->Reset();
   }
 }