PPSModifySingularModes Class Reference

Modifies the alignment modes unconstrained by the track-based alignment. More...

Inheritance diagram for PPSModifySingularModes:

Public Member Functions
	PPSModifySingularModes (const edm::ParameterSet &ps)
Public Member Functions inherited from edm::stream::EDAnalyzer<>
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
	EDAnalyzerBase (const EDAnalyzerBase &)=delete
ModuleDescription const &	moduleDescription () const
const EDAnalyzerBase &	operator= (const EDAnalyzerBase &)=delete
	~EDAnalyzerBase () 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 const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
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	analyze (const edm::Event &e, const edm::EventSetup &es) override
void	beginRun (edm::Run const &, edm::EventSetup const &) override

Private Attributes
edm::ParameterSet	ps_
edm::ESGetToken< CTPPSGeometry, VeryForwardRealGeometryRecord >	tokenRealGeometry_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
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)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Modifies the alignment modes unconstrained by the track-based alignment.

Definition at line 27 of file PPSModifySingularModes.cc.

Constructor & Destructor Documentation

PPSModifySingularModes()

PPSModifySingularModes::PPSModifySingularModes

(

const edm::ParameterSet &

ps

)

Definition at line 46 of file PPSModifySingularModes.cc.

References edm::BeginRun.

    : ps_(ps), tokenRealGeometry_(esConsumes<edm::Transition::BeginRun>()) {}

Member Function Documentation

analyze()

void PPSModifySingularModes::analyze ( const edm::Event &  e, const edm::EventSetup &  es  )

inlineoverrideprivatevirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 38 of file PPSModifySingularModes.cc.

{}

beginRun()

void PPSModifySingularModes::beginRun ( edm::Run const &  , edm::EventSetup const &  es  )

overrideprivatevirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 51 of file PPSModifySingularModes.cc.

References CTPPSDetId::arm(), AlignmentTask::buildGeometry(), HltBtagPostValidation_cff::c, ztail::d, Exception, factorRPFromSensorCorrections(), edm::EventSetup::getData(), edm::ParameterSet::getUntrackedParameter(), input, makeListRunsInFiles::inputFile, CTPPSRPAlignmentCorrectionsMethods::loadFromXML(), convertSQLitetoXML_cfg::output, download_sqlite_cfg::outputFile, ps_, CTPPSDetId::rp(), CTPPSDetId::station(), tokenRealGeometry_, CTPPSRPAlignmentCorrectionsMethods::writeToXML(), z, and testProducerWithPsetDescEmpty_cfi::z2.

                                                                             {
  // get config parameters
  const double z1 = ps_.getUntrackedParameter<double>("z1");
  const double z2 = ps_.getUntrackedParameter<double>("z2");
  const double de_x1 = ps_.getUntrackedParameter<double>("de_x1");
  const double de_x2 = ps_.getUntrackedParameter<double>("de_x2");
  const double de_y1 = ps_.getUntrackedParameter<double>("de_y1");
  const double de_y2 = ps_.getUntrackedParameter<double>("de_y2");
  const double de_rho1 = ps_.getUntrackedParameter<double>("de_rho1");
  const double de_rho2 = ps_.getUntrackedParameter<double>("de_rho2");

  FileInPath inputFileInPath(ps_.getUntrackedParameter<string>("inputFile"));
  const string inputFile = inputFileInPath.fullPath();
  const string outputFile = ps_.getUntrackedParameter<string>("outputFile");

  // validate config parameters
  if (z1 == z2)
    throw cms::Exception("PPS") << "z1 equals z2";

  // calculate slopes and intercepts
  const double a_x = (de_x2 - de_x1) / (z2 - z1), b_x = de_x1 - a_x * z1;
  const double a_y = (de_y2 - de_y1) / (z2 - z1), b_y = de_y1 - a_y * z1;
  const double a_rho = (de_rho2 - de_rho1) / (z2 - z1), b_rho = de_rho1 - a_rho * z1;

  // get geometry
  const auto &geometry = es.getData(tokenRealGeometry_);

  // get input alignments
  CTPPSRPAlignmentCorrectionsDataSequence inputSequence = CTPPSRPAlignmentCorrectionsMethods::loadFromXML(inputFile);
  const auto &input = inputSequence.begin()->second;

  // modify the singular modes
  CTPPSRPAlignmentCorrectionsData output = input;

  for (auto &it : input.getSensorMap()) {
    const auto &sensorId = it.first;

    const auto &c = geometry.sensorTranslation(sensorId);

    // pixels cannot be described by one single value of z, but no better approxiamtion can be made
    const double z = c.z();

    double de_ShX = a_x * z + b_x;
    double de_ShY = a_y * z + b_y;
    const double de_RotZ = a_rho * z + b_rho;

    // add the effect of global rotation (about origin, not sensor centre)
    de_ShX -= +de_RotZ * (c.y() + de_ShY);
    de_ShY -= -de_RotZ * (c.x() + de_ShX);

    CTPPSRPAlignmentCorrectionData d = it.second;
    d.setShX(d.getShX() + de_ShX);
    d.setShY(d.getShY() + de_ShY);
    d.setRotZ(d.getRotZ() + de_RotZ);

    output.setSensorCorrection(sensorId, d);
  }

  // build list of RPs
  vector<unsigned int> rps;
  unsigned int last_rp = 123456;
  for (auto &it : input.getSensorMap()) {
    CTPPSDetId senId(it.first);
    unsigned int rpDecId = senId.arm() * 100 + senId.station() * 10 + senId.rp();

    if (last_rp != rpDecId) {
      rps.push_back(rpDecId);
      last_rp = rpDecId;
    }
  }

  // build alignment geometry (needed for the factorisation below)
  AlignmentGeometry alignmentGeometry;
  vector<unsigned int> excludePlanes;
  AlignmentTask::buildGeometry(rps, excludePlanes, &geometry, 0., alignmentGeometry);

  // factorise output
  CTPPSRPAlignmentCorrectionsData outputExpanded;
  CTPPSRPAlignmentCorrectionsData outputFactored;

  const bool equalWeights = false;
  factorRPFromSensorCorrections(output, outputExpanded, outputFactored, alignmentGeometry, equalWeights, 1);

  // save output
  CTPPSRPAlignmentCorrectionsMethods::writeToXML(outputFactored, outputFile, false, false, true, true, true, true);
}

Member Data Documentation

ps_

edm::ParameterSet PPSModifySingularModes::ps_

private

Definition at line 32 of file PPSModifySingularModes.cc.

Referenced by beginRun().

tokenRealGeometry_

edm::ESGetToken<CTPPSGeometry, VeryForwardRealGeometryRecord> PPSModifySingularModes::tokenRealGeometry_

private

Definition at line 34 of file PPSModifySingularModes.cc.

Referenced by beginRun().