PixelSeeding/src/ThirdHitRZPredictionBase.cc

0001 #include "RecoTracker/PixelSeeding/interface/ThirdHitRZPredictionBase.h"
0002 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
0003 #include "TrackingTools/DetLayers/interface/BarrelDetLayer.h"
0004 #include "TrackingTools/DetLayers/interface/ForwardDetLayer.h"
0005 #include "TrackingTools/DetLayers/interface/DetLayer.h"
0006 #include "DataFormats/GeometrySurface/interface/SimpleDiskBounds.h"
0007 #include "RecoTracker/TkMSParametrization/interface/MultipleScatteringParametrisation.h"
0008
0009 ThirdHitRZPredictionBase::ThirdHitRZPredictionBase() : theBarrel(false), theForward(false), theTolerance(0., 0.) {}
0010
0011 ThirdHitRZPredictionBase::ThirdHitRZPredictionBase(float tolerance, const DetLayer* layer)
0012     : theBarrel(false), theForward(false), theTolerance(tolerance, tolerance) {
0013   if (layer)
0014     initLayer(layer);
0015 }
0016
0017 void ThirdHitRZPredictionBase::initLayer(const DetLayer* layer) {
0018   if (layer->location() == GeomDetEnumerators::barrel) {
0019     theBarrel = true;
0020     theForward = false;
0021     const BarrelDetLayer& bl = reinterpret_cast<const BarrelDetLayer&>(*layer);
0022     float halfThickness = bl.surface().bounds().thickness() / 2;
0023     float radius = bl.specificSurface().radius();
0024     theDetRange = Range(radius - halfThickness, radius + halfThickness);
0025     float maxZ = bl.surface().bounds().length() / 2;
0026     theDetSize = Range(-maxZ, maxZ);
0027   } else if (layer->location() == GeomDetEnumerators::endcap) {
0028     theBarrel = false;
0029     theForward = true;
0030     const ForwardDetLayer& fl = reinterpret_cast<const ForwardDetLayer&>(*layer);
0031     float halfThickness = fl.surface().bounds().thickness() / 2;
0032     float zLayer = fl.position().z();
0033     theDetRange = Range(zLayer - halfThickness, zLayer + halfThickness);
0034     const SimpleDiskBounds& diskRadialBounds = static_cast<const SimpleDiskBounds&>(fl.surface().bounds());
0035     theDetSize = Range(diskRadialBounds.innerRadius(), diskRadialBounds.outerRadius());
0036   } else {
0037     theBarrel = theForward = false;
0038   }
0039 }