SiTrackerMultiRecHit.cc

Go to the documentation of this file.

#include "DataFormats/TrackerRecHit2D/interface/SiTrackerMultiRecHit.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

using namespace std;
using namespace edm;

SiTrackerMultiRecHit::SiTrackerMultiRecHit(const LocalPoint& pos, const LocalError& err, GeomDet const & idet,
                       const std::vector< std::pair<const TrackingRecHit*, float> >& aHitMap, double annealing):
  BaseTrackerRecHit(pos,err, idet,trackerHitRTTI::multi)    
{
  for(std::vector<std::pair<const TrackingRecHit*, float> >::const_iterator ihit = aHitMap.begin(); ihit != aHitMap.end(); ihit++){
    theHits.push_back(ihit->first->clone());
    theWeights.push_back(ihit->second);
  }
  annealing_ = annealing;
}


bool SiTrackerMultiRecHit::sharesInput(const TrackingRecHit* other,
                       SharedInputType what) const
{
  if(geographicalId() != other->geographicalId()&& what==all ) return false;
  vector<const TrackingRecHit*> otherhits=other->recHits();
  if(what==all){
    if(theHits.size()!=other->recHits().size())return false;
    for(vector<const TrackingRecHit*>::iterator otherhit=otherhits.begin();otherhit!=otherhits.end();++otherhit){
      bool found=false;
      for(OwnVector<TrackingRecHit>::const_iterator hit=theHits.begin();hit!=theHits.end();++hit){
    if((hit)->sharesInput(*otherhit,all)){
      found=true;
      break;
    }
      }
      if(found==false){
    return false;
      }
    }
    return true;
  }
  else{
    for(OwnVector<TrackingRecHit>::const_iterator hit=theHits.begin();hit!=theHits.end();++hit){
      if(otherhits.size()!=0){ 
    for(vector<const TrackingRecHit*>::iterator otherhit=otherhits.begin();otherhit!=otherhits.end();++otherhit){
      if((hit)->sharesInput(*otherhit,some))return true;
    }
      }
      else{//otherwise it should be a single rechit
    if((hit)->sharesInput(other,some))return true;
      } 
    }
    return false;
  }
}


vector<const TrackingRecHit*> SiTrackerMultiRecHit::recHits() const{
  vector<const TrackingRecHit*> myhits;
  for(edm::OwnVector<TrackingRecHit>::const_iterator ihit = theHits.begin(); ihit != theHits.end(); ihit++) {
    myhits.push_back(&*ihit);
  }
  return myhits;
}

vector<TrackingRecHit*> SiTrackerMultiRecHit::recHits() {
  return theHits.data();
}


int SiTrackerMultiRecHit::dimension() const{
  //supposing all the hits inside of a MRH have the same id == same type
  int randomComponent = 0; 
  if(theHits[randomComponent].dimension() == 1 ){  return 1;  }
  else if(theHits[randomComponent].dimension() == 2 ){  return 2;  }
  else {  return 0;  }
}

void SiTrackerMultiRecHit::getKfComponents( KfComponentsHolder & holder ) const { 
  if (dimension() == 1)  getKfComponents1D(holder); 
  if (dimension() == 2)  getKfComponents2D(holder); 
}