#include <MuonSeedFromRecHits.h>

Inheritance diagram for MuonSeedFromRecHits:

Public Member Functions
void	add (MuonTransientTrackingRecHit::MuonRecHitPointer hit)

void	clear ()

TrajectorySeed	createSeed (float ptmean, float sptmean, MuonTransientTrackingRecHit::ConstMuonRecHitPointer last) const

MuonTransientTrackingRecHit::ConstMuonRecHitPointer	firstRecHit () const

	MuonSeedFromRecHits ()

unsigned int	nrhit () const

void	setBField (const MagneticField *field)

void	setPtExtractor (const MuonSeedPtExtractor *extractor)

virtual	~MuonSeedFromRecHits ()

Protected Types
typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer	ConstMuonRecHitPointer

typedef MuonTransientTrackingRecHit::MuonRecHitContainer	MuonRecHitContainer

typedef MuonTransientTrackingRecHit::MuonRecHitPointer	MuonRecHitPointer

Protected Attributes
const MagneticField *	theField

const MuonSeedPtExtractor *	thePtExtractor

MuonTransientTrackingRecHit::MuonRecHitContainer	theRhits

Detailed Description

Author: A. Vitelli - INFN Torino; porting R.Bellan - INFN Torino

Generate a seed starting from a list of RecHits

Definition at line 19 of file MuonSeedFromRecHits.h.

Member Typedef Documentation

◆ ConstMuonRecHitPointer

typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer MuonSeedFromRecHits::ConstMuonRecHitPointer

protected

Definition at line 39 of file MuonSeedFromRecHits.h.

◆ MuonRecHitContainer

typedef MuonTransientTrackingRecHit::MuonRecHitContainer MuonSeedFromRecHits::MuonRecHitContainer

protected

Definition at line 37 of file MuonSeedFromRecHits.h.

◆ MuonRecHitPointer

typedef MuonTransientTrackingRecHit::MuonRecHitPointer MuonSeedFromRecHits::MuonRecHitPointer

protected

Definition at line 38 of file MuonSeedFromRecHits.h.

Constructor & Destructor Documentation

◆ MuonSeedFromRecHits()

MuonSeedFromRecHits::MuonSeedFromRecHits ( )

Definition at line 36 of file MuonSeedFromRecHits.cc.

36 : theField(nullptr) {}

◆ ~MuonSeedFromRecHits()

virtual MuonSeedFromRecHits::~MuonSeedFromRecHits ( )

inlinevirtual

Definition at line 22 of file MuonSeedFromRecHits.h.

22 {}

Member Function Documentation

◆ add()

void MuonSeedFromRecHits::add ( MuonTransientTrackingRecHit::MuonRecHitPointer hit )

inline

Definition at line 27 of file MuonSeedFromRecHits.h.

27 { theRhits.push_back(hit); }

References theRhits.

Referenced by counter.Counter::register(), and MuonSeedFinder::seeds().

◆ clear()

void MuonSeedFromRecHits::clear ( void )

inline

Definition at line 30 of file MuonSeedFromRecHits.h.

30 { theRhits.clear(); }

References theRhits.

Referenced by MuonSeedFinder::seeds().

◆ createSeed()

TrajectorySeed MuonSeedFromRecHits::createSeed	(	float	ptmean,
		float	sptmean,
		MuonTransientTrackingRecHit::ConstMuonRecHitPointer	last
	)		const

Definition at line 38 of file MuonSeedFromRecHits.cc.

                                                                                                              {
   const std::string metname = "Muon|RecoMuon|MuonSeedFromRecHits";
  
   MuonPatternRecoDumper debug;
  
   // FIXME: put it into a parameter set!
   double theMinMomentum = 3.0;
   int charge = std::copysign(1, ptmean);
  
   // Minimal pt
   if (fabs(ptmean) < theMinMomentum)
     ptmean = theMinMomentum * charge;
  
   AlgebraicVector t(4);
   AlgebraicSymMatrix mat(5, 0);
  
   // Fill the LocalTrajectoryParameters
   LocalPoint segPos = last->localPosition();
   GlobalVector mom = last->globalPosition() - GlobalPoint();
   GlobalVector polar(GlobalVector::Spherical(mom.theta(), last->globalDirection().phi(), 1.));
   polar *= fabs(ptmean) / polar.perp();
   LocalVector segDirFromPos = last->det()->toLocal(polar);
  
   LocalTrajectoryParameters param(segPos, segDirFromPos, charge);
  
   // this perform H.T() * parErr * H, which is the projection of the
   // the measurement error (rechit rf) to the state error (TSOS rf)
   // Legenda:
   // H => is the 4x5 projection matrix
   // parError the 4x4 parameter error matrix of the RecHit
  
   // LogTrace(metname) << "Projection matrix:\n" << last->projectionMatrix();
   // LogTrace(metname) << "Error matrix:\n" << last->parametersError();
  
   mat = last->parametersError().similarityT(last->projectionMatrix());
  
   float p_err = sqr(sptmean / (ptmean * ptmean));
   mat[0][0] = p_err;
  
   LocalTrajectoryError error(asSMatrix<5>(mat));
  
   // Create the TrajectoryStateOnSurface
   TrajectoryStateOnSurface tsos(param, error, last->det()->surface(), theField);
  
   // The following LogTraces must be moved somewhere else (StandAloneTrajectoryBuilder)
   // Here the TSOS does not have the magnetic field set, so dumpTSOS causes a crash
   // (when LogTrace/LogDebug is activated)
   //LogTrace(metname) << "Trajectory State on Surface before the extrapolation"<<endl;
   //LogTrace(metname) << debug.dumpTSOS(tsos);
  
   // Take the DetLayer on which relies the rechit
   DetId id = last->geographicalId();
   // Segment layer
   //LogTrace(metname) << "The RecSegment relies on: "<<endl;
   //LogTrace(metname) << debug.dumpMuonId(id);
   //LogTrace(metname) << debug.dumpTSOS(tsos);
  
   // Transform it in a TrajectoryStateOnSurface
  
   PTrajectoryStateOnDet const& seedTSOS = trajectoryStateTransform::persistentState(tsos, id.rawId());
  
   edm::OwnVector<TrackingRecHit> container;
   for (unsigned l = 0; l < theRhits.size(); l++) {
     container.push_back(theRhits[l]->hit()->clone());
   }
  
   TrajectorySeed theSeed(seedTSOS, container, alongMomentum);
  
   return theSeed;
 }