GEMPadDigiProducer Class Reference

Inheritance diagram for GEMPadDigiProducer:

Public Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override
	GEMPadDigiProducer (const edm::ParameterSet &ps)
void	produce (edm::Event &, const edm::EventSetup &) override
	~GEMPadDigiProducer () override
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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	buildPads (const GEMDigiCollection &digis, GEMPadDigiCollection &out_pads) const
void	buildPads16GE21 (const GEMDigiCollection &digis, GEMPadDigiCollection &out_pads) const
void	checkGeometry () const
void	checkValid (const GEMPadDigi &pad, const GEMDetId &id) const

Private Attributes
edm::EDGetTokenT< GEMDigiCollection >	digi_token_
	Name of input digi Collection. More...
edm::InputTag	digis_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geom_token_
const GEMGeometry *	geometry_
bool	use16GE21_

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

Definition at line 33 of file GEMPadDigiProducer.cc.

Constructor & Destructor Documentation

GEMPadDigiProducer()

GEMPadDigiProducer::GEMPadDigiProducer ( const edm::ParameterSet &  ps )

explicit

Definition at line 60 of file GEMPadDigiProducer.cc.

References digi_token_, digis_, geom_token_, and edm::ParameterSet::getParameter().

                                                                : use16GE21_(false), geometry_(nullptr) {
  digis_ = ps.getParameter<edm::InputTag>("InputCollection");

  digi_token_ = consumes<GEMDigiCollection>(digis_);
  geom_token_ = esConsumes<GEMGeometry, MuonGeometryRecord, edm::Transition::BeginRun>();

  produces<GEMPadDigiCollection>();
  consumes<GEMDigiCollection>(digis_);
}

~GEMPadDigiProducer()

GEMPadDigiProducer::~GEMPadDigiProducer ( )

override

Definition at line 70 of file GEMPadDigiProducer.cc.

{}

Member Function Documentation

beginRun()

void GEMPadDigiProducer::beginRun ( const edm::Run &  run, const edm::EventSetup &  eventSetup  )

override

Definition at line 79 of file GEMPadDigiProducer.cc.

References GEMGeometry::chamber(), checkGeometry(), options_cfi::eventSetup, geom_token_, geometry_, GEMGeometry::hasGE21(), GEMChamber::nEtaPartitions(), and use16GE21_.

                                                                                    {
  edm::ESHandle<GEMGeometry> hGeom = eventSetup.getHandle(geom_token_);
  geometry_ = &*hGeom;
  // check the number of parititions
  if (geometry_->hasGE21()) {
    use16GE21_ = (geometry_->chamber(GEMDetId(1, 1, 2, 2, 16, 0))->nEtaPartitions() ==
                  GEMPadDigi::NumberPartitions::GE21SplitStrip);
  }

  checkGeometry();
}

buildPads()

void GEMPadDigiProducer::buildPads ( const GEMDigiCollection &  digis, GEMPadDigiCollection &  out_pads  ) const

private

Definition at line 107 of file GEMPadDigiProducer.cc.

References checkValid(), ztail::d, GEMGeometry::etaPartitions(), GEMPadDigi::GE11InValid, GEMPadDigi::GE21InValid, geometry_, l1t::mtf7::ME0, GEMPadDigi::ME0InValid, nPart(), or, AlCaHLTBitMon_ParallelJobs::p, and use16GE21_.

Referenced by produce().

                                                                                                           {
  for (const auto& p : geometry_->etaPartitions()) {
    // when using the GE2/1 geometry with 16 eta partitions
    // ->ignore GE2/1
    if (use16GE21_ and p->isGE21())
      continue;

    // set of <pad, bx, part> pairs, sorted first by pad then by bx
    std::set<std::tuple<int, int, unsigned> > proto_pads;

    // walk over digis in this partition,
    // and stuff them into a set of unique pads (equivalent of OR operation)
    auto digis = det_digis.get(p->id());
    for (auto d = digis.first; d != digis.second; ++d) {
      unsigned pad_num = static_cast<unsigned>(p->padOfStrip(d->strip()));

      auto nPart = GEMPadDigi::NumberPartitions::GE11;
      if (p->isME0()) {
        nPart = GEMPadDigi::NumberPartitions::ME0;
      } else if (p->isGE21()) {
        nPart = GEMPadDigi::NumberPartitions::GE21;
      }
      // check that the input digi is valid
      if ((p->isGE11() and pad_num == GEMPadDigi::GE11InValid) or
          (p->isGE21() and pad_num == GEMPadDigi::GE21InValid) or (p->isME0() and pad_num == GEMPadDigi::ME0InValid)) {
        edm::LogWarning("GEMPadDigiProducer") << "Invalid " << pad_num << " from  " << *d << " in " << p->id();
      }
      proto_pads.emplace(pad_num, d->bx(), nPart);
    }

    // fill the output collections
    for (const auto& d : proto_pads) {
      GEMPadDigi pad_digi(std::get<0>(d), std::get<1>(d), p->subsystem(), std::get<2>(d));
      checkValid(pad_digi, p->id());
      out_pads.insertDigi(p->id(), pad_digi);
    }
  }
}

buildPads16GE21()

void GEMPadDigiProducer::buildPads16GE21 ( const GEMDigiCollection &  digis, GEMPadDigiCollection &  out_pads  ) const

private

Definition at line 146 of file GEMPadDigiProducer.cc.

References checkValid(), ztail::d, GEMGeometry::etaPartitions(), geometry_, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by produce().

                                                                                                                 {
  for (const auto& p : geometry_->etaPartitions()) {
    // when using the GE2/1 geometry with 16 eta partitions
    // ->ignore GE1/1
    if (!p->isGE21())
      continue;

    // ignore eta partition with even numbers
    // these are included in the odd numbered pads
    if (p->id().roll() % 2 == 0)
      continue;

    // set of <pad, bx> pairs, sorted first by pad then by bx
    std::set<std::pair<int, int> > proto_pads;

    // walk over digis in the first partition,
    // and stuff them into a set of unique pads (equivalent of OR operation)
    auto digis = det_digis.get(p->id());

    // proto pads for the odd partitions
    for (auto d = digis.first; d != digis.second; ++d) {
      proto_pads.emplace(d->strip(), d->bx());
    }

    GEMDetId gemId2(
        p->id().region(), p->id().ring(), p->id().station(), p->id().layer(), p->id().chamber(), p->id().roll() + 1);
    auto digis2 = det_digis.get(gemId2);

    // proto pads for the even partitions
    for (auto d = digis2.first; d != digis2.second; ++d) {
      proto_pads.emplace(d->strip(), d->bx());
    }

    // fill the output collections
    for (const auto& d : proto_pads) {
      GEMPadDigi pad_digi(d.first, d.second, p->subsystem(), GEMPadDigi::NumberPartitions::GE21SplitStrip);
      checkValid(pad_digi, p->id());
      out_pads.insertDigi(p->id(), pad_digi);
    }
  }
}

checkGeometry()

void GEMPadDigiProducer::checkGeometry ( ) const

private

Definition at line 196 of file GEMPadDigiProducer.cc.

References GEMGeometry::chamber(), geometry_, GEMGeometry::hasGE11(), GEMGeometry::hasGE21(), GEMGeometry::hasME0(), l1t::mtf7::ME0, GEMChamber::nEtaPartitions(), and use16GE21_.

Referenced by beginRun().

                                             {
  // check that ME0 has 8-eta partitions
  if (geometry_->hasME0()) {
    if (geometry_->chamber(GEMDetId(1, 1, 0, 1, 1, 0))->nEtaPartitions() != GEMPadDigi::NumberPartitions::ME0) {
      edm::LogError("GEMPadDigiProducer") << "ME0 geometry appears corrupted";
    }
  }

  // check that GE1/1 has 8-eta partitions
  if (geometry_->hasGE11()) {
    if (geometry_->chamber(GEMDetId(1, 1, 1, 1, 1, 0))->nEtaPartitions() != GEMPadDigi::NumberPartitions::GE11) {
      edm::LogError("GEMPadDigiProducer") << "GE1/1 geometry appears corrupted";
    }
  }

  if (geometry_->hasGE21()) {
    if (!use16GE21_) {
      // check that GE2/1 has 8-eta partitions
      if (geometry_->chamber(GEMDetId(1, 1, 2, 2, 16, 0))->nEtaPartitions() != GEMPadDigi::NumberPartitions::GE21) {
        edm::LogError("GEMPadDigiProducer") << "GE2/1 geometry (8 partition) appears corrupted";
      }
    } else {
      // check that GE2/1 has 16-eta partitions
      if (geometry_->chamber(GEMDetId(1, 1, 2, 2, 16, 0))->nEtaPartitions() !=
          GEMPadDigi::NumberPartitions::GE21SplitStrip) {
        edm::LogError("GEMPadDigiProducer") << "GE2/1 geometry (16 partition) appears corrupted";
      }
    }
  }
}

checkValid()

void GEMPadDigiProducer::checkValid ( const GEMPadDigi &  pad, const GEMDetId &  id  ) const

private

Definition at line 188 of file GEMPadDigiProducer.cc.

References EcalPhiSymFlatTableProducers_cfi::id, and GEMPadDigi::isValid().

Referenced by buildPads(), and buildPads16GE21().

                                                                                   {
  // check if the pad is valid
  // in principle, invalid pads can appear in the CMS raw data
  if (!pad.isValid()) {
    edm::LogWarning("GEMPadDigiProducer") << "Invalid " << pad << " in " << id;
  }
}

fillDescriptions()

void GEMPadDigiProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 72 of file GEMPadDigiProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                    {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("InputCollection", edm::InputTag("simMuonGEMDigis"));

  descriptions.add("simMuonGEMPadDigisDef", desc);
}

produce()

void GEMPadDigiProducer::produce ( edm::Event &  e, const edm::EventSetup &  eventSetup  )

override

Definition at line 91 of file GEMPadDigiProducer.cc.

References buildPads(), buildPads16GE21(), digi_token_, MillePedeFileConverter_cfg::e, eostools::move(), edm::Handle< T >::product(), and use16GE21_.

                                                                             {
  edm::Handle<GEMDigiCollection> hdigis;
  e.getByToken(digi_token_, hdigis);

  // Create empty output
  std::unique_ptr<GEMPadDigiCollection> pPads(new GEMPadDigiCollection());

  // build the pads
  buildPads(*(hdigis.product()), *pPads);
  if (use16GE21_)
    buildPads16GE21(*(hdigis.product()), *pPads);

  // store them in the event
  e.put(std::move(pPads));
}

Member Data Documentation

digi_token_

edm::EDGetTokenT<GEMDigiCollection> GEMPadDigiProducer::digi_token_

private

Name of input digi Collection.

Definition at line 52 of file GEMPadDigiProducer.cc.

Referenced by GEMPadDigiProducer(), and produce().

digis_

edm::InputTag GEMPadDigiProducer::digis_

private

Definition at line 54 of file GEMPadDigiProducer.cc.

Referenced by GEMPadDigiProducer().

geom_token_

edm::ESGetToken<GEMGeometry, MuonGeometryRecord> GEMPadDigiProducer::geom_token_

private

Definition at line 53 of file GEMPadDigiProducer.cc.

Referenced by beginRun(), and GEMPadDigiProducer().

geometry_

const GEMGeometry* GEMPadDigiProducer::geometry_

private

Definition at line 57 of file GEMPadDigiProducer.cc.

Referenced by beginRun(), buildPads(), buildPads16GE21(), and checkGeometry().

use16GE21_

bool GEMPadDigiProducer::use16GE21_

private

Definition at line 55 of file GEMPadDigiProducer.cc.

Referenced by beginRun(), buildPads(), checkGeometry(), and produce().