WhiteChannelNoiseService Class Reference

#include <WhiteChannelNoiseService.h>

Inheritance diagram for WhiteChannelNoiseService:

Public Member Functions
	WhiteChannelNoiseService (fhicl::ParameterSet const &pset, art::ActivityRegistry &)
	~WhiteChannelNoiseService ()
int	addNoise (detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, Channel chan, AdcSignalVector &sigs) const
std::ostream &	print (std::ostream &out=std::cout, std::string prefix="") const
Public Member Functions inherited from ChannelNoiseService
virtual	~ChannelNoiseService ()=default
virtual void	newEvent ()

Private Attributes
int	m_RandomSeed
int	m_LogLevel
TH1 *	fNoiseHist
	distribution of noise counts More...
CLHEP::HepRandomEngine *	m_pran

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

Detailed Description

Definition at line 23 of file WhiteChannelNoiseService.h.

Constructor & Destructor Documentation

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

Definition at line 23 of file WhiteChannelNoiseService_service.cc.

: m_RandomSeed(0), m_LogLevel(1),
  fNoiseHist(nullptr),
  m_pran(nullptr) {
  const string myname = "WhiteChannelNoiseService::ctor: ";
  // Fetch parameters.
  bool haveSeed = pset.get_if_present<int>("RandomSeed", m_RandomSeed);
  pset.get_if_present<int>("LogLevel", m_LogLevel);
  if ( m_RandomSeed == 0 ) haveSeed = false;
  // Create random number engine.
  if ( haveSeed ) {
    if ( m_LogLevel > 0 ) cout << myname << "WARNING: Using hardwired seed." << endl;
    m_pran = new HepJamesRandom(m_RandomSeed);
  } else {
    string rname = "WhiteChannelNoiseService";
    if ( m_LogLevel > 0 ) cout << myname << "Using NuRandomService." << endl;
    art::ServiceHandle<NuRandomService> seedSvc;
    m_pran = new HepJamesRandom;
    if ( m_LogLevel > 0 ) cout << myname << "    Initial seed: " << m_pran->getSeed() << endl;
    seedSvc->registerEngine(NuRandomService::CLHEPengineSeeder(m_pran), rname);
  }
  if ( m_LogLevel > 0 ) cout << myname << "  Registered seed: " << m_pran->getSeed() << endl;
  art::ServiceHandle<art::TFileService> tfs;
  fNoiseHist     = tfs->make<TH1F>("Noise", ";Noise  (ADC);", 1000,   -10., 10.);
}

WhiteChannelNoiseService::~WhiteChannelNoiseService

(

)

Definition at line 51 of file WhiteChannelNoiseService_service.cc.

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

Member Function Documentation

int WhiteChannelNoiseService::addNoise ( detinfo::DetectorClocksData const &  clockData, detinfo::DetectorPropertiesData const &  detProp, Channel  chan, AdcSignalVector &  sigs  ) const

virtual

Implements ChannelNoiseService.

Definition at line 61 of file WhiteChannelNoiseService_service.cc.

                                                                                  {
  art::ServiceHandle<util::SignalShapingServiceDUNE> sss;
  float fASICGain      = sss->GetASICGain(chan);
  double fShapingTime   = sss->GetShapingTime(chan);
  std::map< double, int > fShapingTimeOrder;
  fShapingTimeOrder = { {0.5, 0}, {1.0, 1}, {2.0, 2}, {3.0, 3} };
  DoubleVec fNoiseFactVec;
  auto tempNoiseVec = sss->GetNoiseFactVec();
  if ( fShapingTimeOrder.find(fShapingTime) != fShapingTimeOrder.end() ) {
    size_t i = 0;
    fNoiseFactVec.resize(2);
    for ( auto& item : tempNoiseVec ) {
      fNoiseFactVec[i] = item.at(fShapingTimeOrder.find( fShapingTime )->second);
      fNoiseFactVec[i] *= fASICGain/4.7;
      ++i;
    }
  } else {
    throw cet::exception("WhiteChannelNoiseService")
      << "\033[93m"
      << "Shaping Time received from signalservices_dune.fcl is not one of allowed values"
      << std::endl
      << "Allowed values: 0.5, 1.0, 2.0, 3.0 usec"
      << "\033[00m"
      << std::endl;
  }
  art::ServiceHandle<geo::Geometry> geo;
  const geo::View_t view = geo->View(chan);
  CLHEP::HepRandomEngine& engine = *m_pran;
  CLHEP::RandGaussQ rGauss_Ind(engine, 0.0, fNoiseFactVec[0]);
  CLHEP::RandGaussQ rGauss_Col(engine, 0.0, fNoiseFactVec[1]);
  for ( AdcSignal& sig : sigs ) {
    double tnoise = 0.0;
    if      ( view==geo::kU ) tnoise = rGauss_Ind.fire();
    else if ( view==geo::kV ) tnoise = rGauss_Ind.fire();
    else                      tnoise = rGauss_Col.fire();
    sig += tnoise;
  }
  return 0;
}

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

virtual

Implements ChannelNoiseService.

Definition at line 105 of file WhiteChannelNoiseService_service.cc.

                                                                          {
  string myprefix = prefix + "  ";
  out << myprefix << "      Random seed: " << m_pran->getSeed() << endl;
  return out;
}

Member Data Documentation

TH1* WhiteChannelNoiseService::fNoiseHist

private

distribution of noise counts

Definition at line 48 of file WhiteChannelNoiseService.h.

int WhiteChannelNoiseService::m_LogLevel

private

Definition at line 45 of file WhiteChannelNoiseService.h.

CLHEP::HepRandomEngine* WhiteChannelNoiseService::m_pran

private

Definition at line 50 of file WhiteChannelNoiseService.h.

int WhiteChannelNoiseService::m_RandomSeed

private

Definition at line 44 of file WhiteChannelNoiseService.h.

---

The documentation for this class was generated from the following files:

• dunesim/dunesim/DetSim/Service/WhiteChannelNoiseService.h
• dunesim/dunesim/DetSim/Service/WhiteChannelNoiseService_service.cc