#include <PixelTripletNoTipGenerator.h>

Inheritance diagram for PixelTripletNoTipGenerator:

Public Member Functions
void	hitTriplets (const TrackingRegion &region, OrderedHitTriplets &trs, const edm::Event &ev, const edm::EventSetup &es, const SeedingLayerSetsHits::SeedingLayerSet &pairLayers, const std::vector< SeedingLayerSetsHits::SeedingLayer > &thirdLayers) override

void	hitTriplets (const TrackingRegion &region, OrderedHitTriplets &result, const edm::EventSetup &es, const HitDoublets &doublets, const RecHitsSortedInPhi *thirdHitMap, const std::vector< const DetLayer > &thirdLayerDetLayer, const int nThirdLayers) override

	PixelTripletNoTipGenerator (const edm::ParameterSet &cfg, edm::ConsumesCollector &iC)

	~PixelTripletNoTipGenerator () override

Public Member Functions inherited from HitTripletGeneratorFromPairAndLayers
	HitTripletGeneratorFromPairAndLayers (unsigned int maxElement=0)

	HitTripletGeneratorFromPairAndLayers (const edm::ParameterSet &pset)

void	init (std::unique_ptr< HitPairGeneratorFromLayerPair > &&pairs, LayerCacheType *layerCache)

const HitPairGeneratorFromLayerPair &	pairGenerator () const

virtual	~HitTripletGeneratorFromPairAndLayers ()

Private Types
typedef CombinedHitTripletGenerator::LayerCacheType	LayerCacheType

Private Attributes
float	extraHitPhiToleranceForPreFiltering

float	extraHitRPhitolerance

float	extraHitRZtolerance

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	fieldToken_

edm::ESGetToken < MultipleScatteringParametrisationMaker, TrackerMultipleScatteringRecord >	msmakerToken_

double	theChi2Cut

double	theNSigma

Additional Inherited Members
Public Types inherited from HitTripletGeneratorFromPairAndLayers
typedef LayerHitMapCache	LayerCacheType

Static Public Member Functions inherited from HitTripletGeneratorFromPairAndLayers
static void	fillDescriptions (edm::ParameterSetDescription &desc)

Protected Attributes inherited from HitTripletGeneratorFromPairAndLayers
LayerCacheType *	theLayerCache

const unsigned int	theMaxElement

std::unique_ptr < HitPairGeneratorFromLayerPair >	thePairGenerator

Detailed Description

Definition at line 15 of file PixelTripletNoTipGenerator.h.

Member Typedef Documentation

typedef CombinedHitTripletGenerator::LayerCacheType PixelTripletNoTipGenerator::LayerCacheType

private

Definition at line 16 of file PixelTripletNoTipGenerator.h.

Constructor & Destructor Documentation

PixelTripletNoTipGenerator::PixelTripletNoTipGenerator	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &	iC
	)

Definition at line 28 of file PixelTripletNoTipGenerator.cc.

     : HitTripletGeneratorFromPairAndLayers(cfg),
       fieldToken_(iC.esConsumes()),
       msmakerToken_(iC.esConsumes()),
       extraHitRZtolerance(cfg.getParameter<double>("extraHitRZtolerance")),
       extraHitRPhitolerance(cfg.getParameter<double>("extraHitRPhitolerance")),
       extraHitPhiToleranceForPreFiltering(cfg.getParameter<double>("extraHitPhiToleranceForPreFiltering")),
       theNSigma(cfg.getParameter<double>("nSigma")),
       theChi2Cut(cfg.getParameter<double>("chi2Cut")) {}

PixelTripletNoTipGenerator::~PixelTripletNoTipGenerator ( )

override

Definition at line 38 of file PixelTripletNoTipGenerator.cc.

38 {}

Member Function Documentation

void PixelTripletNoTipGenerator::hitTriplets	(	const TrackingRegion &	region,
		OrderedHitTriplets &	trs,
		const edm::Event &	ev,
		const edm::EventSetup &	es,
		const SeedingLayerSetsHits::SeedingLayerSet &	pairLayers,
		const std::vector< SeedingLayerSetsHits::SeedingLayer > &	thirdLayers
	)

overridevirtual

Implements HitTripletGeneratorFromPairAndLayers.

Definition at line 40 of file PixelTripletNoTipGenerator.cc.

References GeomDetEnumerators::barrel, Surface::bounds(), PixelRecoUtilities::curvature(), runTauDisplay::dr, extraHitPhiToleranceForPreFiltering, extraHitRZtolerance, fieldToken_, edm::EventSetup::getData(), RecHitsSortedInPhi::hits(), PixelRecoUtilities::inversePt(), phase1PixelTopology::layer, geometryCSVtoXML::line, DetLayer::location(), M_PI, msmakerToken_, TrackingRegion::origin(), TrackingRegion::originRBound(), fireworks::p1, fireworks::p2, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GeometricSearchDet::position(), TrackingRegion::ptMin(), CosmicsPD_Skims::radius, HLT_FULL_cff::region, edm::shift, OrderedHitTriplets::size(), findQualityFiles::size, PixelRecoRange< T >::sort(), BarrelDetLayer::specificSurface(), sqr(), mathSSE::sqrt(), BarrelDetLayer::surface(), ForwardDetLayer::surface(), edm::swap(), theChi2Cut, HitTripletGeneratorFromPairAndLayers::theMaxElement, theNSigma, HitTripletGeneratorFromPairAndLayers::thePairGenerator, Bounds::thickness(), PV3DBase< T, PVType, FrameType >::z(), and beampixel_dqm_sourceclient-live_cfg::zRange.

                                                                                                              {
   //
   //edm::Handle<reco::BeamSpot> bsHandle;
   //ev.getByLabel( theBeamSpotTag, bsHandle);
   //if(!bsHandle.isValid()) return;
   //const reco::BeamSpot & bs = *bsHandle;
   //double errorXY = sqrt( sqr(bs.BeamWidthX()) + sqr(bs.BeamWidthY()) );
   //
 
   const GlobalPoint& bsPoint = region.origin();
   double errorXY = region.originRBound();
   GlobalVector shift = bsPoint - GlobalPoint(0., 0., 0.);
 
   OrderedHitPairs pairs;
   pairs.reserve(30000);
   OrderedHitPairs::const_iterator ip;
   thePairGenerator->hitPairs(region, pairs, ev, es, pairLayers);
 
   if (pairs.empty())
     return;
 
   const DetLayer* firstLayer = thePairGenerator->innerLayer(pairLayers).detLayer();
   const DetLayer* secondLayer = thePairGenerator->outerLayer(pairLayers).detLayer();
   if (!firstLayer || !secondLayer)
     return;
 
   int size = thirdLayers.size();
 
   const RecHitsSortedInPhi** thirdHitMap = new const RecHitsSortedInPhi*[size];
   for (int il = 0; il <= size - 1; il++) {
     thirdHitMap[il] = &(*theLayerCache)(thirdLayers[il], region);
   }
 
   const auto& field = es.getData(fieldToken_);
   const auto& msmaker = es.getData(msmakerToken_);
 
   MultipleScatteringParametrisation sigma1RPhi = msmaker.parametrisation(firstLayer);
   MultipleScatteringParametrisation sigma2RPhi = msmaker.parametrisation(secondLayer);
 
   typedef RecHitsSortedInPhi::Hit Hit;
   for (ip = pairs.begin(); ip != pairs.end(); ip++) {
     GlobalPoint p1((*ip).inner()->globalPosition() - shift);
     GlobalPoint p2((*ip).outer()->globalPosition() - shift);
 
     ThirdHitPredictionFromInvLine predictionRPhiTMP(p1, p2);
     double pt_p1p2 = 1. / PixelRecoUtilities::inversePt(predictionRPhiTMP.curvature(), field);
 
     PixelRecoPointRZ point1(p1.perp(), p1.z());
     PixelRecoPointRZ point2(p2.perp(), p2.z());
 
     PixelRecoLineRZ line(point1, point2);
     double msRPhi1 = sigma1RPhi(pt_p1p2, line.cotLine(), 0.f);
     double msRPhi2 = sigma2RPhi(pt_p1p2, line.cotLine(), point1);
     double sinTheta = 1 / sqrt(1 + sqr(line.cotLine()));
     double cosTheta = fabs(line.cotLine()) / sqrt(1 + sqr(line.cotLine()));
 
     double p1_errorRPhi = sqrt(sqr((*ip).inner()->errorGlobalRPhi()) + sqr(msRPhi1) + sqr(errorXY));
     double p2_errorRPhi = sqrt(sqr((*ip).outer()->errorGlobalRPhi()) + sqr(msRPhi2) + sqr(errorXY));
 
     ThirdHitPredictionFromInvLine predictionRPhi(p1, p2, p1_errorRPhi, p2_errorRPhi);
 
     for (int il = 0; il <= size - 1; il++) {
       const DetLayer* layer = thirdLayers[il].detLayer();
       bool barrelLayer = (layer->location() == GeomDetEnumerators::barrel);
       MultipleScatteringParametrisation sigma3RPhi = msmaker.parametrisation(layer);
       double msRPhi3 = sigma3RPhi(pt_p1p2, line.cotLine(), point2);
 
       Range rRange;
       if (barrelLayer) {
         const BarrelDetLayer& bl = dynamic_cast<const BarrelDetLayer&>(*layer);
         float halfThickness = bl.surface().bounds().thickness() / 2;
         float radius = bl.specificSurface().radius();
         rRange = Range(radius - halfThickness, radius + halfThickness);
       } else {
         const ForwardDetLayer& fl = dynamic_cast<const ForwardDetLayer&>(*layer);
         float halfThickness = fl.surface().bounds().thickness() / 2;
         float zLayer = fl.position().z();
         float zMin = zLayer - halfThickness;
         float zMax = zLayer + halfThickness;
         GlobalVector dr = p2 - p1;
         GlobalPoint p3_a = p2 + dr * (zMin - p2.z()) / dr.z();
         GlobalPoint p3_b = p2 + dr * (zMax - p2.z()) / dr.z();
         if (zLayer * p3_a.z() < 0)
           continue;
         rRange = Range(p3_a.perp(), p3_b.perp());
         rRange.sort();
       }
       double displacment = shift.perp();
       GlobalPoint crossing1 = predictionRPhi.crossing(rRange.min() - displacment) + shift;
       GlobalPoint crossing2 = predictionRPhi.crossing(rRange.max() + displacment) + shift;
       float c1_phi = crossing1.phi();
       float c2_phi = crossing2.phi();
       if (c2_phi < c1_phi)
         swap(c1_phi, c2_phi);
       if (c2_phi - c1_phi > M_PI) {
         c2_phi -= 2 * M_PI;
         swap(c1_phi, c2_phi);
       }
       double extraAngle = (displacment + theNSigma * msRPhi3) / rRange.min() + extraHitPhiToleranceForPreFiltering;
       c1_phi -= extraAngle;
       c2_phi += extraAngle;
       vector<Hit> thirdHits = thirdHitMap[il]->hits(c1_phi, c2_phi);
 
       typedef vector<Hit>::const_iterator IH;
       for (IH th = thirdHits.begin(), eh = thirdHits.end(); th < eh; ++th) {
         GlobalPoint p3((*th)->globalPosition() - shift);
         double p3_errorRPhi = sqrt(sqr((*th)->errorGlobalRPhi()) + sqr(msRPhi3) + sqr(errorXY));
 
         predictionRPhi.add(p3, p3_errorRPhi);
 
         double curvature = predictionRPhi.curvature();
 
         ThirdHitZPrediction zPrediction((*ip).inner()->globalPosition(),
                                         sqrt(sqr((*ip).inner()->errorGlobalR()) + sqr(msRPhi1 / cosTheta)),
                                         sqrt(sqr((*ip).inner()->errorGlobalZ()) + sqr(msRPhi1 / sinTheta)),
                                         (*ip).outer()->globalPosition(),
                                         sqrt(sqr((*ip).outer()->errorGlobalR()) + sqr(msRPhi2 / cosTheta)),
                                         sqrt(sqr((*ip).outer()->errorGlobalZ()) + sqr(msRPhi2 / sinTheta)),
                                         curvature,
                                         theNSigma);
         ThirdHitZPrediction::Range zRange =
             zPrediction((*th)->globalPosition(), sqrt(sqr((*th)->errorGlobalR())) + sqr(msRPhi3 / cosTheta));
 
         double z3Hit = (*th)->globalPosition().z();
         double z3HitError =
             theNSigma * (sqrt(sqr((*th)->errorGlobalZ()) + sqr(msRPhi3 / sinTheta))) + extraHitRZtolerance;
         Range hitZRange(z3Hit - z3HitError, z3Hit + z3HitError);
         bool inside = hitZRange.hasIntersection(zRange);
 
         double curvatureMS = PixelRecoUtilities::curvature(1. / region.ptMin(), field);
         bool ptCut = (predictionRPhi.curvature() - theNSigma * predictionRPhi.errorCurvature() < curvatureMS);
         bool chi2Cut = (predictionRPhi.chi2() < theChi2Cut);
         if (inside && ptCut && chi2Cut) {
           if (theMaxElement != 0 && result.size() >= theMaxElement) {
             result.clear();
             edm::LogError("TooManyTriplets") << " number of triples exceed maximum. no triplets produced.";
             delete[] thirdHitMap;
             return;
           }
           result.push_back(OrderedHitTriplet((*ip).inner(), (*ip).outer(), *th));
         }
         predictionRPhi.remove(p3, p3_errorRPhi);
       }
     }
   }
   delete[] thirdHitMap;
 }

void PixelTripletNoTipGenerator::hitTriplets	(	const TrackingRegion &	region,
		OrderedHitTriplets &	result,
		const edm::EventSetup &	es,
		const HitDoublets &	doublets,
		const RecHitsSortedInPhi **	thirdHitMap,
		const std::vector< const DetLayer * > &	thirdLayerDetLayer,
		const int	nThirdLayers
	)

overridevirtual

Implements HitTripletGeneratorFromPairAndLayers.

Definition at line 192 of file PixelTripletNoTipGenerator.cc.

References Exception.

                                                                      {
   throw cms::Exception("Error") << "PixelTripletNoTipGenerator::hitTriplets is not implemented \n";
 }