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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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

Private Attributes
edm::EDGetTokenT< std::vector< T > >	src_
	Labels for input collections. More...

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

Detailed Description

template<class T>
class FakeTrackProducer< T >

Definition at line 38 of file FakeTrackProducers.cc.

Constructor & Destructor Documentation

template<typename T >

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

explicit

Definition at line 59 of file FakeTrackProducers.cc.

    : src_(consumes<std::vector<T>>(iConfig.getParameter<edm::InputTag>("src")))
//,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>>();
}

template<class T >

FakeTrackProducer< T >::~FakeTrackProducer ( )

inlineoverride

Definition at line 41 of file FakeTrackProducers.cc.

References iEvent, and FakeTrackProducer< T >::produce().

{}

Member Function Documentation

template<class T >

TrajectorySeed::range FakeTrackProducer< T >::getHits ( const TrajectorySeed &  seed ) const

inlineprivate

Definition at line 54 of file FakeTrackProducers.cc.

References TrajectorySeed::recHits().

Referenced by FakeTrackProducer< T >::produce().

{ return seed.recHits(); }

template<class T >

TrajectorySeed::range FakeTrackProducer< T >::getHits ( const TrackCandidate &  seed ) const

inlineprivate

Definition at line 55 of file FakeTrackProducers.cc.

References TrackCandidate::recHits().

{ return seed.recHits(); }

template<class T >

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

inlineprivate

Muon selection.

Definition at line 52 of file FakeTrackProducers.cc.

References TrajectorySeed::startingState().

Referenced by FakeTrackProducer< T >::produce().

{ return seed.startingState(); }

template<class T >

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

inlineprivate

Definition at line 53 of file FakeTrackProducers.cc.

References TrackCandidate::trajectoryStateOnDet().

{ return seed.trajectoryStateOnDet(); }

template<typename T >

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

override

Definition at line 69 of file FakeTrackProducers.cc.

References alongMomentum, beam_dqm_sourceclient-live_cfg::cerr, ALCARECOTkAlJpsiMuMu_cff::charge, LocalTrajectoryParameters::charge(), PTrajectoryStateOnDet::detId(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByToken(), FakeTrackProducer< T >::getHits(), edm::Event::getRefBeforePut(), FakeTrackProducer< T >::getState(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), runTauDisplay::gp, hfClusterShapes_cfi::hits, TrackerGeometry::idToDet(), TrajectoryStateOnSurface::localParameters(), TrackingRecHit::localPosition(), eostools::move(), amptDefaultParameters_cff::mu, MillePedeFileConverter_cfg::out, AlCaHLTBitMon_ParallelJobs::p, edm::Event::put(), DetId::rawId(), reco::TrackExtraBase::setHits(), edm::OwnVector< T, P >::size(), TrackRefitter_38T_cff::src, FakeTrackProducer< T >::src_, GeomDet::surface(), std::swap(), GeomDet::toGlobal(), trajectoryStateTransform::transientState(), x, PV3DBase< T, PVType, FrameType >::x(), testProducerWithPsetDescEmpty_cfi::x1, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by FakeTrackProducer< T >::~FakeTrackProducer().

                                                                                {
  using namespace edm;
  using namespace std;

  edm::ESHandle<TrackerGeometry> theGeometry;
  iSetup.get<TrackerDigiGeometryRecord>().get(theGeometry);
  edm::ESHandle<MagneticField> theMagField;
  iSetup.get<IdealMagneticFieldRecord>().get(theMagField);
  edm::ESHandle<TrackerTopology> httopo;
  iSetup.get<TrackerTopologyRcd>().get(httopo);
  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()));
    TrajectorySeed::range hits = getHits(mu);
    out->back().appendHits(hits.first, hits.second, ttopo);
    // Now Track Extra
    const TrackingRecHit *hit0 = &*hits.first;
    const TrackingRecHit *hit1 = &*(hits.second - 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 (OwnVector<TrackingRecHit>::const_iterator it2 = hits.first; it2 != hits.second; ++it2) {
      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

template<class T >

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

private

Labels for input collections.

Definition at line 47 of file FakeTrackProducers.cc.

Referenced by FakeTrackProducer< T >::produce().