31 bool seedfromoverlaps = ((
std::abs(radius1 - radius2) < 0.1) || (
std::abs(radius3 - radius2) < 0.1)) ?
true :
false;
32 std::vector<const TrackingRecHit *>::const_iterator ohh;
33 std::vector<const TrackingRecHit *>::const_iterator mhh;
34 std::vector<const TrackingRecHit *>::const_iterator ihh;
36 if (!seedfromoverlaps) {
41 float z_diff = mh->globalPosition().z() - oh->globalPosition().z();
42 float midy = mh->globalPosition().y();
43 float outy = oh->globalPosition().y();
44 float midx = mh->globalPosition().x();
45 float outx = oh->globalPosition().x();
46 float dxdy =
std::abs((outx - midx) / (outy - midy));
47 if ((
std::abs(z_diff) < 30) && (outy * midy > 0) && (dxdy < 2)) {
50 float z_diff = mh->globalPosition().z() - ih->globalPosition().z();
51 float inny = ih->globalPosition().y();
52 float innx = ih->globalPosition().x();
53 float dxdy =
std::abs((innx - midx) / (inny - midy));
54 if ((
std::abs(z_diff) < 30) && (inny * midy > 0) && (dxdy < 2) && (!seedfromoverlaps)) {
66 float z_diff = mh->globalPosition().z() - oh->globalPosition().z();
67 float midy = mh->globalPosition().y();
68 float outy = oh->globalPosition().y();
69 float midx = mh->globalPosition().x();
70 float outx = oh->globalPosition().x();
71 float dxdy =
std::abs((outx - midx) / (outy - midy));
72 float DeltaR = oh->globalPosition().perp() - mh->globalPosition().perp();
73 if ((
std::abs(z_diff) < 18) && (
std::abs(oh->globalPosition().phi() - mh->globalPosition().phi()) < 0.05) &&
74 (DeltaR < 0) && (dxdy < 2)) {
77 float z_diff = mh->globalPosition().z() - ih->globalPosition().z();
78 float inny = ih->globalPosition().y();
79 float innx = ih->globalPosition().x();
80 float dxdy =
std::abs((innx - midx) / (inny - midy));
81 if ((
std::abs(z_diff) < 30) && (inny * midy > 0) && (dxdy < 2)) {
const DetLayer * layer() const
const LayerWithHits * theMiddleLayer
Abs< T >::type abs(const T &t)
const LayerWithHits * theInnerLayer
const std::vector< const TrackingRecHit * > & recHits() const
return the recHits of the Layer
const LayerWithHits * theOuterLayer