PixelRod.cc

Go to the documentation of this file.

#include "PixelRod.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "TrackingTools/DetLayers/interface/MeasurementEstimator.h"
#include "TrackingTools/DetLayers/interface/DetLayerException.h"


using namespace std;

typedef GeometricSearchDet::DetWithState DetWithState;

PixelRod::PixelRod(vector<const GeomDet*>& theInputDets):
  DetRodOneR(theInputDets.begin(),theInputDets.end())
{
  theBinFinder = BinFinderType(theDets.begin(),theDets.end());

  //--------- DEBUG INFO --------------
  LogDebug("TkDetLayers") << "==== DEBUG PixelRod ====="; 
  for (vector<const GeomDet*>::const_iterator i=theDets.begin();
       i != theDets.end(); i++){
    LogDebug("TkDetLayers") << "PixelRod's Det pos z,perp,eta,phi: " 
                << (**i).position().z() << " , " 
                << (**i).position().perp() << " , " 
                << (**i).position().eta() << " , " 
                << (**i).position().phi() ;
  }
  LogDebug("TkDetLayers") << "==== end DEBUG PixelRod ====="; 
  //--------------------------------------


}

PixelRod::~PixelRod(){

} 

const vector<const GeometricSearchDet*>& 
PixelRod::components() const{
  throw DetLayerException("PixelRod doesn't have GeometricSearchDet components");
}
 
 
pair<bool, TrajectoryStateOnSurface>
PixelRod::compatible( const TrajectoryStateOnSurface& ts, const Propagator&, 
              const MeasurementEstimator&) const{
  edm::LogError("TkDetLayers") << "temporary dummy implementation of PixelRod::compatible()!!" ;
  return pair<bool,TrajectoryStateOnSurface>();
}

void
PixelRod::compatibleDetsV( const TrajectoryStateOnSurface& startingState,
              const Propagator& prop, 
              const MeasurementEstimator& est,
              std::vector<DetWithState> & result ) const
{  
  typedef MeasurementEstimator::Local2DVector Local2DVector;
  TrajectoryStateOnSurface ts = prop.propagate( startingState, specificSurface());
  if (!ts.isValid()) return;  

  GlobalPoint startPos = ts.globalPosition();

  int closest = theBinFinder.binIndex(startPos.z());
  pair<bool,TrajectoryStateOnSurface> closestCompat = 
    theCompatibilityChecker.isCompatible(theDets[closest],startingState, prop, est);

  if ( closestCompat.first) {
    result.push_back( DetWithState( theDets[closest], closestCompat.second));
  }else{
    if(!closestCompat.second.isValid()) return;  // to investigate why this happens
  }

  const Plane& closestPlane( theDets[closest]->specificSurface() );

  Local2DVector maxDistance = 
    est.maximalLocalDisplacement( closestCompat.second, closestPlane);
  
  float detHalfLen = theDets[closest]->surface().bounds().length()/2.;
  
  // explore neighbours
  for (size_t idet=closest+1; idet < theDets.size(); idet++) {
    LocalPoint nextPos( theDets[idet]->surface().toLocal( closestCompat.second.globalPosition()));
    if (fabs(nextPos.y()) < detHalfLen + maxDistance.y()) {
      if ( !add(idet, result, startingState, prop, est)) break;
    } else {
      break;
    }
  }

  for (int idet=closest-1; idet >= 0; idet--) {
    LocalPoint nextPos( theDets[idet]->surface().toLocal( closestCompat.second.globalPosition()));
    if (fabs(nextPos.y()) < detHalfLen + maxDistance.y()) {
      if ( !add(idet, result, startingState, prop, est)) break;
    } else {
      break;
    }
  }
}


void
PixelRod::groupedCompatibleDetsV( const TrajectoryStateOnSurface&,
                 const Propagator&,
                 const MeasurementEstimator&,
                 std::vector<DetGroup> &) const
{
  LogDebug("TkDetLayers") << "dummy implementation of PixelRod::groupedCompatibleDets()" ;
}