Inheritance diagram for SiStripLorentzAngleFakeESSource:

Public Types
typedef std::unique_ptr< SiStripLorentzAngle >	ReturnType

Public Member Functions
ReturnType	produce (const SiStripLorentzAngleRcd &)

void	setIntervalFor (const edm::eventsetup::EventSetupRecordKey &, const edm::IOVSyncValue &iov, edm::ValidityInterval &iValidity) override

	SiStripLorentzAngleFakeESSource (const edm::ParameterSet &)

	~SiStripLorentzAngleFakeESSource () override

Public Member Functions inherited from edm::ESProducer
	ESProducer ()

	ESProducer (const ESProducer &)=delete

	ESProducer (ESProducer &&)=delete

ESResolverIndex const *	getTokenIndices (unsigned int iIndex) const

ESRecordIndex const *	getTokenRecordIndices (unsigned int iIndex) const

bool	hasMayConsumes () const noexcept

size_t	numberOfTokenIndices (unsigned int iIndex) const

ESProducer &	operator= (const ESProducer &)=delete

ESProducer &	operator= (ESProducer &&)=delete

SerialTaskQueueChain &	queue ()

template<typename Record >
std::optional< std::vector< ESResolverIndex > >	updateFromMayConsumes (unsigned int iIndex, const Record &iRecord) const

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &) final

	~ESProducer () noexcept(false) override

Public Member Functions inherited from edm::ESProductResolverFactoryProducer
	ESProductResolverFactoryProducer ()

	ESProductResolverFactoryProducer (const ESProductResolverFactoryProducer &)=delete

const ESProductResolverFactoryProducer &	operator= (const ESProductResolverFactoryProducer &)=delete

	~ESProductResolverFactoryProducer () noexcept(false) override

Public Member Functions inherited from edm::eventsetup::ESProductResolverProvider
void	createKeyedResolvers (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)

const ComponentDescription &	description () const

	ESProductResolverProvider ()

	ESProductResolverProvider (const ESProductResolverProvider &)=delete

void	fillRecordsNotAllowingConcurrentIOVs (std::set< EventSetupRecordKey > &recordsNotAllowingConcurrentIOVs) const

virtual void	initConcurrentIOVs (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)

bool	isUsingRecord (const EventSetupRecordKey &key) const

KeyedResolvers &	keyedResolvers (const EventSetupRecordKey &iRecordKey, unsigned int iovIndex=0)

const ESProductResolverProvider &	operator= (const ESProductResolverProvider &)=delete

void	setAppendToDataLabel (const edm::ParameterSet &)

void	setDescription (const ComponentDescription &iDescription)

std::set< EventSetupRecordKey >	usingRecords () const

virtual	~ESProductResolverProvider () noexcept(false)

Public Member Functions inherited from edm::EventSetupRecordIntervalFinder
bool	concurrentFinder () const

const eventsetup::ComponentDescription &	descriptionForFinder () const

	EventSetupRecordIntervalFinder ()

	EventSetupRecordIntervalFinder (const EventSetupRecordIntervalFinder &)=delete

std::set< eventsetup::EventSetupRecordKey >	findingForRecords () const

const ValidityInterval &	findIntervalFor (const eventsetup::EventSetupRecordKey &, const IOVSyncValue &)

bool	nonconcurrentAndIOVNeedsUpdate (const eventsetup::EventSetupRecordKey &key, const IOVSyncValue &syncValue) const

const EventSetupRecordIntervalFinder &	operator= (const EventSetupRecordIntervalFinder &)=delete

void	resetInterval (const eventsetup::EventSetupRecordKey &)

void	setDescriptionForFinder (const eventsetup::ComponentDescription &iDescription)

virtual	~EventSetupRecordIntervalFinder () noexcept(false)

Private Attributes
edm::ESGetToken< GeometricDet, IdealGeometryRecord >	m_geomDetToken

std::vector< double >	m_StdDevs_TIB

std::vector< double >	m_StdDevs_TOB

std::vector< double >	m_TIB_EstimatedValuesMax

std::vector< double >	m_TIB_EstimatedValuesMin

std::vector< double >	m_TIB_PerCent_Errs

double	m_TIBmeanPerCentError

double	m_TIBmeanStdDev

double	m_TIBmeanValueMax

double	m_TIBmeanValueMin

std::vector< double >	m_TOB_EstimatedValuesMax

std::vector< double >	m_TOB_EstimatedValuesMin

std::vector< double >	m_TOB_PerCent_Errs

double	m_TOBmeanPerCentError

double	m_TOBmeanStdDev

double	m_TOBmeanValueMax

double	m_TOBmeanValueMin

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	m_tTopoToken

std::vector< bool >	m_uniformTIB

std::vector< bool >	m_uniformTOB

Additional Inherited Members
Static Public Member Functions inherited from edm::eventsetup::ESProductResolverProvider
static void	prevalidate (ConfigurationDescriptions &)

Protected Types inherited from edm::ESProductResolverFactoryProducer
using	EventSetupRecordKey = eventsetup::EventSetupRecordKey

Protected Types inherited from edm::eventsetup::ESProductResolverProvider
using	KeyedResolversVector = std::vector< std::pair< DataKey, std::shared_ptr< ESProductResolver > >>

Protected Member Functions inherited from edm::ESProducer
ESConsumesInfo *	consumesInfoPushBackNew ()

unsigned int	consumesInfoSize () const

template<typename CallbackT , typename TList , typename TRecord >
void	registerProducts (std::shared_ptr< std::pair< unsigned int, std::shared_ptr< CallbackT >>> iCallback, const TList , const TRecord iRecord, const es::Label &iLabel)

template<typename CallbackT , typename TRecord >
void	registerProducts (std::shared_ptr< std::pair< unsigned int, std::shared_ptr< CallbackT >>>, const eventsetup::produce::Null , const TRecord , const es::Label &)

template<typename T >
auto	setWhatProduced (T *iThis, const es::Label &iLabel={})

template<typename T >
auto	setWhatProduced (T iThis, const char iLabel)

template<typename T >
auto	setWhatProduced (T *iThis, const std::string &iLabel)

template<typename T , typename TDecorator >
auto	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel={})

template<typename T , typename TReturn , typename TRecord >
auto	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const es::Label &iLabel={})

template<typename T , typename TReturn , typename TRecord , typename TDecorator >
auto	setWhatProduced (T iThis, TReturn(T ::iMethod)(const TRecord &), const TDecorator &iDec, const es::Label &iLabel={})

template<typename TFunc >
auto	setWhatProduced (TFunc &&func, const es::Label &iLabel={})

template<typename TReturn , typename TRecord , typename TFunc , typename TDecorator >
ESConsumesCollectorT< TRecord >	setWhatProduced (TFunc &&func, TDecorator &&iDec, const es::Label &iLabel={})

void	usesResources (std::vector< std::string > const &)

Protected Member Functions inherited from edm::ESProductResolverFactoryProducer
template<class TFactory >
void	registerFactory (std::unique_ptr< TFactory > iFactory, const std::string &iLabel=std::string())

virtual void	registerFactoryWithKey (const EventSetupRecordKey &iRecord, std::unique_ptr< eventsetup::ESProductResolverFactoryBase > iFactory, const std::string &iLabel=std::string())

KeyedResolversVector	registerResolvers (const EventSetupRecordKey &, unsigned int iovIndex) override

Protected Member Functions inherited from edm::eventsetup::ESProductResolverProvider
template<class T >
void	usingRecord ()

void	usingRecordWithKey (const EventSetupRecordKey &key)

Protected Member Functions inherited from edm::EventSetupRecordIntervalFinder
template<class T >
void	findingRecord ()

void	findingRecordWithKey (const eventsetup::EventSetupRecordKey &)

Detailed Description

Description: Generator of the ideal/fake conditions for the LorentzAngle.

It receives input values with layer granularity and it is able to perform gaussian smearing or use a uniform distribution at the module level. Depending on the parameters passed via cfg, it is able to generate the values per DetId with a gaussian distribution and a uniform distribution. When setting the sigma of the gaussian to 0 and passing a single value the generated values are fixed. For TID and TEC the decision to generate with a uniform distribution comes from the setting for the first layers of TIB and TOB.

Implementation: Port of SiStripLorentzAngleGenerator and templated fake ESSource to an edm::ESProducer

Definition at line 38 of file SiStripLorentzAngleFakeESSource.cc.

Member Typedef Documentation

◆ ReturnType

typedef std::unique_ptr<SiStripLorentzAngle> SiStripLorentzAngleFakeESSource::ReturnType

Definition at line 47 of file SiStripLorentzAngleFakeESSource.cc.

Constructor & Destructor Documentation

◆ SiStripLorentzAngleFakeESSource()

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

Definition at line 114 of file SiStripLorentzAngleFakeESSource.cc.

References gpuPixelDoublets::cc, gather_cfg::cout, edm::ParameterSet::getParameter(), m_geomDetToken, m_StdDevs_TIB, m_StdDevs_TOB, m_TIB_EstimatedValuesMax, m_TIB_EstimatedValuesMin, m_TIB_PerCent_Errs, m_TIBmeanPerCentError, m_TIBmeanStdDev, m_TIBmeanValueMax, m_TIBmeanValueMin, m_TOB_EstimatedValuesMax, m_TOB_EstimatedValuesMin, m_TOB_PerCent_Errs, m_TOBmeanPerCentError, m_TOBmeanStdDev, m_TOBmeanValueMax, m_TOBmeanValueMin, m_tTopoToken, m_uniformTIB, m_uniformTOB, edm::ESProducer::setWhatProduced(), and HcalDetIdTransform::transform().

                                                                                                {
   auto cc = setWhatProduced(this);
   m_tTopoToken = cc.consumes();
   m_geomDetToken = cc.consumes();
   findingRecord<SiStripLorentzAngleRcd>();
 
   m_TIB_EstimatedValuesMin = iConfig.getParameter<std::vector<double>>("TIB_EstimatedValuesMin");
   m_TIB_EstimatedValuesMax = iConfig.getParameter<std::vector<double>>("TIB_EstimatedValuesMax");
   m_TOB_EstimatedValuesMin = iConfig.getParameter<std::vector<double>>("TOB_EstimatedValuesMin");
   m_TOB_EstimatedValuesMax = iConfig.getParameter<std::vector<double>>("TOB_EstimatedValuesMax");
   m_TIB_PerCent_Errs = iConfig.getParameter<std::vector<double>>("TIB_PerCent_Errs");
   m_TOB_PerCent_Errs = iConfig.getParameter<std::vector<double>>("TOB_PerCent_Errs");
 
   // If max values are passed they must be equal in number to the min values.
   if (((!m_TIB_EstimatedValuesMax.empty()) && (m_TIB_EstimatedValuesMin.size() != m_TIB_EstimatedValuesMax.size())) ||
       ((!m_TOB_EstimatedValuesMax.empty()) && (m_TOB_EstimatedValuesMin.size() != m_TOB_EstimatedValuesMax.size()))) {
     std::cout << "ERROR: size of min and max values is different" << std::endl;
     std::cout << "TIB_EstimatedValuesMin.size() = " << m_TIB_EstimatedValuesMin.size()
               << ", TIB_EstimatedValuesMax.size() " << m_TIB_EstimatedValuesMax.size() << std::endl;
     std::cout << "TOB_EstimatedValuesMin.size() = " << m_TOB_EstimatedValuesMin.size()
               << ", TOB_EstimatedValuesMax.size() " << m_TOB_EstimatedValuesMax.size() << std::endl;
   }
 
   m_uniformTIB = std::vector<bool>(m_TIB_EstimatedValuesMin.size(), false);
   m_uniformTOB = std::vector<bool>(m_TOB_EstimatedValuesMin.size(), false);
   setUniform(m_TIB_EstimatedValuesMin, m_TIB_EstimatedValuesMax, m_uniformTIB);
   setUniform(m_TOB_EstimatedValuesMin, m_TOB_EstimatedValuesMax, m_uniformTOB);
 
   // Compute standard deviations
   m_StdDevs_TIB = std::vector<double>(m_TIB_EstimatedValuesMin.size(), 0);
   m_StdDevs_TOB = std::vector<double>(m_TOB_EstimatedValuesMin.size(), 0);
   transform(m_TIB_EstimatedValuesMin.begin(),
             m_TIB_EstimatedValuesMin.end(),
             m_TIB_PerCent_Errs.begin(),
             m_StdDevs_TIB.begin(),
             computeSigma);
   transform(m_TOB_EstimatedValuesMin.begin(),
             m_TOB_EstimatedValuesMin.end(),
             m_TOB_PerCent_Errs.begin(),
             m_StdDevs_TOB.begin(),
             computeSigma);
 
   // Compute mean values to be used with TID and TEC
   m_TIBmeanValueMin = std::accumulate(m_TIB_EstimatedValuesMin.begin(), m_TIB_EstimatedValuesMin.end(), 0.) /
                       double(m_TIB_EstimatedValuesMin.size());
   m_TIBmeanValueMax = std::accumulate(m_TIB_EstimatedValuesMax.begin(), m_TIB_EstimatedValuesMax.end(), 0.) /
                       double(m_TIB_EstimatedValuesMax.size());
   m_TOBmeanValueMin = std::accumulate(m_TOB_EstimatedValuesMin.begin(), m_TOB_EstimatedValuesMin.end(), 0.) /
                       double(m_TOB_EstimatedValuesMin.size());
   m_TOBmeanValueMax = std::accumulate(m_TOB_EstimatedValuesMax.begin(), m_TOB_EstimatedValuesMax.end(), 0.) /
                       double(m_TOB_EstimatedValuesMax.size());
   m_TIBmeanPerCentError =
       std::accumulate(m_TIB_PerCent_Errs.begin(), m_TIB_PerCent_Errs.end(), 0.) / double(m_TIB_PerCent_Errs.size());
   m_TOBmeanPerCentError =
       std::accumulate(m_TOB_PerCent_Errs.begin(), m_TOB_PerCent_Errs.end(), 0.) / double(m_TOB_PerCent_Errs.size());
   m_TIBmeanStdDev = (m_TIBmeanPerCentError / 100) * m_TIBmeanValueMin;
   m_TOBmeanStdDev = (m_TOBmeanPerCentError / 100) * m_TOBmeanValueMin;
 }

◆ ~SiStripLorentzAngleFakeESSource()

SiStripLorentzAngleFakeESSource::~SiStripLorentzAngleFakeESSource ( )

override

Definition at line 173 of file SiStripLorentzAngleFakeESSource.cc.

173 {}

Member Function Documentation

◆ produce()

SiStripLorentzAngleFakeESSource::ReturnType SiStripLorentzAngleFakeESSource::produce ( const SiStripLorentzAngleRcd & iRecord )

Definition at line 182 of file SiStripLorentzAngleFakeESSource.cc.

References hcalRecHitTable_cff::detId, edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::get(), edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::getRecord(), TrackerGeometryUtils::getSiStripDetIds(), m_geomDetToken, m_StdDevs_TIB, m_StdDevs_TOB, m_TIB_EstimatedValuesMax, m_TIB_EstimatedValuesMin, m_TIBmeanStdDev, m_TIBmeanValueMax, m_TIBmeanValueMin, m_TOB_EstimatedValuesMax, m_TOB_EstimatedValuesMin, m_TOBmeanStdDev, m_TOBmeanValueMax, m_TOBmeanValueMin, m_tTopoToken, m_uniformTIB, m_uniformTOB, DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, and StripSubdetector::TOB.

                                            {
   using namespace edm::es;
 
   const auto& geomDet = iRecord.getRecord<TrackerTopologyRcd>().get(m_geomDetToken);
   const auto& tTopo = iRecord.get(m_tTopoToken);
 
   auto lorentzAngle = std::make_unique<SiStripLorentzAngle>();
 
   for (const auto detId : TrackerGeometryUtils::getSiStripDetIds(geomDet)) {
     const DetId detectorId = DetId(detId);
     const int subDet = detectorId.subdetId();
 
     float mobi{0.};
 
     if (subDet == int(StripSubdetector::TIB)) {
       const int layerId = tTopo.tibLayer(detectorId) - 1;
       mobi = hallMobility(m_TIB_EstimatedValuesMin[layerId],
                           m_TIB_EstimatedValuesMax[layerId],
                           m_StdDevs_TIB[layerId],
                           m_uniformTIB[layerId]);
     } else if (subDet == int(StripSubdetector::TOB)) {
       const int layerId = tTopo.tobLayer(detectorId) - 1;
       mobi = hallMobility(m_TOB_EstimatedValuesMin[layerId],
                           m_TOB_EstimatedValuesMax[layerId],
                           m_StdDevs_TOB[layerId],
                           m_uniformTOB[layerId]);
     } else if (subDet == int(StripSubdetector::TID)) {
       // ATTENTION: as of now the uniform generation for TID is decided by the setting for layer 0 of TIB
       mobi = hallMobility(m_TIBmeanValueMin, m_TIBmeanValueMax, m_TIBmeanStdDev, m_uniformTIB[0]);
     }
     if (subDet == int(StripSubdetector::TEC)) {
       if (tTopo.tecRing(detectorId) < 5) {
         // ATTENTION: as of now the uniform generation for TEC is decided by the setting for layer 0 of TIB
         mobi = hallMobility(m_TIBmeanValueMin, m_TIBmeanValueMax, m_TIBmeanStdDev, m_uniformTIB[0]);
       } else {
         // ATTENTION: as of now the uniform generation for TEC is decided by the setting for layer 0 of TOB
         mobi = hallMobility(m_TOBmeanValueMin, m_TOBmeanValueMax, m_TOBmeanStdDev, m_uniformTOB[0]);
       }
     }
 
     if (!lorentzAngle->putLorentzAngle(detId, mobi)) {
       edm::LogError("SiStripLorentzAngleFakeESSource::produce ") << " detid already exists";
     }
   }
 
   return lorentzAngle;
 }

◆ setIntervalFor()

void SiStripLorentzAngleFakeESSource::setIntervalFor	(	const edm::eventsetup::EventSetupRecordKey &	,
		const edm::IOVSyncValue &	iov,
		edm::ValidityInterval &	iValidity
	)

overridevirtual

Implements edm::EventSetupRecordIntervalFinder.

Definition at line 175 of file SiStripLorentzAngleFakeESSource.cc.

References edm::IOVSyncValue::beginOfTime(), and edm::IOVSyncValue::endOfTime().

                                                                                      {
   iValidity = edm::ValidityInterval{iov.beginOfTime(), iov.endOfTime()};
 }