#include <SBNDThermalNoiseServiceInTime.h>

Inheritance diagram for SBNDThermalNoiseServiceInTime:

Public Member Functions
	SBNDThermalNoiseServiceInTime (fhicl::ParameterSet const &pset)

	SBNDThermalNoiseServiceInTime (fhicl::ParameterSet const &pset, art::ActivityRegistry &)

	~SBNDThermalNoiseServiceInTime ()

int	addNoise (detinfo::DetectorClocksData const &clockData, Channel chan, AdcSignalVector &sigs) const override

std::ostream &	print (std::ostream &out=std::cout, std::string prefix="") const override

Public Member Functions inherited from ChannelNoiseService
virtual	~ChannelNoiseService ()=default

virtual void	generateNoise (detinfo::DetectorClocksData const &)

Private Member Functions
void	InitialiseProducerDeps (art::EDProducer *EDProdPointer, fhicl::ParameterSet const &pset) override

Private Attributes
unsigned int	fNoiseArrayPoints
	number of points in randomly generated noise array More...

int	fRandomSeed
	Seed for random number service. If absent or zero, use SeedSvc. More...

int	fLogLevel
	Log message level: 0=quiet, 1=init only, 2+=every event. More...

std::map< double, int >	fShapingTimeOrder

CLHEP::HepRandomEngine *	m_pran

CLHEP::HepRandomEngine *	fNoiseEngine

Additional Inherited Members
Public Types inherited from ChannelNoiseService
typedef unsigned int	Channel

Detailed Description

Definition at line 46 of file SBNDThermalNoiseServiceInTime.h.

Constructor & Destructor Documentation

SBNDThermalNoiseServiceInTime::SBNDThermalNoiseServiceInTime ( fhicl::ParameterSet const & pset )

Definition at line 21 of file SBNDThermalNoiseServiceInTime_service.cc.

                                                             : fRandomSeed(0), fLogLevel(1),
                                                                 m_pran(nullptr), fNoiseEngine(nullptr)
  {
    
   const string myname = "SBNDThermalNoiseServiceInTime::ctor: ";
   fNoiseArrayPoints  = pset.get<unsigned int>("NoiseArrayPoints",1);
   bool haveSeed      = pset.get_if_present<int>("RandomSeed", fRandomSeed);
 
   fShapingTimeOrder = { {0.5, 0 }, {1.0, 1}, {2.0, 2}, {3.0, 3} };
 
   if ( fRandomSeed == 0 ) haveSeed = false;
   pset.get_if_present<int>("LogLevel", fLogLevel);
   int seed = fRandomSeed;
   
   string rname = "SBNDThermalNoiseServiceInTime";
   if ( haveSeed ) {
     if ( fLogLevel > 0 ) cout << myname << "WARNING: Using hardwired seed." << endl;
     m_pran = new HepJamesRandom(seed);
   } else {
     if ( fLogLevel > 0 ) cout << myname << "Using NuRandomService." << endl;
     art::ServiceHandle<NuRandomService> seedSvc;
     m_pran = new HepJamesRandom;
     if ( fLogLevel > 0 ) cout << myname << "    Initial seed: " << m_pran->getSeed() << endl;
     seedSvc->registerEngine(NuRandomService::CLHEPengineSeeder(m_pran), rname);
   }
   if ( fLogLevel > 0 ) cout << myname << "  Registered seed: " << m_pran->getSeed() << endl;
   if ( fLogLevel > 1 ) print() << endl;
 }

SBNDThermalNoiseServiceInTime::SBNDThermalNoiseServiceInTime	(	fhicl::ParameterSet const &	pset,
		art::ActivityRegistry &
	)

Definition at line 53 of file SBNDThermalNoiseServiceInTime_service.cc.

54 : SBNDThermalNoiseServiceInTime(pset) { }

SBNDThermalNoiseServiceInTime::SBNDThermalNoiseServiceInTime

SBNDThermalNoiseServiceInTime(fhicl::ParameterSet const &pset)

Definition: SBNDThermalNoiseServiceInTime_service.cc:21

SBNDThermalNoiseServiceInTime::~SBNDThermalNoiseServiceInTime ( )

Definition at line 58 of file SBNDThermalNoiseServiceInTime_service.cc.

                                                               {
   const string myname = "SBNDThermalNoiseServiceInTime::dtor: ";
   if ( fLogLevel > 0 ) {
     cout << myname << "Deleting random engine with seed " << m_pran->getSeed() << endl;
   }
   delete m_pran;
 }

Member Function Documentation

int SBNDThermalNoiseServiceInTime::addNoise	(	detinfo::DetectorClocksData const &	clockData,
		Channel	chan,
		AdcSignalVector &	sigs
	)		const

overridevirtual

Implements ChannelNoiseService.

Definition at line 68 of file SBNDThermalNoiseServiceInTime_service.cc.

                                                                                        {
 
   //Get services.
   art::ServiceHandle<geo::Geometry> geo;
   art::ServiceHandle<util::SignalShapingServiceSBND> sss;
   
   //Generate Noise:
   size_t view = (size_t)geo->View(chan);
   
   double noise_factor;
   auto tempNoiseVec = sss->GetNoiseFactVec();
   double shapingTime = 2.0; //sss->GetShapingTime(chan);
   double asicGain = sss->GetASICGain(chan);
   
   if (fShapingTimeOrder.find( shapingTime ) != fShapingTimeOrder.end() ) {
     noise_factor = tempNoiseVec[view].at( fShapingTimeOrder.find( shapingTime )->second );
     noise_factor *= asicGain/4.7;
   }
   else {
     throw cet::exception("SBNDThermalNoiseServiceInTime_service")
       << "\033[93m"
       << "Shaping Time recieved from signalshapingservices_sbnd.fcl is not one of the allowed values"
       << std::endl
       << "Allowed values: 0.5, 1.0, 2.0, 3.0 us"
       << "\033[00m"
       << std::endl;
   }
 
   CLHEP::RandGaussQ rGauss(*fNoiseEngine, 0.0, noise_factor);
     
 
   //In this case fNoiseFact is a value in ADC counts
   //It is going to be the Noise RMS
   //loop over all bins in "noise" vector
   //and insert random noise value
   
   for (unsigned int i = 0; i < sigs.size(); i++){
     sigs.at(i) = rGauss.fire();
   }
 
   return 0;
 }

void SBNDThermalNoiseServiceInTime::InitialiseProducerDeps	(	art::EDProducer *	EDProdPointer,
		fhicl::ParameterSet const &	pset
	)

inlineoverrideprivatevirtual

Reimplemented from ChannelNoiseService.

Definition at line 79 of file SBNDThermalNoiseServiceInTime.h.

                                                                                                     {
     
     CLHEP::HepRandomEngine& NoiseEngine((art::ServiceHandle<rndm::NuRandomService>{}->createEngine(*EDProdPointer,"HepJamesRandom","noise",pset,"Seed")));
     fNoiseEngine = &NoiseEngine;
     return; 
   } 

ostream & SBNDThermalNoiseServiceInTime::print	(	std::ostream &	out = `std::cout`,
		std::string	prefix = `""`
	)		const

overridevirtual

Implements ChannelNoiseService.

Definition at line 115 of file SBNDThermalNoiseServiceInTime_service.cc.

                                                                                {
   out << prefix << "SBNDThermalNoiseServiceInTime: " << endl;
   
   out << prefix << "          LogLevel: " <<  fLogLevel << endl;
   out << prefix << "        RandomSeed: " <<  fRandomSeed << endl;
   out << prefix << "  NoiseArrayPoints: " << fNoiseArrayPoints << endl;
   
   return out;
 }