SeedTransformer Class Reference

#include <SeedTransformer.h>

Public Member Functions
std::vector< Trajectory >	seedTransform (const TrajectorySeed &) const
	Performs the fit. More...
	SeedTransformer (const edm::ParameterSet &, edm::ConsumesCollector)
	Constructor. More...
TrajectoryStateOnSurface	seedTransientState (const TrajectorySeed &) const
void	setServices (const edm::EventSetup &)
	Set the services needed by the SeedTransformer. More...
virtual	~SeedTransformer ()
	Destructor. More...

Private Attributes
double	errorRescale
	Error rescale factor. More...
unsigned int	nMinRecHits
	Minimum number of RecHits required to perform the fit. More...
const TrajectoryFitter *	theFitter
edm::ESGetToken< TrajectoryFitter, TrajectoryFitter::Record >	theFitterToken
const MagneticField *	theMagneticField
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	theMagneticFieldToken
const TransientTrackingRecHitBuilder *	theMuonRecHitBuilder
edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	theMuonRecHitBuilderToken
const Propagator *	thePropagator
edm::ESGetToken< Propagator, TrackingComponentsRecord >	thePropagatorToken
const GlobalTrackingGeometry *	theTrackingGeometry
edm::ESGetToken< GlobalTrackingGeometry, GlobalTrackingGeometryRecord >	theTrackingGeometryToken
bool	useSubRecHits

Detailed Description

Description: this class takes a TrajectorySeed, fits its RecHits and returns a vector of Trajectories. If the fit fails, the returned vector is empty.

Author

D. Trocino - University and INFN Torino

Definition at line 35 of file SeedTransformer.h.

Constructor & Destructor Documentation

SeedTransformer()

SeedTransformer::SeedTransformer	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector	iC
	)

Constructor.

Definition at line 47 of file SeedTransformer.cc.

References edm::ConsumesCollector::esConsumes(), edm::ParameterSet::getParameter(), LogTrace, and HLT_2023v12_cff::propagatorName.

                                                                                  {
  LogTrace("Reco|TrackingTools|SeedTransformer") << "SeedTransformer constructor called." << endl << endl;

  theTrackingGeometryToken = iC.esConsumes();
  theMagneticFieldToken = iC.esConsumes();
  auto fitterName = iConfig.getParameter<string>("Fitter");
  theFitterToken = iC.esConsumes(edm::ESInputTag("", fitterName));
  auto muonRecHitBuilderName = iConfig.getParameter<string>("MuonRecHitBuilder");
  theMuonRecHitBuilderToken = iC.esConsumes(edm::ESInputTag("", muonRecHitBuilderName));
  auto propagatorName = iConfig.getParameter<string>("Propagator");
  thePropagatorToken = iC.esConsumes(edm::ESInputTag("", propagatorName));

  nMinRecHits = iConfig.getParameter<unsigned int>("NMinRecHits");
  errorRescale = iConfig.getParameter<double>("RescaleError");
  useSubRecHits = iConfig.getParameter<bool>("UseSubRecHits");
}

~SeedTransformer()

SeedTransformer::~SeedTransformer ( )

virtual

Destructor.

Definition at line 64 of file SeedTransformer.cc.

References LogTrace.

                                  {
  LogTrace("Reco|TrackingTools|SeedTransformer") << "SeedTransformer destructor called." << endl << endl;
}

Member Function Documentation

seedTransform()

vector< Trajectory > SeedTransformer::seedTransform

(

const TrajectorySeed &

aSeed

)

const

Performs the fit.

Definition at line 76 of file SeedTransformer.cc.

References MuonTransientTrackingRecHitBreaker::breakInSubRecHits(), filterCSVwithJSON::copy, TrajectoryStateOnSurface::isValid(), LogTrace, metname, TrajectorySeed::nHits(), rpcPointValidation_cfi::recHit, FastTrackerRecHitMaskProducer_cfi::recHits, TrajectorySeed::recHits(), and TrajectoryStateOnSurface::rescaleError().

Referenced by DirectMuonTrajectoryBuilder::trajectories().

                                                                                   {
  const string metname = "Reco|TrackingTools|SeedTransformer";

  LogTrace(metname) << " Number of valid RecHits:      " << aSeed.nHits() << endl;

  if (aSeed.nHits() < nMinRecHits) {
    LogTrace(metname) << "    --- Too few RecHits, no refit performed! ---" << endl;
    return vector<Trajectory>();
  }

  TrajectoryStateOnSurface aTSOS(seedTransientState(aSeed));

  // Rescale errors before refit, not to bias the result
  aTSOS.rescaleError(errorRescale);

  vector<TransientTrackingRecHit::ConstRecHitPointer> recHits;

  for (auto const& recHit : aSeed.recHits()) {
    if (recHit.isValid()) {
      TransientTrackingRecHit::ConstRecHitPointer ttrh(theMuonRecHitBuilder->build(&recHit));

      if (useSubRecHits) {
        TransientTrackingRecHit::ConstRecHitContainer subHits =
            MuonTransientTrackingRecHitBreaker::breakInSubRecHits(ttrh, 2);
        copy(subHits.begin(), subHits.end(), back_inserter(recHits));
      } else {
        recHits.push_back(ttrh);
      }
    }
  }

  TrajectoryStateOnSurface aInitTSOS = thePropagator->propagate(aTSOS, recHits.front()->det()->surface());

  if (!aInitTSOS.isValid()) {
    LogTrace(metname) << "    --- Initial state for refit not valid! ---" << endl;
    return vector<Trajectory>();
  }

  vector<Trajectory> refittedSeed = theFitter->fit(aSeed, recHits, aInitTSOS);

  if (refittedSeed.empty()) {
    LogTrace(metname) << "    --- Seed fit failed! ---" << endl;
    return vector<Trajectory>();
  }

  else if (!refittedSeed.front().isValid()) {
    LogTrace(metname) << "    --- Seed fitted, but trajectory not valid! ---" << endl;
    return vector<Trajectory>();
  }

  else
    LogTrace(metname) << "    +++ Seed fit succeded! +++" << endl;

  return refittedSeed;
}

seedTransientState()

TrajectoryStateOnSurface SeedTransformer::seedTransientState

(

const TrajectorySeed &

tmpSeed

)

const

Definition at line 132 of file SeedTransformer.cc.

References PTrajectoryStateOnDet::detId(), TrajectorySeed::startingState(), and trajectoryStateTransform::transientState().

                                                                                                {
  PTrajectoryStateOnDet tmpTSOD = tmpSeed.startingState();
  DetId tmpDetId(tmpTSOD.detId());
  const GeomDet* tmpGeomDet = theTrackingGeometry->idToDet(tmpDetId);

  TrajectoryStateOnSurface tmpTSOS =
      trajectoryStateTransform::transientState(tmpTSOD, &(tmpGeomDet->surface()), &(*theMagneticField));

  return tmpTSOS;
}

setServices()

void SeedTransformer::setServices

(

const edm::EventSetup &

iSetup

)

Set the services needed by the SeedTransformer.

Definition at line 68 of file SeedTransformer.cc.

References edm::EventSetup::getData().

Referenced by DirectMuonTrajectoryBuilder::trajectories().

                                                          {
  theTrackingGeometry = &iSetup.getData(theTrackingGeometryToken);
  theMagneticField = &iSetup.getData(theMagneticFieldToken);
  theFitter = &iSetup.getData(theFitterToken);
  theMuonRecHitBuilder = &iSetup.getData(theMuonRecHitBuilderToken);
  thePropagator = &iSetup.getData(thePropagatorToken);
}

Member Data Documentation

errorRescale

double SeedTransformer::errorRescale

private

Error rescale factor.

Definition at line 69 of file SeedTransformer.h.

nMinRecHits

unsigned int SeedTransformer::nMinRecHits

private

Minimum number of RecHits required to perform the fit.

Definition at line 66 of file SeedTransformer.h.

theFitter

const TrajectoryFitter* SeedTransformer::theFitter

private

Definition at line 55 of file SeedTransformer.h.

theFitterToken

edm::ESGetToken<TrajectoryFitter, TrajectoryFitter::Record> SeedTransformer::theFitterToken

private

Definition at line 61 of file SeedTransformer.h.

theMagneticField

const MagneticField* SeedTransformer::theMagneticField

private

Definition at line 54 of file SeedTransformer.h.

theMagneticFieldToken

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> SeedTransformer::theMagneticFieldToken

private

Definition at line 60 of file SeedTransformer.h.

theMuonRecHitBuilder

const TransientTrackingRecHitBuilder* SeedTransformer::theMuonRecHitBuilder

private

Definition at line 56 of file SeedTransformer.h.

theMuonRecHitBuilderToken

edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> SeedTransformer::theMuonRecHitBuilderToken

private

Definition at line 62 of file SeedTransformer.h.

thePropagator

const Propagator* SeedTransformer::thePropagator

private

Definition at line 57 of file SeedTransformer.h.

thePropagatorToken

edm::ESGetToken<Propagator, TrackingComponentsRecord> SeedTransformer::thePropagatorToken

private

Definition at line 63 of file SeedTransformer.h.

theTrackingGeometry

const GlobalTrackingGeometry* SeedTransformer::theTrackingGeometry

private

Definition at line 53 of file SeedTransformer.h.

theTrackingGeometryToken

edm::ESGetToken<GlobalTrackingGeometry, GlobalTrackingGeometryRecord> SeedTransformer::theTrackingGeometryToken

private

Definition at line 59 of file SeedTransformer.h.

useSubRecHits

bool SeedTransformer::useSubRecHits

private

Definition at line 71 of file SeedTransformer.h.