Inheritance diagram for ReferenceTrajectoryFactory:

Public Member Functions
ReferenceTrajectoryFactory *	clone () const override

	ReferenceTrajectoryFactory (const edm::ParameterSet &config)

const ReferenceTrajectoryCollection	trajectories (const edm::EventSetup &setup, const ConstTrajTrackPairCollection &tracks, const reco::BeamSpot &beamSpot) const override
	Produce the reference trajectories. More...

const ReferenceTrajectoryCollection	trajectories (const edm::EventSetup &setup, const ConstTrajTrackPairCollection &tracks, const ExternalPredictionCollection &external, const reco::BeamSpot &beamSpot) const override

	~ReferenceTrajectoryFactory () override

Public Member Functions inherited from TrajectoryFactoryBase
const edm::ParameterSet &	configuration () const

MaterialEffects	materialEffects (void) const

PropagationDirection	propagationDirection (void) const

unsigned int	tracksPerTrajectory () const

	TrajectoryFactoryBase (const edm::ParameterSet &config)

	TrajectoryFactoryBase (const edm::ParameterSet &config, unsigned int tracksPerTrajectory)

virtual	~TrajectoryFactoryBase (void)

Protected Member Functions
const TrajectoryFactoryBase *	bzeroFactory () const

	ReferenceTrajectoryFactory (const ReferenceTrajectoryFactory &other)

Protected Member Functions inherited from TrajectoryFactoryBase
virtual const TrajectoryInput	innermostStateAndRecHits (const ConstTrajTrackPair &track) const

virtual const Trajectory::DataContainer	orderedTrajectoryMeasurements (const Trajectory &trajectory) const

bool	sameSurface (const Surface &s1, const Surface &s2) const

bool	useRecHit (const TransientTrackingRecHit::ConstRecHitPointer &hitPtr) const

Protected Attributes
const TrajectoryFactoryBase *	theBzeroFactory

double	theMass

bool	theUseBzeroIfFieldOff

Protected Attributes inherited from TrajectoryFactoryBase
const bool	allowZeroMaterial_

const bool	includeAPEs_

const bool	useBeamSpot_

Additional Inherited Members
Public Types inherited from TrajectoryFactoryBase
typedef AlignmentAlgorithmBase::ConstTrajTrackPair	ConstTrajTrackPair

typedef AlignmentAlgorithmBase::ConstTrajTrackPairCollection	ConstTrajTrackPairCollection

typedef std::vector< TrajectoryStateOnSurface >	ExternalPredictionCollection

typedef ReferenceTrajectoryBase::MaterialEffects	MaterialEffects

typedef std::vector< ReferenceTrajectoryPtr >	ReferenceTrajectoryCollection

typedef ReferenceTrajectoryBase::ReferenceTrajectoryPtr	ReferenceTrajectoryPtr

typedef std::pair< TrajectoryStateOnSurface, TransientTrackingRecHit::ConstRecHitContainer >	TrajectoryInput

Detailed Description

A factory that produces instances of class ReferenceTrajectory from a given TrajTrackPairCollection. If |B| = 0 T and configuration parameter UseBzeroIfFieldOff is True, hand-over to the BzeroReferenceTrajectoryFactory.

Definition at line 20 of file ReferenceTrajectoryFactory.cc.

Constructor & Destructor Documentation

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

Definition at line 50 of file ReferenceTrajectoryFactory.cc.

References theMass, and theUseBzeroIfFieldOff.

Referenced by clone().

     : TrajectoryFactoryBase(config),
       theMass(config.getParameter<double>("ParticleMass")),
       theUseBzeroIfFieldOff(config.getParameter<bool>("UseBzeroIfFieldOff")),
       theBzeroFactory(nullptr) {
   edm::LogInfo("Alignment") << "@SUB=ReferenceTrajectoryFactory"
                             << "mass: " << theMass
                             << "\nusing Bzero if |B| = 0: " << (theUseBzeroIfFieldOff ? "yes" : "no");
 }

ReferenceTrajectoryFactory::~ReferenceTrajectoryFactory ( void )

override

Definition at line 67 of file ReferenceTrajectoryFactory.cc.

References theBzeroFactory.

67 { delete theBzeroFactory; }

ReferenceTrajectoryFactory::theBzeroFactory

const TrajectoryFactoryBase * theBzeroFactory

Definition: ReferenceTrajectoryFactory.cc:43

ReferenceTrajectoryFactory::ReferenceTrajectoryFactory ( const ReferenceTrajectoryFactory & other )

protected

Definition at line 60 of file ReferenceTrajectoryFactory.cc.

     : TrajectoryFactoryBase(other),
       theMass(other.theMass),
       theUseBzeroIfFieldOff(other.theUseBzeroIfFieldOff),
       theBzeroFactory(nullptr)  // copy data members, but no double pointing to same Bzero factory...
 {}

Member Function Documentation

const TrajectoryFactoryBase * ReferenceTrajectoryFactory::bzeroFactory ( ) const

protected

Definition at line 161 of file ReferenceTrajectoryFactory.cc.

References edm::ParameterSet::addParameter(), TrajectoryFactories_cff::BzeroReferenceTrajectoryFactory, TrajectoryFactoryBase::configuration(), edm::ParameterSet::copyForModify(), DEFINE_EDM_PLUGIN, edm::ParameterSet::eraseSimpleParameter(), edm::ParameterSet::getParameter(), muonDTDigis_cfi::pset, AlCaHLTBitMon_QueryRunRegistry::string, and theBzeroFactory.

Referenced by clone(), and trajectories().

                                                                             {
   if (!theBzeroFactory) {
     const edm::ParameterSet &myPset = this->configuration();
     edm::LogInfo("Alignment") << "@SUB=ReferenceTrajectoryFactory::bzeroFactory"
                               << "Using BzeroReferenceTrajectoryFactory for some (all?) events.";
     // We take the config of this factory, copy it, replace its name and add
     // the momentum parameter as expected by BzeroReferenceTrajectoryFactory and create it:
     edm::ParameterSet pset;
     pset.copyForModify(myPset);
     // next two lines not needed, but may help to better understand log file:
     pset.eraseSimpleParameter("TrajectoryFactoryName");
     pset.addParameter("TrajectoryFactoryName", std::string("BzeroReferenceTrajectoryFactory"));
     pset.addParameter("MomentumEstimate", myPset.getParameter<double>("MomentumEstimateFieldOff"));
     theBzeroFactory = new BzeroReferenceTrajectoryFactory(pset);
   }
   return theBzeroFactory;
 }

ReferenceTrajectoryFactory* ReferenceTrajectoryFactory::clone ( void ) const

inlineoverridevirtual

Implements TrajectoryFactoryBase.

Definition at line 35 of file ReferenceTrajectoryFactory.cc.

References bzeroFactory(), trackingPlots::other, and ReferenceTrajectoryFactory().

35 { return new ReferenceTrajectoryFactory(*this); }

ReferenceTrajectoryFactory::ReferenceTrajectoryFactory

ReferenceTrajectoryFactory(const edm::ParameterSet &config)

Definition: ReferenceTrajectoryFactory.cc:50

const ReferenceTrajectoryFactory::ReferenceTrajectoryCollection ReferenceTrajectoryFactory::trajectories	(	const edm::EventSetup &	setup,
		const ConstTrajTrackPairCollection &	tracks,
		const reco::BeamSpot &	beamSpot
	)		const

overridevirtual

Produce the reference trajectories.

Implements TrajectoryFactoryBase.

Definition at line 69 of file ReferenceTrajectoryFactory.cc.

References ReferenceTrajectoryBase::Config::allowZeroMaterial, TrajectoryFactoryBase::allowZeroMaterial_, pwdgSkimBPark_cfi::beamSpot, bzeroFactory(), looper::config, edm::EventSetup::get(), ReferenceTrajectoryBase::Config::hitsAreReverse, ReferenceTrajectoryBase::Config::includeAPEs, TrajectoryFactoryBase::includeAPEs_, TrajectoryFactoryBase::innermostStateAndRecHits(), input, MagneticField::inTesla(), PV3DBase< T, PVType, FrameType >::mag2(), HLT_2018_cff::magneticField, TrajectoryFactoryBase::materialEffects(), edm::ESHandle< T >::product(), TrajectoryFactoryBase::propagationDirection(), theMass, theUseBzeroIfFieldOff, TrajectoryFactoryBase::trajectories(), ReferenceTrajectoryBase::Config::useBeamSpot, and TrajectoryFactoryBase::useBeamSpot_.

Referenced by trajectories().

                                                                                                                 {
   edm::ESHandle<MagneticField> magneticField;
   setup.get<IdealMagneticFieldRecord>().get(magneticField);
   if (theUseBzeroIfFieldOff && magneticField->inTesla(GlobalPoint(0., 0., 0.)).mag2() < 1.e-6) {
     return this->bzeroFactory()->trajectories(setup, tracks, beamSpot);
   }
 
   ReferenceTrajectoryCollection trajectories;
 
   ConstTrajTrackPairCollection::const_iterator itTracks = tracks.begin();
 
   while (itTracks != tracks.end()) {
     TrajectoryInput input = this->innermostStateAndRecHits(*itTracks);
 
     // Check input: If all hits were rejected, the TSOS is initialized as invalid.
     if (input.first.isValid()) {
       ReferenceTrajectoryBase::Config config(materialEffects(), propagationDirection(), theMass);
       config.useBeamSpot = useBeamSpot_;
       config.includeAPEs = includeAPEs_;
       config.allowZeroMaterial = allowZeroMaterial_;
       // set the flag for reversing the RecHits to false, since they are already in the correct order.
       config.hitsAreReverse = false;
       trajectories.push_back(ReferenceTrajectoryPtr(
           new ReferenceTrajectory(input.first, input.second, magneticField.product(), beamSpot, config)));
     }
 
     ++itTracks;
   }
 
   return trajectories;
 }

const ReferenceTrajectoryFactory::ReferenceTrajectoryCollection ReferenceTrajectoryFactory::trajectories	(	const edm::EventSetup &	setup,
		const ConstTrajTrackPairCollection &	tracks,
		const ExternalPredictionCollection &	external,
		const reco::BeamSpot &	beamSpot
	)		const

overridevirtual

Implements TrajectoryFactoryBase.

Definition at line 102 of file ReferenceTrajectoryFactory.cc.

                                         {
   ReferenceTrajectoryCollection trajectories;
 
   if (tracks.size() != external.size()) {
     edm::LogInfo("ReferenceTrajectories")
         << "@SUB=ReferenceTrajectoryFactory::trajectories"
         << "Inconsistent input:\n"
         << "\tnumber of tracks = " << tracks.size() << "\tnumber of external predictions = " << external.size();
     return trajectories;
   }
 
   edm::ESHandle<MagneticField> magneticField;
   setup.get<IdealMagneticFieldRecord>().get(magneticField);
   if (theUseBzeroIfFieldOff && magneticField->inTesla(GlobalPoint(0., 0., 0.)).mag2() < 1.e-6) {
     return this->bzeroFactory()->trajectories(setup, tracks, external, beamSpot);
   }
 
   ConstTrajTrackPairCollection::const_iterator itTracks = tracks.begin();
   ExternalPredictionCollection::const_iterator itExternal = external.begin();
 
   while (itTracks != tracks.end()) {
     TrajectoryInput input = innermostStateAndRecHits(*itTracks);
     // Check input: If all hits were rejected, the TSOS is initialized as invalid.
     if (input.first.isValid()) {
       if ((*itExternal).isValid() && sameSurface((*itExternal).surface(), input.first.surface())) {
         ReferenceTrajectoryBase::Config config(materialEffects(), propagationDirection(), theMass);
         config.useBeamSpot = useBeamSpot_;
         config.includeAPEs = includeAPEs_;
         config.allowZeroMaterial = allowZeroMaterial_;
         // set the flag for reversing the RecHits to false, since they are already in the correct order.
         config.hitsAreReverse = false;
         ReferenceTrajectoryPtr refTraj(
             new ReferenceTrajectory(*itExternal, input.second, magneticField.product(), beamSpot, config));
 
         AlgebraicSymMatrix externalParamErrors(asHepMatrix<5>((*itExternal).localError().matrix()));
         refTraj->setParameterErrors(externalParamErrors);
         trajectories.push_back(refTraj);
       } else {
         ReferenceTrajectoryBase::Config config(materialEffects(), propagationDirection(), theMass);
         config.useBeamSpot = useBeamSpot_;
         config.includeAPEs = includeAPEs_;
         config.allowZeroMaterial = allowZeroMaterial_;
         config.hitsAreReverse = false;
         trajectories.push_back(ReferenceTrajectoryPtr(
             new ReferenceTrajectory(input.first, input.second, magneticField.product(), beamSpot, config)));
       }
     }
 
     ++itTracks;
     ++itExternal;
   }
 
   return trajectories;
 }