SiTrackerMultiRecHit.h

Go to the documentation of this file.

#ifndef SiTrackerMultiRecHit_H
#define SiTrackerMultiRecHit_H

#include "DataFormats/TrackerRecHit2D/interface/BaseTrackerRecHit.h"
#include "DataFormats/TrackingRecHit/interface/TrackingRecHitFwd.h"
#include <vector>
#include <map>
/*
A rechit type suitable for tracking algorithm that use soft hit-to-track assignement, 
such as the Deterministic Annealing Filter (DAF) or the Multi Track Filter (MTF).
it contains an OwnVector with the component rechits and a vector of weights 
*/
class SiTrackerMultiRecHit : public BaseTrackerRecHit
{
public:
  typedef BaseTrackerRecHit Base;
  SiTrackerMultiRecHit():
    theHits(),
    theWeights(),
    annealing_(0){}
  virtual ~SiTrackerMultiRecHit(){} 
  
  
  SiTrackerMultiRecHit(const LocalPoint&, const LocalError&, GeomDet const & idet,
               const std::vector< std::pair<const TrackingRecHit*, float> >&, double);
  
  virtual SiTrackerMultiRecHit* clone() const {return new SiTrackerMultiRecHit(*this);}
#ifdef NO_DICT
  virtual RecHitPointer cloneSH() const { return std::make_shared<SiTrackerMultiRecHit>(*this);}
#endif
  
//  virtual int dimension() const {return 2;}
  virtual int dimension() const; 
  virtual void getKfComponents( KfComponentsHolder & holder ) const;

  // at the momement nobody care of MultiHit!!!
  // used by trackMerger (to be improved)
  virtual OmniClusterRef const & firstClusterRef() const { return static_cast<BaseTrackerRecHit const *>(&theHits.front())->firstClusterRef();}

  virtual std::vector<const TrackingRecHit*> recHits() const;
   
  virtual std::vector<TrackingRecHit*> recHits() ;
  
  //vector of weights
  std::vector<float> const & weights() const {return theWeights;}
  std::vector<float>  & weights() {return theWeights;}

  //returns the weight for the i component
  using TrackingRecHit::weight;
  float  weight(unsigned int i) const {return theWeights[i];}
  float  & weight(unsigned int i) {return theWeights[i];}

  //get the annealing
  virtual double getAnnealingFactor() const { return annealing_; }
    
  bool sharesInput(const TrackingRecHit* other,
           SharedInputType what) const;

private:
  
  edm::OwnVector<TrackingRecHit> theHits;
  std::vector<float> theWeights;
  double annealing_;    

};

// Comparison operators
inline bool operator<( const SiTrackerMultiRecHit& one, const SiTrackerMultiRecHit& other) {
  if ( one.geographicalId() < other.geographicalId() ) {
    return true;
  } else {
    return false;
  }
}

#endif