#include <SiStripNoiseNormalizedWithApvGainBuilder.h>

Inheritance diagram for SiStripNoiseNormalizedWithApvGainBuilder:

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

	SiStripNoiseNormalizedWithApvGainBuilder (const edm::ParameterSet &iConfig)

	~SiStripNoiseNormalizedWithApvGainBuilder ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

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

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_

double	minimumPosValue_

bool	printdebug_

uint32_t	printDebug_

edm::ParameterSet	pset_

bool	stripLengthMode_

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

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

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 28 of file SiStripNoiseNormalizedWithApvGainBuilder.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 13 of file SiStripNoiseNormalizedWithApvGainBuilder.cc.

                                                                                                                   :
   fp_(iConfig.getUntrackedParameter<edm::FileInPath>("file",edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat"))),
   printdebug_(iConfig.getUntrackedParameter<uint32_t>("printDebug",1)),
   pset_(iConfig),
   electronsPerADC_(0.),
   minimumPosValue_(0.),
   stripLengthMode_(true),
   printDebug_(0)
 {}

SiStripNoiseNormalizedWithApvGainBuilder::~SiStripNoiseNormalizedWithApvGainBuilder ( )

inline

Definition at line 34 of file SiStripNoiseNormalizedWithApvGainBuilder.h.

References analyze(), fillParameters(), fillSubDetParameter(), LayerTriplets::layers(), genParticles_cff::map, AlCaHLTBitMon_QueryRunRegistry::string, and findQualityFiles::v.

34 {};

Member Function Documentation

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

virtual

Definition at line 23 of file SiStripNoiseNormalizedWithApvGainBuilder.cc.

References cond::service::PoolDBOutputService::appendSinceTime(), cond::service::PoolDBOutputService::beginOfTime(), KineDebug3::count(), cond::service::PoolDBOutputService::createNewIOV(), cond::service::PoolDBOutputService::currentTime(), electronsPerADC_, cond::service::PoolDBOutputService::endOfTime(), fp_, edm::FileInPath::fullPath(), muonCSCDigis_cfi::gain, edm::EventSetup::get(), SiStripApvGain::getApvGain(), SiStripApvGain::getDetIds(), SiStripFakeAPVParameters::getIndex(), edm::ParameterSet::getParameter(), SiStripApvGain::getRange(), edm::ParameterSet::getUntrackedParameter(), estimatePileup::inputRange, edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), minimumPosValue_, MuonAssociatorByHits_cfi::obj, printDebug_, printLog(), edm::ESHandle< T >::product(), pset_, SiStripNoises::put(), matplotRender::reader, SiStripNoises::setData(), and stripLengthMode_.

Referenced by ~SiStripNoiseNormalizedWithApvGainBuilder().

 {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<IdealGeometryRecord>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
  
   // Read the gain from the given tag
   edm::ESHandle<SiStripApvGain> inputApvGain;
   iSetup.get<SiStripApvGainRcd>().get( inputApvGain );
   std::vector<uint32_t> inputDetIds;
   inputApvGain->getDetIds(inputDetIds);
 
   // Prepare the new object
   SiStripNoises* obj = new SiStripNoises();
 
   SiStripDetInfoFileReader reader(fp_.fullPath());
 
 
 
   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;
   const std::map<uint32_t, SiStripDetInfoFileReader::DetInfo > DetInfos = reader.getAllData();
   for(std::map<uint32_t, SiStripDetInfoFileReader::DetInfo >::const_iterator it = DetInfos.begin(); it != DetInfos.end(); it++) {
 
     // 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(it->first));
 
     //Generate Noises for det detid
     SiStripNoises::InputVector theSiStripVector;
     float noise = 0.;
     uint32_t detId = it->first;
     SiStripFakeAPVParameters::index sl =  SiStripFakeAPVParameters::getIndex(tTopo, detId);
     unsigned short nApvs = it->second.nApvs;
 
     if(stripLengthMode_) {
       // Use strip length
       double linearSlope = noiseStripLengthLinearSlope.get(sl);
       double linearQuote = noiseStripLengthLinearQuote.get(sl);
       double stripLength = it->second.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(it->first,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;
   }
 }

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.

Referenced by ~SiStripNoiseNormalizedWithApvGainBuilder().

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.

Referenced by ~SiStripNoiseNormalizedWithApvGainBuilder().

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

inlineprivate

Definition at line 50 of file SiStripNoiseNormalizedWithApvGainBuilder.h.

Referenced by analyze().

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