GEMDigiProducer Class Reference

Inheritance diagram for GEMDigiProducer:

Public Types
typedef edm::DetSetVector< GEMDigiSimLink >	GEMDigiSimLinks
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 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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

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 38 of file GEMDigiProducer.cc.

Member Typedef Documentation

GEMDigiSimLinks

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

Definition at line 40 of file GEMDigiProducer.cc.

Constructor & Destructor Documentation

GEMDigiProducer()

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

explicit

Definition at line 62 of file GEMDigiProducer.cc.

                                                          : 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>();
}

References cf_token, Exception, geom_token_, edm::ParameterSet::getParameter(), HLT_FULL_cff::InputTag, edm::Service< T >::isAvailable(), and AlCaHLTBitMon_QueryRunRegistry::string.

~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.

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

References gemDigiModule_, geom_token_, geometry_, and edm::EventSetup::getHandle().

fillDescriptions()

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

static

Definition at line 82 of file GEMDigiProducer.cc.

                                                                                 {
  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<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);
}

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

produce()

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

override

Definition at line 144 of file GEMDigiProducer.cc.

                                                                          {
  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");
}

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

Member Data Documentation

cf_token

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

private

Definition at line 54 of file GEMDigiProducer.cc.

Referenced by GEMDigiProducer(), and produce().

gemDigiModule_

std::unique_ptr<GEMDigiModule> GEMDigiProducer::gemDigiModule_

private

Definition at line 59 of file GEMDigiProducer.cc.

Referenced by beginRun(), and produce().

geom_token_

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

private

Definition at line 55 of file GEMDigiProducer.cc.

Referenced by beginRun(), and GEMDigiProducer().

geometry_

const GEMGeometry* GEMDigiProducer::geometry_

private

Definition at line 57 of file GEMDigiProducer.cc.

Referenced by beginRun(), and produce().