MuonTransientTrackingRecHit Class Reference

#include <MuonTransientTrackingRecHit.h>

Inheritance diagram for MuonTransientTrackingRecHit:

Public Types
typedef std::vector < ConstMuonRecHitPointer >	ConstMuonRecHitContainer
using	ConstMuonRecHitPointer = std::shared_ptr< MuonTransientTrackingRecHit const >
typedef std::vector < MuonRecHitPointer >	MuonRecHitContainer
using	MuonRecHitPointer = std::shared_ptr< MuonTransientTrackingRecHit >
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	isGEM () const
	if this rec hit is a GEM rec hit More...
bool	isME0 () const
	if this rec hit is a ME0 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 final
float	errorGlobalRPhi () const final
float	errorGlobalZ () const final
virtual GlobalPoint	globalPosition () const final
GlobalError	globalPositionError () const final
TrackingRecHitGlobalState	globalState () const
virtual const Surface *	surface () const final
	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

by C. Calabria INFN & Universita  Bari

Definition at line 21 of file MuonTransientTrackingRecHit.h.

Member Typedef Documentation

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

Definition at line 29 of file MuonTransientTrackingRecHit.h.

using MuonTransientTrackingRecHit::ConstMuonRecHitPointer = std::shared_ptr<MuonTransientTrackingRecHit const>

Definition at line 24 of file MuonTransientTrackingRecHit.h.

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

Definition at line 28 of file MuonTransientTrackingRecHit.h.

using MuonTransientTrackingRecHit::MuonRecHitPointer = std::shared_ptr<MuonTransientTrackingRecHit>

Definition at line 23 of file MuonTransientTrackingRecHit.h.

Constructor & Destructor Documentation

virtual MuonTransientTrackingRecHit::~MuonTransientTrackingRecHit ( )

inlinevirtual

Definition at line 31 of file MuonTransientTrackingRecHit.h.

{}

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().

                                                                                                      :
  GenericTransientTrackingRecHit(*geom,*rh){}

MuonTransientTrackingRecHit::MuonTransientTrackingRecHit ( const MuonTransientTrackingRecHit &  other )

private

Copy ctor.

Definition at line 24 of file MuonTransientTrackingRecHit.cc.

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

Member Function Documentation

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

inlinestatic

FIXME virtual ConstMuonRecHitContainer specificTransientHits() const;.

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

{
  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 94 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().

{
  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().

{
  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 223 of file MuonTransientTrackingRecHit.cc.

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

                                               { 
  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 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 invalidateHit(), and transientHits().

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

bool MuonTransientTrackingRecHit::isGEM

(

)

const

if this rec hit is a GEM rec hit

Definition at line 148 of file MuonTransientTrackingRecHit.cc.

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

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

bool MuonTransientTrackingRecHit::isME0

(

)

const

if this rec hit is a ME0 rec hit

Definition at line 152 of file MuonTransientTrackingRecHit.cc.

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

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

bool MuonTransientTrackingRecHit::isRPC

(

)

const

if this rec hit is a RPC rec hit

Definition at line 157 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.

References GeomDet::alignmentPositionError(), LocalErrorBaseExtended< T, ErrorWeightType >::cphixphix(), LocalErrorBaseExtended< T, ErrorWeightType >::cphixx(), LocalErrorBaseExtended< T, ErrorWeightType >::cphixy(), LocalErrorBaseExtended< T, ErrorWeightType >::cphiyphix(), LocalErrorBaseExtended< T, ErrorWeightType >::cphiyphiy(), LocalErrorBaseExtended< T, ErrorWeightType >::cphiyx(), LocalErrorBaseExtended< T, ErrorWeightType >::cphiyy(), LocalErrorBaseExtended< T, ErrorWeightType >::cxx(), LocalErrorBaseExtended< T, ErrorWeightType >::cyx(), LocalErrorBaseExtended< T, ErrorWeightType >::cyy(), TrackingRecHit::det(), AlignmentPositionError::globalError(), GenericTransientTrackingRecHit::localPosition(), NULL, GenericTransientTrackingRecHit::parameters(), GenericTransientTrackingRecHit::parametersError(), ErrorFrameTransformer::transform46(), and PV3DBase< T, PVType, FrameType >::x().

                                                                      {
 
  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 77 of file MuonTransientTrackingRecHit.h.

References LogDebug, and 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().

                                                                                         {
    LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "Getting specificBuild"<<std::endl;
    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 163 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 86 of file MuonTransientTrackingRecHit.h.