Inheritance diagram for CTPPSModifiedOpticalFunctionsESSource:

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

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

	~CTPPSModifiedOpticalFunctionsESSource () 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
double	factor_

edm::ESGetToken< LHCInterpolatedOpticalFunctionsSetCollection, CTPPSInterpolatedOpticsRcd >	inputOpticsToken_

unsigned int	rpDecId_45_F_

unsigned int	rpDecId_45_N_

unsigned int	rpDecId_56_F_

unsigned int	rpDecId_56_N_

std::string	scenario_

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 16 of file CTPPSModifiedOpticalFunctionsESSource.cc.

Constructor & Destructor Documentation

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

Definition at line 38 of file CTPPSModifiedOpticalFunctionsESSource.cc.

References edm::ParameterSet::getParameter(), inputOpticsToken_, edm::ESProducer::setWhatProduced(), and AlCaHLTBitMon_QueryRunRegistry::string.

     : scenario_(iConfig.getParameter<std::string>("scenario")),
       factor_(iConfig.getParameter<double>("factor")),
       rpDecId_45_N_(iConfig.getParameter<unsigned int>("rpId_45_N")),
       rpDecId_45_F_(iConfig.getParameter<unsigned int>("rpId_45_F")),
       rpDecId_56_N_(iConfig.getParameter<unsigned int>("rpId_56_N")),
       rpDecId_56_F_(iConfig.getParameter<unsigned int>("rpId_56_F"))
 {
   auto cc = setWhatProduced(this, iConfig.getParameter<std::string>("outputOpticsLabel"));
 
   inputOpticsToken_ = cc.consumesFrom<LHCInterpolatedOpticalFunctionsSetCollection, CTPPSInterpolatedOpticsRcd>
     (edm::ESInputTag("", iConfig.getParameter<std::string>("inputOpticsLabel")));
 }

CTPPSModifiedOpticalFunctionsESSource::~CTPPSModifiedOpticalFunctionsESSource ( )

inlineoverride

Definition at line 19 of file CTPPSModifiedOpticalFunctionsESSource.cc.

References fillDescriptions(), and produce().

19 {}

Member Function Documentation

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

static

Definition at line 54 of file CTPPSModifiedOpticalFunctionsESSource.cc.

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

Referenced by ~CTPPSModifiedOpticalFunctionsESSource().

                                                                                                        {
   edm::ParameterSetDescription desc;
 
   desc.add<std::string>("inputOpticsLabel", "")->setComment("label of the input optics records");
   desc.add<std::string>("outputOpticsLabel", "modified")->setComment("label of the output optics records");
 
   desc.add<std::string>("scenario", "none")->setComment("name of modification scenario");
 
   desc.add<double>("factor", 0.)->setComment("size of modification (number of sigmas)");
 
   desc.add<unsigned int>("rpId_45_N", 0)->setComment("decimal RP id for 45 near");
   desc.add<unsigned int>("rpId_45_F", 0)->setComment("decimal RP id for 45 far");
   desc.add<unsigned int>("rpId_56_N", 0)->setComment("decimal RP id for 56 near");
   desc.add<unsigned int>("rpId_56_F", 0)->setComment("decimal RP id for 56 far");
 
   descriptions.add("ctppsModifiedOpticalFunctionsESSource", desc);
 }

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

Definition at line 74 of file CTPPSModifiedOpticalFunctionsESSource.cc.

References a, protons_cff::arm, CTPPSDetId::arm(), b, edmIntegrityCheck::d, DEFINE_FWK_EVENTSETUP_MODULE, LHCOpticalFunctionsSet::eLpx, LHCOpticalFunctionsSet::eLx, LHCOpticalFunctionsSet::exd, factor_, edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::get(), mps_fire::i, input, inputOpticsToken_, convertSQLitetoXML_cfg::output, CTPPSDetId::rp(), rpDecId_45_F_, rpDecId_45_N_, rpDecId_56_F_, rpDecId_56_N_, year_2016_postTS2_cff::rpId, scenario_, CTPPSDetId::station(), and protons_cff::xi.

Referenced by ~CTPPSModifiedOpticalFunctionsESSource().

                                                {
   // get input
   LHCInterpolatedOpticalFunctionsSetCollection const &input = iRecord.get(inputOpticsToken_);
 
   // prepare output
   std::shared_ptr<LHCInterpolatedOpticalFunctionsSetCollection> output =
       std::make_shared<LHCInterpolatedOpticalFunctionsSetCollection>(input);
 
   // premare arm/RP id map
   struct ArmInfo {
     unsigned int rpId_N = 0, rpId_F = 0;
   };
 
   std::map<unsigned int, ArmInfo> armInfo;
 
   for (const auto &fsp : *output) {
     CTPPSDetId rpId(fsp.first);
     unsigned int rpDecId = 100 * rpId.arm() + 10 * rpId.station() + rpId.rp();
 
     if (rpDecId == rpDecId_45_N_)
       armInfo[0].rpId_N = fsp.first;
     if (rpDecId == rpDecId_45_F_)
       armInfo[0].rpId_F = fsp.first;
     if (rpDecId == rpDecId_56_N_)
       armInfo[1].rpId_N = fsp.first;
     if (rpDecId == rpDecId_56_F_)
       armInfo[1].rpId_F = fsp.first;
   }
 
   // loop over arms
   bool applied = false;
 
   for (const auto &ap : armInfo) {
     const auto &arm = ap.first;
 
     //printf("* arm %u\n", arm);
 
     auto &of_N = output->find(ap.second.rpId_N)->second;
     auto &of_F = output->find(ap.second.rpId_F)->second;
 
     const double z_N = of_N.getScoringPlaneZ();
     const double z_F = of_F.getScoringPlaneZ();
     const double de_z = (arm == 0) ? z_N - z_F : z_F - z_N;
 
     //printf("  z_N = %.3f m, z_F = %.3f m\n", z_N*1E-2, z_F*1E-2);
 
     if (of_N.m_xi_values.size() != of_N.m_xi_values.size())
       throw cms::Exception("CTPPSModifiedOpticalFunctionsESSource")
           << "Different xi sampling of optical functions in near and far RP.";
 
     // loop over sampling points (in xi)
     for (unsigned int i = 0; i < of_N.m_xi_values.size(); ++i) {
       const double xi = of_N.m_xi_values[i];
 
       double x_d_N = of_N.m_fcn_values[LHCOpticalFunctionsSet::exd][i];
       double x_d_F = of_F.m_fcn_values[LHCOpticalFunctionsSet::exd][i];
 
       double L_x_N = of_N.m_fcn_values[LHCOpticalFunctionsSet::eLx][i];
       double L_x_F = of_F.m_fcn_values[LHCOpticalFunctionsSet::eLx][i];
       double Lp_x = of_N.m_fcn_values[LHCOpticalFunctionsSet::eLpx][i];
 
       //printf("  xi = %.3f, Lp_x = %.3f, %.3f\n", xi, Lp_x, (L_x_F - L_x_N) / de_z);
 
       // apply modification scenario
       if (scenario_ == "none")
         applied = true;
 
       if (scenario_ == "Lx") {
         const double a = 3180., b = 40.;  // cm
         const double de_L_x = factor_ * (a * xi + b);
         L_x_N += de_L_x;
         L_x_F += de_L_x;
         applied = true;
       }
 
       if (scenario_ == "Lpx") {
         const double a = 0.42, b = 0.015;  // dimensionless
         const double de_Lp_x = factor_ * (a * xi + b) * Lp_x;
         Lp_x += de_Lp_x;
         L_x_N -= de_Lp_x * de_z / 2.;
         L_x_F += de_Lp_x * de_z / 2.;
         applied = true;
       }
 
       if (scenario_ == "xd") {
         const double d = 0.08;  // dimensionless
         x_d_N += x_d_N * d * factor_;
         x_d_F += x_d_F * d * factor_;
         applied = true;
       }
 
       // TODO: for test only
       if (scenario_ == "Lx-scale") {
         L_x_N *= factor_;
         L_x_F *= factor_;
         applied = true;
       }
 
       // store updated values
       of_N.m_fcn_values[LHCOpticalFunctionsSet::exd][i] = x_d_N;
       of_F.m_fcn_values[LHCOpticalFunctionsSet::exd][i] = x_d_F;
 
       of_N.m_fcn_values[LHCOpticalFunctionsSet::eLx][i] = L_x_N;
       of_F.m_fcn_values[LHCOpticalFunctionsSet::eLx][i] = L_x_F;
 
       of_N.m_fcn_values[LHCOpticalFunctionsSet::eLpx][i] = Lp_x;
       of_F.m_fcn_values[LHCOpticalFunctionsSet::eLpx][i] = Lp_x;
     }
 
     // re-initialise splines
     of_N.initializeSplines();
     of_F.initializeSplines();
   }
 
   // modification applied?
   if (!applied)
     edm::LogError("CTPPSModifiedOpticalFunctionsESSource") << "Could not apply scenario `" + scenario_ + "'.";
 
   // save modified output
   return output;
 }