SeedFromNuclearInteraction Class Reference

#include <SeedFromNuclearInteraction.h>

Public Member Functions
PropagationDirection	direction () const
recHitContainer	hits () const
const TSOS &	initialTSOS () const
bool	isValid () const
ConstRecHitPointer	outerHit () const
DetId	outerHitDetId () const
GlobalPoint	outerHitPosition () const
AlgebraicMatrix33	rotationMatrix (const GlobalVector &perp) const
	SeedFromNuclearInteraction (const Propagator *prop, const TrackerGeometry *geom, double ptMin)
void	setMeasurements (const TSOS &tsosAtInteractionPoint, ConstRecHitPointer ihit, ConstRecHitPointer ohit)
	Fill all data members from 2 TM's where the first one is supposed to be at the interaction point. More...
void	setMeasurements (TangentHelix &primHelix, const TSOS &inner_TSOS, ConstRecHitPointer ihit, ConstRecHitPointer ohit)
	Fill all data members from 1 TSOS and 2 rec Hits and using the circle associated to the primary track as constraint. More...
FreeTrajectoryState *	stateWithError () const
FreeTrajectoryState *	stateWithError (TangentHelix &helix) const
PTrajectoryStateOnDet const &	trajectoryState () const
TrajectorySeed	TrajSeed () const
const TSOS &	updatedTSOS () const
virtual	~SeedFromNuclearInteraction ()

Private Types
typedef std::vector < ConstRecHitPointer >	ConstRecHitContainer
typedef TransientTrackingRecHit::ConstRecHitPointer	ConstRecHitPointer
typedef edm::OwnVector < TrackingRecHit >	recHitContainer
typedef TrajectoryMeasurement	TM
typedef TrajectoryStateOnSurface	TSOS

Private Member Functions
bool	construct ()

Private Attributes
std::shared_ptr < FreeTrajectoryState >	freeTS_
std::shared_ptr< TSOS >	initialTSOS_
ConstRecHitPointer	innerHit_
bool	isValid_
ConstRecHitPointer	outerHit_
double	ptMin
PTrajectoryStateOnDet	pTraj
ConstRecHitContainer	theHits
const Propagator *	thePropagator
const TrackerGeometry *	theTrackerGeom
std::shared_ptr< TSOS >	updatedTSOS_

Detailed Description

Definition at line 17 of file SeedFromNuclearInteraction.h.

Member Typedef Documentation

typedef std::vector<ConstRecHitPointer> SeedFromNuclearInteraction::ConstRecHitContainer

private

Definition at line 23 of file SeedFromNuclearInteraction.h.

typedef TransientTrackingRecHit::ConstRecHitPointer SeedFromNuclearInteraction::ConstRecHitPointer

private

Definition at line 22 of file SeedFromNuclearInteraction.h.

typedef edm::OwnVector<TrackingRecHit> SeedFromNuclearInteraction::recHitContainer

private

Definition at line 21 of file SeedFromNuclearInteraction.h.

typedef TrajectoryMeasurement SeedFromNuclearInteraction::TM

private

Definition at line 19 of file SeedFromNuclearInteraction.h.

typedef TrajectoryStateOnSurface SeedFromNuclearInteraction::TSOS

private

Definition at line 20 of file SeedFromNuclearInteraction.h.

Constructor & Destructor Documentation

SeedFromNuclearInteraction::SeedFromNuclearInteraction	(	const Propagator *	prop,
		const TrackerGeometry *	geom,
		double	ptMin
	)

Definition at line 14 of file SeedFromNuclearInteraction.cc.

References freeTS_, initialTSOS_, isValid_, and updatedTSOS_.

    : ptMin(iPtMin), thePropagator(prop), theTrackerGeom(geom) {
  isValid_ = true;
  initialTSOS_ = std::make_shared<TrajectoryStateOnSurface>();
  updatedTSOS_ = std::make_shared<TrajectoryStateOnSurface>();
  freeTS_ = std::make_shared<FreeTrajectoryState>();
}

virtual SeedFromNuclearInteraction::~SeedFromNuclearInteraction ( )

inlinevirtual

Definition at line 28 of file SeedFromNuclearInteraction.h.

{}

Member Function Documentation

bool SeedFromNuclearInteraction::construct ( )

private

Definition at line 133 of file SeedFromNuclearInteraction.cc.

References freeTS_, TrackingRecHit::geographicalId(), TrackerGeometry::idToDet(), TrajectoryStateOnSurface::isValid(), LogDebug, outerHitDetId(), trajectoryStateTransform::persistentState(), Propagator::propagate(), pTraj, theHits, thePropagator, theTrackerGeom, KFUpdator::update(), and updatedTSOS_.

Referenced by setMeasurements().

                                           {
  // loop on all hits in theHits
  KFUpdator theUpdator;

  const TrackingRecHit* hit = nullptr;

  LogDebug("NuclearSeedGenerator") << "Seed ** initial state " << freeTS_->cartesianError().matrix();

  for (unsigned int iHit = 0; iHit < theHits.size(); iHit++) {
    hit = theHits[iHit]->hit();
    TrajectoryStateOnSurface state =
        (iHit == 0)
            ? thePropagator->propagate(*freeTS_, theTrackerGeom->idToDet(hit->geographicalId())->surface())
            : thePropagator->propagate(*updatedTSOS_, theTrackerGeom->idToDet(hit->geographicalId())->surface());

    if (!state.isValid())
      return false;

    const TransientTrackingRecHit::ConstRecHitPointer& tth = theHits[iHit];
    updatedTSOS_.reset(new TrajectoryStateOnSurface(theUpdator.update(state, *tth)));
  }

  LogDebug("NuclearSeedGenerator") << "Seed ** updated state " << updatedTSOS_->cartesianError().matrix();

  pTraj = trajectoryStateTransform::persistentState(*updatedTSOS_, outerHitDetId().rawId());
  return true;
}

PropagationDirection SeedFromNuclearInteraction::direction ( ) const

inline

Definition at line 45 of file SeedFromNuclearInteraction.h.

References alongMomentum.

Referenced by stateWithError(), and TrajSeed().

{ return alongMomentum; }

edm::OwnVector< TrackingRecHit > SeedFromNuclearInteraction::hits

(

void

)

const

Definition at line 162 of file SeedFromNuclearInteraction.cc.

References edm::OwnVector< T, P >::push_back(), and theHits.

Referenced by TrajSeed().

                                                                    {
  recHitContainer _hits;
  for (ConstRecHitContainer::const_iterator it = theHits.begin(); it != theHits.end(); it++) {
    _hits.push_back(it->get()->hit()->clone());
  }
  return _hits;
}

const TSOS& SeedFromNuclearInteraction::initialTSOS ( ) const

inline

Definition at line 55 of file SeedFromNuclearInteraction.h.

References initialTSOS_.

{ return *initialTSOS_; }

bool SeedFromNuclearInteraction::isValid ( void  ) const

inline

Definition at line 51 of file SeedFromNuclearInteraction.h.

References isValid_.

Referenced by ntupleDataFormat._Object::_checkIsValid(), and core.AutoHandle.AutoHandle::ReallyLoad().

{ return isValid_; }

ConstRecHitPointer SeedFromNuclearInteraction::outerHit ( ) const

inline

Definition at line 63 of file SeedFromNuclearInteraction.h.

References outerHit_.

{ return outerHit_; }

DetId SeedFromNuclearInteraction::outerHitDetId ( ) const

inline

Definition at line 61 of file SeedFromNuclearInteraction.h.

References outerHit_.

Referenced by construct(), and outerHitPosition().

{ return outerHit_->geographicalId(); }

GlobalPoint SeedFromNuclearInteraction::outerHitPosition ( ) const

inline

Definition at line 57 of file SeedFromNuclearInteraction.h.

References TrackerGeometry::idToDet(), outerHit_, outerHitDetId(), GeomDet::surface(), theTrackerGeom, and Surface::toGlobal().

Referenced by stateWithError().

                                       {
    return theTrackerGeom->idToDet(outerHitDetId())->surface().toGlobal(outerHit_->localPosition());
  }

AlgebraicMatrix33 SeedFromNuclearInteraction::rotationMatrix

(

const GlobalVector &

perp

)

const

Return the rotation matrix to be applied to get parameters in a framework where the z direction is along perp

Definition at line 170 of file SeedFromNuclearInteraction.cc.

References Vector3DBase< T, FrameTag >::cross(), mps_fire::result, Vector3DBase< T, FrameTag >::unit(), unit(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and MetAnalyzer::zAxis.

Referenced by stateWithError().

                                                                                           {
  AlgebraicMatrix33 result;

  // z axis coincides with perp
  GlobalVector zAxis = perp.unit();

  // x axis has no global Z component
  GlobalVector xAxis;
  if (zAxis.x() != 0 || zAxis.y() != 0) {
    // precision is not an issue here, just protect against divizion by zero
    xAxis = GlobalVector(-zAxis.y(), zAxis.x(), 0).unit();
  } else {  // perp coincides with global Z
    xAxis = GlobalVector(1, 0, 0);
  }

  // y axis obtained by cross product
  GlobalVector yAxis(zAxis.cross(xAxis));

  result(0, 0) = xAxis.x();
  result(0, 1) = xAxis.y();
  result(0, 2) = xAxis.z();
  result(1, 0) = yAxis.x();
  result(1, 1) = yAxis.y();
  result(1, 2) = yAxis.z();
  result(2, 0) = zAxis.x();
  result(2, 1) = zAxis.y();
  result(2, 2) = zAxis.z();
  return result;
}

void SeedFromNuclearInteraction::setMeasurements	(	const TSOS &	tsosAtInteractionPoint,
		ConstRecHitPointer	ihit,
		ConstRecHitPointer	ohit
	)

Fill all data members from 2 TM's where the first one is supposed to be at the interaction point.

Definition at line 25 of file SeedFromNuclearInteraction.cc.

References construct(), freeTS_, initialTSOS_, innerHit_, isValid_, outerHit_, ptMin, stateWithError(), and theHits.

                                                                          {
  // delete pointer to TrackingRecHits
  theHits.clear();

  // get the inner and outer transient TrackingRecHits
  innerHit_ = ihit;
  outerHit_ = ohit;

  //theHits.push_back(  inner_TM.recHit() ); // put temporarily - TODO: remove this line
  theHits.push_back(outerHit_);

  initialTSOS_.reset(new TrajectoryStateOnSurface(inner_TSOS));

  // calculate the initial FreeTrajectoryState.
  freeTS_.reset(stateWithError());

  // check transverse momentum
  if (freeTS_->momentum().perp() < ptMin) {
    isValid_ = false;
  } else {
    // convert freeTS_ into a persistent TSOS on the outer surface
    isValid_ = construct();
  }
}

void SeedFromNuclearInteraction::setMeasurements	(	TangentHelix &	primHelix,
		const TSOS &	inner_TSOS,
		ConstRecHitPointer	ihit,
		ConstRecHitPointer	ohit
	)

Fill all data members from 1 TSOS and 2 rec Hits and using the circle associated to the primary track as constraint.

Definition at line 52 of file SeedFromNuclearInteraction.cc.

References construct(), freeTS_, TrackerGeometry::idToDet(), initialTSOS_, innerHit_, isValid_, outerHit_, ptMin, stateWithError(), theHits, theTrackerGeom, and GeomDet::toGlobal().

                                                                          {
  // delete pointer to TrackingRecHits
  theHits.clear();

  // get the inner and outer transient TrackingRecHits
  innerHit_ = ihit;
  outerHit_ = ohit;

  GlobalPoint innerPos =
      theTrackerGeom->idToDet(innerHit_->geographicalId())->surface().toGlobal(innerHit_->localPosition());
  GlobalPoint outerPos =
      theTrackerGeom->idToDet(outerHit_->geographicalId())->surface().toGlobal(outerHit_->localPosition());

  TangentHelix helix(thePrimaryHelix, outerPos, innerPos);

  theHits.push_back(innerHit_);
  theHits.push_back(outerHit_);

  initialTSOS_.reset(new TrajectoryStateOnSurface(inner_TSOS));

  // calculate the initial FreeTrajectoryState from the inner and outer TM assuming that the helix equation is already known.
  freeTS_.reset(stateWithError(helix));

  if (freeTS_->momentum().perp() < ptMin) {
    isValid_ = false;
  } else {
    // convert freeTS_ into a persistent TSOS on the outer surface
    isValid_ = construct();
  }
}

FreeTrajectoryState * SeedFromNuclearInteraction::stateWithError

(

)

const

Definition at line 86 of file SeedFromNuclearInteraction.cc.

References direction(), initialTSOS_, SurfaceOrientation::inner, and outerHitPosition().

Referenced by setMeasurements().

                                                                      {
  // Calculation of the helix assuming that the secondary track has the same direction
  // than the primary track and pass through the inner and outer hits.
  GlobalVector direction = initialTSOS_->globalDirection();
  GlobalPoint inner = initialTSOS_->globalPosition();
  TangentHelix helix(direction, inner, outerHitPosition());

  return stateWithError(helix);
}

FreeTrajectoryState * SeedFromNuclearInteraction::stateWithError

(

TangentHelix &

helix

)

const

Definition at line 96 of file SeedFromNuclearInteraction.cc.

References TangentHelix::charge(), TangentHelix::directionAtVertex(), initialTSOS_, MagneticField::inTesla(), mag(), TangentHelix::rho(), makeMuonMisalignmentScenario::rot, rotationMatrix(), TangentHelix::vertexError(), TangentHelix::vertexPoint(), and z.

                                                                                         {
  //   typedef TkRotation<float> Rotation;

  GlobalVector dirAtVtx = helix.directionAtVertex();
  const MagneticField& mag = initialTSOS_->globalParameters().magneticField();

  // Get the global parameters of the trajectory
  // we assume that the magnetic field at the vertex is equal to the magnetic field at the inner TM.
  GlobalTrajectoryParameters gtp(
      helix.vertexPoint(), dirAtVtx, helix.charge(mag.inTesla(helix.vertexPoint()).z()) / helix.rho(), 0, &mag);

  // Error matrix in a frame where z is in the direction of the track at the vertex
  AlgebraicSymMatrix66 primaryError(initialTSOS_->cartesianError().matrix());
  double p_max = initialTSOS_->globalParameters().momentum().mag();
  AlgebraicMatrix33 rot = this->rotationMatrix(dirAtVtx);

  AlgebraicMatrix66 globalRotation;
  globalRotation.Place_at(rot, 0, 0);
  globalRotation.Place_at(rot, 3, 3);
  AlgebraicSymMatrix66 primaryErrorInNewFrame = ROOT::Math::Similarity(globalRotation, primaryError);

  AlgebraicSymMatrix66 secondaryErrorInNewFrame = AlgebraicMatrixID();
  double p_perp_max = 2;  // energy max of a secondary track emited perpendicularly to the
                          // primary track is +/- 2 GeV
  secondaryErrorInNewFrame(0, 0) = primaryErrorInNewFrame(0, 0) + helix.vertexError() * p_perp_max / p_max;
  secondaryErrorInNewFrame(1, 1) = primaryErrorInNewFrame(1, 1) + helix.vertexError() * p_perp_max / p_max;
  secondaryErrorInNewFrame(2, 2) = helix.vertexError() * helix.vertexError();
  secondaryErrorInNewFrame(3, 3) = p_perp_max * p_perp_max;
  secondaryErrorInNewFrame(4, 4) = p_perp_max * p_perp_max;
  secondaryErrorInNewFrame(5, 5) = p_max * p_max;

  AlgebraicSymMatrix66 secondaryError = ROOT::Math::SimilarityT(globalRotation, secondaryErrorInNewFrame);

  return new FreeTrajectoryState(gtp, CartesianTrajectoryError(secondaryError));
}

PTrajectoryStateOnDet const& SeedFromNuclearInteraction::trajectoryState ( ) const

inline

Definition at line 39 of file SeedFromNuclearInteraction.h.

References pTraj.

Referenced by TrajSeed().

{ return pTraj; }

TrajectorySeed SeedFromNuclearInteraction::TrajSeed ( ) const

inline

Definition at line 49 of file SeedFromNuclearInteraction.h.

References direction(), hits(), and trajectoryState().

{ return TrajectorySeed(trajectoryState(), hits(), direction()); }

const TSOS& SeedFromNuclearInteraction::updatedTSOS ( ) const

inline

Definition at line 53 of file SeedFromNuclearInteraction.h.

References updatedTSOS_.

{ return *updatedTSOS_; }

Member Data Documentation

std::shared_ptr<FreeTrajectoryState> SeedFromNuclearInteraction::freeTS_

private

Initial FreeTrajectoryState

Definition at line 82 of file SeedFromNuclearInteraction.h.

Referenced by construct(), SeedFromNuclearInteraction(), and setMeasurements().

std::shared_ptr<TSOS> SeedFromNuclearInteraction::initialTSOS_

private

Initial TSOS used as input

Definition at line 80 of file SeedFromNuclearInteraction.h.

Referenced by initialTSOS(), SeedFromNuclearInteraction(), setMeasurements(), and stateWithError().

ConstRecHitPointer SeedFromNuclearInteraction::innerHit_

private

Pointer to the hit of the inner TM

Definition at line 75 of file SeedFromNuclearInteraction.h.

Referenced by setMeasurements().

bool SeedFromNuclearInteraction::isValid_

private

check if the seed is valid

Definition at line 70 of file SeedFromNuclearInteraction.h.

Referenced by isValid(), SeedFromNuclearInteraction(), and setMeasurements().

ConstRecHitPointer SeedFromNuclearInteraction::outerHit_

private

Pointer to the outer hit

Definition at line 76 of file SeedFromNuclearInteraction.h.

Referenced by outerHit(), outerHitDetId(), outerHitPosition(), and setMeasurements().

double SeedFromNuclearInteraction::ptMin

private

Minimum transverse momentum of the seed

Definition at line 88 of file SeedFromNuclearInteraction.h.

Referenced by setMeasurements().

PTrajectoryStateOnDet SeedFromNuclearInteraction::pTraj

private

the final persistent TSOS

Definition at line 84 of file SeedFromNuclearInteraction.h.

Referenced by construct(), and trajectoryState().

ConstRecHitContainer SeedFromNuclearInteraction::theHits

private

all the hits to be used to update the

Definition at line 72 of file SeedFromNuclearInteraction.h.

Referenced by construct(), hits(), and setMeasurements().

const Propagator* SeedFromNuclearInteraction::thePropagator

private

Definition at line 90 of file SeedFromNuclearInteraction.h.

Referenced by construct().

const TrackerGeometry* SeedFromNuclearInteraction::theTrackerGeom

private

Definition at line 91 of file SeedFromNuclearInteraction.h.

Referenced by construct(), outerHitPosition(), and setMeasurements().

std::shared_ptr<TSOS> SeedFromNuclearInteraction::updatedTSOS_

private

Final TSOS

Definition at line 78 of file SeedFromNuclearInteraction.h.

Referenced by construct(), SeedFromNuclearInteraction(), and updatedTSOS().