#include <SeedFromConsecutiveHitsCreator.h>

Inheritance diagram for SeedFromConsecutiveHitsCreator:

Public Member Functions
virtual void	init (const TrackingRegion &region, const edm::EventSetup &es, const SeedComparitor *filter) final

virtual void	makeSeed (TrajectorySeedCollection &seedCollection, const SeedingHitSet &hits) final

	SeedFromConsecutiveHitsCreator (const edm::ParameterSet &cfg)

virtual	~SeedFromConsecutiveHitsCreator ()

Public Member Functions inherited from SeedCreator
virtual	~SeedCreator ()

Static Public Member Functions
static void	fillDescriptions (edm::ParameterSetDescription &desc)

static const char *	fillDescriptionsLabel ()

Protected Attributes
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

float	theBOFFMomentum

float	theMinOneOverPtError

float	theOriginTransverseErrorMultiplier

std::string	thePropagatorLabel

edm::ESHandle< TrackerGeometry >	tracker

std::string	TTRHBuilder

Private Member Functions
void	buildSeed (TrajectorySeedCollection &seedCollection, const SeedingHitSet &hits, const FreeTrajectoryState &fts) const

bool	checkHit (const TrajectoryStateOnSurface &tsos, SeedingHitSet::ConstRecHitPointer hit) const

CurvilinearTrajectoryError	initialError (float sin2Theta) const

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

SeedingHitSet::RecHitPointer	refitHit (SeedingHitSet::ConstRecHitPointer hit, const TrajectoryStateOnSurface &state) const

Detailed Description

Definition at line 21 of file SeedFromConsecutiveHitsCreator.h.

Constructor & Destructor Documentation

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

inline

Definition at line 24 of file SeedFromConsecutiveHitsCreator.h.

       : thePropagatorLabel                (cfg.getParameter<std::string>("propagator"))
       , theBOFFMomentum                   (cfg.getParameter<double>("SeedMomentumForBOFF"))
       , theOriginTransverseErrorMultiplier(cfg.getParameter<double>("OriginTransverseErrorMultiplier"))
       , theMinOneOverPtError              (cfg.getParameter<double>("MinOneOverPtError"))
       , TTRHBuilder                       (cfg.getParameter<std::string>("TTRHBuilder"))
       , mfName_(cfg.getParameter<std::string>("magneticField"))
       , forceKinematicWithRegionDirection_(cfg.getParameter<bool>("forceKinematicWithRegionDirection"))
       {}

SeedFromConsecutiveHitsCreator::~SeedFromConsecutiveHitsCreator ( )

virtual

Definition at line 26 of file SeedFromConsecutiveHitsCreator.cc.

26 {}

Member Function Documentation

void SeedFromConsecutiveHitsCreator::buildSeed	(	TrajectorySeedCollection &	seedCollection,
		const SeedingHitSet &	hits,
		const FreeTrajectoryState &	fts
	)		const

private

Definition at line 158 of file SeedFromConsecutiveHitsCreator.cc.

References alongMomentum, checkHit(), TrackingRecHit::geographicalId(), TrackerGeometry::idToDet(), TrajectoryStateOnSurface::isValid(), eostools::move(), trajectoryStateTransform::persistentState(), Propagator::propagate(), PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, edm::OwnVector< T, P >::push_back(), DetId::rawId(), refitHit(), SeedingHitSet::size(), tracker, KFUpdator::update(), TrackInfoProducer_cfi::updatedState, and gsfElectronCkfTrackCandidateMaker_cff::updator.

Referenced by makeSeed().

 {
   const Propagator*  propagator = &(*propagatorHandle);
   
   // get updator
   KFUpdator  updator;
   
   // Now update initial state track using information from seed hits.
   
   TrajectoryStateOnSurface updatedState;
   edm::OwnVector<TrackingRecHit> seedHits;
   
   const TrackingRecHit* hit = 0;
   for ( unsigned int iHit = 0; iHit < hits.size(); iHit++) {
     hit = hits[iHit]->hit();
     TrajectoryStateOnSurface state = (iHit==0) ? 
       propagator->propagate(fts,tracker->idToDet(hit->geographicalId())->surface())
       : propagator->propagate(updatedState, tracker->idToDet(hit->geographicalId())->surface());
     if (!state.isValid()) return;
     
     SeedingHitSet::ConstRecHitPointer   tth = hits[iHit]; 
     
     std::unique_ptr<BaseTrackerRecHit> newtth(refitHit( tth, state));
     
     if (!checkHit(state,&*newtth)) return;
 
     updatedState =  updator.update(state, *newtth);
     if (!updatedState.isValid()) return;
     
     seedHits.push_back(newtth.release());
 
   } 
 
   if(!hit) return;
   
   PTrajectoryStateOnDet const & PTraj = 
     trajectoryStateTransform::persistentState(updatedState, hit->geographicalId().rawId());
   seedCollection.emplace_back(PTraj,std::move(seedHits),alongMomentum);
 
 }

bool SeedFromConsecutiveHitsCreator::checkHit	(	const TrajectoryStateOnSurface &	tsos,
		SeedingHitSet::ConstRecHitPointer	hit
	)		const

private

Definition at line 210 of file SeedFromConsecutiveHitsCreator.cc.

References SeedComparitor::compatible(), and filter.

Referenced by buildSeed().

 { 
     return (filter ? filter->compatible(tsos,hit) : true);
 }

void SeedFromConsecutiveHitsCreator::fillDescriptions ( edm::ParameterSetDescription & desc )

static

Definition at line 28 of file SeedFromConsecutiveHitsCreator.cc.

References edm::ParameterSetDescription::add(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by SeedFromConsecutiveHitsTripletOnlyCreator::fillDescriptions().

                                                                                       {
   desc.add<std::string>("propagator", "PropagatorWithMaterialParabolicMf");
   desc.add<double>("SeedMomentumForBOFF", 5.0);
   desc.add<double>("OriginTransverseErrorMultiplier", 1.0);
   desc.add<double>("MinOneOverPtError", 1.0);
   desc.add<std::string>("TTRHBuilder", "WithTrackAngle");
   desc.add<std::string>("magneticField", "ParabolicMf");
   desc.add<bool>("forceKinematicWithRegionDirection", false);
 }

static const char* SeedFromConsecutiveHitsCreator::fillDescriptionsLabel ( )

inlinestatic

Definition at line 38 of file SeedFromConsecutiveHitsCreator.h.

References dso_hidden, ALCARECOTkAlBeamHalo_cff::filter, hfClusterShapes_cfi::hits, and init.

38 { return "ConsecutiveHits"; }

void SeedFromConsecutiveHitsCreator::init	(	const TrackingRegion &	region,
		const edm::EventSetup &	es,
		const SeedComparitor *	filter
	)

finalvirtual

Implements SeedCreator.

Definition at line 38 of file SeedFromConsecutiveHitsCreator.cc.

References bfield, cloner, filter, edm::EventSetup::get(), isBOFF, mfName_, nomField, MagneticField::nominalValue(), edm::ESHandle< T >::product(), propagatorHandle, region, thePropagatorLabel, tracker, and TTRHBuilder.

                                      {
   region = &iregion;
   filter = ifilter;
   // get tracker
   es.get<TrackerDigiGeometryRecord>().get(tracker);
   // get propagator
   es.get<TrackingComponentsRecord>().get(thePropagatorLabel, propagatorHandle);
   // mag field
   es.get<IdealMagneticFieldRecord>().get(mfName_, bfield);  
   //  edm::ESInputTag mfESInputTag(mfName_);
   //  es.get<IdealMagneticFieldRecord>().get(mfESInputTag, bfield);  
   nomField = bfield->nominalValue();
   isBOFF = (0==nomField);
 
   edm::ESHandle<TransientTrackingRecHitBuilder> builderH;
   es.get<TransientRecHitRecord>().get(TTRHBuilder, builderH);
   auto builder = (TkTransientTrackingRecHitBuilder const *)(builderH.product());
   cloner = (*builder).cloner();
 
 }

CurvilinearTrajectoryError SeedFromConsecutiveHitsCreator::initialError ( float sin2Theta ) const

private

Definition at line 135 of file SeedFromConsecutiveHitsCreator.cc.

References patCaloMETCorrections_cff::C, CurvilinearTrajectoryError::matrix(), hpstanc_transforms::max, TrackingRegion::originRBound(), TrackingRegion::originZBound(), TrackingRegion::ptMin(), region, funct::sqr(), theMinOneOverPtError, and theOriginTransverseErrorMultiplier.

Referenced by makeSeed().

 {
   // Set initial uncertainty on track parameters, using only P.V. constraint and no hit
   // information.
   CurvilinearTrajectoryError newError;  // zeroed
   auto & C = newError.matrix();
 
 // FIXME: minC00. Prevent apriori uncertainty in 1/P from being too small, 
 // to avoid instabilities.
 // N.B. This parameter needs optimising ...
   // Probably OK based on quick study: KS 22/11/12.
   auto sin2th = sin2Theta;
   auto minC00 = sqr(theMinOneOverPtError);
   C[0][0] = std::max(sin2th/sqr(region->ptMin()), minC00);
   auto zErr = sqr(region->originZBound());
   auto transverseErr = sqr(theOriginTransverseErrorMultiplier*region->originRBound());
   C[1][1] = C[2][2] = 1.;  // no good reason. no bad reason....
   C[3][3] = transverseErr;
   C[4][4] = zErr*sin2th + transverseErr*(1.f-sin2th);
 
   return newError;
 }

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

privatevirtual

Reimplemented in SeedFromConsecutiveHitsTripletOnlyCreator, and SeedFromConsecutiveHitsStraightLineCreator.

Definition at line 106 of file SeedFromConsecutiveHitsCreator.cc.

References bfield, GlobalTrajectoryParameters::charge(), SeedComparitor::compatible(), filter, isBOFF, GlobalTrajectoryParameters::momentum(), nomField, TrackingRegion::origin(), GlobalTrajectoryParameters::position(), region, theBOFFMomentum, Vector3DBase< T, FrameTag >::unit(), and unlikely.

Referenced by makeSeed(), SeedFromConsecutiveHitsStraightLineCreator::~SeedFromConsecutiveHitsStraightLineCreator(), and SeedFromConsecutiveHitsTripletOnlyCreator::~SeedFromConsecutiveHitsTripletOnlyCreator().

                                                                {
 
   SeedingHitSet::ConstRecHitPointer tth1 = hits[0];
   SeedingHitSet::ConstRecHitPointer tth2 = hits[1];
   
   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.f/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) : true); 
 }

void SeedFromConsecutiveHitsCreator::makeSeed	(	TrajectorySeedCollection &	seedCollection,
		const SeedingHitSet &	hits
	)

finalvirtual

Implements SeedCreator.

Definition at line 61 of file SeedFromConsecutiveHitsCreator.cc.

References buildSeed(), patCaloMETCorrections_cff::C, FreeTrajectoryState::charge(), constexpr, TrackingRegion::direction(), relativeConstraints::error, forceKinematicWithRegionDirection_, initialError(), initialKinematic(), TkTrackingRegionsMargin< T >::left(), PV3DBase< T, PVType, FrameType >::mag(), PV3DBase< T, PVType, FrameType >::mag2(), GlobalTrajectoryParameters::magneticField(), CurvilinearTrajectoryError::matrix(), hpstanc_transforms::max, GlobalTrajectoryParameters::momentum(), FreeTrajectoryState::momentum(), TrackingRegion::origin(), TrackingRegion::originRBound(), TrackingRegion::originZBound(), FreeTrajectoryState::parameters(), PV3DBase< T, PVType, FrameType >::perp2(), RectangularEtaPhiTrackingRegion::phiMargin(), position, ALCARECOTkAlBeamHalo_cff::ptMin, TrackingRegion::ptMin(), region, TkTrackingRegionsMargin< T >::right(), SeedingHitSet::size(), funct::sqr(), RectangularEtaPhiTrackingRegion::tanLambdaRange(), and PV3DBase< T, PVType, FrameType >::x().

                                                       {
   if ( hits.size() < 2) return;
 
   GlobalTrajectoryParameters kine;
   if (!initialKinematic(kine, hits)) return;
 
   auto sin2Theta = kine.momentum().perp2()/kine.momentum().mag2();
 
   CurvilinearTrajectoryError error = initialError(sin2Theta);
   FreeTrajectoryState fts(kine, error);
 
   if(region->direction().x()!=0 && forceKinematicWithRegionDirection_) // a direction was given, check if it is an etaPhi region
   {
      const RectangularEtaPhiTrackingRegion * etaPhiRegion =  dynamic_cast<const RectangularEtaPhiTrackingRegion *>(region);
      if(etaPhiRegion) {  
     //the following completely reset the kinematics, perhaps it makes no sense and newKine=kine would do better 
     GlobalVector direction=region->direction()/region->direction().mag();
     GlobalVector momentum=direction*fts.momentum().mag();
     GlobalPoint position=region->origin()+5*direction;  
     GlobalTrajectoryParameters newKine(position,momentum,fts.charge(),&fts.parameters().magneticField());
 
     auto ptMin = region->ptMin();
     CurvilinearTrajectoryError newError; //zeroed
         auto & C = newError.matrix();
     constexpr float minC00 = 0.4f;
     C[0][0] = std::max(sin2Theta/sqr(ptMin), minC00);
     auto zErr = sqr(region->originZBound());
     auto transverseErr = sqr(region->originRBound()); // assume equal cxx cyy
         auto twiceDeltaLambda =  std::atan(etaPhiRegion->tanLambdaRange().first)-std::atan(etaPhiRegion->tanLambdaRange().second);
         auto twiceDeltaPhi = etaPhiRegion->phiMargin().right()+etaPhiRegion->phiMargin().left();
     C[1][1] = twiceDeltaLambda*twiceDeltaLambda; //2 sigma of what given in input
     C[2][2] = twiceDeltaPhi*twiceDeltaPhi;
     C[3][3] = transverseErr;
     C[4][4] = zErr*sin2Theta + transverseErr*(1.f-sin2Theta);
     fts =  FreeTrajectoryState(newKine,newError);
     }
   }
 
 
   buildSeed(seedCollection,hits,fts); 
 }

SeedingHitSet::RecHitPointer SeedFromConsecutiveHitsCreator::refitHit	(	SeedingHitSet::ConstRecHitPointer	hit,
		const TrajectoryStateOnSurface &	state
	)		const

private

Definition at line 203 of file SeedFromConsecutiveHitsCreator.cc.

References cloner.

Referenced by buildSeed().

 {
   return (SeedingHitSet::RecHitPointer)(cloner(*hit,state));
 }

Member Data Documentation

edm::ESHandle<MagneticField> SeedFromConsecutiveHitsCreator::bfield

protected

Definition at line 83 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init(), SeedFromConsecutiveHitsStraightLineCreator::initialKinematic(), SeedFromConsecutiveHitsTripletOnlyCreator::initialKinematic(), and initialKinematic().

TkClonerImpl SeedFromConsecutiveHitsCreator::cloner

protected

Definition at line 90 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init(), and refitHit().

const SeedComparitor* SeedFromConsecutiveHitsCreator::filter = 0

protected

Definition at line 80 of file SeedFromConsecutiveHitsCreator.h.

Referenced by checkHit(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), init(), SeedFromConsecutiveHitsTripletOnlyCreator::initialKinematic(), initialKinematic(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

bool SeedFromConsecutiveHitsCreator::forceKinematicWithRegionDirection_

protected

Definition at line 88 of file SeedFromConsecutiveHitsCreator.h.

Referenced by makeSeed().

bool SeedFromConsecutiveHitsCreator::isBOFF = false

protected

Definition at line 85 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init(), SeedFromConsecutiveHitsTripletOnlyCreator::initialKinematic(), and initialKinematic().

std::string SeedFromConsecutiveHitsCreator::mfName_

protected

Definition at line 87 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init().

float SeedFromConsecutiveHitsCreator::nomField

protected

Definition at line 84 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init(), SeedFromConsecutiveHitsTripletOnlyCreator::initialKinematic(), and initialKinematic().

edm::ESHandle<Propagator> SeedFromConsecutiveHitsCreator::propagatorHandle

protected

Definition at line 82 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init().

const TrackingRegion* SeedFromConsecutiveHitsCreator::region = 0

protected

Definition at line 79 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init(), initialError(), SeedFromConsecutiveHitsStraightLineCreator::initialKinematic(), SeedFromConsecutiveHitsTripletOnlyCreator::initialKinematic(), initialKinematic(), and makeSeed().

float SeedFromConsecutiveHitsCreator::theBOFFMomentum

protected

Definition at line 75 of file SeedFromConsecutiveHitsCreator.h.

Referenced by SeedFromConsecutiveHitsTripletOnlyCreator::initialKinematic(), and initialKinematic().

float SeedFromConsecutiveHitsCreator::theMinOneOverPtError

protected

Definition at line 77 of file SeedFromConsecutiveHitsCreator.h.

Referenced by initialError().

float SeedFromConsecutiveHitsCreator::theOriginTransverseErrorMultiplier

protected

Definition at line 76 of file SeedFromConsecutiveHitsCreator.h.

Referenced by initialError().

std::string SeedFromConsecutiveHitsCreator::thePropagatorLabel

protected

Definition at line 74 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init().

edm::ESHandle<TrackerGeometry> SeedFromConsecutiveHitsCreator::tracker

protected

Definition at line 81 of file SeedFromConsecutiveHitsCreator.h.

Referenced by buildSeed(), and init().

std::string SeedFromConsecutiveHitsCreator::TTRHBuilder

protected

Definition at line 86 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init().

Public Member Functions

Static Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation