#include <CondTools/SiStrip/plugins/SiStripNoisesFromDBMiscalibrator.cc>

Inheritance diagram for SiStripNoisesFromDBMiscalibrator:

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

	~SiStripNoisesFromDBMiscalibrator () override=default

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

	~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

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

std::vector< ESResolverIndex > 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

void	endJob () override

SiStripNoises	getNewObject (const std::map< std::pair< uint32_t, int >, float > &theMap)

SiStripNoises	getNewObject_withDefaults (const std::map< std::pair< uint32_t, int >, float > &theMap, const float theDefault)

Private Attributes
const edm::FileInPath	fp_

const bool	m_fillDefaults

const edm::ESGetToken< SiStripNoises, SiStripNoisesRcd >	m_noiseToken

const std::vector< edm::ParameterSet >	m_parameters

const bool	m_perDetIDdebug

const uint32_t	m_printdebug

const bool	m_saveMaps

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	m_tTopoToken

std::unique_ptr< TrackerMap >	missing_map

std::unique_ptr< TrackerMap >	new_payload_map

std::unique_ptr< TrackerMap >	old_payload_map

std::unique_ptr< TrackerMap >	ratio_map

std::unique_ptr< TrackerMap >	scale_map

std::unique_ptr< TrackerMap >	smear_map

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 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 ()

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

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

Description: Class to miscalibrate a SiStrip Noise payload from Database

Implementation: Read a SiStrip Noise payload from DB (either central DB or sqlite file) and apply a miscalibration (either an offset / gaussian smearing or both) returns a local sqlite file with the same since of the original payload

Definition at line 51 of file SiStripNoisesFromDBMiscalibrator.cc.

Constructor & Destructor Documentation

◆ SiStripNoisesFromDBMiscalibrator()

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

explicit

Definition at line 86 of file SiStripNoisesFromDBMiscalibrator.cc.

References edm::ParameterSet::getUntrackedParameter().

     : m_printdebug{iConfig.getUntrackedParameter<uint32_t>("printDebug", 10)},
       m_perDetIDdebug{iConfig.getUntrackedParameter<bool>("perDetIDdebug", false)},
       m_fillDefaults{iConfig.getUntrackedParameter<bool>("fillDefaults", false)},
       m_saveMaps{iConfig.getUntrackedParameter<bool>("saveMaps", true)},
       m_parameters{iConfig.getParameter<std::vector<edm::ParameterSet> >("params")},
       fp_{iConfig.getUntrackedParameter<edm::FileInPath>("file",
                                                          edm::FileInPath(SiStripDetInfoFileReader::kDefaultFile))},
       m_tTopoToken(esConsumes()),
       m_noiseToken(esConsumes()) {
   //now do what ever initialization is needed
 
   scale_map = std::make_unique<TrackerMap>("scale");
   scale_map->setTitle("Tracker Map of Scale factor averaged by module");
   scale_map->setPalette(1);
 
   smear_map = std::make_unique<TrackerMap>("smear");
   smear_map->setTitle("Tracker Map of Smear factor averaged by module");
   smear_map->setPalette(1);
 
   old_payload_map = std::make_unique<TrackerMap>("old_payload");
   old_payload_map->setTitle("Tracker Map of Starting Noise Payload averaged by module");
   old_payload_map->setPalette(1);
 
   new_payload_map = std::make_unique<TrackerMap>("new_payload");
   new_payload_map->setTitle("Tracker Map of Modified Noise Payload averaged by module");
   new_payload_map->setPalette(1);
 
   ratio_map = std::make_unique<TrackerMap>("ratio");
   ratio_map->setTitle("Tracker Map of Average by module of the payload ratio (new/old)");
   ratio_map->setPalette(1);
 
   if (m_fillDefaults) {
     missing_map = std::make_unique<TrackerMap>("uncabled");
     missing_map->setTitle("Tracker Map of uncabled modules");
     missing_map->setPalette(1);
   }
 }

◆ ~SiStripNoisesFromDBMiscalibrator()

SiStripNoisesFromDBMiscalibrator::~SiStripNoisesFromDBMiscalibrator ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 130 of file SiStripNoisesFromDBMiscalibrator.cc.

References SiStripMiscalibrate::Entry::add(), cond::service::PoolDBOutputService::appendOneIOV(), cond::service::PoolDBOutputService::createOneIOV(), cond::service::PoolDBOutputService::currentTime(), ztail::d, spr::find(), edm::EventSetup::getData(), getNewObject(), getNewObject_withDefaults(), SiStripMiscalibrate::getRegionFromString(), SiStripMiscalibrate::getRegionsFromDetId(), edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmiolumiharvest::j, m_fillDefaults, m_noiseToken, m_parameters, m_perDetIDdebug, m_printdebug, m_tTopoToken, SiStripMiscalibrate::Entry::mean(), new_payload_map, hgchebackDigitizer_cfi::noise, me0TriggerPseudoDigis_cff::nStrips, old_payload_map, submitPVValidationJobs::params, l1ctLayer1_patternWriters_cff::partition, FastTimerService_cff::range, ratio_map, nano_mu_digi_cff::region, SiStripMiscalibrate::Entry::reset(), scale_map, smear_map, jetUpdater_cfi::sort, AlCaHLTBitMon_QueryRunRegistry::string, and parallelization::uint.

                                                                                                   {
   using namespace edm;
 
   const auto tTopo = &iSetup.getData(m_tTopoToken);
 
   std::vector<std::string> partitions;
 
   // fill the list of partitions
   for (auto& thePSet : m_parameters) {
     const std::string partition(thePSet.getParameter<std::string>("partition"));
     // only if it is not yet in the list
     if (std::find(partitions.begin(), partitions.end(), partition) == partitions.end()) {
       partitions.push_back(partition);
     }
   }
 
   std::map<sistripsummary::TrackerRegion, SiStripMiscalibrate::Smearings> mapOfSmearings;
 
   for (auto& thePSet : m_parameters) {
     const std::string partition(thePSet.getParameter<std::string>("partition"));
     sistripsummary::TrackerRegion region = SiStripMiscalibrate::getRegionFromString(partition);
 
     bool m_doScale(thePSet.getParameter<bool>("doScale"));
     bool m_doSmear(thePSet.getParameter<bool>("doSmear"));
     double m_scaleFactor(thePSet.getParameter<double>("scaleFactor"));
     double m_smearFactor(thePSet.getParameter<double>("smearFactor"));
 
     SiStripMiscalibrate::Smearings params = SiStripMiscalibrate::Smearings();
     params.setSmearing(m_doScale, m_doSmear, m_scaleFactor, m_smearFactor);
     mapOfSmearings[region] = params;
   }
 
   const auto& stripNoises = iSetup.getData(m_noiseToken);
 
   std::map<std::pair<uint32_t, int>, float> theMap, oldPayloadMap;
 
   std::vector<uint32_t> detid;
   stripNoises.getDetIds(detid);
   for (const auto& d : detid) {
     SiStripNoises::Range range = stripNoises.getRange(d);
 
     auto regions = SiStripMiscalibrate::getRegionsFromDetId(tTopo, d);
 
     // sort by largest to smallest
     std::sort(regions.rbegin(), regions.rend());
 
     SiStripMiscalibrate::Smearings params = SiStripMiscalibrate::Smearings();
 
     for (unsigned int j = 0; j < regions.size(); j++) {
       bool checkRegion = (mapOfSmearings.count(regions[j]) != 0);
 
       if (!checkRegion) {
         // if the subdetector is not in the list and there's no indication for the whole tracker, just use the default
         // i.e. no change
         continue;
       } else {
         params = mapOfSmearings[regions[j]];
         break;
       }
     }
 
     scale_map->fill(d, params.m_scaleFactor);
     smear_map->fill(d, params.m_smearFactor);
 
     int nStrips = 0;
     for (int it = 0; it < (range.second - range.first) * 8 / 9; ++it) {
       auto noise = stripNoises.getNoise(it, range);
       std::pair<uint32_t, int> index = std::make_pair(d, nStrips);
 
       oldPayloadMap[index] = noise;
 
       if (params.m_doScale) {
         noise *= params.m_scaleFactor;
       }
 
       if (params.m_doSmear) {
         float smearedNoise = CLHEP::RandGauss::shoot(noise, params.m_smearFactor);
         noise = smearedNoise;
       }
 
       theMap[index] = noise;
 
       nStrips += 1;
 
     }  // loop over APVs
   }    // loop over DetIds
 
   SiStripNoises theSiStripNoises{};
   if (!m_fillDefaults) {
     theSiStripNoises = this->getNewObject(theMap);
   } else {
     theSiStripNoises = this->getNewObject_withDefaults(theMap, -1.);
   }
 
   // make the payload ratio map
   uint32_t cachedId(0);
   SiStripMiscalibrate::Entry noise_ratio;
   SiStripMiscalibrate::Entry o_noise;
   SiStripMiscalibrate::Entry n_noise;
   uint countDetIds(0);  // count DetIds to print
   uint countStrips(0);
   for (const auto& element : theMap) {
     countStrips++;
     uint32_t DetId = element.first.first;
     int nstrip = element.first.second;
     float new_noise = element.second;
     float old_noise = oldPayloadMap[std::make_pair(DetId, nstrip)];
 
     // flush the counters
     if (cachedId != 0 && DetId != cachedId) {
       ratio_map->fill(cachedId, noise_ratio.mean());
       old_payload_map->fill(cachedId, o_noise.mean());
       new_payload_map->fill(cachedId, n_noise.mean());
 
       //auto test = new_payload_map.get()->smoduleMap;
 
       noise_ratio.reset();
       o_noise.reset();
       n_noise.reset();
       countDetIds++;
 
       if (m_perDetIDdebug && (countDetIds < m_printdebug)) {
         edm::LogPrint("SiStripNoisesFromDBMiscalibrator")
             << "SiStripNoisesFromDBMiscalibrator"
             << "::" << __FUNCTION__ << " detid " << DetId << " \t"
             << " strip " << nstrip << " \t new <noise>: " << std::setw(5) << std::setprecision(2) << n_noise.mean()
             << " \t old <noise>: " << o_noise.mean() << " \t" << std::endl;
       }
     }
 
     // printout for debug
     if ((countStrips < m_printdebug) && !m_perDetIDdebug) {
       edm::LogPrint("SiStripNoisesFromDBMiscalibrator")
           << "SiStripNoisesFromDBMiscalibrator"
           << "::" << __FUNCTION__ << " detid " << DetId << " \t"
           << " strip " << nstrip << " \t new noise: " << std::setw(5) << std::setprecision(2) << new_noise
           << " \t old noise: " << old_noise << " \t" << std::endl;
     }
 
     cachedId = DetId;
     noise_ratio.add(new_noise / old_noise);
     o_noise.add(old_noise);
     n_noise.add(new_noise);
   }
 
   // write out the SiStripNoises record
   edm::Service<cond::service::PoolDBOutputService> poolDbService;
 
   if (poolDbService.isAvailable()) {
     if (poolDbService->isNewTagRequest("SiStripNoisesRcd")) {
       poolDbService->createOneIOV(theSiStripNoises, poolDbService->currentTime(), "SiStripNoisesRcd");
     } else {
       poolDbService->appendOneIOV(theSiStripNoises, poolDbService->currentTime(), "SiStripNoisesRcd");
     }
   } else {
     throw std::runtime_error("PoolDBService required.");
   }
 }

◆ endJob()

void SiStripNoisesFromDBMiscalibrator::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 290 of file SiStripNoisesFromDBMiscalibrator.cc.

References SiStripMiscalibrate::getTruncatedRange(), m_fillDefaults, m_saveMaps, missing_map, new_payload_map, old_payload_map, FastTimerService_cff::range, ratio_map, scale_map, and smear_map.

                                               {
   if (m_saveMaps) {
     scale_map->save(true, 0, 0, "noise_scale_map.pdf");
     scale_map->save(true, 0, 0, "noise_scale_map.png");
 
     smear_map->save(true, 0, 0, "noise_smear_map.pdf");
     smear_map->save(true, 0, 0, "noise_smear_map.png");
 
     ratio_map->save(true, 0, 0, "noise_ratio_map.pdf");
     ratio_map->save(true, 0, 0, "noise_ratio_map.png");
 
     auto range = SiStripMiscalibrate::getTruncatedRange(old_payload_map.get());
 
     old_payload_map->save(true, range.first, range.second, "starting_noise_payload_map.pdf");
     old_payload_map->save(true, range.first, range.second, "starting_noise_payload_map.png");
 
     range = SiStripMiscalibrate::getTruncatedRange(new_payload_map.get());
 
     new_payload_map->save(true, range.first, range.second, "new_noise_payload_map.pdf");
     new_payload_map->save(true, range.first, range.second, "new_noise_payload_map.png");
 
     if (m_fillDefaults) {
       missing_map->save(true, 0, 0, "missing_map.pdf");
       missing_map->save(true, 0, 0, "missing_map.png");
     }
   }
 }

◆ fillDescriptions()

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

static

Definition at line 406 of file SiStripNoisesFromDBMiscalibrator.cc.

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

                                                                                                   {
   edm::ParameterSetDescription desc;
 
   desc.setComment(
       "Creates rescaled / smeared SiStrip Noise payload."
       "PoolDBOutputService must be set up for 'SiSiStripNoisesRcd'.");
 
   edm::ParameterSetDescription descScaler;
   descScaler.setComment(
       "ParameterSet specifying the Strip tracker partition to be scaled / smeared "
       "by a given factor.");
 
   descScaler.add<std::string>("partition", "Tracker");
   descScaler.add<bool>("doScale", true);
   descScaler.add<bool>("doSmear", true);
   descScaler.add<double>("scaleFactor", 1.0);
   descScaler.add<double>("smearFactor", 1.0);
   desc.addVPSet("params", descScaler, std::vector<edm::ParameterSet>(1));
 
   desc.addUntracked<unsigned int>("printDebug", 10);
   desc.addUntracked<bool>("perDetIDdebug", false);
   desc.addUntracked<bool>("fillDefaults", false);
   desc.addUntracked<bool>("saveMaps", true);
 
   descriptions.add("scaleAndSmearSiStripNoises", desc);
 }

◆ getNewObject()

SiStripNoises SiStripNoisesFromDBMiscalibrator::getNewObject ( const std::map< std::pair< uint32_t, int >, float > & theMap )

private

Definition at line 380 of file SiStripNoisesFromDBMiscalibrator.cc.

References hgchebackDigitizer_cfi::noise, and getGTfromDQMFile::obj.

Referenced by analyze().

                                                                                                                 {
   SiStripNoises obj{};
 
   uint32_t PreviousDetId = 0;
   SiStripNoises::InputVector theSiStripVector;
   for (const auto& element : theMap) {
     uint32_t DetId = element.first.first;
     float noise = element.second;
 
     if (DetId != PreviousDetId) {
       if (!theSiStripVector.empty()) {
         if (!obj.put(PreviousDetId, theSiStripVector)) {
           edm::LogError("SiStripNoisesFromDBMiscalibrator")
               << "[SiStripNoisesFromDBMiscalibrator::analyze] detid already exists" << std::endl;
         }
       }
 
       theSiStripVector.clear();
       PreviousDetId = DetId;
     }
     obj.setData(noise, theSiStripVector);
   }
   return obj;
 }

◆ getNewObject_withDefaults()

SiStripNoises SiStripNoisesFromDBMiscalibrator::getNewObject_withDefaults	(	const std::map< std::pair< uint32_t, int >, float > &	theMap,
		const float	theDefault
	)

private

Definition at line 319 of file SiStripNoisesFromDBMiscalibrator.cc.

References fp_, edm::FileInPath::fullPath(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, LogDebug, combine::missing, missing_map, Skims_PA_cff::name, hgchebackDigitizer_cfi::noise, getGTfromDQMFile::obj, SiStripDetInfoFileReader::read(), DQM::reader, and pileupReCalc_HLTpaths::trunc.

Referenced by analyze().

                                                                                  {
   SiStripNoises obj{};
 
   std::vector<uint32_t> missingDetIds;
 
   const auto& reader = SiStripDetInfoFileReader::read(fp_.fullPath());
   const auto& DetInfos = reader.getAllData();
 
   for (const auto& it : DetInfos) {
     const auto& nAPVs = it.second.nApvs;
     //Generate Noise for det detid
     bool isMissing(false);
     SiStripNoises::InputVector theSiStripVector;
     for (int t_strip = 0; t_strip < 128 * nAPVs; ++t_strip) {
       std::pair<uint32_t, int> index = std::make_pair(it.first, t_strip);
 
       if (theMap.find(index) == theMap.end()) {
         LogDebug("SiStripNoisesFromDBMiscalibrator") << "detid " << it.first << " \t"
                                                      << " strip " << t_strip << " \t"
                                                      << " not found" << std::endl;
 
         isMissing = true;
         obj.setData(theDefault, theSiStripVector);
 
       } else {
         float noise = theMap.at(index);
         obj.setData(noise, theSiStripVector);
       }
     }
 
     if (isMissing)
       missingDetIds.push_back(it.first);
 
     if (!obj.put(it.first, theSiStripVector)) {
       edm::LogError("SiStripNoisesFromDBMiscalibrator")
           << "[SiStripNoisesFromDBMiscalibrator::analyze] detid already exists" << std::endl;
     }
   }
 
   if (!missingDetIds.empty()) {
     // open output file
     std::stringstream name;
     name << "missing_modules.txt";
     std::ofstream* ofile = new std::ofstream(name.str(), std::ofstream::trunc);
     if (!ofile->is_open())
       throw "cannot open output file!";
     for (const auto& missing : missingDetIds) {
       edm::LogVerbatim("SiStripNoisesFromDBMiscalibrator") << missing << "  " << 1 << std::endl;
       (*ofile) << missing << " " << 1 << std::endl;
       missing_map->fill(missing, 1);
     }
 
     ofile->close();
     delete ofile;
   }
 
   return obj;
 }