/data/refman/pasoursint/CMSSW_5_3_0/src/Alignment/ReferenceTrajectories/interface/DualReferenceTrajectory.h

Go to the documentation of this file.

#ifndef Alignment_ReferenceTrajectories_DualReferenceTrajectory_H
#define Alignment_ReferenceTrajectories_DualReferenceTrajectory_H

#include "Alignment/ReferenceTrajectories/interface/ReferenceTrajectoryBase.h"
#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"

class ReferenceTrajectory;
namespace reco { class BeamSpot;}

class DualReferenceTrajectory : public ReferenceTrajectoryBase
{

public:

  typedef TransientTrackingRecHit::ConstRecHitContainer ConstRecHitContainer;

  DualReferenceTrajectory( const TrajectoryStateOnSurface &referenceTsos,
                           const ConstRecHitContainer &forwardRecHits,
                           const ConstRecHitContainer &backwardRecHits,
                           const MagneticField *magField,
                           MaterialEffects materialEffects,
                           PropagationDirection propDir,
                           double mass,
                           bool useBeamSpot,
                           const reco::BeamSpot &beamSpot);

  virtual ~DualReferenceTrajectory() {}

  virtual DualReferenceTrajectory* clone() const { return new DualReferenceTrajectory(*this); }

protected:

  DualReferenceTrajectory( unsigned int nPar, unsigned int nHits);

  virtual bool construct(const TrajectoryStateOnSurface &referenceTsos, 
                         const ConstRecHitContainer &forwardRecHits,
                         const ConstRecHitContainer &backwardRecHits,
                         double mass, MaterialEffects materialEffects,
                         const PropagationDirection propDir,
                         const MagneticField *magField,
                         bool useBeamSpot,
                         const reco::BeamSpot &beamSpot);

  virtual ReferenceTrajectory* construct(const TrajectoryStateOnSurface &referenceTsos, 
                                         const ConstRecHitContainer &recHits,
                                         double mass, MaterialEffects materialEffects,
                                         const PropagationDirection propDir,
                                         const MagneticField *magField,
                                         bool useBeamSpot,
                                         const reco::BeamSpot &beamSpot) const;

  virtual AlgebraicVector extractParameters(const TrajectoryStateOnSurface &referenceTsos) const;

  inline const PropagationDirection oppositeDirection( const PropagationDirection propDir ) const
  { return ( propDir == anyDirection ) ? anyDirection : ( ( propDir == alongMomentum ) ? oppositeToMomentum : alongMomentum ); }

};

#endif