#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

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

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::vector< ModuleDescription const * > &modules, 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
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>
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)

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

Constructor & Destructor Documentation

◆ SiStripApvGainRescaler()

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

explicit

Definition at line 66 of file SiStripApvGainRescaler.cc.

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

◆ ~SiStripApvGainRescaler()

SiStripApvGainRescaler::~SiStripApvGainRescaler ( )

override

Definition at line 71 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 81 of file SiStripApvGainRescaler.cc.

                                                                                         {
   using namespace edm;
  
   // take G2_old and G1_old from the regular gain handle
   edm::ESHandle<SiStripGain> g1g2Handle_;
   iSetup.get<SiStripGainRcd>().get(g1g2Handle_);
  
   // take the additional G1_new from the Gain3Rcd (dirty trick)
   edm::ESHandle<SiStripApvGain> g3Handle_;
   iSetup.get<SiStripApvGain3Rcd>().get(g3Handle_);
  
   std::map<std::pair<uint32_t, int>, float> theMap;
  
   std::vector<uint32_t> detid;
   g1g2Handle_->getDetIds(detid);
   for (const auto& d : detid) {
     SiStripApvGain::Range rangeG1_old = g1g2Handle_->getRange(d, 0);
     SiStripApvGain::Range rangeG2_old = g1g2Handle_->getRange(d, 1);
     SiStripApvGain::Range rangeG1_new = g3Handle_->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 = g1g2Handle_->getApvGain(it, rangeG1_old);
       float G2_old = g1g2Handle_->getApvGain(it, rangeG2_old);
       float G1G2_old = G1_old * G2_old;
       float G1_new = g3Handle_->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, edm::EventSetup::get(), get, SiStripGain::getApvGain(), SiStripApvGain::getApvGain(), SiStripGain::getDetIds(), getNewObject(), SiStripApvGain::getRange(), SiStripGain::getRange(), 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 138 of file SiStripApvGainRescaler.cc.

138 {}

◆ endJob()

void SiStripApvGainRescaler::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 141 of file SiStripApvGainRescaler.cc.

141 {}

◆ fillDescriptions()

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

static

Definition at line 177 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(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::setComment(), 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 144 of file SiStripApvGainRescaler.cc.

                                                          {
   std::unique_ptr<SiStripApvGain> obj = std::unique_ptr<SiStripApvGain>(new 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

◆ m_Record

const std::string SiStripApvGainRescaler::m_Record

private

Definition at line 60 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

◆ m_Record