SiStripQualityHotStripIdentifier Class Reference

#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::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () 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 const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
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::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
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_
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::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
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<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
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...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 23 of file SiStripQualityHotStripIdentifier.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 20 of file SiStripQualityHotStripIdentifier.cc.

References bookHistos(), edm::ParameterSet::getUntrackedParameter(), MaxClusterWidth_, MinClusterWidth_, and TrackValidation_cff::pset.

    : ConditionDBWriter<SiStripBadStrip>(iConfig),
      dataLabel_(iConfig.getUntrackedParameter<std::string>("dataLabel", "")),
      conf_(iConfig),
      fp_(iConfig.getUntrackedParameter<edm::FileInPath>("file",
                                                         edm::FileInPath(SiStripDetInfoFileReader::kDefaultFile))),
      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>()) {
  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();
}

SiStripQualityHotStripIdentifier::~SiStripQualityHotStripIdentifier ( )

override

Definition at line 38 of file SiStripQualityHotStripIdentifier.cc.

{}

Member Function Documentation

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

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 140 of file SiStripQualityHotStripIdentifier.cc.

References SiStripRecHit2D::cluster(), Cluster_src_, fillHisto(), spr::find(), edm::Event::getByLabel(), mps_fire::i, edm::isDebugEnabled(), edm::HandleBase::isValid(), LogTrace, MaxClusterWidth_, MinClusterWidth_, SiStripMatchedRecHit2D::monoCluster(), ProjectedSiStripRecHit2D::originalHit(), edm::Handle< T >::product(), contentValuesCheck::ss, SiStripMatchedRecHit2D::stereoCluster(), HLT_FULL_cff::track, Track_src_, HLT_FULL_cff::trackCollection, and tracksCollection_in_EventTree.

                                                                                                 {
  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) {
    int i = 0;
    for (const auto& track : *(trackCollection.product())) {
      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
  for (const auto& dSet : *dsv_SiStripCluster) {
    for (const auto& clus : dSet.data) {
      if (MinClusterWidth_ <= clus.amplitudes().size() && clus.amplitudes().size() <= MaxClusterWidth_) {
        if (std::find(vPSiStripCluster.begin(), vPSiStripCluster.end(), (void*)&clus) == vPSiStripCluster.end()) {
          if (edm::isDebugEnabled())
            ss << " adding cluster to histo for detid " << dSet.id << " with barycenter " << clus.barycenter()
               << std::endl;
          fillHisto(dSet.id, clus.barycenter());
        }
      }
    }
  }
  LogTrace("SiStripQualityHotStripIdentifier") << ss.str();
}

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

inlineoverrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 30 of file SiStripQualityHotStripIdentifier.h.

{}

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

inlineoverrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 34 of file SiStripQualityHotStripIdentifier.h.

References resetHistos().

{ resetHistos(); }

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

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 92 of file SiStripQualityHotStripIdentifier.cc.

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

                                                                                                  {
  tTopo = &iSetup.getData(tTopoToken_);

  resetHistos();

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

void SiStripQualityHotStripIdentifier::algoEndJob ( )

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 102 of file SiStripQualityHotStripIdentifier.cc.

References ClusterPositionHistoMap.

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

void SiStripQualityHotStripIdentifier::bookHistos ( )

private

Definition at line 114 of file SiStripQualityHotStripIdentifier.cc.

References ClusterPositionHistoMap, fp_, edm::FileInPath::fullPath(), and SiStripDetInfoFileReader::read().

Referenced by SiStripQualityHotStripIdentifier().

                                                  {
  edm::LogInfo("SiStripQualityHotStripIdentifier") << " [SiStripQualityHotStripIdentifier::bookHistos] " << std::endl;
  char hname[1024];
  for (const auto& it : SiStripDetInfoFileReader::read(fp_.fullPath()).getAllData()) {
    sprintf(hname, "h_%d", it.first);
    auto 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;
  }
}

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

private

Definition at line 130 of file SiStripQualityHotStripIdentifier.cc.

References ClusterPositionHistoMap.

Referenced by algoAnalyze().

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

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

overrideprivatevirtual

Implements ConditionDBWriter< SiStripBadStrip >.

Definition at line 40 of file SiStripQualityHotStripIdentifier.cc.

References a, HLT_FULL_cff::AlgoName, ClusterPositionHistoMap, conf_, SiStripHotStripAlgorithmFromClusterOccupancy::extractBadStrips(), edm::ParameterSet::getParameter(), SiStripDetInfoFileReader::kDefaultFile, getGTfromDQMFile::obj, sistrip::SpyUtilities::range(), SiStripDetInfoFileReader::read(), SiStripHotStripAlgorithmFromClusterOccupancy::setMinNumEntries(), SiStripHotStripAlgorithmFromClusterOccupancy::setMinNumEntriesPerStrip(), SiStripHotStripAlgorithmFromClusterOccupancy::setProbabilityThreshold(), AlCaHarvesting_cff::SiStripQuality, contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, stripQuality_, and tTopo.

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

    const auto detInfo =
        SiStripDetInfoFileReader::read(edm::FileInPath{SiStripDetInfoFileReader::kDefaultFile}.fullPath());
    SiStripQuality* qobj = new SiStripQuality(detInfo);
    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;
}

void SiStripQualityHotStripIdentifier::resetHistos ( )

private

Definition at line 107 of file SiStripQualityHotStripIdentifier.cc.

References ClusterPositionHistoMap.

Referenced by algoBeginLuminosityBlock(), and algoBeginRun().

                                                   {
  edm::LogInfo("SiStripQualityHotStripIdentifier") << " [SiStripQualityHotStripIdentifier::resetHistos] " << std::endl;
  for (const auto& it : ClusterPositionHistoMap) {
    it.second->Reset();
  }
}

Member Data Documentation

edm::InputTag SiStripQualityHotStripIdentifier::Cluster_src_

private

Definition at line 52 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoAnalyze().

SiStrip::QualityHistosMap SiStripQualityHotStripIdentifier::ClusterPositionHistoMap

private

Definition at line 59 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoEndJob(), bookHistos(), fillHisto(), getNewObject(), and resetHistos().

const edm::ParameterSet SiStripQualityHotStripIdentifier::conf_

private

Definition at line 50 of file SiStripQualityHotStripIdentifier.h.

Referenced by getNewObject().

std::string SiStripQualityHotStripIdentifier::dataLabel_

private

Definition at line 48 of file SiStripQualityHotStripIdentifier.h.

edm::FileInPath SiStripQualityHotStripIdentifier::fp_

private

Definition at line 51 of file SiStripQualityHotStripIdentifier.h.

Referenced by bookHistos().

unsigned short SiStripQualityHotStripIdentifier::MaxClusterWidth_

private

Definition at line 57 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoAnalyze(), and SiStripQualityHotStripIdentifier().

unsigned short SiStripQualityHotStripIdentifier::MinClusterWidth_

private

Definition at line 57 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoAnalyze(), and SiStripQualityHotStripIdentifier().

const SiStripQuality* SiStripQualityHotStripIdentifier::stripQuality_ = nullptr

private

Definition at line 49 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoBeginRun(), and getNewObject().

edm::ESGetToken<SiStripQuality, SiStripQualityRcd> SiStripQualityHotStripIdentifier::stripQualityToken_

private

Definition at line 62 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoBeginRun().

edm::ESWatcher<SiStripQualityRcd> SiStripQualityHotStripIdentifier::stripQualityWatcher_

private

Definition at line 63 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoBeginRun().

edm::InputTag SiStripQualityHotStripIdentifier::Track_src_

private

Definition at line 53 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoAnalyze().

bool SiStripQualityHotStripIdentifier::tracksCollection_in_EventTree

private

Definition at line 54 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoAnalyze().

const TrackerTopology* SiStripQualityHotStripIdentifier::tTopo = nullptr

private

Definition at line 55 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoBeginRun(), and getNewObject().

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiStripQualityHotStripIdentifier::tTopoToken_

private

Definition at line 61 of file SiStripQualityHotStripIdentifier.h.

Referenced by algoBeginRun().