SBNDNoiseServiceFromHist.h

Go to the documentation of this file.

// SBNDNoiseServiceFromHist.h
// Andrew Scarff (University of Sheffield)
// Aug 2019

// Based upon SPhaseChannelNoiseService.h created by Jingbo Wang (UC Davis) for ProtoDUNE.9

// Implementation of a general TPC channel noise model with:
// (1) white noise
// (2) Inherent Gaussian noise in frequency
// (3) MicroBooNE noise in frequency
// (4) Coherent noise (exponential + Gaussian) in frequency 
//     (Note a: phase at each frequency bin is randamized at the moment. Will be updated soon
//      Note b: Currently, consecutive offline channels (configurable) are grouped together and 
//              the same coherent noise waveform is assigned to channels within the same group. )
//
// The default parameters are obtained from the ProtoDUNE-SP data (run 4096)
// fcl file: sbndcode/DetectorSim/Services/SBND_detsim_data_driven_noise.fcl
//

#ifndef SBNDNoiseServiceFromHist_H
#define SBNDNoiseServiceFromHist_H

#include "sbndcode/DetectorSim/Services/ChannelNoiseService.h"

#include "art/Framework/Services/Registry/ServiceHandle.h"
#include "art/Framework/Services/Registry/ServiceDeclarationMacros.h"
#include "art_root_io/TFileService.h"
#include "lardata/DetectorInfoServices/DetectorPropertiesService.h"
#include "lardata/Utilities/LArFFT.h"
#include "larcore/Geometry/Geometry.h"
#include "nurandom/RandomUtils/NuRandomService.h"
#include "sbndcode/Utilities/SignalShapingServiceSBND.h"

#include "CLHEP/Random/JamesRandom.h"
#include "CLHEP/Random/RandFlat.h"

#include "TH1F.h"
#include "TRandom3.h"
#include "TF1.h"
#include "TMath.h"

#include <vector>
#include <iostream>
#include <sstream>

class TH1;
namespace CLHEP {
class HepRandomEngine;
}

class SBNDNoiseServiceFromHist : public ChannelNoiseService {

public:

  // Ctor.
  SBNDNoiseServiceFromHist(fhicl::ParameterSet const& pset);

  // Ctor.
  SBNDNoiseServiceFromHist(fhicl::ParameterSet const& pset, art::ActivityRegistry&);

  // Dtor.
  ~SBNDNoiseServiceFromHist();

  // Add noise to a signal array.
  int addNoise(Channel chan, AdcSignalVector& sigs) const;

  // Print the configuration.
  std::ostream& print(std::ostream& out =std::cout, std::string prefix ="") const;

private:
 
  // General parameters
  unsigned int fNoiseArrayPoints;  ///< number of points in randomly generated noise array
  int          fRandomSeed;        ///< Seed for random number service. If absent or zero, use SeedSvc.
  int          fLogLevel;          ///< Log message level: 0=quiet, 1=init only, 2+=every event

  CLHEP::HepRandomEngine* m_pran;


};

DECLARE_ART_SERVICE_INTERFACE_IMPL(SBNDNoiseServiceFromHist, ChannelNoiseService, LEGACY)

#endif