107 if (
positiveYOnly && (outrhit->globalPosition().y() < 0 || innrhit->globalPosition().y() < 0 ||
108 outrhit->globalPosition().y() < innrhit->globalPosition().y()))
110 if (
negativeYOnly && (outrhit->globalPosition().y() > 0 || innrhit->globalPosition().y() > 0 ||
111 outrhit->globalPosition().y() > innrhit->globalPosition().y()))
119 float ch = helix.stateAtVertex().charge();
120 float mom2 = gv.
x() * gv.x() + gv.y() * gv.y() + gv.z() * gv.z();
135 LogDebug(
"CosmicSeedFinder") <<
"outerState " << outerState;
138 LogDebug(
"CosmicSeedFinder") <<
"outerUpdated " << outerUpdated;
146 edm::LogError(
"TooManySeeds") <<
"Found too many seeds, bailing out.\n";
156 LogDebug(
"CosmicSeedFinder") <<
"outerState " << outerState;
159 LogDebug(
"CosmicSeedFinder") <<
"outerUpdated " << outerUpdated;
167 edm::LogError(
"TooManySeeds") <<
"Found too many seeds, bailing out.\n";
184 LogDebug(
"CosmicSeedFinder") <<
"inner point of the seed " <<
inner <<
" outer point of the seed " <<
outer;
188 if (
positiveYOnly && (outrhit->globalPosition().y() < 0 || innrhit->globalPosition().y() < 0 ||
189 outrhit->globalPosition().y() < innrhit->globalPosition().y()))
191 if (
negativeYOnly && (outrhit->globalPosition().y() > 0 || innrhit->globalPosition().y() > 0 ||
192 outrhit->globalPosition().y() > innrhit->globalPosition().y()))
200 for (
int i = 0;
i < 2;
i++) {
202 int predsign = (2 *
i) - 1;
209 LogDebug(
"CosmicSeedFinder") <<
" FirstTSOS " << CosmicSeed;
214 LogDebug(
"CosmicSeedFinder") <<
"outerState " << outerState;
217 LogDebug(
"CosmicSeedFinder") <<
"outerUpdated " << outerUpdated;
220 outerUpdated, (*(
HitPairs[is].
outer())).geographicalId().rawId());
225 edm::LogError(
"TooManySeeds") <<
"Found too many seeds, bailing out.\n";
231 edm::LogWarning(
"CosmicSeedFinder") <<
" SeedForCosmics first update failed ";
233 edm::LogWarning(
"CosmicSeedFinder") <<
" SeedForCosmics first propagation failed ";
239 LogDebug(
"CosmicSeedFinder") <<
" FirstTSOS " << CosmicSeed;
244 LogDebug(
"CosmicSeedFinder") <<
"outerState " << outerState;
247 LogDebug(
"CosmicSeedFinder") <<
"outerUpdated " << outerUpdated;
250 outerUpdated, (*(
HitPairs[is].
outer())).geographicalId().rawId());
255 edm::LogError(
"TooManySeeds") <<
"Found too many seeds, bailing out.\n";
261 edm::LogWarning(
"CosmicSeedFinder") <<
" SeedForCosmics first update failed ";
263 edm::LogWarning(
"CosmicSeedFinder") <<
" SeedForCosmics first propagation failed ";
edm::ESHandle< TrackerGeometry > tracker
constexpr bool isNotFinite(T x)
ROOT::Math::SMatrixIdentity AlgebraicMatrixID
Log< level::Error, false > LogError
PropagatorWithMaterial * thePropagatorOp
OrderedHitTriplets HitTriplets
BaseTrackerRecHit const * ConstRecHitPointer
PropagatorWithMaterial * thePropagatorAl
TrajectoryStateOnSurface update(const TrajectoryStateOnSurface &, const TrackingRecHit &) const override
const TrackerGeomDet * idToDet(DetId) const override
unsigned int size() const override
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
const Plane & surface() const
The nominal surface of the GeomDet.
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
ROOT::Math::SMatrix< double, 5, 5, ROOT::Math::MatRepSym< double, 5 > > AlgebraicSymMatrix55
TEveGeoShape * clone(const TEveElement *element, TEveElement *parent)
TrajectoryStateOnSurface propagate(STA const &state, SUR const &surface) const
int nominalValue() const
The nominal field value for this map in kGauss.
Log< level::Warning, false > LogWarning
unsigned int size() const override
edm::ESHandle< MagneticField > magfield