SeedFromConsecutiveHitsTripletOnlyCreator Class Reference

#include <SeedFromConsecutiveHitsTripletOnlyCreator.h>

Inheritance diagram for SeedFromConsecutiveHitsTripletOnlyCreator:

Public Member Functions
	SeedFromConsecutiveHitsTripletOnlyCreator (const edm::ParameterSet &cfg)
virtual	~SeedFromConsecutiveHitsTripletOnlyCreator ()
Public Member Functions inherited from SeedFromConsecutiveHitsCreator
virtual void	init (const TrackingRegion &region, const edm::EventSetup &es, const SeedComparitor *filter)
virtual void	makeSeed (TrajectorySeedCollection &seedCollection, const SeedingHitSet &hits)
	SeedFromConsecutiveHitsCreator (const edm::ParameterSet &cfg)
	SeedFromConsecutiveHitsCreator (const std::string &propagator="PropagatorWithMaterial", double seedMomentumForBOFF=-5.0, double aOriginTransverseErrorMultiplier=1.0, double aMinOneOverPtError=1.0, const std::string &bname="WithTrackAngle")
virtual	~SeedFromConsecutiveHitsCreator ()
Public Member Functions inherited from SeedCreator
virtual	~SeedCreator ()

Private Member Functions
virtual bool	initialKinematic (GlobalTrajectoryParameters &kine, const SeedingHitSet &hits) const

Additional Inherited Members
Protected Attributes inherited from SeedFromConsecutiveHitsCreator
edm::ESHandle< MagneticField >	bfield
TkClonerImpl	cloner
const SeedComparitor *	filter = 0
bool	forceKinematicWithRegionDirection_
bool	isBOFF = false
std::string	mfName_
float	nomField
edm::ESHandle< Propagator >	propagatorHandle
const TrackingRegion *	region = 0
double	theBOFFMomentum
double	theMinOneOverPtError
double	theOriginTransverseErrorMultiplier
std::string	thePropagatorLabel
edm::ESHandle< TrackerGeometry >	tracker
std::string	TTRHBuilder
bool	useSimpleMF_

Detailed Description

Definition at line 6 of file SeedFromConsecutiveHitsTripletOnlyCreator.h.

Constructor & Destructor Documentation

SeedFromConsecutiveHitsTripletOnlyCreator::SeedFromConsecutiveHitsTripletOnlyCreator ( const edm::ParameterSet &  cfg )

inline

Definition at line 9 of file SeedFromConsecutiveHitsTripletOnlyCreator.h.

                                                                         :
    SeedFromConsecutiveHitsCreator(cfg) { }

virtual SeedFromConsecutiveHitsTripletOnlyCreator::~SeedFromConsecutiveHitsTripletOnlyCreator ( )

inlinevirtual

Definition at line 12 of file SeedFromConsecutiveHitsTripletOnlyCreator.h.

{}

Member Function Documentation

bool SeedFromConsecutiveHitsTripletOnlyCreator::initialKinematic ( GlobalTrajectoryParameters &  kine, const SeedingHitSet &  hits  ) const

privatevirtual

Reimplemented from SeedFromConsecutiveHitsCreator.

Definition at line 12 of file SeedFromConsecutiveHitsTripletOnlyCreator.cc.

References SeedFromConsecutiveHitsCreator::bfield, GlobalTrajectoryParameters::charge(), SeedComparitor::compatible(), SeedFromConsecutiveHitsCreator::filter, SeedFromConsecutiveHitsCreator::isBOFF, trackerHitRTTI::isSingleType(), GlobalTrajectoryParameters::momentum(), SeedFromConsecutiveHitsCreator::nomField, TrackingRegion::origin(), GlobalTrajectoryParameters::position(), SeedFromConsecutiveHitsCreator::region, SeedingHitSet::size(), SiStripDetId::TEC, SeedFromConsecutiveHitsCreator::theBOFFMomentum, SiStripDetId::TID, Vector3DBase< T, FrameTag >::unit(), and unlikely.

                                                                   {

  SeedingHitSet::ConstRecHitPointer tth1 = hits[0];
  SeedingHitSet::ConstRecHitPointer tth2 = hits[1];
  


  if (hits.size()==3 && !( trackerHitRTTI::isSingleType(*hits[2]) && (hits[2]->geographicalId().subdetId()==SiStripDetId::TID || 
                              hits[2]->geographicalId().subdetId()==SiStripDetId::TEC ) 
              ) ) {
    //if 3rd hit is mono and endcap pT is not well defined so take initial state from pair
    SeedingHitSet::ConstRecHitPointer tth3 = hits[2];
    FastHelix helix(tth3->globalPosition(), tth2->globalPosition(), tth1->globalPosition(), nomField, &*bfield, tth1->globalPosition());
    kine = helix.stateAtVertex();
    if unlikely(isBOFF && (theBOFFMomentum > 0)) {
      kine = GlobalTrajectoryParameters(kine.position(),
                    kine.momentum().unit() * theBOFFMomentum,
                    kine.charge(),
                    &*bfield);
    }
    return (filter ? filter->compatible(hits, kine, helix, *region) : true); 
  } 
  
  const GlobalPoint& vertexPos = region->origin();

  FastHelix helix(tth2->globalPosition(), tth1->globalPosition(), vertexPos, nomField,&*bfield);
  if (helix.isValid()) {
    kine = helix.stateAtVertex();
  } else {
    GlobalVector initMomentum(tth2->globalPosition() - vertexPos);
    initMomentum *= (100./initMomentum.perp()); 
    kine = GlobalTrajectoryParameters(vertexPos, initMomentum, 1, &*bfield);
  } 

  if unlikely(isBOFF && (theBOFFMomentum > 0)) {
      kine = GlobalTrajectoryParameters(kine.position(),
                    kine.momentum().unit() * theBOFFMomentum,
                    kine.charge(),
                    &*bfield);
  }
  return (filter ? filter->compatible(hits, kine, helix, *region) : true); 
}