#include <RPDetDigitizer.h>

Public Member Functions
	RPDetDigitizer (const edm::ParameterSet &params, CLHEP::HepRandomEngine &eng, RPDetId det_id, const CTPPSRPAlignmentCorrectionsData *alignments, const CTPPSGeometry &geom)

void	run (const std::vector< PSimHit > &input, const std::vector< int > &input_links, std::vector< TotemRPDigi > &output_digi, simromanpot::DigiPrimaryMapType &output_digi_links)

Private Attributes
RPDetId	det_id_

bool	links_persistence_

bool	misalignment_simulation_on_

bool	noNoise_

int	numStrips_

double	theNoiseInElectrons

std::unique_ptr< RPDisplacementGenerator >	theRPDisplacementGenerator

std::unique_ptr< RPGaussianTailNoiseAdder >	theRPGaussianTailNoiseAdder

std::unique_ptr< RPHitChargeConverter >	theRPHitChargeConverter

std::unique_ptr< RPPileUpSignals >	theRPPileUpSignals

std::unique_ptr< RPVFATSimulator >	theRPVFATSimulator

double	theStripThresholdInE

int	verbosity_

Detailed Description

Definition at line 26 of file RPDetDigitizer.h.

Constructor & Destructor Documentation

◆ RPDetDigitizer()

RPDetDigitizer::RPDetDigitizer	(	const edm::ParameterSet &	params,
		CLHEP::HepRandomEngine &	eng,
		RPDetId	det_id,
		const CTPPSRPAlignmentCorrectionsData *	alignments,
		const CTPPSGeometry &	geom
	)

Definition at line 9 of file RPDetDigitizer.cc.

References PixelBaryCentreAnalyzer_cfg::alignments, det_id_, RPTopology::DetStripNo(), relativeConstraints::geom, links_persistence_, misalignment_simulation_on_, noNoise_, numStrips_, submitPVValidationJobs::params, theNoiseInElectrons, theRPDisplacementGenerator, theRPGaussianTailNoiseAdder, theRPHitChargeConverter, theRPPileUpSignals, theRPVFATSimulator, theStripThresholdInE, and verbosity_.

     : det_id_(det_id) {
   verbosity_ = params.getParameter<int>("RPVerbosity");
   numStrips_ = RPTopology().DetStripNo();
   theNoiseInElectrons = params.getParameter<double>("RPEquivalentNoiseCharge300um");
   theStripThresholdInE = params.getParameter<double>("RPVFATThreshold");
   noNoise_ = params.getParameter<bool>("RPNoNoise");
   misalignment_simulation_on_ = params.getParameter<bool>("RPDisplacementOn");
   links_persistence_ = params.getParameter<bool>("RPDigiSimHitRelationsPresistence");
 
   theRPGaussianTailNoiseAdder = std::make_unique<RPGaussianTailNoiseAdder>(
       numStrips_, theNoiseInElectrons, theStripThresholdInE, eng, verbosity_);
   theRPPileUpSignals = std::make_unique<RPPileUpSignals>(params, det_id_);
   theRPVFATSimulator = std::make_unique<RPVFATSimulator>(params, det_id_);
   theRPHitChargeConverter = std::make_unique<RPHitChargeConverter>(params, eng, det_id_);
   theRPDisplacementGenerator = std::make_unique<RPDisplacementGenerator>(
       params.getParameter<bool>("RPDisplacementOn"), det_id_, alignments, geom);
 }

Member Function Documentation

◆ run()

void RPDetDigitizer::run	(	const std::vector< PSimHit > &	input,
		const std::vector< int > &	input_links,
		std::vector< TotemRPDigi > &	output_digi,
		simromanpot::DigiPrimaryMapType &	output_digi_links
	)

Definition at line 33 of file RPDetDigitizer.cc.

References det_id_, mps_fire::i, input, links_persistence_, LogDebug, misalignment_simulation_on_, noNoise_, theRPDisplacementGenerator, theRPGaussianTailNoiseAdder, theRPHitChargeConverter, theRPPileUpSignals, theRPVFATSimulator, and verbosity_.

                                                                            {
   if (verbosity_)
     LogDebug("RPDetDigitizer ") << det_id_ << " received input.size()=" << input.size() << "\n";
   theRPPileUpSignals->reset();
 
   bool links_persistence_checked = links_persistence_ && input_links.size() == input.size();
 
   int input_size = input.size();
   for (int i = 0; i < input_size; ++i) {
     simromanpot::strip_charge_map the_strip_charge_map;
     if (misalignment_simulation_on_)
       the_strip_charge_map = theRPHitChargeConverter->processHit(theRPDisplacementGenerator->displace(input[i]));
     else
       the_strip_charge_map = theRPHitChargeConverter->processHit(input[i]);
 
     if (verbosity_)
       LogDebug("RPHitChargeConverter ") << det_id_ << " returned hits=" << the_strip_charge_map.size() << "\n";
     if (links_persistence_checked)
       theRPPileUpSignals->add(the_strip_charge_map, input_links[i]);
     else
       theRPPileUpSignals->add(the_strip_charge_map, 0);
   }
 
   const simromanpot::strip_charge_map &theSignal = theRPPileUpSignals->dumpSignal();
   simromanpot::strip_charge_map_links_type &theSignalProvenance = theRPPileUpSignals->dumpLinks();
   simromanpot::strip_charge_map afterNoise;
   if (noNoise_)
     afterNoise = theSignal;
   else
     afterNoise = theRPGaussianTailNoiseAdder->addNoise(theSignal);
 
   theRPVFATSimulator->ConvertChargeToHits(afterNoise, theSignalProvenance, output_digi, output_digi_links);
 }