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 ExternalPredictionCollection &external, const reco::BeamSpot &beamSpot) const override

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

	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

◆ ConstRecHitCollection

typedef TwoBodyDecayTrajectory::ConstRecHitCollection TwoBodyDecayTrajectoryFactory::ConstRecHitCollection

Definition at line 32 of file TwoBodyDecayTrajectoryFactory.cc.

◆ TsosContainer

typedef TwoBodyDecayTrajectoryState::TsosContainer TwoBodyDecayTrajectoryFactory::TsosContainer

Definition at line 31 of file TwoBodyDecayTrajectoryFactory.cc.

◆ VirtualMeasurement

typedef TwoBodyDecayVirtualMeasurement TwoBodyDecayTrajectoryFactory::VirtualMeasurement

Definition at line 30 of file TwoBodyDecayTrajectoryFactory.cc.

Constructor & Destructor Documentation

◆ TwoBodyDecayTrajectoryFactory()

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

Definition at line 74 of file TwoBodyDecayTrajectoryFactory.cc.

     : 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");
 }

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

Referenced by clone().

◆ ~TwoBodyDecayTrajectoryFactory()

TwoBodyDecayTrajectoryFactory::~TwoBodyDecayTrajectoryFactory ( )

override

Definition at line 87 of file TwoBodyDecayTrajectoryFactory.cc.

87 {}

Member Function Documentation

◆ clone()

TwoBodyDecayTrajectoryFactory* TwoBodyDecayTrajectoryFactory::clone ( void ) const

inlineoverridevirtual

Implements TrajectoryFactoryBase.

Definition at line 47 of file TwoBodyDecayTrajectoryFactory.cc.

47 { return new TwoBodyDecayTrajectoryFactory(*this); }

References TwoBodyDecayTrajectoryFactory().

◆ constructTrajectories()

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.

                                    {
   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;
 }

References TrajectoryFactoryBase::allowZeroMaterial_, pwdgSkimBPark_cfi::beamSpot, looper::config, TwoBodyDecay::covariance(), TwoBodyDecay::hasError(), TrajectoryFactoryBase::includeAPEs_, TrajectoryFactoryBase::innermostStateAndRecHits(), input1, input2, TwoBodyDecayTrajectoryState::isValid(), match(), TrajectoryFactoryBase::materialEffects(), TrajectoryFactoryBase::propagationDirection(), FastTrackerRecHitMaskProducer_cfi::recHits, mps_fire::result, theConstructTsosWithErrorsFlag, theSecondaryMass, theUseRefittedStateFlag, trajectories(), TwoBodyDecayTrajectoryState::trajectoryStates(), and TrajectoryFactoryBase::useBeamSpot_.

Referenced by trajectories().

◆ match()

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

protected

Definition at line 227 of file TwoBodyDecayTrajectoryFactory.cc.

                                                                                                          {
   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));
 }

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

Referenced by constructTrajectories().

◆ trajectories() [1/2]

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.

                                         {
   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;
 }

References pwdgSkimBPark_cfi::beamSpot, TwoBodyDecay::chi2(), constructTrajectories(), TwoBodyDecayFitter::estimate(), get, TwoBodyDecay::isValid(), HLT_2018_cff::magneticField, edm::second(), reco::TransientTrack::setES(), singleTopDQM_cfi::setup, theChi2CutValue, theFitter, thePrimaryMass, thePrimaryWidth, theSecondaryMass, PDWG_EXOHSCP_cff::tracks, and trajectories().

◆ trajectories() [2/2]

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.

                                                                                                                 {
   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;
 }