CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Types | Private Attributes
PixelTripletNoTipGenerator Class Reference

#include <PixelTripletNoTipGenerator.h>

Inheritance diagram for PixelTripletNoTipGenerator:
HitTripletGeneratorFromPairAndLayers

Public Member Functions

virtual 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)
 
virtual ~PixelTripletNoTipGenerator ()
 
- 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
 
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 13 of file PixelTripletNoTipGenerator.h.

Member Typedef Documentation

typedef CombinedHitTripletGenerator::LayerCacheType PixelTripletNoTipGenerator::LayerCacheType
private

Definition at line 14 of file PixelTripletNoTipGenerator.h.

Constructor & Destructor Documentation

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

Definition at line 24 of file PixelTripletNoTipGenerator.cc.

25  : HitTripletGeneratorFromPairAndLayers(cfg),
26  extraHitRZtolerance(cfg.getParameter<double>("extraHitRZtolerance")),
27  extraHitRPhitolerance(cfg.getParameter<double>("extraHitRPhitolerance")),
28  extraHitPhiToleranceForPreFiltering(cfg.getParameter<double>("extraHitPhiToleranceForPreFiltering")),
29  theNSigma(cfg.getParameter<double>("nSigma")),
30  theChi2Cut(cfg.getParameter<double>("chi2Cut"))
31 { }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
PixelTripletNoTipGenerator::~PixelTripletNoTipGenerator ( )
virtual

Definition at line 33 of file PixelTripletNoTipGenerator.cc.

33 {}

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 35 of file PixelTripletNoTipGenerator.cc.

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

42 {
43 
44 //
45 //edm::Handle<reco::BeamSpot> bsHandle;
46 //ev.getByLabel( theBeamSpotTag, bsHandle);
47 //if(!bsHandle.isValid()) return;
48 //const reco::BeamSpot & bs = *bsHandle;
49 //double errorXY = sqrt( sqr(bs.BeamWidthX()) + sqr(bs.BeamWidthY()) );
50 //
51 
52  GlobalPoint bsPoint = region.origin();
53  double errorXY = region.originRBound();
54  GlobalVector shift = bsPoint - GlobalPoint(0.,0.,0.);
55 
56  OrderedHitPairs pairs; pairs.reserve(30000);
57  OrderedHitPairs::const_iterator ip;
58  thePairGenerator->hitPairs(region,pairs,ev,es, pairLayers);
59 
60  if (pairs.size() ==0) return;
61 
62  int size = thirdLayers.size();
63 
64  const RecHitsSortedInPhi **thirdHitMap = new const RecHitsSortedInPhi*[size];
65  for (int il=0; il <=size-1; il++) {
66  thirdHitMap[il] = &(*theLayerCache)(thirdLayers[il], region, es);
67  }
68 
69  const DetLayer * firstLayer = thePairGenerator->innerLayer(pairLayers).detLayer();
70  const DetLayer * secondLayer = thePairGenerator->outerLayer(pairLayers).detLayer();
71  if (!firstLayer || !secondLayer) return;
72 
73  MultipleScatteringParametrisation sigma1RPhi( firstLayer, es);
74  MultipleScatteringParametrisation sigma2RPhi( secondLayer, es);
75 
76  typedef RecHitsSortedInPhi::Hit Hit;
77  for (ip = pairs.begin(); ip != pairs.end(); ip++) {
78 
79  GlobalPoint p1((*ip).inner()->globalPosition()-shift);
80  GlobalPoint p2((*ip).outer()->globalPosition()-shift);
81 
82  ThirdHitPredictionFromInvLine predictionRPhiTMP(p1, p2 );
83  double pt_p1p2 = 1./PixelRecoUtilities::inversePt(predictionRPhiTMP.curvature(),es);
84 
85  PixelRecoPointRZ point1(p1.perp(), p1.z());
86  PixelRecoPointRZ point2(p2.perp(), p2.z());
87 
88  PixelRecoLineRZ line(point1, point2);
89  double msRPhi1 = sigma1RPhi(pt_p1p2, line.cotLine(), 0.f);
90  double msRPhi2 = sigma2RPhi(pt_p1p2, line.cotLine(),point1);
91  double sinTheta = 1/sqrt(1+sqr(line.cotLine()));
92  double cosTheta = fabs(line.cotLine())/sqrt(1+sqr(line.cotLine()));
93 
94  double p1_errorRPhi = sqrt(sqr((*ip).inner()->errorGlobalRPhi())+sqr(msRPhi1) +sqr(errorXY));
95  double p2_errorRPhi = sqrt(sqr((*ip).outer()->errorGlobalRPhi())+sqr(msRPhi2) +sqr(errorXY));
96 
97  ThirdHitPredictionFromInvLine predictionRPhi(p1, p2, p1_errorRPhi, p2_errorRPhi );
98 
99  for (int il=0; il <=size-1; il++) {
100 
101  const DetLayer * layer = thirdLayers[il].detLayer();
102  bool barrelLayer = (layer->location() == GeomDetEnumerators::barrel);
103  MultipleScatteringParametrisation sigma3RPhi( layer, es);
104  double msRPhi3 = sigma3RPhi(pt_p1p2, line.cotLine(),point2);
105 
106  Range rRange;
107  if (barrelLayer) {
108  const BarrelDetLayer& bl = dynamic_cast<const BarrelDetLayer&>(*layer);
109  float halfThickness = bl.surface().bounds().thickness()/2;
110  float radius = bl.specificSurface().radius();
111  rRange = Range(radius-halfThickness, radius+halfThickness);
112  } else {
113  const ForwardDetLayer& fl = dynamic_cast<const ForwardDetLayer&>(*layer);
114  float halfThickness = fl.surface().bounds().thickness()/2;
115  float zLayer = fl.position().z() ;
116  float zMin = zLayer-halfThickness;
117  float zMax = zLayer+halfThickness;
118  GlobalVector dr = p2-p1;
119  GlobalPoint p3_a = p2+dr*(zMin-p2.z())/dr.z();
120  GlobalPoint p3_b = p2+dr*(zMax-p2.z())/dr.z();
121  if (zLayer * p3_a.z() < 0) continue;
122  rRange = Range(p3_a.perp(), p3_b.perp());
123  rRange.sort();
124  }
125  double displacment = shift.perp();
126  GlobalPoint crossing1 = predictionRPhi.crossing(rRange.min()-displacment)+shift;
127  GlobalPoint crossing2 = predictionRPhi.crossing(rRange.max()+displacment)+shift;
128  float c1_phi= crossing1.phi();
129  float c2_phi= crossing2.phi();
130  if (c2_phi < c1_phi) swap(c1_phi,c2_phi);
131  if (c2_phi-c1_phi > M_PI) { c2_phi -= 2*M_PI; swap(c1_phi,c2_phi); }
132  double extraAngle = (displacment+theNSigma*msRPhi3)/rRange.min()+extraHitPhiToleranceForPreFiltering;
133  c1_phi -= extraAngle;
134  c2_phi += extraAngle;
135  vector<Hit> thirdHits = thirdHitMap[il]->hits(c1_phi, c2_phi) ;
136 
137  typedef vector<Hit>::const_iterator IH;
138  for (IH th=thirdHits.begin(), eh=thirdHits.end(); th < eh; ++th) {
139  GlobalPoint p3((*th)->globalPosition()-shift);
140  double p3_errorRPhi = sqrt(sqr((*th)->errorGlobalRPhi()) +sqr(msRPhi3) + sqr(errorXY));
141 
142  predictionRPhi.add(p3,p3_errorRPhi);
143 
144  double curvature = predictionRPhi.curvature();
145 
146  ThirdHitZPrediction zPrediction( (*ip).inner()->globalPosition(), sqrt(sqr((*ip).inner()->errorGlobalR())+sqr(msRPhi1/cosTheta)), sqrt( sqr((*ip).inner()->errorGlobalZ())+ sqr(msRPhi1/sinTheta)),
147  (*ip).outer()->globalPosition(), sqrt(sqr((*ip).outer()->errorGlobalR())+sqr(msRPhi2/cosTheta)), sqrt( sqr((*ip).outer()->errorGlobalZ())+sqr(msRPhi2/sinTheta)),
148  curvature, theNSigma);
149  ThirdHitZPrediction::Range zRange = zPrediction((*th)->globalPosition(), sqrt(sqr((*th)->errorGlobalR()))+sqr(msRPhi3/cosTheta));
150 
151  double z3Hit = (*th)->globalPosition().z();
152  double z3HitError = theNSigma*(sqrt(sqr((*th)->errorGlobalZ()) + sqr(msRPhi3/sinTheta) ))+extraHitRZtolerance;
153  Range hitZRange(z3Hit-z3HitError, z3Hit+z3HitError);
154  bool inside = hitZRange.hasIntersection(zRange);
155 
156  double curvatureMS = PixelRecoUtilities::curvature(1./region.ptMin(),es);
157  bool ptCut = (predictionRPhi.curvature()-theNSigma*predictionRPhi.errorCurvature() < curvatureMS);
158  bool chi2Cut = (predictionRPhi.chi2() < theChi2Cut);
159  if (inside && ptCut && chi2Cut) {
160  if (theMaxElement!=0 && result.size() >= theMaxElement){
161  result.clear();
162  edm::LogError("TooManyTriplets")<<" number of triples exceed maximum. no triplets produced.";
163  delete [] thirdHitMap;
164  return;
165  }
166  result.push_back( OrderedHitTriplet( (*ip).inner(), (*ip).outer(), *th));
167  }
168  predictionRPhi.remove(p3,p3_errorRPhi);
169  }
170  }
171  }
172  delete [] thirdHitMap;
173 }
TrackingRegion::originRBound
float originRBound() const
bounds the particle vertex in the transverse plane
Definition: TrackingRegion.h:77
RecHitsSortedInPhi::hits
std::vector< Hit > hits(float phiMin, float phiMax) const
Definition: RecHitsSortedInPhi.cc:98
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:72
TrackingRegion::origin
GlobalPoint const & origin() const
Definition: TrackingRegion.h:74
PixelRecoRange::sort
PixelRecoRange< T > & sort()
Definition: PixelRecoRange.h:60
DetLayer::location
virtual Location location() const =0
Which part of the detector (barrel, endcap)
BarrelDetLayer::surface
virtual const BoundSurface & surface() const final
GeometricSearchDet interface.
Definition: BarrelDetLayer.h:32
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
Surface::bounds
const Bounds & bounds() const
Definition: Surface.h:120
OrderedHitPairs::size
virtual unsigned int size() const
Definition: OrderedHitPairs.h:13
edm::swap
void swap(Association< C > &lhs, Association< C > &rhs)
Definition: Association.h:116
beampixel_dqm_sourceclient-live_cfg.zRange
tuple zRange
Definition: beampixel_dqm_sourceclient-live_cfg.py:102
HLT_FULL_cff.region
tuple region
Definition: HLT_FULL_cff.py:87245
mps_fire.result
tuple result
Definition: mps_fire.py:84
Bounds::thickness
virtual float thickness() const =0
PixelRecoUtilities::inversePt
T inversePt(T curvature, const edm::EventSetup &iSetup)
Definition: PixelRecoUtilities.h:45
PixelRecoUtilities::curvature
T curvature(T InversePt, const edm::EventSetup &iSetup)
Definition: PixelRecoUtilities.h:41
Range
PixelRecoRange< float > Range
Definition: PixelQuadrupletGenerator.cc:31
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
p2
double p2[4]
Definition: TauolaWrapper.h:90
Hit
SeedingHitSet::ConstRecHitPointer Hit
Definition: SeedGeneratorFromProtoTracksEDProducer.cc:32
RecHitsSortedInPhi::Hit
BaseTrackerRecHit const * Hit
Definition: RecHitsSortedInPhi.h:18
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:51
HitTripletGeneratorFromPairAndLayers::thePairGenerator
std::unique_ptr< HitPairGeneratorFromLayerPair > thePairGenerator
Definition: HitTripletGeneratorFromPairAndLayers.h:49
GeometricSearchDet::position
virtual const Surface::PositionType & position() const
Returns position of the surface.
Definition: GeometricSearchDet.h:33
BarrelDetLayer::specificSurface
virtual const BoundCylinder & specificSurface() const final
Extension of the interface.
Definition: BarrelDetLayer.h:43
TrackingRegion::ptMin
float ptMin() const
minimal pt of interest
Definition: TrackingRegion.h:83
ForwardDetLayer::surface
virtual const BoundSurface & surface() const final
The surface of the GeometricSearchDet.
Definition: ForwardDetLayer.h:31
p1
double p1[4]
Definition: TauolaWrapper.h:89
funct::sqr
Square< F >::type sqr(const F &f)
Definition: Square.h:13
edm::shift
static unsigned int const shift
Definition: LuminosityBlockID.cc:9
CosmicsPD_Skims.radius
tuple radius
Definition: CosmicsPD_Skims.py:135
findQualityFiles.size
tuple size
Write out results.
Definition: findQualityFiles.py:442
p3
double p3[4]
Definition: TauolaWrapper.h:91
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 174 of file PixelTripletNoTipGenerator.cc.

References Exception.

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

Member Data Documentation

float PixelTripletNoTipGenerator::extraHitPhiToleranceForPreFiltering
private

Definition at line 36 of file PixelTripletNoTipGenerator.h.

Referenced by hitTriplets().

float PixelTripletNoTipGenerator::extraHitRPhitolerance
private

Definition at line 35 of file PixelTripletNoTipGenerator.h.

float PixelTripletNoTipGenerator::extraHitRZtolerance
private

Definition at line 34 of file PixelTripletNoTipGenerator.h.

Referenced by hitTriplets().

double PixelTripletNoTipGenerator::theChi2Cut
private

Definition at line 38 of file PixelTripletNoTipGenerator.h.

Referenced by hitTriplets().

double PixelTripletNoTipGenerator::theNSigma
private

Definition at line 37 of file PixelTripletNoTipGenerator.h.

Referenced by hitTriplets().

Generated for CMSSW Reference Manual by doxygen 1.8.5