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

Inheritance diagram for SiStripApvGainRescaler:

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

	~SiStripApvGainRescaler () override

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

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

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

Private Attributes
const edm::ESGetToken< SiStripGain, SiStripGainRcd >	g1g2Token_

const edm::ESGetToken< SiStripApvGain, SiStripApvGain3Rcd >	g3Token_

const uint32_t	m_printdebug

const std::string	m_Record

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: Utility class to rescale the values of SiStrip G2 by the ratio of G1_old/G1_new: this is useful in the case in which a Gain2 payload needs to recycled after a G1 update to keep the G1*G2 product constant

Implementation: [Notes on implementation]

Definition at line 45 of file SiStripApvGainRescaler.cc.

Constructor & Destructor Documentation

◆ SiStripApvGainRescaler()

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

explicit

Definition at line 69 of file SiStripApvGainRescaler.cc.

References edm::ParameterSet::getUntrackedParameter().

     : m_printdebug{iConfig.getUntrackedParameter<uint32_t>("printDebug", 1)},
       m_Record(iConfig.getParameter<std::string>("Record")),
       g1g2Token_(esConsumes()),
       g3Token_(esConsumes()) {}

◆ ~SiStripApvGainRescaler()

SiStripApvGainRescaler::~SiStripApvGainRescaler ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 81 of file SiStripApvGainRescaler.cc.

References cond::service::PoolDBOutputService::currentTime(), ztail::d, g1g2Token_, g3Token_, edm::EventSetup::getData(), getNewObject(), edm::Service< T >::isAvailable(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, m_Record, and cond::service::PoolDBOutputService::writeOneIOV().

                                                                                         {
   using namespace edm;
 
   const auto& g1g2 = iSetup.getData(g1g2Token_);
   const auto& g3 = iSetup.getData(g3Token_);
 
   std::map<std::pair<uint32_t, int>, float> theMap;
 
   std::vector<uint32_t> detid;
   g1g2.getDetIds(detid);
   for (const auto& d : detid) {
     SiStripApvGain::Range rangeG1_old = g1g2.getRange(d, 0);
     SiStripApvGain::Range rangeG2_old = g1g2.getRange(d, 1);
     SiStripApvGain::Range rangeG1_new = g3.getRange(d);
 
     int nAPV = 0;
     for (int it = 0; it < rangeG1_old.second - rangeG1_old.first; it++) {
       nAPV++;
 
       std::pair<uint32_t, int> index = std::make_pair(d, nAPV);
 
       float G1_old = g1g2.getApvGain(it, rangeG1_old);
       float G2_old = g1g2.getApvGain(it, rangeG2_old);
       float G1G2_old = G1_old * G2_old;
       float G1_new = g3.getApvGain(it, rangeG1_new);
 
       // this is based on G1_old*G2_old = G1_new * G2_new ==> G2_new = (G1_old*G2_old)/G1_new
 
       float NewGain = G1G2_old / G1_new;
 
       // DO NOT RESCALE APVs set to the default value
       if (G2_old != 1.) {
         theMap[index] = NewGain;
       } else {
         theMap[index] = 1.;
       }
 
     }  // loop over APVs
   }    // loop over DetIds
 
   std::unique_ptr<SiStripApvGain> theAPVGains = this->getNewObject(theMap);
 
   // write out the APVGains record
   edm::Service<cond::service::PoolDBOutputService> poolDbService;
 
   if (poolDbService.isAvailable())
     poolDbService->writeOneIOV(*theAPVGains, poolDbService->currentTime(), m_Record);
   else
     throw std::runtime_error("PoolDBService required.");
 }

◆ fillDescriptions()

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

static

Definition at line 170 of file SiStripApvGainRescaler.cc.

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

                                                                                         {
   edm::ParameterSetDescription desc;
 
   desc.setComment(
       " Utility class to rescale the values of SiStrip G2 by the ratio of G1_old/G1_new: this is useful in the case in "
       "which a Gain2 payload needs to recycled after a G1 update to keep the G1*G2 product constant."
       "PoolDBOutputService must be set up for 'SiStripApvGainRcd'.");
 
   desc.add<std::string>("Record", "SiStripApvGainRcd");
   desc.addUntracked<unsigned int>("printDebug", 1);
   descriptions.add("rescaleGain2byGain1", desc);
 }

◆ getNewObject()

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

private

Definition at line 133 of file SiStripApvGainRescaler.cc.

References m_printdebug, getGTfromDQMFile::obj, and FastTimerService_cff::range.

Referenced by analyze().

                                                          {
   std::unique_ptr<SiStripApvGain> obj = std::make_unique<SiStripApvGain>();
 
   std::vector<float> theSiStripVector;
   uint32_t PreviousDetId = 0;
   unsigned int countDetIds(0);  // count DetIds to print
   for (const auto& element : theMap) {
     uint32_t DetId = element.first.first;
     if (DetId != PreviousDetId) {
       if (!theSiStripVector.empty()) {
         SiStripApvGain::Range range(theSiStripVector.begin(), theSiStripVector.end());
         if (!obj->put(PreviousDetId, range))
           edm::LogError("SiStripApvGainRescaler") << "Bug to put detId = " << PreviousDetId << "\n";
       }
       theSiStripVector.clear();
       PreviousDetId = DetId;
       countDetIds++;
     }
     theSiStripVector.push_back(element.second);
 
     if (countDetIds <= m_printdebug) {
       edm::LogInfo("SiStripApvGainRescaler")
           << __FUNCTION__ << " DetId: " << DetId << " APV:   " << element.first.second << " Gain:  " << element.second;
     }
   }
 
   if (!theSiStripVector.empty()) {
     SiStripApvGain::Range range(theSiStripVector.begin(), theSiStripVector.end());
     if (!obj->put(PreviousDetId, range))
       edm::LogError("SiStripApvGainRescaler") << "Bug to put detId = " << PreviousDetId << "\n";
   }
 
   return obj;
 }

Member Data Documentation

◆ g1g2Token_

const edm::ESGetToken<SiStripGain, SiStripGainRcd> SiStripApvGainRescaler::g1g2Token_

private

Definition at line 61 of file SiStripApvGainRescaler.cc.

Referenced by analyze().

◆ g3Token_

const edm::ESGetToken<SiStripApvGain, SiStripApvGain3Rcd> SiStripApvGainRescaler::g3Token_

private

Definition at line 63 of file SiStripApvGainRescaler.cc.

Referenced by analyze().

◆ m_printdebug

const uint32_t SiStripApvGainRescaler::m_printdebug

private

Definition at line 57 of file SiStripApvGainRescaler.cc.

Referenced by getNewObject().

◆ m_Record

const std::string SiStripApvGainRescaler::m_Record

private

Definition at line 58 of file SiStripApvGainRescaler.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SiStripApvGainRescaler()

◆ ~SiStripApvGainRescaler()

Member Function Documentation

◆ analyze()

◆ fillDescriptions()

◆ getNewObject()

Member Data Documentation

◆ g1g2Token_

◆ g3Token_

◆ m_printdebug

◆ m_Record