|
|
|
by Edmund Widl, see CMS NOTE-2007/032. extension for BreakPoints or BrokenLines by Claus Kleinwort
Definition at line 31 of file TwoBodyDecayTrajectoryFactory.cc.
| typedef TwoBodyDecayTrajectory::ConstRecHitCollection TwoBodyDecayTrajectoryFactory::ConstRecHitCollection |
Definition at line 37 of file TwoBodyDecayTrajectoryFactory.cc.
Definition at line 36 of file TwoBodyDecayTrajectoryFactory.cc.
Definition at line 35 of file TwoBodyDecayTrajectoryFactory.cc.
| 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.
{}
| virtual TwoBodyDecayTrajectoryFactory* TwoBodyDecayTrajectoryFactory::clone | ( | void | ) | const [inline, virtual] |
Implements TrajectoryFactoryBase.
Definition at line 52 of file TwoBodyDecayTrajectoryFactory.cc.
References TwoBodyDecayTrajectoryFactory().
{ return new TwoBodyDecayTrajectoryFactory(*this); }
| 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 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< T >::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;
}
| 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< T >::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;
}
bool TwoBodyDecayTrajectoryFactory::theConstructTsosWithErrorsFlag [protected] |
Definition at line 72 of file TwoBodyDecayTrajectoryFactory.cc.
Referenced by constructTrajectories().
Definition at line 68 of file TwoBodyDecayTrajectoryFactory.cc.
Referenced by trajectories().
double TwoBodyDecayTrajectoryFactory::theNSigmaCutValue [protected] |
Definition at line 70 of file TwoBodyDecayTrajectoryFactory.cc.
Referenced by match().
bool TwoBodyDecayTrajectoryFactory::theUseRefittedStateFlag [protected] |
Definition at line 71 of file TwoBodyDecayTrajectoryFactory.cc.
Referenced by constructTrajectories().
Definition at line 66 of file TwoBodyDecayTrajectoryFactory.cc.
Referenced by constructTrajectories(), produceVirtualMeasurement(), and trajectories().
1.7.3