SeedGeneratorForCRack Class Reference

#include <SeedGeneratorForCRack.h>

Public Types
typedef TrajectoryStateOnSurface	TSOS

Public Member Functions
void	init (const SiStripRecHit2DCollection &collstereo, const SiStripRecHit2DCollection &collrphi, const SiStripMatchedRecHit2DCollection &collmatched, const edm::EventSetup &c)
void	run (TrajectorySeedCollection &, const edm::EventSetup &c)
	SeedGeneratorForCRack (const edm::ParameterSet &conf, edm::ConsumesCollector)
void	seeds (TrajectorySeedCollection &output, const edm::EventSetup &c, const TrackingRegion &region)
virtual	~SeedGeneratorForCRack ()

Private Attributes
std::string	geometry
OrderedHitPairs	HitPairs
edm::ESHandle< MagneticField >	magfield
float	multipleScatteringFactor
GlobalTrackingRegion	region
double	seedMomentum
float	seedpt
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	theMagfieldToken
CosmicHitPairGenerator *	thePairGenerator
PropagatorWithMaterial *	thePropagatorAl
PropagatorWithMaterial *	thePropagatorOp
const edm::ESGetToken< GeometricSearchTracker, TrackerRecoGeometryRecord >	theSearchTrackerToken
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	theTrackerToken
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	theTTopoToken
const edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	theTTRHToken
KFUpdator *	theUpdator
edm::ESHandle< TrackerGeometry >	tracker
const TransientTrackingRecHitBuilder *	TTTRHBuilder

Detailed Description

Definition at line 25 of file SeedGeneratorForCRack.h.

Member Typedef Documentation

TSOS

typedef TrajectoryStateOnSurface SeedGeneratorForCRack::TSOS

Definition at line 27 of file SeedGeneratorForCRack.h.

Constructor & Destructor Documentation

SeedGeneratorForCRack()

SeedGeneratorForCRack::SeedGeneratorForCRack	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector	iCC
	)

Definition at line 29 of file SeedGeneratorForCRack.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), multipleScatteringFactor, seedMomentum, seedpt, and AlCaHLTBitMon_QueryRunRegistry::string.

    : theMagfieldToken(iCC.esConsumes()),
      theTrackerToken(iCC.esConsumes()),
      theSearchTrackerToken(iCC.esConsumes()),
      theTTopoToken(iCC.esConsumes()),
      theTTRHToken(iCC.esConsumes(edm::ESInputTag("", conf.getParameter<std::string>("TTRHBuilder")))),
      region(conf.getParameter<double>("ptMin"),
             conf.getParameter<double>("originRadius"),
             conf.getParameter<double>("originHalfLength"),
             conf.getParameter<double>("originZPosition")) {
  seedpt = conf.getParameter<double>("SeedPt");
  geometry = conf.getUntrackedParameter<std::string>("GeometricStructure", "STANDARD");
  multipleScatteringFactor = conf.getUntrackedParameter<double>("multipleScatteringFactor", 1.0);
  seedMomentum = conf.getUntrackedParameter<double>("SeedMomentum", 1);
}

~SeedGeneratorForCRack()

virtual SeedGeneratorForCRack::~SeedGeneratorForCRack ( )

inlinevirtual

Definition at line 29 of file SeedGeneratorForCRack.h.

{};

Member Function Documentation

init()

void SeedGeneratorForCRack::init	(	const SiStripRecHit2DCollection &	collstereo,
		const SiStripRecHit2DCollection &	collrphi,
		const SiStripMatchedRecHit2DCollection &	collmatched,
		const edm::EventSetup &	c
	)

Definition at line 6 of file SeedGeneratorForCRack.cc.

References alongMomentum, edm::EventSetup::getData(), edm::EventSetup::getHandle(), CosmicHitPairGenerator::hitPairs(), HitPairs, LogDebug, magfield, oppositeToMomentum, region, OrderedHitPairs::size(), theMagfieldToken, thePairGenerator, thePropagatorAl, thePropagatorOp, theSearchTrackerToken, theTrackerToken, theTTopoToken, theTTRHToken, theUpdator, HLT_2023v12_cff::track, tracker, and TTTRHBuilder.

Referenced by CRackSeedGenerator::produce().

                                                              {
  magfield = iSetup.getHandle(theMagfieldToken);
  tracker = iSetup.getHandle(theTrackerToken);
  thePropagatorAl = new PropagatorWithMaterial(alongMomentum, 0.1057, &(*magfield));
  thePropagatorOp = new PropagatorWithMaterial(oppositeToMomentum, 0.1057, &(*magfield));
  theUpdator = new KFUpdator();

  // get the transient builder
  TTTRHBuilder = &iSetup.getData(theTTRHToken);

  GeometricSearchTracker const& track = iSetup.getData(theSearchTrackerToken);
  TrackerTopology const& ttopo = iSetup.getData(theTTopoToken);

  CosmicLayerPairs cosmiclayers(geometry, collrphi, collmatched, track, ttopo);
  thePairGenerator = new CosmicHitPairGenerator(cosmiclayers, *tracker);
  HitPairs.clear();
  thePairGenerator->hitPairs(region, HitPairs);
  LogDebug("CosmicSeedFinder") << "Initialized with " << HitPairs.size() << " hit pairs" << std::endl;
}

run()

void SeedGeneratorForCRack::run	(	TrajectorySeedCollection &	output,
		const edm::EventSetup &	c
	)

Definition at line 45 of file SeedGeneratorForCRack.cc.

References region, seeds(), thePairGenerator, thePropagatorAl, thePropagatorOp, and theUpdator.

Referenced by CRackSeedGenerator::produce().

                                                                                             {
  seeds(output, iSetup, region);
  delete thePairGenerator;
  delete thePropagatorAl;
  delete thePropagatorOp;
  delete theUpdator;
}

seeds()

void SeedGeneratorForCRack::seeds	(	TrajectorySeedCollection &	output,
		const edm::EventSetup &	c,
		const TrackingRegion &	region
	)

Definition at line 53 of file SeedGeneratorForCRack.cc.

References alongMomentum, clone(), HitPairs, hfClusterShapes_cfi::hits, mps_fire::i, TrackerGeometry::idToDet(), SurfaceOrientation::inner, TrajectoryStateOnSurface::isValid(), LogDebug, magfield, multipleScatteringFactor, oppositeToMomentum, SurfaceOrientation::outer, trajectoryStateTransform::persistentState(), PropagatorWithMaterial::propagate(), TrajectoryStateOnSurface::rescaleError(), seedMomentum, OrderedHitPairs::size(), GeomDet::surface(), thePropagatorOp, theUpdator, Surface::toGlobal(), tracker, Vector3DBase< T, FrameTag >::unit(), and KFUpdator::update().

Referenced by run().

                                                                {
  for (unsigned int is = 0; is < HitPairs.size(); is++) {
    GlobalPoint inner = tracker->idToDet((*(HitPairs[is].inner())).geographicalId())
                            ->surface()
                            .toGlobal((*(HitPairs[is].inner())).localPosition());
    GlobalPoint outer = tracker->idToDet((*(HitPairs[is].outer())).geographicalId())
                            ->surface()
                            .toGlobal((*(HitPairs[is].outer())).localPosition());

    LogDebug("CosmicSeedFinder") << "inner point of the seed " << inner << " outer point of the seed " << outer;
    SeedingHitSet::ConstRecHitPointer inrhit = HitPairs[is].inner();
    SeedingHitSet::ConstRecHitPointer outrhit = HitPairs[is].outer();

    edm::OwnVector<TrackingRecHit> hits;
    hits.push_back(HitPairs[is].outer()->hit()->clone());

    for (int i = 0; i < 2; i++) {
      //FIRST STATE IS CALCULATED CONSIDERING THAT THE CHARGE CAN BE POSITIVE OR NEGATIVE
      int predsign = (2 * i) - 1;
      if ((outer.y() - inner.y()) > 0) {
        GlobalVector momentum = GlobalVector(inner - outer);
        momentum = momentum.unit() * seedMomentum;
        GlobalTrajectoryParameters Gtp(inner, momentum, predsign, &(*magfield));
        AlgebraicSymMatrix55 errMatrix = ROOT::Math::SMatrixIdentity();
        TSOS innerState = TSOS(Gtp,
                               CurvilinearTrajectoryError(errMatrix),
                               tracker->idToDet((HitPairs[is].inner()->hit())->geographicalId())->surface());
        const TSOS innerUpdated = theUpdator->update(innerState, *inrhit);
        //Cosmic Seed update inner...
        LogDebug("CosmicSeedFinder") << " FirstTSOS " << innerUpdated;

        //First propagation
        const TSOS outerState = thePropagatorOp->propagate(
            innerUpdated, tracker->idToDet((*(HitPairs[is].outer())).geographicalId())->surface());
        if (outerState.isValid()) {
          LogDebug("CosmicSeedFinder") << "outerState " << outerState;
          TSOS outerUpdated = theUpdator->update(outerState, *outrhit);
          //fudge factor for multiple scattering
          outerUpdated.rescaleError(multipleScatteringFactor);
          if (outerUpdated.isValid()) {
            LogDebug("CosmicSeedFinder") << "outerUpdated " << outerUpdated;

            PTrajectoryStateOnDet PTraj = trajectoryStateTransform::persistentState(
                outerUpdated, (*(HitPairs[is].outer())).geographicalId().rawId());
            output.push_back(TrajectorySeed(PTraj, hits, alongMomentum));

          } else
            edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first update failed ";
        } else
          edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first propagation failed ";

      } else {
        GlobalVector momentum = GlobalVector(outer - inner);
        momentum = momentum.unit() * seedMomentum;
        GlobalTrajectoryParameters Gtp(inner, momentum, predsign, &(*magfield));
        AlgebraicSymMatrix55 errMatrix = ROOT::Math::SMatrixIdentity();
        TSOS innerState = TSOS(Gtp,
                               CurvilinearTrajectoryError(errMatrix),
                               tracker->idToDet((HitPairs[is].inner()->hit())->geographicalId())->surface());
        const TSOS innerUpdated = theUpdator->update(innerState, *inrhit);
        LogDebug("CosmicSeedFinder") << " FirstTSOS " << innerState;

        //First propagation
        const TSOS outerState = thePropagatorOp->propagate(
            innerUpdated, tracker->idToDet((*(HitPairs[is].outer())).geographicalId())->surface());
        if (outerState.isValid()) {
          LogDebug("CosmicSeedFinder") << "outerState " << outerState;
          TSOS outerUpdated = theUpdator->update(outerState, *outrhit);
          //fudge factor for multiple scattering
          outerUpdated.rescaleError(multipleScatteringFactor);
          if (outerUpdated.isValid()) {
            LogDebug("CosmicSeedFinder") << "outerUpdated " << outerUpdated;
            PTrajectoryStateOnDet PTraj = trajectoryStateTransform::persistentState(
                outerUpdated, (*(HitPairs[is].outer())).geographicalId().rawId());
            output.push_back(TrajectorySeed(PTraj, hits, oppositeToMomentum));

          } else
            edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first update failed ";
        } else
          edm::LogWarning("CosmicSeedFinder") << " SeedForCosmics first propagation failed ";
      }
    }
  }
}

Member Data Documentation

geometry

std::string SeedGeneratorForCRack::geometry

private

Definition at line 55 of file SeedGeneratorForCRack.h.

HitPairs

OrderedHitPairs SeedGeneratorForCRack::HitPairs

private

Definition at line 57 of file SeedGeneratorForCRack.h.

Referenced by init(), and seeds().

magfield

edm::ESHandle<MagneticField> SeedGeneratorForCRack::magfield

private

Definition at line 48 of file SeedGeneratorForCRack.h.

Referenced by init(), and seeds().

multipleScatteringFactor

float SeedGeneratorForCRack::multipleScatteringFactor

private

Definition at line 58 of file SeedGeneratorForCRack.h.

Referenced by SeedGeneratorForCRack(), and seeds().

region

GlobalTrackingRegion SeedGeneratorForCRack::region

private

Definition at line 46 of file SeedGeneratorForCRack.h.

Referenced by init(), and run().

seedMomentum

double SeedGeneratorForCRack::seedMomentum

private

Definition at line 59 of file SeedGeneratorForCRack.h.

Referenced by SeedGeneratorForCRack(), and seeds().

seedpt

float SeedGeneratorForCRack::seedpt

private

Definition at line 56 of file SeedGeneratorForCRack.h.

Referenced by SeedGeneratorForCRack().

theMagfieldToken

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> SeedGeneratorForCRack::theMagfieldToken

private

Definition at line 40 of file SeedGeneratorForCRack.h.

Referenced by init().

thePairGenerator

CosmicHitPairGenerator* SeedGeneratorForCRack::thePairGenerator

private

Definition at line 47 of file SeedGeneratorForCRack.h.

Referenced by init(), and run().

thePropagatorAl

PropagatorWithMaterial* SeedGeneratorForCRack::thePropagatorAl

private

Definition at line 52 of file SeedGeneratorForCRack.h.

Referenced by init(), and run().

thePropagatorOp

PropagatorWithMaterial* SeedGeneratorForCRack::thePropagatorOp

private

Definition at line 53 of file SeedGeneratorForCRack.h.

Referenced by init(), run(), and seeds().

theSearchTrackerToken

const edm::ESGetToken<GeometricSearchTracker, TrackerRecoGeometryRecord> SeedGeneratorForCRack::theSearchTrackerToken

private

Definition at line 42 of file SeedGeneratorForCRack.h.

Referenced by init().

theTrackerToken

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SeedGeneratorForCRack::theTrackerToken

private

Definition at line 41 of file SeedGeneratorForCRack.h.

Referenced by init().

theTTopoToken

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SeedGeneratorForCRack::theTTopoToken

private

Definition at line 43 of file SeedGeneratorForCRack.h.

Referenced by init().

theTTRHToken

const edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> SeedGeneratorForCRack::theTTRHToken

private

Definition at line 44 of file SeedGeneratorForCRack.h.

Referenced by init().

theUpdator

KFUpdator* SeedGeneratorForCRack::theUpdator

private

Definition at line 51 of file SeedGeneratorForCRack.h.

Referenced by init(), run(), and seeds().

tracker

edm::ESHandle<TrackerGeometry> SeedGeneratorForCRack::tracker

private

Definition at line 49 of file SeedGeneratorForCRack.h.

Referenced by init(), and seeds().

TTTRHBuilder

const TransientTrackingRecHitBuilder* SeedGeneratorForCRack::TTTRHBuilder

private

Definition at line 54 of file SeedGeneratorForCRack.h.

Referenced by init().