#include <DTRecSegment2DExtendedProducer.h>

Inheritance diagram for DTRecSegment2DExtendedProducer:

Public Member Functions
	DTRecSegment2DExtendedProducer (const edm::ParameterSet &)
	Constructor. More...

void	produce (edm::Event &event, const edm::EventSetup &setup) override
	The method which produces the 2D-segments. More...

	~DTRecSegment2DExtendedProducer () override
	Destructor. More...

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
bool	debug

edm::EDGetTokenT< DTRecClusterCollection >	recClusToken_

edm::EDGetTokenT< DTRecHitCollection >	recHits1DToken_

DTCombinatorialExtendedPatternReco *	theAlgo

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Producer for DT segment in one projection.

Author: Stefano Lacaprara - INFN Legnaro stefa.nosp@m.no.l.nosp@m.acapr.nosp@m.ara@.nosp@m.pd.in.nosp@m.fn.i.nosp@m.t; Riccardo Bellan - INFN TO ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Definition at line 32 of file DTRecSegment2DExtendedProducer.h.

Constructor & Destructor Documentation

DTRecSegment2DExtendedProducer::DTRecSegment2DExtendedProducer ( const edm::ParameterSet & pset )

Constructor.

Definition at line 35 of file DTRecSegment2DExtendedProducer.cc.

References gather_cfg::cout, debug, edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

                                                                                           {
   // Set verbose output
   debug = pset.getUntrackedParameter<bool>("debug"); 
 
   // the name of the 1D rec hits collection
   recHits1DToken_ = consumes<DTRecHitCollection>(pset.getParameter<InputTag>("recHits1DLabel"));
   recClusToken_ = consumes<DTRecClusterCollection>(pset.getParameter<InputTag>("recClusLabel"));
 
   if(debug)
     cout << "[DTRecSegment2DExtendedProducer] Constructor called" << endl;
 
   produces<DTRecSegment2DCollection>();
 
   // Get the concrete reconstruction algo from the factory
   theAlgo = new DTCombinatorialExtendedPatternReco(pset.getParameter<ParameterSet>("Reco2DAlgoConfig"));
 }

DTRecSegment2DExtendedProducer::~DTRecSegment2DExtendedProducer ( )

override

Destructor.

Definition at line 53 of file DTRecSegment2DExtendedProducer.cc.

References gather_cfg::cout, and debug.

                                                                 {
   if(debug)
     cout << "[DTRecSegment2DExtendedProducer] Destructor called" << endl;
   delete theAlgo;
 }

Member Function Documentation

void DTRecSegment2DExtendedProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	setup
	)

override

The method which produces the 2D-segments.

Definition at line 60 of file DTRecSegment2DExtendedProducer.cc.

References edm::OwnVector< T, P >::begin(), popcon2dropbox::copy(), gather_cfg::cout, debug, edm::OwnVector< T, P >::empty(), edm::OwnVector< T, P >::end(), edm::EventSetup::get(), DTRangeMapAccessor::layersBySuperLayer(), eostools::move(), edm::OwnVector< T, P >::size(), DTGeometry::superLayer(), and DTLayerId::superlayerId().

                                                            {
   if(debug)
     cout << "[DTRecSegment2DExtendedProducer] produce called" << endl;
   // Get the DT Geometry
   ESHandle<DTGeometry> dtGeom;
   setup.get<MuonGeometryRecord>().get(dtGeom);
 
   theAlgo->setES(setup);
   
   // Get the 1D rechits from the event
   Handle<DTRecHitCollection> allHits; 
   event.getByToken(recHits1DToken_, allHits);
 
   // Get the 1D clusters from the event
   Handle<DTRecClusterCollection> dtClusters; 
   event.getByToken(recClusToken_, dtClusters);
   theAlgo->setClusters(vector<DTSLRecCluster>(dtClusters->begin(),
                                               dtClusters->end()));
 
   // Create the pointer to the collection which will store the rechits
   auto segments = std::make_unique<DTRecSegment2DCollection>();
 
   // Iterate through all hit collections ordered by LayerId
   DTRecHitCollection::id_iterator dtLayerIt;
   DTSuperLayerId oldSlId;
   for (dtLayerIt = allHits->id_begin(); dtLayerIt != allHits->id_end(); ++dtLayerIt){
     // The layerId
     DTLayerId layerId = (*dtLayerIt);
     const DTSuperLayerId SLId = layerId.superlayerId();
     if (SLId==oldSlId) continue; // I'm on the same SL as before
     oldSlId = SLId;
 
     if(debug) cout <<"Reconstructing the 2D segments in "<< SLId << endl;
 
     const DTSuperLayer* sl = dtGeom->superLayer(SLId);
 
     // Get all the rec hit in the same superLayer in which layerId relies 
     DTRecHitCollection::range range =
       allHits->get(DTRangeMapAccessor::layersBySuperLayer(SLId));
 
     // Fill the vector with the 1D RecHit
     vector<DTRecHit1DPair> pairs(range.first,range.second);
 
     if(debug) cout << "Number of 1D-RecHit pairs " << pairs.size() << endl;
 
     //if(debug) cout << "Start the 2D-segments Reco "<< endl;
     OwnVector<DTSLRecSegment2D> segs = theAlgo->reconstruct(sl, pairs);
     if(debug) {
       cout << "Number of Reconstructed segments: " << segs.size() << endl;
       copy(segs.begin(), segs.end(),
            ostream_iterator<DTSLRecSegment2D>(cout, "\n"));
     }
 
     if (!segs.empty() )
       segments->put(SLId, segs.begin(),segs.end());
   }
   event.put(std::move(segments));
 }