#include <PixelTripletNoTipGenerator.h>

Inheritance diagram for PixelTripletNoTipGenerator:

Public Member Functions
virtual void	hitTriplets (const TrackingRegion &region, OrderedHitTriplets &trs, const edm::Event &ev, const edm::EventSetup &es)

void	init (const HitPairGenerator &pairs, LayerCacheType *layerCache) override

const HitPairGenerator &	pairGenerator () const

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

void	setSeedingLayers (SeedingLayerSetsHits::SeedingLayerSet pairLayers, std::vector< SeedingLayerSetsHits::SeedingLayer > thirdLayers) override

virtual	~PixelTripletNoTipGenerator ()

Public Member Functions inherited from HitTripletGeneratorFromPairAndLayers
virtual	~HitTripletGeneratorFromPairAndLayers ()

Public Member Functions inherited from HitTripletGenerator
virtual void	clear ()

	HitTripletGenerator (unsigned int size=500)

virtual void	hitTriplets (const TrackingRegion &reg, OrderedHitTriplets &prs, const edm::EventSetup &es)

virtual const OrderedHitTriplets &	run (const TrackingRegion &region, const edm::Event &ev, const edm::EventSetup &es)

virtual	~HitTripletGenerator ()

Public Member Functions inherited from OrderedHitsGenerator
	OrderedHitsGenerator ()

virtual	~OrderedHitsGenerator ()

Private Types
typedef CombinedHitTripletGenerator::LayerCacheType	LayerCacheType

Private Attributes
float	extraHitPhiToleranceForPreFiltering

float	extraHitRPhitolerance

float	extraHitRZtolerance

double	theChi2Cut

LayerCacheType *	theLayerCache

std::vector < SeedingLayerSetsHits::SeedingLayer >	theLayers

double	theNSigma

HitPairGenerator *	thePairGenerator

Additional Inherited Members
Public Types inherited from HitTripletGeneratorFromPairAndLayers
typedef LayerHitMapCache	LayerCacheType

Public Attributes inherited from OrderedHitsGenerator
unsigned int	theMaxElement

Detailed Description

Definition at line 14 of file PixelTripletNoTipGenerator.h.

Member Typedef Documentation

typedef CombinedHitTripletGenerator::LayerCacheType PixelTripletNoTipGenerator::LayerCacheType

private

Definition at line 15 of file PixelTripletNoTipGenerator.h.

Constructor & Destructor Documentation

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

Definition at line 24 of file PixelTripletNoTipGenerator.cc.

     : thePairGenerator(0),
       theLayerCache(0),
       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"))
 { }

virtual PixelTripletNoTipGenerator::~PixelTripletNoTipGenerator ( )

inlinevirtual

Definition at line 19 of file PixelTripletNoTipGenerator.h.

References thePairGenerator.

19 { delete thePairGenerator; }

PixelTripletNoTipGenerator::thePairGenerator

HitPairGenerator * thePairGenerator

Definition: PixelTripletNoTipGenerator.h:32

Member Function Documentation

void PixelTripletNoTipGenerator::hitTriplets	(	const TrackingRegion &	region,
		OrderedHitTriplets &	trs,
		const edm::Event &	ev,
		const edm::EventSetup &	es
	)

virtual

Implements HitTripletGenerator.

Definition at line 47 of file PixelTripletNoTipGenerator.cc.

References GeomDetEnumerators::barrel, Surface::bounds(), PixelRecoUtilities::curvature(), SeedingLayerSetsHits::SeedingLayer::detLayer(), extraHitPhiToleranceForPreFiltering, extraHitRZtolerance, HitPairGenerator::hitPairs(), RecHitsSortedInPhi::hits(), HitPairGeneratorFromLayerPair::innerLayer(), PixelRecoUtilities::inversePt(), geometryCSVtoXML::line, DetLayer::location(), M_PI, TrackingRegion::origin(), TrackingRegion::originRBound(), HitPairGeneratorFromLayerPair::outerLayer(), p1, p2, p3, filesave_online::pairs, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GeometricSearchDet::position(), TrackingRegion::ptMin(), CosmicsPD_Skims::radius, edm::shift, OrderedHitPairs::size(), OrderedHitTriplets::size(), findQualityFiles::size, PixelRecoRange< T >::sort(), BarrelDetLayer::specificSurface(), funct::sqr(), mathSSE::sqrt(), ForwardDetLayer::surface(), BarrelDetLayer::surface(), swap(), theChi2Cut, theLayers, OrderedHitsGenerator::theMaxElement, theNSigma, 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()) );
 //
 
   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);
 
   if (pairs.size() ==0) return;
 
   int size = theLayers.size();
 
   const RecHitsSortedInPhi **thirdHitMap = new const RecHitsSortedInPhi*[size];
   for (int il=0; il <=size-1; il++) {
      thirdHitMap[il] = &(*theLayerCache)(theLayers[il], region, ev, es);
   }
 
   const HitPairGeneratorFromLayerPair * pairGen = dynamic_cast<const HitPairGeneratorFromLayerPair *>(thePairGenerator);
   const DetLayer * firstLayer = pairGen->innerLayer().detLayer();
   const DetLayer * secondLayer = pairGen->outerLayer().detLayer();
   if (!firstLayer || !secondLayer) return;
 
   MultipleScatteringParametrisation sigma1RPhi( firstLayer, es);
   MultipleScatteringParametrisation sigma2RPhi( secondLayer, es);
 
   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(),es);
 
     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 = theLayers[il].detLayer();
       bool barrelLayer = (layer->location() == GeomDetEnumerators::barrel);
       MultipleScatteringParametrisation sigma3RPhi( layer, es);
       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(),es);
          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::init	(	const HitPairGenerator &	pairs,
		LayerCacheType *	layerCache
	)

overridevirtual

Implements HitTripletGeneratorFromPairAndLayers.

Definition at line 34 of file PixelTripletNoTipGenerator.cc.

References HitPairGenerator::clone(), theLayerCache, and thePairGenerator.

 {
   thePairGenerator = pairs.clone();
   theLayerCache = layerCache;
 }

const HitPairGenerator& PixelTripletNoTipGenerator::pairGenerator ( ) const

inline

Definition at line 29 of file PixelTripletNoTipGenerator.h.

References thePairGenerator.

29 { return *thePairGenerator; }

PixelTripletNoTipGenerator::thePairGenerator

HitPairGenerator * thePairGenerator

Definition: PixelTripletNoTipGenerator.h:32

void PixelTripletNoTipGenerator::setSeedingLayers	(	SeedingLayerSetsHits::SeedingLayerSet	pairLayers,
		std::vector< SeedingLayerSetsHits::SeedingLayer >	thirdLayers
	)

overridevirtual

Implements HitTripletGeneratorFromPairAndLayers.

Definition at line 41 of file PixelTripletNoTipGenerator.cc.

References HitPairGenerator::setSeedingLayers(), theLayers, and thePairGenerator.

                                                                                                            {
   thePairGenerator->setSeedingLayers(pairLayers);
   theLayers = thirdLayers;
 }