#include <GEMRecHitProducer.h>

Inheritance diagram for GEMRecHitProducer:

Public Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override

	GEMRecHitProducer (const edm::ParameterSet &config)
	Constructor. More...

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

	~GEMRecHitProducer () override
	Destructor. More...

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () 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
GEMRecHitBaseAlgo *	theAlgo

edm::EDGetTokenT< GEMDigiCollection >	theGEMDigiToken

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
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::stream::EDProducerBase
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

Module for GEMRecHit production.

Author: M. Maggim – INFN Bari

Definition at line 40 of file GEMRecHitProducer.h.

Constructor & Destructor Documentation

GEMRecHitProducer::GEMRecHitProducer ( const edm::ParameterSet & config )

Constructor.

Definition at line 37 of file GEMRecHitProducer.cc.

References reco::get(), and edm::ParameterSet::getParameter().

                                                               {
 
   // Set verbose output
 
   produces<GEMRecHitCollection>();
 
   theGEMDigiToken = consumes<GEMDigiCollection>(config.getParameter<edm::InputTag>("gemDigiLabel"));  
 
   // Get the concrete reconstruction algo from the factory
 
   string theAlgoName = config.getParameter<string>("recAlgo");
   theAlgo = GEMRecHitAlgoFactory::get()->create(theAlgoName,
                         config.getParameter<ParameterSet>("recAlgoConfig"));
 
   // Get masked- and dead-strip information
 
   /* GEMMaskedStripsObj = new GEMMaskedStrips();
 
   GEMDeadStripsObj = new GEMDeadStrips();
 
   maskSource = config.getParameter<std::string>("maskSource");
 
   if (maskSource == "File") {
     edm::FileInPath fp = config.getParameter<edm::FileInPath>("maskvecfile");
     std::ifstream inputFile(fp.fullPath().c_str(), std::ios::in);
     if ( !inputFile ) {
       std::cerr << "Masked Strips File cannot not be opened" << std::endl;
       exit(1);
     }
     while ( inputFile.good() ) {
       GEMMaskedStrips::MaskItem Item;
       inputFile >> Item.rawId >> Item.strip;
       if ( inputFile.good() ) MaskVec.push_back(Item);
     }
     inputFile.close();
   }
 
   deadSource = config.getParameter<std::string>("deadSource");
 
   if (deadSource == "File") {
     edm::FileInPath fp = config.getParameter<edm::FileInPath>("deadvecfile");
     std::ifstream inputFile(fp.fullPath().c_str(), std::ios::in);
     if ( !inputFile ) {
       std::cerr << "Dead Strips File cannot not be opened" << std::endl;
       exit(1);
     }
     while ( inputFile.good() ) {
       GEMDeadStrips::DeadItem Item;
       inputFile >> Item.rawId >> Item.strip;
       if ( inputFile.good() ) DeadVec.push_back(Item);
     }
     inputFile.close();
   }
   */
 }

GEMRecHitProducer::~GEMRecHitProducer ( )

override

Destructor.

Definition at line 94 of file GEMRecHitProducer.cc.

                                      {
 
   delete theAlgo;
   // delete GEMMaskedStripsObj;
   // delete GEMDeadStripsObj;
 
 }

Member Function Documentation

void GEMRecHitProducer::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	setup
	)

override

Definition at line 104 of file GEMRecHitProducer.cc.

                                                                            {
 
   // Getting the masked-strip information
   /*
   if ( maskSource == "EventSetup" ) {
     edm::ESHandle<GEMMaskedStrips> readoutMaskedStrips;
     setup.get<GEMMaskedStripsRcd>().get(readoutMaskedStrips);
     const GEMMaskedStrips* tmp_obj = readoutMaskedStrips.product();
     GEMMaskedStripsObj->MaskVec = tmp_obj->MaskVec;
     delete tmp_obj;
   }
   else if ( maskSource == "File" ) {
     std::vector<GEMMaskedStrips::MaskItem>::iterator posVec;
     for ( posVec = MaskVec.begin(); posVec != MaskVec.end(); ++posVec ) {
       GEMMaskedStrips::MaskItem Item; 
       Item.rawId = (*posVec).rawId;
       Item.strip = (*posVec).strip;
       GEMMaskedStripsObj->MaskVec.push_back(Item);
     }
   }
   */
   // Getting the dead-strip information
   /*
   if ( deadSource == "EventSetup" ) {
     edm::ESHandle<GEMDeadStrips> readoutDeadStrips;
     setup.get<GEMDeadStripsRcd>().get(readoutDeadStrips);
     const GEMDeadStrips* tmp_obj = readoutDeadStrips.product();
     GEMDeadStripsObj->DeadVec = tmp_obj->DeadVec;
     delete tmp_obj;
   }
   else if ( deadSource == "File" ) {
     std::vector<GEMDeadStrips::DeadItem>::iterator posVec;
     for ( posVec = DeadVec.begin(); posVec != DeadVec.end(); ++posVec ) {
       GEMDeadStrips::DeadItem Item;
       Item.rawId = (*posVec).rawId;
       Item.strip = (*posVec).strip;
       GEMDeadStripsObj->DeadVec.push_back(Item);
     }
   }
   */
 }

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

override

The method which produces the rechits.

Definition at line 148 of file GEMRecHitProducer.cc.

References edm::OwnVector< T, P >::begin(), edm::OwnVector< T, P >::empty(), edm::OwnVector< T, P >::end(), GEMGeometry::etaPartition(), edm::EventSetup::get(), RecoTauDiscriminantConfiguration::mask, and eostools::move().

                                                                      {
 
   // Get the GEM Geometry
 
   ESHandle<GEMGeometry> gemGeom;
   setup.get<MuonGeometryRecord>().get(gemGeom);
 
   // Get the digis from the event
 
   Handle<GEMDigiCollection> digis; 
   event.getByToken(theGEMDigiToken,digis);
 
   // Pass the EventSetup to the algo
 
   theAlgo->setES(setup);
 
   // Create the pointer to the collection which will store the rechits
 
   auto recHitCollection = std::make_unique<GEMRecHitCollection>();
 
   // Iterate through all digi collections ordered by LayerId   
 
   GEMDigiCollection::DigiRangeIterator gemdgIt;
   for (gemdgIt = digis->begin(); gemdgIt != digis->end();
        ++gemdgIt){
        
     // The layerId
     const GEMDetId& gemId = (*gemdgIt).first;
 
     // Get the GeomDet from the setup
     const GEMEtaPartition* roll = gemGeom->etaPartition(gemId);
 
     // Get the iterators over the digis associated with this LayerId
     const GEMDigiCollection::Range& range = (*gemdgIt).second;
 
 
     // Getting the roll mask, that includes dead strips, for the given GEMDet
     EtaPartitionMask mask;
     /*
     int rawId = gemId.rawId();
     int Size = GEMMaskedStripsObj->MaskVec.size();
     for (int i = 0; i < Size; i++ ) {
       if ( GEMMaskedStripsObj->MaskVec[i].rawId == rawId ) {
     int bit = GEMMaskedStripsObj->MaskVec[i].strip;
     mask.set(bit-1);
       }
     }
 
     Size = GEMDeadStripsObj->DeadVec.size();
     for (int i = 0; i < Size; i++ ) {
       if ( GEMDeadStripsObj->DeadVec[i].rawId == rawId ) {
     int bit = GEMDeadStripsObj->DeadVec[i].strip;
     mask.set(bit-1);
       }
     }
     */
     // Call the reconstruction algorithm    
 
     OwnVector<GEMRecHit> recHits =
       theAlgo->reconstruct(*roll, gemId, range, mask);
     
     if(!recHits.empty()) //FIXME: is it really needed?
       recHitCollection->put(gemId, recHits.begin(), recHits.end());
   }
 
   event.put(std::move(recHitCollection));
 
 }

Member Data Documentation

GEMRecHitBaseAlgo* GEMRecHitProducer::theAlgo

private

Definition at line 61 of file GEMRecHitProducer.h.

edm::EDGetTokenT<GEMDigiCollection> GEMRecHitProducer::theGEMDigiToken

private

Definition at line 58 of file GEMRecHitProducer.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation