test
CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Types | Public Member Functions | Private Member Functions | Private Attributes
AlignmentProducer Class Reference

#include <AlignmentProducer.h>

Inheritance diagram for AlignmentProducer:
edm::ESProducerLooper edm::ESProducer edm::EventSetupRecordIntervalFinder edm::EDLooper edm::ESProxyFactoryProducer edm::EDLooperBase edm::eventsetup::DataProxyProvider

Public Types

typedef std::vector< Alignable * > Alignables
 
typedef std::pair< const
Trajectory *, const
reco::Track * > 
ConstTrajTrackPair
 
typedef std::vector
< ConstTrajTrackPair
ConstTrajTrackPairCollection
 
typedef
AlignmentAlgorithmBase::RunNumber 
RunNumber
 
typedef
AlignmentAlgorithmBase::RunRange 
RunRange
 
typedef std::vector< RunRangeRunRanges
 
- Public Types inherited from edm::eventsetup::DataProxyProvider
typedef std::vector< std::pair
< DataKey,
edm::propagate_const
< std::shared_ptr< DataProxy > > > > 
KeyedProxies
 
typedef std::vector
< EventSetupRecordKey
Keys
 
typedef std::map
< EventSetupRecordKey,
KeyedProxies
RecordProxies
 
- Public Types inherited from edm::EDLooperBase
enum  Status { kContinue, kStop }
 

Public Member Functions

 AlignmentProducer (const edm::ParameterSet &iConfig)
 Constructor. More...
 
virtual void beginLuminosityBlock (const edm::LuminosityBlock &lumiBlock, const edm::EventSetup &setup)
 Called at lumi block start, calling algorithm's beginLuminosityBlock. More...
 
virtual void beginOfJob (const edm::EventSetup &)
 Called at beginning of job. More...
 
virtual void beginRun (const edm::Run &run, const edm::EventSetup &setup)
 Called at run start and calling algorithms beginRun. More...
 
virtual Status duringLoop (const edm::Event &event, const edm::EventSetup &setup)
 Called at each event. More...
 
virtual void endLuminosityBlock (const edm::LuminosityBlock &lumiBlock, const edm::EventSetup &setup)
 Called at lumi block end, calling algorithm's endLuminosityBlock. More...
 
virtual void endOfJob ()
 Called at end of job. More...
 
virtual Status endOfLoop (const edm::EventSetup &, unsigned int iLoop)
 Called at end of loop. More...
 
virtual void endRun (const edm::Run &run, const edm::EventSetup &setup)
 Called at run end - currently reading TkFittedLasBeam if an InpuTag is given for that. More...
 
virtual boost::shared_ptr
< CSCGeometry
produceCSC (const MuonGeometryRecord &iRecord)
 Produce the muon CSC geometry. More...
 
virtual boost::shared_ptr
< DTGeometry
produceDT (const MuonGeometryRecord &iRecord)
 Produce the muon DT geometry. More...
 
virtual boost::shared_ptr
< TrackerGeometry
produceTracker (const TrackerDigiGeometryRecord &iRecord)
 Produce the tracker geometry. More...
 
virtual void startingNewLoop (unsigned int iLoop)
 Called at beginning of loop. More...
 
 ~AlignmentProducer ()
 Destructor. More...
 
- Public Member Functions inherited from edm::ESProducerLooper
 ESProducerLooper ()
 
virtual std::set
< eventsetup::EventSetupRecordKey
modifyingRecords () const
 
- Public Member Functions inherited from edm::ESProducer
 ESProducer ()
 
virtual ~ESProducer ()
 
- Public Member Functions inherited from edm::ESProxyFactoryProducer
 ESProxyFactoryProducer ()
 
virtual void newInterval (const eventsetup::EventSetupRecordKey &iRecordType, const ValidityInterval &iInterval)
 overrides DataProxyProvider method More...
 
virtual ~ESProxyFactoryProducer ()
 
- Public Member Functions inherited from edm::eventsetup::DataProxyProvider
 DataProxyProvider ()
 
const ComponentDescriptiondescription () const
 
bool isUsingRecord (const EventSetupRecordKey &) const
 
const KeyedProxieskeyedProxies (const EventSetupRecordKey &iRecordKey) const
 
void resetProxies (const EventSetupRecordKey &iRecordType)
 
void resetProxiesIfTransient (const EventSetupRecordKey &iRecordType)
 
void setAppendToDataLabel (const edm::ParameterSet &)
 
void setDescription (const ComponentDescription &iDescription)
 
std::set< EventSetupRecordKeyusingRecords () const
 
virtual ~DataProxyProvider ()
 
- Public Member Functions inherited from edm::EventSetupRecordIntervalFinder
const
eventsetup::ComponentDescription
descriptionForFinder () const
 
 EventSetupRecordIntervalFinder ()
 
std::set
< eventsetup::EventSetupRecordKey
findingForRecords () const
 
const ValidityIntervalfindIntervalFor (const eventsetup::EventSetupRecordKey &, const IOVSyncValue &)
 
void setDescriptionForFinder (const eventsetup::ComponentDescription &iDescription)
 
virtual ~EventSetupRecordIntervalFinder ()
 
- Public Member Functions inherited from edm::EDLooper
 EDLooper ()
 
 EDLooper (EDLooper const &)=delete
 
EDLooperoperator= (EDLooper const &)=delete
 
virtual ~EDLooper ()
 
- Public Member Functions inherited from edm::EDLooperBase
virtual void attachTo (ActivityRegistry &)
 Override this method if you need to monitor the state of the processing. More...
 
virtual void beginOfJob ()
 
void copyInfo (ScheduleInfo const &)
 
void doBeginLuminosityBlock (LuminosityBlockPrincipal &, EventSetup const &, ProcessContext *)
 
void doBeginRun (RunPrincipal &, EventSetup const &, ProcessContext *)
 
Status doDuringLoop (EventPrincipal &eventPrincipal, EventSetup const &es, ProcessingController &, StreamContext *)
 
void doEndLuminosityBlock (LuminosityBlockPrincipal &, EventSetup const &, ProcessContext *)
 
Status doEndOfLoop (EventSetup const &es)
 
void doEndRun (RunPrincipal &, EventSetup const &, ProcessContext *)
 
void doStartingNewLoop ()
 
 EDLooperBase ()
 
 EDLooperBase (EDLooperBase const &)=delete
 
EDLooperBaseoperator= (EDLooperBase const &)=delete
 
void prepareForNextLoop (eventsetup::EventSetupProvider *esp)
 
void setActionTable (ExceptionToActionTable const *actionTable)
 
void setModuleChanger (ModuleChanger *)
 
virtual ~EDLooperBase ()
 

Private Member Functions

void addSurveyInfo_ (Alignable *)
 Add survey info to an alignable. More...
 
template<class G , class Rcd , class ErrRcd >
void applyDB (G *geometry, const edm::EventSetup &iSetup, const AlignTransform &globalPosition) const
 
template<class G , class DeformationRcd >
void applyDB (G *geometry, const edm::EventSetup &iSetup) const
 Apply DB constants for surface deformations. More...
 
void createGeometries_ (const edm::EventSetup &)
 Create tracker and muon geometries. More...
 
RunRanges makeNonOverlappingRunRanges (const edm::VParameterSet &RunRangeSelectionVPSet)
 
void readInSurveyRcds (const edm::EventSetup &)
 read in survey records More...
 
void simpleMisalignment_ (const Alignables &alivec, const std::string &selection, float shift, float rot, bool local)
 Apply random shifts and rotations to selected alignables, according to configuration. More...
 
void writeDB (Alignments *alignments, const std::string &alignRcd, AlignmentErrorsExtended *alignmentErrors, const std::string &errRcd, const AlignTransform *globalCoordinates, cond::Time_t time) const
 
void writeDB (AlignmentSurfaceDeformations *alignmentSurfaceDeformations, const std::string &surfaceDeformationRcd, cond::Time_t time) const
 
void writeForRunRange (cond::Time_t time)
 

Private Attributes

const bool applyDbAlignment_
 
const edm::InputTag beamSpotTag_
 
const bool checkDbAlignmentValidity_
 
const edm::InputTag clusterValueMapTag_
 
const bool doMisalignmentScenario_
 
const bool doMuon_
 
const bool doTracker_
 
const AlignmentsglobalPositions_
 GlobalPositions that might be read from DB, NULL otherwise. More...
 
int nevent_
 
const bool saveApeToDB_
 
const bool saveDeformationsToDB_
 
const bool saveToDB_
 
const int stNFixAlignables_
 
const double stRandomRotation_
 
const double stRandomShift_
 
AlignableExtrastheAlignableExtras
 
AlignableMuontheAlignableMuon
 
AlignableTrackertheAlignableTracker
 
AlignmentAlgorithmBasetheAlignmentAlgo
 
AlignmentParameterStoretheAlignmentParameterStore
 
std::vector
< IntegratedCalibrationBase * > 
theCalibrations
 
const unsigned int theMaxLoops
 
std::vector
< AlignmentMonitorBase * > 
theMonitors
 
boost::shared_ptr< CSCGeometrytheMuonCSC
 
boost::shared_ptr< DTGeometrytheMuonDT
 
edm::ParameterSet theParameterSet
 
const SurveyErrorstheSurveyErrors
 
unsigned int theSurveyIndex
 
const AlignmentstheSurveyValues
 
boost::shared_ptr
< TrackerGeometry
theTracker
 
const edm::InputTag tjTkAssociationMapTag_
 
const edm::InputTag tkLasBeamTag_
 
const bool useExtras_
 
const bool useSurvey_
 
edm::ESWatcher
< CSCSurveyErrorExtendedRcd
watchCSCSurveyErrRcd_
 
edm::ESWatcher< CSCSurveyRcdwatchCSCSurveyRcd_
 
edm::ESWatcher
< DTSurveyErrorExtendedRcd
watchDTSurveyErrRcd_
 
edm::ESWatcher< DTSurveyRcdwatchDTSurveyRcd_
 
edm::ESWatcher
< TrackerSurveyErrorExtendedRcd
watchTkSurveyErrRcd_
 
edm::ESWatcher< TrackerSurveyRcdwatchTkSurveyRcd_
 

Additional Inherited Members

- Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void prevalidate (ConfigurationDescriptions &)
 
- Protected Member Functions inherited from edm::ESProducerLooper
virtual void registerFactoryWithKey (const eventsetup::EventSetupRecordKey &iRecord, std::auto_ptr< eventsetup::ProxyFactoryBase > &iFactory, const std::string &iLabel=std::string())
 
void setIntervalFor (const eventsetup::EventSetupRecordKey &iKey, const IOVSyncValue &iTime, ValidityInterval &oInterval)
 
- Protected Member Functions inherited from edm::ESProducer
template<typename T >
void setWhatProduced (T *iThis, const es::Label &iLabel=es::Label())
 
template<typename T >
void setWhatProduced (T *iThis, const char *iLabel)
 
template<typename T >
void setWhatProduced (T *iThis, const std::string &iLabel)
 
template<typename T , typename TDecorator >
void setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel=es::Label())
 
template<typename T , typename TReturn , typename TRecord >
void setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const es::Label &iLabel=es::Label())
 
template<typename T , typename TReturn , typename TRecord , typename TArg >
void setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel=es::Label())
 
- Protected Member Functions inherited from edm::ESProxyFactoryProducer
template<class TFactory >
void registerFactory (std::auto_ptr< TFactory > iFactory, const std::string &iLabel=std::string())
 
virtual void registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList)
 override DataProxyProvider method More...
 
- Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
void eraseAll (const EventSetupRecordKey &iRecordKey)
 deletes all the Proxies in aStream More...
 
void invalidateProxies (const EventSetupRecordKey &iRecordKey)
 
template<class T >
void usingRecord ()
 
void usingRecordWithKey (const EventSetupRecordKey &)
 
- Protected Member Functions inherited from edm::EventSetupRecordIntervalFinder
template<class T >
void findingRecord ()
 
void findingRecordWithKey (const eventsetup::EventSetupRecordKey &)
 
- Protected Member Functions inherited from edm::EDLooperBase
ModuleChangermoduleChanger ()
 This only returns a non-zero value during the call to endOfLoop. More...
 
ScheduleInfo const * scheduleInfo () const
 This returns a non-zero value after the constructor has been called. More...
 

Detailed Description

Package : Alignment/CommonAlignmentProducer Description : calls alignment algorithms

Author
: Frederic Ronga Revision :
Revision:
1.27

last update:

Date:
2012/06/13 16:23:30

by :

Author:
yana

Definition at line 55 of file AlignmentProducer.h.

Member Typedef Documentation

Definition at line 59 of file AlignmentProducer.h.

Definition at line 60 of file AlignmentProducer.h.

Definition at line 61 of file AlignmentProducer.h.

Definition at line 63 of file AlignmentProducer.h.

Definition at line 64 of file AlignmentProducer.h.

typedef std::vector<RunRange> AlignmentProducer::RunRanges

Definition at line 65 of file AlignmentProducer.h.

Constructor & Destructor Documentation

AlignmentProducer::AlignmentProducer ( const edm::ParameterSet iConfig)

Constructor.

Definition at line 80 of file AlignmentProducer.cc.

References edm::ParameterSet::addUntrackedParameter(), beamerCreator::create(), doMuon_, doTracker_, Exception, reco::get(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mergeVDriftHistosByStation::name, produceCSC(), produceDT(), produceTracker(), edm::ESProducer::setWhatProduced(), AlCaHLTBitMon_QueryRunRegistry::string, theAlignmentAlgo, theCalibrations, and theMonitors.

80  :
84  nevent_(0), theParameterSet(iConfig),
85  theMaxLoops( iConfig.getUntrackedParameter<unsigned int>("maxLoops") ),
86  stNFixAlignables_(iConfig.getParameter<int>("nFixAlignables") ),
87  stRandomShift_(iConfig.getParameter<double>("randomShift")),
88  stRandomRotation_(iConfig.getParameter<double>("randomRotation")),
89  applyDbAlignment_( iConfig.getUntrackedParameter<bool>("applyDbAlignment")),
90  checkDbAlignmentValidity_( iConfig.getUntrackedParameter<bool>("checkDbAlignmentValidity")),
91  doMisalignmentScenario_(iConfig.getParameter<bool>("doMisalignmentScenario")),
92  saveToDB_(iConfig.getParameter<bool>("saveToDB")),
93  saveApeToDB_(iConfig.getParameter<bool>("saveApeToDB")),
94  saveDeformationsToDB_(iConfig.getParameter<bool>("saveDeformationsToDB")),
95  doTracker_( iConfig.getUntrackedParameter<bool>("doTracker") ),
96  doMuon_( iConfig.getUntrackedParameter<bool>("doMuon") ),
97  useExtras_( iConfig.getUntrackedParameter<bool>("useExtras") ),
98  useSurvey_( iConfig.getParameter<bool>("useSurvey") ),
99  tjTkAssociationMapTag_(iConfig.getParameter<edm::InputTag>("tjTkAssociationMapTag")),
100  beamSpotTag_(iConfig.getParameter<edm::InputTag>("beamSpotTag")),
101  tkLasBeamTag_(iConfig.getParameter<edm::InputTag>("tkLasBeamTag")),
102  clusterValueMapTag_(iConfig.getParameter<edm::InputTag>("hitPrescaleMapTag"))
103 {
104  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::AlignmentProducer";
105 
106  // Tell the framework what data is being produced
107  if (doTracker_) {
109  }
110  if (doMuon_) {
113  }
114 
115  // Create the alignment algorithm
116  edm::ParameterSet algoConfig = iConfig.getParameter<edm::ParameterSet>( "algoConfig" );
117  edm::VParameterSet iovSelection = iConfig.getParameter<edm::VParameterSet>( "RunRangeSelection" );
118  algoConfig.addUntrackedParameter<edm::VParameterSet>( "RunRangeSelection", iovSelection );
119  std::string algoName = algoConfig.getParameter<std::string>( "algoName" );
120  theAlignmentAlgo = AlignmentAlgorithmPluginFactory::get( )->create( algoName, algoConfig );
121 
122  // Check if found
123  if ( !theAlignmentAlgo )
124  throw cms::Exception("BadConfig") << "Couldn't find algorithm called " << algoName;
125 
126  // Now create monitors:
127  edm::ParameterSet monitorConfig = iConfig.getParameter<edm::ParameterSet>( "monitorConfig" );
128  std::vector<std::string> monitors = monitorConfig.getUntrackedParameter<std::vector<std::string> >( "monitors" );
129  for (std::vector<std::string>::const_iterator miter = monitors.begin(); miter != monitors.end(); ++miter) {
130  AlignmentMonitorBase* newMonitor = AlignmentMonitorPluginFactory::get()->create(*miter, monitorConfig.getUntrackedParameter<edm::ParameterSet>(*miter));
131 
132  if (!newMonitor) throw cms::Exception("BadConfig") << "Couldn't find monitor named " << *miter;
133 
134  theMonitors.push_back(newMonitor);
135  }
136 
137  // Finally create integrated calibrations:
138  edm::VParameterSet calibrations = iConfig.getParameter<edm::VParameterSet>("calibrations");
139  for (auto iCalib = calibrations.begin(); iCalib != calibrations.end(); ++iCalib) {
140  const std::string name(iCalib->getParameter<std::string>("calibrationName"));
142  // exception comes from line before: if (!theCalibrations.back()) throw cms::Exception(..) << ..;
143  }
144 
145 }
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
const bool doMisalignmentScenario_
AlignmentAlgorithmBase * theAlignmentAlgo
virtual boost::shared_ptr< TrackerGeometry > produceTracker(const TrackerDigiGeometryRecord &iRecord)
Produce the tracker geometry.
std::vector< ParameterSet > VParameterSet
Definition: ParameterSet.h:33
virtual boost::shared_ptr< DTGeometry > produceDT(const MuonGeometryRecord &iRecord)
Produce the muon DT geometry.
std::vector< IntegratedCalibrationBase * > theCalibrations
edm::ParameterSet theParameterSet
const edm::InputTag tkLasBeamTag_
void setWhatProduced(T *iThis, const es::Label &iLabel=es::Label())
Definition: ESProducer.h:115
const bool applyDbAlignment_
const unsigned int theMaxLoops
AlignableExtras * theAlignableExtras
const bool checkDbAlignmentValidity_
const edm::InputTag clusterValueMapTag_
std::vector< AlignmentMonitorBase * > theMonitors
AlignmentParameterStore * theAlignmentParameterStore
const edm::InputTag beamSpotTag_
const edm::InputTag tjTkAssociationMapTag_
AlignableMuon * theAlignableMuon
const bool saveDeformationsToDB_
const double stRandomShift_
void addUntrackedParameter(std::string const &name, T const &value)
Definition: ParameterSet.h:208
virtual boost::shared_ptr< CSCGeometry > produceCSC(const MuonGeometryRecord &iRecord)
Produce the muon CSC geometry.
AlignableTracker * theAlignableTracker
const Alignments * globalPositions_
GlobalPositions that might be read from DB, NULL otherwise.
T get(const Candidate &c)
Definition: component.h:55
const double stRandomRotation_
AlignmentProducer::~AlignmentProducer ( )

Destructor.

Definition at line 150 of file AlignmentProducer.cc.

References globalPositions_, theAlignableExtras, theAlignableMuon, theAlignableTracker, theAlignmentAlgo, theAlignmentParameterStore, and theCalibrations.

151 {
152  delete theAlignmentAlgo;
153 
154  // Delete monitors as well??
155 
156  for (auto iCal = theCalibrations.begin(); iCal != theCalibrations.end(); ++iCal) {
157  delete *iCal; // delete integrated calibration pointed to by (*iCal)
158  }
159 
161  delete theAlignableExtras;
162  delete theAlignableTracker;
163  delete theAlignableMuon;
164 
165  delete globalPositions_;
166 }
AlignmentAlgorithmBase * theAlignmentAlgo
std::vector< IntegratedCalibrationBase * > theCalibrations
AlignableExtras * theAlignableExtras
AlignmentParameterStore * theAlignmentParameterStore
AlignableMuon * theAlignableMuon
AlignableTracker * theAlignableTracker
const Alignments * globalPositions_
GlobalPositions that might be read from DB, NULL otherwise.

Member Function Documentation

void AlignmentProducer::addSurveyInfo_ ( Alignable ali)
private

Add survey info to an alignable.

Definition at line 676 of file AlignmentProducer.cc.

References Alignable::alignableObjectId(), AlCaHLTBitMon_QueryRunRegistry::comp, Alignable::components(), relativeConstraints::error, Exception, i, 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(), theSurveyErrors, theSurveyIndex, theSurveyValues, and AlignableSurface::width().

Referenced by readInSurveyRcds().

677 {
678  const std::vector<Alignable*>& comp = ali->components();
679 
680  unsigned int nComp = comp.size();
681 
682  for (unsigned int i = 0; i < nComp; ++i) addSurveyInfo_(comp[i]);
683 
685 
686  if ( ali->id() != error.rawId() ||
687  ali->alignableObjectId() != error.structureType() )
688  {
689  throw cms::Exception("DatabaseError")
690  << "Error reading survey info from DB. Mismatched id!";
691  }
692 
693  const CLHEP::Hep3Vector& pos = theSurveyValues->m_align[theSurveyIndex].translation();
694  const CLHEP::HepRotation& rot = theSurveyValues->m_align[theSurveyIndex].rotation();
695 
696  AlignableSurface surf( align::PositionType( pos.x(), pos.y(), pos.z() ),
697  align::RotationType( rot.xx(), rot.xy(), rot.xz(),
698  rot.yx(), rot.yy(), rot.yz(),
699  rot.zx(), rot.zy(), rot.zz() ) );
700 
701  surf.setWidth( ali->surface().width() );
702  surf.setLength( ali->surface().length() );
703 
704  ali->setSurvey( new SurveyDet( surf, error.matrix() ) );
705 
706  ++theSurveyIndex;
707 }
align::Scalar width() const
align::ID id() const
Return the ID of Alignable, i.e. DetId of &#39;first&#39; component GeomDet(Unit).
Definition: Alignable.h:185
int i
Definition: DBlmapReader.cc:9
ErrorMatrix matrix() const
Definition: SurveyError.h:76
const Alignments * theSurveyValues
void addSurveyInfo_(Alignable *)
Add survey info to an alignable.
uint8_t structureType() const
Definition: SurveyError.h:66
virtual Alignables components() const =0
Return vector of all direct components.
std::vector< AlignTransform > m_align
Definition: Alignments.h:19
unsigned int theSurveyIndex
virtual StructureType alignableObjectId() const =0
Return the alignable type identifier.
align::ID rawId() const
Definition: SurveyError.h:71
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:131
const SurveyErrors * theSurveyErrors
align::Scalar length() const
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:306
std::vector< SurveyError > m_surveyErrors
Definition: SurveyErrors.h:23
template<class G , class Rcd , class ErrRcd >
void AlignmentProducer::applyDB ( G *  geometry,
const edm::EventSetup iSetup,
const AlignTransform globalPosition 
) const
private

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

Definition at line 775 of file AlignmentProducer.cc.

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().

777 {
778  // 'G' is the geometry class for that DB should be applied,
779  // 'Rcd' is the record class for its Alignments
780  // 'ErrRcd' is the record class for its AlignmentErrorsExtended
781  // 'globalCoordinates' are global transformation for this geometry
782 
783  const Rcd & record = iSetup.get<Rcd>();
785  const edm::ValidityInterval & validity = record.validityInterval();
786  const edm::IOVSyncValue first = validity.first();
787  const edm::IOVSyncValue last = validity.last();
788  if (first!=edm::IOVSyncValue::beginOfTime() ||
790  throw cms::Exception("DatabaseError")
791  << "@SUB=AlignmentProducer::applyDB"
792  << "\nTrying to apply "
793  << record.key().name()
794  << " with multiple IOVs in tag.\n"
795  << "Validity range is "
796  << first.eventID().run() << " - " << last.eventID().run();
797  }
798  }
799 
800  edm::ESHandle<Alignments> alignments;
801  record.get(alignments);
802 
804  iSetup.get<ErrRcd>().get(alignmentErrors);
805 
806  GeometryAligner aligner;
807  aligner.applyAlignments<G>(geometry, &(*alignments), &(*alignmentErrors),
808  globalCoordinates);
809 }
RunNumber_t run() const
Definition: EventID.h:39
const EventID & eventID() const
Definition: IOVSyncValue.h:42
Class to update a given geometry with a set of alignments.
JetCorrectorParameters::Record record
Definition: classes.h:7
static const IOVSyncValue & endOfTime()
Definition: IOVSyncValue.cc:97
const IOVSyncValue & last() const
const bool checkDbAlignmentValidity_
static const IOVSyncValue & beginOfTime()
tuple G
Definition: callgraph.py:12
void applyAlignments(C *geometry, const Alignments *alignments, const AlignmentErrorsExtended *alignmentErrors, const AlignTransform &globalCoordinates)
const T & get() const
Definition: EventSetup.h:56
ESHandle< TrackerGeometry > geometry
const IOVSyncValue & first() const
template<class G , class DeformationRcd >
void AlignmentProducer::applyDB ( G *  geometry,
const edm::EventSetup iSetup 
) const
private

Apply DB constants for surface deformations.

Definition at line 816 of file AlignmentProducer.cc.

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().

817 {
818  // 'G' is the geometry class for that DB should be applied,
819  // 'DeformationRcd' is the record class for its surface deformations
820 
821  const DeformationRcd & record = iSetup.get<DeformationRcd>();
823  const edm::ValidityInterval & validity = record.validityInterval();
824  const edm::IOVSyncValue first = validity.first();
825  const edm::IOVSyncValue last = validity.last();
826  if (first!=edm::IOVSyncValue::beginOfTime() ||
828  throw cms::Exception("DatabaseError")
829  << "@SUB=AlignmentProducer::applyDB"
830  << "\nTrying to apply "
831  << record.key().name()
832  << " with multiple IOVs in tag.\n"
833  << "Validity range is "
834  << first.eventID().run() << " - " << last.eventID().run();
835  }
836  }
838  record.get(surfaceDeformations);
839 
840  GeometryAligner aligner;
841  aligner.attachSurfaceDeformations<G>(geometry, &(*surfaceDeformations));
842 }
RunNumber_t run() const
Definition: EventID.h:39
const EventID & eventID() const
Definition: IOVSyncValue.h:42
Class to update a given geometry with a set of alignments.
JetCorrectorParameters::Record record
Definition: classes.h:7
static const IOVSyncValue & endOfTime()
Definition: IOVSyncValue.cc:97
const IOVSyncValue & last() const
const bool checkDbAlignmentValidity_
static const IOVSyncValue & beginOfTime()
void attachSurfaceDeformations(C *geometry, const AlignmentSurfaceDeformations *surfaceDeformations)
tuple G
Definition: callgraph.py:12
const T & get() const
Definition: EventSetup.h:56
ESHandle< TrackerGeometry > geometry
const IOVSyncValue & first() const
void AlignmentProducer::beginLuminosityBlock ( const edm::LuminosityBlock lumiBlock,
const edm::EventSetup setup 
)
virtual

Called at lumi block start, calling algorithm's beginLuminosityBlock.

Reimplemented from edm::EDLooperBase.

Definition at line 560 of file AlignmentProducer.cc.

References AlignmentAlgorithmBase::beginLuminosityBlock(), and theAlignmentAlgo.

562 {
563  theAlignmentAlgo->beginLuminosityBlock(setup); // do not forward edm::LuminosityBlock
564 }
virtual void beginLuminosityBlock(const edm::EventSetup &setup)
called at begin of luminosity block (no lumi block info passed yet)
AlignmentAlgorithmBase * theAlignmentAlgo
void AlignmentProducer::beginOfJob ( const edm::EventSetup iSetup)
virtual

Called at beginning of job.

Reimplemented from edm::EDLooperBase.

Definition at line 199 of file AlignmentProducer.cc.

References AlignmentAlgorithmBase::addCalibrations(), align::AlignableMuon, AlignmentParameterStore_cfi::AlignmentParameterStore, applyDbAlignment_, MuonScenarioBuilder::applyScenario(), TrackerScenarioBuilder::applyScenario(), createGeometries_(), align::DetectorGlobalPosition(), doMisalignmentScenario_, doMuon_, doTracker_, Exception, edm::EventSetup::get(), edm::ParameterSet::getParameter(), globalPositions_, AlignmentAlgorithmBase::initialize(), DetId::Muon, edm::ESHandle< class >::product(), simpleMisalignment_(), stNFixAlignables_, stRandomRotation_, stRandomShift_, AlCaHLTBitMon_QueryRunRegistry::string, theAlignableExtras, theAlignableMuon, theAlignableTracker, theAlignmentAlgo, theAlignmentParameterStore, theCalibrations, theMonitors, theMuonCSC, theMuonDT, theParameterSet, theTracker, DetId::Tracker, and useExtras_.

200 {
201  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob";
202 
203  //Retrieve tracker topology from geometry
204  edm::ESHandle<TrackerTopology> tTopoHandle;
205  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
206  const TrackerTopology* const tTopo = tTopoHandle.product();
207 
208  // Create the geometries from the ideal geometries (first time only)
209  this->createGeometries_( iSetup );
210 
211  // Retrieve and apply alignments, if requested (requires DB setup)
212  if ( applyDbAlignment_ ) {
213  // we need GlobalPositionRcd - and have to keep track for later removal
214  // before writing again to DB...
215  edm::ESHandle<Alignments> globalPositionRcd;
216  iSetup.get<GlobalPositionRcd>().get(globalPositionRcd);
217  globalPositions_ = new Alignments(*globalPositionRcd);
218 
219  if ( doTracker_ ) { // apply to tracker
220  this->applyDB<TrackerGeometry,TrackerAlignmentRcd,TrackerAlignmentErrorExtendedRcd>
221  (&(*theTracker), iSetup,
223  this->applyDB<TrackerGeometry,TrackerSurfaceDeformationRcd>(&(*theTracker), iSetup);
224  }
225 
226  if ( doMuon_ ) { // apply to tracker
227  this->applyDB<DTGeometry,DTAlignmentRcd,DTAlignmentErrorExtendedRcd>
228  (&(*theMuonDT), iSetup,
230  this->applyDB<CSCGeometry,CSCAlignmentRcd,CSCAlignmentErrorExtendedRcd>
231  (&(*theMuonCSC), iSetup,
233  }
234  }
235 
236  // Create alignable tracker and muon
237  if (doTracker_) {
238  theAlignableTracker = new AlignableTracker( &(*theTracker), tTopo );
239  }
240 
241  if (doMuon_) {
242  theAlignableMuon = new AlignableMuon( &(*theMuonDT), &(*theMuonCSC) );
243  }
244 
245  if (useExtras_) {
247  }
248 
249  // Create alignment parameter builder
250  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob"
251  << "Creating AlignmentParameterBuilder";
252  edm::ParameterSet aliParamBuildCfg =
253  theParameterSet.getParameter<edm::ParameterSet>("ParameterBuilder");
254  AlignmentParameterBuilder alignmentParameterBuilder(theAlignableTracker,
257  aliParamBuildCfg );
258  // Fix alignables if requested
259  if (stNFixAlignables_>0) alignmentParameterBuilder.fixAlignables(stNFixAlignables_);
260 
261  // Get list of alignables
262  Alignables theAlignables = alignmentParameterBuilder.alignables();
263  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob"
264  << "got " << theAlignables.size() << " alignables";
265 
266  // Create AlignmentParameterStore
267  edm::ParameterSet aliParamStoreCfg =
269  theAlignmentParameterStore = new AlignmentParameterStore(theAlignables, aliParamStoreCfg);
270  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob"
271  << "AlignmentParameterStore created!";
272 
273  // Apply misalignment scenario to alignable tracker and muon if requested
274  // WARNING: this assumes scenarioConfig can be passed to both muon and tracker
276  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob"
277  << "Applying misalignment scenario to "
278  << (doTracker_ ? "tracker" : "")
279  << (doMuon_ ? (doTracker_ ? " and muon" : "muon") : ".");
280  edm::ParameterSet scenarioConfig
281  = theParameterSet.getParameter<edm::ParameterSet>( "MisalignmentScenario" );
282  if (doTracker_) {
283  TrackerScenarioBuilder scenarioBuilder( theAlignableTracker );
284  scenarioBuilder.applyScenario( scenarioConfig );
285  }
286  if (doMuon_) {
287  MuonScenarioBuilder muonScenarioBuilder( theAlignableMuon );
288  muonScenarioBuilder.applyScenario( scenarioConfig );
289  }
290  } else {
291  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob"
292  << "NOT applying misalignment scenario!";
293  }
294 
295  // Apply simple misalignment
296  const std::string sParSel(theParameterSet.getParameter<std::string>("parameterSelectorSimple"));
297  this->simpleMisalignment_(theAlignables, sParSel, stRandomShift_, stRandomRotation_, true);
298 
299  // Initialize alignment algorithm and integrated calibration and pass the latter to algorithm
300  theAlignmentAlgo->initialize( iSetup,
302  theAlignmentParameterStore );
303  for (auto iCal = theCalibrations.begin(); iCal != theCalibrations.end(); ++iCal) {
305  }
306  // Not all algorithms support calibrations - so do not pass empty vector
307  // and throw if non-empty and not supported:
309  throw cms::Exception("BadConfig") << "[AlignmentProducer::beginOfJob]\n"
310  << "Configured " << theCalibrations.size() << " calibration(s) "
311  << "for algorithm not supporting it.";
312  }
313 
314  for (std::vector<AlignmentMonitorBase*>::const_iterator monitor = theMonitors.begin();
315  monitor != theMonitors.end(); ++monitor) {
316  (*monitor)->beginOfJob(theAlignableTracker, theAlignableMuon, theAlignmentParameterStore);
317  }
318 }
T getParameter(std::string const &) const
const bool doMisalignmentScenario_
Builds a scenario from configuration and applies it to the alignable Muon.
std::vector< Alignable * > Alignables
Definition: Alignable.h:251
AlignmentAlgorithmBase * theAlignmentAlgo
void simpleMisalignment_(const Alignables &alivec, const std::string &selection, float shift, float rot, bool local)
Apply random shifts and rotations to selected alignables, according to configuration.
std::vector< IntegratedCalibrationBase * > theCalibrations
edm::ParameterSet theParameterSet
void createGeometries_(const edm::EventSetup &)
Create tracker and muon geometries.
const bool applyDbAlignment_
virtual bool addCalibrations(const Calibrations &)
AlignableExtras * theAlignableExtras
boost::shared_ptr< CSCGeometry > theMuonCSC
std::vector< AlignmentMonitorBase * > theMonitors
AlignmentParameterStore * theAlignmentParameterStore
virtual void initialize(const edm::EventSetup &setup, AlignableTracker *tracker, AlignableMuon *muon, AlignableExtras *extras, AlignmentParameterStore *store)=0
Call at beginning of job (must be implemented in derived class)
Definition: DetId.h:18
AlignableMuon * theAlignableMuon
const double stRandomShift_
const T & get() const
Definition: EventSetup.h:56
T const * product() const
Definition: ESHandle.h:86
AlignableTracker * theAlignableTracker
const AlignTransform & DetectorGlobalPosition(const Alignments &allGlobals, const DetId &id)
Builds a scenario from configuration and applies it to the alignable tracker.
const Alignments * globalPositions_
GlobalPositions that might be read from DB, NULL otherwise.
const double stRandomRotation_
void AlignmentProducer::beginRun ( const edm::Run run,
const edm::EventSetup setup 
)
virtual

Called at run start and calling algorithms beginRun.

Reimplemented from edm::EDLooperBase.

Definition at line 535 of file AlignmentProducer.cc.

References AlignmentAlgorithmBase::beginRun(), and theAlignmentAlgo.

536 {
537  theAlignmentAlgo->beginRun(setup); // do not forward edm::Run...
538 }
virtual void beginRun(const edm::EventSetup &setup)
called at begin of run
AlignmentAlgorithmBase * theAlignmentAlgo
void AlignmentProducer::createGeometries_ ( const edm::EventSetup iSetup)
private

Create tracker and muon geometries.

Definition at line 644 of file AlignmentProducer.cc.

References TrackerGeomBuilderFromGeometricDet::build(), CSCGeometryBuilderFromDDD::build(), DTGeometryBuilderFromDDD::build(), doMuon_, doTracker_, edm::EventSetup::get(), edm::ESHandle< class >::product(), theMuonCSC, theMuonDT, and theTracker.

Referenced by beginOfJob().

645 {
647  iSetup.get<IdealGeometryRecord>().get( cpv );
648 
649  if (doTracker_) {
650  edm::ESHandle<GeometricDet> geometricDet;
651  iSetup.get<IdealGeometryRecord>().get( geometricDet );
652 
654  iSetup.get<PTrackerParametersRcd>().get( ptp );
655 
656  edm::ESHandle<TrackerTopology> tTopoHandle;
657  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
658  const TrackerTopology* const tTopo = tTopoHandle.product();
659 
660  TrackerGeomBuilderFromGeometricDet trackerBuilder;
661  theTracker = boost::shared_ptr<TrackerGeometry>( trackerBuilder.build(&(*geometricDet), *ptp, tTopo ));
662  }
663 
664  if (doMuon_) {
666  iSetup.get<MuonNumberingRecord>().get(mdc);
667  DTGeometryBuilderFromDDD DTGeometryBuilder;
669  theMuonDT = boost::shared_ptr<DTGeometry>(new DTGeometry );
670  DTGeometryBuilder.build( theMuonDT, &(*cpv), *mdc);
671  theMuonCSC = boost::shared_ptr<CSCGeometry>( new CSCGeometry );
672  CSCGeometryBuilder.build( theMuonCSC, &(*cpv), *mdc );
673  }
674 }
boost::shared_ptr< TrackerGeometry > theTracker
TrackerGeometry * build(const GeometricDet *gd, const PTrackerParameters &ptp, const TrackerTopology *tTopo)
void build(boost::shared_ptr< CSCGeometry > geom, const DDCompactView *fv, const MuonDDDConstants &muonConstants)
Build the geometry.
boost::shared_ptr< CSCGeometry > theMuonCSC
boost::shared_ptr< DTGeometry > theMuonDT
void build(boost::shared_ptr< DTGeometry > theGeometry, const DDCompactView *cview, const MuonDDDConstants &muonConstants)
const T & get() const
Definition: EventSetup.h:56
T const * product() const
Definition: ESHandle.h:86
edm::EDLooper::Status AlignmentProducer::duringLoop ( const edm::Event event,
const edm::EventSetup setup 
)
virtual

Called at each event.

Implements edm::EDLooper.

Definition at line 473 of file AlignmentProducer.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, beamSpotTag_, clusterValueMapTag_, edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, edm::InputTag::encode(), ZMuMuCategoriesSequences_cff::eventInfo, edm::Event::getByLabel(), i, edm::EventBase::id(), AlignableExtras::initializeBeamSpot(), edm::EDLooperBase::kContinue, nevent_, readInSurveyRcds(), AlignmentAlgorithmBase::run(), theAlignableExtras, theAlignmentAlgo, theMonitors, and tjTkAssociationMapTag_.

475 {
476  ++nevent_;
477 
478  // reading in survey records
479  this->readInSurveyRcds(setup);
480 
481  // Printout event number
482  for ( int i=10; i<10000000; i*=10 )
483  if ( nevent_<10*i && (nevent_%i)==0 )
484  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::duringLoop"
485  << "Events processed: " << nevent_;
486 
487  // Retrieve trajectories and tracks from the event
488  // -> merely skip if collection is empty
490  if (event.getByLabel(tjTkAssociationMapTag_, m_TrajTracksMap)) {
491 
492  // Form pairs of trajectories and tracks
493  ConstTrajTrackPairCollection trajTracks;
494  for ( TrajTrackAssociationCollection::const_iterator iPair = m_TrajTracksMap->begin();
495  iPair != m_TrajTracksMap->end(); ++iPair) {
496  trajTracks.push_back( ConstTrajTrackPair( &(*(*iPair).key), &(*(*iPair).val) ) );
497  }
499  event.getByLabel(beamSpotTag_, beamSpot);
500 
501  if (nevent_==1 && theAlignableExtras) {
502  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::duringLoop"
503  << "initializing AlignableBeamSpot" << std::endl;
504  theAlignableExtras->initializeBeamSpot(beamSpot->x0(), beamSpot->y0(), beamSpot->z0(),
505  beamSpot->dxdz(), beamSpot->dydz());
506  }
507 
508  // Run the alignment algorithm with its input
509  const AliClusterValueMap *clusterValueMapPtr = 0;
510  if(clusterValueMapTag_.encode().size()){//check that the input tag is not empty
511  edm::Handle<AliClusterValueMap> clusterValueMap;
512  event.getByLabel(clusterValueMapTag_, clusterValueMap);
513  clusterValueMapPtr = &(*clusterValueMap);
514  }
515 
516  const AlignmentAlgorithmBase::EventInfo eventInfo(event.id(), trajTracks, *beamSpot,
517  clusterValueMapPtr);
518  theAlignmentAlgo->run(setup, eventInfo);
519 
520 
521  for (std::vector<AlignmentMonitorBase*>::const_iterator monitor = theMonitors.begin();
522  monitor != theMonitors.end(); ++monitor) {
523  (*monitor)->duringLoop(event, setup, trajTracks); // forward eventInfo?
524  }
525  } else {
526  edm::LogError("Alignment") << "@SUB=AlignmentProducer::duringLoop"
527  << "No track collection found: skipping event";
528  }
529 
530 
531  return kContinue;
532 }
int i
Definition: DBlmapReader.cc:9
virtual void run(const edm::EventSetup &setup, const EventInfo &eventInfo)=0
Run the algorithm (must be implemented in derived class)
AlignmentAlgorithmBase * theAlignmentAlgo
void readInSurveyRcds(const edm::EventSetup &)
read in survey records
std::string encode() const
Definition: InputTag.cc:164
define event information passed to algorithms
AlignableExtras * theAlignableExtras
const edm::InputTag clusterValueMapTag_
std::vector< AlignmentMonitorBase * > theMonitors
void initializeBeamSpot(double x, double y, double z, double dxdz, double dydz)
Initialize the alignable beam spot with the given parameters.
const edm::InputTag beamSpotTag_
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:418
const edm::InputTag tjTkAssociationMapTag_
std::vector< ConstTrajTrackPair > ConstTrajTrackPairCollection
edm::EventID id() const
Definition: EventBase.h:59
std::pair< const Trajectory *, const reco::Track * > ConstTrajTrackPair
void AlignmentProducer::endLuminosityBlock ( const edm::LuminosityBlock lumiBlock,
const edm::EventSetup setup 
)
virtual

Called at lumi block end, calling algorithm's endLuminosityBlock.

Reimplemented from edm::EDLooperBase.

Definition at line 567 of file AlignmentProducer.cc.

References AlignmentAlgorithmBase::endLuminosityBlock(), and theAlignmentAlgo.

569 {
570  theAlignmentAlgo->endLuminosityBlock(setup); // do not forward edm::LuminosityBlock
571 }
AlignmentAlgorithmBase * theAlignmentAlgo
virtual void endLuminosityBlock(const edm::EventSetup &setup)
called at end of luminosity block (no lumi block info passed yet)
void AlignmentProducer::endOfJob ( )
virtual

Called at end of job.

Reimplemented from edm::EDLooperBase.

Definition at line 322 of file AlignmentProducer.cc.

References AlignableExtras::beamSpot(), edm::ParameterSet::getParameter(), makeNonOverlappingRunRanges(), nevent_, AlignmentParameters::parameters(), cond::runnumber, saveApeToDB_, saveDeformationsToDB_, saveToDB_, AlignmentAlgorithmBase::setParametersForRunRange(), theAlignableExtras, theAlignmentAlgo, theCalibrations, theMonitors, theParameterSet, cond::timeTypeSpecs, and writeForRunRange().

323 {
324  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::endOfJob";
325 
326  for (std::vector<AlignmentMonitorBase*>::const_iterator monitor = theMonitors.begin(); monitor != theMonitors.end(); ++monitor) {
327  (*monitor)->endOfJob();
328  }
329 
330  if (0 == nevent_) {
331  edm::LogError("Alignment") << "@SUB=AlignmentProducer::endOfJob" << "Did not process any "
332  << "events in last loop, do not dare to store to DB.";
333  } else {
334 
335  // Expand run ranges and make them unique
336  edm::VParameterSet runRangeSelectionVPSet(theParameterSet.getParameter<edm::VParameterSet>("RunRangeSelection"));
337  RunRanges uniqueRunRanges(this->makeNonOverlappingRunRanges(runRangeSelectionVPSet));
338  if (uniqueRunRanges.empty()) { // create dummy IOV
339  const RunRange runRange(cond::timeTypeSpecs[cond::runnumber].beginValue,
341  uniqueRunRanges.push_back(runRange);
342  }
343 
344  std::vector<AlgebraicVector> beamSpotParameters;
345 
346  for (RunRanges::const_iterator iRunRange = uniqueRunRanges.begin();
347  iRunRange != uniqueRunRanges.end();
348  ++iRunRange) {
349 
351 
352  // Save alignments to database
354  this->writeForRunRange((*iRunRange).first);
355 
356  // Deal with extra alignables, e.g. beam spot
357  if (theAlignableExtras) {
359  if (!alis.empty()) {
360  BeamSpotAlignmentParameters *beamSpotAliPars = dynamic_cast<BeamSpotAlignmentParameters*>(alis[0]->alignmentParameters());
361  beamSpotParameters.push_back(beamSpotAliPars->parameters());
362  }
363  }
364  }
365 
366  if (theAlignableExtras) {
367  std::ostringstream bsOutput;
368 
369  std::vector<AlgebraicVector>::const_iterator itPar = beamSpotParameters.begin();
370  for (RunRanges::const_iterator iRunRange = uniqueRunRanges.begin();
371  iRunRange != uniqueRunRanges.end();
372  ++iRunRange, ++itPar) {
373  bsOutput << "Run range: " << (*iRunRange).first << " - " << (*iRunRange).second << "\n";
374  bsOutput << " Displacement: x=" << (*itPar)[0] << ", y=" << (*itPar)[1] << "\n";
375  bsOutput << " Slope: dx/dz=" << (*itPar)[2] << ", dy/dz=" << (*itPar)[3] << "\n";
376  }
377 
378  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::endOfJob"
379  << "Parameters for alignable beamspot:\n"
380  << bsOutput.str();
381  }
382 
383  for (auto iCal = theCalibrations.begin(); iCal != theCalibrations.end(); ++iCal) {
384  (*iCal)->endOfJob();
385  }
386 
387  }
388 }
const TimeTypeSpecs timeTypeSpecs[]
Definition: Time.cc:22
T getParameter(std::string const &) const
std::vector< Alignable * > Alignables
Definition: Alignable.h:251
AlignmentAlgorithmBase * theAlignmentAlgo
std::vector< ParameterSet > VParameterSet
Definition: ParameterSet.h:33
std::vector< IntegratedCalibrationBase * > theCalibrations
edm::ParameterSet theParameterSet
const AlgebraicVector & parameters(void) const
Get alignment parameters.
AlignableExtras * theAlignableExtras
virtual bool setParametersForRunRange(const RunRange &rr)
std::pair< RunNumber, RunNumber > RunRange
std::vector< AlignmentMonitorBase * > theMonitors
RunRanges makeNonOverlappingRunRanges(const edm::VParameterSet &RunRangeSelectionVPSet)
std::vector< RunRange > RunRanges
const bool saveDeformationsToDB_
void writeForRunRange(cond::Time_t time)
Alignables & beamSpot()
Return beam spot alignable as a vector with one element.
edm::EDLooper::Status AlignmentProducer::endOfLoop ( const edm::EventSetup iSetup,
unsigned int  iLoop 
)
virtual

Called at end of loop.

Implements edm::EDLooperBase.

Definition at line 438 of file AlignmentProducer.cc.

References edm::EDLooperBase::kContinue, edm::EDLooperBase::kStop, nevent_, AlignmentAlgorithmBase::terminate(), theAlignmentAlgo, theCalibrations, theMaxLoops, and theMonitors.

439 {
440 
441  if (0 == nevent_) {
442  // beginOfJob is usually called by the framework in the first event of the first loop
443  // (a hack: beginOfJob needs the EventSetup that is not well defined without an event)
444  // and the algorithms rely on the initialisations done in beginOfJob. We cannot call
445  // this->beginOfJob(iSetup); here either since that will access the EventSetup to get
446  // some geometry information that is not defined either without having seen an event.
447  edm::LogError("Alignment") << "@SUB=AlignmentProducer::endOfLoop"
448  << "Did not process any events in loop " << iLoop
449  << ", stop processing without terminating algorithm.";
450  return kStop;
451  }
452 
453  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::endOfLoop"
454  << "Ending loop " << iLoop << ", terminating algorithm.";
455 
456  theAlignmentAlgo->terminate(iSetup);
457  // FIXME: Should this be done in algorithm::terminate(const edm::EventSetup& iSetup)??
458  for (auto iCal = theCalibrations.begin(); iCal != theCalibrations.end(); ++iCal) {
459  (*iCal)->endOfLoop();
460  }
461 
462  for (std::vector<AlignmentMonitorBase*>::const_iterator monitor = theMonitors.begin(); monitor != theMonitors.end(); ++monitor) {
463  (*monitor)->endOfLoop(iSetup);
464  }
465 
466  if ( iLoop == theMaxLoops-1 || iLoop >= theMaxLoops ) return kStop;
467  else return kContinue;
468 }
AlignmentAlgorithmBase * theAlignmentAlgo
virtual void terminate(const edm::EventSetup &iSetup)=0
Call at end of each loop (must be implemented in derived class)
std::vector< IntegratedCalibrationBase * > theCalibrations
const unsigned int theMaxLoops
std::vector< AlignmentMonitorBase * > theMonitors
void AlignmentProducer::endRun ( const edm::Run run,
const edm::EventSetup setup 
)
virtual

Called at run end - currently reading TkFittedLasBeam if an InpuTag is given for that.

Reimplemented from edm::EDLooperBase.

Definition at line 541 of file AlignmentProducer.cc.

References edm::InputTag::encode(), AlignmentAlgorithmBase::endRun(), edm::Run::getByLabel(), edm::RunBase::id(), theAlignmentAlgo, and tkLasBeamTag_.

542 {
543  // call with or without las beam info...
545  if (tkLasBeamTag_.encode().size()) { // non-empty InputTag
548  run.getByLabel(tkLasBeamTag_, lasBeams);
549  run.getByLabel(tkLasBeamTag_, tsoses);
550 
551  theAlignmentAlgo->endRun(EndRunInfo(run.id(), &(*lasBeams), &(*tsoses)), setup);
552  } else {
553  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::endRun"
554  << "No Tk LAS beams to forward to algorithm.";
555  theAlignmentAlgo->endRun(EndRunInfo(run.id(), 0, 0), setup);
556  }
557 }
bool getByLabel(std::string const &label, Handle< PROD > &result) const
Definition: Run.h:219
RunID const & id() const
Definition: RunBase.h:39
AlignmentAlgorithmBase * theAlignmentAlgo
const edm::InputTag tkLasBeamTag_
std::string encode() const
Definition: InputTag.cc:164
virtual void endRun(const EndRunInfo &runInfo, const edm::EventSetup &setup)
called at end of run - order of arguments like in EDProducer etc.
AlignmentAlgorithmBase::EndRunInfo EndRunInfo
define run information passed to algorithms (in endRun)
AlignmentProducer::RunRanges AlignmentProducer::makeNonOverlappingRunRanges ( const edm::VParameterSet RunRangeSelectionVPSet)
private

Definition at line 968 of file AlignmentProducer.cc.

References cond::TimeTypeSpecs::beginValue, cond::TimeTypeSpecs::endValue, plotBeamSpotDB::first, i, cond::runnumber, groupFilesInBlocks::temp, cond::timeTypeSpecs, and edm::tokenize().

Referenced by endOfJob().

969 {
970  static bool oldRunRangeSelectionWarning = false;
971 
974 
975  RunRanges uniqueRunRanges;
976  if (!RunRangeSelectionVPSet.empty()) {
977 
978  std::map<RunNumber,RunNumber> uniqueFirstRunNumbers;
979 
980  for (std::vector<edm::ParameterSet>::const_iterator ipset = RunRangeSelectionVPSet.begin();
981  ipset != RunRangeSelectionVPSet.end();
982  ++ipset) {
983  const std::vector<std::string> RunRangeStrings = (*ipset).getParameter<std::vector<std::string> >("RunRanges");
984  for (std::vector<std::string>::const_iterator irange = RunRangeStrings.begin();
985  irange != RunRangeStrings.end();
986  ++irange) {
987 
988  if ((*irange).find(':')==std::string::npos) {
989 
990  RunNumber first = beginValue;
991  long int temp = strtol((*irange).c_str(), 0, 0);
992  if (temp!=-1) first = temp;
993  uniqueFirstRunNumbers[first] = first;
994 
995  } else {
996 
997  if (!oldRunRangeSelectionWarning) {
998  edm::LogWarning("BadConfig") << "@SUB=AlignmentProducer::makeNonOverlappingRunRanges"
999  << "Config file contains old format for 'RunRangeSelection'. Only the start run\n"
1000  << "number is used internally. The number of the last run is ignored and can be\n"
1001  << "safely removed from the config file.\n";
1002  oldRunRangeSelectionWarning = true;
1003  }
1004 
1005  std::vector<std::string> tokens = edm::tokenize(*irange, ":");
1006  long int temp;
1007  RunNumber first = beginValue;
1008  temp = strtol(tokens[0].c_str(), 0, 0);
1009  if (temp!=-1) first = temp;
1010  uniqueFirstRunNumbers[first] = first;
1011  }
1012  }
1013  }
1014 
1015  for (std::map<RunNumber,RunNumber>::iterator iFirst = uniqueFirstRunNumbers.begin();
1016  iFirst!=uniqueFirstRunNumbers.end();
1017  ++iFirst) {
1018  uniqueRunRanges.push_back(std::pair<RunNumber,RunNumber>((*iFirst).first, endValue));
1019  }
1020  for (unsigned int i = 0;i<uniqueRunRanges.size()-1;++i) {
1021  uniqueRunRanges[i].second = uniqueRunRanges[i+1].first - 1;
1022  }
1023 
1024  } else {
1025 
1026  uniqueRunRanges.push_back(std::pair<RunNumber,RunNumber>(beginValue, endValue));
1027 
1028  }
1029 
1030  return uniqueRunRanges;
1031 }
const TimeTypeSpecs timeTypeSpecs[]
Definition: Time.cc:22
int i
Definition: DBlmapReader.cc:9
Time_t beginValue
Definition: Time.h:45
std::vector< RunRange > RunRanges
std::vector< std::string > tokenize(std::string const &input, std::string const &separator)
breaks the input string into tokens, delimited by the separator
Definition: Parse.cc:57
Time_t endValue
Definition: Time.h:46
cond::RealTimeType< cond::runnumber >::type RunNumber
boost::shared_ptr< CSCGeometry > AlignmentProducer::produceCSC ( const MuonGeometryRecord iRecord)
virtual

Produce the muon CSC geometry.

Definition at line 190 of file AlignmentProducer.cc.

References theMuonCSC.

Referenced by AlignmentProducer().

191 {
192  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::produceCSC";
193  return theMuonCSC;
194 }
boost::shared_ptr< CSCGeometry > theMuonCSC
boost::shared_ptr< DTGeometry > AlignmentProducer::produceDT ( const MuonGeometryRecord iRecord)
virtual

Produce the muon DT geometry.

Definition at line 181 of file AlignmentProducer.cc.

References theMuonDT.

Referenced by AlignmentProducer().

182 {
183  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::produceDT";
184  return theMuonDT;
185 }
boost::shared_ptr< DTGeometry > theMuonDT
boost::shared_ptr< TrackerGeometry > AlignmentProducer::produceTracker ( const TrackerDigiGeometryRecord iRecord)
virtual

Produce the tracker geometry.

Definition at line 172 of file AlignmentProducer.cc.

References theTracker.

Referenced by AlignmentProducer().

173 {
174  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::produceTracker";
175  return theTracker;
176 }
boost::shared_ptr< TrackerGeometry > theTracker
void AlignmentProducer::readInSurveyRcds ( const edm::EventSetup iSetup)
private

read in survey records

Definition at line 709 of file AlignmentProducer.cc.

References addSurveyInfo_(), edm::ESWatcher< T >::check(), AlignableMuon::CSCEndcaps(), doMuon_, doTracker_, AlignableMuon::DTBarrel(), edm::EventSetup::get(), theAlignableMuon, theAlignableTracker, theSurveyErrors, theSurveyIndex, theSurveyValues, useSurvey_, watchTkSurveyErrRcd_, and watchTkSurveyRcd_.

Referenced by duringLoop().

709  {
710 
711  // Get Survey Rcds and add Survey Info
712  if ( doTracker_ && useSurvey_ ){
713  bool tkSurveyBool = watchTkSurveyRcd_.check(iSetup);
714  bool tkSurveyErrBool = watchTkSurveyErrRcd_.check(iSetup);
715  edm::LogInfo("Alignment") << "watcher tksurveyrcd: " << tkSurveyBool;
716  edm::LogInfo("Alignment") << "watcher tksurveyerrrcd: " << tkSurveyErrBool;
717  if ( tkSurveyBool || tkSurveyErrBool){
718 
719  edm::LogInfo("Alignment") << "ADDING THE SURVEY INFORMATION";
721  edm::ESHandle<SurveyErrors> surveyErrors;
722 
723  iSetup.get<TrackerSurveyRcd>().get(surveys);
724  iSetup.get<TrackerSurveyErrorExtendedRcd>().get(surveyErrors);
725 
726  theSurveyIndex = 0;
727  theSurveyValues = &*surveys;
728  theSurveyErrors = &*surveyErrors;
730  }
731  }
732 
733  if ( doMuon_ && useSurvey_) {
734  bool DTSurveyBool = watchTkSurveyRcd_.check(iSetup);
735  bool DTSurveyErrBool = watchTkSurveyErrRcd_.check(iSetup);
736  bool CSCSurveyBool = watchTkSurveyRcd_.check(iSetup);
737  bool CSCSurveyErrBool = watchTkSurveyErrRcd_.check(iSetup);
738 
739  if ( DTSurveyBool || DTSurveyErrBool || CSCSurveyBool || CSCSurveyErrBool ){
740  edm::ESHandle<Alignments> dtSurveys;
741  edm::ESHandle<SurveyErrors> dtSurveyErrors;
742  edm::ESHandle<Alignments> cscSurveys;
743  edm::ESHandle<SurveyErrors> cscSurveyErrors;
744 
745  iSetup.get<DTSurveyRcd>().get(dtSurveys);
746  iSetup.get<DTSurveyErrorExtendedRcd>().get(dtSurveyErrors);
747  iSetup.get<CSCSurveyRcd>().get(cscSurveys);
748  iSetup.get<CSCSurveyErrorExtendedRcd>().get(cscSurveyErrors);
749 
750  theSurveyIndex = 0;
751  theSurveyValues = &*dtSurveys;
752  theSurveyErrors = &*dtSurveyErrors;
753  std::vector<Alignable*> barrels = theAlignableMuon->DTBarrel();
754  for (std::vector<Alignable*>::const_iterator iter = barrels.begin(); iter != barrels.end(); ++iter) {
755  addSurveyInfo_(*iter);
756  }
757 
758  theSurveyIndex = 0;
759  theSurveyValues = &*cscSurveys;
760  theSurveyErrors = &*cscSurveyErrors;
761  std::vector<Alignable*> endcaps = theAlignableMuon->CSCEndcaps();
762  for (std::vector<Alignable*>::const_iterator iter = endcaps.begin(); iter != endcaps.end(); ++iter) {
763  addSurveyInfo_(*iter);
764  }
765  }
766  }
767 
768 }
align::Alignables DTBarrel()
const Alignments * theSurveyValues
void addSurveyInfo_(Alignable *)
Add survey info to an alignable.
align::Alignables CSCEndcaps()
unsigned int theSurveyIndex
const SurveyErrors * theSurveyErrors
AlignableMuon * theAlignableMuon
const T & get() const
Definition: EventSetup.h:56
bool check(const edm::EventSetup &iSetup)
Definition: ESWatcher.h:57
edm::ESWatcher< TrackerSurveyRcd > watchTkSurveyRcd_
AlignableTracker * theAlignableTracker
edm::ESWatcher< TrackerSurveyErrorExtendedRcd > watchTkSurveyErrRcd_
void AlignmentProducer::simpleMisalignment_ ( const Alignables alivec,
const std::string &  selection,
float  shift,
float  rot,
bool  local 
)
private

Apply random shifts and rotations to selected alignables, according to configuration.

Definition at line 575 of file AlignmentProducer.cc.

References funct::abs(), Alignable::alignmentParameters(), AlignmentParameterSelector::convertParamSel(), RigidBodyAlignmentParameters::dalpha, RigidBodyAlignmentParameters::dbeta, RigidBodyAlignmentParameters::dgamma, RigidBodyAlignmentParameters::dx, RigidBodyAlignmentParameters::dy, RigidBodyAlignmentParameters::dz, Exception, Alignable::move(), RigidBodyAlignmentParameters::N_PARAM, convertSQLitetoXML_cfg::output, alignCSCRings::r, random, Alignable::rotateInGlobalFrame(), Alignable::rotateInLocalFrame(), indexGen::s2, AlignmentParameters::selector(), Alignable::surface(), AlignableSurface::toGlobal(), and align::toMatrix().

Referenced by beginOfJob().

577 {
578 
579  std::ostringstream output; // collecting output
580 
581  if (shift > 0. || rot > 0.) {
582  output << "Adding random flat shift of max size " << shift
583  << " and adding random flat rotation of max size " << rot <<" to ";
584 
585  std::vector<bool> commSel(0);
586  if (selection != "-1") {
587  AlignmentParameterSelector aSelector(0,0); // no alignable needed here...
588  const std::vector<char> cSel(aSelector.convertParamSel(selection));
589  if (cSel.size() < RigidBodyAlignmentParameters::N_PARAM) {
590  throw cms::Exception("BadConfig")
591  << "[AlignmentProducer::simpleMisalignment_]\n"
592  << "Expect selection string '" << selection << "' to be at least of length "
593  << RigidBodyAlignmentParameters::N_PARAM << " or to be '-1'.\n"
594  << "(Most probably you have to adjust the parameter 'parameterSelectorSimple'.)";
595  }
596  for (std::vector<char>::const_iterator cIter = cSel.begin(); cIter != cSel.end(); ++cIter) {
597  commSel.push_back(*cIter == '0' ? false : true);
598  }
599  output << "parameters defined by (" << selection
600  << "), representing (x,y,z,alpha,beta,gamma),";
601  } else {
602  output << "the active parameters of each alignable,";
603  }
604  output << " in " << (local ? "local" : "global") << " frame.";
605 
606  for (std::vector<Alignable*>::const_iterator it = alivec.begin(); it != alivec.end(); ++it) {
607  Alignable* ali=(*it);
608  std::vector<bool> mysel(commSel.empty() ? ali->alignmentParameters()->selector() : commSel);
609 
610  if (std::abs(shift)>0.00001) {
611  double s0 = 0., s1 = 0., s2 = 0.;
612  if (mysel[RigidBodyAlignmentParameters::dx]) s0 = shift * double(random()%1000-500)/500.;
613  if (mysel[RigidBodyAlignmentParameters::dy]) s1 = shift * double(random()%1000-500)/500.;
614  if (mysel[RigidBodyAlignmentParameters::dz]) s2 = shift * double(random()%1000-500)/500.;
615 
616  if (local) ali->move( ali->surface().toGlobal(align::LocalVector(s0,s1,s2)) );
617  else ali->move( align::GlobalVector(s0,s1,s2) );
618 
619  //AlignmentPositionError ape(dx,dy,dz);
620  //ali->addAlignmentPositionError(ape);
621  }
622 
623  if (std::abs(rot)>0.00001) {
625  if (mysel[RigidBodyAlignmentParameters::dalpha]) r(1)=rot*double(random()%1000-500)/500.;
626  if (mysel[RigidBodyAlignmentParameters::dbeta]) r(2)=rot*double(random()%1000-500)/500.;
627  if (mysel[RigidBodyAlignmentParameters::dgamma]) r(3)=rot*double(random()%1000-500)/500.;
628 
629  const align::RotationType mrot = align::toMatrix(r);
630  if (local) ali->rotateInLocalFrame(mrot);
631  else ali->rotateInGlobalFrame(mrot);
632 
633  //ali->addAlignmentPositionErrorFromRotation(mrot);
634  }
635  } // end loop on alignables
636  } else {
637  output << "No simple misalignment added!";
638  }
639  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::simpleMisalignment_" << output.str();
640 }
virtual void rotateInLocalFrame(const RotationType &rotation)
Rotation intepreted in the local reference frame.
Definition: Alignable.cc:120
selection
main part
Definition: corrVsCorr.py:98
virtual void move(const GlobalVector &displacement)=0
Movement with respect to the global reference frame.
const std::vector< bool > & selector(void) const
Get alignment parameter selector vector.
TRandom random
Definition: MVATrainer.cc:138
AlignmentParameters * alignmentParameters() const
Get the AlignmentParameters.
Definition: Alignable.h:57
tuple s2
Definition: indexGen.py:106
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:131
AlgebraicVector EulerAngles
Definition: Definitions.h:36
virtual void rotateInGlobalFrame(const RotationType &rotation)=0
align::GlobalPoints toGlobal(const align::LocalPoints &) const
Return in global coord given a set of local points.
RotationType toMatrix(const EulerAngles &)
Convert rotation angles about x-, y-, z-axes to matrix.
Definition: Utilities.cc:40
static unsigned int const shift
void AlignmentProducer::startingNewLoop ( unsigned int  iLoop)
virtual

Called at beginning of loop.

Implements edm::EDLooperBase.

Definition at line 392 of file AlignmentProducer.cc.

References AlignableTracker::alignmentErrors(), AlignableTracker::alignments(), GeometryAligner::applyAlignments(), GeometryAligner::attachSurfaceDeformations(), AlignableMuon::cscAlignmentErrorsExtended(), AlignableMuon::cscAlignments(), doMuon_, doTracker_, AlignableMuon::dtAlignmentErrorsExtended(), AlignableMuon::dtAlignments(), nevent_, AlignmentAlgorithmBase::startNewLoop(), Alignable::surfaceDeformations(), theAlignableMuon, theAlignableTracker, theAlignmentAlgo, theCalibrations, and theMonitors.

393 {
394  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::startingNewLoop"
395  << "Starting loop number " << iLoop;
396 
397  nevent_ = 0;
398 
400  // FIXME: Should this be done in algorithm::startNewLoop()??
401  for (auto iCal = theCalibrations.begin(); iCal != theCalibrations.end(); ++iCal) {
402  (*iCal)->startNewLoop();
403  }
404 
405  for (std::vector<AlignmentMonitorBase*>::const_iterator monitor = theMonitors.begin(); monitor != theMonitors.end(); ++monitor) {
406  (*monitor)->startingNewLoop();
407  }
408 
409  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::startingNewLoop"
410  << "Now physically apply alignments to geometry...";
411 
412 
413  // Propagate changes to reconstruction geometry (from initialisation or iteration)
414  GeometryAligner aligner;
415  if ( doTracker_ ) {
416  std::auto_ptr<Alignments> alignments(theAlignableTracker->alignments());
417  std::auto_ptr<AlignmentErrorsExtended> alignmentErrors(theAlignableTracker->alignmentErrors());
418  aligner.applyAlignments<TrackerGeometry>( &(*theTracker),&(*alignments),&(*alignmentErrors), AlignTransform() ); // don't apply global a second time!
419  std::auto_ptr<AlignmentSurfaceDeformations> aliDeforms(theAlignableTracker->surfaceDeformations());
420  aligner.attachSurfaceDeformations<TrackerGeometry>(&(*theTracker), &(*aliDeforms));
421 
422  }
423  if ( doMuon_ ) {
424  std::auto_ptr<Alignments> dtAlignments( theAlignableMuon->dtAlignments());
425  std::auto_ptr<AlignmentErrorsExtended> dtAlignmentErrorsExtended( theAlignableMuon->dtAlignmentErrorsExtended());
426  std::auto_ptr<Alignments> cscAlignments( theAlignableMuon->cscAlignments());
427  std::auto_ptr<AlignmentErrorsExtended> cscAlignmentErrorsExtended( theAlignableMuon->cscAlignmentErrorsExtended());
428 
429  aligner.applyAlignments<DTGeometry>( &(*theMuonDT), &(*dtAlignments), &(*dtAlignmentErrorsExtended), AlignTransform() ); // don't apply global a second time!
430  aligner.applyAlignments<CSCGeometry>( &(*theMuonCSC), &(*cscAlignments), &(*cscAlignmentErrorsExtended), AlignTransform() ); // nope!
431  }
432 }
AlignmentAlgorithmBase * theAlignmentAlgo
AlignmentErrorsExtended * alignmentErrors() const
Return alignment errors, sorted by DetId.
Class to update a given geometry with a set of alignments.
std::vector< IntegratedCalibrationBase * > theCalibrations
Alignments * dtAlignments()
void attachSurfaceDeformations(C *geometry, const AlignmentSurfaceDeformations *surfaceDeformations)
std::vector< AlignmentMonitorBase * > theMonitors
AlignmentSurfaceDeformations * surfaceDeformations() const
Return surface deformations, sorted by DetId.
Definition: Alignable.cc:245
AlignmentErrorsExtended * cscAlignmentErrorsExtended()
AlignableMuon * theAlignableMuon
void applyAlignments(C *geometry, const Alignments *alignments, const AlignmentErrorsExtended *alignmentErrors, const AlignTransform &globalCoordinates)
Alignments * cscAlignments()
AlignmentErrorsExtended * dtAlignmentErrorsExtended()
AlignableTracker * theAlignableTracker
Alignments * alignments() const
Return alignments, sorted by DetId.
void AlignmentProducer::writeDB ( Alignments alignments,
const std::string &  alignRcd,
AlignmentErrorsExtended alignmentErrors,
const std::string &  errRcd,
const AlignTransform globalCoordinates,
cond::Time_t  time 
) const
private

Write alignment and/or errors to DB for record names (removes *globalCoordinates before writing if non-null...). Takes over ownership of alignments and alignmentErrrors.

Definition at line 890 of file AlignmentProducer.cc.

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

Referenced by writeForRunRange().

896 {
897  Alignments * tempAlignments = alignments;
898  AlignmentErrorsExtended * tempAlignmentErrorsExtended = alignmentErrors;
899 
900  // Call service
902  if (!poolDb.isAvailable()) { // Die if not available
903  delete tempAlignments; // promised to take over ownership...
904  delete tempAlignmentErrorsExtended; // dito
905  throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
906  }
907 
908  if (globalCoordinates // happens only if (applyDbAlignment_ == true)
909  && globalCoordinates->transform() != AlignTransform::Transform::Identity) {
910 
911  tempAlignments = new Alignments(); // temporary storage for
912  tempAlignmentErrorsExtended = new AlignmentErrorsExtended(); // final alignments and errors
913 
914  GeometryAligner aligner;
915  aligner.removeGlobalTransform(alignments, alignmentErrors,
916  *globalCoordinates,
917  tempAlignments, tempAlignmentErrorsExtended);
918 
919  delete alignments; // have to delete original alignments
920  delete alignmentErrors; // same thing for the errors
921 
922  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::writeDB"
923  << "globalCoordinates removed from alignments (" << alignRcd
924  << ") and errors (" << alignRcd << ").";
925  }
926 
927  if (saveToDB_) {
928  edm::LogInfo("Alignment") << "Writing Alignments for run " << time
929  << " to " << alignRcd << ".";
930  poolDb->writeOne<Alignments>(tempAlignments, time, alignRcd);
931  } else { // poolDb->writeOne(..) takes over 'alignments' ownership,...
932  delete tempAlignments; // ...otherwise we have to delete, as promised!
933  }
934 
935  if (saveApeToDB_) {
936  edm::LogInfo("Alignment") << "Writing AlignmentErrorsExtended for run " << time
937  << " to " << errRcd << ".";
938  poolDb->writeOne<AlignmentErrorsExtended>(tempAlignmentErrorsExtended, time, errRcd);
939  } else { // poolDb->writeOne(..) takes over 'alignmentErrors' ownership,...
940  delete tempAlignmentErrorsExtended; // ...otherwise we have to delete, as promised!
941  }
942 }
Class to update a given geometry with a set of alignments.
void removeGlobalTransform(const Alignments *alignments, const AlignmentErrorsExtended *alignmentErrors, const AlignTransform &globalCoordinates, Alignments *newAlignments, AlignmentErrorsExtended *newAlignmentErrorsExtended)
bool isAvailable() const
Definition: Service.h:46
void writeOne(T *payload, Time_t time, const std::string &recordName, bool withlogging=false)
Transform transform() const
void AlignmentProducer::writeDB ( AlignmentSurfaceDeformations alignmentSurfaceDeformations,
const std::string &  surfaceDeformationRcd,
cond::Time_t  time 
) const
private

Write surface deformations (bows & kinks) to DB for given record name Takes over ownership of alignmentsurfaceDeformations.

Definition at line 946 of file AlignmentProducer.cc.

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

949 {
950  // Call service
952  if (!poolDb.isAvailable()) { // Die if not available
953  delete alignmentSurfaceDeformations; // promised to take over ownership...
954  throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
955  }
956 
957  if (saveDeformationsToDB_) {
958  edm::LogInfo("Alignment") << "Writing AlignmentSurfaceDeformations for run " << time
959  << " to " << surfaceDeformationRcd << ".";
960  poolDb->writeOne<AlignmentSurfaceDeformations>(alignmentSurfaceDeformations, time,
961  surfaceDeformationRcd);
962  } else { // poolDb->writeOne(..) takes over 'surfaceDeformation' ownership,...
963  delete alignmentSurfaceDeformations; // ...otherwise we have to delete, as promised!
964  }
965 }
bool isAvailable() const
Definition: Service.h:46
void writeOne(T *payload, Time_t time, const std::string &recordName, bool withlogging=false)
const bool saveDeformationsToDB_
void AlignmentProducer::writeForRunRange ( cond::Time_t  time)
private

Definition at line 845 of file AlignmentProducer.cc.

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

Referenced by endOfJob().

846 {
847  if ( doTracker_ ) { // first tracker
848  const AlignTransform *trackerGlobal = 0; // will be 'removed' from constants
849  if (globalPositions_) { // i.e. applied before in applyDB
852  }
853 
854  Alignments *alignments = theAlignableTracker->alignments();
856  this->writeDB(alignments, "TrackerAlignmentRcd",
857  alignmentErrors, "TrackerAlignmentErrorExtendedRcd", trackerGlobal,
858  time);
859  }
860 
861  if ( doMuon_ ) { // now muon
862  const AlignTransform *muonGlobal = 0; // will be 'removed' from constants
863  if (globalPositions_) { // i.e. applied before in applyDB
865  DetId(DetId::Muon));
866  }
867  // Get alignments+errors, first DT - ownership taken over by writeDB(..), so no delete
868  Alignments *alignments = theAlignableMuon->dtAlignments();
870  this->writeDB(alignments, "DTAlignmentRcd",
871  alignmentErrors, "DTAlignmentErrorExtendedRcd", muonGlobal,
872  time);
873 
874  // Get alignments+errors, now CSC - ownership taken over by writeDB(..), so no delete
875  alignments = theAlignableMuon->cscAlignments();
876  alignmentErrors = theAlignableMuon->cscAlignmentErrorsExtended();
877  this->writeDB(alignments, "CSCAlignmentRcd",
878  alignmentErrors, "CSCAlignmentErrorExtendedRcd", muonGlobal,
879  time);
880  }
881 
882  // Save surface deformations to database
884  AlignmentSurfaceDeformations *alignmentSurfaceDeformations = theAlignableTracker->surfaceDeformations();
885  this->writeDB(alignmentSurfaceDeformations, "TrackerSurfaceDeformationRcd", time);
886  }
887 }
AlignmentErrorsExtended * alignmentErrors() const
Return alignment errors, sorted by DetId.
void writeDB(Alignments *alignments, const std::string &alignRcd, AlignmentErrorsExtended *alignmentErrors, const std::string &errRcd, const AlignTransform *globalCoordinates, cond::Time_t time) const
Alignments * dtAlignments()
AlignmentSurfaceDeformations * surfaceDeformations() const
Return surface deformations, sorted by DetId.
Definition: Alignable.cc:245
AlignmentErrorsExtended * cscAlignmentErrorsExtended()
Definition: DetId.h:18
AlignableMuon * theAlignableMuon
const bool saveDeformationsToDB_
Alignments * cscAlignments()
AlignmentErrorsExtended * dtAlignmentErrorsExtended()
AlignableTracker * theAlignableTracker
const AlignTransform & DetectorGlobalPosition(const Alignments &allGlobals, const DetId &id)
Alignments * alignments() const
Return alignments, sorted by DetId.
const Alignments * globalPositions_
GlobalPositions that might be read from DB, NULL otherwise.

Member Data Documentation

const bool AlignmentProducer::applyDbAlignment_
private

Definition at line 181 of file AlignmentProducer.h.

Referenced by beginOfJob().

const edm::InputTag AlignmentProducer::beamSpotTag_
private

Definition at line 189 of file AlignmentProducer.h.

Referenced by duringLoop().

const bool AlignmentProducer::checkDbAlignmentValidity_
private

Definition at line 181 of file AlignmentProducer.h.

Referenced by applyDB().

const edm::InputTag AlignmentProducer::clusterValueMapTag_
private

Definition at line 191 of file AlignmentProducer.h.

Referenced by duringLoop().

const bool AlignmentProducer::doMisalignmentScenario_
private

Definition at line 182 of file AlignmentProducer.h.

Referenced by beginOfJob().

const bool AlignmentProducer::doMuon_
private
const bool AlignmentProducer::doTracker_
private
const Alignments* AlignmentProducer::globalPositions_
private

GlobalPositions that might be read from DB, NULL otherwise.

Definition at line 170 of file AlignmentProducer.h.

Referenced by beginOfJob(), writeForRunRange(), and ~AlignmentProducer().

int AlignmentProducer::nevent_
private

Definition at line 172 of file AlignmentProducer.h.

Referenced by duringLoop(), endOfJob(), endOfLoop(), and startingNewLoop().

const bool AlignmentProducer::saveApeToDB_
private

Definition at line 183 of file AlignmentProducer.h.

Referenced by endOfJob(), and writeDB().

const bool AlignmentProducer::saveDeformationsToDB_
private

Definition at line 183 of file AlignmentProducer.h.

Referenced by endOfJob(), writeDB(), and writeForRunRange().

const bool AlignmentProducer::saveToDB_
private

Definition at line 183 of file AlignmentProducer.h.

Referenced by endOfJob(), and writeDB().

const int AlignmentProducer::stNFixAlignables_
private

Definition at line 179 of file AlignmentProducer.h.

Referenced by beginOfJob().

const double AlignmentProducer::stRandomRotation_
private

Definition at line 180 of file AlignmentProducer.h.

Referenced by beginOfJob().

const double AlignmentProducer::stRandomShift_
private

Definition at line 180 of file AlignmentProducer.h.

Referenced by beginOfJob().

AlignableExtras* AlignmentProducer::theAlignableExtras
private

Definition at line 162 of file AlignmentProducer.h.

Referenced by beginOfJob(), duringLoop(), endOfJob(), and ~AlignmentProducer().

AlignableMuon* AlignmentProducer::theAlignableMuon
private
AlignableTracker* AlignmentProducer::theAlignableTracker
private
AlignmentAlgorithmBase* AlignmentProducer::theAlignmentAlgo
private
AlignmentParameterStore* AlignmentProducer::theAlignmentParameterStore
private

Definition at line 158 of file AlignmentProducer.h.

Referenced by beginOfJob(), and ~AlignmentProducer().

std::vector<IntegratedCalibrationBase*> AlignmentProducer::theCalibrations
private
const unsigned int AlignmentProducer::theMaxLoops
private

Definition at line 177 of file AlignmentProducer.h.

Referenced by endOfLoop().

std::vector<AlignmentMonitorBase*> AlignmentProducer::theMonitors
private
boost::shared_ptr<CSCGeometry> AlignmentProducer::theMuonCSC
private

Definition at line 168 of file AlignmentProducer.h.

Referenced by beginOfJob(), createGeometries_(), and produceCSC().

boost::shared_ptr<DTGeometry> AlignmentProducer::theMuonDT
private

Definition at line 167 of file AlignmentProducer.h.

Referenced by beginOfJob(), createGeometries_(), and produceDT().

edm::ParameterSet AlignmentProducer::theParameterSet
private

Definition at line 173 of file AlignmentProducer.h.

Referenced by beginOfJob(), and endOfJob().

const SurveyErrors* AlignmentProducer::theSurveyErrors
private

Definition at line 155 of file AlignmentProducer.h.

Referenced by addSurveyInfo_(), and readInSurveyRcds().

unsigned int AlignmentProducer::theSurveyIndex
private

Definition at line 153 of file AlignmentProducer.h.

Referenced by addSurveyInfo_(), and readInSurveyRcds().

const Alignments* AlignmentProducer::theSurveyValues
private

Definition at line 154 of file AlignmentProducer.h.

Referenced by addSurveyInfo_(), and readInSurveyRcds().

boost::shared_ptr<TrackerGeometry> AlignmentProducer::theTracker
private

Definition at line 166 of file AlignmentProducer.h.

Referenced by beginOfJob(), createGeometries_(), and produceTracker().

const edm::InputTag AlignmentProducer::tjTkAssociationMapTag_
private

Definition at line 188 of file AlignmentProducer.h.

Referenced by duringLoop().

const edm::InputTag AlignmentProducer::tkLasBeamTag_
private

Definition at line 190 of file AlignmentProducer.h.

Referenced by endRun().

const bool AlignmentProducer::useExtras_
private

Definition at line 184 of file AlignmentProducer.h.

Referenced by beginOfJob().

const bool AlignmentProducer::useSurvey_
private

Definition at line 185 of file AlignmentProducer.h.

Referenced by readInSurveyRcds().

edm::ESWatcher<CSCSurveyErrorExtendedRcd> AlignmentProducer::watchCSCSurveyErrRcd_
private

Definition at line 199 of file AlignmentProducer.h.

edm::ESWatcher<CSCSurveyRcd> AlignmentProducer::watchCSCSurveyRcd_
private

Definition at line 198 of file AlignmentProducer.h.

edm::ESWatcher<DTSurveyErrorExtendedRcd> AlignmentProducer::watchDTSurveyErrRcd_
private

Definition at line 197 of file AlignmentProducer.h.

edm::ESWatcher<DTSurveyRcd> AlignmentProducer::watchDTSurveyRcd_
private

Definition at line 196 of file AlignmentProducer.h.

edm::ESWatcher<TrackerSurveyErrorExtendedRcd> AlignmentProducer::watchTkSurveyErrRcd_
private

Definition at line 195 of file AlignmentProducer.h.

Referenced by readInSurveyRcds().

edm::ESWatcher<TrackerSurveyRcd> AlignmentProducer::watchTkSurveyRcd_
private

Definition at line 194 of file AlignmentProducer.h.

Referenced by readInSurveyRcds().