FakeTrackProducer< T > Class Template Reference

Inheritance diagram for FakeTrackProducer< T >:

Public Member Functions
	FakeTrackProducer (const edm::ParameterSet &iConfig)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	~FakeTrackProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
const PTrajectoryStateOnDet &	getState (const TrajectorySeed &seed) const
	Muon selection. More...
const PTrajectoryStateOnDet &	getState (const TrackCandidate &seed) const

Private Attributes
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geometryToken_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_
edm::EDGetTokenT< std::vector< T > >	src_
	Labels for input collections. More...
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	trackerTopoToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

template<class T>
class FakeTrackProducer< T >

Definition at line 39 of file FakeTrackProducers.cc.

Constructor & Destructor Documentation

FakeTrackProducer()

template<typename T >

FakeTrackProducer< T >::FakeTrackProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 61 of file FakeTrackProducers.cc.

    : src_(consumes<std::vector<T>>(iConfig.getParameter<edm::InputTag>("src"))),
      geometryToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>()),
      magFieldToken_(esConsumes<MagneticField, IdealMagneticFieldRecord>()),
      trackerTopoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd>())
//,selector_(iConfig.existsAs<std::string>("cut") ? iConfig.getParameter<std::string>("cut") : "", true)
{
  produces<std::vector<reco::Track>>();
  produces<std::vector<reco::TrackExtra>>();
  produces<edm::OwnVector<TrackingRecHit>>();
}

~FakeTrackProducer()

template<class T >

FakeTrackProducer< T >::~FakeTrackProducer ( )

inlineoverride

Definition at line 42 of file FakeTrackProducers.cc.

{}

Member Function Documentation

getState() [1/2]

template<class T >

const PTrajectoryStateOnDet& FakeTrackProducer< T >::getState ( const TrajectorySeed &  seed ) const

inlineprivate

Muon selection.

Definition at line 56 of file FakeTrackProducers.cc.

References fileCollector::seed.

{ return seed.startingState(); }

getState() [2/2]

template<class T >

const PTrajectoryStateOnDet& FakeTrackProducer< T >::getState ( const TrackCandidate &  seed ) const

inlineprivate

Definition at line 57 of file FakeTrackProducers.cc.

References fileCollector::seed.

{ return seed.trajectoryStateOnDet(); }

produce()

template<typename T >

void FakeTrackProducer< T >::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 74 of file FakeTrackProducers.cc.

References alongMomentum, DMR_cfg::cerr, ALCARECOTkAlJpsiMuMu_cff::charge, PTrajectoryStateOnDet::detId(), TrackingRecHit::geographicalId(), edm::EventSetup::getHandle(), hfClusterShapes_cfi::hits, TrackerGeometry::idToDet(), iEvent, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, TrackingRecHit::localPosition(), eostools::move(), amptDefaultParameters_cff::mu, MillePedeFileConverter_cfg::out, AlCaHLTBitMon_ParallelJobs::p, DetId::rawId(), reco::TrackExtraBase::setHits(), edm::OwnVector< T, P >::size(), TrackRefitter_38T_cff::src, GeomDet::surface(), edm::swap(), GeomDet::toGlobal(), trajectoryStateTransform::transientState(), x, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                {
  using namespace edm;
  using namespace std;

  edm::ESHandle<TrackerGeometry> theGeometry = iSetup.getHandle(geometryToken_);
  edm::ESHandle<MagneticField> theMagField = iSetup.getHandle(magFieldToken_);
  edm::ESHandle<TrackerTopology> httopo = iSetup.getHandle(trackerTopoToken_);
  const TrackerTopology &ttopo = *httopo;

  Handle<vector<T>> src;
  iEvent.getByToken(src_, src);

  unique_ptr<vector<reco::Track>> out(new vector<reco::Track>());
  out->reserve(src->size());
  unique_ptr<vector<reco::TrackExtra>> outEx(new vector<reco::TrackExtra>());
  outEx->reserve(src->size());
  unique_ptr<OwnVector<TrackingRecHit>> outHits(new OwnVector<TrackingRecHit>());

  TrackingRecHitRefProd rHits = iEvent.getRefBeforePut<TrackingRecHitCollection>();
  reco::TrackExtraRefProd rTrackExtras = iEvent.getRefBeforePut<reco::TrackExtraCollection>();
  for (typename vector<T>::const_iterator it = src->begin(), ed = src->end(); it != ed; ++it) {
    const T &mu = *it;
    //if (!selector_(mu)) continue;
    const PTrajectoryStateOnDet &pstate = getState(mu);
    const GeomDet *det = theGeometry->idToDet(DetId(pstate.detId()));
    if (det == nullptr) {
      std::cerr << "ERROR:  bogus detid " << pstate.detId() << std::endl;
      continue;
    }
    TrajectoryStateOnSurface state = trajectoryStateTransform::transientState(pstate, &det->surface(), &*theMagField);
    GlobalPoint gx = state.globalPosition();
    GlobalVector gp = state.globalMomentum();
    reco::Track::Point x(gx.x(), gx.y(), gx.z());
    reco::Track::Vector p(gp.x(), gp.y(), gp.z());
    int charge = state.localParameters().charge();
    out->push_back(reco::Track(1.0, 1.0, x, p, charge, reco::Track::CovarianceMatrix()));
    auto hits = mu.recHits();
    out->back().appendHits(hits.begin(), hits.end(), ttopo);
    // Now Track Extra
    const TrackingRecHit *hit0 = &*hits.begin();
    const TrackingRecHit *hit1 = &*(hits.end() - 1);
    const GeomDet *det0 = theGeometry->idToDet(hit0->geographicalId());
    const GeomDet *det1 = theGeometry->idToDet(hit1->geographicalId());
    if (det0 == nullptr || det1 == nullptr) {
      std::cerr << "ERROR:  bogus detids at beginning or end of range" << std::endl;
      continue;
    }
    GlobalPoint gx0 = det0->toGlobal(hit0->localPosition());
    GlobalPoint gx1 = det1->toGlobal(hit1->localPosition());
    reco::Track::Point x0(gx0.x(), gx0.y(), gx0.z());
    reco::Track::Point x1(gx1.x(), gx1.y(), gx1.z());
    if (x0.R() > x1.R())
      std::swap(x0, x1);
    outEx->push_back(reco::TrackExtra(x1,
                                      p,
                                      true,
                                      x0,
                                      p,
                                      true,
                                      reco::Track::CovarianceMatrix(),
                                      hit0->geographicalId().rawId(),
                                      reco::Track::CovarianceMatrix(),
                                      hit1->geographicalId().rawId(),
                                      alongMomentum));
    out->back().setExtra(reco::TrackExtraRef(rTrackExtras, outEx->size() - 1));
    reco::TrackExtra &ex = outEx->back();
    auto const firstHitIndex = outHits->size();
    for (auto const &it2 : hits) {
      outHits->push_back(it2);
    }
    ex.setHits(rHits, firstHitIndex, outHits->size() - firstHitIndex);
  }

  iEvent.put(std::move(out));
  iEvent.put(std::move(outEx));
  iEvent.put(std::move(outHits));
}

Member Data Documentation

geometryToken_

template<class T >

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> FakeTrackProducer< T >::geometryToken_

private

Definition at line 49 of file FakeTrackProducers.cc.

magFieldToken_

template<class T >

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> FakeTrackProducer< T >::magFieldToken_

private

Definition at line 50 of file FakeTrackProducers.cc.

src_

template<class T >

edm::EDGetTokenT<std::vector<T> > FakeTrackProducer< T >::src_

private

Labels for input collections.

Definition at line 48 of file FakeTrackProducers.cc.

trackerTopoToken_

template<class T >

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> FakeTrackProducer< T >::trackerTopoToken_

private

Definition at line 51 of file FakeTrackProducers.cc.