#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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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 &&)=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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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	beginJob () override

void	endJob () override

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

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

edm::ESGetToken< SiStripApvGain, SiStripApvGain3Rcd >	g3Token_

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

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 70 of file SiStripApvGainRescaler.cc.

     : m_Record(iConfig.getParameter<std::string>("Record")), g1g2Token_(esConsumes()), g3Token_(esConsumes()) {
   //now do what ever initialization is needed
 }

◆ ~SiStripApvGainRescaler()

SiStripApvGainRescaler::~SiStripApvGainRescaler ( )

override

Definition at line 75 of file SiStripApvGainRescaler.cc.

                                                 {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 85 of file SiStripApvGainRescaler.cc.

                                                                                         {
   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->writeOne(theAPVGains.get(), poolDbService->currentTime(), m_Record);
   else
     throw std::runtime_error("PoolDBService required.");
 }

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

◆ beginJob()

void SiStripApvGainRescaler::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 137 of file SiStripApvGainRescaler.cc.

137 {}

◆ endJob()

void SiStripApvGainRescaler::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 140 of file SiStripApvGainRescaler.cc.

140 {}

◆ fillDescriptions()

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

static

Definition at line 176 of file SiStripApvGainRescaler.cc.

                                                                                         {
   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");
   descriptions.add("rescaleGain2byGain1", desc);
 }

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

◆ getNewObject()

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

private

Definition at line 143 of file SiStripApvGainRescaler.cc.

                                                          {
   std::unique_ptr<SiStripApvGain> obj = std::make_unique<SiStripApvGain>();
  
   std::vector<float> theSiStripVector;
   uint32_t PreviousDetId = 0;
   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))
           printf("Bug to put detId = %i\n", PreviousDetId);
       }
       theSiStripVector.clear();
       PreviousDetId = DetId;
     }
     theSiStripVector.push_back(element.second);
  
     edm::LogInfo("SiStripApvGainRescaler")
         << " DetId: " << DetId << " APV:   " << element.first.second << " Gain:  " << element.second << std::endl;
   }
  
   if (!theSiStripVector.empty()) {
     SiStripApvGain::Range range(theSiStripVector.begin(), theSiStripVector.end());
     if (!obj->put(PreviousDetId, range))
       printf("Bug to put detId = %i\n", PreviousDetId);
   }
  
   return obj;
 }

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

Referenced by analyze().

Member Data Documentation

◆ g1g2Token_

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

private

Definition at line 62 of file SiStripApvGainRescaler.cc.

Referenced by analyze().

◆ g3Token_

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

private

Definition at line 64 of file SiStripApvGainRescaler.cc.

Referenced by analyze().

◆ m_Record

const std::string SiStripApvGainRescaler::m_Record

private

Definition at line 59 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()

◆ beginJob()

◆ endJob()

◆ fillDescriptions()

◆ getNewObject()

Member Data Documentation

◆ g1g2Token_

◆ g3Token_

◆ m_Record