#include <ShapedCohProtoDUNENoiseService.h>

Inheritance diagram for ShapedCohProtoDUNENoiseService:

Public Member Functions
	ShapedCohProtoDUNENoiseService (fhicl::ParameterSet const &)

	ShapedCohProtoDUNENoiseService (fhicl::ParameterSet const &, art::ActivityRegistry &)

	~ShapedCohProtoDUNENoiseService ()

int	addNoise (detinfo::DetectorClocksData const &, detinfo::DetectorPropertiesData const &, Channel, AdcSignalVector &) const

void	newEvent ()

std::ostream &	print (std::ostream &, std::string) const

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

Protected Member Functions
int	addFEMBNoise (const Channel, AdcSignalVector &, detinfo::DetectorClocksData const &) const

int	addHVNoise (const Channel, AdcSignalVector &, detinfo::DetectorClocksData const &) const

int	addShapedNoise (const Channel, AdcSignalVector &, detinfo::DetectorClocksData const &) const

Private Member Functions
void	init ()

Private Attributes
CLHEP::HepRandomEngine *	m_pran

art::ServiceHandle< geo::Geometry >	m_geo

double	m_collection_plane_noise

double	m_induction_plane_noise

double	m_collection_plane_noise_rms

double	m_induction_plane_noise_rms

std::map< size_t, size_t >	m_channel_femb

bool	m_random_phase

double	m_frequency_rms

double	m_amplitude_rms

std::vector< std::vector< double > >	m_FEMBCo_Frq_nominal

std::vector< std::vector< double > >	m_FEMBCo_Amp_nominal

std::vector< std::vector< double > >	m_FEMBCo_Phs_nominal

double	m_FEMBCo_Wfm_this_event_vec [m_n_apa][m_n_femb][m_n_tick]

std::vector< double >	m_HV1_Frq_nominal

std::vector< double >	m_HV1_Amp_nominal

std::vector< double >	m_HV1_Phs_nominal

double	m_HV1_Wfm_this_event_vec [m_n_tick]

bool	m_init = 0

Static Private Attributes
static const size_t	m_n_apa = 12

static const size_t	m_n_tick = 8192

static const size_t	m_n_wire_per_apa = 2560

static const size_t	m_n_femb = 20

static const size_t	m_n_max_coh_noise = 9

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

Detailed Description

Definition at line 15 of file ShapedCohProtoDUNENoiseService.h.

Constructor & Destructor Documentation

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

Definition at line 19 of file ShapedCohProtoDUNENoiseService_service.cc.

                                                                                             {
   
   art::ServiceHandle<rndm::NuRandomService> seedSvc;
   m_pran = new CLHEP::HepJamesRandom();
   seedSvc->registerEngine(rndm::NuRandomService::CLHEPengineSeeder(m_pran), "ShapedCohProtoDUNENoiseService_rand");
   
   std::cout << "Initialising ShapedCohProtoDUNENoiseService\n";
   double amplitude_col = pset.get<double>("amplitude_multiplicator_collection");
   double amplitude_ind = pset.get<double>("amplitude_multiplicator_induction");
   m_collection_plane_noise = pset.get<double>("collection_plane_noise") * amplitude_col;
   m_induction_plane_noise  = pset.get<double>("induction_plane_noise" ) * amplitude_ind;
   m_collection_plane_noise_rms = pset.get<double>("collection_plane_noise_rms") * amplitude_col;
   m_induction_plane_noise_rms  = pset.get<double>("induction_plane_noise_rms" ) * amplitude_ind;
   
   m_FEMBCo_Frq_nominal = pset.get<std::vector<std::vector<double>>>("FEMBCo_Frq");
   m_FEMBCo_Amp_nominal = pset.get<std::vector<std::vector<double>>>("FEMBCo_Amp");
   m_FEMBCo_Phs_nominal = pset.get<std::vector<std::vector<double>>>("FEMBCo_Phs");
   m_random_phase = pset.get<bool>("random_phase_noise");
   
   assert(m_FEMBCo_Frq_nominal.size() == m_n_femb);
   assert(m_FEMBCo_Amp_nominal.size() == m_n_femb);
   assert(m_FEMBCo_Phs_nominal.size() == m_n_femb);
 
   for (size_t femb=0; femb<m_n_femb; ++femb) {
     assert(m_FEMBCo_Frq_nominal[femb].size() <= m_n_max_coh_noise);
     assert(m_FEMBCo_Amp_nominal[femb].size() <= m_n_max_coh_noise);
     assert(m_FEMBCo_Phs_nominal[femb].size() <= m_n_max_coh_noise);
   }
   
   // HV noise
   m_HV1_Frq_nominal = pset.get<std::vector<double>>("HV1_Frq");
   m_HV1_Amp_nominal = pset.get<std::vector<double>>("HV1_Amp");
   m_HV1_Phs_nominal = pset.get<std::vector<double>>("HV1_Phs");
   assert(m_HV1_Frq_nominal.size() <= m_n_max_coh_noise);
   assert(m_HV1_Amp_nominal.size() <= m_n_max_coh_noise);
   assert(m_HV1_Phs_nominal.size() <= m_n_max_coh_noise);
   m_init = false;
 }

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

Definition at line 74 of file ShapedCohProtoDUNENoiseService_service.cc.

74 : ShapedCohProtoDUNENoiseService(pset) { }

ShapedCohProtoDUNENoiseService::ShapedCohProtoDUNENoiseService

ShapedCohProtoDUNENoiseService(fhicl::ParameterSet const &)

Definition: ShapedCohProtoDUNENoiseService_service.cc:19

ShapedCohProtoDUNENoiseService::~ShapedCohProtoDUNENoiseService ( )

Definition at line 76 of file ShapedCohProtoDUNENoiseService_service.cc.

76 { }

Member Function Documentation

int ShapedCohProtoDUNENoiseService::addFEMBNoise	(	const Channel	c,
		AdcSignalVector &	adcs,
		detinfo::DetectorClocksData const &	clock
	)		const

protected

Definition at line 123 of file ShapedCohProtoDUNENoiseService_service.cc.

                                                                                                                                      {
   const size_t femb = m_channel_femb.at(c%m_n_wire_per_apa);
   const size_t apa = c/m_n_wire_per_apa;
 
   for (size_t sample=0; sample<adcs.size(); ++sample)
     adcs[sample] += m_FEMBCo_Wfm_this_event_vec[apa][femb][sample];
   
   return 0;
 }

int ShapedCohProtoDUNENoiseService::addHVNoise	(	const Channel	c,
		AdcSignalVector &	adcs,
		detinfo::DetectorClocksData const &	clock
	)		const

protected

Definition at line 133 of file ShapedCohProtoDUNENoiseService_service.cc.

                                                                                                                                    {
 
   for (size_t sample=0; sample<adcs.size(); ++sample)
     adcs[sample] += m_HV1_Wfm_this_event_vec[sample];
   
   return 0;
 }

int ShapedCohProtoDUNENoiseService::addNoise	(	detinfo::DetectorClocksData const &	clock,
		detinfo::DetectorPropertiesData const &	,
		Channel	c,
		AdcSignalVector &	adcs
	)		const

virtual

Implements ChannelNoiseService.

Definition at line 80 of file ShapedCohProtoDUNENoiseService_service.cc.

                                                                                      {
 
   assert(adcs.size() <= m_n_tick);
   int ret = addFEMBNoise(c, adcs, clock) + addHVNoise(c, adcs, clock) + addShapedNoise(c, adcs, clock);
   return ret;
 }

int ShapedCohProtoDUNENoiseService::addShapedNoise	(	const Channel	c,
		AdcSignalVector &	adcs,
		detinfo::DetectorClocksData const &	clock
	)		const

protected

Definition at line 89 of file ShapedCohProtoDUNENoiseService_service.cc.

                                                                                                                                        {
   
   art::ServiceHandle<geo::Geometry> geo;
   const geo::View_t view = geo->View(c);
   CLHEP::HepRandomEngine& engine = *m_pran;
   CLHEP::RandGauss g_ind(engine, m_induction_plane_noise , m_induction_plane_noise_rms );
   CLHEP::RandGauss g_col(engine, m_collection_plane_noise, m_collection_plane_noise_rms);
   
   double ampl=-1;
   while (ampl<0) {
     switch (view) {
     case geo::kU:
     case geo::kV:
       ampl = g_ind.fire();
       break;
     default:
       ampl = g_col.fire();
       break;
     }
   }
   CLHEP::RandGauss rgauss(engine, 0, ampl);
   
   std::vector<double> noise_vec;
   for (size_t i=0; i<m_n_tick; ++i)
     noise_vec.push_back(rgauss.fire());
   
   art::ServiceHandle<util::SignalShapingServiceDUNE> sss;
   sss->ConvoluteElectronicResponse(clock, c, noise_vec);
   
   std::transform(adcs.begin(), adcs.end(), noise_vec.begin(), adcs.begin(), std::plus<double>());
 
   return 0;
 }

void ShapedCohProtoDUNENoiseService::init ( )

private

Definition at line 58 of file ShapedCohProtoDUNENoiseService_service.cc.

                                           {
   art::ServiceHandle<dune::PdspChannelMapService> channel_map;
   
   std::map<size_t,std::vector<size_t>> femb_se;
   for (size_t chan=0; chan<m_n_wire_per_apa; ++chan) {
     // int n_wire = chan%m_n_wire_per_apa; // that is one little useless statement
     int femb = channel_map->FEMBFromOfflineChannel  (chan)-1; // so that it starts from 0
     int slot = channel_map->SlotIdFromOfflineChannel(chan); // if we ever need this...
     // int femb_asic = channel_map->FEMBChannelFromOfflineChannel(chan) / 32; or this
     femb_se[femb+slot*4].push_back(chan);
     m_channel_femb[chan] = femb+slot*4;
   }
   m_init = true;
   std::cout << "ShapedCohProtoDUNENoiseService initialised!\n";
 }

void ShapedCohProtoDUNENoiseService::newEvent ( )

virtual

Reimplemented from ChannelNoiseService.

Definition at line 143 of file ShapedCohProtoDUNENoiseService_service.cc.

                                               {
   if (not m_init) init();
 
   CLHEP::RandGauss gaus(*m_pran);
   CLHEP::RandFlat flat(*m_pran);
   const double two_pi_tick_period = 6.2831853072/2000000;
   for (size_t apa=0; apa<m_n_apa; ++apa) {
     for (size_t femb=0; femb<m_FEMBCo_Frq_nominal.size(); ++femb) {
       for (size_t noise=0; noise<m_FEMBCo_Frq_nominal[femb].size(); ++noise) {
       
         double frq = -1;
         double amp = -1;
         double phs = -1;
         
         while (frq<0)
           frq = gaus.fire(m_FEMBCo_Frq_nominal[femb][noise], m_frequency_rms);
      
         while (amp<0)
           amp = gaus.fire(m_FEMBCo_Amp_nominal[femb][noise], m_amplitude_rms);
 
         if (m_random_phase) {
           while (phs<0)
             phs = flat.fire(0, 6.2831853072);
         } else {
           throw std::runtime_error("ShapedCohProtoDUNENoiseService::newEvent(): Non random phases are not implemented yet");
         }
         
         if (noise == 0) {
           for (size_t sample=0; sample<m_n_tick; ++sample){
             double two_pi_t_f = two_pi_tick_period*sample*frq;
             m_FEMBCo_Wfm_this_event_vec[apa][femb][sample] = amp * sin(phs+two_pi_t_f);
           }
         } else {
           for (size_t sample=0; sample<m_n_tick; ++sample){
             double two_pi_t_f = two_pi_tick_period*sample*frq;
             m_FEMBCo_Wfm_this_event_vec[apa][femb][sample] += amp * sin(phs+two_pi_t_f);
           }
         }
       }
     }
   }
   
   for (size_t noise=0; noise<m_HV1_Frq_nominal.size(); ++noise) {
     double frq = -1;
     double amp = -1;
     double phs = -1;
       
     while (frq<0)
       frq = gaus.fire(m_HV1_Frq_nominal[noise], m_frequency_rms);
      
     while (amp<0)
       amp = gaus.fire(m_HV1_Amp_nominal[noise], m_amplitude_rms);
 
     if (m_random_phase) {
       while (phs<0)
         phs = flat.fire(0, 6.2831853072);
     } else {
       throw std::runtime_error("ShaoedCohProtoDUNENoise::newEvent(): Non random phases are not implemented yet");
     }
         
     if (noise == 0) {
       for (size_t sample=0; sample<m_n_tick; ++sample){
         double two_pi_t_f = two_pi_tick_period*sample*frq;
         m_HV1_Wfm_this_event_vec[sample] = amp * sin(phs+two_pi_t_f);
       }
     } else {
       for (size_t sample=0; sample<m_n_tick; ++sample){
         double two_pi_t_f = two_pi_tick_period*sample*frq;
         m_HV1_Wfm_this_event_vec[sample] += amp * sin(phs+two_pi_t_f);
       }
     }
   }
 }

std::ostream & ShapedCohProtoDUNENoiseService::print	(	std::ostream &	o,
		std::string
	)		const

virtual

Implements ChannelNoiseService.

Definition at line 217 of file ShapedCohProtoDUNENoiseService_service.cc.

                                                                                 {
   return o;
 }

Member Data Documentation

double ShapedCohProtoDUNENoiseService::m_amplitude_rms

private

Definition at line 55 of file ShapedCohProtoDUNENoiseService.h.

std::map<size_t, size_t> ShapedCohProtoDUNENoiseService::m_channel_femb

private

Definition at line 51 of file ShapedCohProtoDUNENoiseService.h.

double ShapedCohProtoDUNENoiseService::m_collection_plane_noise

private

Definition at line 46 of file ShapedCohProtoDUNENoiseService.h.

double ShapedCohProtoDUNENoiseService::m_collection_plane_noise_rms

private

Definition at line 48 of file ShapedCohProtoDUNENoiseService.h.

std::vector<std::vector<double> > ShapedCohProtoDUNENoiseService::m_FEMBCo_Amp_nominal

private

Definition at line 58 of file ShapedCohProtoDUNENoiseService.h.

std::vector<std::vector<double> > ShapedCohProtoDUNENoiseService::m_FEMBCo_Frq_nominal

private

Definition at line 57 of file ShapedCohProtoDUNENoiseService.h.

std::vector<std::vector<double> > ShapedCohProtoDUNENoiseService::m_FEMBCo_Phs_nominal

private

Definition at line 59 of file ShapedCohProtoDUNENoiseService.h.

double ShapedCohProtoDUNENoiseService::m_FEMBCo_Wfm_this_event_vec[m_n_apa][m_n_femb][m_n_tick]

private

Definition at line 62 of file ShapedCohProtoDUNENoiseService.h.

double ShapedCohProtoDUNENoiseService::m_frequency_rms

private

Definition at line 54 of file ShapedCohProtoDUNENoiseService.h.

art::ServiceHandle<geo::Geometry> ShapedCohProtoDUNENoiseService::m_geo

private

Definition at line 44 of file ShapedCohProtoDUNENoiseService.h.

std::vector<double> ShapedCohProtoDUNENoiseService::m_HV1_Amp_nominal

private

Definition at line 66 of file ShapedCohProtoDUNENoiseService.h.

std::vector<double> ShapedCohProtoDUNENoiseService::m_HV1_Frq_nominal

private

Definition at line 65 of file ShapedCohProtoDUNENoiseService.h.

std::vector<double> ShapedCohProtoDUNENoiseService::m_HV1_Phs_nominal

private

Definition at line 67 of file ShapedCohProtoDUNENoiseService.h.

double ShapedCohProtoDUNENoiseService::m_HV1_Wfm_this_event_vec[m_n_tick]

private

Definition at line 70 of file ShapedCohProtoDUNENoiseService.h.

double ShapedCohProtoDUNENoiseService::m_induction_plane_noise

private

Definition at line 47 of file ShapedCohProtoDUNENoiseService.h.

double ShapedCohProtoDUNENoiseService::m_induction_plane_noise_rms

private

Definition at line 49 of file ShapedCohProtoDUNENoiseService.h.

bool ShapedCohProtoDUNENoiseService::m_init = 0

private

Definition at line 72 of file ShapedCohProtoDUNENoiseService.h.

const size_t ShapedCohProtoDUNENoiseService::m_n_apa = 12

staticprivate

Definition at line 37 of file ShapedCohProtoDUNENoiseService.h.

const size_t ShapedCohProtoDUNENoiseService::m_n_femb = 20

staticprivate

Definition at line 40 of file ShapedCohProtoDUNENoiseService.h.

const size_t ShapedCohProtoDUNENoiseService::m_n_max_coh_noise = 9

staticprivate

Definition at line 41 of file ShapedCohProtoDUNENoiseService.h.

const size_t ShapedCohProtoDUNENoiseService::m_n_tick = 8192

staticprivate

Definition at line 38 of file ShapedCohProtoDUNENoiseService.h.

const size_t ShapedCohProtoDUNENoiseService::m_n_wire_per_apa = 2560

staticprivate

Definition at line 39 of file ShapedCohProtoDUNENoiseService.h.

CLHEP::HepRandomEngine* ShapedCohProtoDUNENoiseService::m_pran

private

Definition at line 43 of file ShapedCohProtoDUNENoiseService.h.

bool ShapedCohProtoDUNENoiseService::m_random_phase

private

Definition at line 53 of file ShapedCohProtoDUNENoiseService.h.

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

dunesim/dunesim/DetSim/Service/ShapedCohProtoDUNENoiseService.h
dunesim/dunesim/DetSim/Service/ShapedCohProtoDUNENoiseService_service.cc

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation