#include <SiStripGainRandomCalculator.h>

Inheritance diagram for SiStripGainRandomCalculator:

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

	~SiStripGainRandomCalculator () override

Public Member Functions inherited from ConditionDBWriter< SiStripApvGain >
	ConditionDBWriter (const edm::ParameterSet &iConfig)

	~ConditionDBWriter () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources >
	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

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	algoAnalyze (const edm::Event &, const edm::EventSetup &) override

std::unique_ptr< SiStripApvGain >	getNewObject () override

Private Attributes
std::vector< std::pair< uint32_t, unsigned short > >	detid_apvs_

double	meanGain_

double	minimumPosValue_

bool	printdebug_

double	sigmaGain_

edm::ESWatcher< TrackerDigiGeometryRecord >	tkDigiGeomRcdWatcher_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from ConditionDBWriter< SiStripApvGain >
static void	fillPSetDescription (edm::ParameterSetDescription &desc)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from ConditionDBWriter< SiStripApvGain >
void	setDoStore (const bool doStore)
	When set to false the payload will not be written to the db. More...

void	storeOnDbNow ()

cond::Time_t	timeOfLastIOV ()

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

Definition at line 31 of file SiStripGainRandomCalculator.h.

Constructor & Destructor Documentation

◆ SiStripGainRandomCalculator()

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

explicit

Definition at line 20 of file SiStripGainRandomCalculator.cc.

References detid_apvs_, edm::EDConsumerBase::esConsumes(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), meanGain_, minimumPosValue_, printdebug_, sigmaGain_, and tkGeomToken_.

     : ConditionDBWriter<SiStripApvGain>(iConfig) {
   edm::LogInfo("SiStripGainRandomCalculator::SiStripGainRandomCalculator");
 
   //   std::string Mode=iConfig.getParameter<std::string>("Mode");
   //   if (Mode==std::string("Gaussian")) GaussianMode_=true;
   //   else if (IOVMode==std::string("Constant")) ConstantMode_=true;
   //   else  edm::LogError("SiStripGainRandomCalculator::SiStripGainRandomCalculator(): ERROR - unknown generation mode...will not store anything on the DB") << std::endl;
 
   detid_apvs_.clear();
 
   meanGain_ = iConfig.getParameter<double>("MeanGain");
   sigmaGain_ = iConfig.getParameter<double>("SigmaGain");
   minimumPosValue_ = iConfig.getParameter<double>("MinPositiveGain");
   printdebug_ = iConfig.getUntrackedParameter<bool>("printDebug", false);
 
   tkGeomToken_ = esConsumes();
 }

◆ ~SiStripGainRandomCalculator()

SiStripGainRandomCalculator::~SiStripGainRandomCalculator ( )

override

Definition at line 39 of file SiStripGainRandomCalculator.cc.

                                                           {
   edm::LogInfo("SiStripGainRandomCalculator::~SiStripGainRandomCalculator");
 }

Member Function Documentation

◆ algoAnalyze()

void SiStripGainRandomCalculator::algoAnalyze	(	const edm::Event &	event,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 43 of file SiStripGainRandomCalculator.cc.

References edm::ESWatcher< T >::check(), createTree::dd, detid_apvs_, beamvalidation::exit(), edm::EventSetup::getData(), AlCaHLTBitMon_ParallelJobs::p, printdebug_, tkDigiGeomRcdWatcher_, and tkGeomToken_.

                                                                                                 {
   if (tkDigiGeomRcdWatcher_.check(iSetup)) {
     const auto& dd = iSetup.getData(tkGeomToken_);
     edm::LogInfo("SiStripGainRandomCalculator::algoAnalyze - got new geometry  ") << std::endl;
 
     detid_apvs_.clear();
 
     edm::LogInfo("SiStripGainCalculator") << " There are " << dd.detUnits().size() << " detectors" << std::endl;
 
     for (const auto& it : dd.detUnits()) {
       if (dynamic_cast<const StripGeomDetUnit*>(it) != nullptr) {
         uint32_t detid = (it->geographicalId()).rawId();
         const StripTopology& p = dynamic_cast<const StripGeomDetUnit*>(it)->specificTopology();
         unsigned short NAPVs = p.nstrips() / 128;
         if (NAPVs < 1 || NAPVs > 6) {
           edm::LogError("SiStripGainCalculator")
               << " Problem with Number of strips in detector.. " << p.nstrips() << " Exiting program" << endl;
           exit(1);
         }
         detid_apvs_.push_back(pair<uint32_t, unsigned short>(detid, NAPVs));
         if (printdebug_)
           edm::LogInfo("SiStripGainCalculator") << "detid " << detid << " apvs " << NAPVs;
       }
     }
   }
 }

◆ getNewObject()

std::unique_ptr< SiStripApvGain > SiStripGainRandomCalculator::getNewObject ( )

overrideprivatevirtual

Implements ConditionDBWriter< SiStripApvGain >.

Definition at line 70 of file SiStripGainRandomCalculator.cc.

References gather_cfg::cout, detid_apvs_, PedestalClient_cfi::gain, dqmiolumiharvest::j, meanGain_, minimumPosValue_, getGTfromDQMFile::obj, printdebug_, FastTimerService_cff::range, sigmaGain_, and trackerHitRTTI::vector.

                                                                         {
   std::cout << "SiStripGainRandomCalculator::getNewObject called" << std::endl;
 
   auto obj = std::make_unique<SiStripApvGain>();
 
   for (std::vector<pair<uint32_t, unsigned short> >::const_iterator it = detid_apvs_.begin(); it != detid_apvs_.end();
        it++) {
     //Generate Gain for det detid
     std::vector<float> theSiStripVector;
     for (unsigned short j = 0; j < it->second; j++) {
       float gain;
 
       //      if(sigmaGain_/meanGain_ < 0.00001) gain = meanGain_;
       //      else{
       gain = CLHEP::RandGauss::shoot(meanGain_, sigmaGain_);
       if (gain <= minimumPosValue_)
         gain = minimumPosValue_;
       //      }
 
       if (printdebug_)
         edm::LogInfo("SiStripGainCalculator") << "detid " << it->first << " \t"
                                               << " apv " << j << " \t" << gain << " \t" << std::endl;
       theSiStripVector.push_back(gain);
     }
 
     SiStripApvGain::Range range(theSiStripVector.begin(), theSiStripVector.end());
     if (!obj->put(it->first, range))
       edm::LogError("SiStripGainCalculator") << "[SiStripGainCalculator::beginJob] detid already exists" << std::endl;
   }
 
   return obj;
 }