|
|
#include <SeedGeneratorForCosmics.h>
Definition at line 26 of file SeedGeneratorForCosmics.h.
◆ TSOS
◆ SeedGeneratorForCosmics()
◆ ~SeedGeneratorForCosmics()
virtual SeedGeneratorForCosmics::~SeedGeneratorForCosmics |
( |
| ) |
|
|
inlinevirtual |
◆ init()
Definition at line 7 of file SeedGeneratorForCosmics.cc.
28 cosmiclayers.init(collstereo, collrphi, collmatched, iSetup);
36 cosmiclayers2.
init(collstereo, collrphi, collmatched,
geometry, iSetup);
References alongMomentum, builderName, edm::EventSetup::get(), get, CosmicHitPairGenerator::hitPairs(), HitPairs, hitsforseeds, CosmicHitTripletGenerator::hitTriplets(), HitTriplets, CosmicLayerPairs::init(), CosmicLayerTriplets::init(), LogDebug, magfield, oppositeToMomentum, edm::ESHandle< T >::product(), region, thePairGenerator, thePropagatorAl, thePropagatorOp, theTripletGenerator, theUpdator, tracker, and TTTRHBuilder.
Referenced by CosmicSeedGenerator::produce().
◆ run()
◆ seeds()
Definition at line 74 of file SeedGeneratorForCosmics.cc.
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 Mom =
sqrt(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;
189 if (
positiveYOnly && (outrhit->globalPosition().y() < 0 || innrhit->globalPosition().y() < 0 ||
190 outrhit->globalPosition().y() < innrhit->globalPosition().y()))
192 if (
negativeYOnly && (outrhit->globalPosition().y() > 0 || innrhit->globalPosition().y() > 0 ||
193 outrhit->globalPosition().y() > innrhit->globalPosition().y()))
201 for (
int i = 0;
i < 2;
i++) {
203 int predsign = (2 *
i) - 1;
210 LogDebug(
"CosmicSeedFinder") <<
" FirstTSOS " << CosmicSeed;
215 LogDebug(
"CosmicSeedFinder") <<
"outerState " << outerState;
218 LogDebug(
"CosmicSeedFinder") <<
"outerUpdated " << outerUpdated;
221 outerUpdated, (*(
HitPairs[is].
outer())).geographicalId().rawId());
226 edm::LogError(
"TooManySeeds") <<
"Found too many seeds, bailing out.\n";
232 edm::LogWarning(
"CosmicSeedFinder") <<
" SeedForCosmics first update failed ";
234 edm::LogWarning(
"CosmicSeedFinder") <<
" SeedForCosmics first propagation failed ";
240 LogDebug(
"CosmicSeedFinder") <<
" FirstTSOS " << CosmicSeed;
245 LogDebug(
"CosmicSeedFinder") <<
"outerState " << outerState;
248 LogDebug(
"CosmicSeedFinder") <<
"outerUpdated " << outerUpdated;
251 outerUpdated, (*(
HitPairs[is].
outer())).geographicalId().rawId());
256 edm::LogError(
"TooManySeeds") <<
"Found too many seeds, bailing out.\n";
262 edm::LogWarning(
"CosmicSeedFinder") <<
" SeedForCosmics first update failed ";
264 edm::LogWarning(
"CosmicSeedFinder") <<
" SeedForCosmics first propagation failed ";
References alongMomentum, clone(), HitPairs, hfClusterShapes_cfi::hits, HitTriplets, mps_fire::i, TrackerGeometry::idToDet(), SurfaceOrientation::inner, CommonMethods::isnan(), TrajectoryStateOnSurface::isValid(), LogDebug, mag(), magfield, maxSeeds_, negativeYOnly, MagneticField::nominalValue(), oppositeToMomentum, SurfaceOrientation::outer, convertSQLitetoXML_cfg::output, trajectoryStateTransform::persistentState(), positiveYOnly, PropagatorWithMaterial::propagate(), seedpt, OrderedHitPairs::size(), OrderedHitTriplets::size(), mathSSE::sqrt(), GeomDet::surface(), thePropagatorAl, thePropagatorOp, theUpdator, Surface::toGlobal(), tracker, and KFUpdator::update().
Referenced by run().
◆ builderName
std::string SeedGeneratorForCosmics::builderName |
|
private |
◆ geometry
std::string SeedGeneratorForCosmics::geometry |
|
private |
◆ HitPairs
◆ hitsforseeds
std::string SeedGeneratorForCosmics::hitsforseeds |
|
private |
◆ HitTriplets
◆ magfield
◆ maxSeeds_
int32_t SeedGeneratorForCosmics::maxSeeds_ |
|
private |
◆ negativeYOnly
bool SeedGeneratorForCosmics::negativeYOnly |
|
private |
◆ positiveYOnly
bool SeedGeneratorForCosmics::positiveYOnly |
|
private |
◆ region
◆ seedpt
float SeedGeneratorForCosmics::seedpt |
|
private |
◆ thePairGenerator
◆ thePropagatorAl
◆ thePropagatorOp
◆ theTripletGenerator
◆ theUpdator
KFUpdator* SeedGeneratorForCosmics::theUpdator |
|
private |
◆ tracker
◆ TTTRHBuilder
const TrackerGeomDet * idToDet(DetId) const override
T const * product() const
edm::ESHandle< MagneticField > magfield
TrajectoryStateOnSurface update(const TrajectoryStateOnSurface &, const TrackingRecHit &) const override
TrajectoryStateOnSurface propagate(STA const &state, SUR const &surface) const
edm::ESHandle< TrackerGeometry > tracker
BaseTrackerRecHit const * ConstRecHitPointer
T getUntrackedParameter(std::string const &, T const &) const
Log< level::Info, false > LogInfo
Log< level::Warning, false > LogWarning
ROOT::Math::SMatrixIdentity AlgebraicMatrixID
CosmicHitPairGenerator * thePairGenerator
const Plane & surface() const
The nominal surface of the GeomDet.
PropagatorWithMaterial * thePropagatorAl
GlobalTrackingRegion region
OrderedHitTriplets HitTriplets
PropagatorWithMaterial * thePropagatorOp
TEveGeoShape * clone(const TEveElement *element, TEveElement *parent)
CosmicHitTripletGenerator * theTripletGenerator
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
int nominalValue() const
The nominal field value for this map in kGauss.
bool seeds(TrajectorySeedCollection &output, const edm::EventSetup &c, const TrackingRegion ®ion)
void init(const SiStripRecHit2DCollection &collstereo, const SiStripRecHit2DCollection &collrphi, const SiStripMatchedRecHit2DCollection &collmatched, std::string geometry, const edm::EventSetup &iSetup)
unsigned int size() const override
Log< level::Error, false > LogError
void hitTriplets(const TrackingRegion ®, OrderedHitTriplets &prs, const edm::EventSetup &iSetup)
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
T getParameter(std::string const &) const
unsigned int size() const override
ROOT::Math::SMatrix< double, 5, 5, ROOT::Math::MatRepSym< double, 5 > > AlgebraicSymMatrix55
void hitPairs(const TrackingRegion ®, OrderedHitPairs &prs, const edm::EventSetup &iSetup)
form base class
const TransientTrackingRecHitBuilder * TTTRHBuilder