#include <Alignment/TrackerAlignment/plugins/MCMisalignmentScaler.cc>

Inheritance diagram for MCMisalignmentScaler:

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

	~MCMisalignmentScaler () override=default

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

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

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

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Types
using	ScalerMap = std::unordered_map< unsigned int, std::unordered_map< int, double > >

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

ScalerMap	decodeSubDetectors (const edm::VParameterSet &)

Private Attributes
const edm::ESGetToken< Alignments, TrackerAlignmentRcd >	aliToken_

bool	firstEvent_ {true}

const edm::ESGetToken< GeometricDet, IdealGeometryRecord >	geomDetToken_

const double	outlierPullToIdealCut_

const edm::ESGetToken< SiPixelQuality, SiPixelQualityRcd >	pixelQualityToken_

const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd >	ptitpToken_

const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd >	ptpToken_

const bool	pullBadModulesToIdeal_

const ScalerMap	scalers_

const edm::ESGetToken< SiStripQuality, SiStripQualityRcd >	stripQualityToken_

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	topoToken_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Description: Plugin to rescale misalignment wrt. ideal geometry

Implementation:

The plugin takes the ideal geometry and the alignment object and rescales the position difference by the scaling factor provided by the user.

Definition at line 62 of file MCMisalignmentScaler.cc.

Member Typedef Documentation

◆ ScalerMap

using MCMisalignmentScaler::ScalerMap = std::unordered_map<unsigned int, std::unordered_map<int, double> >

private

Definition at line 77 of file MCMisalignmentScaler.cc.

Constructor & Destructor Documentation

◆ MCMisalignmentScaler()

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

explicit

Definition at line 92 of file MCMisalignmentScaler.cc.

References decodeSubDetectors(), and edm::ParameterSet::getParameter().

     : pixelQualityToken_(esConsumes()),
       stripQualityToken_(esConsumes()),
       geomDetToken_(esConsumes()),
       ptpToken_(esConsumes()),
       ptitpToken_(esConsumes()),
       topoToken_(esConsumes()),
       aliToken_(esConsumes()),
       scalers_{decodeSubDetectors(iConfig.getParameter<edm::VParameterSet>("scalers"))},
       pullBadModulesToIdeal_{iConfig.getUntrackedParameter<bool>("pullBadModulesToIdeal")},
       outlierPullToIdealCut_{iConfig.getUntrackedParameter<double>("outlierPullToIdealCut")} {}

◆ ~MCMisalignmentScaler()

MCMisalignmentScaler::~MCMisalignmentScaler ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 109 of file MCMisalignmentScaler.cc.

References a, aliToken_, b, cond::TimeTypeSpecs::beginValue, TrackerGeomBuilderFromGeometricDet::build(), Exception, firstEvent_, geomDetToken_, edm::EventSetup::getData(), edm::Service< T >::isAvailable(), SiStripQuality::IsModuleBad(), SiPixelQuality::IsModuleBad(), Alignments::m_align, outlierPullToIdealCut_, PixelSubdetector::PixelBarrel, pixelQualityToken_, ptitpToken_, ptpToken_, pullBadModulesToIdeal_, AlignTransform::rawId(), cond::runnumber, scaleFactor, scalers_, TrackerTopology::side(), SiStripBadComponentsDQMServiceTemplate_cfg::since, jetsAK4_CHS_cff::sort, stripQualityToken_, GeomDetEnumerators::subDetId, StripSubdetector::TIB, TrackerTopology::tibSide(), cond::timeTypeSpecs, StripSubdetector::TOB, TrackerTopology::tobSide(), topoToken_, PbPb_ZMuSkimMuonDPG_cff::tracker, and cond::service::PoolDBOutputService::writeOneIOV().

                                                                                {
   if (!firstEvent_)
     return;
   firstEvent_ = false;
 
   // get handle on bad modules
   const SiPixelQuality* pixelModules = &iSetup.getData(pixelQualityToken_);
   const SiStripQuality* stripModules = &iSetup.getData(stripQualityToken_);
 
   // get the tracker geometry
   const GeometricDet* geometricDet = &iSetup.getData(geomDetToken_);
   const PTrackerParameters& ptp = iSetup.getData(ptpToken_);
   const PTrackerAdditionalParametersPerDet* ptitp = &iSetup.getData(ptitpToken_);
   const TrackerTopology* topology = &iSetup.getData(topoToken_);
 
   TrackerGeomBuilderFromGeometricDet trackerBuilder;
   auto tracker = std::unique_ptr<TrackerGeometry>{trackerBuilder.build(geometricDet, ptitp, ptp, topology)};
 
   auto dets = tracker->dets();
   std::sort(dets.begin(), dets.end(), [](const auto& a, const auto& b) {
     return a->geographicalId().rawId() < b->geographicalId().rawId();
   });
 
   // get the input alignment
   const Alignments* alignments = &iSetup.getData(aliToken_);
 
   if (dets.size() != alignments->m_align.size()) {
     throw cms::Exception("GeometryMismatch") << "Size mismatch between alignments (size=" << alignments->m_align.size()
                                              << ") and ideal geometry (size=" << dets.size() << ")";
   }
 
   Alignments rescaledAlignments{};
   {
     auto outlierCounter{0};
     auto ideal = dets.cbegin();
     const auto& ideal_end = dets.cend();
     auto misaligned = alignments->m_align.cbegin();
     for (; ideal != ideal_end; ++ideal, ++misaligned) {
       if ((*ideal)->geographicalId().rawId() != misaligned->rawId()) {
         throw cms::Exception("GeometryMismatch") << "Order differs between Dets in alignments ideal geometry.";
       }
 
       // determine scale factor
       const auto& subDetId = (*ideal)->geographicalId().subdetId();
       auto side = topology->side((*ideal)->geographicalId());
       if (side == 0) {
         switch (subDetId) {
           case PixelSubdetector::PixelBarrel:
             side = 1;  // both sides are treated identical -> pick one of them
             break;
           case StripSubdetector::TIB:
             side = topology->tibSide((*ideal)->geographicalId());
             break;
           case StripSubdetector::TOB:
             side = topology->tobSide((*ideal)->geographicalId());
             break;
           default:
             break;
         }
       }
       auto scaleFactor = scalers_.find(subDetId)->second.find(side)->second;
 
       if (pullBadModulesToIdeal_ &&
           (pixelModules->IsModuleBad(misaligned->rawId()) || stripModules->IsModuleBad(misaligned->rawId()))) {
         scaleFactor = 0.0;
       }
 
       auto x_diff = misaligned->translation().x() - (*ideal)->position().x();
       auto y_diff = misaligned->translation().y() - (*ideal)->position().y();
       auto z_diff = misaligned->translation().z() - (*ideal)->position().z();
 
       auto xx_diff = misaligned->rotation().xx() - (*ideal)->rotation().xx();
       auto xy_diff = misaligned->rotation().xy() - (*ideal)->rotation().xy();
       auto xz_diff = misaligned->rotation().xz() - (*ideal)->rotation().xz();
       auto yx_diff = misaligned->rotation().yx() - (*ideal)->rotation().yx();
       auto yy_diff = misaligned->rotation().yy() - (*ideal)->rotation().yy();
       auto yz_diff = misaligned->rotation().yz() - (*ideal)->rotation().yz();
       auto zx_diff = misaligned->rotation().zx() - (*ideal)->rotation().zx();
       auto zy_diff = misaligned->rotation().zy() - (*ideal)->rotation().zy();
       auto zz_diff = misaligned->rotation().zz() - (*ideal)->rotation().zz();
 
       if (outlierPullToIdealCut_ > 0.0 &&
           (x_diff * x_diff + y_diff * y_diff + z_diff * z_diff) > outlierPullToIdealCut_ * outlierPullToIdealCut_) {
         ++outlierCounter;
         edm::LogInfo("Alignment") << outlierCounter << ") Outlier found in subdetector " << subDetId
                                   << ":  delta x: " << x_diff << ",  delta y: " << y_diff << ",  delta z: " << z_diff
                                   << ",  delta xx: " << xx_diff << ",  delta xy: " << xy_diff
                                   << ",  delta xz: " << xz_diff << ",  delta yx: " << yx_diff
                                   << ",  delta yx: " << yy_diff << ",  delta yy: " << yz_diff
                                   << ",  delta zz: " << zx_diff << ",  delta zy: " << zy_diff
                                   << ",  delta zz: " << zz_diff << "\n";
         scaleFactor = 0.0;
       }
 
       const AlignTransform::Translation rescaledTranslation{(*ideal)->position().x() + scaleFactor * x_diff,
                                                             (*ideal)->position().y() + scaleFactor * y_diff,
                                                             (*ideal)->position().z() + scaleFactor * z_diff};
 
       const AlignTransform::Rotation rescaledRotation{
           CLHEP::HepRep3x3{(*ideal)->rotation().xx() + scaleFactor * xx_diff,
                            (*ideal)->rotation().xy() + scaleFactor * xy_diff,
                            (*ideal)->rotation().xz() + scaleFactor * xz_diff,
                            (*ideal)->rotation().yx() + scaleFactor * yx_diff,
                            (*ideal)->rotation().yy() + scaleFactor * yy_diff,
                            (*ideal)->rotation().yz() + scaleFactor * yz_diff,
                            (*ideal)->rotation().zx() + scaleFactor * zx_diff,
                            (*ideal)->rotation().zy() + scaleFactor * zy_diff,
                            (*ideal)->rotation().zz() + scaleFactor * zz_diff}};
 
       const AlignTransform rescaledTransform{rescaledTranslation, rescaledRotation, misaligned->rawId()};
       rescaledAlignments.m_align.emplace_back(rescaledTransform);
     }
   }
 
   edm::Service<cond::service::PoolDBOutputService> poolDb;
   if (!poolDb.isAvailable()) {
     throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
   }
   edm::LogInfo("Alignment") << "Writing rescaled tracker-alignment record.";
   const auto& since = cond::timeTypeSpecs[cond::runnumber].beginValue;
   poolDb->writeOneIOV(rescaledAlignments, since, "TrackerAlignmentRcd");
 }

◆ decodeSubDetectors()

MCMisalignmentScaler::ScalerMap MCMisalignmentScaler::decodeSubDetectors ( const edm::VParameterSet & psets )

private

Definition at line 232 of file MCMisalignmentScaler.cc.

References funct::abs(), Exception, DQMScaleToClient_cfi::factor, logInfo, Skims_PA_cff::name, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, muonDTDigis_cfi::pset, AlCaHLTBitMon_QueryRunRegistry::string, GeomDetEnumerators::subDetId, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, and StripSubdetector::TOB.

Referenced by MCMisalignmentScaler().

                                                                                                     {
   // initialize scaler map
   ScalerMap subDetMap;
   for (unsigned int subDetId = 1; subDetId <= 6; ++subDetId) {
     subDetMap[subDetId][1] = 1.0;
     subDetMap[subDetId][2] = 1.0;
   }
 
   // apply scale factors from configuration
   for (const auto& pset : psets) {
     const auto& name = pset.getUntrackedParameter<std::string>("subDetector");
     const auto& factor = pset.getUntrackedParameter<double>("factor");
 
     std::vector<int> sides;
     if (name.find('-') != std::string::npos)
       sides.push_back(1);
     if (name.find('+') != std::string::npos)
       sides.push_back(2);
     if (sides.empty()) {  // -> use both sides
       sides.push_back(1);
       sides.push_back(2);
     }
 
     if (name.find("Tracker") != std::string::npos) {
       for (unsigned int subDetId = 1; subDetId <= 6; ++subDetId) {
         for (const auto& side : sides)
           subDetMap[subDetId][side] *= factor;
       }
       if (sides.size() == 1) {
         // if only one side to be scaled
         // -> scale also the other side for PXB (subdetid = 1)
         subDetMap[PixelSubdetector::PixelBarrel][std::abs(sides[0] - 2) + 1] *= factor;
       }
     } else if (name.find("PXB") != std::string::npos) {
       // ignore sides for PXB
       subDetMap[PixelSubdetector::PixelBarrel][1] *= factor;
       subDetMap[PixelSubdetector::PixelBarrel][2] *= factor;
     } else if (name.find("PXF") != std::string::npos) {
       for (const auto& side : sides)
         subDetMap[PixelSubdetector::PixelEndcap][side] *= factor;
     } else if (name.find("TIB") != std::string::npos) {
       for (const auto& side : sides)
         subDetMap[StripSubdetector::TIB][side] *= factor;
     } else if (name.find("TOB") != std::string::npos) {
       for (const auto& side : sides)
         subDetMap[StripSubdetector::TOB][side] *= factor;
     } else if (name.find("TID") != std::string::npos) {
       for (const auto& side : sides)
         subDetMap[StripSubdetector::TID][side] *= factor;
     } else if (name.find("TEC") != std::string::npos) {
       for (const auto& side : sides)
         subDetMap[StripSubdetector::TEC][side] *= factor;
     } else {
       throw cms::Exception("BadConfig") << "@SUB=MCMisalignmentScaler::decodeSubDetectors\n"
                                         << "Unknown tracker subdetector: " << name
                                         << "\nSupported options: Tracker, PXB, PXF, TIB, TOB, TID, TEC "
                                         << "(possibly decorated with '+' or '-')";
     }
   }
 
   std::stringstream logInfo;
   logInfo << "MC misalignment scale factors:\n";
   for (const auto& subdet : subDetMap) {
     logInfo << "  Subdet " << subdet.first << "\n";
     for (const auto& side : subdet.second) {
       logInfo << "    side " << side.first << ": " << side.second << "\n";
     }
     logInfo << "\n";
   }
   edm::LogInfo("Alignment") << logInfo.str();
 
   return subDetMap;
 }

◆ fillDescriptions()

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

static

Definition at line 307 of file MCMisalignmentScaler.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::addUntracked(), submitPVResolutionJobs::desc, edm::ParameterSetDescription::setComment(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                       {
   edm::ParameterSetDescription desc;
   desc.setComment(
       "Creates rescaled MC misalignment scenario. "
       "PoolDBOutputService must be set up for 'TrackerAlignmentRcd'.");
   edm::ParameterSetDescription descScaler;
   descScaler.setComment(
       "ParameterSet specifying the tracker part to be scaled "
       "by a given factor.");
   descScaler.addUntracked<std::string>("subDetector", "Tracker");
   descScaler.addUntracked<double>("factor", 1.0);
   desc.addVPSet("scalers", descScaler, std::vector<edm::ParameterSet>(1));
   desc.addUntracked<bool>("pullBadModulesToIdeal", false);
   desc.addUntracked<double>("outlierPullToIdealCut", -1.0);
   descriptions.add("mcMisalignmentScaler", desc);
 }

Member Data Documentation

◆ aliToken_

const edm::ESGetToken<Alignments, TrackerAlignmentRcd> MCMisalignmentScaler::aliToken_

private

Definition at line 76 of file MCMisalignmentScaler.cc.

Referenced by analyze().

◆ firstEvent_

bool MCMisalignmentScaler::firstEvent_ {true}

private

Definition at line 86 of file MCMisalignmentScaler.cc.

Referenced by analyze().

◆ geomDetToken_

const edm::ESGetToken<GeometricDet, IdealGeometryRecord> MCMisalignmentScaler::geomDetToken_

private

Definition at line 72 of file MCMisalignmentScaler.cc.

Referenced by analyze().

◆ outlierPullToIdealCut_

const double MCMisalignmentScaler::outlierPullToIdealCut_

private

Definition at line 85 of file MCMisalignmentScaler.cc.

Referenced by analyze().

◆ pixelQualityToken_

const edm::ESGetToken<SiPixelQuality, SiPixelQualityRcd> MCMisalignmentScaler::pixelQualityToken_

private

Definition at line 70 of file MCMisalignmentScaler.cc.

Referenced by analyze().

◆ ptitpToken_

const edm::ESGetToken<PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd> MCMisalignmentScaler::ptitpToken_

private

Definition at line 74 of file MCMisalignmentScaler.cc.

Referenced by analyze().

◆ ptpToken_

const edm::ESGetToken<PTrackerParameters, PTrackerParametersRcd> MCMisalignmentScaler::ptpToken_

private

Definition at line 73 of file MCMisalignmentScaler.cc.

Referenced by analyze().

◆ pullBadModulesToIdeal_

const bool MCMisalignmentScaler::pullBadModulesToIdeal_

private

Definition at line 84 of file MCMisalignmentScaler.cc.

Referenced by analyze().

◆ scalers_

const ScalerMap MCMisalignmentScaler::scalers_

private

Definition at line 83 of file MCMisalignmentScaler.cc.

Referenced by analyze().

◆ stripQualityToken_

const edm::ESGetToken<SiStripQuality, SiStripQualityRcd> MCMisalignmentScaler::stripQualityToken_

private

Definition at line 71 of file MCMisalignmentScaler.cc.

Referenced by analyze().

◆ topoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> MCMisalignmentScaler::topoToken_

private

Definition at line 75 of file MCMisalignmentScaler.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ ScalerMap

Constructor & Destructor Documentation

◆ MCMisalignmentScaler()

◆ ~MCMisalignmentScaler()

Member Function Documentation

◆ analyze()

◆ decodeSubDetectors()

◆ fillDescriptions()

Member Data Documentation

◆ aliToken_

◆ firstEvent_

◆ geomDetToken_

◆ outlierPullToIdealCut_

◆ pixelQualityToken_

◆ ptitpToken_

◆ ptpToken_

◆ pullBadModulesToIdeal_

◆ scalers_

◆ stripQualityToken_

◆ topoToken_