SiStripApvGainRescaler Class Reference

#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
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 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)
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
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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(), SiStripApvGain::getRange(), SiStripGain::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.

{
}

void SiStripApvGainRescaler::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 156 of file SiStripApvGainRescaler.cc.

{
}

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 hgcalPlots::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().