Inheritance diagram for RPDigiProducer:

Public Member Functions
	RPDigiProducer (const edm::ParameterSet &)

	~RPDigiProducer () override

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () 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)

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

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Private Types
typedef std::map< unsigned int, std::vector< PSimHit > >	simhit_map

typedef simhit_map::iterator	simhit_map_iterator

Private Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override

edm::DetSet< TotemRPDigi >	convertRPStripDetSet (const edm::DetSet< TotemRPDigi > &)

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::ParameterSet	conf_

DeadChannelsManager	deadChannelsManager

CLHEP::HepRandomEngine *	rndEngine_ = 0

std::vector< std::string >	RP_hit_containers_

bool	simulateDeadChannels

std::map< RPDetId, std::unique_ptr< RPDetDigitizer > >	theAlgoMap

edm::EDGetTokenT< CrossingFrame< PSimHit > >	tokenCrossingFrameTotemRP

int	verbosity_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	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

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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 53 of file RPDigiProducer.cc.

Member Typedef Documentation

typedef std::map<unsigned int, std::vector<PSimHit> > RPDigiProducer::simhit_map

private

Definition at line 68 of file RPDigiProducer.cc.

typedef simhit_map::iterator RPDigiProducer::simhit_map_iterator

private

Definition at line 69 of file RPDigiProducer.cc.

Constructor & Destructor Documentation

RPDigiProducer::RPDigiProducer ( const edm::ParameterSet & conf )

explicit

Definition at line 90 of file RPDigiProducer.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), HLT_2018_cff::InputTag, RP_hit_containers_, simulateDeadChannels, tokenCrossingFrameTotemRP, and verbosity_.

                                                           : conf_(conf) {
   //now do what ever other initialization is needed
   produces<edm::DetSetVector<TotemRPDigi>>();
 
   // register data to consume
   tokenCrossingFrameTotemRP = consumes<CrossingFrame<PSimHit>>(edm::InputTag("mix", "g4SimHitsTotemHitsRP", ""));
 
   RP_hit_containers_ = conf.getParameter<std::vector<std::string>>("ROUList");
   verbosity_ = conf.getParameter<int>("RPVerbosity");
 
   simulateDeadChannels = false;
   if (conf.exists(
           "simulateDeadChannels")) {  //check if "simulateDeadChannels" variable is defined in configuration file
     simulateDeadChannels = conf.getParameter<bool>("simulateDeadChannels");
   }
 }

RPDigiProducer::~RPDigiProducer ( )

override

Definition at line 107 of file RPDigiProducer.cc.

                                 {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

void RPDigiProducer::beginRun	(	const edm::Run &	beginrun,
		const edm::EventSetup &	es
	)

overrideprivate

Definition at line 202 of file RPDigiProducer.cc.

References deadChannelsManager, edm::EventSetup::get(), and simulateDeadChannels.

                                                                              {
   // get analysis mask to mask channels
   if (simulateDeadChannels) {
     edm::ESHandle<TotemAnalysisMask> analysisMask;
     es.get<TotemReadoutRcd>().get(analysisMask);
     deadChannelsManager = DeadChannelsManager(analysisMask);  //set analysisMask in deadChannelsManager
   }
 }

edm::DetSet< TotemRPDigi > RPDigiProducer::convertRPStripDetSet ( const edm::DetSet< TotemRPDigi > & rpstrip_detset )

private

Definition at line 211 of file RPDigiProducer.cc.

References edm::DetSet< T >::data, edm::DetSet< T >::detId(), edm::DetSet< T >::reserve(), and edm::DetSet< T >::size().

Referenced by produce().

                                                                                                         {
   edm::DetSet<TotemRPDigi> rpdigi_detset(rpstrip_detset.detId());
   rpdigi_detset.reserve(rpstrip_detset.size());
 
   for (std::vector<TotemRPDigi>::const_iterator stripIterator = rpstrip_detset.data.begin();
        stripIterator < rpstrip_detset.data.end();
        ++stripIterator) {
     rpdigi_detset.push_back(TotemRPDigi(stripIterator->stripNumber()));
   }
 
   return rpdigi_detset;
 }

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

static

Definition at line 224 of file RPDigiProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                 {
   //RPSiDetDigitizer
   //all distances in [mm]
   edm::ParameterSetDescription desc;
   desc.add<bool>("RPLandauFluctuations", true);
   desc.add<bool>("RPDisplacementOn", false);
   desc.add<int>("RPVerbosity", 0);
   desc.add<double>("RPVFATThreshold", 9000.0);
   desc.add<double>("RPTopEdgePosition", 1.5);
   desc.add<double>("RPActiveEdgeSmearing", 0.013);
   desc.add<double>("RPEquivalentNoiseCharge300um", 1000.0);
   desc.add<int>("RPVFATTriggerMode", 2);
   desc.add<std::vector<double>>("RPInterStripSmearing",
                                 {
                                     0.011,
                                 });
   desc.add<double>("RPSharingSigmas", 5.0);  //how many sigmas taken into account for the edges and inter strips
   desc.add<double>("RPGeVPerElectron", 3.61e-09);
   desc.add<double>("RPActiveEdgePosition", 0.034);  //from the physical edge
   desc.add<bool>("RPDeadStripSimulationOn", false);
   desc.add<std::vector<std::string>>("ROUList",
                                      {
                                          "TotemHitsRP",
                                      });
   desc.add<bool>("RPNoNoise", false);
   desc.add<bool>("RPDigiSimHitRelationsPresistence", false);  //save links betweend digi, clusters and OSCAR/Geant4 hits
   desc.add<std::string>("mixLabel", "mix");
   desc.add<int>("RPChargeDivisionsPerThickness", 5);
   desc.add<double>("RPDeltaProductionCut", 0.120425);  //[MeV]
   desc.add<double>("RPBottomEdgePosition", 1.5);
   desc.add<double>("RPBottomEdgeSmearing", 0.011);
   desc.add<double>("RPTopEdgeSmearing", 0.011);
   desc.add<std::string>("InputCollection", "g4SimHitsTotemHitsRP");
   desc.add<double>("RPInterStripCoupling",
                    1.0);  //fraction of charge going to the strip, the missing part is taken by its neighbours
   desc.add<double>("RPDeadStripProbability", 0.001);
   desc.add<int>("RPChargeDivisionsPerStrip", 15);
   descriptions.add("RPSiDetDigitizer", desc);
 }

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

overrideprivate

Definition at line 116 of file RPDigiProducer.cc.

References MixCollection< T >::MixItr::bunch(), cuy::col, conf_, convertRPStripDetSet(), KineDebug3::count(), edm::DetSet< T >::data, Exception, edm::Event::getByLabel(), edm::RandomNumberGenerator::getEngine(), MixCollection< T >::MixItr::getTrigger(), edm::Service< T >::isAvailable(), eostools::move(), edm::Handle< T >::product(), edm::Event::put(), rndEngine_, writedatasetfile::run, edm::Event::streamID(), theAlgoMap, and verbosity_.

                                                                           {
   using namespace edm;
 
   // initialize random engine
   if (!rndEngine_) {
     Service<RandomNumberGenerator> rng;
     if (!rng.isAvailable()) {
       throw cms::Exception("Configuration")
           << "This class requires the RandomNumberGeneratorService\n"
              "which is not present in the configuration file.  You must add the service\n"
              "in the configuration file or remove the modules that require it.";
     }
     rndEngine_ = &(rng->getEngine(iEvent.streamID()));
   }
 
   // Step A: Get Inputs
   edm::Handle<CrossingFrame<PSimHit>> cf;
   iEvent.getByLabel("mix", "g4SimHitsTotemHitsRP", cf);
 
   if (verbosity_) {
     edm::LogInfo("RPDigiProducer") << "\n\n=================== Starting SimHit access"
                                    << "  ==================="
                                    << "\n";
 
     MixCollection<PSimHit> col{cf.product(), std::pair(-0, 0)};
     MixCollection<PSimHit>::iterator cfi;
     int count = 0;
     for (cfi = col.begin(); cfi != col.end(); cfi++) {
       edm::LogInfo("RPDigiProducer") << " Hit " << count << " has tof " << cfi->timeOfFlight() << " trackid "
                                      << cfi->trackId() << " bunchcr " << cfi.bunch() << " trigger " << cfi.getTrigger()
                                      << ", from EncodedEventId: " << cfi->eventId().bunchCrossing() << " "
                                      << cfi->eventId().event() << " bcr from MixCol " << cfi.bunch() << "\n";
       edm::LogInfo("RPDigiProducer") << " Hit: " << (*cfi) << "\n";
       count++;
     }
   }
 
   MixCollection<PSimHit> allRPHits{cf.product(), std::pair(0, 0)};
 
   if (verbosity_)
     edm::LogInfo("RPDigiProducer") << "Input MixCollection size = " << allRPHits.size() << "\n";
 
   //Loop on PSimHit
   simhit_map simHitMap_;
   simHitMap_.clear();
 
   MixCollection<PSimHit>::iterator isim;
   for (isim = allRPHits.begin(); isim != allRPHits.end(); ++isim) {
     simHitMap_[(*isim).detUnitId()].push_back((*isim));
   }
 
   // Step B: LOOP on hits in event
   std::vector<edm::DetSet<TotemRPDigi>> DigiVector;
   DigiVector.reserve(400);
   DigiVector.clear();
 
   for (simhit_map_iterator it = simHitMap_.begin(); it != simHitMap_.end(); ++it) {
     edm::DetSet<TotemRPDigi> digi_collector(it->first);
 
     if (theAlgoMap.find(it->first) == theAlgoMap.end()) {
       theAlgoMap[it->first] =
           std::unique_ptr<RPDetDigitizer>(new RPDetDigitizer(conf_, *rndEngine_, it->first, iSetup));
     }
 
     std::vector<int> input_links;
     simromanpot::DigiPrimaryMapType output_digi_links;
 
     (theAlgoMap.find(it->first)->second)
         ->run(simHitMap_[it->first], input_links, digi_collector.data, output_digi_links);
 
     if (!digi_collector.data.empty()) {
       DigiVector.push_back(convertRPStripDetSet(digi_collector));
     }
   }
 
   // Step C: create empty output collection
   std::unique_ptr<edm::DetSetVector<TotemRPDigi>> digi_output(new edm::DetSetVector<TotemRPDigi>(DigiVector));
 
   if (verbosity_) {
     edm::LogInfo("RPDigiProducer") << "digi_output->size()=" << digi_output->size() << "\n";
   }
   // Step D: write output to file
   iEvent.put(std::move(digi_output));
 }

Member Data Documentation

edm::ParameterSet RPDigiProducer::conf_

private

Definition at line 71 of file RPDigiProducer.cc.

Referenced by produce().

DeadChannelsManager RPDigiProducer::deadChannelsManager

private

this variable answers the question whether given channel is dead or not

Definition at line 80 of file RPDigiProducer.cc.

Referenced by beginRun().

CLHEP::HepRandomEngine* RPDigiProducer::rndEngine_ = 0

private

Definition at line 74 of file RPDigiProducer.cc.

Referenced by produce().

std::vector<std::string> RPDigiProducer::RP_hit_containers_

private

Definition at line 67 of file RPDigiProducer.cc.

Referenced by RPDigiProducer().

bool RPDigiProducer::simulateDeadChannels

private

this variable indicates whether we take into account dead channels or simulate as if all channels work ok (by default we do not simulate dead channels)

Definition at line 85 of file RPDigiProducer.cc.

Referenced by beginRun(), and RPDigiProducer().

std::map<RPDetId, std::unique_ptr<RPDetDigitizer> > RPDigiProducer::theAlgoMap

private

Definition at line 72 of file RPDigiProducer.cc.

Referenced by produce().

edm::EDGetTokenT<CrossingFrame<PSimHit> > RPDigiProducer::tokenCrossingFrameTotemRP

private

Definition at line 87 of file RPDigiProducer.cc.

Referenced by RPDigiProducer().

int RPDigiProducer::verbosity_

private