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, boost::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 const *)
 
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, AlignmentErrors *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< CSCSurveyErrorRcdwatchCSCSurveyErrRcd_
 
edm::ESWatcher< CSCSurveyRcdwatchCSCSurveyRcd_
 
edm::ESWatcher< DTSurveyErrorRcdwatchDTSurveyErrRcd_
 
edm::ESWatcher< DTSurveyRcdwatchDTSurveyRcd_
 
edm::ESWatcher
< TrackerSurveyErrorRcd
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
ModuleChanger const * moduleChanger () const
 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 56 of file AlignmentProducer.h.

Member Typedef Documentation

Definition at line 60 of file AlignmentProducer.h.

Definition at line 61 of file AlignmentProducer.h.

Definition at line 62 of file AlignmentProducer.h.

Definition at line 64 of file AlignmentProducer.h.

Definition at line 65 of file AlignmentProducer.h.

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

Definition at line 66 of file AlignmentProducer.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 77 of file AlignmentProducer.cc.

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

77  :
81  nevent_(0), theParameterSet(iConfig),
82  theMaxLoops( iConfig.getUntrackedParameter<unsigned int>("maxLoops") ),
83  stNFixAlignables_(iConfig.getParameter<int>("nFixAlignables") ),
84  stRandomShift_(iConfig.getParameter<double>("randomShift")),
85  stRandomRotation_(iConfig.getParameter<double>("randomRotation")),
86  applyDbAlignment_( iConfig.getUntrackedParameter<bool>("applyDbAlignment")),
87  checkDbAlignmentValidity_( iConfig.getUntrackedParameter<bool>("checkDbAlignmentValidity")),
88  doMisalignmentScenario_(iConfig.getParameter<bool>("doMisalignmentScenario")),
89  saveToDB_(iConfig.getParameter<bool>("saveToDB")),
90  saveApeToDB_(iConfig.getParameter<bool>("saveApeToDB")),
91  saveDeformationsToDB_(iConfig.getParameter<bool>("saveDeformationsToDB")),
92  doTracker_( iConfig.getUntrackedParameter<bool>("doTracker") ),
93  doMuon_( iConfig.getUntrackedParameter<bool>("doMuon") ),
94  useExtras_( iConfig.getUntrackedParameter<bool>("useExtras") ),
95  useSurvey_( iConfig.getParameter<bool>("useSurvey") ),
96  tjTkAssociationMapTag_(iConfig.getParameter<edm::InputTag>("tjTkAssociationMapTag")),
97  beamSpotTag_(iConfig.getParameter<edm::InputTag>("beamSpotTag")),
98  tkLasBeamTag_(iConfig.getParameter<edm::InputTag>("tkLasBeamTag")),
99  clusterValueMapTag_(iConfig.getParameter<edm::InputTag>("hitPrescaleMapTag"))
100 {
101  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::AlignmentProducer";
102 
103  // Tell the framework what data is being produced
104  if (doTracker_) {
106  }
107  if (doMuon_) {
110  }
111 
112  // Create the alignment algorithm
113  edm::ParameterSet algoConfig = iConfig.getParameter<edm::ParameterSet>( "algoConfig" );
114  edm::VParameterSet iovSelection = iConfig.getParameter<edm::VParameterSet>( "RunRangeSelection" );
115  algoConfig.addUntrackedParameter<edm::VParameterSet>( "RunRangeSelection", iovSelection );
116  std::string algoName = algoConfig.getParameter<std::string>( "algoName" );
117  theAlignmentAlgo = AlignmentAlgorithmPluginFactory::get( )->create( algoName, algoConfig );
118 
119  // Check if found
120  if ( !theAlignmentAlgo )
121  throw cms::Exception("BadConfig") << "Couldn't find algorithm called " << algoName;
122 
123  // Now create monitors:
124  edm::ParameterSet monitorConfig = iConfig.getParameter<edm::ParameterSet>( "monitorConfig" );
125  std::vector<std::string> monitors = monitorConfig.getUntrackedParameter<std::vector<std::string> >( "monitors" );
126  for (std::vector<std::string>::const_iterator miter = monitors.begin(); miter != monitors.end(); ++miter) {
127  AlignmentMonitorBase* newMonitor = AlignmentMonitorPluginFactory::get()->create(*miter, monitorConfig.getUntrackedParameter<edm::ParameterSet>(*miter));
128 
129  if (!newMonitor) throw cms::Exception("BadConfig") << "Couldn't find monitor named " << *miter;
130 
131  theMonitors.push_back(newMonitor);
132  }
133 
134  // Finally create integrated calibrations:
135  edm::VParameterSet calibrations = iConfig.getParameter<edm::VParameterSet>("calibrations");
136  for (auto iCalib = calibrations.begin(); iCalib != calibrations.end(); ++iCalib) {
137  const std::string name(iCalib->getParameter<std::string>("calibrationName"));
139  // exception comes from line before: if (!theCalibrations.back()) throw cms::Exception(..) << ..;
140  }
141 
142 }
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:206
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.
SurfaceDeformation * create(int type, const std::vector< double > &params)
T get(const Candidate &c)
Definition: component.h:55
const double stRandomRotation_
AlignmentProducer::~AlignmentProducer ( )

Destructor.

Definition at line 147 of file AlignmentProducer.cc.

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

148 {
149  delete theAlignmentAlgo;
150 
151  // Delete monitors as well??
152 
153  for (auto iCal = theCalibrations.begin(); iCal != theCalibrations.end(); ++iCal) {
154  delete *iCal; // delete integrated calibration pointed to by (*iCal)
155  }
156 
158  delete theAlignableExtras;
159  delete theAlignableTracker;
160  delete theAlignableMuon;
161 
162  delete globalPositions_;
163 }
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 665 of file AlignmentProducer.cc.

References Alignable::alignableObjectId(), AlCaHLTBitMon_QueryRunRegistry::comp, Alignable::components(), relativeConstraints::error, edm::hlt::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().

666 {
667  const std::vector<Alignable*>& comp = ali->components();
668 
669  unsigned int nComp = comp.size();
670 
671  for (unsigned int i = 0; i < nComp; ++i) addSurveyInfo_(comp[i]);
672 
674 
675  if ( ali->id() != error.rawId() ||
676  ali->alignableObjectId() != error.structureType() )
677  {
678  throw cms::Exception("DatabaseError")
679  << "Error reading survey info from DB. Mismatched id!";
680  }
681 
682  const CLHEP::Hep3Vector& pos = theSurveyValues->m_align[theSurveyIndex].translation();
683  const CLHEP::HepRotation& rot = theSurveyValues->m_align[theSurveyIndex].rotation();
684 
685  AlignableSurface surf( align::PositionType( pos.x(), pos.y(), pos.z() ),
686  align::RotationType( rot.xx(), rot.xy(), rot.xz(),
687  rot.yx(), rot.yy(), rot.yz(),
688  rot.zx(), rot.zy(), rot.zz() ) );
689 
690  surf.setWidth( ali->surface().width() );
691  surf.setLength( ali->surface().length() );
692 
693  ali->setSurvey( new SurveyDet( surf, error.matrix() ) );
694 
695  ++theSurveyIndex;
696 }
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:180
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:126
const SurveyErrors * theSurveyErrors
align::Scalar length() const
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:268
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 764 of file AlignmentProducer.cc.

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

766 {
767  // 'G' is the geometry class for that DB should be applied,
768  // 'Rcd' is the record class for its Alignments
769  // 'ErrRcd' is the record class for its AlignmentErrors
770  // 'globalCoordinates' are global transformation for this geometry
771 
772  const Rcd & record = iSetup.get<Rcd>();
774  const edm::ValidityInterval & validity = record.validityInterval();
775  const edm::IOVSyncValue first = validity.first();
776  const edm::IOVSyncValue last = validity.last();
777  if (first!=edm::IOVSyncValue::beginOfTime() ||
779  throw cms::Exception("DatabaseError")
780  << "@SUB=AlignmentProducer::applyDB"
781  << "\nTrying to apply "
782  << record.key().name()
783  << " with multiple IOVs in tag.\n"
784  << "Validity range is "
785  << first.eventID().run() << " - " << last.eventID().run();
786  }
787  }
788 
789  edm::ESHandle<Alignments> alignments;
790  record.get(alignments);
791 
792  edm::ESHandle<AlignmentErrors> alignmentErrors;
793  iSetup.get<ErrRcd>().get(alignmentErrors);
794 
795  GeometryAligner aligner;
796  aligner.applyAlignments<G>(geometry, &(*alignments), &(*alignmentErrors),
797  globalCoordinates);
798 }
RunNumber_t run() const
Definition: EventID.h:42
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 applyAlignments(C *geometry, const Alignments *alignments, const AlignmentErrors *alignmentErrors, const AlignTransform &globalCoordinates)
bool first
Definition: L1TdeRCT.cc:75
tuple G
Definition: callgraph.py:12
const T & get() const
Definition: EventSetup.h:55
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 805 of file AlignmentProducer.cc.

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

806 {
807  // 'G' is the geometry class for that DB should be applied,
808  // 'DeformationRcd' is the record class for its surface deformations
809 
810  const DeformationRcd & record = iSetup.get<DeformationRcd>();
812  const edm::ValidityInterval & validity = record.validityInterval();
813  const edm::IOVSyncValue first = validity.first();
814  const edm::IOVSyncValue last = validity.last();
815  if (first!=edm::IOVSyncValue::beginOfTime() ||
817  throw cms::Exception("DatabaseError")
818  << "@SUB=AlignmentProducer::applyDB"
819  << "\nTrying to apply "
820  << record.key().name()
821  << " with multiple IOVs in tag.\n"
822  << "Validity range is "
823  << first.eventID().run() << " - " << last.eventID().run();
824  }
825  }
827  record.get(surfaceDeformations);
828 
829  GeometryAligner aligner;
830  aligner.attachSurfaceDeformations<G>(geometry, &(*surfaceDeformations));
831 }
RunNumber_t run() const
Definition: EventID.h:42
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)
bool first
Definition: L1TdeRCT.cc:75
tuple G
Definition: callgraph.py:12
const T & get() const
Definition: EventSetup.h:55
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 557 of file AlignmentProducer.cc.

References AlignmentAlgorithmBase::beginLuminosityBlock(), and theAlignmentAlgo.

559 {
560  theAlignmentAlgo->beginLuminosityBlock(setup); // do not forward edm::LuminosityBlock
561 }
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 196 of file AlignmentProducer.cc.

References AlignmentAlgorithmBase::addCalibrations(), align::AlignableMuon, AlignmentParameterStore_cfi::AlignmentParameterStore, applyDbAlignment_, MuonScenarioBuilder::applyScenario(), TrackerScenarioBuilder::applyScenario(), createGeometries_(), align::DetectorGlobalPosition(), doMisalignmentScenario_, doMuon_, doTracker_, edm::hlt::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_.

197 {
198  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob";
199 
200  //Retrieve tracker topology from geometry
201  edm::ESHandle<TrackerTopology> tTopoHandle;
202  iSetup.get<IdealGeometryRecord>().get(tTopoHandle);
203  const TrackerTopology* const tTopo = tTopoHandle.product();
204 
205  // Create the geometries from the ideal geometries (first time only)
206  this->createGeometries_( iSetup );
207 
208  // Retrieve and apply alignments, if requested (requires DB setup)
209  if ( applyDbAlignment_ ) {
210  // we need GlobalPositionRcd - and have to keep track for later removal
211  // before writing again to DB...
212  edm::ESHandle<Alignments> globalPositionRcd;
213  iSetup.get<GlobalPositionRcd>().get(globalPositionRcd);
214  globalPositions_ = new Alignments(*globalPositionRcd);
215 
216  if ( doTracker_ ) { // apply to tracker
217  this->applyDB<TrackerGeometry,TrackerAlignmentRcd,TrackerAlignmentErrorRcd>
218  (&(*theTracker), iSetup,
220  this->applyDB<TrackerGeometry,TrackerSurfaceDeformationRcd>(&(*theTracker), iSetup);
221  }
222 
223  if ( doMuon_ ) { // apply to tracker
224  this->applyDB<DTGeometry,DTAlignmentRcd,DTAlignmentErrorRcd>
225  (&(*theMuonDT), iSetup,
227  this->applyDB<CSCGeometry,CSCAlignmentRcd,CSCAlignmentErrorRcd>
228  (&(*theMuonCSC), iSetup,
230  }
231  }
232 
233  // Create alignable tracker and muon
234  if (doTracker_) {
235  theAlignableTracker = new AlignableTracker( &(*theTracker), tTopo );
236  }
237 
238  if (doMuon_) {
239  theAlignableMuon = new AlignableMuon( &(*theMuonDT), &(*theMuonCSC) );
240  }
241 
242  if (useExtras_) {
244  }
245 
246  // Create alignment parameter builder
247  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob"
248  << "Creating AlignmentParameterBuilder";
249  edm::ParameterSet aliParamBuildCfg =
250  theParameterSet.getParameter<edm::ParameterSet>("ParameterBuilder");
251  AlignmentParameterBuilder alignmentParameterBuilder(theAlignableTracker,
254  aliParamBuildCfg );
255  // Fix alignables if requested
256  if (stNFixAlignables_>0) alignmentParameterBuilder.fixAlignables(stNFixAlignables_);
257 
258  // Get list of alignables
259  Alignables theAlignables = alignmentParameterBuilder.alignables();
260  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob"
261  << "got " << theAlignables.size() << " alignables";
262 
263  // Create AlignmentParameterStore
264  edm::ParameterSet aliParamStoreCfg =
266  theAlignmentParameterStore = new AlignmentParameterStore(theAlignables, aliParamStoreCfg);
267  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob"
268  << "AlignmentParameterStore created!";
269 
270  // Apply misalignment scenario to alignable tracker and muon if requested
271  // WARNING: this assumes scenarioConfig can be passed to both muon and tracker
273  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob"
274  << "Applying misalignment scenario to "
275  << (doTracker_ ? "tracker" : "")
276  << (doMuon_ ? (doTracker_ ? " and muon" : "muon") : ".");
277  edm::ParameterSet scenarioConfig
278  = theParameterSet.getParameter<edm::ParameterSet>( "MisalignmentScenario" );
279  if (doTracker_) {
280  TrackerScenarioBuilder scenarioBuilder( theAlignableTracker );
281  scenarioBuilder.applyScenario( scenarioConfig );
282  }
283  if (doMuon_) {
284  MuonScenarioBuilder muonScenarioBuilder( theAlignableMuon );
285  muonScenarioBuilder.applyScenario( scenarioConfig );
286  }
287  } else {
288  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::beginOfJob"
289  << "NOT applying misalignment scenario!";
290  }
291 
292  // Apply simple misalignment
293  const std::string sParSel(theParameterSet.getParameter<std::string>("parameterSelectorSimple"));
294  this->simpleMisalignment_(theAlignables, sParSel, stRandomShift_, stRandomRotation_, true);
295 
296  // Initialize alignment algorithm and integrated calibration and pass the latter to algorithm
297  theAlignmentAlgo->initialize( iSetup,
299  theAlignmentParameterStore );
300  for (auto iCal = theCalibrations.begin(); iCal != theCalibrations.end(); ++iCal) {
302  }
303  // Not all algorithms support calibrations - so do not pass empty vector
304  // and throw if non-empty and not supported:
306  throw cms::Exception("BadConfig") << "[AlignmentProducer::beginOfJob]\n"
307  << "Configured " << theCalibrations.size() << " calibration(s) "
308  << "for algorithm not supporting it.";
309  }
310 
311  for (std::vector<AlignmentMonitorBase*>::const_iterator monitor = theMonitors.begin();
312  monitor != theMonitors.end(); ++monitor) {
313  (*monitor)->beginOfJob(theAlignableTracker, theAlignableMuon, theAlignmentParameterStore);
314  }
315 }
std::vector< Alignable * > Alignables
T getParameter(std::string const &) const
const bool doMisalignmentScenario_
Builds a scenario from configuration and applies it to the alignable Muon.
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.
virtual bool addCalibrations(const std::vector< IntegratedCalibrationBase * > &iCals)
const bool applyDbAlignment_
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:55
T const * product() const
Definition: ESHandle.h:62
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 532 of file AlignmentProducer.cc.

References AlignmentAlgorithmBase::beginRun(), and theAlignmentAlgo.

533 {
534  theAlignmentAlgo->beginRun(setup); // do not forward edm::Run...
535 }
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 641 of file AlignmentProducer.cc.

References TrackerGeomBuilderFromGeometricDet::build(), CSCGeometryBuilderFromDDD::build(), DTGeometryBuilderFromDDD::build(), doMuon_, doTracker_, edm::EventSetup::get(), theMuonCSC, theMuonDT, theParameterSet, and theTracker.

Referenced by beginOfJob().

642 {
644  iSetup.get<IdealGeometryRecord>().get( cpv );
645 
646  if (doTracker_) {
647  edm::ESHandle<GeometricDet> geometricDet;
648  iSetup.get<IdealGeometryRecord>().get( geometricDet );
649  TrackerGeomBuilderFromGeometricDet trackerBuilder;
650  theTracker = boost::shared_ptr<TrackerGeometry>( trackerBuilder.build(&(*geometricDet), theParameterSet ));
651  }
652 
653  if (doMuon_) {
655  iSetup.get<MuonNumberingRecord>().get(mdc);
656  DTGeometryBuilderFromDDD DTGeometryBuilder;
658  theMuonDT = boost::shared_ptr<DTGeometry>(new DTGeometry );
659  DTGeometryBuilder.build( theMuonDT, &(*cpv), *mdc);
660  theMuonCSC = boost::shared_ptr<CSCGeometry>( new CSCGeometry );
661  CSCGeometryBuilder.build( theMuonCSC, &(*cpv), *mdc );
662  }
663 }
boost::shared_ptr< TrackerGeometry > theTracker
TrackerGeometry * build(const GeometricDet *gd, const edm::ParameterSet &pSet)
edm::ParameterSet theParameterSet
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:55
edm::EDLooper::Status AlignmentProducer::duringLoop ( const edm::Event event,
const edm::EventSetup setup 
)
virtual

Called at each event.

Implements edm::EDLooper.

Definition at line 470 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_.

472 {
473  ++nevent_;
474 
475  // reading in survey records
476  this->readInSurveyRcds(setup);
477 
478  // Printout event number
479  for ( int i=10; i<10000000; i*=10 )
480  if ( nevent_<10*i && (nevent_%i)==0 )
481  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::duringLoop"
482  << "Events processed: " << nevent_;
483 
484  // Retrieve trajectories and tracks from the event
485  // -> merely skip if collection is empty
487  if (event.getByLabel(tjTkAssociationMapTag_, m_TrajTracksMap)) {
488 
489  // Form pairs of trajectories and tracks
490  ConstTrajTrackPairCollection trajTracks;
491  for ( TrajTrackAssociationCollection::const_iterator iPair = m_TrajTracksMap->begin();
492  iPair != m_TrajTracksMap->end(); ++iPair) {
493  trajTracks.push_back( ConstTrajTrackPair( &(*(*iPair).key), &(*(*iPair).val) ) );
494  }
496  event.getByLabel(beamSpotTag_, beamSpot);
497 
498  if (nevent_==1 && theAlignableExtras) {
499  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::duringLoop"
500  << "initializing AlignableBeamSpot" << std::endl;
501  theAlignableExtras->initializeBeamSpot(beamSpot->x0(), beamSpot->y0(), beamSpot->z0(),
502  beamSpot->dxdz(), beamSpot->dydz());
503  }
504 
505  // Run the alignment algorithm with its input
506  const AliClusterValueMap *clusterValueMapPtr = 0;
507  if(clusterValueMapTag_.encode().size()){//check that the input tag is not empty
508  edm::Handle<AliClusterValueMap> clusterValueMap;
509  event.getByLabel(clusterValueMapTag_, clusterValueMap);
510  clusterValueMapPtr = &(*clusterValueMap);
511  }
512 
513  const AlignmentAlgorithmBase::EventInfo eventInfo(event.id(), trajTracks, *beamSpot,
514  clusterValueMapPtr);
515  theAlignmentAlgo->run(setup, eventInfo);
516 
517 
518  for (std::vector<AlignmentMonitorBase*>::const_iterator monitor = theMonitors.begin();
519  monitor != theMonitors.end(); ++monitor) {
520  (*monitor)->duringLoop(event, setup, trajTracks); // forward eventInfo?
521  }
522  } else {
523  edm::LogError("Alignment") << "@SUB=AlignmentProducer::duringLoop"
524  << "No track collection found: skipping event";
525  }
526 
527 
528  return kContinue;
529 }
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
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:390
const edm::InputTag tjTkAssociationMapTag_
std::vector< ConstTrajTrackPair > ConstTrajTrackPairCollection
edm::EventID id() const
Definition: EventBase.h:56
std::pair< const Trajectory *, const reco::Track * > ConstTrajTrackPair
define event information passed to algorithms
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 564 of file AlignmentProducer.cc.

References AlignmentAlgorithmBase::endLuminosityBlock(), and theAlignmentAlgo.

566 {
567  theAlignmentAlgo->endLuminosityBlock(setup); // do not forward edm::LuminosityBlock
568 }
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 319 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().

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

Called at end of loop.

Implements edm::EDLooperBase.

Definition at line 435 of file AlignmentProducer.cc.

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

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

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

539 {
540  // call with or without las beam info...
541  typedef AlignmentAlgorithmBase::EndRunInfo EndRunInfo;
542  if (tkLasBeamTag_.encode().size()) { // non-empty InputTag
545  run.getByLabel(tkLasBeamTag_, lasBeams);
546  run.getByLabel(tkLasBeamTag_, tsoses);
547 
548  theAlignmentAlgo->endRun(EndRunInfo(run.id(), &(*lasBeams), &(*tsoses)), setup);
549  } else {
550  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::endRun"
551  << "No Tk LAS beams to forward to algorithm.";
552  theAlignmentAlgo->endRun(EndRunInfo(run.id(), 0, 0), setup);
553  }
554 }
bool getByLabel(std::string const &label, Handle< PROD > &result) const
Definition: Run.h:200
RunID const & id() const
Definition: RunBase.h:41
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.
define run information passed to algorithms (in endRun)
AlignmentProducer::RunRanges AlignmentProducer::makeNonOverlappingRunRanges ( const edm::VParameterSet RunRangeSelectionVPSet)
private

Definition at line 957 of file AlignmentProducer.cc.

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

Referenced by endOfJob().

958 {
959  static bool oldRunRangeSelectionWarning = false;
960 
963 
964  RunRanges uniqueRunRanges;
965  if (!RunRangeSelectionVPSet.empty()) {
966 
967  std::map<RunNumber,RunNumber> uniqueFirstRunNumbers;
968 
969  for (std::vector<edm::ParameterSet>::const_iterator ipset = RunRangeSelectionVPSet.begin();
970  ipset != RunRangeSelectionVPSet.end();
971  ++ipset) {
972  const std::vector<std::string> RunRangeStrings = (*ipset).getParameter<std::vector<std::string> >("RunRanges");
973  for (std::vector<std::string>::const_iterator irange = RunRangeStrings.begin();
974  irange != RunRangeStrings.end();
975  ++irange) {
976 
977  if ((*irange).find(':')==std::string::npos) {
978 
979  RunNumber first = beginValue;
980  long int temp = strtol((*irange).c_str(), 0, 0);
981  if (temp!=-1) first = temp;
982  uniqueFirstRunNumbers[first] = first;
983 
984  } else {
985 
986  if (!oldRunRangeSelectionWarning) {
987  edm::LogWarning("BadConfig") << "@SUB=AlignmentProducer::makeNonOverlappingRunRanges"
988  << "Config file contains old format for 'RunRangeSelection'. Only the start run\n"
989  << "number is used internally. The number of the last run is ignored and can be\n"
990  << "safely removed from the config file.\n";
991  oldRunRangeSelectionWarning = true;
992  }
993 
994  std::vector<std::string> tokens = edm::tokenize(*irange, ":");
995  long int temp;
996  RunNumber first = beginValue;
997  temp = strtol(tokens[0].c_str(), 0, 0);
998  if (temp!=-1) first = temp;
999  uniqueFirstRunNumbers[first] = first;
1000  }
1001  }
1002  }
1003 
1004  for (std::map<RunNumber,RunNumber>::iterator iFirst = uniqueFirstRunNumbers.begin();
1005  iFirst!=uniqueFirstRunNumbers.end();
1006  ++iFirst) {
1007  uniqueRunRanges.push_back(std::pair<RunNumber,RunNumber>((*iFirst).first, endValue));
1008  }
1009  for (unsigned int i = 0;i<uniqueRunRanges.size()-1;++i) {
1010  uniqueRunRanges[i].second = uniqueRunRanges[i+1].first - 1;
1011  }
1012 
1013  } else {
1014 
1015  uniqueRunRanges.push_back(std::pair<RunNumber,RunNumber>(beginValue, endValue));
1016 
1017  }
1018 
1019  return uniqueRunRanges;
1020 }
const TimeTypeSpecs timeTypeSpecs[]
Definition: Time.cc:22
int i
Definition: DBlmapReader.cc:9
Time_t beginValue
Definition: Time.h:45
AlignmentAlgorithmBase::RunNumber RunNumber
bool first
Definition: L1TdeRCT.cc:75
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
std::vector< RunRange > RunRanges
boost::shared_ptr< CSCGeometry > AlignmentProducer::produceCSC ( const MuonGeometryRecord iRecord)
virtual

Produce the muon CSC geometry.

Definition at line 187 of file AlignmentProducer.cc.

References theMuonCSC.

Referenced by AlignmentProducer().

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

Produce the muon DT geometry.

Definition at line 178 of file AlignmentProducer.cc.

References theMuonDT.

Referenced by AlignmentProducer().

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

Produce the tracker geometry.

Definition at line 169 of file AlignmentProducer.cc.

References theTracker.

Referenced by AlignmentProducer().

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

read in survey records

Definition at line 698 of file AlignmentProducer.cc.

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

Referenced by duringLoop().

698  {
699 
700  // Get Survey Rcds and add Survey Info
701  if ( doTracker_ && useSurvey_ ){
702  bool tkSurveyBool = watchTkSurveyRcd_.check(iSetup);
703  bool tkSurveyErrBool = watchTkSurveyErrRcd_.check(iSetup);
704  edm::LogInfo("Alignment") << "watcher tksurveyrcd: " << tkSurveyBool;
705  edm::LogInfo("Alignment") << "watcher tksurveyerrrcd: " << tkSurveyErrBool;
706  if ( tkSurveyBool || tkSurveyErrBool){
707 
708  edm::LogInfo("Alignment") << "ADDING THE SURVEY INFORMATION";
710  edm::ESHandle<SurveyErrors> surveyErrors;
711 
712  iSetup.get<TrackerSurveyRcd>().get(surveys);
713  iSetup.get<TrackerSurveyErrorRcd>().get(surveyErrors);
714 
715  theSurveyIndex = 0;
716  theSurveyValues = &*surveys;
717  theSurveyErrors = &*surveyErrors;
719  }
720  }
721 
722  if ( doMuon_ && useSurvey_) {
723  bool DTSurveyBool = watchTkSurveyRcd_.check(iSetup);
724  bool DTSurveyErrBool = watchTkSurveyErrRcd_.check(iSetup);
725  bool CSCSurveyBool = watchTkSurveyRcd_.check(iSetup);
726  bool CSCSurveyErrBool = watchTkSurveyErrRcd_.check(iSetup);
727 
728  if ( DTSurveyBool || DTSurveyErrBool || CSCSurveyBool || CSCSurveyErrBool ){
729  edm::ESHandle<Alignments> dtSurveys;
730  edm::ESHandle<SurveyErrors> dtSurveyErrors;
731  edm::ESHandle<Alignments> cscSurveys;
732  edm::ESHandle<SurveyErrors> cscSurveyErrors;
733 
734  iSetup.get<DTSurveyRcd>().get(dtSurveys);
735  iSetup.get<DTSurveyErrorRcd>().get(dtSurveyErrors);
736  iSetup.get<CSCSurveyRcd>().get(cscSurveys);
737  iSetup.get<CSCSurveyErrorRcd>().get(cscSurveyErrors);
738 
739  theSurveyIndex = 0;
740  theSurveyValues = &*dtSurveys;
741  theSurveyErrors = &*dtSurveyErrors;
742  std::vector<Alignable*> barrels = theAlignableMuon->DTBarrel();
743  for (std::vector<Alignable*>::const_iterator iter = barrels.begin(); iter != barrels.end(); ++iter) {
745  }
746 
747  theSurveyIndex = 0;
748  theSurveyValues = &*cscSurveys;
749  theSurveyErrors = &*cscSurveyErrors;
750  std::vector<Alignable*> endcaps = theAlignableMuon->CSCEndcaps();
751  for (std::vector<Alignable*>::const_iterator iter = endcaps.begin(); iter != endcaps.end(); ++iter) {
753  }
754  }
755  }
756 
757 }
align::Alignables DTBarrel()
const Alignments * theSurveyValues
void addSurveyInfo_(Alignable *)
Add survey info to an alignable.
align::Alignables CSCEndcaps()
unsigned int theSurveyIndex
const SurveyErrors * theSurveyErrors
edm::ESWatcher< TrackerSurveyErrorRcd > watchTkSurveyErrRcd_
AlignableMuon * theAlignableMuon
const T & get() const
Definition: EventSetup.h:55
bool check(const edm::EventSetup &iSetup)
Definition: ESWatcher.h:58
edm::ESWatcher< TrackerSurveyRcd > watchTkSurveyRcd_
AlignableTracker * theAlignableTracker
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 572 of file AlignmentProducer.cc.

References funct::abs(), Alignable::alignmentParameters(), AlignmentParameterSelector::convertParamSel(), RigidBodyAlignmentParameters::dalpha, RigidBodyAlignmentParameters::dbeta, RigidBodyAlignmentParameters::dgamma, RigidBodyAlignmentParameters::dx, RigidBodyAlignmentParameters::dy, RigidBodyAlignmentParameters::dz, edm::hlt::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().

574 {
575 
576  std::ostringstream output; // collecting output
577 
578  if (shift > 0. || rot > 0.) {
579  output << "Adding random flat shift of max size " << shift
580  << " and adding random flat rotation of max size " << rot <<" to ";
581 
582  std::vector<bool> commSel(0);
583  if (selection != "-1") {
584  AlignmentParameterSelector aSelector(0,0); // no alignable needed here...
585  const std::vector<char> cSel(aSelector.convertParamSel(selection));
586  if (cSel.size() < RigidBodyAlignmentParameters::N_PARAM) {
587  throw cms::Exception("BadConfig")
588  << "[AlignmentProducer::simpleMisalignment_]\n"
589  << "Expect selection string '" << selection << "' to be at least of length "
590  << RigidBodyAlignmentParameters::N_PARAM << " or to be '-1'.\n"
591  << "(Most probably you have to adjust the parameter 'parameterSelectorSimple'.)";
592  }
593  for (std::vector<char>::const_iterator cIter = cSel.begin(); cIter != cSel.end(); ++cIter) {
594  commSel.push_back(*cIter == '0' ? false : true);
595  }
596  output << "parameters defined by (" << selection
597  << "), representing (x,y,z,alpha,beta,gamma),";
598  } else {
599  output << "the active parameters of each alignable,";
600  }
601  output << " in " << (local ? "local" : "global") << " frame.";
602 
603  for (std::vector<Alignable*>::const_iterator it = alivec.begin(); it != alivec.end(); ++it) {
604  Alignable* ali=(*it);
605  std::vector<bool> mysel(commSel.empty() ? ali->alignmentParameters()->selector() : commSel);
606 
607  if (std::abs(shift)>0.00001) {
608  double s0 = 0., s1 = 0., s2 = 0.;
609  if (mysel[RigidBodyAlignmentParameters::dx]) s0 = shift * double(random()%1000-500)/500.;
610  if (mysel[RigidBodyAlignmentParameters::dy]) s1 = shift * double(random()%1000-500)/500.;
611  if (mysel[RigidBodyAlignmentParameters::dz]) s2 = shift * double(random()%1000-500)/500.;
612 
613  if (local) ali->move( ali->surface().toGlobal(align::LocalVector(s0,s1,s2)) );
614  else ali->move( align::GlobalVector(s0,s1,s2) );
615 
616  //AlignmentPositionError ape(dx,dy,dz);
617  //ali->addAlignmentPositionError(ape);
618  }
619 
620  if (std::abs(rot)>0.00001) {
622  if (mysel[RigidBodyAlignmentParameters::dalpha]) r(1)=rot*double(random()%1000-500)/500.;
623  if (mysel[RigidBodyAlignmentParameters::dbeta]) r(2)=rot*double(random()%1000-500)/500.;
624  if (mysel[RigidBodyAlignmentParameters::dgamma]) r(3)=rot*double(random()%1000-500)/500.;
625 
626  const align::RotationType mrot = align::toMatrix(r);
627  if (local) ali->rotateInLocalFrame(mrot);
628  else ali->rotateInGlobalFrame(mrot);
629 
630  //ali->addAlignmentPositionErrorFromRotation(mrot);
631  }
632  } // end loop on alignables
633  } else {
634  output << "No simple misalignment added!";
635  }
636  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::simpleMisalignment_" << output.str();
637 }
virtual void rotateInLocalFrame(const RotationType &rotation)
Rotation intepreted in the local reference frame.
Definition: Alignable.cc:91
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:126
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 389 of file AlignmentProducer.cc.

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

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

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

Referenced by writeForRunRange().

885 {
886  Alignments * tempAlignments = alignments;
887  AlignmentErrors * tempAlignmentErrors = alignmentErrors;
888 
889  // Call service
891  if (!poolDb.isAvailable()) { // Die if not available
892  delete tempAlignments; // promised to take over ownership...
893  delete tempAlignmentErrors; // dito
894  throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
895  }
896 
897  if (globalCoordinates // happens only if (applyDbAlignment_ == true)
898  && globalCoordinates->transform() != AlignTransform::Transform::Identity) {
899 
900  tempAlignments = new Alignments(); // temporary storage for
901  tempAlignmentErrors = new AlignmentErrors(); // final alignments and errors
902 
903  GeometryAligner aligner;
904  aligner.removeGlobalTransform(alignments, alignmentErrors,
905  *globalCoordinates,
906  tempAlignments, tempAlignmentErrors);
907 
908  delete alignments; // have to delete original alignments
909  delete alignmentErrors; // same thing for the errors
910 
911  edm::LogInfo("Alignment") << "@SUB=AlignmentProducer::writeDB"
912  << "globalCoordinates removed from alignments (" << alignRcd
913  << ") and errors (" << alignRcd << ").";
914  }
915 
916  if (saveToDB_) {
917  edm::LogInfo("Alignment") << "Writing Alignments for run " << time
918  << " to " << alignRcd << ".";
919  poolDb->writeOne<Alignments>(tempAlignments, time, alignRcd);
920  } else { // poolDb->writeOne(..) takes over 'alignments' ownership,...
921  delete tempAlignments; // ...otherwise we have to delete, as promised!
922  }
923 
924  if (saveApeToDB_) {
925  edm::LogInfo("Alignment") << "Writing AlignmentErrors for run " << time
926  << " to " << errRcd << ".";
927  poolDb->writeOne<AlignmentErrors>(tempAlignmentErrors, time, errRcd);
928  } else { // poolDb->writeOne(..) takes over 'alignmentErrors' ownership,...
929  delete tempAlignmentErrors; // ...otherwise we have to delete, as promised!
930  }
931 }
Class to update a given geometry with a set of alignments.
void removeGlobalTransform(const Alignments *alignments, const AlignmentErrors *alignmentErrors, const AlignTransform &globalCoordinates, Alignments *newAlignments, AlignmentErrors *newAlignmentErrors)
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 935 of file AlignmentProducer.cc.

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

938 {
939  // Call service
941  if (!poolDb.isAvailable()) { // Die if not available
942  delete alignmentSurfaceDeformations; // promised to take over ownership...
943  throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
944  }
945 
946  if (saveDeformationsToDB_) {
947  edm::LogInfo("Alignment") << "Writing AlignmentSurfaceDeformations for run " << time
948  << " to " << surfaceDeformationRcd << ".";
949  poolDb->writeOne<AlignmentSurfaceDeformations>(alignmentSurfaceDeformations, time,
950  surfaceDeformationRcd);
951  } else { // poolDb->writeOne(..) takes over 'surfaceDeformation' ownership,...
952  delete alignmentSurfaceDeformations; // ...otherwise we have to delete, as promised!
953  }
954 }
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 834 of file AlignmentProducer.cc.

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

Referenced by endOfJob().

835 {
836  if ( doTracker_ ) { // first tracker
837  const AlignTransform *trackerGlobal = 0; // will be 'removed' from constants
838  if (globalPositions_) { // i.e. applied before in applyDB
841  }
842 
843  Alignments *alignments = theAlignableTracker->alignments();
845  this->writeDB(alignments, "TrackerAlignmentRcd",
846  alignmentErrors, "TrackerAlignmentErrorRcd", trackerGlobal,
847  time);
848  }
849 
850  if ( doMuon_ ) { // now muon
851  const AlignTransform *muonGlobal = 0; // will be 'removed' from constants
852  if (globalPositions_) { // i.e. applied before in applyDB
854  DetId(DetId::Muon));
855  }
856  // Get alignments+errors, first DT - ownership taken over by writeDB(..), so no delete
857  Alignments *alignments = theAlignableMuon->dtAlignments();
858  AlignmentErrors *alignmentErrors = theAlignableMuon->dtAlignmentErrors();
859  this->writeDB(alignments, "DTAlignmentRcd",
860  alignmentErrors, "DTAlignmentErrorRcd", muonGlobal,
861  time);
862 
863  // Get alignments+errors, now CSC - ownership taken over by writeDB(..), so no delete
864  alignments = theAlignableMuon->cscAlignments();
865  alignmentErrors = theAlignableMuon->cscAlignmentErrors();
866  this->writeDB(alignments, "CSCAlignmentRcd",
867  alignmentErrors, "CSCAlignmentErrorRcd", muonGlobal,
868  time);
869  }
870 
871  // Save surface deformations to database
873  AlignmentSurfaceDeformations *alignmentSurfaceDeformations = theAlignableTracker->surfaceDeformations();
874  this->writeDB(alignmentSurfaceDeformations, "TrackerSurfaceDeformationRcd", time);
875  }
876 }
AlignmentErrors * dtAlignmentErrors()
Alignments * dtAlignments()
AlignmentErrors * alignmentErrors() const
Return alignment errors, sorted by DetId.
AlignmentErrors * cscAlignmentErrors()
void writeDB(Alignments *alignments, const std::string &alignRcd, AlignmentErrors *alignmentErrors, const std::string &errRcd, const AlignTransform *globalCoordinates, cond::Time_t time) const
AlignmentSurfaceDeformations * surfaceDeformations() const
Return surface deformations, sorted by DetId.
Definition: Alignable.cc:216
Definition: DetId.h:18
AlignableMuon * theAlignableMuon
const bool saveDeformationsToDB_
Alignments * cscAlignments()
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 182 of file AlignmentProducer.h.

Referenced by beginOfJob().

const edm::InputTag AlignmentProducer::beamSpotTag_
private

Definition at line 190 of file AlignmentProducer.h.

Referenced by duringLoop().

const bool AlignmentProducer::checkDbAlignmentValidity_
private

Definition at line 182 of file AlignmentProducer.h.

Referenced by applyDB().

const edm::InputTag AlignmentProducer::clusterValueMapTag_
private

Definition at line 192 of file AlignmentProducer.h.

Referenced by duringLoop().

const bool AlignmentProducer::doMisalignmentScenario_
private

Definition at line 183 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 171 of file AlignmentProducer.h.

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

int AlignmentProducer::nevent_
private

Definition at line 173 of file AlignmentProducer.h.

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

const bool AlignmentProducer::saveApeToDB_
private

Definition at line 184 of file AlignmentProducer.h.

Referenced by endOfJob(), and writeDB().

const bool AlignmentProducer::saveDeformationsToDB_
private

Definition at line 184 of file AlignmentProducer.h.

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

const bool AlignmentProducer::saveToDB_
private

Definition at line 184 of file AlignmentProducer.h.

Referenced by endOfJob(), and writeDB().

const int AlignmentProducer::stNFixAlignables_
private

Definition at line 180 of file AlignmentProducer.h.

Referenced by beginOfJob().

const double AlignmentProducer::stRandomRotation_
private

Definition at line 181 of file AlignmentProducer.h.

Referenced by beginOfJob().

const double AlignmentProducer::stRandomShift_
private

Definition at line 181 of file AlignmentProducer.h.

Referenced by beginOfJob().

AlignableExtras* AlignmentProducer::theAlignableExtras
private

Definition at line 163 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 159 of file AlignmentProducer.h.

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

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

Definition at line 178 of file AlignmentProducer.h.

Referenced by endOfLoop().

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

Definition at line 169 of file AlignmentProducer.h.

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

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

Definition at line 168 of file AlignmentProducer.h.

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

edm::ParameterSet AlignmentProducer::theParameterSet
private

Definition at line 174 of file AlignmentProducer.h.

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

const SurveyErrors* AlignmentProducer::theSurveyErrors
private

Definition at line 156 of file AlignmentProducer.h.

Referenced by addSurveyInfo_(), and readInSurveyRcds().

unsigned int AlignmentProducer::theSurveyIndex
private

Definition at line 154 of file AlignmentProducer.h.

Referenced by addSurveyInfo_(), and readInSurveyRcds().

const Alignments* AlignmentProducer::theSurveyValues
private

Definition at line 155 of file AlignmentProducer.h.

Referenced by addSurveyInfo_(), and readInSurveyRcds().

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

Definition at line 167 of file AlignmentProducer.h.

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

const edm::InputTag AlignmentProducer::tjTkAssociationMapTag_
private

Definition at line 189 of file AlignmentProducer.h.

Referenced by duringLoop().

const edm::InputTag AlignmentProducer::tkLasBeamTag_
private

Definition at line 191 of file AlignmentProducer.h.

Referenced by endRun().

const bool AlignmentProducer::useExtras_
private

Definition at line 185 of file AlignmentProducer.h.

Referenced by beginOfJob().

const bool AlignmentProducer::useSurvey_
private

Definition at line 186 of file AlignmentProducer.h.

Referenced by readInSurveyRcds().

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

Definition at line 200 of file AlignmentProducer.h.

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

Definition at line 199 of file AlignmentProducer.h.

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

Definition at line 198 of file AlignmentProducer.h.

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

Definition at line 197 of file AlignmentProducer.h.

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

Definition at line 196 of file AlignmentProducer.h.

Referenced by readInSurveyRcds().

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

Definition at line 195 of file AlignmentProducer.h.

Referenced by readInSurveyRcds().