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
std::vector< ConsumesInfo >	consumesInfo () const

	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

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) 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.

34 { }

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

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation