DTRecSegment2DProducer Class Reference

#include <DTRecSegment2DProducer.h>

Inheritance diagram for DTRecSegment2DProducer:

Public Member Functions
	DTRecSegment2DProducer (const edm::ParameterSet &)
	Constructor. More...
void	produce (edm::Event &event, const edm::EventSetup &setup) override
	The method which produces the 2D-segments. More...
	~DTRecSegment2DProducer () override
	Destructor. More...
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final

Private Attributes
bool	debug
edm::EDGetTokenT< DTRecHitCollection >	recHits1DToken_
DTRecSegment2DBaseAlgo *	theAlgo

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

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 31 of file DTRecSegment2DProducer.h.

Constructor & Destructor Documentation

DTRecSegment2DProducer::DTRecSegment2DProducer

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 36 of file DTRecSegment2DProducer.cc.

References gather_cfg::cout, debug, reco::get(), 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"));

  if(debug)
    cout << "[DTRecSegment2DProducer] Constructor called" << endl;

  produces<DTRecSegment2DCollection>();

  // Get the concrete reconstruction algo from the factory
  string theAlgoName = pset.getParameter<string>("Reco2DAlgoName");
  if(debug) cout << "the Reco2D AlgoName is " << theAlgoName << endl;
  theAlgo = DTRecSegment2DAlgoFactory::get()->create(theAlgoName,
                                                     pset.getParameter<ParameterSet>("Reco2DAlgoConfig"));
}

DTRecSegment2DProducer::~DTRecSegment2DProducer ( )

override

Destructor.

Definition at line 56 of file DTRecSegment2DProducer.cc.

References gather_cfg::cout, and debug.

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

Member Function Documentation

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

override

The method which produces the 2D-segments.

Definition at line 63 of file DTRecSegment2DProducer.cc.

References edm::OwnVector< T, P >::begin(), 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 << "[DTRecSegment2DProducer] 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);

  // 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;

    if (!segs.empty() )
      segments->put(SLId, segs.begin(),segs.end());
  }
  event.put(std::move(segments));
}

Member Data Documentation

bool DTRecSegment2DProducer::debug

private

Definition at line 50 of file DTRecSegment2DProducer.h.

Referenced by runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), and runTauIdMVA.TauIDEmbedder::runTauID().

edm::EDGetTokenT<DTRecHitCollection> DTRecSegment2DProducer::recHits1DToken_

private

Definition at line 56 of file DTRecSegment2DProducer.h.

DTRecSegment2DBaseAlgo* DTRecSegment2DProducer::theAlgo

private

Definition at line 53 of file DTRecSegment2DProducer.h.