MPulseGenerator Class Reference
[Readers]

Inheritance diagram for MPulseGenerator:

List of all members.

---

Detailed Description

Pulse generator.

Author:

Marco Vignati

Pulse generator can generate pulses with noise in different channels. This allows to have at the same time many channels with different features to test algorithms. Random generation involves time interval between pulses, amplitude and noise. See details in the configuration file description below.
Some material, together with examples can be found here: http://corleone.roma1.infn.it/~vignatim/cuore/generator/

• Config example:

      reader PulseGenerator
      verbosity = info
      enable = true
      RunNumber = 1
      EventNumberPrintStep = 1000
      Seed = 12345
      ConfigurationFile = genconf.txt
      RoundToADC = true
      NoiseFile = avg_ps.root
      SpectrumFile = spectrum.root
      GenerateFromEnergy = false
      CalibrationFile = cal_ds25.txt
      endmod
      

• Parameters
  • RunNumber [ unsigned int = 1]: run number (the generator will add 700000 to it)
  • Seed [unsigned int = 12345]: seed of the random number generator
  • RoundToADC [bool = true]: round samples to the Old-DAQ ADC resolution
  • ConfigurationFile [string = GeneratorConfiguration.txt]: configuration file for the generator, see description below
  • NoiseFile [string = avg_ps.root] if some channel needs noise generation from an input PS, that is fetched from this file
  • SpectrumFile [string = spectrum.root] if some channel needs amplitude sampled by an input spectrum, that is fetched from this file
  • GenerateFromEnergy [bool = false] generate amplitude from energy distribution that will be then converted into mV through the calibration function and the coefficients contained in the CalibrationFile
  • CalibrationFile [string = ""] file containing the calibration function and its coefficents for each channel. This file is the same produced by MCalCoefficients

• Configuration file:
  This is a sample config file for some Cuoricino channel.

      # begin BoloGeneratorConfig
      # QMatrix
      # 4 22
      # ch   N   f_s  Bes_c R_S     R_L c_p     Alpha    tau_1    tau_2    tau_r       base    pre_trig a_distr a_d_p1 a_d_p2 t_distr t_d_p1 t_d_p2 n_distr n_d_p1 n_d_p2
         1   512 125  12    167.22  54  329.317 0.748827 0.304307 2.23273  0.00617952  711.14  0.958763    4   975.566    0      2     0.02    0       1      -1     30
        12   512 125  12    35.4399 54  320.072 0.879846 0.23286  0.850443 0.019176    961.961 0.953571    4   1701.61    0      2     0.02    0       1      -1     30
        51   512 125  12    38.4842 54  844.168 0.995524 0.340698 5.64396  0.0319884   829.161 0.959332    4   1486.05    0      2     0.02    0       1      -1     30
        63   512 125  12    45.0786 54  625.471 0.853619 0.180573 0.521086 0.0279043   924.812 0.955147    4   1111.93    0      2     0.02    0       1      -1     30
      # end BoloGeneratorConfig
      

where:

• ch: channel number
• N: number of samples per generated window
• f_s: sampling frequency (Hz)
• Bes_c: Bessel filter cutoff frequency(Hz)
• R_S: work resistance (MOhm)
• R_L: load resistance (GOhm)
• c_p: wires capacity (pF)
• Alpha: ratio between the amplitudes of two decay constants (range: 0-1)
• tau_1: first decay constant (s)
• tau_2: second decay constant (s)
• tau_r: rise time (s)
• base: baseline (mV)
• pre_trig: pre trigger region size (s)
• a_distr: amplitude distribution, types are:
  • 0. no pulse generation
  • 1. uniform distribution between a_d_p1 and a_d_p2 (mV or keV)
  • 2. exponential distribution with average a_d_p1 (mV or keV)
  • 3. fixed amplitudes of values a_d_p1 (mV or keV)
  • 4. generation from input histogram, needs Spectrum_xx TH1D's int the input file SpectrumFile. note that if the input spectrum is in keV you must use GenerateFromEnergy = true and provide a calibration file to get the correct amplitude in mV.
• t_distr: distribution of time interval between events, types are:
  • 1. uniform distribution between t_d_p1 and t_d_p2, both in seconds.
  • 2. exponential distribution with average 1/t_d_p1 (this generate a poisson process with rate t_d_p1 (Hz))
  • 3. fixed interval at t_d_p1(s)
• n_distr: noise distribution, types are:
  • 0. no noise generation
  • 1. generation from input Power Spectrum with randomization parameter n_d_p2 (default to 30, do not change it unless you now what you are doing), and average rms n_d_p1(mV), if n_d_p1 < 0 the measured noise RMS is used: this needs NoiseAvgPS_xx QVectors in the input file "NoiseFile",
  • 2. white noise with rms n_d_p1(mV)
  • 3. sinusoidal noise with amplitude n_d_p1(mV) and frequency n_d_p2(Hz)
  • 4. white noise with rms n_d_p1(mV) filtered with the bessel, plus a random noise after the bessel with rms n_d_p2(mV)
  • 6. pink noise (1/f) with rms n_d_p1(mV)
  • 7. pink noise (1/f) with rms n_d_p1(mV) filtered with the bessel

• Reads:
  • NOTHING

• Writes
  • QRawEvent There is no trigger in the generator, thus the pulse(s) start position is saved in the trigger information (QSampleInfo) of the Event.
  • Main Pulse MC truth is saved in the following AuxData:
    • PG_Max: Pulse Maximum (mV)
    • PG_MaxPos: Pulse Maximum Position (s)
    • PG_Energy: Pulse Energy (keV) (if GenerateFromEnergy = true)
    • PG_DeltaT: Time since last event on this channel (s)
  • If other pulses are present in the window, the MC truth variables of each of them (except the main pulse) are saved in the following AuxData (QVector):
    • PG_OE_Max
    • PG_OE_MaxPos
    • PG_OE_Energy
    • PG_OE_DeltaT
    • PG_OE_Time: absolute time (s), the time of the event corresponding to this pulse.

Public Member Functions
	MPulseGenerator (QSequence *s)
	constructor
	~MPulseGenerator ()
	destructor
void	Init ()
QEvent *	Do (QEvent *ev)
	Doit method is called for each event, getting the event as argument.
void	Done ()
	Done method is called after event loop.
Static Private Member Functions
static bool	ComparePulseGenerator (QPulseGenerator *one, QPulseGenerator *two)
Private Attributes
std::vector< QPulseGenerator * >	fChannelGenerators
std::string	fConfigurationFile
int	fSeed
time_t	fStartTime
gsl_rng *	fRandomGen
bool	fRoundToADC
bool	fGenerateFromEnergy
int	fThisRun
int	fRunNumber
double	fApolloADC2mV
double	fOldDAQADC2mV
double *	fADC2mV
QRunData *	fRunData
QVector	fTimes
QVector	fDeltaTs
QVector	fEnergies
QVector	fAmplitudes
QVector	fMax
QVector	fMaxPos

---

Member Function Documentation

void MPulseGenerator::Init

(

)

[virtual]

Init method is called before event loop

Implements QReader.

---

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

• /usr/local/cuore/pkg/modrespfunc/MPulseGenerator.hh
• /usr/local/cuore/pkg/modrespfunc/MPulseGenerator.cc

---