Inheritance diagram for TrackerSystematicMisalignments:

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

void	beginJob () override
	Read ideal tracker geometry from DB. More...

	TrackerSystematicMisalignments (const edm::ParameterSet &)

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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

	~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

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	applySystematicMisalignment (Alignable *)

align::GlobalVector	findSystematicMis (const align::PositionType &, const bool blindToZ, const bool blindToR)

Private Attributes
const edm::ESGetToken< AlignmentErrorsExtended, TrackerAlignmentErrorExtendedRcd >	aliErrorToken_

const edm::ESGetToken< Alignments, TrackerAlignmentRcd >	aliToken_

const edm::ESGetToken< GeometricDet, IdealGeometryRecord >	geomDetToken_

const edm::ESGetToken< Alignments, GlobalPositionRcd >	gprToken_

double	m_bowingEpsilon

double	m_ellipticalDelta

double	m_ellipticalEpsilon

bool	m_fromDBGeom

double	m_layerRotEpsilon

double	m_radialEpsilon

double	m_sagittaDelta

double	m_sagittaEpsilon

double	m_skewDelta

double	m_skewEpsilon

double	m_telescopeEpsilon

double	m_twistEpsilon

double	m_zExpEpsilon

bool	oldMinusZconvention

const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd >	ptitpToken_

const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd >	ptpToken_

bool	suppressBlindMvmts

AlignableTracker *	theAlignableTracker

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

Static Public Member Functions inherited from edm::one::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 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

Class to misaligned tracker from DB.

Date: 2012/06/13 09:24:50

Revision: 1.5

Author: Chung Khim Lae

Definition at line 47 of file TrackerSystematicMisalignments.cc.

Constructor & Destructor Documentation

◆ TrackerSystematicMisalignments()

TrackerSystematicMisalignments::TrackerSystematicMisalignments ( const edm::ParameterSet & cfg )

Definition at line 100 of file TrackerSystematicMisalignments.cc.

References looper::cfg, m_bowingEpsilon, m_ellipticalDelta, m_ellipticalEpsilon, m_fromDBGeom, m_layerRotEpsilon, m_radialEpsilon, m_sagittaDelta, m_sagittaEpsilon, m_skewDelta, m_skewEpsilon, m_telescopeEpsilon, m_twistEpsilon, m_zExpEpsilon, oldMinusZconvention, and suppressBlindMvmts.

     : geomDetToken_(esConsumes()),
       ptpToken_(esConsumes()),
       ptitpToken_(esConsumes()),
       topoToken_(esConsumes()),
       aliToken_(esConsumes()),
       aliErrorToken_(esConsumes()),
       gprToken_(esConsumes()),
       theAlignableTracker(nullptr) {
   // use existing geometry
   m_fromDBGeom = cfg.getUntrackedParameter<bool>("fromDBGeom");
 
   // constants
   m_radialEpsilon = cfg.getUntrackedParameter<double>("radialEpsilon");
   m_telescopeEpsilon = cfg.getUntrackedParameter<double>("telescopeEpsilon");
   m_layerRotEpsilon = cfg.getUntrackedParameter<double>("layerRotEpsilon");
   m_bowingEpsilon = cfg.getUntrackedParameter<double>("bowingEpsilon");
   m_zExpEpsilon = cfg.getUntrackedParameter<double>("zExpEpsilon");
   m_twistEpsilon = cfg.getUntrackedParameter<double>("twistEpsilon");
   m_ellipticalEpsilon = cfg.getUntrackedParameter<double>("ellipticalEpsilon");
   m_skewEpsilon = cfg.getUntrackedParameter<double>("skewEpsilon");
   m_sagittaEpsilon = cfg.getUntrackedParameter<double>("sagittaEpsilon");
 
   m_ellipticalDelta = cfg.getUntrackedParameter<double>("ellipticalDelta");
   m_skewDelta = cfg.getUntrackedParameter<double>("skewDelta");
   m_sagittaDelta = cfg.getUntrackedParameter<double>("sagittaDelta");
 
   if (m_radialEpsilon > -990.0) {
     edm::LogWarning("MisalignedTracker") << "Applying radial ...";
   }
   if (m_telescopeEpsilon > -990.0) {
     edm::LogWarning("MisalignedTracker") << "Applying telescope ...";
   }
   if (m_layerRotEpsilon > -990.0) {
     edm::LogWarning("MisalignedTracker") << "Applying layer rotation ...";
   }
   if (m_bowingEpsilon > -990.0) {
     edm::LogWarning("MisalignedTracker") << "Applying bowing ...";
   }
   if (m_zExpEpsilon > -990.0) {
     edm::LogWarning("MisalignedTracker") << "Applying z-expansion ...";
   }
   if (m_twistEpsilon > -990.0) {
     edm::LogWarning("MisalignedTracker") << "Applying twist ...";
   }
   if (m_ellipticalEpsilon > -990.0) {
     edm::LogWarning("MisalignedTracker") << "Applying elliptical ...";
   }
   if (m_skewEpsilon > -990.0) {
     edm::LogWarning("MisalignedTracker") << "Applying skew ...";
   }
   if (m_sagittaEpsilon > -990.0) {
     edm::LogWarning("MisalignedTracker") << "Applying sagitta ...";
   }
 
   // get flag for suppression of blind movements
   suppressBlindMvmts = cfg.getUntrackedParameter<bool>("suppressBlindMvmts");
   if (suppressBlindMvmts) {
     edm::LogWarning("MisalignedTracker") << "Blind movements suppressed (TIB/TOB in z, TID/TEC in r)";
   }
 
   // compatibility with old (weird) z convention
   oldMinusZconvention = cfg.getUntrackedParameter<bool>("oldMinusZconvention");
   if (oldMinusZconvention) {
     edm::LogWarning("MisalignedTracker") << "Old z convention: dz --> -dz";
   } else {
     edm::LogWarning("MisalignedTracker") << "New z convention: dz --> dz";
   }
 }

Member Function Documentation

◆ analyze()

void TrackerSystematicMisalignments::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 172 of file TrackerSystematicMisalignments.cc.

References aliErrorToken_, AlignableTracker::alignmentErrors(), AlignableTracker::alignments(), aliToken_, GeometryAligner::applyAlignments(), applySystematicMisalignment(), cond::service::PoolDBOutputService::beginOfTime(), TrackerGeomBuilderFromGeometricDet::build(), align::DetectorGlobalPosition(), relativeConstraints::geom, geomDetToken_, gprToken_, edm::Service< T >::isAvailable(), m_fromDBGeom, ptitpToken_, ptpToken_, singleTopDQM_cfi::setup, AlCaHLTBitMon_QueryRunRegistry::string, theAlignableTracker, topoToken_, DetId::Tracker, PbPb_ZMuSkimMuonDPG_cff::tracker, and cond::service::PoolDBOutputService::writeOneIOV().

                                                                                               {
   //Retrieve tracker topology from geometry
   const GeometricDet* geom = &setup.getData(geomDetToken_);
   const PTrackerParameters& ptp = setup.getData(ptpToken_);
   const PTrackerAdditionalParametersPerDet* ptitp = &setup.getData(ptitpToken_);
   const TrackerTopology* tTopo = &setup.getData(topoToken_);
 
   TrackerGeometry* tracker = TrackerGeomBuilderFromGeometricDet().build(geom, ptitp, ptp, tTopo);
 
   //take geometry from DB or randomly generate geometry
   if (m_fromDBGeom) {
     //build the tracker
     const Alignments* alignments = &setup.getData(aliToken_);
     const AlignmentErrorsExtended* alignmentErrors = &setup.getData(aliErrorToken_);
     const Alignments* globalPositionRcd = &setup.getData(gprToken_);
 
     //apply the latest alignments
     GeometryAligner aligner;
     aligner.applyAlignments<TrackerGeometry>(&(*tracker),
                                              &(*alignments),
                                              &(*alignmentErrors),
                                              align::DetectorGlobalPosition(*globalPositionRcd, DetId(DetId::Tracker)));
   }
 
   theAlignableTracker = new AlignableTracker(&(*tracker), tTopo);
 
   applySystematicMisalignment(&(*theAlignableTracker));
 
   // -------------- writing out to alignment record --------------
   Alignments myAlignments = *(theAlignableTracker->alignments());
   AlignmentErrorsExtended myAlignmentErrorsExtended = *(theAlignableTracker->alignmentErrors());
 
   // Store alignment[Error]s to DB
   edm::Service<cond::service::PoolDBOutputService> poolDbService;
   std::string theAlignRecordName = "TrackerAlignmentRcd";
   std::string theErrorRecordName = "TrackerAlignmentErrorExtendedRcd";
 
   // Call service
   if (!poolDbService.isAvailable())  // Die if not available
     throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
 
   poolDbService->writeOneIOV<Alignments>(myAlignments, poolDbService->beginOfTime(), theAlignRecordName);
   poolDbService->writeOneIOV<AlignmentErrorsExtended>(
       myAlignmentErrorsExtended, poolDbService->beginOfTime(), theErrorRecordName);
 }

◆ applySystematicMisalignment()

void TrackerSystematicMisalignments::applySystematicMisalignment ( Alignable * ali )

private

Definition at line 218 of file TrackerSystematicMisalignments.cc.

References Alignable::alignableObjectId(), AlCaHLTBitMon_QueryRunRegistry::comp, Alignable::components(), findSystematicMis(), Alignable::geomDetId(), Alignable::globalPosition(), mps_fire::i, personalPlayback::level, Alignable::move(), DetId::subdetId(), suppressBlindMvmts, SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, and SiStripDetId::TOB.

Referenced by analyze().

                                                                                {
   const align::Alignables& comp = ali->components();
   unsigned int nComp = comp.size();
   //move then do for lower level object
   //for issue of det vs detunit
   bool usecomps = true;
   if ((ali->alignableObjectId() == 2) && (nComp >= 1))
     usecomps = false;
   for (unsigned int i = 0; i < nComp; ++i) {
     if (usecomps)
       applySystematicMisalignment(comp[i]);
   }
 
   // if suppression of blind mvmts: check if subdet is blind to a certain mode
   bool blindToZ(false), blindToR(false);
   if (suppressBlindMvmts) {
     const int subdetid = ali->geomDetId().subdetId();
     switch (subdetid) {
       // TIB/TON blind to z
       case SiStripDetId::TIB:
       case SiStripDetId::TOB:
         blindToZ = true;
         break;
       // TID/TEC blind to R
       case SiStripDetId::TID:
       case SiStripDetId::TEC:
         blindToR = true;
         break;
       default:
         break;
     }
   }
 
   const int level = ali->alignableObjectId();
   if ((level == 1) || (level == 2)) {
     const align::PositionType gP = ali->globalPosition();
     const align::GlobalVector gVec = findSystematicMis(gP, blindToZ, blindToR);
     ali->move(gVec);
   }
 }

◆ beginJob()

void TrackerSystematicMisalignments::beginJob ( )

overridevirtual

Read ideal tracker geometry from DB.

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 170 of file TrackerSystematicMisalignments.cc.

170 {}

◆ findSystematicMis()

align::GlobalVector TrackerSystematicMisalignments::findSystematicMis	(	const align::PositionType &	globalPos,
		const bool	blindToZ,
		const bool	blindToR
	)

private

Definition at line 259 of file TrackerSystematicMisalignments.cc.

References funct::cos(), l1tTrackerHTMiss_cfi::deltaZ, m_bowingEpsilon, m_ellipticalDelta, m_ellipticalEpsilon, m_layerRotEpsilon, m_radialEpsilon, m_sagittaDelta, m_sagittaEpsilon, m_skewDelta, m_skewEpsilon, m_telescopeEpsilon, m_twistEpsilon, m_zExpEpsilon, oldMinusZconvention, PV3DBase< T, PVType, FrameType >::phi(), funct::sin(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by applySystematicMisalignment().

                                                                                            {
   //align::GlobalVector TrackerSystematicMisalignments::findSystematicMis( align::PositionType globalPos ){
   // calculates shift for the current alignable
   // all corrections are calculated w.r.t. the original geometry
   double deltaX = 0.0;
   double deltaY = 0.0;
   double deltaZ = 0.0;
   const double oldX = globalPos.x();
   const double oldY = globalPos.y();
   const double oldZ = globalPos.z();
   const double oldPhi = globalPos.phi();
   const double oldR = sqrt(globalPos.x() * globalPos.x() + globalPos.y() * globalPos.y());
 
   if (m_radialEpsilon > -990.0 && !blindToR) {
     deltaX += m_radialEpsilon * oldX;
     deltaY += m_radialEpsilon * oldY;
   }
   if (m_telescopeEpsilon > -990.0 && !blindToZ) {
     deltaZ += m_telescopeEpsilon * oldR;
   }
   if (m_layerRotEpsilon > -990.0) {
     // The following number was chosen such that the Layer Rotation systematic
     // misalignment would not cause an overall rotation of the tracker.
     const double Roffset = 57.0;
     const double xP = oldR * cos(oldPhi + m_layerRotEpsilon * (oldR - Roffset));
     const double yP = oldR * sin(oldPhi + m_layerRotEpsilon * (oldR - Roffset));
     deltaX += (xP - oldX);
     deltaY += (yP - oldY);
   }
   if (m_bowingEpsilon > -990.0 && !blindToR) {
     const double trackeredgePlusZ = 271.846;
     const double bowfactor = m_bowingEpsilon * (trackeredgePlusZ * trackeredgePlusZ - oldZ * oldZ);
     deltaX += oldX * bowfactor;
     deltaY += oldY * bowfactor;
   }
   if (m_zExpEpsilon > -990.0 && !blindToZ) {
     deltaZ += oldZ * m_zExpEpsilon;
   }
   if (m_twistEpsilon > -990.0) {
     const double xP = oldR * cos(oldPhi + m_twistEpsilon * oldZ);
     const double yP = oldR * sin(oldPhi + m_twistEpsilon * oldZ);
     deltaX += (xP - oldX);
     deltaY += (yP - oldY);
   }
   if (m_ellipticalEpsilon > -990.0 && !blindToR) {
     deltaX += oldX * m_ellipticalEpsilon * cos(2.0 * oldPhi + m_ellipticalDelta);
     deltaY += oldY * m_ellipticalEpsilon * cos(2.0 * oldPhi + m_ellipticalDelta);
   }
   if (m_skewEpsilon > -990.0 && !blindToZ) {
     deltaZ += m_skewEpsilon * cos(oldPhi + m_skewDelta);
   }
   if (m_sagittaEpsilon > -990.0) {
     // deltaX += oldX/fabs(oldX)*m_sagittaEpsilon; // old one...
     deltaX += oldR * m_sagittaEpsilon * sin(m_sagittaDelta);
     deltaY += oldR * m_sagittaEpsilon * cos(m_sagittaDelta);  //Delta y is cos so that delta=0 reflects the old behavior
   }
 
   // Compatibility with old version <= 1.5
   if (oldMinusZconvention)
     deltaZ = -deltaZ;
 
   align::GlobalVector gV(deltaX, deltaY, deltaZ);
   return gV;
 }