MuonReSeeder Class Reference

Inheritance diagram for MuonReSeeder:

Public Member Functions
	MuonReSeeder (const edm::ParameterSet &iConfig)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	~MuonReSeeder () 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 Attributes
bool	debug_
	Dump deug information. More...
bool	insideOut_
	Do inside-out. More...
int	layersToKeep_
	How many hits to keep from the muon trajectory. More...
TrackTransformer	refitter_
	Track Transformer. More...
StringCutObjectSelector < reco::Muon >	selector_
	Muon selection. More...
edm::EDGetTokenT< edm::View < reco::Muon > >	src_
	Labels for input collections. More...
const edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	tTopoToken_

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

Definition at line 30 of file MuonReSeeder.cc.

Constructor & Destructor Documentation

MuonReSeeder::MuonReSeeder ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 58 of file MuonReSeeder.cc.

    : src_(consumes<edm::View<reco::Muon>>(iConfig.getParameter<edm::InputTag>("src"))),
      tTopoToken_(esConsumes()),
      selector_(iConfig.existsAs<std::string>("cut") ? iConfig.getParameter<std::string>("cut") : "", true),
      layersToKeep_(iConfig.getParameter<int32_t>("layersToKeep")),
      insideOut_(iConfig.getParameter<bool>("insideOut")),
      debug_(iConfig.getUntrackedParameter<bool>("debug", false)),
      refitter_(iConfig, consumesCollector()) {
  produces<std::vector<TrajectorySeed>>();
}

MuonReSeeder::~MuonReSeeder ( )

inlineoverride

Definition at line 33 of file MuonReSeeder.cc.

{}

Member Function Documentation

void MuonReSeeder::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 69 of file MuonReSeeder.cc.

References funct::abs(), alongMomentum, gather_cfg::cout, debug_, dataset::end, reco::LeafCandidate::eta(), edm::OwnVector< T, P >::front(), TrackingRecHit::geographicalId(), edm::Event::getByToken(), edm::EventSetup::getData(), TrajectoryStateOnSurface::globalPosition(), TrackingRecHit::hit(), mps_fire::i, insideOut_, edm::Ref< C, T, F >::isNull(), TrajectoryStateOnSurface::isValid(), TrackingRecHit::isValid(), layersToKeep_, eostools::move(), RPCpg::mu, oppositeToMomentum, submitPVResolutionJobs::out, PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), edm::OwnVector< T, P >::push_back(), edm::Event::put(), DetId::rawId(), TrajectoryMeasurement::recHit(), refitter_, fileCollector::seed, selector_, TrackTransformer::setServices(), alcazmumu_cfi::src, src_, DetId::subdetId(), reco::Muon::track(), TrackTransformer::transform(), tTopoToken_, TrajectoryMeasurement::updatedState(), and PV3DBase< T, PVType, FrameType >::z().

                                                                        {
  using namespace edm;
  using namespace std;

  refitter_.setServices(iSetup);

  Handle<View<reco::Muon>> src;
  iEvent.getByToken(src_, src);

  //Retrieve tracker topology from geometry
  const TrackerTopology &tTopo = iSetup.getData(tTopoToken_);

  auto out = std::make_unique<std::vector<TrajectorySeed>>();
  unsigned int nsrc = src->size();
  out->reserve(nsrc);

  for (View<reco::Muon>::const_iterator it = src->begin(), ed = src->end(); it != ed; ++it) {
    const reco::Muon &mu = *it;
    if (mu.track().isNull() || !selector_(mu))
      continue;
    std::vector<Trajectory> traj = refitter_.transform(*mu.track());
    if (traj.size() != 1)
      continue;
    edm::OwnVector<TrackingRecHit> seedHits;
    const std::vector<TrajectoryMeasurement> &tms = traj.front().measurements();
    TrajectoryStateOnSurface tsos;
    const TrackingRecHit *hit = nullptr;
    bool fromInside = (insideOut_ == (traj.front().direction() == alongMomentum));
    if (debug_) {
      std::cout << "Considering muon of pt " << mu.pt() << ", eta = " << mu.eta() << ", phi = " << mu.phi()
                << std::endl;
      std::cout << "Trajectory is " << (traj.front().direction() == alongMomentum ? "along" : "opposite")
                << " to momentum, so will start from " << (fromInside ? "inside" : "outside") << std::endl;
    }
    const TrajectoryMeasurement &tin = (fromInside ? tms.front() : tms.back());
    const TrajectoryMeasurement &tou = (fromInside ? tms.front() : tms.back());
    if (debug_) {
      std::cout << "IN state: subdetector   = " << tin.recHit()->geographicalId().subdetId() << std::endl;
      std::cout << "          global pos Rho  " << tin.updatedState().globalPosition().perp() << ", Z   "
                << tin.updatedState().globalPosition().z() << std::endl;
      std::cout << "OU state: subdetector   = " << tou.recHit()->geographicalId().subdetId() << std::endl;
      std::cout << "          global pos Rho  " << tou.updatedState().globalPosition().perp() << ", Z   "
                << tou.updatedState().globalPosition().z() << std::endl;
    }
    int lastSubdet = 0, lastLayer = -1;
    for (int i = (fromInside ? 0 : tms.size() - 1),
             end = (fromInside ? tms.size() : -1),
             step = (fromInside ? +1 : -1),
             taken = 0;
         (end - i) * step > 0;
         i += step) {
      const TrackingRecHit *lastHit = hit;
      hit = tms[i].recHit()->hit();
      if (debug_)
        std::cout << "  considering hit " << i << ": rechit on " << (hit ? hit->geographicalId().rawId() : -1)
                  << std::endl;
      if (!hit)
        continue;
      int subdet = hit->geographicalId().subdetId();
      int lay = tTopo.layer(hit->geographicalId());
      if (subdet != lastSubdet || lay != lastLayer) {
        // I'm on a new layer
        if (lastHit != nullptr && taken == layersToKeep_) {
          // I've had enough layers, I can stop here
          hit = lastHit;
          break;
        }
        lastSubdet = subdet;
        lastLayer = lay;
        taken++;
      }
      seedHits.push_back(*hit);
      tsos = tms[i].updatedState().isValid()
                 ? tms[i].updatedState()
                 : (abs(i - end) < abs(i) ? tms[i].forwardPredictedState() : tms[i].backwardPredictedState());
      if (debug_) {
        std::cout << "     hit  : subdetector   = " << tms[i].recHit()->geographicalId().subdetId() << std::endl;
        if (hit->isValid()) {
          std::cout << "            global pos Rho  " << tms[i].recHit()->globalPosition().perp() << ", Z   "
                    << tms[i].recHit()->globalPosition().z() << std::endl;
        } else {
          std::cout << "            invalid tracking rec hit, so no global position" << std::endl;
        }
        std::cout << "     state: global pos Rho  " << tsos.globalPosition().perp() << ", Z   "
                  << tsos.globalPosition().z() << std::endl;
      }
    }
    if (!tsos.isValid())
      continue;
    PTrajectoryStateOnDet const &PTraj = trajectoryStateTransform::persistentState(tsos, hit->geographicalId().rawId());
    TrajectorySeed seed(PTraj, std::move(seedHits), insideOut_ ? alongMomentum : oppositeToMomentum);
    out->push_back(seed);
  }

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

Member Data Documentation

bool MuonReSeeder::debug_

private

Dump deug information.

Definition at line 52 of file MuonReSeeder.cc.

Referenced by produce().

bool MuonReSeeder::insideOut_

private

Do inside-out.

Definition at line 49 of file MuonReSeeder.cc.

Referenced by produce().

int MuonReSeeder::layersToKeep_

private

How many hits to keep from the muon trajectory.

Definition at line 46 of file MuonReSeeder.cc.

Referenced by produce().

TrackTransformer MuonReSeeder::refitter_

private

Track Transformer.

Definition at line 55 of file MuonReSeeder.cc.

Referenced by produce().

StringCutObjectSelector<reco::Muon> MuonReSeeder::selector_

private

Muon selection.

Definition at line 43 of file MuonReSeeder.cc.

Referenced by produce().

edm::EDGetTokenT<edm::View<reco::Muon> > MuonReSeeder::src_

private

Labels for input collections.

Definition at line 39 of file MuonReSeeder.cc.

Referenced by produce().

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> MuonReSeeder::tTopoToken_

private

Definition at line 40 of file MuonReSeeder.cc.

Referenced by produce().