#include <SiStripNoiseNormalizedWithApvGainBuilder.h>

Inheritance diagram for SiStripNoiseNormalizedWithApvGainBuilder:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

	SiStripNoiseNormalizedWithApvGainBuilder (const edm::ParameterSet &iConfig)

	~SiStripNoiseNormalizedWithApvGainBuilder () 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 &&)=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)

Private Member Functions
void	fillParameters (std::map< int, std::vector< double > > &mapToFill, const std::string &parameterName) const
	Fills the parameters read from cfg and matching the name in the given map. More...

void	fillSubDetParameter (std::map< int, std::vector< double > > &mapToFill, const std::vector< double > &v, const int subDet, const unsigned short layers) const

void	printLog (const uint32_t detId, const unsigned short strip, const double &noise) const

Private Attributes
double	electronsPerADC_

edm::FileInPath	fp_

edm::ESGetToken< SiStripApvGain, SiStripApvGainRcd >	inputApvGainToken_

double	minimumPosValue_

bool	printdebug_

uint32_t	printDebug_

edm::ParameterSet	pset_

bool	stripLengthMode_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tGeomToken_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

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

Produces a noise tag using the same settings as the service used in the DummyDBWriter, but it receives a SiStripApvGain tag from the EventSetup and uses the gain values (per apv) to rescale the noise (per strip).

Definition at line 32 of file SiStripNoiseNormalizedWithApvGainBuilder.h.

Constructor & Destructor Documentation

◆ SiStripNoiseNormalizedWithApvGainBuilder()

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

explicit

Definition at line 12 of file SiStripNoiseNormalizedWithApvGainBuilder.cc.

     : printdebug_(iConfig.getUntrackedParameter<uint32_t>("printDebug", 1)),
       pset_(iConfig),
       electronsPerADC_(0.),
       minimumPosValue_(0.),
       stripLengthMode_(true),
       printDebug_(0) {
   tTopoToken_ = esConsumes<TrackerTopology, TrackerTopologyRcd>();
   tGeomToken_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>();
   inputApvGainToken_ = esConsumes<SiStripApvGain, SiStripApvGainRcd>();
 }

References inputApvGainToken_, tGeomToken_, and tTopoToken_.

◆ ~SiStripNoiseNormalizedWithApvGainBuilder()

SiStripNoiseNormalizedWithApvGainBuilder::~SiStripNoiseNormalizedWithApvGainBuilder ( )

inlineoverride

Definition at line 36 of file SiStripNoiseNormalizedWithApvGainBuilder.h.

36 {};

Member Function Documentation

◆ analyze()

void SiStripNoiseNormalizedWithApvGainBuilder::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 24 of file SiStripNoiseNormalizedWithApvGainBuilder.cc.

                                                                                                        {
   const auto& tTopo = iSetup.getData(tTopoToken_);
   const auto& tGeom = iSetup.getData(tGeomToken_);
   // Read the gain from the given tag
   const auto& inputApvGain = iSetup.getData(inputApvGainToken_);
   std::vector<uint32_t> inputDetIds;
   inputApvGain.getDetIds(inputDetIds);
  
   // Prepare the new object
   SiStripNoises* obj = new SiStripNoises();
  
   stripLengthMode_ = pset_.getParameter<bool>("StripLengthMode");
  
   //parameters for random noise generation. not used if Strip length mode is chosen
   SiStripFakeAPVParameters meanNoise{pset_, "MeanNoise"};
   SiStripFakeAPVParameters sigmaNoise{pset_, "SigmaNoise"};
   minimumPosValue_ = pset_.getParameter<double>("MinPositiveNoise");
  
   //parameters for strip length proportional noise generation. not used if random mode is chosen
   SiStripFakeAPVParameters noiseStripLengthLinearSlope{pset_, "NoiseStripLengthSlope"};
   SiStripFakeAPVParameters noiseStripLengthLinearQuote{pset_, "NoiseStripLengthQuote"};
   electronsPerADC_ = pset_.getParameter<double>("electronPerAdc");
  
   printDebug_ = pset_.getUntrackedParameter<uint32_t>("printDebug", 5);
  
   unsigned int count = 0;
   for (const auto det : tGeom.detUnits()) {
     const StripGeomDetUnit* stripDet = dynamic_cast<const StripGeomDetUnit*>(det);
     if (stripDet != nullptr) {
       const DetId detId = stripDet->geographicalId();
       // Find if this DetId is in the input tag and if so how many are the Apvs for which it contains information
       SiStripApvGain::Range inputRange(inputApvGain.getRange(detId));
  
       //Generate Noises for det detid
       SiStripNoises::InputVector theSiStripVector;
       float noise = 0.;
       SiStripFakeAPVParameters::index sl = SiStripFakeAPVParameters::getIndex(&tTopo, detId);
       unsigned short nApvs = stripDet->specificTopology().nstrips() / 128;
  
       if (stripLengthMode_) {
         // Use strip length
         double linearSlope = noiseStripLengthLinearSlope.get(sl);
         double linearQuote = noiseStripLengthLinearQuote.get(sl);
         double stripLength = stripDet->specificTopology().stripLength();
         for (unsigned short j = 0; j < nApvs; ++j) {
           double gain = inputApvGain.getApvGain(j, inputRange);
  
           for (unsigned short stripId = 0; stripId < 128; ++stripId) {
             noise = ((linearSlope * stripLength + linearQuote) / electronsPerADC_) * gain;
             if (count < printDebug_)
               printLog(detId, stripId + 128 * j, noise);
             obj->setData(noise, theSiStripVector);
           }
         }
       } else {
         // Use random generator
         double meanN = meanNoise.get(sl);
         double sigmaN = sigmaNoise.get(sl);
         for (unsigned short j = 0; j < nApvs; ++j) {
           double gain = inputApvGain.getApvGain(j, inputRange);
  
           for (unsigned short stripId = 0; stripId < 128; ++stripId) {
             noise = (CLHEP::RandGauss::shoot(meanN, sigmaN)) * gain;
             if (noise <= minimumPosValue_)
               noise = minimumPosValue_;
             if (count < printDebug_)
               printLog(detId, stripId + 128 * j, noise);
             obj->setData(noise, theSiStripVector);
           }
         }
       }
       ++count;
  
       if (!obj->put(detId, theSiStripVector)) {
         edm::LogError("SiStripNoisesFakeESSource::produce ") << " detid already exists" << std::endl;
       }
     }
   }
  
   //End now write data in DB
   edm::Service<cond::service::PoolDBOutputService> mydbservice;
  
   if (mydbservice.isAvailable()) {
     if (mydbservice->isNewTagRequest("SiStripNoisesRcd")) {
       mydbservice->createNewIOV<SiStripNoises>(
           obj, mydbservice->beginOfTime(), mydbservice->endOfTime(), "SiStripNoisesRcd");
     } else {
       mydbservice->appendSinceTime<SiStripNoises>(obj, mydbservice->currentTime(), "SiStripNoisesRcd");
     }
   } else {
     edm::LogError("SiStripNoiseNormalizedWithApvGainBuilder") << "Service is unavailable" << std::endl;
   }
 }

References cond::service::PoolDBOutputService::appendSinceTime(), cond::service::PoolDBOutputService::beginOfTime(), KineDebug3::count(), cond::service::PoolDBOutputService::createNewIOV(), cond::service::PoolDBOutputService::currentTime(), electronsPerADC_, cond::service::PoolDBOutputService::endOfTime(), PedestalClient_cfi::gain, GeomDet::geographicalId(), edm::EventSetup::getData(), SiStripFakeAPVParameters::getIndex(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), inputApvGainToken_, pileupCalc::inputRange, edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), dqmiolumiharvest::j, minimumPosValue_, hgchebackDigitizer_cfi::noise, StripTopology::nstrips(), getGTfromDQMFile::obj, printDebug_, printLog(), pset_, StripGeomDetUnit::specificTopology(), StripTopology::stripLength(), stripLengthMode_, tGeomToken_, and tTopoToken_.

◆ fillParameters()

void SiStripNoiseNormalizedWithApvGainBuilder::fillParameters	(	std::map< int, std::vector< double > > &	mapToFill,
		const std::string &	parameterName
	)		const

private

Fills the parameters read from cfg and matching the name in the given map.

◆ fillSubDetParameter()

void SiStripNoiseNormalizedWithApvGainBuilder::fillSubDetParameter	(	std::map< int, std::vector< double > > &	mapToFill,
		const std::vector< double > &	v,
		const int	subDet,
		const unsigned short	layers
	)		const

private

Fills the map with the paramters for the given subdetector.
Each vector "v" holds the parameters for the layers/rings, if the vector has only one parameter all the layers/rings get that parameter.
The only other possibility is that the number of parameters equals the number of layers, otherwise an exception of type "Configuration" will be thrown.

◆ printLog()

void SiStripNoiseNormalizedWithApvGainBuilder::printLog	(	const uint32_t	detId,
		const unsigned short	strip,
		const double &	noise
	)		const

inlineprivate

Definition at line 55 of file SiStripNoiseNormalizedWithApvGainBuilder.h.

                                                                                                     {
     edm::LogInfo("SiStripNoisesDummyCalculator")
         << "detid: " << detId << " strip: " << strip << " noise: " << noise << " \t" << std::endl;
   }