#include <AlignmentProducerBase.h>

Inheritance diagram for AlignmentProducerBase:

Protected Member Functions
	AlignmentProducerBase (const edm::ParameterSet &)

void	beginLuminosityBlockImpl (const edm::LuminosityBlock &, const edm::EventSetup &)
	begin lumi block More...

void	beginRunImpl (const edm::Run &, const edm::EventSetup &)
	begin run More...

void	endLuminosityBlockImpl (const edm::LuminosityBlock &, const edm::EventSetup &)
	end lumi block More...

void	endRunImpl (const edm::Run &, const edm::EventSetup &)
	end run More...

bool	finish ()

virtual bool	getAliClusterValueMap (const edm::Event &, edm::Handle< AliClusterValueMap > &)=0

virtual bool	getBeamSpot (const edm::Event &, edm::Handle< reco::BeamSpot > &)=0

virtual bool	getTkFittedLasBeamCollection (const edm::Run &, edm::Handle< TkFittedLasBeamCollection > &)=0

virtual bool	getTrajTrackAssociationCollection (const edm::Event &, edm::Handle< TrajTrackAssociationCollection > &)=0

virtual bool	getTsosVectorCollection (const edm::Run &, edm::Handle< TsosVectorCollection > &)=0

void	initAlignmentAlgorithm (const edm::EventSetup &, bool update=false)

int	nEvent () const

bool	processEvent (const edm::Event &, const edm::EventSetup &)
	Process event. More...

void	startProcessing ()
	Start processing of events. More...

void	terminateProcessing (const edm::EventSetup *=0)
	Terminate processing of events. More...

virtual	~AlignmentProducerBase ()(false)

Protected Attributes
const edm::InputTag	beamSpotTag_
	BeamSpot. More...

const edm::InputTag	clusterValueMapTag_
	ValueMap containing associtaion cluster-flag. More...

const bool	doMuon_

const bool	doTracker_

std::shared_ptr< CSCGeometry >	muonCSCGeometry_

std::shared_ptr< DTGeometry >	muonDTGeometry_

const edm::InputTag	tjTkAssociationMapTag_
	Map with tracks/trajectories. More...

const edm::InputTag	tkLasBeamTag_
	LAS beams in edm::Run (ignore if empty) More...

std::shared_ptr< TrackerGeometry >	trackerGeometry_

const bool	useExtras_

Private Member Functions
void	addSurveyInfo (Alignable *)
	Adds survey info to an Alignable. More...

void	applyAlignmentsToDB (const edm::EventSetup &)

void	applyAlignmentsToGeometry ()

template<class G , class Rcd , class ErrRcd >
void	applyDB (G *, const edm::EventSetup &, const AlignTransform &) const

template<class G , class DeformationRcd >
void	applyDB (G *, const edm::EventSetup &) const
	Applies DB constants for SurfaceDeformations. More...

void	applyMisalignment ()
	Applies misalignment scenario to . More...

void	buildParameterStore ()
	Creates the , which manages all Alignables. More...

void	createAlignables (const TrackerTopology *, bool update=false)

void	createAlignmentAlgorithm ()
	Creates the choosen alignment algorithm. More...

void	createCalibrations ()
	Creates the calibrations. More...

void	createGeometries (const edm::EventSetup &, const TrackerTopology *)
	Creates ideal geometry from IdealGeometryRecord. More...

void	createMonitors ()
	Creates the monitors. More...

template<typename T >
bool	hasParameter (const edm::ParameterSet &, const std::string &name)

void	initBeamSpot (const edm::Event &)
	Initializes Beamspot of Alignables . More...

void	readInSurveyRcds (const edm::EventSetup &)
	Reads in survey records. More...

bool	setupChanged (const edm::EventSetup &)
	Checks if one of the EventSetup-Records has changed. More...

void	simpleMisalignment (const align::Alignables &, const std::string &, float, float, bool)
	Applies misalignment scenario to . More...

void	storeAlignmentsToDB ()
	Writes Alignments (i.e. Records) to database-file. More...

void	writeDB (Alignments , const std::string &, AlignmentErrorsExtended , const std::string &, const AlignTransform *, cond::Time_t) const

void	writeDB (AlignmentSurfaceDeformations *, const std::string &, cond::Time_t) const

void	writeForRunRange (cond::Time_t)

Private Attributes
AlignableExtras *	alignableExtras_ { 0 }

AlignableMuon *	alignableMuon_ { 0 }

AlignableTracker *	alignableTracker_ { 0 }

std::unique_ptr< AlignmentAlgorithmBase >	alignmentAlgo_

AlignmentParameterStore *	alignmentParameterStore_ { 0 }

const bool	applyDbAlignment_

edm::Handle< reco::BeamSpot >	beamSpot_

Calibrations	calibrations_

const bool	checkDbAlignmentValidity_

edm::ParameterSet	config_

const bool	doMisalignmentScenario_

const bool	enableAlignableUpdates_

cond::Time_t	firstRun_ {cond::timeTypeSpecs[cond::runnumber].endValue}

std::unique_ptr< const Alignments >	globalPositions_
	GlobalPositions that might be read from DB, nullptr otherwise. More...

bool	isAlgoInitialized_ {false}

bool	isDuringLoop_ {false}

AlignmentMonitors	monitors_

int	nevent_ {0}

bool	runAtPCL_ {false}

const bool	saveApeToDB_

const bool	saveDeformationsToDB_

const bool	saveToDB_

const int	stNFixAlignables_

const double	stRandomRotation_

const double	stRandomShift_

const SurveyErrors *	surveyErrors_ { 0 }

size_t	surveyIndex_ {0}

const Alignments *	surveyValues_ { 0 }

const align::RunRanges	uniqueRunRanges_

const bool	useSurvey_

edm::ESWatcher< CSCAlignmentErrorExtendedRcd >	watchCSCAlErrExtRcd_

edm::ESWatcher< CSCAlignmentRcd >	watchCSCAlRcd_

edm::ESWatcher< CSCSurveyErrorExtendedRcd >	watchCSCSurveyErrExtRcd_

edm::ESWatcher< CSCSurveyRcd >	watchCSCSurveyRcd_

edm::ESWatcher< DTAlignmentErrorExtendedRcd >	watchDTAlErrExtRcd_

edm::ESWatcher< DTAlignmentRcd >	watchDTAlRcd_

edm::ESWatcher< DTSurveyErrorExtendedRcd >	watchDTSurveyErrExtRcd_

edm::ESWatcher< DTSurveyRcd >	watchDTSurveyRcd_

edm::ESWatcher< GlobalPositionRcd >	watchGlobalPositionRcd_

edm::ESWatcher< IdealGeometryRecord >	watchIdealGeometryRcd_

edm::ESWatcher< TrackerSurveyErrorExtendedRcd >	watchTkSurveyErrExtRcd_

edm::ESWatcher< TrackerSurveyRcd >	watchTkSurveyRcd_

edm::ESWatcher< TrackerAlignmentErrorExtendedRcd >	watchTrackerAlErrorExtRcd_

edm::ESWatcher< TrackerAlignmentRcd >	watchTrackerAlRcd_

edm::ESWatcher< TrackerSurfaceDeformationRcd >	watchTrackerSurDeRcd_

Detailed Description

Definition at line 71 of file AlignmentProducerBase.h.

Constructor & Destructor Documentation

AlignmentProducerBase::AlignmentProducerBase ( const edm::ParameterSet & config )

protected

Definition at line 35 of file AlignmentProducerBase.cc.

References ALCARECOPromptCalibProdSiPixelAli0T_cff::algoConfig, applyDbAlignment_, beamSpotTag_, cond::TimeTypeSpecs::beginValue, checkDbAlignmentValidity_, clusterValueMapTag_, config_, createAlignmentAlgorithm(), createCalibrations(), createMonitors(), doMisalignmentScenario_, doMuon_, enableAlignableUpdates_, Exception, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), edm::ParameterSet::getUntrackedParameter(), align::makeUniqueRunRanges(), runAtPCL_, cond::runnumber, saveApeToDB_, saveDeformationsToDB_, saveToDB_, stNFixAlignables_, stRandomRotation_, stRandomShift_, cond::timeTypeSpecs, tjTkAssociationMapTag_, tkLasBeamTag_, uniqueRunRanges_, useExtras_, and useSurvey_.

                                                                           :
   doTracker_{config.getUntrackedParameter<bool>("doTracker")},
   doMuon_{config.getUntrackedParameter<bool>("doMuon")},
   useExtras_{config.getUntrackedParameter<bool>("useExtras")},
   tjTkAssociationMapTag_{config.getParameter<edm::InputTag>("tjTkAssociationMapTag")},
   beamSpotTag_{config.getParameter<edm::InputTag>("beamSpotTag")},
   tkLasBeamTag_{config.getParameter<edm::InputTag>("tkLasBeamTag")},
   clusterValueMapTag_{config.getParameter<edm::InputTag>("hitPrescaleMapTag")},
   uniqueRunRanges_
   {align::makeUniqueRunRanges(config.getParameter<edm::VParameterSet>("RunRangeSelection"),
                               cond::timeTypeSpecs[cond::runnumber].beginValue)},
   config_{config},
   stNFixAlignables_{config.getParameter<int>("nFixAlignables")},
   stRandomShift_{config.getParameter<double>("randomShift")},
   stRandomRotation_{config.getParameter<double>("randomRotation")},
   applyDbAlignment_{config.getUntrackedParameter<bool>("applyDbAlignment")},
   checkDbAlignmentValidity_{config.getUntrackedParameter<bool>("checkDbAlignmentValidity")},
   doMisalignmentScenario_{config.getParameter<bool>("doMisalignmentScenario")},
   saveToDB_{config.getParameter<bool>("saveToDB")},
   saveApeToDB_{config.getParameter<bool>("saveApeToDB")},
   saveDeformationsToDB_{config.getParameter<bool>("saveDeformationsToDB")},
   useSurvey_{config.getParameter<bool>("useSurvey")},
   enableAlignableUpdates_{config.getParameter<bool>("enableAlignableUpdates")}
 {
   edm::LogInfo("Alignment") << "@SUB=AlignmentProducerBase::AlignmentProducerBase";
 
   const auto& algoConfig = config_.getParameterSet("algoConfig");
   if (hasParameter<bool>(config_, "runAtPCL")) {
     // configured in main config?
     runAtPCL_ = config_.getParameter<bool>("runAtPCL");
 
     if (hasParameter<bool>(algoConfig, "runAtPCL") &&
         (runAtPCL_ != algoConfig.getParameter<bool>("runAtPCL"))) {
       throw cms::Exception("BadConfig")
         << "Inconsistent settings for 'runAtPCL' in configuration of the "
         << "alignment producer and the alignment algorithm.";
     }
 
   } else if (hasParameter<bool>(algoConfig, "runAtPCL")) {
     // configured in algo config?
     runAtPCL_ = algoConfig.getParameter<bool>("runAtPCL");
 
   } else {
     // assume 'false' if it was not configured
     runAtPCL_ = false;
   }
 
   createAlignmentAlgorithm();
   createMonitors();
   createCalibrations();
 }

AlignmentProducerBase::~AlignmentProducerBase ( )

protectedvirtual

Definition at line 90 of file AlignmentProducerBase.cc.

References alignableExtras_, alignableMuon_, alignableTracker_, alignmentParameterStore_, and calibrations_.

 {
   for (auto& iCal: calibrations_) delete iCal;
 
   delete alignmentParameterStore_;
   delete alignableExtras_;
   delete alignableTracker_;
   delete alignableMuon_;
 }

Member Function Documentation

void AlignmentProducerBase::addSurveyInfo ( Alignable * ali )

private

Adds survey info to an Alignable.

Definition at line 885 of file AlignmentProducerBase.cc.

References Alignable::alignableObjectId(), AlCaHLTBitMon_QueryRunRegistry::comp, Alignable::components(), relativeConstraints::error, Exception, Alignable::id(), AlignableSurface::length(), Alignments::m_align, SurveyErrors::m_surveyErrors, SurveyError::matrix(), SurveyError::rawId(), makeMuonMisalignmentScenario::rot, AlignableSurface::setLength(), Alignable::setSurvey(), AlignableSurface::setWidth(), SurveyError::structureType(), Alignable::surface(), surveyErrors_, surveyIndex_, surveyValues_, and AlignableSurface::width().

Referenced by readInSurveyRcds().

 {
   const auto& comps = ali->components();
 
   for (const auto& comp: comps) addSurveyInfo(comp);
 
   const SurveyError& error = surveyErrors_->m_surveyErrors[surveyIndex_];
 
   if ( ali->id() != error.rawId() ||
        ali->alignableObjectId() != error.structureType() )
     {
       throw cms::Exception("DatabaseError")
         << "Error reading survey info from DB. Mismatched id!";
     }
 
   const auto& pos = surveyValues_->m_align[surveyIndex_].translation();
   const auto& rot = surveyValues_->m_align[surveyIndex_].rotation();
 
   AlignableSurface surf(align::PositionType(pos.x(), pos.y(), pos.z()),
                         align::RotationType(rot.xx(), rot.xy(), rot.xz(),
                                             rot.yx(), rot.yy(), rot.yz(),
                                             rot.zx(), rot.zy(), rot.zz()));
 
   surf.setWidth(ali->surface().width());
   surf.setLength(ali->surface().length());
 
   ali->setSurvey(new SurveyDet(surf, error.matrix()));
 
   ++surveyIndex_;
 }

void AlignmentProducerBase::applyAlignmentsToDB ( const edm::EventSetup & setup )

private

Applies Alignments from Database (GlobalPositionRcd) to Geometry

Definition at line 548 of file AlignmentProducerBase.cc.

References applyDB(), applyDbAlignment_, align::DetectorGlobalPosition(), doMuon_, doTracker_, edm::EventSetup::get(), globalPositions_, DetId::Muon, muonCSCGeometry_, muonDTGeometry_, GeneralSetup::setup(), DetId::Tracker, and trackerGeometry_.

Referenced by initAlignmentAlgorithm().

 {
   // Retrieve and apply alignments, if requested (requires z setup)
   if (applyDbAlignment_) {
     // we need GlobalPositionRcd - and have to keep track for later removal
     // before writing again to DB...
 
     edm::ESHandle<Alignments> globalAlignments;
     setup.get<GlobalPositionRcd>().get(globalAlignments);
     globalPositions_ = std::make_unique<Alignments>(*globalAlignments);
 
     if (doTracker_) {
       applyDB<TrackerGeometry,
               TrackerAlignmentRcd,
               TrackerAlignmentErrorExtendedRcd>
         (trackerGeometry_.get(), setup,
          align::DetectorGlobalPosition(*globalPositions_, DetId(DetId::Tracker))
          );
 
       applyDB<TrackerGeometry,
               TrackerSurfaceDeformationRcd>(trackerGeometry_.get(), setup);
     }
 
     if (doMuon_) {
       applyDB<DTGeometry,
               DTAlignmentRcd,
               DTAlignmentErrorExtendedRcd>
         (muonDTGeometry_.get(), setup,
          align::DetectorGlobalPosition(*globalPositions_, DetId(DetId::Muon))
          );
 
       applyDB<CSCGeometry,
               CSCAlignmentRcd,
               CSCAlignmentErrorExtendedRcd>
         (muonCSCGeometry_.get(), setup,
          align::DetectorGlobalPosition(*globalPositions_, DetId(DetId::Muon))
          );
     }
   }
 }

void AlignmentProducerBase::applyAlignmentsToGeometry ( )

private

Applies Alignments, AlignmentErrors and SurfaceDeformations to

Definition at line 771 of file AlignmentProducerBase.cc.

References alignableMuon_, alignableTracker_, AlignableTracker::alignmentErrors(), AlignableTracker::alignments(), GeometryAligner::applyAlignments(), GeometryAligner::attachSurfaceDeformations(), AlignableMuon::cscAlignmentErrorsExtended(), AlignableMuon::cscAlignments(), doMuon_, doTracker_, AlignableMuon::dtAlignmentErrorsExtended(), AlignableMuon::dtAlignments(), Exception, muonCSCGeometry_, muonDTGeometry_, Alignable::surfaceDeformations(), and trackerGeometry_.

Referenced by initAlignmentAlgorithm(), and startProcessing().

 {
   edm::LogInfo("Alignment")
     << "@SUB=AlignmentProducerBase::applyAlignmentsToGeometry"
     << "Now physically apply alignments to  geometry...";
 
   // Propagate changes to reconstruction geometry (from initialisation or iteration)
   GeometryAligner aligner;
 
   if (doTracker_) {
     if (!alignableTracker_) {
       throw cms::Exception("LogicError")
         << "@SUB=AlignmentProducerBase::applyAlignmentsToGeometry\n"
         << "Trying to apply tracker alignment before creating it.";
     }
 
     std::unique_ptr<Alignments> alignments{alignableTracker_->alignments()};
     std::unique_ptr<AlignmentErrorsExtended>
       alignmentErrExt{alignableTracker_->alignmentErrors()};
     std::unique_ptr<AlignmentSurfaceDeformations>
       aliDeforms{alignableTracker_->surfaceDeformations()};
 
     aligner.applyAlignments(trackerGeometry_.get(), alignments.get(),
                             alignmentErrExt.get(), AlignTransform());
     aligner.attachSurfaceDeformations(trackerGeometry_.get(), aliDeforms.get());
   }
 
   if (doMuon_) {
     if (!alignableMuon_) {
       throw cms::Exception("LogicError")
         << "@SUB=AlignmentProducerBase::applyAlignmentsToGeometry\n"
         << "Trying to apply muon alignment before creating it.";
     }
 
     std::unique_ptr<Alignments> dtAlignments{alignableMuon_->dtAlignments()};
     std::unique_ptr<Alignments> cscAlignments{alignableMuon_->cscAlignments()};
 
     std::unique_ptr<AlignmentErrorsExtended>
       dtAlignmentErrExt{alignableMuon_->dtAlignmentErrorsExtended()};
     std::unique_ptr<AlignmentErrorsExtended>
       cscAlignmentErrExt{alignableMuon_->cscAlignmentErrorsExtended()};
 
     aligner.applyAlignments(muonDTGeometry_.get(), dtAlignments.get(),
                             dtAlignmentErrExt.get(), AlignTransform());
     aligner.applyAlignments(muonCSCGeometry_.get(), cscAlignments.get(),
                             cscAlignmentErrExt.get(), AlignTransform());
   }
 }

template<class G , class Rcd , class ErrRcd >

void AlignmentProducerBase::applyDB	(	G *	geometry,
		const edm::EventSetup &	iSetup,
		const AlignTransform &	globalCoordinates
	)		const

private

Applies DB constants belonging to (Err)Rcd to Geometry, taking into account 'globalPosition' correction.

Definition at line 296 of file AlignmentProducerBase.h.

References GeometryAligner::applyAlignments(), edm::IOVSyncValue::beginOfTime(), checkDbAlignmentValidity_, edm::IOVSyncValue::endOfTime(), edm::IOVSyncValue::eventID(), Exception, edm::ValidityInterval::first(), plotBeamSpotDB::first, callgraph::G, geometry, edm::EventSetup::get(), edm::ValidityInterval::last(), plotBeamSpotDB::last, record, and edm::EventID::run().

Referenced by applyAlignmentsToDB().

 {
   // 'G' is the geometry class for that DB should be applied,
   // 'Rcd' is the record class for its Alignments
   // 'ErrRcd' is the record class for its AlignmentErrorsExtended
   // 'globalCoordinates' are global transformation for this geometry
 
   const Rcd & record = iSetup.get<Rcd>();
   if (checkDbAlignmentValidity_) {
     const edm::ValidityInterval & validity = record.validityInterval();
     const edm::IOVSyncValue first = validity.first();
     const edm::IOVSyncValue last = validity.last();
     if (first!=edm::IOVSyncValue::beginOfTime() ||
         last!=edm::IOVSyncValue::endOfTime()) {
       throw cms::Exception("DatabaseError")
         << "@SUB=AlignmentProducerBase::applyDB"
         << "\nTrying to apply "
         << record.key().name()
         << " with multiple IOVs in tag.\n"
         << "Validity range is "
         << first.eventID().run() << " - " << last.eventID().run();
     }
   }
 
   edm::ESHandle<Alignments> alignments;
   record.get(alignments);
 
   edm::ESHandle<AlignmentErrorsExtended> alignmentErrors;
   iSetup.get<ErrRcd>().get(alignmentErrors);
 
   GeometryAligner aligner;
   aligner.applyAlignments<G>(geometry, &(*alignments), &(*alignmentErrors),
                              globalCoordinates);
 }

template<class G , class DeformationRcd >

void AlignmentProducerBase::applyDB	(	G *	geometry,
		const edm::EventSetup &	iSetup
	)		const

private

Applies DB constants for SurfaceDeformations.

Definition at line 335 of file AlignmentProducerBase.h.

References GeometryAligner::attachSurfaceDeformations(), edm::IOVSyncValue::beginOfTime(), checkDbAlignmentValidity_, edm::IOVSyncValue::endOfTime(), edm::IOVSyncValue::eventID(), Exception, edm::ValidityInterval::first(), plotBeamSpotDB::first, callgraph::G, geometry, edm::EventSetup::get(), edm::ValidityInterval::last(), plotBeamSpotDB::last, record, and edm::EventID::run().

 {
   // 'G' is the geometry class for that DB should be applied,
   // 'DeformationRcd' is the record class for its surface deformations
 
   const DeformationRcd & record = iSetup.get<DeformationRcd>();
   if (checkDbAlignmentValidity_) {
     const edm::ValidityInterval & validity = record.validityInterval();
     const edm::IOVSyncValue first = validity.first();
     const edm::IOVSyncValue last = validity.last();
     if (first!=edm::IOVSyncValue::beginOfTime() ||
         last!=edm::IOVSyncValue::endOfTime()) {
       throw cms::Exception("DatabaseError")
         << "@SUB=AlignmentProducerBase::applyDB"
         << "\nTrying to apply "
         << record.key().name()
         << " with multiple IOVs in tag.\n"
         << "Validity range is "
         << first.eventID().run() << " - " << last.eventID().run();
     }
   }
   edm::ESHandle<AlignmentSurfaceDeformations> surfaceDeformations;
   record.get(surfaceDeformations);
 
   GeometryAligner aligner;
   aligner.attachSurfaceDeformations<G>(geometry, &(*surfaceDeformations));
 }

void AlignmentProducerBase::applyMisalignment ( )

private

Applies misalignment scenario to .

Definition at line 661 of file AlignmentProducerBase.cc.

References alignableMuon_, AlignmentParameterStore::alignables(), alignableTracker_, alignmentParameterStore_, MuonScenarioBuilder::applyScenario(), TrackerScenarioBuilder::applyScenario(), config_, doMisalignmentScenario_, doMuon_, doTracker_, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), simpleMisalignment(), stRandomRotation_, stRandomShift_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by initAlignmentAlgorithm().

 {
   // Apply misalignment scenario to alignable tracker and muon if requested
   // WARNING: this assumes scenarioConfig can be passed to both muon and tracker
 
   if (doMisalignmentScenario_ && (doTracker_ || doMuon_)) {
     edm::LogInfo("Alignment")
       << "@SUB=AlignmentProducerBase::applyMisalignment"
       << "Applying misalignment scenario to "
       << (doTracker_ ? "tracker" : "")
       << (doMuon_    ? (doTracker_ ? " and muon" : "muon") : ".");
 
     const auto& scenarioConfig = config_.getParameterSet("MisalignmentScenario");
 
     if (doTracker_) {
       TrackerScenarioBuilder scenarioBuilder(alignableTracker_);
       scenarioBuilder.applyScenario(scenarioConfig);
     }
     if (doMuon_) {
       MuonScenarioBuilder muonScenarioBuilder(alignableMuon_);
       muonScenarioBuilder.applyScenario(scenarioConfig);
     }
 
   } else {
     edm::LogInfo("Alignment")
       << "@SUB=AlignmentProducerBase::applyMisalignment"
       << "NOT applying misalignment scenario!";
   }
 
   // Apply simple misalignment
   const auto& sParSel =
     config_.getParameter<std::string>("parameterSelectorSimple");
   simpleMisalignment(alignmentParameterStore_->alignables(), sParSel,
                      stRandomShift_, stRandomRotation_, true);
 }

void AlignmentProducerBase::beginLuminosityBlockImpl	(	const edm::LuminosityBlock &	,
		const edm::EventSetup &	setup
	)

protected

begin lumi block

Definition at line 293 of file AlignmentProducerBase.cc.

References alignmentAlgo_.

Referenced by AlignmentProducerAsAnalyzer::beginLuminosityBlock(), and AlignmentProducer::beginLuminosityBlock().

 {
   // Do not forward edm::LuminosityBlock
   alignmentAlgo_->beginLuminosityBlock(setup);
 }

void AlignmentProducerBase::beginRunImpl	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

protected

begin run

Definition at line 243 of file AlignmentProducerBase.cc.

References alignmentAlgo_, enableAlignableUpdates_, firstRun_, edm::RunBase::id(), initAlignmentAlgorithm(), edm::RunID::run(), runAtPCL_, and setupChanged().

Referenced by AlignmentProducerAsAnalyzer::beginRun(), and AlignmentProducer::beginRun().

 {
   const bool changed{setupChanged(setup)};
   if (changed) {
     edm::LogInfo("Alignment")
       << "@SUB=AlignmentProducerBase::beginRunImpl"
       << "EventSetup-Record changed.";
 
     // updatable alignables are currently not used at PCL, but event setup
     // changes require a complete re-initialization
     if (runAtPCL_) {
       initAlignmentAlgorithm(setup, /* update = */ false);
     } else if (enableAlignableUpdates_) {
       initAlignmentAlgorithm(setup, /* update = */ true);
     }
   }
 
   alignmentAlgo_->beginRun(run, setup,
                            changed && (runAtPCL_ || enableAlignableUpdates_));
 
   //store the first run analyzed to be used for setting the IOV (for PCL)
   if (firstRun_ > static_cast<cond::Time_t>(run.id().run())) {
     firstRun_ = static_cast<cond::Time_t>(run.id().run());
   }
 }

void AlignmentProducerBase::buildParameterStore ( )

private

Creates the , which manages all Alignables.

Definition at line 622 of file AlignmentProducerBase.cc.

References alignableExtras_, alignableMuon_, alignableTracker_, AlignmentParameterStore_cfi::AlignmentParameterStore, alignmentParameterStore_, config_, edm::ParameterSet::getParameter(), and stNFixAlignables_.

Referenced by initAlignmentAlgorithm().

 {
   // Create alignment parameter builder
   edm::LogInfo("Alignment")
     << "@SUB=AlignmentProducerBase::buildParameterStore"
     << "Creating AlignmentParameterBuilder";
 
   const auto& alParamBuildCfg =
     config_.getParameter<edm::ParameterSet>("ParameterBuilder");
   const auto& alParamStoreCfg =
     config_.getParameter<edm::ParameterSet>("ParameterStore");
 
   AlignmentParameterBuilder alignmentParameterBuilder{alignableTracker_,
                                                       alignableMuon_,
                                                       alignableExtras_,
                                                       alParamBuildCfg};
 
   // Fix alignables if requested
   if (stNFixAlignables_ > 0) {
     alignmentParameterBuilder.fixAlignables(stNFixAlignables_);
   }
 
   // Get list of alignables
   const auto& alignables = alignmentParameterBuilder.alignables();
   edm::LogInfo("Alignment")
     << "@SUB=AlignmentProducerBase::buildParameterStore"
     << "got " << alignables.size() << " alignables";
 
   // Create AlignmentParameterStore
   alignmentParameterStore_ =
     new AlignmentParameterStore(alignables, alParamStoreCfg);
   edm::LogInfo("Alignment")
     << "@SUB=AlignmentProducerBase::buildParameterStore"
     << "AlignmentParameterStore created!";
 }

void AlignmentProducerBase::createAlignables	(	const TrackerTopology *	tTopo,
		bool	update = `false`
	)

private

Creates Alignables from the previously loaded Geometry

Definition at line 592 of file AlignmentProducerBase.cc.

References alignableExtras_, align::AlignableMuon, alignableMuon_, alignableTracker_, doMuon_, doTracker_, muonCSCGeometry_, muonDTGeometry_, trackerGeometry_, AlignableTracker::update(), AlignableMuon::update(), and useExtras_.

Referenced by initAlignmentAlgorithm().

 {
   if (doTracker_) {
     if (update) {
       alignableTracker_->update(trackerGeometry_.get(), tTopo);
     } else {
       alignableTracker_ = new AlignableTracker(trackerGeometry_.get(), tTopo);
     }
   }
 
   if (doMuon_) {
     if (update) {
       alignableMuon_->update(muonDTGeometry_.get(), muonCSCGeometry_.get());
     } else {
       alignableMuon_ = new AlignableMuon(muonDTGeometry_.get(), muonCSCGeometry_.get());
     }
   }
 
   if (useExtras_) {
     if (update) {
       // FIXME: Requires further code changes to track beam spot condition changes
     } else {
       alignableExtras_ = new AlignableExtras();
     }
   }
 }

void AlignmentProducerBase::createAlignmentAlgorithm ( )

private

Creates the choosen alignment algorithm.

Definition at line 314 of file AlignmentProducerBase.cc.

References ALCARECOPromptCalibProdSiPixelAli0T_cff::algoConfig, HIPAlignmentAlgorithm_cfi::algoName, alignmentAlgo_, config_, enableAlignableUpdates_, Exception, reco::get(), edm::ParameterSet::getParameter(), runAtPCL_, AlCaHLTBitMon_QueryRunRegistry::string, and uniqueRunRanges_.

Referenced by AlignmentProducerBase().

 {
   auto algoConfig = config_.getParameter<edm::ParameterSet>("algoConfig");
   algoConfig.addUntrackedParameter("RunRangeSelection",
                                    config_.getParameter<edm::VParameterSet>("RunRangeSelection"));
   algoConfig.addUntrackedParameter<align::RunNumber>("firstIOV",
                                                      runAtPCL_ ?
                                                      1 :
                                                      uniqueRunRanges_.front().first);
   algoConfig.addUntrackedParameter("enableAlignableUpdates",
                                    enableAlignableUpdates_);
 
   const auto& algoName = algoConfig.getParameter<std::string>("algoName");
   alignmentAlgo_ = std::unique_ptr<AlignmentAlgorithmBase>
     {AlignmentAlgorithmPluginFactory::get()->create(algoName, algoConfig)};
 
   if (!alignmentAlgo_) {
     throw cms::Exception("BadConfig")
       << "Couldn't find the called alignment algorithm: " << algoName;
   }
 }

void AlignmentProducerBase::createCalibrations ( )

private

Creates the calibrations.

Definition at line 358 of file AlignmentProducerBase.cc.

References alignmentAlgo_, AlignmentProducer_cff::calibrations, calibrations_, config_, beamerCreator::create(), Exception, reco::get(), edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by AlignmentProducerBase().

 {
   const auto& calibrations = config_.getParameter<edm::VParameterSet>("calibrations");
   for (const auto& iCalib: calibrations) {
     calibrations_.push_back(IntegratedCalibrationPluginFactory::get()
                             ->create(iCalib.getParameter<std::string>("calibrationName"),
                                      iCalib));
   }
 
   // Not all algorithms support calibrations - so do not pass empty vector
   // and throw if non-empty and not supported:
   if (!calibrations_.empty()) {
     if (alignmentAlgo_->supportsCalibrations()) {
       alignmentAlgo_->addCalibrations(calibrations_);
 
     } else {
       throw cms::Exception("BadConfig")
         << "@SUB=AlignmentProducerBase::createCalibrations\n"
         << "Configured " << calibrations_.size() << " calibration(s) "
         << "for algorithm not supporting it.";
     }
   }
 }

void AlignmentProducerBase::createGeometries	(	const edm::EventSetup &	iSetup,
		const TrackerTopology *	tTopo
	)

private

Creates ideal geometry from IdealGeometryRecord.

Definition at line 515 of file AlignmentProducerBase.cc.

References TrackerGeomBuilderFromGeometricDet::build(), CSCGeometryBuilderFromDDD::build(), DTGeometryBuilderFromDDD::build(), doMuon_, doTracker_, edm::EventSetup::get(), muonCSCGeometry_, muonDTGeometry_, and trackerGeometry_.

Referenced by initAlignmentAlgorithm().

 {
   if (doTracker_) {
     edm::ESHandle<GeometricDet> geometricDet;
     iSetup.get<IdealGeometryRecord>().get(geometricDet);
 
     edm::ESHandle<PTrackerParameters> ptp;
     iSetup.get<PTrackerParametersRcd>().get(ptp);
 
     TrackerGeomBuilderFromGeometricDet trackerBuilder;
 
     trackerGeometry_ = std::shared_ptr<TrackerGeometry>
       (trackerBuilder.build(&(*geometricDet), *ptp, tTopo));
   }
 
   if (doMuon_) {
     edm::ESTransientHandle<DDCompactView> cpv;
     iSetup.get<IdealGeometryRecord>().get(cpv);
     edm::ESHandle<MuonDDDConstants> mdc;
     iSetup.get<MuonNumberingRecord>().get(mdc);
     DTGeometryBuilderFromDDD DTGeometryBuilder;
     CSCGeometryBuilderFromDDD CSCGeometryBuilder;
     muonDTGeometry_ = std::make_shared<DTGeometry>();
     DTGeometryBuilder.build(muonDTGeometry_, &(*cpv), *mdc);
     muonCSCGeometry_ = std::make_shared<CSCGeometry>();
     CSCGeometryBuilder.build(muonCSCGeometry_, &(*cpv), *mdc );
   }
 }

void AlignmentProducerBase::createMonitors ( )

private

Creates the monitors.

Definition at line 339 of file AlignmentProducerBase.cc.

References config_, Exception, reco::get(), edm::ParameterSet::getParameter(), AlignmentProducer_cff::monitorConfig, gather_cfg::monitors, monitors_, and eostools::move().

Referenced by AlignmentProducerBase().

 {
   const auto& monitorConfig = config_.getParameter<edm::ParameterSet>("monitorConfig");
   auto monitors = monitorConfig.getUntrackedParameter<std::vector<std::string> >("monitors");
   for (const auto& miter: monitors) {
     std::unique_ptr<AlignmentMonitorBase> newMonitor
       {AlignmentMonitorPluginFactory::get()
        ->create(miter, monitorConfig.getUntrackedParameterSet(miter))};
 
     if (!newMonitor) {
       throw cms::Exception("BadConfig") << "Couldn't find monitor named " << miter;
     }
     monitors_.emplace_back(std::move(newMonitor));
   }
 }

void AlignmentProducerBase::endLuminosityBlockImpl	(	const edm::LuminosityBlock &	,
		const edm::EventSetup &	setup
	)

protected

end lumi block

Definition at line 303 of file AlignmentProducerBase.cc.

References alignmentAlgo_.

Referenced by AlignmentProducerAsAnalyzer::endLuminosityBlock(), and AlignmentProducer::endLuminosityBlock().

 {
   // Do not forward edm::LuminosityBlock
   alignmentAlgo_->endLuminosityBlock(setup);
 }

void AlignmentProducerBase::endRunImpl	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

protected

end run

Definition at line 272 of file AlignmentProducerBase.cc.

References alignmentAlgo_, edm::InputTag::encode(), getTkFittedLasBeamCollection(), getTsosVectorCollection(), edm::RunBase::id(), and tkLasBeamTag_.

Referenced by AlignmentProducerAsAnalyzer::endRun(), and AlignmentProducer::endRun().

 {
   if (!tkLasBeamTag_.encode().empty()) {
     edm::Handle<TkFittedLasBeamCollection> lasBeams;
     edm::Handle<TsosVectorCollection> tsoses;
     getTkFittedLasBeamCollection(run, lasBeams);
     getTsosVectorCollection(run, tsoses);
 
     alignmentAlgo_->endRun(EndRunInfo(run.id(), &(*lasBeams), &(*tsoses)), setup);
   } else {
     edm::LogInfo("Alignment")
       << "@SUB=AlignmentProducerBase::endRunImpl"
       << "No Tk LAS beams to forward to algorithm.";
     alignmentAlgo_->endRun(EndRunInfo(run.id(), nullptr, nullptr), setup);
   }
 
 }

bool AlignmentProducerBase::finish ( )

protected

Steers activities after end of job, terminates the AlignmentAlgorithm

Definition at line 919 of file AlignmentProducerBase.cc.

References alignmentAlgo_, calibrations_, monitors_, nevent_, saveApeToDB_, saveDeformationsToDB_, saveToDB_, and storeAlignmentsToDB().

Referenced by progressbar.ProgressBar::__next__(), AlignmentProducerAsAnalyzer::endJob(), AlignmentProducer::endOfJob(), and nEvent().

 {
 
   for (const auto& monitor: monitors_) monitor->endOfJob();
 
   if (alignmentAlgo_->processesEvents() && nevent_ == 0) {
     return false;
   }
 
   for (const auto& iCal:calibrations_) iCal->endOfJob();
 
   if (saveToDB_ || saveApeToDB_ || saveDeformationsToDB_) {
     if (alignmentAlgo_->storeAlignments()) storeAlignmentsToDB();
   } else {
     edm::LogInfo("Alignment")
       << "@SUB=AlignmentProducerBase::finish"
       << "No payload to be stored!";
   }
 
   return true;
 }

virtual bool AlignmentProducerBase::getAliClusterValueMap	(	const edm::Event &	,
		edm::Handle< AliClusterValueMap > &
	)

protectedpure virtual

Implemented in AlignmentProducer, and AlignmentProducerAsAnalyzer.

Referenced by nEvent(), and processEvent().

virtual bool AlignmentProducerBase::getBeamSpot	(	const edm::Event &	,
		edm::Handle< reco::BeamSpot > &
	)

protectedpure virtual

Implemented in AlignmentProducer, and AlignmentProducerAsAnalyzer.

Referenced by initBeamSpot(), and nEvent().

virtual bool AlignmentProducerBase::getTkFittedLasBeamCollection	(	const edm::Run &	,
		edm::Handle< TkFittedLasBeamCollection > &
	)

protectedpure virtual

Implemented in AlignmentProducer, and AlignmentProducerAsAnalyzer.

Referenced by endRunImpl(), and nEvent().

virtual bool AlignmentProducerBase::getTrajTrackAssociationCollection	(	const edm::Event &	,
		edm::Handle< TrajTrackAssociationCollection > &
	)

protectedpure virtual

Implemented in AlignmentProducer, and AlignmentProducerAsAnalyzer.

Referenced by nEvent(), and processEvent().

virtual bool AlignmentProducerBase::getTsosVectorCollection	(	const edm::Run &	,
		edm::Handle< TsosVectorCollection > &
	)

protectedpure virtual

Implemented in AlignmentProducer, and AlignmentProducerAsAnalyzer.

Referenced by endRunImpl(), and nEvent().

template<typename T >

bool AlignmentProducerBase::hasParameter	(	const edm::ParameterSet &	config,
		const std::string &	name
	)

private

Definition at line 1133 of file AlignmentProducerBase.cc.

References edm::Exception::categoryCode(), edm::errors::Configuration, MillePedeFileConverter_cfg::e, edm::ParameterSet::getParameter(), cms::Exception::message(), and dataset::name.

 {
   try {
     config.getParameter<T>(name);
   } catch (const edm::Exception& e) {
     if (e.categoryCode() == edm::errors::Configuration) {
       if (e.message().find("MissingParameter") != std::string::npos) {
         return false;
       } else {
         throw;
       }
     } else {
       throw;
     }
   }
 
   return true;
 }

void AlignmentProducerBase::initAlignmentAlgorithm	(	const edm::EventSetup &	setup,
		bool	update = `false`
	)

protected

Creates Geometry and Alignables of the Tracker and initializes the AlignmentAlgorithm

Definition at line 440 of file AlignmentProducerBase.cc.

References alignableExtras_, alignableMuon_, alignableTracker_, alignmentAlgo_, alignmentParameterStore_, applyAlignmentsToDB(), applyAlignmentsToGeometry(), applyMisalignment(), buildParameterStore(), calibrations_, createAlignables(), createGeometries(), edm::EventSetup::get(), isAlgoInitialized_, monitors_, edm::ESHandle< T >::product(), and startProcessing().

Referenced by AlignmentProducer::beginOfJob(), beginRunImpl(), nEvent(), and processEvent().

 {
   edm::LogInfo("Alignment")
     << "@SUB=AlignmentProducerBase::initAlignmentAlgorithm"
     << "Bwgin";
 
   auto isTrueUpdate = update && isAlgoInitialized_;
 
   // Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   setup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   // Create the geometries from the ideal geometries
   createGeometries(setup, tTopo);
 
   applyAlignmentsToDB(setup);
   createAlignables(tTopo, isTrueUpdate);
   buildParameterStore();
   applyMisalignment();
 
   // Initialize alignment algorithm and integrated calibration and pass the
   // latter to algorithm
   edm::LogInfo("Alignment")
     << "@SUB=AlignmentProducerBase::initAlignmentAlgorithm"
     << "Initializing alignment algorithm.";
   alignmentAlgo_->initialize(setup,
                              alignableTracker_,
                              alignableMuon_,
                              alignableExtras_,
                              alignmentParameterStore_);
   isAlgoInitialized_ = true;
 
   applyAlignmentsToGeometry();
 
   if (!isTrueUpdate) {                // only needed the first time
     for (const auto& iCal: calibrations_) {
       iCal->beginOfJob(alignableTracker_, alignableMuon_, alignableExtras_);
     }
     for (const auto& monitor: monitors_) {
      monitor->beginOfJob(alignableTracker_, alignableMuon_, alignmentParameterStore_);
     }
   }
   startProcessing();            // needed if derived class is non-EDLooper-based
                                 // has no effect, if called during loop
 
   edm::LogInfo("Alignment")
     << "@SUB=AlignmentProducerBase::initAlignmentAlgorithm"
     << "End";
 }

void AlignmentProducerBase::initBeamSpot ( const edm::Event & event )

private

Initializes Beamspot of Alignables .

Definition at line 495 of file AlignmentProducerBase.cc.

References alignableExtras_, beamSpot_, reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), getBeamSpot(), AlignableExtras::initializeBeamSpot(), nevent_, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

Referenced by processEvent().

 {
   getBeamSpot(event, beamSpot_);
 
   if (nevent_== 0 && alignableExtras_) {
     edm::LogInfo("Alignment")
       << "@SUB=AlignmentProducerBase::initBeamSpot"
       << "Initializing AlignableBeamSpot";
 
     alignableExtras_->initializeBeamSpot(beamSpot_->x0(),
                                          beamSpot_->y0(),
                                          beamSpot_->z0(),
                                          beamSpot_->dxdz(),
                                          beamSpot_->dydz());
   }
 }

int AlignmentProducerBase::nEvent ( ) const

inlineprotected

Definition at line 101 of file AlignmentProducerBase.h.

References finish(), getAliClusterValueMap(), getBeamSpot(), getTkFittedLasBeamCollection(), getTrajTrackAssociationCollection(), getTsosVectorCollection(), initAlignmentAlgorithm(), nevent_, and update.

Referenced by AlignmentProducer::endOfLoop().

101 { return nevent_; }

AlignmentProducerBase::nevent_

int nevent_

Definition: AlignmentProducerBase.h:237

bool AlignmentProducerBase::processEvent	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

protected

Process event.

Definition at line 156 of file AlignmentProducerBase.cc.

References alignmentAlgo_, beamSpot_, edm::AssociationMap< Tag >::begin(), clusterValueMapTag_, enableAlignableUpdates_, edm::InputTag::encode(), edm::AssociationMap< Tag >::end(), dumpTauVariables_cfi::eventInfo, getAliClusterValueMap(), getTrajTrackAssociationCollection(), mps_fire::i, initAlignmentAlgorithm(), initBeamSpot(), monitors_, nevent_, readInSurveyRcds(), runAtPCL_, and setupChanged().

Referenced by AlignmentProducerAsAnalyzer::analyze(), and AlignmentProducer::duringLoop().

 {
   if (setupChanged(setup)) {
     edm::LogInfo("Alignment")
       << "@SUB=AlignmentProducerBase::processEvent"
       << "EventSetup-Record changed.";
 
     // updatable alignables are currently not used at PCL, but event setup
     // changes require a complete re-initialization
     if (runAtPCL_) {
       initAlignmentAlgorithm(setup, /* update = */ false);
     } else if (enableAlignableUpdates_) {
       initAlignmentAlgorithm(setup, /* update = */ true);
     }
   }
 
   initBeamSpot(event); // must happen every event and before incrementing 'nevent_'
 
   ++nevent_;  // must happen before the check below;
               // otherwise subsequent checks fail for "EmptySource"
 
   if (!alignmentAlgo_->processesEvents()) {
     edm::LogInfo("Alignment")
       << "@SUB=AlignmentProducerBase::processEvent"
       << "Skipping event. The current configuration of the alignment algorithm "
       << "does not need to process any events.";
     return false;
   }
 
   // reading in survey records
   readInSurveyRcds(setup);
 
   // Printout event number
   for ( int i=10; i<10000000; i*=10 ) {
     if ( nevent_<10*i && (nevent_%i)==0 ) {
       edm::LogInfo("Alignment")
         << "@SUB=AlignmentProducerBase::processEvent"
         << "Events processed: " << nevent_;
     }
   }
 
   // Retrieve trajectories and tracks from the event
   // -> merely skip if collection is empty
   edm::Handle<TrajTrackAssociationCollection> handleTrajTracksCollection;
 
   if (getTrajTrackAssociationCollection(event, handleTrajTracksCollection)) {
 
     // Form pairs of trajectories and tracks
     ConstTrajTrackPairs trajTracks;
     for (auto iter  = handleTrajTracksCollection->begin();
               iter != handleTrajTracksCollection->end();
             ++iter) {
       trajTracks.push_back(ConstTrajTrackPair(&(*(*iter).key), &(*(*iter).val)));
     }
 
     // Run the alignment algorithm with its input
     const AliClusterValueMap* clusterValueMapPtr{nullptr};
     if (!clusterValueMapTag_.encode().empty()) {
       edm::Handle<AliClusterValueMap> clusterValueMap;
       getAliClusterValueMap(event, clusterValueMap);
       clusterValueMapPtr = &(*clusterValueMap);
     }
 
 
     const AlignmentAlgorithmBase::EventInfo eventInfo{event.id(),
                                                       trajTracks,
                                                       *beamSpot_,
                                                       clusterValueMapPtr};
     alignmentAlgo_->run(setup, eventInfo);
 
 
     for (const auto& monitor: monitors_) {
       monitor->duringLoop(event, setup, trajTracks); // forward eventInfo?
     }
   } else {
     edm::LogError("Alignment")
       << "@SUB=AlignmentProducerBase::processEvent"
       << "No track collection found: skipping event";
   }
 
   return true;
 }

void AlignmentProducerBase::readInSurveyRcds ( const edm::EventSetup & iSetup )

private

Reads in survey records.

Definition at line 823 of file AlignmentProducerBase.cc.

References addSurveyInfo(), alignableMuon_, alignableTracker_, Reference_intrackfit_cff::barrel, edm::ESWatcher< T >::check(), AlignableMuon::CSCEndcaps(), doMuon_, doTracker_, AlignableMuon::DTBarrel(), makeMuonMisalignmentScenario::endcap, edm::EventSetup::get(), surveyErrors_, surveyIndex_, surveyValues_, useSurvey_, watchTkSurveyErrExtRcd_, and watchTkSurveyRcd_.

Referenced by processEvent().

 {
 
   // Get Survey Rcds and add Survey Info
   if ( doTracker_ && useSurvey_ ){
     bool tkSurveyBool = watchTkSurveyRcd_.check(iSetup);
     bool tkSurveyErrBool = watchTkSurveyErrExtRcd_.check(iSetup);
     edm::LogInfo("Alignment") << "watcher tksurveyrcd: " << tkSurveyBool;
     edm::LogInfo("Alignment") << "watcher tksurveyerrrcd: " << tkSurveyErrBool;
     if ( tkSurveyBool || tkSurveyErrBool){
 
       edm::LogInfo("Alignment") << "ADDING THE SURVEY INFORMATION";
       edm::ESHandle<Alignments> surveys;
       edm::ESHandle<SurveyErrors> surveyErrors;
 
       iSetup.get<TrackerSurveyRcd>().get(surveys);
       iSetup.get<TrackerSurveyErrorExtendedRcd>().get(surveyErrors);
 
       surveyIndex_  = 0;
       surveyValues_ = &*surveys;
       surveyErrors_ = &*surveyErrors;
       addSurveyInfo(alignableTracker_);
     }
   }
 
   if ( doMuon_ && useSurvey_) {
     bool DTSurveyBool = watchTkSurveyRcd_.check(iSetup);
     bool DTSurveyErrBool = watchTkSurveyErrExtRcd_.check(iSetup);
     bool CSCSurveyBool = watchTkSurveyRcd_.check(iSetup);
     bool CSCSurveyErrBool = watchTkSurveyErrExtRcd_.check(iSetup);
 
     if ( DTSurveyBool || DTSurveyErrBool || CSCSurveyBool || CSCSurveyErrBool ){
       edm::ESHandle<Alignments> dtSurveys;
       edm::ESHandle<SurveyErrors> dtSurveyErrors;
       edm::ESHandle<Alignments> cscSurveys;
       edm::ESHandle<SurveyErrors> cscSurveyErrors;
 
       iSetup.get<DTSurveyRcd>().get(dtSurveys);
       iSetup.get<DTSurveyErrorExtendedRcd>().get(dtSurveyErrors);
       iSetup.get<CSCSurveyRcd>().get(cscSurveys);
       iSetup.get<CSCSurveyErrorExtendedRcd>().get(cscSurveyErrors);
 
       surveyIndex_  = 0;
       surveyValues_ = &*dtSurveys;
       surveyErrors_ = &*dtSurveyErrors;
       const auto& barrels = alignableMuon_->DTBarrel();
       for (const auto& barrel: barrels) addSurveyInfo(barrel);
 
       surveyIndex_  = 0;
       surveyValues_ = &*cscSurveys;
       surveyErrors_ = &*cscSurveyErrors;
       const auto& endcaps = alignableMuon_->CSCEndcaps();
       for (const auto& endcap: endcaps) addSurveyInfo(endcap);
 
     }
   }
 
 }

bool AlignmentProducerBase::setupChanged ( const edm::EventSetup & setup )

private

Checks if one of the EventSetup-Records has changed.

Definition at line 385 of file AlignmentProducerBase.cc.

References edm::ESWatcher< T >::check(), doMuon_, doTracker_, watchCSCAlErrExtRcd_, watchCSCAlRcd_, watchDTAlErrExtRcd_, watchDTAlRcd_, watchGlobalPositionRcd_, watchIdealGeometryRcd_, watchTrackerAlErrorExtRcd_, watchTrackerAlRcd_, and watchTrackerSurDeRcd_.

Referenced by beginRunImpl(), and processEvent().

 {
   bool changed{false};
 
   if (watchIdealGeometryRcd_.check(setup)) {
     changed = true;
   }
 
   if (watchGlobalPositionRcd_.check(setup)) {
     changed = true;
   }
 
   if (doTracker_) {
     if (watchTrackerAlRcd_.check(setup)) {
       changed = true;
     }
 
     if (watchTrackerAlErrorExtRcd_.check(setup)) {
       changed = true;
     }
 
     if (watchTrackerSurDeRcd_.check(setup)) {
       changed = true;
     }
   }
 
   if (doMuon_) {
     if (watchDTAlRcd_.check(setup)) {
       changed = true;
     }
 
     if (watchDTAlErrExtRcd_.check(setup)) {
       changed = true;
     }
 
     if (watchCSCAlRcd_.check(setup)) {
       changed = true;
     }
 
     if (watchCSCAlErrExtRcd_.check(setup)) {
       changed = true;
     }
   }
 
   /* TODO: ExtraAlignables: Which record(s) to check?
    *
    if (useExtras_) {}
   */
 
   return changed;
 }

void AlignmentProducerBase::simpleMisalignment	(	const align::Alignables &	alivec,
		const std::string &	selection,
		float	shift,
		float	rot,
		bool	local
	)

private

Applies misalignment scenario to .

Definition at line 700 of file AlignmentProducerBase.cc.

References funct::abs(), AlignmentParameterSelector::convertParamSel(), RigidBodyAlignmentParameters::dalpha, RigidBodyAlignmentParameters::dbeta, RigidBodyAlignmentParameters::dgamma, RigidBodyAlignmentParameters::dx, RigidBodyAlignmentParameters::dy, RigidBodyAlignmentParameters::dz, Exception, RigidBodyAlignmentParameters::N_PARAM, convertSQLitetoXML_cfg::output, alignCSCRings::r, random, indexGen::s2, and align::toMatrix().

Referenced by applyMisalignment().

 {
 
   std::ostringstream output; // collecting output
 
   if (shift > 0. || rot > 0.) {
     output << "Adding random flat shift of max size " << shift
            << " and adding random flat rotation of max size " << rot <<" to ";
 
     std::vector<bool> commSel(0);
     if (selection != "-1") {
       AlignmentParameterSelector aSelector(nullptr,nullptr); // no alignable needed here...
       const std::vector<char> cSel(aSelector.convertParamSel(selection));
       if (cSel.size() < RigidBodyAlignmentParameters::N_PARAM) {
         throw cms::Exception("BadConfig")
           << "[AlignmentProducerBase::simpleMisalignment_]\n"
           << "Expect selection string '" << selection << "' to be at least of length "
           << RigidBodyAlignmentParameters::N_PARAM << " or to be '-1'.\n"
           << "(Most probably you have to adjust the parameter 'parameterSelectorSimple'.)";
       }
       for (const auto& cIter: cSel) {
         commSel.push_back(cIter == '0' ? false : true);
       }
       output << "parameters defined by (" << selection
              << "), representing (x,y,z,alpha,beta,gamma),";
     } else {
       output << "the active parameters of each alignable,";
     }
     output << " in " << (local ? "local" : "global") << " frame.";
 
     for (const auto& ali: alivec) {
       std::vector<bool> mysel(commSel.empty() ? ali->alignmentParameters()->selector() : commSel);
 
       if (std::abs(shift)>0.00001) {
         double s0 = 0., s1 = 0., s2 = 0.;
         if (mysel[RigidBodyAlignmentParameters::dx]) s0 = shift * double(random()%1000-500)/500.;
         if (mysel[RigidBodyAlignmentParameters::dy]) s1 = shift * double(random()%1000-500)/500.;
         if (mysel[RigidBodyAlignmentParameters::dz]) s2 = shift * double(random()%1000-500)/500.;
 
         if (local) ali->move(ali->surface().toGlobal(align::LocalVector(s0,s1,s2)));
         else       ali->move(align::GlobalVector(s0,s1,s2));
 
         //AlignmentPositionError ape(dx,dy,dz);
         //ali->addAlignmentPositionError(ape);
       }
 
       if (std::abs(rot)>0.00001) {
         align::EulerAngles r(3);
         if (mysel[RigidBodyAlignmentParameters::dalpha]) r(1)=rot*double(random()%1000-500)/500.;
         if (mysel[RigidBodyAlignmentParameters::dbeta])  r(2)=rot*double(random()%1000-500)/500.;
         if (mysel[RigidBodyAlignmentParameters::dgamma]) r(3)=rot*double(random()%1000-500)/500.;
 
         const align::RotationType mrot = align::toMatrix(r);
         if (local) ali->rotateInLocalFrame(mrot);
         else       ali->rotateInGlobalFrame(mrot);
 
         //ali->addAlignmentPositionErrorFromRotation(mrot);
       }
     } // end loop on alignables
   } else {
     output << "No simple misalignment added!";
   }
   edm::LogInfo("Alignment")
     << "@SUB=AlignmentProducerBase::simpleMisalignment" << output.str();
 }

void AlignmentProducerBase::startProcessing ( )

protected

Start processing of events.

Definition at line 103 of file AlignmentProducerBase.cc.

References alignmentAlgo_, applyAlignmentsToGeometry(), calibrations_, Exception, isAlgoInitialized_, isDuringLoop_, monitors_, and nevent_.

Referenced by initAlignmentAlgorithm(), and AlignmentProducer::startingNewLoop().

 {
   if (isDuringLoop_) return;
 
   edm::LogInfo("Alignment")
     << "@SUB=AlignmentProducerBase::startProcessing"
     << "Begin";
 
   if (!isAlgoInitialized_) {
     throw cms::Exception("LogicError")
       << "@SUB=AlignmentProducerBase::startProcessing\n"
       << "Trying to start event processing before initializing the alignment "
       << "algorithm.";
   }
 
   nevent_ = 0;
 
   alignmentAlgo_->startNewLoop();
 
   // FIXME: Should this be done in algorithm::startNewLoop()??
   for (const auto& iCal: calibrations_) iCal->startNewLoop();
   for (const auto& monitor: monitors_) monitor->startingNewLoop();
 
   applyAlignmentsToGeometry();
   isDuringLoop_ = true;
 }

void AlignmentProducerBase::storeAlignmentsToDB ( )

private

Writes Alignments (i.e. Records) to database-file.

Definition at line 944 of file AlignmentProducerBase.cc.

References alignableExtras_, alignmentAlgo_, AlignableExtras::beamSpot(), config_, Exception, firstRun_, edm::ParameterSet::getParameterSetVector(), align::makeUniqueRunRanges(), runAtPCL_, saveApeToDB_, saveDeformationsToDB_, saveToDB_, uniqueRunRanges_, and writeForRunRange().

Referenced by finish().

 {
   const auto runRangeSelectionVPSet =
     config_.getParameterSetVector("RunRangeSelection");
 
   // handle PCL use case
   const auto& uniqueRunRanges =
     (runAtPCL_ ?
      align::makeUniqueRunRanges(runRangeSelectionVPSet, firstRun_) :
      uniqueRunRanges_);
 
   std::vector<AlgebraicVector> beamSpotParameters;
 
   for (const auto& iRunRange: uniqueRunRanges) {
 
     alignmentAlgo_->setParametersForRunRange(iRunRange);
 
     // Save alignments to database
     if (saveToDB_ || saveApeToDB_ || saveDeformationsToDB_) {
       writeForRunRange(iRunRange.first);
     }
 
     // Deal with extra alignables, e.g. beam spot
     if (alignableExtras_) {
       auto& alis = alignableExtras_->beamSpot();
       if (!alis.empty()) {
         auto beamSpotAliPars =
           dynamic_cast<BeamSpotAlignmentParameters*>(alis[0]->alignmentParameters());
         if (!beamSpotAliPars) {
           throw cms::Exception("LogicError")
             << "@SUB=AlignmentProducerBase::storeAlignmentsToDB\n"
             << "First alignable of alignableExtras_ does not have "
             << "'BeamSpotAlignmentParameters', while it should have.";
         }
 
         beamSpotParameters.push_back(beamSpotAliPars->parameters());
       }
     }
   }
 
   if (alignableExtras_) {
     std::ostringstream bsOutput;
 
     auto itPar = beamSpotParameters.cbegin();
     for (auto iRunRange = uniqueRunRanges.cbegin();
          iRunRange != uniqueRunRanges.cend();
          ++iRunRange, ++itPar) {
       bsOutput << "Run range: " << (*iRunRange).first << " - " << (*iRunRange).second << "\n";
       bsOutput << "  Displacement: x=" << (*itPar)[0] << ", y=" << (*itPar)[1] << "\n";
       bsOutput << "  Slope: dx/dz=" << (*itPar)[2] << ", dy/dz=" << (*itPar)[3] << "\n";
     }
 
     edm::LogInfo("Alignment")
       << "@SUB=AlignmentProducerBase::storeAlignmentsToDB"
       << "Parameters for alignable beamspot:\n" << bsOutput.str();
   }
 }

void AlignmentProducerBase::terminateProcessing ( const edm::EventSetup * setup = 0 )

protected

Terminate processing of events.

Definition at line 133 of file AlignmentProducerBase.cc.

References alignmentAlgo_, calibrations_, isDuringLoop_, and monitors_.

Referenced by AlignmentProducerAsAnalyzer::endJob(), and AlignmentProducer::endOfLoop().

 {
   if (!isDuringLoop_) return;
 
   edm::LogInfo("Alignment")
     << "@SUB=AlignmentProducerBase::terminateProcessing"
     << "Terminating algorithm.";
   if (setup) {
     alignmentAlgo_->terminate(*setup);
   } else {
     alignmentAlgo_->terminate();
   }
 
   // FIXME: Should this be done in algorithm::terminate()??
   for (const auto& iCal:calibrations_) iCal->endOfLoop();
   for (const auto& monitor: monitors_) monitor->endOfLoop();
 
   isDuringLoop_ = false;
 }

void AlignmentProducerBase::writeDB	(	Alignments *	alignments,
		const std::string &	alignRcd,
		AlignmentErrorsExtended *	alignmentErrors,
		const std::string &	errRcd,
		const AlignTransform *	globalCoordinates,
		cond::Time_t	time
	)		const

private

Writes Alignments and/or AlignmentErrors to DB for record names (removes *globalCoordinates before writing if non-null...). Takes over ownership of Alignments and AlignmentErrors.

Definition at line 1051 of file AlignmentProducerBase.cc.

References Exception, edm::Service< T >::isAvailable(), GeometryAligner::removeGlobalTransform(), saveApeToDB_, saveToDB_, ntuplemaker::time, AlignTransform::transform(), and cond::service::PoolDBOutputService::writeOne().

Referenced by writeForRunRange().

 {
   Alignments* tempAlignments = alignments;
   AlignmentErrorsExtended* tempAlignmentErrorsExtended = alignmentErrors;
 
   // Call service
   edm::Service<cond::service::PoolDBOutputService> poolDb;
   if (!poolDb.isAvailable()) { // Die if not available
     delete tempAlignments;      // promised to take over ownership...
     delete tempAlignmentErrorsExtended; // dito
     throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
   }
 
   if (globalCoordinates  // happens only if (applyDbAlignment_ == true)
       && globalCoordinates->transform() != AlignTransform::Transform::Identity) {
 
     tempAlignments = new Alignments();            // temporary storage for
     tempAlignmentErrorsExtended = new AlignmentErrorsExtended();  // final alignments and errors
 
     GeometryAligner aligner;
     aligner.removeGlobalTransform(alignments, alignmentErrors,
                                   *globalCoordinates,
                                   tempAlignments, tempAlignmentErrorsExtended);
 
     delete alignments;       // have to delete original alignments
     delete alignmentErrors;  // same thing for the errors
 
     edm::LogInfo("Alignment")
       << "@SUB=AlignmentProducerBase::writeDB"
       << "globalCoordinates removed from alignments (" << alignRcd
       << ") and errors (" << alignRcd << ").";
   }
 
   if (saveToDB_) {
     edm::LogInfo("Alignment") << "Writing Alignments for run " << time
                               << " to " << alignRcd << ".";
     poolDb->writeOne<Alignments>(tempAlignments, time, alignRcd);
   } else { // poolDb->writeOne(..) takes over 'alignments' ownership,...
     delete tempAlignments; // ...otherwise we have to delete, as promised!
   }
 
   if (saveApeToDB_) {
     edm::LogInfo("Alignment") << "Writing AlignmentErrorsExtended for run " << time
                               << " to " << errRcd << ".";
     poolDb->writeOne<AlignmentErrorsExtended>(tempAlignmentErrorsExtended, time, errRcd);
   } else { // poolDb->writeOne(..) takes over 'alignmentErrors' ownership,...
     delete tempAlignmentErrorsExtended; // ...otherwise we have to delete, as promised!
   }
 }

void AlignmentProducerBase::writeDB	(	AlignmentSurfaceDeformations *	alignmentSurfaceDeformations,
		const std::string &	surfaceDeformationRcd,
		cond::Time_t	time
	)		const

private

Writes SurfaceDeformations (bows & kinks) to DB for given record name Takes over ownership of AlignmentSurfaceDeformations.

Definition at line 1108 of file AlignmentProducerBase.cc.

References Exception, edm::Service< T >::isAvailable(), saveDeformationsToDB_, ntuplemaker::time, and cond::service::PoolDBOutputService::writeOne().

 {
   // Call service
   edm::Service<cond::service::PoolDBOutputService> poolDb;
   if (!poolDb.isAvailable()) { // Die if not available
     delete alignmentSurfaceDeformations; // promised to take over ownership...
     throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
   }
 
   if (saveDeformationsToDB_) {
     edm::LogInfo("Alignment") << "Writing AlignmentSurfaceDeformations for run " << time
                               << " to " << surfaceDeformationRcd  << ".";
     poolDb->writeOne<AlignmentSurfaceDeformations>(alignmentSurfaceDeformations, time,
                                                    surfaceDeformationRcd);
   } else { // poolDb->writeOne(..) takes over 'surfaceDeformation' ownership,...
     delete alignmentSurfaceDeformations; // ...otherwise we have to delete, as promised!
   }
 }

void AlignmentProducerBase::writeForRunRange ( cond::Time_t time )

private

Writes Alignments and AlignmentErrors for all sub detectors and the given run number

Definition at line 1005 of file AlignmentProducerBase.cc.

References alignableMuon_, alignableTracker_, AlignableTracker::alignmentErrors(), AlignableTracker::alignments(), AlignableMuon::cscAlignmentErrorsExtended(), AlignableMuon::cscAlignments(), align::DetectorGlobalPosition(), doMuon_, doTracker_, AlignableMuon::dtAlignmentErrorsExtended(), AlignableMuon::dtAlignments(), globalPositions_, DetId::Muon, saveDeformationsToDB_, Alignable::surfaceDeformations(), DetId::Tracker, and writeDB().

Referenced by storeAlignmentsToDB().

 {
   if (doTracker_) { // first tracker
     const AlignTransform* trackerGlobal{nullptr}; // will be 'removed' from constants
     if (globalPositions_) { // i.e. applied before in applyDB
       trackerGlobal = &align::DetectorGlobalPosition(*globalPositions_,
                                                      DetId(DetId::Tracker));
     }
 
     auto alignments = alignableTracker_->alignments();
     auto alignmentErrors = alignableTracker_->alignmentErrors();
     this->writeDB(alignments, "TrackerAlignmentRcd",
                   alignmentErrors, "TrackerAlignmentErrorExtendedRcd", trackerGlobal,
                   time);
 
     // Save surface deformations to database
     if (saveDeformationsToDB_) {
       auto alignmentSurfaceDeformations = alignableTracker_->surfaceDeformations();
       this->writeDB(alignmentSurfaceDeformations, "TrackerSurfaceDeformationRcd", time);
     }
   }
 
   if (doMuon_) { // now muon
     const AlignTransform* muonGlobal{nullptr}; // will be 'removed' from constants
     if (globalPositions_) { // i.e. applied before in applyDB
       muonGlobal = &align::DetectorGlobalPosition(*globalPositions_,
                                                   DetId(DetId::Muon));
     }
     // Get alignments+errors, first DT - ownership taken over by writeDB(..), so no delete
     auto alignments = alignableMuon_->dtAlignments();
     auto alignmentErrors = alignableMuon_->dtAlignmentErrorsExtended();
     this->writeDB(alignments, "DTAlignmentRcd",
                   alignmentErrors, "DTAlignmentErrorExtendedRcd", muonGlobal,
                   time);
 
     // Get alignments+errors, now CSC - ownership taken over by writeDB(..), so no delete
     alignments = alignableMuon_->cscAlignments();
     alignmentErrors = alignableMuon_->cscAlignmentErrorsExtended();
     this->writeDB(alignments, "CSCAlignmentRcd",
                   alignmentErrors, "CSCAlignmentErrorExtendedRcd", muonGlobal,
                   time);
   }
 }