#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::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

const edm::ESGetToken< DTGeometry, MuonGeometryRecord >	dtGeomToken_

edm::EDGetTokenT< DTRecClusterCollection >	recClusToken_

edm::EDGetTokenT< DTRecHitCollection >	recHits1DToken_

DTCombinatorialExtendedPatternReco *	theAlgo

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

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 34 of file DTRecSegment2DExtendedProducer.h.

Constructor & Destructor Documentation

◆ DTRecSegment2DExtendedProducer()

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

Constructor.

Definition at line 36 of file DTRecSegment2DExtendedProducer.cc.

References gather_cfg::cout, debug, muonDTDigis_cfi::pset, recClusToken_, recHits1DToken_, and theAlgo.

     : dtGeomToken_(esConsumes()) {
   // 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"), consumesCollector());
 }

◆ ~DTRecSegment2DExtendedProducer()

DTRecSegment2DExtendedProducer::~DTRecSegment2DExtendedProducer ( )

override

Destructor.

Definition at line 56 of file DTRecSegment2DExtendedProducer.cc.

References gather_cfg::cout, debug, and theAlgo.

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

Member Function Documentation

◆ produce()

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

override

The method which produces the 2D-segments.

Definition at line 63 of file DTRecSegment2DExtendedProducer.cc.

References edm::OwnVector< T, P >::begin(), filterCSVwithJSON::copy, gather_cfg::cout, debug, dtGeomToken_, edm::OwnVector< T, P >::empty(), edm::OwnVector< T, P >::end(), DTRangeMapAccessor::layersBySuperLayer(), eostools::move(), isotrackApplyRegressor::range, recClusToken_, recHits1DToken_, DTCombinatorialExtendedPatternReco::reconstruct(), DTCombinatorialExtendedPatternReco::setClusters(), DTCombinatorialExtendedPatternReco::setES(), singleTopDQM_cfi::setup, edm::OwnVector< T, P >::size(), DTGeometry::superLayer(), DTLayerId::superlayerId(), and theAlgo.

                                                                                         {
   if (debug)
     cout << "[DTRecSegment2DExtendedProducer] produce called" << endl;
   // Get the DT Geometry
   const DTGeometry& dtGeom = setup.getData(dtGeomToken_);
 
   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));
 }