MuonReSeeder Class Reference

Inheritance diagram for MuonReSeeder:

Public Member Functions
	MuonReSeeder (const edm::ParameterSet &iConfig)
virtual void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
virtual	~MuonReSeeder ()
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducerBase ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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...

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
Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 31 of file MuonReSeeder.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 59 of file MuonReSeeder.cc.

                                                          :
    src_(consumes<edm::View<reco::Muon> >(iConfig.getParameter<edm::InputTag>("src"))),
    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)
{
    produces<std::vector<TrajectorySeed> >(); 
}

virtual MuonReSeeder::~MuonReSeeder ( )

inlinevirtual

Definition at line 34 of file MuonReSeeder.cc.

{ }

Member Function Documentation

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

overridevirtual

Implements edm::stream::EDProducerBase.

Definition at line 71 of file MuonReSeeder.cc.

References funct::abs(), alongMomentum, gather_cfg::cout, debug_, end, reco::LeafCandidate::eta(), edm::OwnVector< T, P >::front(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByToken(), TrajectoryStateOnSurface::globalPosition(), TrackingRecHit::hit(), i, insideOut_, edm::Ref< C, T, F >::isNull(), TrajectoryStateOnSurface::isValid(), TrackingRecHit::isValid(), layersToKeep_, RPCpg::mu, oppositeToMomentum, dbtoconf::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_, relval_parameters_module::step, DetId::subdetId(), reco::Muon::track(), TrackTransformer::transform(), 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
    edm::ESHandle<TrackerTopology> tTopo;
    iSetup.get<IdealGeometryRecord>().get(tTopo);

    auto_ptr<vector<TrajectorySeed> > out(new 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  = 0;
        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 != 0 && 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(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 40 of file MuonReSeeder.cc.

Referenced by produce().