SBNDNoiseServiceFromHist Class Reference

#include <SBNDNoiseServiceFromHist.h>

Inheritance diagram for SBNDNoiseServiceFromHist:

Public Member Functions
	SBNDNoiseServiceFromHist (fhicl::ParameterSet const &pset)
	SBNDNoiseServiceFromHist (fhicl::ParameterSet const &pset, art::ActivityRegistry &)
	~SBNDNoiseServiceFromHist ()
int	addNoise (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 int	addNoise (detinfo::DetectorClocksData const &, Channel chan, AdcSignalVector &sigs) const =0
virtual void	generateNoise (detinfo::DetectorClocksData const &)
virtual void	InitialiseProducerDeps (art::EDProducer *EDProdPointer, fhicl::ParameterSet const &pset)

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...
CLHEP::HepRandomEngine *	m_pran

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

Detailed Description

Definition at line 51 of file SBNDNoiseServiceFromHist.h.

Constructor & Destructor Documentation

SBNDNoiseServiceFromHist::SBNDNoiseServiceFromHist

(

fhicl::ParameterSet const &

pset

)

Definition at line 20 of file SBNDNoiseServiceFromHist_service.cc.

  : fRandomSeed(0), fLogLevel(1),  m_pran(nullptr), fNoiseEngine(nullptr)
{
  const string myname = "SBNDNoiseServiceFromHist::ctor: ";
  fNoiseArrayPoints  = pset.get<unsigned int>("NoiseArrayPoints");
  bool haveSeed      = pset.get_if_present<int>("RandomSeed", fRandomSeed);


  fShapingTimeOrder = { {0.5, 0 }, {1.0, 1}, {2.0, 2}, {3.0, 3} };
  fNoiseWidth        = pset.get< double              >("NoiseWidth");
  fNoiseRand         = pset.get< double              >("NoiseRand");
  fLowCutoff         = pset.get< double              >("LowCutoff");
  
  //Getting noise histo
  fNoiseHistoName = p.get< std::string         >("NoiseHistoName");
  
  cet::search_path sp("FW_SEARCH_PATH");
  sp.find_file(p.get<std::string>("NoiseFileFname"), fNoiseFileFname);
  
  TFile in(fNoiseFileFname.c_str(), "READ");
  if (!in.IsOpen()) {
    throw art::Exception(art::errors::FileOpenError)
      << "Could not find Root file '" << fNoiseHistoName
      << "' for noise histogram\n";
  }
  fNoiseHist = (TH1D*) in.Get(fNoiseHistoName.c_str());
  
  if (!fNoiseHist) {
    throw art::Exception(art::errors::NotFound)
      << "Could not find noise histogram '" << fNoiseHistoName
      << "' in Root file '" << in.GetPath() << "'\n";
  }
  // release the histogram from its file; now we own it
  fNoiseHist->SetDirectory(nullptr);
  in.Close();


  if ( fRandomSeed == 0 ) haveSeed = false;
  pset.get_if_present<int>("LogLevel", fLogLevel);
  int seed = fRandomSeed;
  
  string rname = "SBNDNoiseServiceFromHist";
  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;
  generateNoise();
  if ( fLogLevel > 1 ) print() << endl;
}

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

Definition at line 80 of file SBNDNoiseServiceFromHist_service.cc.

: SBNDNoiseServiceFromHist(pset) { }

SBNDNoiseServiceFromHist::~SBNDNoiseServiceFromHist

(

)

Definition at line 85 of file SBNDNoiseServiceFromHist_service.cc.

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

Member Function Documentation

int SBNDNoiseServiceFromHist::addNoise	(	Channel	chan,
		AdcSignalVector &	sigs
	)		const

Definition at line 95 of file SBNDNoiseServiceFromHist_service.cc.

                                                                                {

  //Get services.
  art::ServiceHandle<geo::Geometry> geo;
  art::ServiceHandle<util::SignalShapingServiceSBND> sss;
  art::ServiceHandle<util::LArFFT> fFFT;
  
  //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);

  size_t fNTicks = fFFT->FFTSize();

  
  if (fShapingTimeOrder.find( shapingTime ) != fShapingTimeOrder.end() ) {
    noise_factor = tempNoiseVec[view].at( fShapingTimeOrder.find( shapingTime )->second );
    noise_factor *= asicGain/4.7;
  }
  else {
    throw cet::exception("SBNDNoiseServiceFromHist_service.cc")
      << "\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::RandFlat flat(*fNoiseEngine, -1, 1);


  if (sigs.size() != fNTicks)
    throw cet::exception("SBNDNoiseServiceFromHist_service.cc")
        << "\033[93m"
        << "Frequency noise vector length must match fNTicks (FFT size)"
        << " ... " << sigs.size() << " != " << fNTicks
        << "\033[00m"
        << std::endl;

  // noise in frequency space
  std::vector<TComplex> noiseFrequency(fNTicks / 2 + 1, 0.);

  double pval = 0.;
  double lofilter = 0.;
  double phase = 0.;
  double rnd[2] = {0.};

  // width of frequencyBin in kHz
  double binWidth = 1.0 / (fNTicks * fSampleRate * 1.0e-6);

  for (size_t i = 0; i < fNTicks / 2 + 1; ++i) {
    // exponential noise spectrum
    flat.fireArray(2, rnd, 0, 1);
    pval = fNoiseHist->GetBinContent(i) * ((1 - fNoiseRand) + 2 * fNoiseRand * rnd[0]) * noise_factor;
    phase = rnd[1] * 2.*TMath::Pi();
    TComplex tc(pval * cos(phase), pval * sin(phase));
    noiseFrequency.at(i) += tc;
  }


  // inverse FFT MCSignal
  fFFT->DoInvFFT(noiseFrequency, sigs);

  noiseFrequency.clear();

  // multiply each noise value by fNTicks as the InvFFT
  // divides each bin by fNTicks assuming that a forward FFT
  // has already been done.
  for (unsigned int i = 0; i < sigs.size(); ++i) {
    sigs.at(i) *= 1.*fNTicks;
  }

  return 0;
}

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

virtual

Implements ChannelNoiseService.

Definition at line 178 of file SBNDNoiseServiceFromHist_service.cc.

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

Member Data Documentation

int SBNDNoiseServiceFromHist::fLogLevel

private

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

Definition at line 75 of file SBNDNoiseServiceFromHist.h.

unsigned int SBNDNoiseServiceFromHist::fNoiseArrayPoints

private

number of points in randomly generated noise array

Definition at line 73 of file SBNDNoiseServiceFromHist.h.

int SBNDNoiseServiceFromHist::fRandomSeed

private

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

Definition at line 74 of file SBNDNoiseServiceFromHist.h.

CLHEP::HepRandomEngine* SBNDNoiseServiceFromHist::m_pran

private

Definition at line 77 of file SBNDNoiseServiceFromHist.h.

---

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

• SBNDNoiseServiceFromHist.h
• SBNDNoiseServiceFromHist_service.cc