/data/refman/pasoursint/CMSSW_5_3_8_patch3/src/TrackingTools/DetLayers/src/DetRodOneR.cc

Go to the documentation of this file.

#include "TrackingTools/DetLayers/interface/DetRodOneR.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h" 
#include "TrackingTools/DetLayers/interface/RodPlaneBuilderFromDet.h"

#include <Utilities/General/interface/precomputed_value_sort.h>
#include <Geometry/CommonDetUnit/interface/DetSorting.h>

#include <algorithm>
#include <cmath>

using namespace std;

DetRodOneR::~DetRodOneR(){}

DetRodOneR::DetRodOneR(vector<const GeomDet*>::const_iterator first,
                       vector<const GeomDet*>::const_iterator last)
  : theDets(first,last)
{
  initialize();
}

DetRodOneR::DetRodOneR( const vector<const GeomDet*>& dets)
  : theDets(dets) 
{
  initialize();
}


void DetRodOneR::initialize()
{
  // assume the dets ARE in a rod;
  // sort them in Z

  precomputed_value_sort( theDets.begin(), theDets.end(), geomsort::DetZ());
  
  setPlane( RodPlaneBuilderFromDet()( theDets));
  
}


// It needs that the basic component to have the compatible() method
bool DetRodOneR::add( int idet, vector<DetWithState>& result,
                      const TrajectoryStateOnSurface& startingState,
                      const Propagator& prop, 
                      const MeasurementEstimator& est) const
{
  pair<bool,TrajectoryStateOnSurface> compat = 
    theCompatibilityChecker.isCompatible(theDets[idet],startingState, prop, est);
  
  if (compat.first) {
    result.push_back( DetWithState( theDets[idet], compat.second));
  }

  return compat.first;
}