#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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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)

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_

unsigned long long	m_cacheID_

double	meanGain_

double	minimumPosValue_

bool	printdebug_

double	sigmaGain_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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)

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

Definition at line 27 of file SiStripGainRandomCalculator.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 23 of file SiStripGainRandomCalculator.cc.

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

     : ConditionDBWriter<SiStripApvGain>(iConfig), m_cacheID_(0) {
   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);
 }

SiStripGainRandomCalculator::~SiStripGainRandomCalculator ( )

override

Definition at line 40 of file SiStripGainRandomCalculator.cc.

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

Member Function Documentation

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

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 44 of file SiStripGainRandomCalculator.cc.

References detid_apvs_, TrackerGeometry::detUnits(), beamvalidation::exit(), edm::EventSetup::get(), m_cacheID_, StripTopology::nstrips(), AlCaHLTBitMon_ParallelJobs::p, and printdebug_.

                                                                                                 {
   unsigned long long cacheID = iSetup.get<TrackerDigiGeometryRecord>().cacheIdentifier();
 
   if (m_cacheID_ != cacheID) {
     m_cacheID_ = cacheID;
 
     edm::ESHandle<TrackerGeometry> pDD;
 
     iSetup.get<TrackerDigiGeometryRecord>().get(pDD);
     edm::LogInfo("SiStripGainRandomCalculator::algoAnalyze - got new geometry  ") << std::endl;
 
     detid_apvs_.clear();
 
     edm::LogInfo("SiStripGainCalculator") << " There are " << pDD->detUnits().size() << " detectors" << std::endl;
 
     for (const auto& it : pDD->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;
       }
     }
   }
 }

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

overrideprivatevirtual

Implements ConditionDBWriter< SiStripApvGain >.

Definition at line 77 of file SiStripGainRandomCalculator.cc.

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

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