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

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

Definition at line 39 of file GEMDigiProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 61 of file GEMDigiProducer.cc.

References cf_token, Exception, geom_token_, edm::ParameterSet::getParameter(), HLT_FULL_cff::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 ( )

overridedefault

Member Function Documentation

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

override

Definition at line 138 of file GEMDigiProducer.cc.

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

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

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

static

Definition at line 81 of file GEMDigiProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::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);
 }

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

override

Definition at line 144 of file GEMDigiProducer.cc.

References cf_token, GEMGeometry::etaPartitions(), gemDigiModule_, geometry_, edm::Event::getByToken(), edm::RandomNumberGenerator::getEngine(), LogDebug, eostools::move(), edm::Handle< T >::product(), edm::Event::put(), trackerHits::simHits, and edm::Event::streamID().

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