#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

typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer MuonSeedFromRecHits::ConstMuonRecHitPointer

protected

Definition at line 39 of file MuonSeedFromRecHits.h.

typedef MuonTransientTrackingRecHit::MuonRecHitContainer MuonSeedFromRecHits::MuonRecHitContainer

protected

Definition at line 37 of file MuonSeedFromRecHits.h.

typedef MuonTransientTrackingRecHit::MuonRecHitPointer MuonSeedFromRecHits::MuonRecHitPointer

protected

Definition at line 38 of file MuonSeedFromRecHits.h.

Constructor & Destructor Documentation

MuonSeedFromRecHits::MuonSeedFromRecHits ( )

Definition at line 33 of file MuonSeedFromRecHits.cc.

 : theField(0)
 {
 }

virtual MuonSeedFromRecHits::~MuonSeedFromRecHits ( )

inlinevirtual

Definition at line 23 of file MuonSeedFromRecHits.h.

23 {}

Member Function Documentation

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

inline

Definition at line 28 of file MuonSeedFromRecHits.h.

References theRhits.

Referenced by MuonSeedFinder::seeds().

28 { theRhits.push_back(hit); }

MuonSeedFromRecHits::theRhits

MuonTransientTrackingRecHit::MuonRecHitContainer theRhits

Definition: MuonSeedFromRecHits.h:41

void MuonSeedFromRecHits::clear ( void )

inline

Definition at line 31 of file MuonSeedFromRecHits.h.

References theRhits.

Referenced by MuonSeedFinder::seeds().

31 {theRhits.clear();}

MuonSeedFromRecHits::theRhits

MuonTransientTrackingRecHit::MuonRecHitContainer theRhits

Definition: MuonSeedFromRecHits.h:41

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

Definition at line 39 of file MuonSeedFromRecHits.cc.

References alongMomentum, DeDxDiscriminatorTools::charge(), clone(), debug, relativeConstraints::error, prof2calltree::l, metname, PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), edm::OwnVector< T, P >::push_back(), funct::sqr(), AlCaHLTBitMon_QueryRunRegistry::string, edmStreamStallGrapher::t, theField, theRhits, and PV3DBase< T, PVType, FrameType >::theta().

Referenced by MuonCSCSeedFromRecHits::createDefaultEndcapSeed(), MuonOverlapSeedFromRecHits::makeSeed(), MuonCSCSeedFromRecHits::makeSeed(), MuonDTSeedFromRecHits::seed(), and ExhaustiveMuonTrajectoryBuilder::trajectories().

 {
   
   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;
 }