icarus_tool::NoiseFromHist Class Reference

Inheritance diagram for icarus_tool::NoiseFromHist:

Public Member Functions
	NoiseFromHist (const fhicl::ParameterSet &pset)
	~NoiseFromHist ()
void	configure (const fhicl::ParameterSet &pset) override
void	nextEvent () override
void	generateNoise (CLHEP::HepRandomEngine &, CLHEP::HepRandomEngine &, icarusutil::TimeVec &, detinfo::DetectorPropertiesData const &, double, const geo::PlaneID &, unsigned int) override

Private Attributes
double	fNoiseRand
std::string	fInputNoiseHistFileName
std::string	fHistogramName
double	fHistNormFactor
std::vector< double >	fNoiseHistVec
std::unique_ptr < icarus_signal_processing::ICARUSFFT < double > >	fFFT

Additional Inherited Members
Private Member Functions inherited from icarus_tool::IGenNoise
virtual	~IGenNoise () noexcept=default

Detailed Description

Definition at line 36 of file NoiseFromHist_tool.cc.

Constructor & Destructor Documentation

icarus_tool::NoiseFromHist::NoiseFromHist ( const fhicl::ParameterSet &  pset )

explicit

Definition at line 72 of file NoiseFromHist_tool.cc.

{
    configure(pset);
}

icarus_tool::NoiseFromHist::~NoiseFromHist

(

)

Definition at line 77 of file NoiseFromHist_tool.cc.

{
}

Member Function Documentation

void icarus_tool::NoiseFromHist::configure ( const fhicl::ParameterSet &  pset )

overridevirtual

Implements icarus_tool::IGenNoise.

Definition at line 81 of file NoiseFromHist_tool.cc.

{
    // Recover the histogram used for noise generation
    fNoiseRand              = pset.get< double>("NoiseRand");
    fInputNoiseHistFileName = pset.get<std::string>("NoiseHistFileName");
    fHistogramName          = pset.get<std::string>("HistogramName");
    fHistNormFactor         = pset.get<double>("HistNormFactor");
    
    std::string fullFileName;
    cet::search_path searchPath("FW_SEARCH_PATH");
    searchPath.find_file(fInputNoiseHistFileName, fullFileName);
    
    TFile inputFile(fullFileName.c_str(), "READ");
    
    if (!inputFile.IsOpen())
        throw cet::exception("NoiseFromHist::configure") << "Unable to open input file: " << fInputNoiseHistFileName << std::endl;
    
    TH1D* histPtr = (TH1D*)inputFile.Get(fHistogramName.c_str());
    
    if (!histPtr)
        throw cet::exception("NoiseFromHist::configure") << "Unable to recover desired histogram: " << fHistogramName << std::endl;
    
    fNoiseHistVec.resize(histPtr->GetNbinsX(), 0.);
    
    for(size_t histIdx = 0; histIdx < size_t(histPtr->GetNbinsX()); histIdx++)
        fNoiseHistVec[histIdx] = histPtr->GetBinContent(histIdx+1);
    
    // Close the input file
    inputFile.Close();

    // Now set up our plans for doing the convolution
    auto const clockData = art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataForJob();
    int numberTimeSamples = clockData.NumberTimeSamples();

    fFFT = std::make_unique<icarus_signal_processing::ICARUSFFT<double>>(numberTimeSamples);
   
    return;
}

void icarus_tool::NoiseFromHist::generateNoise ( CLHEP::HepRandomEngine &  engine, CLHEP::HepRandomEngine &  , icarusutil::TimeVec &  noise, detinfo::DetectorPropertiesData const &  detProp, double  noise_factor, const geo::PlaneID &  planeID, unsigned int  board  )

overridevirtual

Implements icarus_tool::IGenNoise.

Definition at line 120 of file NoiseFromHist_tool.cc.

{    
    CLHEP::RandFlat flat(engine,-1,1);
    
    size_t nFFTTicks = detProp.NumberTimeSamples();
    
    if(noise.size() != nFFTTicks)
        throw cet::exception("SimWireICARUS")
        << "\033[93m"
        << "Frequency noise vector length must match FFT Ticks (FFT size)"
        << " ... " << noise.size() << " != " << nFFTTicks
        << "\033[00m"
        << std::endl;
    
    // noise in frequency space
    std::vector<std::complex<double>> noiseFrequency(nFFTTicks/2+1, 0.);
    
    double pval        = 0.;
    double phase       = 0.;
    double rnd[2]      = {0.};
    double scaleFactor = fHistNormFactor * noise_factor;
    
    // width of frequencyBin in kHz
    
    for(size_t i=0; i< nFFTTicks/2 + 1; ++i)
    {
        // exponential noise spectrum
        flat.fireArray(2,rnd,0,1);
        
        pval = fNoiseHistVec[i] * ((1-fNoiseRand) + 2 * fNoiseRand*rnd[0]) * scaleFactor;

        phase = rnd[1] * 2. * M_PI;

        std::complex<double> tc(pval*cos(phase),pval*sin(phase));

        noiseFrequency.at(i) += tc;
    }
    
    // inverse FFT MCSignal
    fFFT->inverseFFT(noiseFrequency, noise);
    
    noiseFrequency.clear();

    return;
}

void icarus_tool::NoiseFromHist::nextEvent ( )

inlineoverridevirtual

Implements icarus_tool::IGenNoise.

Definition at line 45 of file NoiseFromHist_tool.cc.

{return;};

Member Data Documentation

std::unique_ptr<icarus_signal_processing::ICARUSFFT<double> > icarus_tool::NoiseFromHist::fFFT

private

Definition at line 67 of file NoiseFromHist_tool.cc.

double icarus_tool::NoiseFromHist::fHistNormFactor

private

Definition at line 61 of file NoiseFromHist_tool.cc.

std::string icarus_tool::NoiseFromHist::fHistogramName

private

Definition at line 59 of file NoiseFromHist_tool.cc.

std::string icarus_tool::NoiseFromHist::fInputNoiseHistFileName

private

Definition at line 58 of file NoiseFromHist_tool.cc.

std::vector<double> icarus_tool::NoiseFromHist::fNoiseHistVec

private

Definition at line 65 of file NoiseFromHist_tool.cc.

double icarus_tool::NoiseFromHist::fNoiseRand

private

Definition at line 57 of file NoiseFromHist_tool.cc.

---

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

• NoiseFromHist_tool.cc