SeedFromConsecutiveHitsCreator Class Reference

#include <SeedFromConsecutiveHitsCreator.h>

Inheritance diagram for SeedFromConsecutiveHitsCreator:

Public Member Functions
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)
virtual	~SeedFromConsecutiveHitsCreator ()
Public Member Functions inherited from SeedCreator
virtual	~SeedCreator ()

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
double	theBOFFMomentum
double	theMinOneOverPtError
double	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 25 of file SeedFromConsecutiveHitsCreator.cc.

{}

Member Function Documentation

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

private

Definition at line 149 of file SeedFromConsecutiveHitsCreator.cc.

References alongMomentum, checkHit(), SeedComparitor::compatible(), filter, TrackingRecHit::geographicalId(), TrajectoryStateOnSurface::isValid(), eostools::move(), trajectoryStateTransform::persistentState(), Propagator::propagate(), HLT_25ns14e33_v1_cff::propagator, edm::OwnVector< T, P >::push_back(), DetId::rawId(), refitHit(), fileCollector::seed, SeedingHitSet::size(), tracker, KFUpdator::update(), and HLT_25ns14e33_v1_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());

  } 

  
  PTrajectoryStateOnDet const & PTraj = 
    trajectoryStateTransform::persistentState(updatedState, hit->geographicalId().rawId());
  TrajectorySeed seed(PTraj,std::move(seedHits),alongMomentum); 
  if ( !filter || filter->compatible(seed)) seedCollection.push_back(std::move(seed));

}

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

private

Definition at line 201 of file SeedFromConsecutiveHitsCreator.cc.

References SeedComparitor::compatible(), and filter.

Referenced by buildSeed().

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

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

virtual

Implements SeedCreator.

Definition at line 27 of file SeedFromConsecutiveHitsCreator.cc.

References bfield, cloner, filter, edm::EventSetup::get(), isBOFF, mfName_, nomField, edm::ESHandle< class >::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 128 of file SeedFromConsecutiveHitsCreator.cc.

References funct::C, bookConverter::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.
  AlgebraicSymMatrix55 C = ROOT::Math::SMatrixIdentity();

// 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.
  float sin2th = sin2Theta;
  float minC00 = sqr(theMinOneOverPtError);
  C[0][0] = std::max(sin2th/sqr(region->ptMin()), minC00);
  float zErr = sqr(region->originZBound());
  float transverseErr = sqr(theOriginTransverseErrorMultiplier*region->originRBound());
  C[3][3] = transverseErr;
  C[4][4] = zErr*sin2th + transverseErr*(1.f-sin2th);

  return CurvilinearTrajectoryError(C);
}

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

privatevirtual

Reimplemented in SeedFromConsecutiveHitsStraightLineCreator, and SeedFromConsecutiveHitsTripletOnlyCreator.

Definition at line 99 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().

                                                               {

  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./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); 
}

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

virtual

Implements SeedCreator.

Definition at line 50 of file SeedFromConsecutiveHitsCreator.cc.

References buildSeed(), funct::C, FreeTrajectoryState::charge(), HLTFastRecoForTau_cff::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, TrackingRegion::direction(), relativeConstraints::error, RectangularEtaPhiTrackingRegion::etaRange(), forceKinematicWithRegionDirection_, initialError(), initialKinematic(), TkTrackingRegionsMargin< T >::left(), PV3DBase< T, PVType, FrameType >::mag(), PV3DBase< T, PVType, FrameType >::mag2(), GlobalTrajectoryParameters::magneticField(), bookConverter::max, GlobalTrajectoryParameters::momentum(), FreeTrajectoryState::momentum(), TrackingRegion::origin(), TrackingRegion::originRBound(), TrackingRegion::originZBound(), FreeTrajectoryState::parameters(), PV3DBase< T, PVType, FrameType >::perp2(), RectangularEtaPhiTrackingRegion::phiMargin(), position, PtMinSelector_cfg::ptMin, TrackingRegion::ptMin(), region, TkTrackingRegionsMargin< T >::right(), SeedingHitSet::size(), funct::sqr(), and PV3DBase< T, PVType, FrameType >::x().

                                                      {
  if ( hits.size() < 2) return;

  GlobalTrajectoryParameters kine;
  if (!initialKinematic(kine, hits)) return;

  float 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());

    GlobalError vertexErr( sqr(region->originRBound()), 0, sqr(region->originRBound()),
                   0, 0, sqr(region->originZBound()));

    float ptMin = region->ptMin();
    AlgebraicSymMatrix55 C = ROOT::Math::SMatrixIdentity();

    float minC00 = 0.4;
    C[0][0] = std::max(sin2Theta/sqr(ptMin), minC00);
    float zErr = vertexErr.czz();
    float transverseErr = vertexErr.cxx(); // assume equal cxx cyy
        float deltaEta = (etaPhiRegion->etaRange().first-etaPhiRegion->etaRange().second)/2.;
        float deltaPhi = (etaPhiRegion->phiMargin().right()+etaPhiRegion->phiMargin().left())/2.;
    C[1][1] = deltaEta*deltaEta*4; //2 sigma of what given in input
    C[2][2] = deltaPhi*deltaPhi*4;
    C[3][3] = transverseErr;
    C[4][4] = zErr*sin2Theta + transverseErr*(1-sin2Theta);
    CurvilinearTrajectoryError newError(C);
    fts =  FreeTrajectoryState(newKine,newError);
    }
  }


  buildSeed(seedCollection,hits,fts); 
}

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

private

Definition at line 194 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 80 of file SeedFromConsecutiveHitsCreator.h.

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

TkClonerImpl SeedFromConsecutiveHitsCreator::cloner

protected

Definition at line 87 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init(), and refitHit().

const SeedComparitor* SeedFromConsecutiveHitsCreator::filter = 0

protected

Definition at line 77 of file SeedFromConsecutiveHitsCreator.h.

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

bool SeedFromConsecutiveHitsCreator::forceKinematicWithRegionDirection_

protected

Definition at line 85 of file SeedFromConsecutiveHitsCreator.h.

Referenced by makeSeed().

bool SeedFromConsecutiveHitsCreator::isBOFF = false

protected

Definition at line 82 of file SeedFromConsecutiveHitsCreator.h.

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

std::string SeedFromConsecutiveHitsCreator::mfName_

protected

Definition at line 84 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init().

float SeedFromConsecutiveHitsCreator::nomField

protected

Definition at line 81 of file SeedFromConsecutiveHitsCreator.h.

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

edm::ESHandle<Propagator> SeedFromConsecutiveHitsCreator::propagatorHandle

protected

Definition at line 79 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init().

const TrackingRegion* SeedFromConsecutiveHitsCreator::region = 0

protected

Definition at line 76 of file SeedFromConsecutiveHitsCreator.h.

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

double SeedFromConsecutiveHitsCreator::theBOFFMomentum

protected

Definition at line 72 of file SeedFromConsecutiveHitsCreator.h.

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

double SeedFromConsecutiveHitsCreator::theMinOneOverPtError

protected

Definition at line 74 of file SeedFromConsecutiveHitsCreator.h.

Referenced by initialError().

double SeedFromConsecutiveHitsCreator::theOriginTransverseErrorMultiplier

protected

Definition at line 73 of file SeedFromConsecutiveHitsCreator.h.

Referenced by initialError().

std::string SeedFromConsecutiveHitsCreator::thePropagatorLabel

protected

Definition at line 71 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init().

edm::ESHandle<TrackerGeometry> SeedFromConsecutiveHitsCreator::tracker

protected

Definition at line 78 of file SeedFromConsecutiveHitsCreator.h.

Referenced by buildSeed(), and init().

std::string SeedFromConsecutiveHitsCreator::TTRHBuilder

protected

Definition at line 83 of file SeedFromConsecutiveHitsCreator.h.

Referenced by init().