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

Inheritance diagram for SiStripApvGainRescaler:

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

	~SiStripApvGainRescaler ()

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

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

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
virtual void	analyze (const edm::Event &, const edm::EventSetup &) override

virtual void	beginJob () override

virtual 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
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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)

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

Constructor & Destructor Documentation

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

explicit

Definition at line 68 of file SiStripApvGainRescaler.cc.

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

SiStripApvGainRescaler::~SiStripApvGainRescaler ( )

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

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

overrideprivatevirtual

Definition at line 88 of file SiStripApvGainRescaler.cc.

References cond::service::PoolDBOutputService::currentTime(), edmIntegrityCheck::d, edm::EventSetup::get(), SiStripGain::getApvGain(), SiStripApvGain::getApvGain(), SiStripGain::getDetIds(), getNewObject(), SiStripGain::getRange(), SiStripApvGain::getRange(), edm::Service< T >::isAvailable(), m_Record, and cond::service::PoolDBOutputService::writeOne().

 {
    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.");
  
 }

void SiStripApvGainRescaler::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 150 of file SiStripApvGainRescaler.cc.

151 {

152 }

void SiStripApvGainRescaler::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 156 of file SiStripApvGainRescaler.cc.

157 {

158 }

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

static

Definition at line 197 of file SiStripApvGainRescaler.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, edm::ParameterSetDescription::setComment(), 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");
   descriptions.add("rescaleGain2byGain1", desc);
 }

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

private

Definition at line 162 of file SiStripApvGainRescaler.cc.

References GetRecoTauVFromDQM_MC_cff::obj.

Referenced by analyze().

 {
   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;
 }

Member Data Documentation

const std::string SiStripApvGainRescaler::m_Record

private

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

Member Function Documentation

Member Data Documentation