GEMDigiProducer Class Reference

Inheritance diagram for GEMDigiProducer:

Public Types
typedef edm::DetSetVector< GEMDigiSimLink >	GEMDigiSimLinks
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

Public Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override
	GEMDigiProducer (const edm::ParameterSet &ps)
void	produce (edm::Event &, const edm::EventSetup &) override
	~GEMDigiProducer () 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 Attributes
edm::EDGetTokenT< CrossingFrame< PSimHit > >	cf_token
std::unique_ptr< GEMDigiModule >	gemDigiModule_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geom_token_
const GEMGeometry *	geometry_

Detailed Description

Definition at line 37 of file GEMDigiProducer.cc.

Member Typedef Documentation

GEMDigiSimLinks

typedef edm::DetSetVector<GEMDigiSimLink> GEMDigiProducer::GEMDigiSimLinks

Definition at line 39 of file GEMDigiProducer.cc.

Constructor & Destructor Documentation

GEMDigiProducer()

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

explicit

Definition at line 61 of file GEMDigiProducer.cc.

References cf_token, Exception, geom_token_, edm::ParameterSet::getParameter(), ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                          : gemDigiModule_(std::make_unique<GEMDigiModule>(ps)) {
  produces<GEMDigiCollection>();
  produces<GEMDigiSimLinks>("GEM");

  edm::Service<edm::RandomNumberGenerator> rng;
  if (!rng.isAvailable()) {
    throw cms::Exception("Configuration")
        << "GEMDigiProducer::GEMDigiProducer() - RandomNumberGeneratorService is not present in configuration file.\n"
        << "Add the service in the configuration file or remove the modules that require it.";
  }

  std::string mix_(ps.getParameter<std::string>("mixLabel"));
  std::string collection_(ps.getParameter<std::string>("inputCollection"));

  cf_token = consumes<CrossingFrame<PSimHit> >(edm::InputTag(mix_, collection_));
  geom_token_ = esConsumes<GEMGeometry, MuonGeometryRecord, edm::Transition::BeginRun>();
}

~GEMDigiProducer()

GEMDigiProducer::~GEMDigiProducer ( )

overridedefault

Member Function Documentation

beginRun()

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

override

Definition at line 138 of file GEMDigiProducer.cc.

References options_cfi::eventSetup, gemDigiModule_, geom_token_, and geometry_.

                                                                             {
  edm::ESHandle<GEMGeometry> hGeom = eventSetup.getHandle(geom_token_);
  gemDigiModule_->setGeometry(&*hGeom);
  geometry_ = &*hGeom;
}

fillDescriptions()

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

static

Definition at line 81 of file GEMDigiProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                 {
  edm::ParameterSetDescription desc;
  desc.add<std::string>("inputCollection", "g4SimHitsMuonGEMHits");
  desc.add<std::string>("mixLabel", "mix");

  desc.add<double>("signalPropagationSpeed", 0.66);
  desc.add<double>("timeResolution", 5.);
  desc.add<double>("timeJitter", 1.0);
  desc.add<double>("averageShapingTime", 50.0);
  desc.add<double>("averageEfficiency", 0.98);
  desc.add<double>("averageNoiseRate", 0.001);
  // intrinsic noise rate (Hz/cm^2)

  // in terms of 25 ns
  desc.add<int>("minBunch", -5);
  desc.add<int>("maxBunch", 3);

  desc.add<bool>("fixedRollRadius", true);
  // Uses fixed radius in the center of the roll
  desc.add<bool>("digitizeOnlyMuons", false);
  desc.add<bool>("simulateBkgNoise", false);
  // false == No background simulation
  desc.add<bool>("simulateNoiseCLS", true);
  desc.add<bool>("simulateElectronBkg", true);
  // flase == simulate only neutral bkg
  desc.add<bool>("simulateIntrinsicNoise", false);
  desc.add<bool>("bx0filter", false);

  desc.add<double>("instLumi", 7.5);
  // in units of 1E34 cm^-2 s^-1. Internally the background is parmetrized from FLUKA+GEANT result at 5E+34 (PU 140). We are adding a 1.5 factor for PU 200
  desc.add<double>("rateFact", 1.0);
  // We are adding also a safety factor of 2 to tak into account the new beam pipe effect (not yet known). Hits can be thrown away later at re-digi step. Parameters are kept in sync with the ones used in the GEM digitizer
  desc.add<double>("bxWidth", 25E-9);
  desc.add<double>("referenceInstLumi", 5.);
  // referecne inst. luminosity 5E+34 cm^-2s^-1
  desc.add<double>("resolutionX", 0.03);

  // The follwing parameters are needed to model the background contribution
  // The parameters have been obtained after the fit of th perdicted by FLUKA
  // By default the backgroundmodeling with these parameters should be disabled with
  // the 9_2_X release setting simulateBkgNoise = false
  desc.add<double>("GE11ModNeuBkgParam0", 5710.23);
  desc.add<double>("GE11ModNeuBkgParam1", -43.3928);
  desc.add<double>("GE11ModNeuBkgParam2", 0.0863681);
  desc.add<double>("GE21ModNeuBkgParam0", 1440.44);
  desc.add<double>("GE21ModNeuBkgParam1", -7.48607);
  desc.add<double>("GE21ModNeuBkgParam2", 0.0103078);
  desc.add<double>("GE11ElecBkgParam0", 406.249);
  desc.add<double>("GE11ElecBkgParam1", -2.90939);
  desc.add<double>("GE11ElecBkgParam2", 0.00548191);
  desc.add<double>("GE21ElecBkgParam0", 97.0505);
  desc.add<double>("GE21ElecBkgParam1", -43.3928);
  desc.add<double>("GE21ElecBkgParam2", 00.000550599);

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

produce()

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

override

Definition at line 144 of file GEMDigiProducer.cc.

References cf_token, hcalRecHitTable_cff::detId, MillePedeFileConverter_cfg::e, GEMGeometry::etaPartitions(), gemDigiModule_, geometry_, edm::RandomNumberGenerator::getEngine(), hfClusterShapes_cfi::hits, LogDebug, eostools::move(), edm::Handle< T >::product(), nano_mu_digi_cff::rawId, nano_mu_digi_cff::roll, and FastTrackerRecHitCombiner_cfi::simHits.

                                                                          {
  edm::Service<edm::RandomNumberGenerator> rng;
  CLHEP::HepRandomEngine* engine = &rng->getEngine(e.streamID());

  edm::Handle<CrossingFrame<PSimHit> > cf;
  e.getByToken(cf_token, cf);

  MixCollection<PSimHit> hits{cf.product()};

  // Create empty output
  auto digis = std::make_unique<GEMDigiCollection>();
  auto gemDigiSimLinks = std::make_unique<GEMDigiSimLinks>();

  // arrange the hits by eta partition
  std::map<uint32_t, edm::PSimHitContainer> hitMap;
  for (const auto& hit : hits) {
    hitMap[GEMDetId(hit.detUnitId()).rawId()].emplace_back(hit);
  }

  // simulate signal and noise for each eta partition
  const auto& etaPartitions(geometry_->etaPartitions());

  for (const auto& roll : etaPartitions) {
    const GEMDetId detId(roll->id());
    const uint32_t rawId(detId.rawId());
    const auto& simHits(hitMap[rawId]);

    LogDebug("GEMDigiProducer") << "GEMDigiProducer: found " << simHits.size() << " hit(s) in eta partition" << rawId;

    gemDigiModule_->simulate(roll, simHits, engine);
    gemDigiModule_->fillDigis(rawId, *digis);
    (*gemDigiSimLinks).insert(gemDigiModule_->gemDigiSimLinks());
  }

  // store them in the event
  e.put(std::move(digis));
  e.put(std::move(gemDigiSimLinks), "GEM");
}

Member Data Documentation

cf_token

edm::EDGetTokenT<CrossingFrame<PSimHit> > GEMDigiProducer::cf_token

private

Definition at line 53 of file GEMDigiProducer.cc.

Referenced by GEMDigiProducer(), and produce().

gemDigiModule_

std::unique_ptr<GEMDigiModule> GEMDigiProducer::gemDigiModule_

private

Definition at line 58 of file GEMDigiProducer.cc.

Referenced by beginRun(), and produce().

geom_token_

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

private

Definition at line 54 of file GEMDigiProducer.cc.

Referenced by beginRun(), and GEMDigiProducer().

geometry_

const GEMGeometry* GEMDigiProducer::geometry_

private

Definition at line 56 of file GEMDigiProducer.cc.

Referenced by beginRun(), and produce().