Inheritance diagram for CTPPSInterpolatedOpticalFunctionsESSource:

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

std::shared_ptr< LHCInterpolatedOpticalFunctionsSetCollection >	produce (const CTPPSInterpolatedOpticsRcd &)

	~CTPPSInterpolatedOpticalFunctionsESSource () override

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

ESProxyIndex const *	getTokenIndices (unsigned int iIndex) const

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &) final

	~ESProducer ()(false) override

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

void	newInterval (const eventsetup::EventSetupRecordKey &iRecordType, const ValidityInterval &iInterval) override
	overrides DataProxyProvider method More...

	~ESProxyFactoryProducer () noexcept(false) override

Public Member Functions inherited from edm::eventsetup::DataProxyProvider
	DataProxyProvider ()

const ComponentDescription &	description () const

bool	isUsingRecord (const EventSetupRecordKey &) const

const KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey) const

void	resetProxies (const EventSetupRecordKey &iRecordType)

void	resetProxiesIfTransient (const EventSetupRecordKey &iRecordType)

void	setAppendToDataLabel (const edm::ParameterSet &)

void	setDescription (const ComponentDescription &iDescription)

std::set< EventSetupRecordKey >	usingRecords () const

virtual	~DataProxyProvider () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)

Private Attributes
float	currentCrossingAngle_

std::shared_ptr< LHCInterpolatedOpticalFunctionsSetCollection >	currentData_

bool	currentDataValid_

std::string	lhcInfoLabel_

std::string	opticsLabel_

Additional Inherited Members
Public Types inherited from edm::eventsetup::DataProxyProvider
typedef std::vector< std::pair< DataKey, edm::propagate_const< std::shared_ptr< DataProxy > > > >	KeyedProxies

typedef std::vector< EventSetupRecordKey >	Keys

typedef std::map< EventSetupRecordKey, KeyedProxies >	RecordProxies

Protected Member Functions inherited from edm::ESProducer
	ESProducer (const ESProducer &)=delete

ESProducer const &	operator= (const ESProducer &)=delete

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 TArg >
ESConsumesCollectorT< TRecord >	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel={})

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

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

void	registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList) override
	override DataProxyProvider method More...

Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
void	eraseAll (const EventSetupRecordKey &iRecordKey)
	deletes all the Proxies in aStream More...

void	invalidateProxies (const EventSetupRecordKey &iRecordKey)

template<class T >
void	usingRecord ()

void	usingRecordWithKey (const EventSetupRecordKey &)

Detailed Description

Definition at line 20 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

Constructor & Destructor Documentation

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

Definition at line 42 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

References opticsLabel_, and edm::ESProducer::setWhatProduced().

                                                                                                                    :
   lhcInfoLabel_(iConfig.getParameter<std::string>("lhcInfoLabel")),
   opticsLabel_(iConfig.getParameter<std::string>("opticsLabel")),
   currentCrossingAngle_(-1.),
   currentDataValid_(false)
 {
   setWhatProduced(this, opticsLabel_);
 }

CTPPSInterpolatedOpticalFunctionsESSource::~CTPPSInterpolatedOpticalFunctionsESSource ( )

inlineoverride

Definition at line 24 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

References fillDescriptions(), and produce().

24 {}

Member Function Documentation

void CTPPSInterpolatedOpticalFunctionsESSource::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 53 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~CTPPSInterpolatedOpticalFunctionsESSource().

 {
   edm::ParameterSetDescription desc;
 
   desc.add<std::string>("lhcInfoLabel", "")
     ->setComment("label of the LHCInfo record");
 
   desc.add<std::string>("opticsLabel", "")
     ->setComment("label of the optics record");
 
   descriptions.add("ctppsInterpolatedOpticalFunctionsESSource", desc);
 }

std::shared_ptr< LHCInterpolatedOpticalFunctionsSetCollection > CTPPSInterpolatedOpticalFunctionsESSource::produce ( const CTPPSInterpolatedOpticsRcd & iRecord )

Definition at line 68 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

References LHCInfo::crossingAngle(), currentCrossingAngle_, currentData_, currentDataValid_, DEFINE_FWK_EVENTSETUP_MODULE, MillePedeFileConverter_cfg::e, Exception, plotBeamSpotDB::first, edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::getRecord(), LHCInterpolatedOpticalFunctionsSet::initializeSplines(), lhcInfoLabel_, LHCOpticalFunctionsSet::m_fcn_values, LHCOpticalFunctionsSet::m_xi_values, LHCOpticalFunctionsSet::m_z, eostools::move(), GetRecoTauVFromDQM_MC_cff::next, LHCOpticalFunctionsSet::nFunctions, opticsLabel_, pi, year_2016_postTS2_cff::rpId, and protons_cff::xi.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), HTMLExport.HTMLExportStatic::export(), and ~CTPPSInterpolatedOpticalFunctionsESSource().

 {
   // get the input data
   edm::ESHandle<LHCOpticalFunctionsSetCollection> hOFColl;
   iRecord.getRecord<CTPPSOpticsRcd>().get(opticsLabel_, hOFColl);
 
   edm::ESHandle<LHCInfo> hLHCInfo;
   iRecord.getRecord<LHCInfoRcd>().get(lhcInfoLabel_, hLHCInfo);
 
   // is there anything to do?
   if (currentDataValid_ && hLHCInfo->crossingAngle() == currentCrossingAngle_)
     return currentData_;
 
   // is crossing angle reasonable (LHCInfo is correctly filled in DB)?
   if (hLHCInfo->crossingAngle() == 0.)
   {
     edm::LogInfo("CTPPSInterpolatedOpticalFunctionsESSource") << "Invalid crossing angle, no optical functions produced.";
 
     currentDataValid_ = false;
     currentCrossingAngle_ = -1;
     currentData_ = std::make_shared<LHCInterpolatedOpticalFunctionsSetCollection>();
 
     return currentData_;
   }
 
   // set new crossing angle
   currentCrossingAngle_ = hLHCInfo->crossingAngle();
   edm::LogInfo("CTPPSInterpolatedOpticalFunctionsESSource") << "Crossing angle has changed to " << currentCrossingAngle_ << ".";
 
   // is input optics available ?
   if (hOFColl->empty())
   {
     edm::LogInfo("CTPPSInterpolatedOpticalFunctionsESSource") << "No input optics available, no optical functions produced.";
 
     currentDataValid_ = false;
     currentData_ = std::make_shared<LHCInterpolatedOpticalFunctionsSetCollection>();
 
     return currentData_;
   }
 
   // regular case with single-xangle input
   if (hOFColl->size() == 1)
   {
     const auto &it = hOFColl->begin();
 
     // does the input xangle correspond to the actual one?
     if (fabs(currentCrossingAngle_ - it->first) > 1e-6)
       throw cms::Exception("CTPPSInterpolatedOpticalFunctionsESSource") << "Cannot interpolate: input given only for xangle "
         << it->first << " while interpolation requested for " << currentCrossingAngle_ << ".";
 
     currentData_ = std::make_shared<LHCInterpolatedOpticalFunctionsSetCollection>();
     for (const auto &rp_p : it->second)
     {
       const auto rpId = rp_p.first;
       LHCInterpolatedOpticalFunctionsSet iof(rp_p.second);
       iof.initializeSplines();
       currentData_->emplace(rpId, std::move(iof));
     }
 
     currentDataValid_ = true;
   }
 
   // regular case with multi-xangle input
   if (hOFColl->size() > 1)
   {
     // find the closest xangle points for interpolation
     auto it1 = hOFColl->begin();
     auto it2 = std::next(it1);
 
     if (currentCrossingAngle_ > it1->first)
     {
       for (; it1 != hOFColl->end(); ++it1)
       {
         it2 = std::next(it1);
 
         if (it2 == hOFColl->end())
         {
           it2 = it1;
           it1 = std::prev(it1);
           break;
         }
 
         if (it1->first <= currentCrossingAngle_ && currentCrossingAngle_ < it2->first)
           break;
       }
     }
 
     const auto &xangle1 = it1->first;
     const auto &xangle2 = it2->first;
 
     const auto &ofs1 = it1->second;
     const auto &ofs2 = it2->second;
 
     // do the interpoaltion RP by RP
     currentData_ = std::make_shared<LHCInterpolatedOpticalFunctionsSetCollection>();
     for (const auto &rp_p : ofs1)
     {
       const auto rpId = rp_p.first;
       const auto &rp_it2 = ofs2.find(rpId);
       if (rp_it2 == ofs2.end())
         throw cms::Exception("CTPPSInterpolatedOpticalFunctionsESSource") << "RP mismatch between ofs1 and ofs2.";
 
       const auto &of1 = rp_p.second;
       const auto &of2 = rp_it2->second;
 
       const size_t num_xi_vals1 = of1.getXiValues().size();
       const size_t num_xi_vals2 = of2.getXiValues().size();
 
       if (num_xi_vals1 != num_xi_vals2)
         throw cms::Exception("CTPPSInterpolatedOpticalFunctionsESSource") << "Size mismatch between ofs1 and ofs2.";
 
       const size_t num_xi_vals = num_xi_vals1;
 
       LHCInterpolatedOpticalFunctionsSet iof;
       iof.m_z = of1.getScoringPlaneZ();
       iof.m_fcn_values.resize(LHCInterpolatedOpticalFunctionsSet::nFunctions);
       iof.m_xi_values.resize(num_xi_vals);
 
       for (size_t fi = 0; fi < of1.getFcnValues().size(); ++fi)
       {
         iof.m_fcn_values[fi].resize(num_xi_vals);
 
         for (size_t pi = 0; pi < num_xi_vals; ++pi)
         {
           double xi = of1.getXiValues()[pi];
           double xi_control = of2.getXiValues()[pi];
 
           if (fabs(xi - xi_control) > 1e-6)
             throw cms::Exception("CTPPSInterpolatedOpticalFunctionsESSource") << "Xi mismatch between ofs1 and ofs2.";
 
           iof.m_xi_values[pi] = xi;
 
           double v1 = of1.getFcnValues()[fi][pi];
           double v2 = of2.getFcnValues()[fi][pi];
           iof.m_fcn_values[fi][pi] = v1 + (v2 - v1) / (xangle2 - xangle1) * (currentCrossingAngle_ - xangle1);
         }
       }
 
       iof.initializeSplines();
 
       currentDataValid_ = true;
       currentData_->emplace(rpId, std::move(iof));
     }
   }
 
   return currentData_;
 }