SiTrackerMultiRecHitUpdator Class Reference

#include <SiTrackerMultiRecHitUpdator.h>

Public Types
typedef std::pair< LocalPoint, LocalError >	LocalParameters

Public Member Functions
const std::vector< double > &	annealingProgram () const
virtual TransientTrackingRecHit::RecHitPointer	buildMultiRecHit (const std::vector< const TrackingRecHit * > &rhv, const TrajectoryStateOnSurface &tsos, MeasurementDetWithData &measDet, float annealing=1.) const
double	ComputeWeight (const TrajectoryStateOnSurface &tsos, const TransientTrackingRecHit &aRecHit, bool CutWeight, double annealing=1.) const
template<unsigned int N>
double	ComputeWeight (const TrajectoryStateOnSurface &tsos, const TransientTrackingRecHit &aRecHit, bool CutWeight, double annealing=1.) const
const std::vector< double > &	getAnnealingProgram () const
const TransientTrackingRecHitBuilder *	getBuilder () const
	SiTrackerMultiRecHitUpdator (const TransientTrackingRecHitBuilder *builder, const TrackingRecHitPropagator *hitpropagator, const float Chi2Cut1D, const float Chi2Cut2D, const std::vector< double > &anAnnealingProgram, bool debug)
virtual TransientTrackingRecHit::RecHitPointer	update (TransientTrackingRecHit::ConstRecHitPointer original, const TrajectoryStateOnSurface &tsos, MeasurementDetWithData &measDet, double annealing=1.) const
TransientTrackingRecHit::RecHitPointer	update (TransientTrackingRecHit::ConstRecHitContainer &tcomponents, const TrajectoryStateOnSurface &tsos, MeasurementDetWithData &measDet, double annealing=1.) const
virtual	~SiTrackerMultiRecHitUpdator ()

Private Member Functions
LocalParameters	calcParameters (const TrajectoryStateOnSurface &tsos, std::vector< std::pair< const TrackingRecHit *, float > > &aHitMap) const
template<unsigned int N>
LocalParameters	calcParameters (const TrajectoryStateOnSurface &tsos, std::vector< std::pair< const TrackingRecHit *, float > > &aHitMap) const
bool	TIDorTEChit (const TrackingRecHit *const &hit) const

Private Attributes
bool	debug_
const std::vector< double >	theAnnealingProgram
const TransientTrackingRecHitBuilder *	theBuilder
double	theChi2Cut1D
double	theChi2Cut2D
TkClonerImpl	theHitCloner
const TrackingRecHitPropagator *	theHitPropagator

Detailed Description

Builds a SiTrackerMultiRecHit out of a vector of TrackingRecHit or updates an existing SiTrackerMultiRecHit given a tsos.

Author

tropiano, genta in May 2014 by brondolin

Definition at line 27 of file SiTrackerMultiRecHitUpdator.h.

Member Typedef Documentation

typedef std::pair<LocalPoint,LocalError> SiTrackerMultiRecHitUpdator::LocalParameters

Definition at line 31 of file SiTrackerMultiRecHitUpdator.h.

Constructor & Destructor Documentation

SiTrackerMultiRecHitUpdator::SiTrackerMultiRecHitUpdator	(	const TransientTrackingRecHitBuilder *	builder,
		const TrackingRecHitPropagator *	hitpropagator,
		const float	Chi2Cut1D,
		const float	Chi2Cut2D,
		const std::vector< double > &	anAnnealingProgram,
		bool	debug
	)

Definition at line 19 of file SiTrackerMultiRecHitUpdator.cc.

References HLT_25ns10e33_v2_cff::cloner, and theHitCloner.

                                        :
  theBuilder(builder),
  theHitPropagator(hitpropagator),
  theChi2Cut1D(Chi2Cut1D),
  theChi2Cut2D(Chi2Cut2D),
  theAnnealingProgram(anAnnealingProgram),
  debug_(debug){
    theHitCloner = static_cast<TkTransientTrackingRecHitBuilder const *>(builder)->cloner();
  }

virtual SiTrackerMultiRecHitUpdator::~SiTrackerMultiRecHitUpdator ( )

inlinevirtual

Definition at line 37 of file SiTrackerMultiRecHitUpdator.h.

{};

Member Function Documentation

const std::vector<double>& SiTrackerMultiRecHitUpdator::annealingProgram ( ) const

inline

Definition at line 65 of file SiTrackerMultiRecHitUpdator.h.

References theAnnealingProgram.

{return theAnnealingProgram;}

TransientTrackingRecHit::RecHitPointer SiTrackerMultiRecHitUpdator::buildMultiRecHit ( const std::vector< const TrackingRecHit * > &  rhv, const TrajectoryStateOnSurface &  tsos, MeasurementDetWithData &  measDet, float  annealing = 1.  ) const

virtual

Definition at line 35 of file SiTrackerMultiRecHitUpdator.cc.

References TransientTrackingRecHitBuilder::build(), LogTrace, theBuilder, and update().

                                                                                                     {

  LogTrace("SiTrackerMultiRecHitUpdator") << "Calling SiTrackerMultiRecHitUpdator::buildMultiRecHit with AnnealingFactor: "  << annealing;

  TransientTrackingRecHit::ConstRecHitContainer tcomponents;    
  for (std::vector<const TrackingRecHit*>::const_iterator iter = rhv.begin(); iter != rhv.end(); iter++){

    TransientTrackingRecHit::RecHitPointer transient = theBuilder->build(*iter);
    if(transient->isValid()) tcomponents.push_back(transient);

  }
  return update(tcomponents, tsos, measDet, annealing); 
  
}

SiTrackerMultiRecHitUpdator::LocalParameters SiTrackerMultiRecHitUpdator::calcParameters ( const TrajectoryStateOnSurface &  tsos, std::vector< std::pair< const TrackingRecHit *, float > > &  aHitMap  ) const

private

Definition at line 291 of file SiTrackerMultiRecHitUpdator.cc.

References Exception, plotBeamSpotDB::first, LogTrace, and mergeVDriftHistosByStation::name.

Referenced by update().

                                                                                                                                                                                            {

  //supposing all the hits inside of a MRH have the same dimension
  LogTrace("SiTrackerMultiRecHitUpdator")<<"SiTrackerMultiRecHitUpdator::LocalParameters: dim first recHit: " << aHitMap[0].first->dimension() <<std::endl;
  switch (aHitMap[0].first->dimension()) {
    case 1: return calcParameters<1>(tsos,aHitMap);
    case 2: return calcParameters<2>(tsos,aHitMap);
  }
  throw cms::Exception("Rec hit of invalid dimension for computing MRH (not 1,2)") <<
        "The value was " << aHitMap[0].first->dimension() <<
        ", type is " << typeid(aHitMap[0].first).name() << "\n";

}

template<unsigned int N>

LocalParameters SiTrackerMultiRecHitUpdator::calcParameters ( const TrajectoryStateOnSurface &  tsos, std::vector< std::pair< const TrackingRecHit *, float > > &  aHitMap  ) const

private

double SiTrackerMultiRecHitUpdator::ComputeWeight	(	const TrajectoryStateOnSurface &	tsos,
		const TransientTrackingRecHit &	aRecHit,
		bool	CutWeight,
		double	annealing = 1.
	)		const

Definition at line 191 of file SiTrackerMultiRecHitUpdator.cc.

References TrackingRecHit::dimension(), Exception, and mergeVDriftHistosByStation::name.

Referenced by update().

                                                                                                           {
     switch (aRecHit.dimension()) {
         case 1: return ComputeWeight<1>(tsos,aRecHit,CutWeight,annealing);
         case 2: return ComputeWeight<2>(tsos,aRecHit,CutWeight,annealing);
         case 3: return ComputeWeight<3>(tsos,aRecHit,CutWeight,annealing);
         case 4: return ComputeWeight<4>(tsos,aRecHit,CutWeight,annealing);
         case 5: return ComputeWeight<5>(tsos,aRecHit,CutWeight,annealing);
     }
     throw cms::Exception("Rec hit of invalid dimension (not 1,2,3,4,5)") <<
          "The value was " << aRecHit.dimension() << 
         ", type is " << typeid(aRecHit).name() << "\n";
}

template<unsigned int N>

double SiTrackerMultiRecHitUpdator::ComputeWeight	(	const TrajectoryStateOnSurface &	tsos,
		const TransientTrackingRecHit &	aRecHit,
		bool	CutWeight,
		double	annealing = 1.
	)		const

const std::vector<double>& SiTrackerMultiRecHitUpdator::getAnnealingProgram ( ) const

inline

Definition at line 66 of file SiTrackerMultiRecHitUpdator.h.

References theAnnealingProgram.

Referenced by DAFTrackProducerAlgorithm::runWithCandidate().

{return theAnnealingProgram;}

const TransientTrackingRecHitBuilder* SiTrackerMultiRecHitUpdator::getBuilder ( ) const

inline

Definition at line 68 of file SiTrackerMultiRecHitUpdator.h.

References theBuilder.

Referenced by SimpleDAFHitCollector::SimpleDAFHitCollector().

{return theBuilder;}

bool SiTrackerMultiRecHitUpdator::TIDorTEChit ( const TrackingRecHit *const &  hit ) const

private

Definition at line 380 of file SiTrackerMultiRecHitUpdator.cc.

References DetId::det(), TrackingRecHit::geographicalId(), DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TID, and DetId::Tracker.

                                                                                   {

  DetId hitId = hit->geographicalId();

  if( hitId.det() == DetId::Tracker && 
      ( hitId.subdetId() == StripSubdetector::TEC || hitId.subdetId() == StripSubdetector::TID) ) {
    return true;
  }    

  return false;
}

TransientTrackingRecHit::RecHitPointer SiTrackerMultiRecHitUpdator::update ( TransientTrackingRecHit::ConstRecHitPointer  original, const TrajectoryStateOnSurface &  tsos, MeasurementDetWithData &  measDet, double  annealing = 1.  ) const

virtual

Definition at line 52 of file SiTrackerMultiRecHitUpdator.cc.

References Exception, TrajectoryStateOnSurface::isValid(), LogTrace, TkClonerImpl::makeShared(), and theHitCloner.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), buildMultiRecHit(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), Vispa.Gui.PortConnection.PointToPointConnection::updateConnection(), and DAFTrackProducerAlgorithm::updateHits().

                                                                {

  LogTrace("SiTrackerMultiRecHitUpdator") << "Calling SiTrackerMultiRecHitUpdator::update with AnnealingFactor: "  << annealing;

  if(!tsos.isValid()) {
    //return original->clone();
    throw cms::Exception("SiTrackerMultiRecHitUpdator") << "!!! MultiRecHitUpdator::update(..): tsos NOT valid!!! ";
  } 

  //check if to clone is the right thing
  if(original->isValid()){
    if (original->transientHits().empty()){
      return theHitCloner.makeShared(original,tsos);
    }
  } else {
    return theHitCloner.makeShared(original,tsos);
  }

  TransientTrackingRecHit::ConstRecHitContainer tcomponents = original->transientHits();    
  return update(tcomponents, tsos, measDet, annealing);
}

TransientTrackingRecHit::RecHitPointer SiTrackerMultiRecHitUpdator::update	(	TransientTrackingRecHit::ConstRecHitContainer &	tcomponents,
		const TrajectoryStateOnSurface &	tsos,
		MeasurementDetWithData &	measDet,
		double	annealing = 1.
	)		const

Definition at line 77 of file SiTrackerMultiRecHitUpdator.cc.

References calcParameters(), cropTnPTrees::cloned, ComputeWeight(), counter, alignCSCRings::e, Exception, plotBeamSpotDB::first, MeasurementDet::geomDet(), align::invalid, TrajectoryStateOnSurface::isValid(), LogTrace, TkClonerImpl::makeShared(), MeasurementDetWithData::mdet(), TrackingRecHit::missing, mergeVDriftHistosByStation::name, lumiCalc2::normmap, AlCaHLTBitMon_ParallelJobs::p, GloballyPositioned< T >::position(), TrackingRecHitPropagator::project(), edm::second(), TrajectoryStateOnSurface::surface(), theBuilder, theHitCloner, and theHitPropagator.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                                                                {

  if (tcomponents.empty()){
    LogTrace("SiTrackerMultiRecHitUpdator") << "Empty components vector passed to SiTrackerMultiRecHitUpdator::update, returning an InvalidTransientRecHit ";
    return std::make_shared<InvalidTrackingRecHit>(measDet.mdet().geomDet(), TrackingRecHit::missing);
  }     
  
  if(!tsos.isValid()) {
    LogTrace("SiTrackerMultiRecHitUpdator")<<"SiTrackerMultiRecHitUpdator::update: tsos NOT valid!!!, returning an InvalidTransientRecHit";
    return std::make_shared<InvalidTrackingRecHit>(measDet.mdet().geomDet(), TrackingRecHit::missing);
  }
  
  std::vector<TransientTrackingRecHit::RecHitPointer> updatedcomponents;
  const GeomDet* geomdet = 0;

  //running on all over the MRH components 
  for (TransientTrackingRecHit::ConstRecHitContainer::const_iterator iter = tcomponents.begin(); iter != tcomponents.end(); iter++){

    //the first rechit must belong to the same surface of TSOS
    if (iter == tcomponents.begin()) {

      if (&((*iter)->det()->surface())!=&(tsos.surface())){
    throw cms::Exception("SiTrackerMultiRecHitUpdator") << "the Trajectory state and the first rechit "
    "passed to the SiTrackerMultiRecHitUpdator lay on different surfaces!: state lays on surface " 
        << tsos.surface().position() << " hit with detid " << (*iter)->det()->geographicalId().rawId() 
        << " lays on surface " << (*iter)->det()->surface().position();
      }

      geomdet = (*iter)->det();

    }

    //if the rechit does not belong to the surface of the tsos
    //GenericProjectedRecHit2D is used to prepagate
    if (&((*iter)->det()->surface())!=&(tsos.surface())){

      TransientTrackingRecHit::RecHitPointer cloned = theHitPropagator->project<GenericProjectedRecHit2D>(*iter,
                                         *geomdet, tsos, theBuilder);
      //if it is used a sensor by sensor grouping this should not appear
      if (cloned->isValid()) updatedcomponents.push_back(cloned);

    } else {
      TransientTrackingRecHit::RecHitPointer cloned = theHitCloner.makeShared(*iter,tsos);
      if (cloned->isValid()){
        updatedcomponents.push_back(cloned);
      }
    }
  } 

  std::vector<std::pair<const TrackingRecHit*, float> > mymap;
  std::vector<std::pair<const TrackingRecHit*, float> > normmap;
  
  double a_sum=0, c_sum=0;


  for(std::vector<TransientTrackingRecHit::RecHitPointer>::iterator ihit = updatedcomponents.begin(); 
    ihit != updatedcomponents.end(); ihit++) {

    double a_i = ComputeWeight(tsos, *(*ihit), false, annealing); //exp(-0.5*Chi2)
    LogTrace("SiTrackerMultiRecHitUpdator")<< "\t\t a_i:" << a_i ;
    //double c_i = ComputeWeight(tsos, *(*ihit), true, annealing);  //exp(-0.5*theChi2Cut/annealing)/(2.*M_PI*sqrt(det));
    //LogTrace("SiTrackerMultiRecHitUpdator")<< "\t\t c_i:" << c_i ;
    mymap.push_back(std::pair<const TrackingRecHit*, float>((*ihit)->hit(), a_i));

    a_sum += a_i;
    //with the new definition, the cut weight is computed only once
    if( ihit == updatedcomponents.begin() ) c_sum = ComputeWeight(tsos, *(*ihit), true, annealing);   //exp(-0.5*theChi2Cut/annealing)
  }
  double total_sum = a_sum + c_sum;    
  LogTrace("SiTrackerMultiRecHitUpdator")<< "\t\t c_sum:" << c_sum ;
  LogTrace("SiTrackerMultiRecHitUpdator")<< "\t\t total sum:" << total_sum ;

  unsigned int counter = 0;
  bool invalid = true;
  for(std::vector<TransientTrackingRecHit::RecHitPointer>::iterator ihit = updatedcomponents.begin(); 
    ihit != updatedcomponents.end(); ihit++) {

    double p = ((mymap[counter].second)/total_sum > 1.e-12 ? (mymap[counter].second)/total_sum : 1.e-12);
    //ORCA: float p = ((mymap[counter].second)/total_sum > 0.01 ? (mymap[counter].second)/total_sum : 1.e-6);


    LogTrace("SiTrackerMultiRecHitUpdator")<< "  Component hit type " << typeid(*mymap[counter].first).name()
                                           << " and dim:" << mymap[counter].first->dimension()
                                           << " position (PRECISE!!!)" << mymap[counter].first->localPosition()
                                           << " error " << mymap[counter].first->localPositionError()
                                           << " with weight " << p ;

    if( p > 10e-6 ){
      invalid = false;
      normmap.push_back(std::pair<const TrackingRecHit*,float>(mymap[counter].first, p));
    }

    counter++;
  }
 
  if(!invalid){

    SiTrackerMultiRecHitUpdator::LocalParameters param = calcParameters(tsos, normmap);
    SiTrackerMultiRecHit updated(param.first, param.second, *normmap.front().first->det(), normmap, annealing);
    LogTrace("SiTrackerMultiRecHitUpdator") << " Updated Hit position " << updated.localPosition() 
                            << " updated error " << updated.localPositionError() << std::endl;

    return std::make_shared<SiTrackerMultiRecHit>(param.first, param.second, *normmap.front().first->det(), normmap, annealing);

  } else {
    LogTrace("SiTrackerMultiRecHitUpdator") << " No hits with weight (> 10e-6) have been found for this MRH." << std::endl;
    return std::make_shared<InvalidTrackingRecHit>(measDet.mdet().geomDet(), TrackingRecHit::missing);
  }
}

Member Data Documentation

bool SiTrackerMultiRecHitUpdator::debug_

private

Definition at line 84 of file SiTrackerMultiRecHitUpdator.h.

const std::vector<double> SiTrackerMultiRecHitUpdator::theAnnealingProgram

private

Definition at line 82 of file SiTrackerMultiRecHitUpdator.h.

Referenced by annealingProgram(), and getAnnealingProgram().

const TransientTrackingRecHitBuilder* SiTrackerMultiRecHitUpdator::theBuilder

private

Definition at line 78 of file SiTrackerMultiRecHitUpdator.h.

Referenced by buildMultiRecHit(), getBuilder(), and update().

double SiTrackerMultiRecHitUpdator::theChi2Cut1D

private

Definition at line 80 of file SiTrackerMultiRecHitUpdator.h.

double SiTrackerMultiRecHitUpdator::theChi2Cut2D

private

Definition at line 81 of file SiTrackerMultiRecHitUpdator.h.

TkClonerImpl SiTrackerMultiRecHitUpdator::theHitCloner

private

Definition at line 83 of file SiTrackerMultiRecHitUpdator.h.

Referenced by SiTrackerMultiRecHitUpdator(), and update().

const TrackingRecHitPropagator* SiTrackerMultiRecHitUpdator::theHitPropagator

private

Definition at line 79 of file SiTrackerMultiRecHitUpdator.h.

Referenced by update().