#include <MuonTransientTrackingRecHit.h>

Inheritance diagram for MuonTransientTrackingRecHit:

Public Types
typedef std::vector < ConstMuonRecHitPointer >	ConstMuonRecHitContainer

typedef MuonTransientTrackingRecHit const *	ConstMuonRecHitPointer

typedef std::vector < MuonRecHitPointer >	MuonRecHitContainer

typedef MuonTransientTrackingRecHit *	MuonRecHitPointer

Public Types inherited from GenericTransientTrackingRecHit
typedef TrackingRecHit::Type	Type

Public Types inherited from TrackingRecHit
typedef std::vector < ConstRecHitPointer >	ConstRecHitContainer

using	ConstRecHitPointer = std::shared_ptr< TrackingRecHit const >

typedef unsigned int	id_type

typedef std::vector < ConstRecHitPointer >	RecHitContainer

using	RecHitPointer = std::shared_ptr< TrackingRecHit const >

enum	SharedInputType { all, some }
	definition of equality via shared input More...

enum	Type { valid = 0, missing = 1, inactive = 2, bad = 3, missing_inner = 4, missing_outer = 5 }

Public Member Functions
virtual double	chi2 () const
	Chi square of the fit for segments, else 0. More...

virtual int	degreesOfFreedom () const
	Degrees of freedom for segments, else 0. More...

virtual GlobalVector	globalDirection () const
	Direction in 3D for segments, otherwise (0,0,0) More...

virtual GlobalError	globalDirectionError () const
	Error on the global direction. More...

void	invalidateHit ()

bool	isCSC () const
	if this rec hit is a CSC rec hit More...

bool	isDT () const
	if this rec hit is a DT rec hit More...

bool	isRPC () const
	if this rec hit is a RPC rec hit More...

virtual LocalVector	localDirection () const
	Direction in 3D for segments, otherwise (0,0,0) More...

virtual LocalError	localDirectionError () const
	Error on the local direction. More...

virtual AlgebraicSymMatrix	parametersError () const override

virtual ConstRecHitContainer	transientHits () const override
	return the sub components of this transient rechit More...

virtual	~MuonTransientTrackingRecHit ()

Public Member Functions inherited from GenericTransientTrackingRecHit
virtual bool	canImproveWithTrack () const override

TrackingRecHit *	cloneHit () const

virtual int	dimension () const

virtual const TrackingRecHit *	hit () const

virtual LocalPoint	localPosition () const override

virtual LocalError	localPositionError () const override

virtual AlgebraicVector	parameters () const

virtual AlgebraicMatrix	projectionMatrix () const

virtual std::vector< const TrackingRecHit * >	recHits () const override
	Access to component RecHits (if any) More...

virtual std::vector < TrackingRecHit * >	recHits () override
	Non-const access to component RecHits (if any) More...

virtual	~GenericTransientTrackingRecHit ()

Public Member Functions inherited from TValidTrackingRecHit
TrackingRecHit *	cloneHit () const

virtual float	clusterProbability () const
	cluster probability, overloaded by pixel rechits. More...

float	errorGlobalR () const

float	errorGlobalRPhi () const

float	errorGlobalZ () const

virtual GlobalPoint	globalPosition () const

GlobalError	globalPositionError () const

TrackingRecHitGlobalState	globalState () const

virtual const Surface *	surface () const

	TValidTrackingRecHit (const GeomDet &geom)

template<typename... Args>
	TValidTrackingRecHit (Args &&...args)

Public Member Functions inherited from TrackingRecHit
RecHitPointer	cloneForFit (const GeomDet &idet) const

virtual RecHitPointer	cloneSH () const

const GeomDet *	det () const

virtual const GeomDetUnit *	detUnit () const

DetId	geographicalId () const

virtual void	getKfComponents (KfComponentsHolder &holder) const

unsigned int	getRTTI () const

Type	getType () const

virtual bool	hasPositionAndError () const
	to be redefined by daughter class More...

bool	isValid () const

id_type	rawId () const

virtual void	recHitsV (std::vector< const TrackingRecHit * > &) const

virtual void	recHitsV (std::vector< TrackingRecHit * > &)

virtual void	setDet (const GeomDet &idet)

virtual bool	sharesInput (const TrackingRecHit *other, SharedInputType what) const

	TrackingRecHit (DetId id, Type type=valid)

	TrackingRecHit (id_type id=0, Type type=valid)

	TrackingRecHit (DetId id, unsigned int rt, Type type=valid)

	TrackingRecHit (const GeomDet &idet, Type type=valid)

	TrackingRecHit (const GeomDet &idet, unsigned int rt, Type type=valid)

	TrackingRecHit (const GeomDet &idet, TrackingRecHit const &rh)

Type	type () const

virtual float	weight () const

virtual	~TrackingRecHit ()

Static Public Member Functions
static RecHitPointer	build (const GeomDet geom, const TrackingRecHit rh)
	FIXME virtual ConstMuonRecHitContainer specificTransientHits() const;. More...

static MuonRecHitPointer	specificBuild (const GeomDet geom, const TrackingRecHit rh)

Static Public Member Functions inherited from GenericTransientTrackingRecHit
static RecHitPointer	build (const GeomDet geom, const TrackingRecHit rh)

Private Member Functions
virtual MuonTransientTrackingRecHit *	clone () const override

	MuonTransientTrackingRecHit (const GeomDet geom, const TrackingRecHit rh)
	Construct from a TrackingRecHit and its GeomDet. More...

	MuonTransientTrackingRecHit (const MuonTransientTrackingRecHit &other)
	Copy ctor. More...

Friends
class	kkkwwwxxxyyyzzz

Additional Inherited Members
Static Public Attributes inherited from TrackingRecHit
static const int	rttiShift = 24

static const int	typeMask = 0xf

Protected Member Functions inherited from GenericTransientTrackingRecHit
	GenericTransientTrackingRecHit (const GeomDet &geom, const TrackingRecHit &rh)

	GenericTransientTrackingRecHit (const GeomDet &geom, TrackingRecHit *rh)
	for derived classes convenience, does not clone! More...

	GenericTransientTrackingRecHit (const GenericTransientTrackingRecHit &other)

Protected Member Functions inherited from TrackingRecHit
void	setId (id_type iid)

void	setRTTI (unsigned int rt)

void	setType (Type ttype)

Protected Attributes inherited from GenericTransientTrackingRecHit
TrackingRecHit *	trackingRecHit_

Detailed Description

A TransientTrackingRecHit for muons.

Author: C. Liu Purdue University

Definition at line 17 of file MuonTransientTrackingRecHit.h.

Member Typedef Documentation

typedef std::vector<ConstMuonRecHitPointer> MuonTransientTrackingRecHit::ConstMuonRecHitContainer

Definition at line 31 of file MuonTransientTrackingRecHit.h.

typedef MuonTransientTrackingRecHit const* MuonTransientTrackingRecHit::ConstMuonRecHitPointer

Definition at line 25 of file MuonTransientTrackingRecHit.h.

typedef std::vector<MuonRecHitPointer> MuonTransientTrackingRecHit::MuonRecHitContainer

Definition at line 30 of file MuonTransientTrackingRecHit.h.

typedef MuonTransientTrackingRecHit* MuonTransientTrackingRecHit::MuonRecHitPointer

Definition at line 24 of file MuonTransientTrackingRecHit.h.

Constructor & Destructor Documentation

virtual MuonTransientTrackingRecHit::~MuonTransientTrackingRecHit ( )

inlinevirtual

Definition at line 33 of file MuonTransientTrackingRecHit.h.

33 {}

MuonTransientTrackingRecHit::MuonTransientTrackingRecHit	(	const GeomDet *	geom,
		const TrackingRecHit *	rh
	)

private

Construct from a TrackingRecHit and its GeomDet.

Definition at line 21 of file MuonTransientTrackingRecHit.cc.

Referenced by build(), clone(), and specificBuild().

21 :

22 GenericTransientTrackingRecHit(*geom,*rh){}

GenericTransientTrackingRecHit::GenericTransientTrackingRecHit

GenericTransientTrackingRecHit(const GeomDet &geom, const TrackingRecHit &rh)

Definition: GenericTransientTrackingRecHit.h:44

MuonTransientTrackingRecHit::MuonTransientTrackingRecHit ( const MuonTransientTrackingRecHit & other )

private

Copy ctor.

Definition at line 24 of file MuonTransientTrackingRecHit.cc.

24 :

25 GenericTransientTrackingRecHit(*other.det(), *(other.hit())) {}

GenericTransientTrackingRecHit::GenericTransientTrackingRecHit

GenericTransientTrackingRecHit(const GeomDet &geom, const TrackingRecHit &rh)

Definition: GenericTransientTrackingRecHit.h:44

GenericTransientTrackingRecHit::hit

virtual const TrackingRecHit * hit() const

Definition: GenericTransientTrackingRecHit.h:27

TrackingRecHit::det

const GeomDet * det() const

Definition: TrackingRecHit.h:118

Member Function Documentation

static RecHitPointer MuonTransientTrackingRecHit::build	(	const GeomDet *	geom,
		const TrackingRecHit *	rh
	)

inlinestatic

FIXME virtual ConstMuonRecHitContainer specificTransientHits() const;.

Definition at line 69 of file MuonTransientTrackingRecHit.h.

References MuonTransientTrackingRecHit().

Referenced by GlobalTrajectoryBuilderBase::fixTEC().

                                                                               {
     return RecHitPointer( new MuonTransientTrackingRecHit(geom, rh));
   }

double MuonTransientTrackingRecHit::chi2 ( void ) const

virtual

Chi square of the fit for segments, else 0.

Definition at line 126 of file MuonTransientTrackingRecHit.cc.

References GenericTransientTrackingRecHit::hit().

Referenced by CosmicMuonSeedGenerator::checkQuality(), and MuonSeedCleaner::NChi2OfSegment().

 {
   if (dynamic_cast<const RecSegment*>(hit()))
     return dynamic_cast<const RecSegment*>(hit())->chi2();
   else return 0.;
 }

virtual MuonTransientTrackingRecHit* MuonTransientTrackingRecHit::clone ( void ) const

inlineoverrideprivatevirtual

Reimplemented from GenericTransientTrackingRecHit.

Definition at line 89 of file MuonTransientTrackingRecHit.h.

References MuonTransientTrackingRecHit().

Referenced by transientHits().

                                                                {
     return new MuonTransientTrackingRecHit(*this);
   }

int MuonTransientTrackingRecHit::degreesOfFreedom ( ) const

virtual

Degrees of freedom for segments, else 0.

Definition at line 133 of file MuonTransientTrackingRecHit.cc.

References GenericTransientTrackingRecHit::hit().

Referenced by CosmicMuonSeedGenerator::leftIsBetter(), and MuonSeedCleaner::NChi2OfSegment().

 {
   if (dynamic_cast<const RecSegment*>(hit()))
     return dynamic_cast<const RecSegment*>(hit())->degreesOfFreedom();
   else return 0;
 }

GlobalVector MuonTransientTrackingRecHit::globalDirection ( ) const

virtual

Direction in 3D for segments, otherwise (0,0,0)

Definition at line 44 of file MuonTransientTrackingRecHit.cc.

References TrackingRecHit::det(), localDirection(), and TValidTrackingRecHit::surface().

Referenced by CosmicMuonSeedGenerator::areCorrelated(), CosmicMuonSeedGenerator::createSeed(), and MuonSeedOrcaPatternRecognition::discriminator().

 {
   return  (det()->surface().toGlobal(localDirection()));
 }

GlobalError MuonTransientTrackingRecHit::globalDirectionError ( ) const

virtual

Error on the global direction.

Definition at line 49 of file MuonTransientTrackingRecHit.cc.

References TrackingRecHit::det(), localDirectionError(), TValidTrackingRecHit::surface(), and ErrorFrameTransformer::transform().

 {
   return ErrorFrameTransformer().transform( localDirectionError(), (det()->surface()));
 }

void MuonTransientTrackingRecHit::invalidateHit ( )

Definition at line 214 of file MuonTransientTrackingRecHit.cc.

References TrackingRecHit::bad, GenericTransientTrackingRecHit::dimension(), isCSC(), isDT(), GenericTransientTrackingRecHit::recHits(), TrackingRecHit::setType(), and GenericTransientTrackingRecHit::trackingRecHit_.

Referenced by StandAloneMuonFilter::createDefaultTrajectory(), and MuonTrajectoryUpdator::update().

                                                { 
   setType(bad); trackingRecHit_->setType(bad); 
 
 
   if (isDT()){
     if(dimension() > 1){ // MB4s have 2D, but formatted in 4D segments 
       std::vector<TrackingRecHit*> seg2D = recHits(); // 4D --> 2D
       // load 1D hits (2D --> 1D)
       for(std::vector<TrackingRecHit*>::iterator it = seg2D.begin(); it != seg2D.end(); ++it){
     std::vector<TrackingRecHit*> hits1D =  (*it)->recHits();
     (*it)->setType(bad);
     for(std::vector<TrackingRecHit*>::iterator it2 = hits1D.begin(); it2 != hits1D.end(); ++it2)
       (*it2)->setType(bad);
       }
     }
   }
   else if(isCSC())
     if(dimension() == 4){
       std::vector<TrackingRecHit*>  hits = recHits(); // load 2D hits (4D --> 1D)
       for(std::vector<TrackingRecHit*>::iterator it = hits.begin(); it != hits.end(); ++it)
     (*it)->setType(bad);
     }
   
   
 }

bool MuonTransientTrackingRecHit::isCSC ( ) const

if this rec hit is a CSC rec hit

Definition at line 144 of file MuonTransientTrackingRecHit.cc.

References MuonSubdetId::CSC, and TrackingRecHit::geographicalId().

Referenced by CosmicMuonSeedGenerator::checkQuality(), and invalidateHit().

                                              {
   return  (geographicalId().subdetId() == MuonSubdetId::CSC);
 }

bool MuonTransientTrackingRecHit::isDT ( ) const

if this rec hit is a DT rec hit

Definition at line 140 of file MuonTransientTrackingRecHit.cc.

References MuonSubdetId::DT, and TrackingRecHit::geographicalId().

Referenced by CosmicMuonSeedGenerator::checkQuality(), MuonSeedOrcaPatternRecognition::discriminator(), invalidateHit(), and transientHits().

                                             {
   return  (geographicalId().subdetId() == MuonSubdetId::DT);
 }

bool MuonTransientTrackingRecHit::isRPC ( ) const

if this rec hit is a RPC rec hit

Definition at line 148 of file MuonTransientTrackingRecHit.cc.

References TrackingRecHit::geographicalId(), and MuonSubdetId::RPC.

                                              {
   return  (geographicalId().subdetId() == MuonSubdetId::RPC);
 }

LocalVector MuonTransientTrackingRecHit::localDirection ( ) const

virtual

Direction in 3D for segments, otherwise (0,0,0)

Definition at line 28 of file MuonTransientTrackingRecHit.cc.

References GenericTransientTrackingRecHit::hit().

Referenced by globalDirection().

                                                               {
 
   if (dynamic_cast<const RecSegment*>(hit()) )
      return dynamic_cast<const RecSegment*>(hit())->localDirection(); 
   else return LocalVector(0.,0.,0.);
 
 }

LocalError MuonTransientTrackingRecHit::localDirectionError ( ) const

virtual

Error on the local direction.

Definition at line 36 of file MuonTransientTrackingRecHit.cc.

References GenericTransientTrackingRecHit::hit().

Referenced by globalDirectionError().

                                                                   {
 
   if (dynamic_cast<const RecSegment*>(hit()))
      return dynamic_cast<const RecSegment*>(hit())->localDirectionError();
   else return LocalError(0.,0.,0.);
 
 }

AlgebraicSymMatrix MuonTransientTrackingRecHit::parametersError ( ) const

overridevirtual

Reimplemented from GenericTransientTrackingRecHit.

Definition at line 55 of file MuonTransientTrackingRecHit.cc.

Referenced by CosmicMuonSeedGenerator::createSeed().

                                                                       {
  
   AlgebraicSymMatrix err = GenericTransientTrackingRecHit::parametersError();
   AlgebraicVector par = GenericTransientTrackingRecHit::parameters();
 
   const AlignmentPositionError* APE = det()->alignmentPositionError();
   if (APE != NULL) {
     AlgebraicVector positions(2,0);
     AlgebraicVector directions(2,0);
 
     if(err.num_row() == 1) {
      positions[0] = 0.;
      positions[1] = 0.;
      directions[0] = 0.;
      directions[1] = 0.;
      LocalErrorExtended lape = ErrorFrameTransformer().transform46(APE->globalError(),positions,directions);
      err[0][0] += lape.cxx();
     } else if (err.num_row() == 2) {
      positions[0] = localPosition().x();
      positions[1] = 0.;
      directions[0] = 0.;
      directions[1] = 0.;
      LocalErrorExtended lape = ErrorFrameTransformer().transform46(APE->globalError(),positions,directions);
 
      AlgebraicSymMatrix lapeMatrix(2,0);
      lapeMatrix[1][1] = lape.cxx();
      lapeMatrix[0][0] = lape.cphixphix();
      lapeMatrix[0][1] = lape.cphixx();
 
      if(err.num_row() != lapeMatrix.num_row())
        throw cms::Exception("MuonTransientTrackingRecHit::parametersError")
         <<"Discrepancy between alignment error matrix and error matrix: APE "
         << lapeMatrix.num_row()
         << ", error matrix " << err.num_row()
         << std::endl;
 
      err += lapeMatrix;
     } else if (err.num_row() == 4) { 
      positions[0] = par[2];
      positions[1] = par[3];
      directions[0] = par[0];
      directions[1] = par[1];
 
      LocalErrorExtended lape = ErrorFrameTransformer().transform46(APE->globalError(),positions,directions);
 
      AlgebraicSymMatrix lapeMatrix(4,0);
      lapeMatrix[2][2] = lape.cxx();
      lapeMatrix[2][3] = lape.cyx();
      lapeMatrix[3][3] = lape.cyy();
      lapeMatrix[0][0] = lape.cphixphix();
      lapeMatrix[0][1] = lape.cphiyphix();
      lapeMatrix[1][1] = lape.cphiyphiy();
 
      lapeMatrix[0][2] = lape.cphixx();
      lapeMatrix[0][3] = lape.cphixy();
      lapeMatrix[1][3] = lape.cphiyy();
      lapeMatrix[1][2] = lape.cphiyx();
 
      if(err.num_row() != lapeMatrix.num_row())
        throw cms::Exception("MuonTransientTrackingRecHit::parametersError")
         <<"Discrepancy between alignment error matrix and error matrix: APE "
         << lapeMatrix.num_row()
         << ", error matrix " << err.num_row()
         << std::endl;
 
      err += lapeMatrix;
     }
   }
   return err;
 }

static MuonRecHitPointer MuonTransientTrackingRecHit::specificBuild	(	const GeomDet *	geom,
		const TrackingRecHit *	rh
	)

inlinestatic

Definition at line 73 of file MuonTransientTrackingRecHit.h.

References MuonTransientTrackingRecHit().

Referenced by SegmentToTrackAssociator::associate(), SegmentsTrackAssociator::associate(), MuonTransientTrackingRecHitBuilder::build(), ResidualRefitting::CollectTrackHits(), StandAloneMuonFilter::createDefaultTrajectory(), MuonShowerInformationFiller::fillHitsByStation(), MuonShowerInformationFiller::hitsFromSegments(), MuonSeedCleaner::NChi2OfSegment(), ResidualRefitting::NewTrackMeasurements(), MuonSeedCleaner::NRecHitsFromSegment(), SETPatternRecognition::produce(), MuonMillepedeTrackRefitter::produce(), MuonDetLayerMeasurements::recHits(), ExhaustiveMuonTrajectoryBuilder::trajectories(), and MuonTrajectoryUpdator::update().

                                                                                          {
     return MuonRecHitPointer(new MuonTransientTrackingRecHit(geom, rh));
   }

TransientTrackingRecHit::ConstRecHitContainer MuonTransientTrackingRecHit::transientHits ( ) const

overridevirtual

return the sub components of this transient rechit

Reimplemented from TrackingRecHit.

Definition at line 154 of file MuonTransientTrackingRecHit.cc.

References clone(), GeomDet::components(), filterCSVwithJSON::copy, TrackingRecHit::det(), GenericTransientTrackingRecHit::dimension(), GeomDet::geographicalId(), isDT(), GenericTransientTrackingRecHit::recHits(), and tmp.

                                                                                             {
 
   ConstRecHitContainer theSubTransientRecHits;
   
   // the sub rec hit of this TransientRecHit
   std::vector<const TrackingRecHit*> ownRecHits = recHits();
 
   if(ownRecHits.size() == 0){
     theSubTransientRecHits.push_back(TransientTrackingRecHit::RecHitPointer(clone()));
     return theSubTransientRecHits;
   }
   
   // the components of the geom det on which reside this rechit
   std::vector<const GeomDet *> geomDets = det()->components();
 
   if(isDT() && dimension() == 2 && ownRecHits.front()->dimension() == 1 
      && (geomDets.size() == 3 || geomDets.size() == 2) ){ // it is a phi segment!!
     
     std::vector<const GeomDet *> subGeomDets;
 
     int sl = 1;
     for(std::vector<const GeomDet *>::const_iterator geoDet = geomDets.begin();
     geoDet != geomDets.end(); ++geoDet){
       if(sl != 3){ // FIXME!! this maybe is not always true
     std::vector<const GeomDet *> tmp = (*geoDet)->components();
     std::copy(tmp.begin(),tmp.end(),back_inserter(subGeomDets));
       }
       ++sl;
     }
     geomDets.clear();
     geomDets = subGeomDets;
   }
   
   // Fill the GeomDet map
   std::map<DetId,const GeomDet*> gemDetMap;
   
   for (std::vector<const GeomDet*>::const_iterator subDet = geomDets.begin(); 
        subDet != geomDets.end(); ++subDet)
     gemDetMap[ (*subDet)->geographicalId() ] = *subDet;
   
   std::map<DetId,const GeomDet*>::iterator gemDetMap_iter;
   
   // Loop in order to check the ids
   for (std::vector<const TrackingRecHit*>::const_iterator rechit = ownRecHits.begin(); 
        rechit != ownRecHits.end(); ++rechit){
     
     gemDetMap_iter = gemDetMap.find( (*rechit)->geographicalId() );
     
     if(gemDetMap_iter != gemDetMap.end() )
       theSubTransientRecHits.push_back(TransientTrackingRecHit::RecHitPointer(new MuonTransientTrackingRecHit(gemDetMap_iter->second, 
                                      *rechit)) );
     else if( (*rechit)->geographicalId() == det()->geographicalId() ) // Phi in DT is on Chamber
       theSubTransientRecHits.push_back(TransientTrackingRecHit::RecHitPointer(new MuonTransientTrackingRecHit(det(), 
                                        *rechit)) );
   }
   return theSubTransientRecHits;
 
 }

Friends And Related Function Documentation

friend class kkkwwwxxxyyyzzz

friend

Definition at line 81 of file MuonTransientTrackingRecHit.h.

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Friends

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation