Inheritance diagram for TwoBodyDecayTrajectoryFactory:

Public Types
typedef TwoBodyDecayTrajectory::ConstRecHitCollection	ConstRecHitCollection

typedef TwoBodyDecayTrajectoryState::TsosContainer	TsosContainer

typedef TwoBodyDecayVirtualMeasurement	VirtualMeasurement

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

Public Member Functions
TwoBodyDecayTrajectoryFactory *	clone () const override

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

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

	TwoBodyDecayTrajectoryFactory (const edm::ParameterSet &config)

	~TwoBodyDecayTrajectoryFactory () 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 ReferenceTrajectoryCollection	constructTrajectories (const ConstTrajTrackPairCollection &tracks, const TwoBodyDecay &tbd, const MagneticField *magField, const reco::BeamSpot &beamSpot, bool setParameterErrors) const

bool	match (const TrajectoryStateOnSurface &state, const TransientTrackingRecHit::ConstRecHitPointer &recHit) const

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
double	theChi2CutValue

bool	theConstructTsosWithErrorsFlag

TwoBodyDecayFitter	theFitter

double	theNSigmaCutValue

double	thePrimaryMass

double	thePrimaryWidth

double	theSecondaryMass

bool	theUseRefittedStateFlag

Protected Attributes inherited from TrajectoryFactoryBase
const bool	allowZeroMaterial_

const bool	includeAPEs_

const bool	useBeamSpot_

Detailed Description

by Edmund Widl, see CMS NOTE-2007/032. extension for BreakPoints or BrokenLines by Claus Kleinwort

Definition at line 28 of file TwoBodyDecayTrajectoryFactory.cc.

Member Typedef Documentation

typedef TwoBodyDecayTrajectory::ConstRecHitCollection TwoBodyDecayTrajectoryFactory::ConstRecHitCollection

Definition at line 32 of file TwoBodyDecayTrajectoryFactory.cc.

typedef TwoBodyDecayTrajectoryState::TsosContainer TwoBodyDecayTrajectoryFactory::TsosContainer

Definition at line 31 of file TwoBodyDecayTrajectoryFactory.cc.

typedef TwoBodyDecayVirtualMeasurement TwoBodyDecayTrajectoryFactory::VirtualMeasurement

Definition at line 30 of file TwoBodyDecayTrajectoryFactory.cc.

Constructor & Destructor Documentation

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

Definition at line 74 of file TwoBodyDecayTrajectoryFactory.cc.

References edm::ParameterSet::getParameter(), theChi2CutValue, theConstructTsosWithErrorsFlag, theNSigmaCutValue, thePrimaryMass, thePrimaryWidth, theSecondaryMass, and theUseRefittedStateFlag.

Referenced by clone().

     : TrajectoryFactoryBase(config, 2), theFitter(config) {
   const edm::ParameterSet ppc = config.getParameter<edm::ParameterSet>("ParticleProperties");
   thePrimaryMass = ppc.getParameter<double>("PrimaryMass");
   thePrimaryWidth = ppc.getParameter<double>("PrimaryWidth");
   theSecondaryMass = ppc.getParameter<double>("SecondaryMass");
 
   theNSigmaCutValue = config.getParameter<double>("NSigmaCut");
   theChi2CutValue = config.getParameter<double>("Chi2Cut");
   theUseRefittedStateFlag = config.getParameter<bool>("UseRefittedState");
   theConstructTsosWithErrorsFlag = config.getParameter<bool>("ConstructTsosWithErrors");
 }

TwoBodyDecayTrajectoryFactory::~TwoBodyDecayTrajectoryFactory ( )

override

Definition at line 87 of file TwoBodyDecayTrajectoryFactory.cc.

87 {}

Member Function Documentation

TwoBodyDecayTrajectoryFactory* TwoBodyDecayTrajectoryFactory::clone ( void ) const

inlineoverridevirtual

Implements TrajectoryFactoryBase.

Definition at line 47 of file TwoBodyDecayTrajectoryFactory.cc.

References pwdgSkimBPark_cfi::beamSpot, constructTrajectories(), match(), rpcPointValidation_cfi::recHit, PDWG_EXOHSCP_cff::tracks, and TwoBodyDecayTrajectoryFactory().

47 { return new TwoBodyDecayTrajectoryFactory(*this); }

TwoBodyDecayTrajectoryFactory::TwoBodyDecayTrajectoryFactory

TwoBodyDecayTrajectoryFactory(const edm::ParameterSet &config)

Definition: TwoBodyDecayTrajectoryFactory.cc:74

const TwoBodyDecayTrajectoryFactory::ReferenceTrajectoryCollection TwoBodyDecayTrajectoryFactory::constructTrajectories	(	const ConstTrajTrackPairCollection &	tracks,
		const TwoBodyDecay &	tbd,
		const MagneticField *	magField,
		const reco::BeamSpot &	beamSpot,
		bool	setParameterErrors
	)		const

protected

Definition at line 170 of file TwoBodyDecayTrajectoryFactory.cc.

References ReferenceTrajectoryBase::Config::allowZeroMaterial, TrajectoryFactoryBase::allowZeroMaterial_, looper::config, ReferenceTrajectoryBase::Config::constructTsosWithErrors, TwoBodyDecay::covariance(), TwoBodyDecay::hasError(), ReferenceTrajectoryBase::Config::hitsAreReverse, ReferenceTrajectoryBase::Config::includeAPEs, TrajectoryFactoryBase::includeAPEs_, TrajectoryFactoryBase::innermostStateAndRecHits(), input1, input2, TwoBodyDecayTrajectoryState::isValid(), match(), TrajectoryFactoryBase::materialEffects(), TrajectoryFactoryBase::propagationDirection(), FastTrackerRecHitMaskProducer_cfi::recHits, mps_fire::result, ReferenceTrajectoryBase::setParameterErrors(), theConstructTsosWithErrorsFlag, theSecondaryMass, theUseRefittedStateFlag, trajectories(), TwoBodyDecayTrajectoryState::trajectoryStates(), ReferenceTrajectoryBase::Config::useBeamSpot, TrajectoryFactoryBase::useBeamSpot_, and ReferenceTrajectoryBase::Config::useRefittedState.

Referenced by clone(), and trajectories().

                                    {
   ReferenceTrajectoryCollection trajectories;
 
   // get innermost valid trajectory state and hits from the tracks
   TrajectoryInput input1 = this->innermostStateAndRecHits(tracks[0]);
   TrajectoryInput input2 = this->innermostStateAndRecHits(tracks[1]);
 
   if (!(input1.first.isValid() && input2.first.isValid()))
     return trajectories;
 
   // produce TwoBodyDecayTrajectoryState (input for TwoBodyDecayTrajectory)
   TsosContainer tsos(input1.first, input2.first);
   ConstRecHitCollection recHits(input1.second, input2.second);
   TwoBodyDecayTrajectoryState trajectoryState(tsos, tbd, theSecondaryMass, magField);
 
   if (!trajectoryState.isValid()) {
     trajectories.push_back(ReferenceTrajectoryPtr(new TwoBodyDecayTrajectory()));
     return trajectories;
   }
 
   // always use the refitted trajectory state for matching
   // FIXME FIXME CLONE
   //TrackingRecHit::ConstRecHitPointer updatedRecHit1( recHits.first.front()->clone( tsos.first ) );
   // TrackingRecHit::ConstRecHitPointer updatedRecHit2( recHits.second.front()->clone( tsos.second ) );
 
   TrackingRecHit::ConstRecHitPointer updatedRecHit1(recHits.first.front());
   TrackingRecHit::ConstRecHitPointer updatedRecHit2(recHits.second.front());
 
   bool valid1 = match(trajectoryState.trajectoryStates(true).first, updatedRecHit1);
 
   bool valid2 = match(trajectoryState.trajectoryStates(true).second, updatedRecHit2);
 
   if (!valid1 || !valid2) {
     trajectories.push_back(ReferenceTrajectoryPtr(new TwoBodyDecayTrajectory()));
     return trajectories;
   }
 
   ReferenceTrajectoryBase::Config config(materialEffects(), propagationDirection());
   config.useBeamSpot = useBeamSpot_;
   config.includeAPEs = includeAPEs_;
   config.allowZeroMaterial = allowZeroMaterial_;
   config.useRefittedState = theUseRefittedStateFlag;
   config.constructTsosWithErrors = theConstructTsosWithErrorsFlag;
   // set the flag for reversing the RecHits to false, since they are already in the correct order.
   config.hitsAreReverse = false;
   TwoBodyDecayTrajectory *result = new TwoBodyDecayTrajectory(trajectoryState, recHits, magField, beamSpot, config);
   if (setParameterErrors && tbd.hasError())
     result->setParameterErrors(tbd.covariance());
   trajectories.push_back(ReferenceTrajectoryPtr(result));
   return trajectories;
 }

bool TwoBodyDecayTrajectoryFactory::match	(	const TrajectoryStateOnSurface &	state,
		const TransientTrackingRecHit::ConstRecHitPointer &	recHit
	)		const

protected

Definition at line 227 of file TwoBodyDecayTrajectoryFactory.cc.

References DEFINE_EDM_PLUGIN, TrajectoryStateOnSurface::localPosition(), mathSSE::sqrt(), theNSigmaCutValue, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by clone(), and constructTrajectories().

                                                                                                          {
   LocalPoint lp1 = state.localPosition();
   LocalPoint lp2 = recHit->localPosition();
 
   double deltaX = lp1.x() - lp2.x();
   double deltaY = lp1.y() - lp2.y();
 
   LocalError le = recHit->localPositionError();
 
   double varX = le.xx();
   double varY = le.yy();
 
   return ((fabs(deltaX) / sqrt(varX) < theNSigmaCutValue) && (fabs(deltaY) / sqrt(varY) < theNSigmaCutValue));
 }

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

overridevirtual

Produce the trajectories.

Implements TrajectoryFactoryBase.

Definition at line 89 of file TwoBodyDecayTrajectoryFactory.cc.

References pwdgSkimBPark_cfi::beamSpot, TwoBodyDecay::chi2(), constructTrajectories(), TwoBodyDecayFitter::estimate(), edm::EventSetup::get(), TwoBodyDecay::isValid(), HLT_2018_cff::magneticField, edm::ESHandle< T >::product(), edm::second(), reco::TransientTrack::setES(), theChi2CutValue, theFitter, thePrimaryMass, thePrimaryWidth, and theSecondaryMass.

Referenced by constructTrajectories(), and trajectories().

                                                                                                                 {
   ReferenceTrajectoryCollection trajectories;
 
   edm::ESHandle<MagneticField> magneticField;
   setup.get<IdealMagneticFieldRecord>().get(magneticField);
 
   if (tracks.size() == 2) {
     // produce transient tracks from persistent tracks
     std::vector<reco::TransientTrack> transientTracks(2);
 
     transientTracks[0] = reco::TransientTrack(*tracks[0].second, magneticField.product());
     transientTracks[0].setES(setup);
 
     transientTracks[1] = reco::TransientTrack(*tracks[1].second, magneticField.product());
     transientTracks[1].setES(setup);
 
     // estimate the decay parameters
     VirtualMeasurement vm(thePrimaryMass, thePrimaryWidth, theSecondaryMass, beamSpot);
     TwoBodyDecay tbd = theFitter.estimate(transientTracks, vm);
 
     if (!tbd.isValid() || (tbd.chi2() > theChi2CutValue)) {
       trajectories.push_back(ReferenceTrajectoryPtr(new TwoBodyDecayTrajectory()));
       return trajectories;
     }
 
     return constructTrajectories(tracks, tbd, magneticField.product(), beamSpot, false);
   } else {
     edm::LogInfo("ReferenceTrajectories") << "@SUB=TwoBodyDecayTrajectoryFactory::trajectories"
                                           << "Need 2 tracks, got " << tracks.size() << ".\n";
   }
 
   return trajectories;
 }

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

overridevirtual

Implements TrajectoryFactoryBase.

Definition at line 124 of file TwoBodyDecayTrajectoryFactory.cc.

References pwdgSkimBPark_cfi::beamSpot, TwoBodyDecay::chi2(), constructTrajectories(), TwoBodyDecayFitter::estimate(), edm::EventSetup::get(), TwoBodyDecay::isValid(), HLT_2018_cff::magneticField, edm::ESHandle< T >::product(), edm::second(), reco::TransientTrack::setES(), theChi2CutValue, theFitter, thePrimaryMass, thePrimaryWidth, theSecondaryMass, and trajectories().

                                         {
   ReferenceTrajectoryCollection trajectories;
 
   edm::ESHandle<MagneticField> magneticField;
   setup.get<IdealMagneticFieldRecord>().get(magneticField);
 
   if (tracks.size() == 2 && external.size() == 2) {
     if (external[0].isValid() && external[1].isValid())  // Include external estimates
     {
       // produce transient tracks from persistent tracks
       std::vector<reco::TransientTrack> transientTracks(2);
 
       transientTracks[0] = reco::TransientTrack(*tracks[0].second, magneticField.product());
       transientTracks[0].setES(setup);
 
       transientTracks[1] = reco::TransientTrack(*tracks[1].second, magneticField.product());
       transientTracks[1].setES(setup);
 
       // estimate the decay parameters. the transient tracks are not really associated to the
       // the external tsos, but this is o.k., because the only information retrieved from them
       // is the magnetic field.
       VirtualMeasurement vm(thePrimaryMass, thePrimaryWidth, theSecondaryMass, beamSpot);
       TwoBodyDecay tbd = theFitter.estimate(transientTracks, external, vm);
 
       if (!tbd.isValid() || (tbd.chi2() > theChi2CutValue)) {
         trajectories.push_back(ReferenceTrajectoryPtr(new TwoBodyDecayTrajectory()));
         return trajectories;
       }
 
       return constructTrajectories(tracks, tbd, magneticField.product(), beamSpot, true);
     } else {
       // Return without external estimate
       trajectories = this->trajectories(setup, tracks, beamSpot);
     }
   } else {
     edm::LogInfo("ReferenceTrajectories") << "@SUB=TwoBodyDecayTrajectoryFactory::trajectories"
                                           << "Need 2 tracks, got " << tracks.size() << ".\n";
   }
 
   return trajectories;
 }