/data/refman/pasoursint/CMSSW_5_3_3/src/Alignment/ReferenceTrajectories/src/DualBzeroReferenceTrajectory.cc

Go to the documentation of this file.

#include "Alignment/ReferenceTrajectories/interface/DualBzeroReferenceTrajectory.h"

#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"

#include "DataFormats/CLHEP/interface/AlgebraicObjects.h" 
#include "DataFormats/TrajectoryState/interface/LocalTrajectoryParameters.h"

#include "Alignment/ReferenceTrajectories/interface/BzeroReferenceTrajectory.h"


DualBzeroReferenceTrajectory::DualBzeroReferenceTrajectory( const TrajectoryStateOnSurface &referenceTsos,
                                                            const ConstRecHitContainer &forwardRecHits,
                                                            const ConstRecHitContainer &backwardRecHits,
                                                            const MagneticField *magField,
                                                            MaterialEffects materialEffects,
                                                            PropagationDirection propDir,
                                                            double mass,
                                                            double momentumEstimate, 
                                                            bool useBeamSpot,
                                                            const reco::BeamSpot &beamSpot)
  : DualReferenceTrajectory(referenceTsos.localParameters().mixedFormatVector().kSize - 1,
                            numberOfUsedRecHits(forwardRecHits) + numberOfUsedRecHits(backwardRecHits) - 1),
    theMomentumEstimate(momentumEstimate)
{
    theValidityFlag = DualReferenceTrajectory::construct(referenceTsos,
                                                         forwardRecHits,
                                                         backwardRecHits,
                                                         mass, materialEffects,
                                                         propDir, magField,
                                                         useBeamSpot, beamSpot);
}


ReferenceTrajectory*
DualBzeroReferenceTrajectory::construct(const TrajectoryStateOnSurface &referenceTsos, 
                                        const ConstRecHitContainer &recHits,
                                        double mass, MaterialEffects materialEffects,
                                        const PropagationDirection propDir,
                                        const MagneticField *magField,
                                        bool useBeamSpot,
                                        const reco::BeamSpot &beamSpot) const
{
  if (materialEffects >= breakPoints)  throw cms::Exception("BadConfig")
    << "[DualBzeroReferenceTrajectory::construct] Wrong MaterialEffects: " << materialEffects;
  
  return new BzeroReferenceTrajectory(referenceTsos, recHits,
                                      false, magField,
                                      materialEffects, propDir,
                                      mass, theMomentumEstimate, useBeamSpot, beamSpot);
}


AlgebraicVector
DualBzeroReferenceTrajectory::extractParameters(const TrajectoryStateOnSurface &referenceTsos) const
{
  AlgebraicVector param = asHepVector<5>( referenceTsos.localParameters().mixedFormatVector() );
  return param.sub( 2, 5 );
}