StripCompactDigiSimLinksProducer Class Reference

Inheritance diagram for StripCompactDigiSimLinksProducer:

Public Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override
	StripCompactDigiSimLinksProducer (const edm::ParameterSet &iConfig)
	~StripCompactDigiSimLinksProducer () override
Public Member Functions inherited from edm::one::EDProducer<>
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () 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)
	~ProducerBase () noexcept(false) override
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
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
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
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
uint32_t	maxHoleSize_
edm::InputTag	src_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	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::one::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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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 23 of file StripCompactDigiSimLinksProducer.cc.

Constructor & Destructor Documentation

StripCompactDigiSimLinksProducer::StripCompactDigiSimLinksProducer

(

const edm::ParameterSet &

iConfig

)

Definition at line 35 of file StripCompactDigiSimLinksProducer.cc.

    : src_(iConfig.getParameter<edm::InputTag>("src")), maxHoleSize_(iConfig.getParameter<uint32_t>("maxHoleSize")) {
  produces<StripCompactDigiSimLinks>();
}

StripCompactDigiSimLinksProducer::~StripCompactDigiSimLinksProducer ( )

override

Definition at line 40 of file StripCompactDigiSimLinksProducer.cc.

{}

Member Function Documentation

void StripCompactDigiSimLinksProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &    )

override

Definition at line 42 of file StripCompactDigiSimLinksProducer.cc.

References edm::DetSet< T >::begin(), MessageLogger_cfi::cerr, DEBUG, DEFINE_FWK_MODULE, spr::find(), edm::Event::getByLabel(), StripCompactDigiSimLinks::Filler::insert(), crabWrapper::key, maxHoleSize_, eostools::move(), convertSQLitetoXML_cfg::output, edm::Event::put(), findQualityFiles::size, TrackRefitter_38T_cff::src, src_, and edm::swap().

                                                                                      {
  using namespace edm;
  Handle<DetSetVector<StripDigiSimLink>> src;
  iEvent.getByLabel(src_, src);

  StripCompactDigiSimLinks::Filler output;

  int previousStrip;      // previous strip with at least one link (might not be the
                          // strip of the previous link there are overlapping
                          // clusters)
  int previousLinkStrip;  // strip of the previous link (can be the same as the
                          // one of this link if there are overlapping clusters)
  std::vector<StripCompactDigiSimLinks::key_type> thisStripSignals;      // particles on this strip
  std::vector<StripCompactDigiSimLinks::key_type> previousStripSignals;  // particles on the previous strip

  for (auto const &det : *src) {
    DEBUG(std::cerr << "\n\nProcessing detset " << det.detId() << ", size = " << det.size() << std::endl;)
    previousStrip = -2;      // previous strip with at least one link (might not be the strip of
                             // the previous link there are overlapping clusters)
    previousLinkStrip = -2;  // strip of the previous link (can be the same as the one of this
                             // link if there are overlapping clusters)
    thisStripSignals.clear();
    previousStripSignals.clear();
    for (DetSet<StripDigiSimLink>::const_iterator it = det.begin(), ed = det.end(); it != ed; ++it) {
      DEBUG(std::cerr << "  processing digiSimLink on strip " << it->channel() << " left by particle "
                      << it->SimTrackId() << ", event " << it->eventId().rawId() << std::endl;)
      if (int(it->channel()) != previousLinkStrip) {
        previousStrip = previousLinkStrip;
        DEBUG(std::cerr << "   strip changed!" << std::endl;)
        swap(thisStripSignals, previousStripSignals);
        thisStripSignals.clear();
      }
      DEBUG(std::cerr << "    previous strip " << previousStrip << ", previous link strip " << previousLinkStrip
                      << std::endl;)
      // DEBUG(std::cerr << "    on previous strip: "; for(auto const& k :
      // previousStripSignals) { std::cerr << "(ev " << k.first.rawId() << ", id
      // " << k.second << ") "; } std::cerr << std::endl;) DEBUG(std::cerr << "
      // on this strip: ";     for(auto const& k :     thisStripSignals) {
      // std::cerr << "(ev " << k.first.rawId() << ", id " << k.second << ") ";
      // } std::cerr << std::endl;)
      StripCompactDigiSimLinks::key_type key(it->eventId(), it->SimTrackId());
      bool alreadyClusterized = false;
      if (int(it->channel()) == previousStrip + 1) {
        DEBUG(std::cerr << "  on next strip" << std::endl;)
        if (std::find(previousStripSignals.begin(), previousStripSignals.end(), key) != previousStripSignals.end()) {
          alreadyClusterized = true;
          DEBUG(std::cerr << "    and part of previous cluster" << std::endl;)
        }
      }
      if (!alreadyClusterized) {
        DEBUG(std::cerr << "   clusterize!" << std::endl;)
        unsigned int size = 1;
        int myLastStrip = it->channel();  // last known strip with signals from this particle
        for (DetSet<StripDigiSimLink>::const_iterator it2 = it + 1; it2 < ed; ++it2) {
          DEBUG(std::cerr << "     digiSimLink on strip " << it2->channel() << " left by particle " << it2->SimTrackId()
                          << ", event " << it2->eventId().rawId() << std::endl;)
          if ((it2->channel() - myLastStrip) > maxHoleSize_ + 1) {
            DEBUG(std::cerr << "        found hole of size " << (it2->channel() - myLastStrip) << ", stopping."
                            << std::endl;)
            break;
          }
          if ((it2->eventId() == key.first) && (it2->SimTrackId() == key.second)) {
            size++;
            DEBUG(std::cerr << "        extending cluster, now size = " << size << std::endl;)
            myLastStrip = it2->channel();
          }
        }
        output.insert(key, StripCompactDigiSimLinks::HitRecord(det.detId(), it->channel(), size));
      }
      if (int(it->channel()) != previousLinkStrip) {
        previousLinkStrip = it->channel();
      }
      thisStripSignals.push_back(key);
      DEBUG(std::cerr << "    ending state " << previousStrip << ", previous link strip " << previousLinkStrip
                      << std::endl;)
      // DEBUG(std::cerr << "    on previous strip: "; for(auto const& k :
      // previousStripSignals) { std::cerr << "(ev " << k.first.rawId() << ", id
      // " << k.second << ") "; } std::cerr << std::endl;) DEBUG(std::cerr << "
      // on this strip: ";     for(auto const& k :     thisStripSignals) {
      // std::cerr << "(ev " << k.first.rawId() << ", id " << k.second << ") ";
      // } std::cerr << std::endl;)
      DEBUG(std::cerr << std::endl;)
    }
  }

  std::unique_ptr<StripCompactDigiSimLinks> ptr(new StripCompactDigiSimLinks(output));
  iEvent.put(std::move(ptr));
}

Member Data Documentation

uint32_t StripCompactDigiSimLinksProducer::maxHoleSize_

private

Definition at line 32 of file StripCompactDigiSimLinksProducer.cc.

Referenced by produce().

edm::InputTag StripCompactDigiSimLinksProducer::src_

private

Definition at line 31 of file StripCompactDigiSimLinksProducer.cc.

Referenced by produce().