SiStripQualityHotStripIdentifierRoot Class Reference

#include <SiStripQualityHotStripIdentifierRoot.h>

Inheritance diagram for SiStripQualityHotStripIdentifierRoot:

Public Types
typedef dqm::legacy::DQMStore	DQMStore
typedef dqm::legacy::MonitorElement	MonitorElement
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	SiStripQualityHotStripIdentifierRoot (const edm::ParameterSet &)
	~SiStripQualityHotStripIdentifierRoot () 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::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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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	algoBeginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) override
void	algoBeginRun (const edm::Run &, const edm::EventSetup &) override
void	algoEndJob () override
void	bookHistos ()
std::unique_ptr< SiStripBadStrip >	getNewObject () override

Private Attributes
uint32_t	calibrationthreshold
SiStrip::QualityHistosMap	ClusterPositionHistoMap
const edm::ParameterSet	conf_
std::string	dirpath
DQMStore *	dqmStore_
TFile *	file0
std::string	filename
unsigned short	MaxClusterWidth_
double	MeanNumberOfCluster
unsigned short	MinClusterWidth_
const SiStripQuality *	SiStripQuality_
edm::ESGetToken< SiStripQuality, SiStripQualityRcd >	stripQualityToken_
edm::ESWatcher< SiStripQualityRcd >	stripQualityWatcher_
SiStripHotStripAlgorithmFromClusterOccupancy *	theIdentifier
SiStripBadAPVAlgorithmFromClusterOccupancy *	theIdentifier2
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy *	theIdentifier3
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_
double	TotNumberOfEvents
const TrackerGeometry *	tracker_
const TrackerTopology *	tTopo
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_
bool	UseInputDB_

Additional Inherited Members
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 35 of file SiStripQualityHotStripIdentifierRoot.h.

Member Typedef Documentation

DQMStore

typedef dqm::legacy::DQMStore SiStripQualityHotStripIdentifierRoot::DQMStore

Definition at line 38 of file SiStripQualityHotStripIdentifierRoot.h.

MonitorElement

typedef dqm::legacy::MonitorElement SiStripQualityHotStripIdentifierRoot::MonitorElement

Definition at line 37 of file SiStripQualityHotStripIdentifierRoot.h.

Constructor & Destructor Documentation

SiStripQualityHotStripIdentifierRoot()

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

explicit

Definition at line 25 of file SiStripQualityHotStripIdentifierRoot.cc.

    : ConditionDBWriter<SiStripBadStrip>(iConfig),
      UseInputDB_(iConfig.getUntrackedParameter<bool>("UseInputDB", false)),
      conf_(iConfig),
      tracker_(nullptr),
      tTopo(nullptr),
      filename(iConfig.getUntrackedParameter<std::string>("rootFilename", "CondDB_TKCC_20X_v3_hlt_50822.root")),
      dirpath(iConfig.getUntrackedParameter<std::string>("rootDirPath", "")),
      TotNumberOfEvents(0),
      MeanNumberOfCluster(0),
      calibrationthreshold(iConfig.getUntrackedParameter<uint32_t>("CalibrationThreshold", 10000)),
      tTopoToken_(esConsumes<edm::Transition::BeginRun>()),
      tkGeomToken_(esConsumes<edm::Transition::BeginRun>()),
      stripQualityToken_(UseInputDB_
                             ? decltype(stripQualityToken_)(esConsumes<edm::Transition::BeginRun>(
                                   edm::ESInputTag{"", iConfig.getUntrackedParameter<std::string>("dataLabel", "")}))
                             : decltype(stripQualityToken_){}) {
  dqmStore_ = edm::Service<DQMStore>().operator->();
 
  if (!filename.empty()) {
    edm::LogInfo("SiStripQualityHotStripIdentifierRoot") << " before opening file " << filename.c_str();
    dqmStore_->open(filename, false);
    edm::LogInfo("SiStripQualityHotStripIdentifierRoot") << " after opening file ";
  }
}

References edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

~SiStripQualityHotStripIdentifierRoot()

SiStripQualityHotStripIdentifierRoot::~SiStripQualityHotStripIdentifierRoot ( )

override

Definition at line 51 of file SiStripQualityHotStripIdentifierRoot.cc.

{}

Member Function Documentation

algoAnalyze()

void SiStripQualityHotStripIdentifierRoot::algoAnalyze ( const edm::Event &  , const edm::EventSetup &    )

inlineoverrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 52 of file SiStripQualityHotStripIdentifierRoot.h.

{};

algoBeginLuminosityBlock()

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

inlineoverrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 47 of file SiStripQualityHotStripIdentifierRoot.h.

{}

algoBeginRun()

void SiStripQualityHotStripIdentifierRoot::algoBeginRun ( const edm::Run &  iRun, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 245 of file SiStripQualityHotStripIdentifierRoot.cc.

                                                                                                       {
  tTopo = &iSetup.getData(tTopoToken_);
  tracker_ = &iSetup.getData(tkGeomToken_);
  if (UseInputDB_ && stripQualityWatcher_.check(iSetup)) {
    SiStripQuality_ = &iSetup.getData(stripQualityToken_);
  }
}

References edm::ESWatcher< T >::check(), edm::EventSetup::getData(), SiStripQuality_, stripQualityToken_, stripQualityWatcher_, tkGeomToken_, tracker_, tTopo, tTopoToken_, and UseInputDB_.

algoEndJob()

void SiStripQualityHotStripIdentifierRoot::algoEndJob ( )

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 253 of file SiStripQualityHotStripIdentifierRoot.cc.

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

References ClusterPositionHistoMap.

bookHistos()

void SiStripQualityHotStripIdentifierRoot::bookHistos ( )

private

Definition at line 258 of file SiStripQualityHotStripIdentifierRoot.cc.

                                                      {
  edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
      << " [SiStripQualityHotStripIdentifierRoot::bookHistos] " << dirpath << std::endl;
  std::vector<MonitorElement*> MEs = dqmStore_->getAllContents(dirpath);
  //"DQMData/Run 50908/SiStrip/MechanicalView");
  //                                   "/DQMData/Run 50908/SiStrip/Run summary/MechanicalView/TID/side_2/wheel_3/ring_2/mono_modules/module_402676874");
  edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
      << " [SiStripQualityHotStripIdentifierRoot::bookHistos] vector size " << MEs.size() << std::endl;
 
  std::vector<MonitorElement*>::const_iterator iter = MEs.begin();
  std::vector<MonitorElement*>::const_iterator iterEnd = MEs.end();
 
  edm::ParameterSet parameters = conf_.getParameter<edm::ParameterSet>("AlgoParameters");
  bool gotNentries = true;
  if (parameters.getUntrackedParameter<uint32_t>("NumberOfEvents", 0) == 0 &&
      parameters.getUntrackedParameter<double>("OccupancyThreshold", 0) != 0)
    gotNentries = false;
  edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
      << "[SiStripQualityHotStripIdentifierRoot::bookHistos]  gotNentries flag " << gotNentries << " number of "
      << parameters.getUntrackedParameter<uint32_t>("NumberOfEvents", 0) << " occup "
      << parameters.getUntrackedParameter<double>("OccupancyThreshold", 0) << " OccupancyHisto"
      << parameters.getUntrackedParameter<std::string>("OccupancyHisto") << std::endl;
 
  // Check Number of Events
  for (; iter != iterEnd; ++iter) {
    std::string me_name = (*iter)->getName();
 
    if (!gotNentries && strstr(me_name.c_str(), "TotalNumberOfCluster__T") != nullptr &&
        strstr(me_name.c_str(), "Profile") == nullptr) {
      TotNumberOfEvents = ((TH1F*)(*iter)->getTH1F())->GetEntries();
      MeanNumberOfCluster = ((TH1F*)(*iter)->getTH1F())->GetMean();
      edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
          << "Total Number of Events: " << TotNumberOfEvents << std::endl;
 
      gotNentries = true;
      edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
          << "[SiStripQualityHotStripIdentifierRoot::bookHistos]  gotNentries flag " << gotNentries << std::endl;
      break;
    }
  }
  if (!gotNentries) {
    edm::LogWarning("SiStripQualityHotStripIdentifierRoot")
        << " [SiStripQualityHotStripIdentifierRoot::bookHistos] :: Histogram with to check # of events missing"
        << std::endl;
  }
  for (; iter != iterEnd; ++iter) {
    std::string me_name = (*iter)->getName();
 
    if (strstr(me_name.c_str(), (parameters.getUntrackedParameter<std::string>("OccupancyHisto")).c_str()) == nullptr)
      continue;
 
    unsigned int detid = 0;
    char title[128];
    sprintf(title, "%s", me_name.c_str());
    char* saveptr;
    char* ptr = strtok_r(title, "__", &saveptr);
    int c = 0;
    while (ptr != nullptr) {
      if (c == 2) {
        detid = atol(ptr);
        break;
      }
      ptr = strtok_r(nullptr, "_", &saveptr);
      c++;
    }
    LogDebug("SiStripQualityHotStripIdentifierRoot")
        << " [SiStripQualityHotStripIdentifierRoot::bookHistos] detid " << detid << std::endl;
 
    ClusterPositionHistoMap[detid] = std::make_shared<TH1F>(*(*iter)->getTH1F());
  }
}

References c, ClusterPositionHistoMap, conf_, dirpath, dqmStore_, dqm::implementation::IGetter::getAllContents(), edm::ParameterSet::getParameter(), LogDebug, MeanNumberOfCluster, EcalMEFormatter_cfi::MEs, AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, and TotNumberOfEvents.

Referenced by getNewObject().

getNewObject()

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

overrideprivatevirtual

Implements ConditionDBWriter< SiStripBadStrip >.

Definition at line 53 of file SiStripQualityHotStripIdentifierRoot.cc.

                                                                                  {
  edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
      << "SiStripQualityHotStripIdentifierRoot::getNewObject called" << std::endl;
 
  auto obj = std::make_unique<SiStripBadStrip>();
 
  edm::ParameterSet parameters = conf_.getParameter<edm::ParameterSet>("AlgoParameters");
  std::string AlgoName = parameters.getParameter<std::string>("AlgoName");
 
  bookHistos();
 
  if (TotNumberOfEvents >= calibrationthreshold) {
    edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
        << " [SiStripQualityHotStripIdentifierRoot::getNewObject] Total number of events is " << TotNumberOfEvents
        << "  which is bigger than  " << calibrationthreshold << ". Calibration is preparing to be launched."
        << std::endl;
    if (MeanNumberOfCluster > 0.01) {
      edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
          << "Mean Number of Cluster is " << MeanNumberOfCluster << ". Calibration is launched." << std::endl;
 
      if (AlgoName == "SiStripHotStripAlgorithmFromClusterOccupancy") {
        edm::LogInfo("SiStripQualityHotStripIdentifierRoot") << " [SiStripQualityHotStripIdentifierRoot::getNewObject] "
                                                                "call to SiStripHotStripAlgorithmFromClusterOccupancy"
                                                             << std::endl;
 
        theIdentifier = new SiStripHotStripAlgorithmFromClusterOccupancy(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));
        theIdentifier->setOccupancyThreshold(parameters.getUntrackedParameter<double>("OccupancyThreshold", 1.E-5));
        theIdentifier->setNumberOfEvents(TotNumberOfEvents);
        theIdentifier->setOutputFileName(
            conf_.getUntrackedParameter<std::string>("OccupancyRootFile", "Occupancy.root"),
            conf_.getUntrackedParameter<bool>("WriteOccupancyRootFile", false));
        theIdentifier->setTrackerGeometry(tracker_);
 
        const auto detInfo =
            SiStripDetInfoFileReader::read(edm::FileInPath{SiStripDetInfoFileReader::kDefaultFile}.fullPath());
        SiStripQuality* qobj = new SiStripQuality(detInfo);
        theIdentifier->extractBadStrips(
            qobj,
            ClusterPositionHistoMap,
            SiStripQuality_);  //here I insert SiStripQuality as input and get qobj as output
 
        //----------
 
        edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
            << " [SiStripQualityHotStripIdentifierRoot::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("SiStripQualityHotStripIdentifierRoot")
                << "[SiStripQualityHotStripIdentifierRoot::getNewObject] detid already exists" << std::endl;
        }
        edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
            << " [SiStripQualityHotStripIdentifierRoot::getNewObject] " << ss.str() << std::endl;
      } else if (AlgoName == "SiStripBadAPVAlgorithmFromClusterOccupancy") {
        edm::LogInfo("SiStripQualityHotStripIdentifierRoot") << " [SiStripQualityHotStripIdentifierRoot::getNewObject] "
                                                                "call to SiStripBadAPVAlgorithmFromClusterOccupancy"
                                                             << std::endl;
 
        theIdentifier2 = new SiStripBadAPVAlgorithmFromClusterOccupancy(conf_, tTopo);
        theIdentifier2->setLowOccupancyThreshold(parameters.getUntrackedParameter<double>("LowOccupancyThreshold", 5));
        theIdentifier2->setHighOccupancyThreshold(
            parameters.getUntrackedParameter<double>("HighOccupancyThreshold", 10));
        theIdentifier2->setAbsoluteLowThreshold(parameters.getUntrackedParameter<double>("AbsoluteLowThreshold", 0));
        theIdentifier2->setNumberIterations(parameters.getUntrackedParameter<uint32_t>("NumberIterations", 2));
        theIdentifier2->setAbsoluteOccupancyThreshold(
            parameters.getUntrackedParameter<double>("OccupancyThreshold", 1.E-5));
        theIdentifier2->setNumberOfEvents(TotNumberOfEvents);
        theIdentifier2->setMinNumOfEvents();
        theIdentifier2->setOutputFileName(
            conf_.getUntrackedParameter<std::string>("OccupancyRootFile", "Occupancy.root"),
            conf_.getUntrackedParameter<bool>("WriteOccupancyRootFile", false));
        theIdentifier2->setTrackerGeometry(tracker_);
 
        const auto detInfo =
            SiStripDetInfoFileReader::read(edm::FileInPath{SiStripDetInfoFileReader::kDefaultFile}.fullPath());
        SiStripQuality* qobj = new SiStripQuality(detInfo);
        theIdentifier2->extractBadAPVs(qobj, ClusterPositionHistoMap, SiStripQuality_);
 
        //----------
 
        edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
            << " [SiStripQualityHotStripIdentifierRoot::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("SiStripQualityHotStripIdentifierRoot")
                << "[SiStripQualityHotStripIdentifierRoot::getNewObject] detid already exists" << std::endl;
        }
        edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
            << " [SiStripQualityHotStripIdentifierRoot::getNewObject] " << ss.str() << std::endl;
 
      } else if (AlgoName == "SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy") {
        edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
            << " [SiStripQualityHotStripIdentifierRoot::getNewObject] call to "
               "SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy"
            << std::endl;
 
        theIdentifier3 = new SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy(conf_, tTopo);
        theIdentifier3->setProbabilityThreshold(
            parameters.getUntrackedParameter<double>("ProbabilityThreshold", 1.E-7));
        theIdentifier3->setMinNumEntries(parameters.getUntrackedParameter<uint32_t>("MinNumEntries", 100));
        theIdentifier3->setMinNumEntriesPerStrip(
            parameters.getUntrackedParameter<uint32_t>("MinNumEntriesPerStrip", 5));
        theIdentifier3->setNumberOfEvents(TotNumberOfEvents);
        theIdentifier3->setOutputFileName(
            conf_.getUntrackedParameter<std::string>("OccupancyRootFile", "Occupancy.root"),
            conf_.getUntrackedParameter<bool>("WriteOccupancyRootFile", false),
            conf_.getUntrackedParameter<std::string>("DQMHistoOutputFile", "DQMHistos.root"),
            conf_.getUntrackedParameter<bool>("WriteDQMHistoOutputFile", false));
        theIdentifier3->setTrackerGeometry(tracker_);
        theIdentifier3->setLowOccupancyThreshold(parameters.getUntrackedParameter<double>("LowOccupancyThreshold", 5));
        theIdentifier3->setHighOccupancyThreshold(
            parameters.getUntrackedParameter<double>("HighOccupancyThreshold", 10));
        theIdentifier3->setAbsoluteLowThreshold(parameters.getUntrackedParameter<double>("AbsoluteLowThreshold", 0));
        theIdentifier3->setNumberIterations(parameters.getUntrackedParameter<uint32_t>("NumberIterations", 2));
        theIdentifier3->setAbsoluteOccupancyThreshold(
            parameters.getUntrackedParameter<double>("OccupancyThreshold", 1.E-5));
        theIdentifier3->setMinNumOfEvents();
 
        const auto detInfo =
            SiStripDetInfoFileReader::read(edm::FileInPath{SiStripDetInfoFileReader::kDefaultFile}.fullPath());
        SiStripQuality* qobj = new SiStripQuality(detInfo);
        theIdentifier3->extractBadAPVSandStrips(
            qobj,
            ClusterPositionHistoMap,
            SiStripQuality_);  //here I insert SiStripQuality as input and get qobj as output
 
        //----------
 
        edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
            << " [SiStripQualityHotStripIdentifierRoot::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("SiStripQualityHotStripIdentifierRoot")
                << "[SiStripQualityHotStripIdentifierRoot::getNewObject] detid already exists" << std::endl;
        }
        edm::LogInfo("SiStripQualityHotStripIdentifierRoot")
            << " [SiStripQualityHotStripIdentifierRoot::getNewObject] " << ss.str() << std::endl;
      } else {
        edm::LogError("SiStripQualityHotStripIdentifierRoot") << " [SiStripQualityHotStripIdentifierRoot::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("SiStripQualityHotStripIdentifierRoot")
              << "[SiStripQualityHotStripIdentifierRoot::getNewObject] detid already exists" << std::endl;
      }
    }
 
    else {
      edm::LogWarning("SiStripQualityHotStripIdentifierRoot")
          << " [SiStripQualityHotStripIdentifierRoot::getNewObject] :: Empty Events :: "
          << "Mean Number of Cluster is  " << MeanNumberOfCluster << " Calibration NOT launched" << std::endl;
      setDoStore(false);  // Don't put anything in the sqlite-file!
    }
  } else {
    edm::LogWarning("SiStripQualityHotStripIdentifierRoot")
        << " [SiStripQualityHotStripIdentifierRoot::getNewObject] :: Not Enough Events :: "
        << "Total number of events is  " << TotNumberOfEvents << " Calibration NOT launched" << std::endl;
    setDoStore(false);  // Don't put anything in the sqlite-file!
  }
 
  return obj;
}

References a, HLTSiStripMonitoring_cff::AlgoName, bookHistos(), calibrationthreshold, ClusterPositionHistoMap, conf_, SiStripBadAPVAlgorithmFromClusterOccupancy::extractBadAPVs(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::extractBadAPVSandStrips(), SiStripHotStripAlgorithmFromClusterOccupancy::extractBadStrips(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), SiStripDetInfoFileReader::kDefaultFile, MeanNumberOfCluster, getGTfromDQMFile::obj, FastTimerService_cff::range, SiStripDetInfoFileReader::read(), SiStripBadAPVAlgorithmFromClusterOccupancy::setAbsoluteLowThreshold(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setAbsoluteLowThreshold(), SiStripBadAPVAlgorithmFromClusterOccupancy::setAbsoluteOccupancyThreshold(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setAbsoluteOccupancyThreshold(), ConditionDBWriter< SiStripBadStrip >::setDoStore(), SiStripBadAPVAlgorithmFromClusterOccupancy::setHighOccupancyThreshold(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setHighOccupancyThreshold(), SiStripBadAPVAlgorithmFromClusterOccupancy::setLowOccupancyThreshold(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setLowOccupancyThreshold(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setMinNumEntries(), SiStripHotStripAlgorithmFromClusterOccupancy::setMinNumEntries(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setMinNumEntriesPerStrip(), SiStripHotStripAlgorithmFromClusterOccupancy::setMinNumEntriesPerStrip(), SiStripBadAPVAlgorithmFromClusterOccupancy::setMinNumOfEvents(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setMinNumOfEvents(), SiStripBadAPVAlgorithmFromClusterOccupancy::setNumberIterations(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setNumberIterations(), SiStripBadAPVAlgorithmFromClusterOccupancy::setNumberOfEvents(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setNumberOfEvents(), SiStripHotStripAlgorithmFromClusterOccupancy::setNumberOfEvents(), SiStripHotStripAlgorithmFromClusterOccupancy::setOccupancyThreshold(), SiStripBadAPVAlgorithmFromClusterOccupancy::setOutputFileName(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setOutputFileName(), SiStripHotStripAlgorithmFromClusterOccupancy::setOutputFileName(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setProbabilityThreshold(), SiStripHotStripAlgorithmFromClusterOccupancy::setProbabilityThreshold(), SiStripBadAPVAlgorithmFromClusterOccupancy::setTrackerGeometry(), SiStripHotStripAlgorithmFromClusterOccupancy::setTrackerGeometry(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setTrackerGeometry(), AlCaHarvesting_cff::SiStripQuality, SiStripQuality_, contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, theIdentifier, theIdentifier2, theIdentifier3, TotNumberOfEvents, tracker_, and tTopo.

Member Data Documentation

calibrationthreshold

uint32_t SiStripQualityHotStripIdentifierRoot::calibrationthreshold

private

Definition at line 73 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by getNewObject().

ClusterPositionHistoMap

SiStrip::QualityHistosMap SiStripQualityHotStripIdentifierRoot::ClusterPositionHistoMap

private

Definition at line 80 of file SiStripQualityHotStripIdentifierRoot.h.

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

conf_

const edm::ParameterSet SiStripQualityHotStripIdentifierRoot::conf_

private

Definition at line 60 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by bookHistos(), and getNewObject().

dirpath

std::string SiStripQualityHotStripIdentifierRoot::dirpath

private

Definition at line 69 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by bookHistos().

dqmStore_

DQMStore* SiStripQualityHotStripIdentifierRoot::dqmStore_

private

Definition at line 66 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by bookHistos().

file0

TFile* SiStripQualityHotStripIdentifierRoot::file0

private

Definition at line 68 of file SiStripQualityHotStripIdentifierRoot.h.

filename

std::string SiStripQualityHotStripIdentifierRoot::filename

private

Definition at line 69 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by python.rootplot.rootmath.Target::__repr__(), and utils.unpickler::run().

MaxClusterWidth_

unsigned short SiStripQualityHotStripIdentifierRoot::MaxClusterWidth_

private

Definition at line 70 of file SiStripQualityHotStripIdentifierRoot.h.

MeanNumberOfCluster

double SiStripQualityHotStripIdentifierRoot::MeanNumberOfCluster

private

Definition at line 72 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by bookHistos(), and getNewObject().

MinClusterWidth_

unsigned short SiStripQualityHotStripIdentifierRoot::MinClusterWidth_

private

Definition at line 70 of file SiStripQualityHotStripIdentifierRoot.h.

SiStripQuality_

const SiStripQuality* SiStripQualityHotStripIdentifierRoot::SiStripQuality_

private

Definition at line 64 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun(), and getNewObject().

stripQualityToken_

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

private

Definition at line 77 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun().

stripQualityWatcher_

edm::ESWatcher<SiStripQualityRcd> SiStripQualityHotStripIdentifierRoot::stripQualityWatcher_

private

Definition at line 78 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun().

theIdentifier

SiStripHotStripAlgorithmFromClusterOccupancy* SiStripQualityHotStripIdentifierRoot::theIdentifier

private

Definition at line 81 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by getNewObject().

theIdentifier2

SiStripBadAPVAlgorithmFromClusterOccupancy* SiStripQualityHotStripIdentifierRoot::theIdentifier2

private

Definition at line 82 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by getNewObject().

theIdentifier3

SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy* SiStripQualityHotStripIdentifierRoot::theIdentifier3

private

Definition at line 83 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by getNewObject().

tkGeomToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiStripQualityHotStripIdentifierRoot::tkGeomToken_

private

Definition at line 76 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun().

TotNumberOfEvents

double SiStripQualityHotStripIdentifierRoot::TotNumberOfEvents

private

Definition at line 71 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by bookHistos(), and getNewObject().

tracker_

const TrackerGeometry* SiStripQualityHotStripIdentifierRoot::tracker_

private

Definition at line 62 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun(), and getNewObject().

tTopo

const TrackerTopology* SiStripQualityHotStripIdentifierRoot::tTopo

private

Definition at line 63 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun(), and getNewObject().

tTopoToken_

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

private

Definition at line 75 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun().

UseInputDB_

bool SiStripQualityHotStripIdentifierRoot::UseInputDB_

private

Definition at line 59 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun().