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
virtual TwoBodyDecayTrajectoryFactory *	clone () const

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

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

	TwoBodyDecayTrajectoryFactory (const edm::ParameterSet &config)

	~TwoBodyDecayTrajectoryFactory ()

Public Member Functions inherited from TrajectoryFactoryBase
const MaterialEffects	materialEffects (void) const

const PropagationDirection	propagationDirection (void) const

	TrajectoryFactoryBase (const edm::ParameterSet &config)

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

void	produceVirtualMeasurement (const edm::ParameterSet &config)

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
bool	theConstructTsosWithErrorsFlag

TwoBodyDecayFitter *	theFitter

double	theNSigmaCutValue

bool	theUseRefittedStateFlag

VirtualMeasurement	theVM

Protected Attributes inherited from TrajectoryFactoryBase
bool	theUseBeamSpot

Detailed Description

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

Definition at line 31 of file TwoBodyDecayTrajectoryFactory.cc.

Member Typedef Documentation

typedef TwoBodyDecayTrajectory::ConstRecHitCollection TwoBodyDecayTrajectoryFactory::ConstRecHitCollection

Definition at line 37 of file TwoBodyDecayTrajectoryFactory.cc.

typedef TwoBodyDecayTrajectoryState::TsosContainer TwoBodyDecayTrajectoryFactory::TsosContainer

Definition at line 36 of file TwoBodyDecayTrajectoryFactory.cc.

typedef TwoBodyDecayVirtualMeasurement TwoBodyDecayTrajectoryFactory::VirtualMeasurement

Definition at line 35 of file TwoBodyDecayTrajectoryFactory.cc.

Constructor & Destructor Documentation

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

Definition at line 80 of file TwoBodyDecayTrajectoryFactory.cc.

References produceVirtualMeasurement().

Referenced by clone().

   : TrajectoryFactoryBase( config ),
     theFitter( new TwoBodyDecayFitter( config ) ),
     theNSigmaCutValue( config.getParameter< double >( "NSigmaCut" ) ),
     theUseRefittedStateFlag( config.getParameter< bool >( "UseRefittedState" ) ),
     theConstructTsosWithErrorsFlag( config.getParameter< bool >( "ConstructTsosWithErrors" ) )
 {
   produceVirtualMeasurement( config );
 }

TwoBodyDecayTrajectoryFactory::~TwoBodyDecayTrajectoryFactory ( )

inline

Definition at line 40 of file TwoBodyDecayTrajectoryFactory.cc.

40 {}

Member Function Documentation

virtual TwoBodyDecayTrajectoryFactory* TwoBodyDecayTrajectoryFactory::clone ( void ) const

inlinevirtual

Implements TrajectoryFactoryBase.

Definition at line 52 of file TwoBodyDecayTrajectoryFactory.cc.

References TwoBodyDecayTrajectoryFactory().

52 { return new TwoBodyDecayTrajectoryFactory(*this); }

TwoBodyDecayTrajectoryFactory::TwoBodyDecayTrajectoryFactory

TwoBodyDecayTrajectoryFactory(const edm::ParameterSet &config)

Definition: TwoBodyDecayTrajectoryFactory.cc:80

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 187 of file TwoBodyDecayTrajectoryFactory.cc.

References TwoBodyDecay::covariance(), TwoBodyDecay::hasError(), TrajectoryFactoryBase::innermostStateAndRecHits(), input1, input2, match(), TrajectoryFactoryBase::materialEffects(), TrajectoryFactoryBase::propagationDirection(), query::result, TwoBodyDecayVirtualMeasurement::secondaryMass(), ReferenceTrajectoryBase::setParameterErrors(), theConstructTsosWithErrorsFlag, theUseRefittedStateFlag, theVM, and trajectories().

Referenced by 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, theVM.secondaryMass(), magField );
 
   if ( !trajectoryState.isValid() )
   {
     trajectories.push_back( ReferenceTrajectoryPtr( new TwoBodyDecayTrajectory() ) );
     return trajectories;
   }
 
   // always use the refitted trajectory state for matching
   TransientTrackingRecHit::ConstRecHitPointer updatedRecHit1( recHits.first.front()->clone( tsos.first ) );
   bool valid1 = match( trajectoryState.trajectoryStates( true ).first,
                updatedRecHit1 );
 
   TransientTrackingRecHit::ConstRecHitPointer updatedRecHit2( recHits.second.front()->clone( tsos.second ) );
   bool valid2 = match( trajectoryState.trajectoryStates( true ).second,
                updatedRecHit2 );
 
   if ( !valid1 || !valid2 )
   {
     trajectories.push_back( ReferenceTrajectoryPtr( new TwoBodyDecayTrajectory() ) );
     return trajectories;
   }
 
   // set the flag for reversing the RecHits to false, since they are already in the correct order.
   TwoBodyDecayTrajectory* result = new TwoBodyDecayTrajectory( trajectoryState, recHits, magField,
                                    materialEffects(), propagationDirection(),
                                    false, beamSpot, theUseRefittedStateFlag,
                                    theConstructTsosWithErrorsFlag );
   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 238 of file TwoBodyDecayTrajectoryFactory.cc.

References AlignmentPositionError::globalError(), asciidump::le, TrajectoryStateOnSurface::localPosition(), mathSSE::sqrt(), theNSigmaCutValue, ErrorFrameTransformer::transform(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by 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();
 
   AlignmentPositionError* gape = recHit->det()->alignmentPositionError();
   if ( gape )
   {
     ErrorFrameTransformer eft;
     LocalError lape = eft.transform( gape->globalError(), recHit->det()->surface() );
 
     varX += lape.xx();
     varY += lape.yy();
   }
 
   return ( ( fabs(deltaX)/sqrt(varX) < theNSigmaCutValue ) && ( fabs(deltaY)/sqrt(varY) < theNSigmaCutValue ) );
 }

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

protected

Definition at line 266 of file TwoBodyDecayTrajectoryFactory.cc.

References edm::ParameterSet::getParameter(), and theVM.

Referenced by TwoBodyDecayTrajectoryFactory().

 {
   const edm::ParameterSet bsc = config.getParameter< edm::ParameterSet >( "BeamSpot" );
   const edm::ParameterSet ppc = config.getParameter< edm::ParameterSet >( "ParticleProperties" );
   // FIXME: Should get 3D beamspot and errors from BeamSpot input from
   //        event with extrapolation of tracks to beam line?
 
   GlobalPoint theBeamSpot( bsc.getParameter< double >( "MeanX" ),
                bsc.getParameter< double >( "MeanY" ),
                bsc.getParameter< double >( "MeanZ" ) );
 
   GlobalError theBeamSpotError( bsc.getParameter< double >( "VarXX" ),
                 bsc.getParameter< double >( "VarXY" ),
                 bsc.getParameter< double >( "VarYY" ),
                 bsc.getParameter< double >( "VarXZ" ),
                 bsc.getParameter< double >( "VarYZ" ),
                 bsc.getParameter< double >( "VarZZ" ) );
 
   theVM = VirtualMeasurement( ppc.getParameter< double >( "PrimaryMass" ),
                   ppc.getParameter< double >( "PrimaryWidth" ),
                   ppc.getParameter< double >( "SecondaryMass" ),
                   theBeamSpot, theBeamSpotError );
   return;
 }

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

virtual

Produce the trajectories.

Implements TrajectoryFactoryBase.

Definition at line 92 of file TwoBodyDecayTrajectoryFactory.cc.

References constructTrajectories(), TwoBodyDecayFitter::estimate(), funct::false, edm::EventSetup::get(), TwoBodyDecay::isValid(), edm::ESHandle< class >::product(), edm::second(), reco::TransientTrack::setES(), theFitter, and theVM.

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
     TwoBodyDecay tbd = theFitter->estimate( transientTracks, theVM );
 
     if ( !tbd.isValid() )
     {
       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

virtual

Implements TrajectoryFactoryBase.

Definition at line 134 of file TwoBodyDecayTrajectoryFactory.cc.

References constructTrajectories(), TwoBodyDecayFitter::estimate(), edm::EventSetup::get(), TwoBodyDecay::isValid(), edm::ESHandle< class >::product(), edm::second(), reco::TransientTrack::setES(), theFitter, theVM, trajectories(), and funct::true.

 {
   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.
       TwoBodyDecay tbd = theFitter->estimate( transientTracks, external, theVM );
 
       if ( !tbd.isValid() )
       {
     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;
 }