Class that handles the channel trigger. More...

#include <ChanTrigger.h>

Public Member Functions
	ChanTrigger ()
	Channel trigger constructur. More...

void	InitializeEachBand (Anita *anita1)
	Initialize trigger bands. More...

void	ApplyAntennaGain (Settings settings1, Anita anita1, Balloon bn1, Screen panel1, int ant, Vector &n_eplane, Vector &n_hplane, Vector &n_normal)
	Apply the antenna gain. More...

void	TriggerPath (Settings settings1, Anita anita1, int ant, Balloon *bn1)
	Apply trigger path. More...

void	DigitizerPath (Settings settings1, Anita anita1, int ant, Balloon *bn1)
	Apply digitizer path. More...

void	TimeShiftAndSignalFluct (Settings settings1, Anita anita1, int ilayer, int ifold, double volts_rx_rfcm_lab_e_all[48][512], double volts_rx_rfcm_lab_h_all[48][512])
	Time shift and fluctuate signal. More...

void	L1Trigger (Anita anita1, double timedomain_output_1[5][Anita::NFOUR], double timedomain_output_2[5][Anita::NFOUR], double powerthreshold[2][5], int channels_passing_e_forglob, int *channels_passing_h_forglob, int &npass)
	The L1 trigger of the Anita trigger scheme. More...

double	getRate ()
	Returns the thisrate variable value (in MHz) More...

double	rateToThreshold (double rate, int band)
	Calculates the trigger threshold for an antenna via a fit from the "singles" rate and band identifier. More...

void	WhichBandsPass (Settings settings1, Anita anita1, GlobalTrigger globaltrig1, Balloon bn1, int ilayer, int ifold, double dangle, double emfrac, double hadfrac, double thresholds[2][5])
	Which bands passes the trigger. More...

void	WhichBandsPassTrigger1 (Settings settings1, Anita anita1, GlobalTrigger globaltrig1, Balloon bn1, int ilayer, int ifold, double thresholds[2][5])
	Which bands passes the trigger (for trigger scheme 0 and 1) More...

void	WhichBandsPassTrigger2 (Settings settings1, Anita anita1, GlobalTrigger globaltrig1, Balloon bn1, int ilayer, int ifold, double dangle, double emfrac, double hadfrac, double thresholds[2][5])
	Which bands passes the trigger (for trigger scheme larger than 2) More...

void	GetThresholds (Settings settings1, Anita anita1, int ilayer, double thresholds[2][5])
	Sets the threshold values based on which payload and where the antenna is located physically. More...

void	DiodeConvolution (Settings settings1, Anita anita1, GlobalTrigger *globaltrig1, int ilayer, int ifold, double mindiodeconvl[5], double onediodeconvl[5], double psignal[5][Anita::NFOUR], double timedomain_output[5][Anita::NFOUR], double timedomain_output_justNoise[5][Anita::NFOUR], int ibinshift, int ipol, double thresholds[2][5])
	Apply the diode convolution. More...

void	addToChannelSums (Settings settings1, Anita anita1, int ibw, int k)
	Increment the volts in each band. More...

void	applyImpulseResponseDigitizer (Settings settings1, Anita anita1, int nPoints, int ant, double *x, double y[512], bool pol)
	Apply impulse response to digitizer path. More...

void	applyImpulseResponseTrigger (Settings settings1, Anita anita1, int ant, double y[512], double *vhz, bool pol)
	Apply impulse response to trigger path. More...

void	getNoiseFromFlight (Settings settings1, Anita anita1, int ant, bool also_digi=true)
	Add noise from ANITA-3 flight to the time domain waveforms. More...

void	injectImpulseAfterAntenna (Anita *anita1, int ant)
	Inject pulse after the antenna (used for trigger efficiency scans) More...

TGraph *	getPulserAtAMPA (Anita *anita1, int ant)
	Get time domain graph of pulse at AMPA (used for trigger efficiency scans) More...

void	saveTriggerWaveforms (Anita *anita1, double sig0[48], double sig1[48], double noise0[48], double noise1[48])
	Save signal and noise waveforms at trigger. More...

void	saveDigitizerWaveforms (Anita *anita1, double sig0[48], double sig1[48], double noise0[48], double noise1[48])
	Save signal and noise waveforms at digitizer. More...

void	injectImpulseAtSurf (Anita *anita1, double volts_triggerPath_e[Anita::HALFNFOUR], double volts_triggerPath_h[Anita::HALFNFOUR], int ant)
	Inject pulse at the surf (used for trigger efficiency scans) More...

void	calculateCW (Anita *anita1, double frequency, double phase, double amplitude)
	Add CW. More...

double	applyButterworthFilter (double ff, double ampl, int notchStatus[3])
	Apply Butterworth Filter. More...

Static Public Member Functions
static void	ConvertEHtoLREfield (double, double, double &, double &)
	Convert E and H to left and right e field. More...

static void	ConvertEHtoLREnergy (double, double, double &, double &)
	Convert E and H to left and right energy. More...

static double	GetNoise (Settings *settings1, double altitude_bn, double geoid, double theta, double bw, double temp)
	Get Noise. More...

static double	FindPeak (double *waveform, int n)
	Find peak voltage of a waveform. More...

static int	IsItUnmasked (unsigned short surfTrigBandMask[9][2], int ibw, int ilayer, int ifold, int ipol)
	Returns whether the indicated antenna and band are "masked". More...

Public Attributes
double	vhz_rx [2][5][Anita::NFREQ]
	Array of amplitudes in the Fourier domain (V/Hz) after the antenna gain. Indeces stand for [ipol][iband][ifreq]. More...

double	volts_rx_forfft [2][5][Anita::HALFNFOUR]
	Array of time domain after the antenna gain. Indeces stand for [ipol][iband][itime]. More...

vector< int >	flag_e [5]
	Which bands pass trigger e. More...

vector< int >	flag_h [5]
	Which bands pass trigger h. More...

double	bwslice_volts_pol0 [5]
	Sum voltage for each slice in bandwidth for the lcp polarization. More...

double	bwslice_volts_pol1 [5]
	Sum voltage for each slice in bandwidth for the rcp polarization. More...

double	bwslice_energy_pol0 [5]
	Square the sum of voltage for each slice in bandwidth for the 0th polarization. More...

double	bwslice_energy_pol1 [5]
	Square the sum of voltage for each slice in bandwidth for the 1st polarization. More...

double	bwslice_volts_pol0_em [5]
	Component of the voltage that comes from the em shower for 0th polarization. More...

double	bwslice_volts_pol1_em [5]
	Component of the voltage that comes from the em shower for 1st polarization. More...

double	bwslice_volts_pole [5]
	Sum voltage for each slice in bandwidth for the e polarization. More...

double	bwslice_energy_pole [5]
	Square the sum of voltage for each slice in bandwidth for e polarization. The 5th element is the full band. More...

double	bwslice_volts_polh [5]
	Sum voltage for each slice in bandwidth for the h polarization. More...

double	bwslice_energy_polh [5]
	Square the sum of voltage for each slice in bandwidth for h polarization. The 5th element is the full band. More...

double	volts_rx_rfcm_lab [2][Anita::HALFNFOUR]
	For digitizer path, time domain voltage vs. time after rx, rfcm's and lab. More...

double	volts_rx_rfcm_lab_all [2][48][Anita::HALFNFOUR]
	For digitizer path, time domain voltage vs. time after rx, rfcm's and lab. More...

double	volts_rx_rfcm [2][Anita::HALFNFOUR]
	For digitizer path, time domain voltage vs. time after rx, rfcm's. More...

double	justNoise_digPath [2][Anita::HALFNFOUR]
	For digitizer path, time domain noise from flight. More...

double	justNoise_trigPath [2][Anita::HALFNFOUR]
	For trigger path, time domain noise from flight. More...

double	cw_digPath [2][Anita::HALFNFOUR]
	For digitizer path, time domain cw. More...

double	justSig_trigPath [2][Anita::HALFNFOUR]
	Just signal in trigger path. More...

double	justSig_digPath [2][Anita::HALFNFOUR]
	Just signal in trigger path. More...

double	signal_eachband [2][Anita::NBANDS_MAX]
	Signal in each band. More...

double	threshold_eachband [2][Anita::NBANDS_MAX]
	Threshold in each band. More...

double	noise_eachband [2][Anita::NBANDS_MAX]
	Noise in each band. More...

int	passes_eachband [2][Anita::NBANDS_MAX]
	Whether the signal passes or not each band. More...

vector< double >	vsignal_eachband [2]
	Signal in each band. More...

vector< double >	vthreshold_eachband [2]
	Threshold in each band. More...

vector< double >	vnoise_eachband [2]
	Noise in each band. More...

vector< int >	vpasses_eachband [2]
	Whether the signal passes or not each band. More...

double	v_banding_rfcm [2][5][Anita::NFREQ]
	This is Volts/m as a function of frequency after rfcm's and banding. More...

double	v_banding_rfcm_forfft [2][5][HALFNFOUR]
	Starts out as V/s vs. freq after banding, rfcm, after fft it is V vs. t. More...

double	vm_banding_rfcm_forfft [2][5][HALFNFOUR]
	Tunnel diode input (signal + noise) More...

double	vm_banding_rfcm_forfft_justNoise [2][5][HALFNFOUR]
	Tunnel diode input (just noise) More...

double	v_banding_rfcm_forfft_temp [2][5][HALFNFOUR]
	Use for the averaging over 10 neighboring bins. More...

double	integral_vmmhz
	Electric field integral. More...

int	unwarned
	Whether we have warned the user about resetting thresholds when they are beyond the measured bounds. More...

Private Member Functions
double	ADCCountstoPowerThreshold (Anita *anita1, int ipol, int iant)

Private Attributes
double	thisrate
	Rate in MHz. More...

double	thispowerthresh
	Relative power threshold. More...

double	e_component
	E comp along polarization. More...

double	h_component
	H comp along polarization. More...

double	n_component
	normal comp along polarization More...

double	e_component_kvector
	component of e-field along the rx e-plane More...

double	h_component_kvector
	component of the e-field along the rx h-plane More...

double	n_component_kvector
	component of the e-field along the normal More...

double	hitangle_e
	angle the ray hits the antenna wrt e-plane More...

double	hitangle_h
	angle the ray hits the antenna wrt h-plane More...

Static Private Attributes
static const int	NSURF =9
	Number of surfs. More...

static const int	NSURFPLUSONE =10
	Number of surfs plus one. More...

static const int	NSURFMINUSONE =8
	Number of surfs minus on1. More...

static const int	NCHANNELS =32
	Number of channells on each surf. More...

static const int	NPOINTS =4073
	Max number of points from surf measurements. More...

static const unsigned	NFOUR = 1024
	Number of points in Fourier space. More...

static const unsigned	HALFNFOUR = 512
	Half of the number of points in the Fourier space. More...

Detailed Description

Class that handles the channel trigger.

Constructor & Destructor Documentation

◆ ChanTrigger()

ChanTrigger::ChanTrigger ( )

Channel trigger constructur.

Member Function Documentation

◆ ADCCountstoPowerThreshold()

double ChanTrigger::ADCCountstoPowerThreshold	(	Anita *	anita1,
		int	ipol,
		int	iant
	)

private

ADC counts to relative power threshold

Function to convert adc thresholds into relative power thresholds that can be handled by icemc So far this only works for ANITA-3 and full band trigger

Parameters

anita1	:: Anita - anita object
ipol	:: int - which polarisation
iant	:: int - which antenna number

Returns: relative power threshold (double)

◆ addToChannelSums()

void ChanTrigger::addToChannelSums	(	Settings *	settings1,
		Anita *	anita1,
		int	ibw,
		int	k
	)

Increment the volts in each band.

This increments the bwslice_volts_...'s at each frequency

Parameters

settings1	:: Settings - simulation settings
anita1	:: Anita - anita payload object
ibw	:: int - which band
k	:: int - frequency bin

◆ ApplyAntennaGain()

void ChanTrigger::ApplyAntennaGain	(	Settings *	settings1,
		Anita *	anita1,
		Balloon *	bn1,
		Screen *	panel1,
		int	ant,
		Vector &	n_eplane,
		Vector &	n_hplane,
		Vector &	n_normal
	)

Apply the antenna gain.

Loop through all points in the screen, and apply the antenna gain to each point All the waveforms are then summed with screen point weights and delay volts_rx_forfft = time domain waveforms after antenna gain vhz_rx = amplitude in Fourier domain after antenna gain

Parameters

settings1	:: Settings - the simulation settings
anita1	:: Anita - anita object
bn1	:: Balloon - balloon object
panel1	:: Screen - screen object
ant	:: int- antenna number
n_eplane	:: Vector
n_hplane	:: Vector
n_normal	:: Vector

◆ applyButterworthFilter()

double ChanTrigger::applyButterworthFilter	(	double	ff,
		double	ampl,
		int	notchStatus[3]
	)

Apply Butterworth Filter.

This is an approximation of the notch filters flown during ANITA4

Parameters

ff	:: frequenzy in Hz
ampl	:: amplitude
filters	:: array of three integers indicating notch status (on/off)

Returns: frequency domain amplitude scaled by the filter

◆ applyImpulseResponseDigitizer()

void ChanTrigger::applyImpulseResponseDigitizer	(	Settings *	settings1,
		Anita *	anita1,
		int	nPoints,
		int	ant,
		double *	x,
		double	y[512],
		bool	pol
	)

Apply impulse response to digitizer path.

This can only work when FFTtools is also linked

Parameters

settings1	:: Settings - simulation settings
anita1	:: Anita - anita payload object
nPoints	:: int - number of points in time domain
ant	:: int - antenna number
x	:: double* - time values
y	:: double[512] - output voltages
pol	:: bool - which polarization

◆ applyImpulseResponseTrigger()

void ChanTrigger::applyImpulseResponseTrigger	(	Settings *	settings1,
		Anita *	anita1,
		int	ant,
		double	y[512],
		double *	vhz,
		bool	pol
	)

Apply impulse response to trigger path.

This can only work when FFTtools is also linked

Parameters

settings1	:: Settings - simulation settings
anita1	:: Anita - anita payload object
ant	:: int - antenna number
y	:: double[512] - output voltages
vhz	:: double* - amplitude in Fourier domain
pol	:: bool - which polarization

◆ calculateCW()

void ChanTrigger::calculateCW	(	Anita *	anita1,
		double	frequency,
		double	phase,
		double	amplitude
	)

Add CW.

Add CW in time domain

Parameters

anita1	:: Anita - anita payload object
frequency	at which to simulate cw
phase
amplitude

◆ ConvertEHtoLREfield()

void ChanTrigger::ConvertEHtoLREfield	(	double	e_component,
		double	h_component,
		double &	lcp_component,
		double &	rcp_component
	)

static

Convert E and H to left and right e field.

Convert E and H to left and right e field

Parameters

e_component	:: double -
h_component	:: double -
lcp_component	:: double - referenced lcp component
rcp_component	:: double - referenced rcp component

◆ ConvertEHtoLREnergy()

void ChanTrigger::ConvertEHtoLREnergy	(	double	e_component,
		double	h_component,
		double &	lcp_component,
		double &	rcp_component
	)

static

Convert E and H to left and right energy.

Convert E and H to left and right energy

Parameters

e_component	:: double -
h_component	:: double -
lcp_component	:: double - referenced lcp component
rcp_component	:: double - referenced rcp component

◆ DigitizerPath()

void ChanTrigger::DigitizerPath	(	Settings *	settings1,
		Anita *	anita1,
		int	ant,
		Balloon *	bn1
	)

Apply digitizer path.

Apply digitizer path to each channel volts_rx_rfcm_lab = time domain waveforms after digitizer path vhz_rx_rfcm_lab_e = amplitude in Fourier domain after digitizer path

Parameters

settings1	:: Settings - the simulation settings
anita1	:: Anita - anita object
ant	:: int - antenna number

◆ DiodeConvolution()

void ChanTrigger::DiodeConvolution	(	Settings *	settings1,
		Anita *	anita1,
		GlobalTrigger *	globaltrig1,
		int	ilayer,
		int	ifold,
		double	mindiodeconvl[5],
		double	onediodeconvl[5],
		double	psignal[5][Anita::NFOUR],
		double	timedomain_output[5][Anita::NFOUR],
		double	timedomain_output_justNoise[5][Anita::NFOUR],
		int	ibinshift,
		int	ipol,
		double	thresholds[2][5]
	)

Apply the diode convolution.

Parameters

settings1	:: Settings - simulation settings
anita1	:: Anita - anita object
globaltrig1	:: GlobalTrigger - global trigger object
ilayer	:: int - layer number
ifold	:: int - phi sector
mindiodeconvl	:: double[5]
onediodeconvl	:: double[5]
psignal	:: double[5][NFOUR]
timedomain_output	:: double[5][NFOUR]
ibinshift	:: int
ipol	:: int - which polarization
thresholds	:: double[2][5] - relative power thresholds for each pol and band

◆ FindPeak()

double ChanTrigger::FindPeak	(	double *	waveform,
		int	n
	)

static

Find peak voltage of a waveform.

Parameters

waveform	:: double* - waveform
n	:: int - number of points in waveform

◆ GetNoise()

double ChanTrigger::GetNoise	(	Settings *	settings1,
		double	altitude_bn,
		double	geoid,
		double	theta,
		double	bw,
		double	temp
	)

static

Get Noise.

Parameters

settings1	:: Settings - the simulation settings
altitude_bn	:: double - altitude of the balloon
geoid	:: double -
bw	:: double -
theta	:: double -
temp	:: double -

◆ getNoiseFromFlight()

void ChanTrigger::getNoiseFromFlight	(	Settings *	settings1,
		Anita *	anita1,
		int	ant,
		bool	also_digi = `true`
	)

Add noise from ANITA-3 flight to the time domain waveforms.

Rayleigh distribution parameters of ANITA-3 thermal noise are read in anita.cc This function generates random noise using amplitude of Rayleigh distributions. This can only work when FFTtools is also linked

Parameters

anita1	:: Anita - anita payload object
pol	:: int - which polarization
ant	:: int - which antennta
also_digi	:: bool - also fill in digitizer

◆ getPulserAtAMPA()

TGraph* ChanTrigger::getPulserAtAMPA	(	Anita *	anita1,
		int	ant
	)

Get time domain graph of pulse at AMPA (used for trigger efficiency scans)

Pulser waveform is read in anita.cc

Parameters

anita1	:: Anita - anita payload object
ant	:: int - which antennta

◆ getRate()

double ChanTrigger::getRate ( )

Returns the thisrate variable value (in MHz)

Todo:: Refactor the class so that this function is either deprecated or "thisrate" is implemented differently.

◆ GetThresholds()

void ChanTrigger::GetThresholds	(	Settings *	settings1,
		Anita *	anita1,
		int	ilayer,
		double	thresholds[2][5]
	)

Sets the threshold values based on which payload and where the antenna is located physically.

The nadir antennas had a separate threshold from the other antennas due to the way that they were "OR"'d between their two neighbors.

Todo:: Deprecate in favor of the more robust boost::multi_array or the more specialized PayloadArray class. Both have multi-index access to the same items.

◆ InitializeEachBand()

void ChanTrigger::InitializeEachBand ( Anita * anita1 )

Initialize trigger bands.

Initialize trigger bands

Parameters

anita1 :: Anita - anita object

◆ injectImpulseAfterAntenna()

void ChanTrigger::injectImpulseAfterAntenna	(	Anita *	anita1,
		int	ant
	)

Inject pulse after the antenna (used for trigger efficiency scans)

Pulser waveform is read in anita.cc THIS FUNCTION IS ONLY USED WITH THE OLD TRIGGER IMPULSE RESPONSE

Parameters

anita1	:: Anita - anita payload object
ant	:: int - which antennta

◆ injectImpulseAtSurf()

void ChanTrigger::injectImpulseAtSurf	(	Anita *	anita1,
		double	volts_triggerPath_e[Anita::HALFNFOUR],
		double	volts_triggerPath_h[Anita::HALFNFOUR],
		int	ant
	)

Inject pulse at the surf (used for trigger efficiency scans)

Pulser waveforms are read in anita.cc volts_triggerPath_e/h is substituded with pulser waveform + noise This can only work when FFTtools is also linked

Parameters

anita1	:: Anita - anita payload object
ant	:: int - which antennta
volts_triggerPath_e	:: double[NFOUR] - time domain waveform at surf
volts_triggerPath_h	:: double[NFOUR] - time domain waveform at surf

◆ IsItUnmasked()

int ChanTrigger::IsItUnmasked	(	unsigned short	surfTrigBandMask[9][2],
		int	ibw,
		int	ilayer,
		int	ifold,
		int	ipol
	)

static

Returns whether the indicated antenna and band are "masked".

Only works for Anita-1 and 2

Parameters

surfTrigBandMask	:: unsigned short [9][2] - surf masks
ibw	:: int - which band
ilayer	:: int - layer number
ifold	:: int - phi sector
ipol	:: int - pol number

Returns: whether the channel and band are masked or not

◆ L1Trigger()

void ChanTrigger::L1Trigger	(	Anita *	anita1,
		double	timedomain_output_1[5][Anita::NFOUR],
		double	timedomain_output_2[5][Anita::NFOUR],
		double	powerthreshold[2][5],
		int *	channels_passing_e_forglob,
		int *	channels_passing_h_forglob,
		int &	npass
	)

The L1 trigger of the Anita trigger scheme.

Parameters

anita1	:: Anita - anita payload object
timedomain_output_1	:: double [4][NFOUR] - time domain output for each band e
timedomain_output_2	:: double [4][NFOUR] - time domain output for each band h
powerthreshold	:: double[2][5] - relative power thresholds for each pol and each band
channels_passing_e_forglob	:: double* - array of channels passing the L1 trigger for e (used in the GlobalTrigger class)
channels_passing_h_forglob	:: double* - array of channels passing the L1 trigger for h (used in the GlobalTrigger class)
npass	:: &int - number of bands passing the L1 trigger

◆ rateToThreshold()

double ChanTrigger::rateToThreshold	(	double	rate,
		int	band
	)

Calculates the trigger threshold for an antenna via a fit from the "singles" rate and band identifier.

This code is stolen from Stephen's AnitaHardwareTrigger in his Aesop software

Parameters

rate	is the desired singles rate, in MHz.
band	ranges from 0 to 3

Returns: trigger threshold, in units of Diode Output / Average Diode Output The output is good for all antennas, and the curves are good fits in the region ~1 MHz to 25 MHz. The curves look like they should continue beyond that range as well.

Todo:: The hard-coded values should be explained and the band dependence should be more explicit, because the actual bands can be modified elsewhere in the code and this function won't be changed unless the user specifically knows to.

◆ saveDigitizerWaveforms()

void ChanTrigger::saveDigitizerWaveforms	(	Anita *	anita1,
		double	sig0[48],
		double	sig1[48],
		double	noise0[48],
		double	noise1[48]
	)

Save signal and noise waveforms at digitizer.

Parameters

anita1	:: Anita - anita payload object
sig	:: double[2][NFOUR/2] - output
noise	:: double[2][NFOUR/2] - output

◆ saveTriggerWaveforms()

void ChanTrigger::saveTriggerWaveforms	(	Anita *	anita1,
		double	sig0[48],
		double	sig1[48],
		double	noise0[48],
		double	noise1[48]
	)

Save signal and noise waveforms at trigger.

Parameters

anita1	:: Anita - anita payload object
sig	:: double[2][NFOUR/2] - output
noise	:: double[2][NFOUR/2] - output

◆ TimeShiftAndSignalFluct()

void ChanTrigger::TimeShiftAndSignalFluct	(	Settings *	settings1,
		Anita *	anita1,
		int	ilayer,
		int	ifold,
		double	volts_rx_rfcm_lab_e_all[48][512],
		double	volts_rx_rfcm_lab_h_all[48][512]
	)

Time shift and fluctuate signal.

Apply time delays to each channel and add thermal noise

Parameters

settings1	:: Settings - the simulation settings
anita1	:: Anita - anita object
ilayer	:: int - payload layer
ifold	:: int - payload phi sector
volts_rx_rfcm_lab_e_all	:: double [48][512] - time domain waveform for each channel (VPOL)
volts_rx_rfcm_lab_h_all	:: double [48][512] - time domain waveform for each channel (HPOL)

◆ TriggerPath()

void ChanTrigger::TriggerPath	(	Settings *	settings1,
		Anita *	anita1,
		int	ant,
		Balloon *	bn1
	)

Apply trigger path.

Apply trigger path to each channel v_banding_rfcm_forfft = time domain waveforms after trigger path v_banding_rfcm = amplitude in Fourier domain after trigger path

Parameters

settings1	:: Settings - the simulation settings
anita1	:: Anita - anita object
ant	:: int - antenna number

◆ WhichBandsPass()

void ChanTrigger::WhichBandsPass	(	Settings *	settings1,
		Anita *	anita1,
		GlobalTrigger *	globaltrig1,
		Balloon *	bn1,
		int	ilayer,
		int	ifold,
		double	dangle,
		double	emfrac,
		double	hadfrac,
		double	thresholds[2][5]
	)

Which bands passes the trigger.

Parameters

settings1	:: Settings - simulation settings
anita1	:: Anita - anita object
globaltrig1	:: GlobalTrigger - global trigger object
bn1	:: Balloon - balloon object
ilayer	:: int - layer number
ifold	:: int - phi sector
dangle	:: int -
emfrac	:: double - electromagnetic fraction of the shower
hadfrac	:: double - hadronic fraction of the shower
thresholds	:: double [2][5] - relative power thresholds for each pol and band

◆ WhichBandsPassTrigger1()

void ChanTrigger::WhichBandsPassTrigger1	(	Settings *	settings1,
		Anita *	anita1,
		GlobalTrigger *	globaltrig1,
		Balloon *	bn1,
		int	ilayer,
		int	ifold,
		double	thresholds[2][5]
	)

Which bands passes the trigger (for trigger scheme 0 and 1)

Parameters

settings1	:: Settings - simulation settings
anita1	:: Anita - anita object
globaltrig1	:: GlobalTrigger - global trigger object
bn1	:: Balloon - balloon object
ilayer	:: int - layer number
ifold	:: int - phi sector
thresholds	:: double [2][5] - relative power thresholds for each pol and band

◆ WhichBandsPassTrigger2()

void ChanTrigger::WhichBandsPassTrigger2	(	Settings *	settings1,
		Anita *	anita1,
		GlobalTrigger *	globaltrig1,
		Balloon *	bn1,
		int	ilayer,
		int	ifold,
		double	dangle,
		double	emfrac,
		double	hadfrac,
		double	thresholds[2][5]
	)

Which bands passes the trigger (for trigger scheme larger than 2)

Parameters

settings1	:: Settings - simulation settings
anita1	:: Anita - anita object
globaltrig1	:: GlobalTrigger - global trigger object
bn1	:: Balloon - balloon object
ilayer	:: int - layer number
ifold	:: int - phi sector
dangle	:: int -
emfrac	:: double - electromagnetic fraction of the shower
hadfrac	:: double - hadronic fraction of the shower
thresholds	:: double [2][5] - relative power thresholds for each pol and band

Member Data Documentation

◆ bwslice_energy_pol0

double ChanTrigger::bwslice_energy_pol0[5]

Square the sum of voltage for each slice in bandwidth for the 0th polarization.

◆ bwslice_energy_pol1

double ChanTrigger::bwslice_energy_pol1[5]

Square the sum of voltage for each slice in bandwidth for the 1st polarization.

◆ bwslice_energy_pole

double ChanTrigger::bwslice_energy_pole[5]

Square the sum of voltage for each slice in bandwidth for e polarization. The 5th element is the full band.

◆ bwslice_energy_polh

double ChanTrigger::bwslice_energy_polh[5]

Square the sum of voltage for each slice in bandwidth for h polarization. The 5th element is the full band.

◆ bwslice_volts_pol0

double ChanTrigger::bwslice_volts_pol0[5]

Sum voltage for each slice in bandwidth for the lcp polarization.

◆ bwslice_volts_pol0_em

double ChanTrigger::bwslice_volts_pol0_em[5]

Component of the voltage that comes from the em shower for 0th polarization.

◆ bwslice_volts_pol1

double ChanTrigger::bwslice_volts_pol1[5]

Sum voltage for each slice in bandwidth for the rcp polarization.

◆ bwslice_volts_pol1_em

double ChanTrigger::bwslice_volts_pol1_em[5]

Component of the voltage that comes from the em shower for 1st polarization.

◆ bwslice_volts_pole

double ChanTrigger::bwslice_volts_pole[5]

Sum voltage for each slice in bandwidth for the e polarization.

◆ bwslice_volts_polh

double ChanTrigger::bwslice_volts_polh[5]

Sum voltage for each slice in bandwidth for the h polarization.

◆ cw_digPath

double ChanTrigger::cw_digPath[2][Anita::HALFNFOUR]

For digitizer path, time domain cw.

◆ e_component

double ChanTrigger::e_component

private

E comp along polarization.

◆ e_component_kvector

double ChanTrigger::e_component_kvector

private

component of e-field along the rx e-plane

◆ flag_e

vector<int> ChanTrigger::flag_e[5]

Which bands pass trigger e.

◆ flag_h

vector<int> ChanTrigger::flag_h[5]

Which bands pass trigger h.

◆ h_component

double ChanTrigger::h_component

private

H comp along polarization.

◆ h_component_kvector

double ChanTrigger::h_component_kvector

private

component of the e-field along the rx h-plane

◆ HALFNFOUR

const unsigned ChanTrigger::HALFNFOUR = 512

staticprivate

Half of the number of points in the Fourier space.

◆ hitangle_e

double ChanTrigger::hitangle_e

private

angle the ray hits the antenna wrt e-plane

◆ hitangle_h

double ChanTrigger::hitangle_h

private

angle the ray hits the antenna wrt h-plane

◆ integral_vmmhz

double ChanTrigger::integral_vmmhz

Electric field integral.

◆ justNoise_digPath

double ChanTrigger::justNoise_digPath[2][Anita::HALFNFOUR]

For digitizer path, time domain noise from flight.

◆ justNoise_trigPath

double ChanTrigger::justNoise_trigPath[2][Anita::HALFNFOUR]

For trigger path, time domain noise from flight.

◆ justSig_digPath

double ChanTrigger::justSig_digPath[2][Anita::HALFNFOUR]

Just signal in trigger path.

◆ justSig_trigPath

double ChanTrigger::justSig_trigPath[2][Anita::HALFNFOUR]

Just signal in trigger path.

◆ n_component

double ChanTrigger::n_component

private

normal comp along polarization

◆ n_component_kvector

double ChanTrigger::n_component_kvector

private

component of the e-field along the normal

◆ NCHANNELS

const int ChanTrigger::NCHANNELS =32

staticprivate

Number of channells on each surf.

◆ NFOUR

const unsigned ChanTrigger::NFOUR = 1024

staticprivate

Number of points in Fourier space.

◆ noise_eachband

double ChanTrigger::noise_eachband[2][Anita::NBANDS_MAX]

Noise in each band.

◆ NPOINTS

const int ChanTrigger::NPOINTS =4073

staticprivate

Max number of points from surf measurements.

◆ NSURF

const int ChanTrigger::NSURF =9

staticprivate

Number of surfs.

◆ NSURFMINUSONE

const int ChanTrigger::NSURFMINUSONE =8

staticprivate

Number of surfs minus on1.

◆ NSURFPLUSONE

const int ChanTrigger::NSURFPLUSONE =10

staticprivate

Number of surfs plus one.

◆ passes_eachband

int ChanTrigger::passes_eachband[2][Anita::NBANDS_MAX]

Whether the signal passes or not each band.

◆ signal_eachband

double ChanTrigger::signal_eachband[2][Anita::NBANDS_MAX]

Signal in each band.

◆ thispowerthresh

double ChanTrigger::thispowerthresh

private

Relative power threshold.

◆ thisrate

double ChanTrigger::thisrate

private

Rate in MHz.

◆ threshold_eachband

double ChanTrigger::threshold_eachband[2][Anita::NBANDS_MAX]

Threshold in each band.

◆ unwarned

int ChanTrigger::unwarned

Whether we have warned the user about resetting thresholds when they are beyond the measured bounds.

◆ v_banding_rfcm

double ChanTrigger::v_banding_rfcm[2][5][Anita::NFREQ]

This is Volts/m as a function of frequency after rfcm's and banding.

◆ v_banding_rfcm_forfft

double ChanTrigger::v_banding_rfcm_forfft[2][5][HALFNFOUR]

Starts out as V/s vs. freq after banding, rfcm, after fft it is V vs. t.

◆ v_banding_rfcm_forfft_temp

double ChanTrigger::v_banding_rfcm_forfft_temp[2][5][HALFNFOUR]

Use for the averaging over 10 neighboring bins.

◆ vhz_rx

double ChanTrigger::vhz_rx[2][5][Anita::NFREQ]

Array of amplitudes in the Fourier domain (V/Hz) after the antenna gain. Indeces stand for [ipol][iband][ifreq].

◆ vm_banding_rfcm_forfft

double ChanTrigger::vm_banding_rfcm_forfft[2][5][HALFNFOUR]

Tunnel diode input (signal + noise)

◆ vm_banding_rfcm_forfft_justNoise

double ChanTrigger::vm_banding_rfcm_forfft_justNoise[2][5][HALFNFOUR]

Tunnel diode input (just noise)

◆ vnoise_eachband

vector<double> ChanTrigger::vnoise_eachband[2]

Noise in each band.

◆ volts_rx_forfft

double ChanTrigger::volts_rx_forfft[2][5][Anita::HALFNFOUR]

Array of time domain after the antenna gain. Indeces stand for [ipol][iband][itime].

◆ volts_rx_rfcm

double ChanTrigger::volts_rx_rfcm[2][Anita::HALFNFOUR]

For digitizer path, time domain voltage vs. time after rx, rfcm's.

◆ volts_rx_rfcm_lab

double ChanTrigger::volts_rx_rfcm_lab[2][Anita::HALFNFOUR]

For digitizer path, time domain voltage vs. time after rx, rfcm's and lab.

◆ volts_rx_rfcm_lab_all

double ChanTrigger::volts_rx_rfcm_lab_all[2][48][Anita::HALFNFOUR]

For digitizer path, time domain voltage vs. time after rx, rfcm's and lab.

◆ vpasses_eachband

vector<int> ChanTrigger::vpasses_eachband[2]

Whether the signal passes or not each band.

◆ vsignal_eachband

vector<double> ChanTrigger::vsignal_eachband[2]

Signal in each band.

◆ vthreshold_eachband

vector<double> ChanTrigger::vthreshold_eachband[2]

Threshold in each band.

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

Public Member Functions

Static Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ChanTrigger()

Member Function Documentation

◆ ADCCountstoPowerThreshold()

◆ addToChannelSums()

◆ ApplyAntennaGain()

◆ applyButterworthFilter()

◆ applyImpulseResponseDigitizer()

◆ applyImpulseResponseTrigger()

◆ calculateCW()

◆ ConvertEHtoLREfield()

◆ ConvertEHtoLREnergy()

◆ DigitizerPath()

◆ DiodeConvolution()

◆ FindPeak()

◆ GetNoise()

◆ getNoiseFromFlight()

◆ getPulserAtAMPA()

◆ getRate()

◆ GetThresholds()

◆ InitializeEachBand()

◆ injectImpulseAfterAntenna()

◆ injectImpulseAtSurf()

◆ IsItUnmasked()

◆ L1Trigger()

◆ rateToThreshold()

◆ saveDigitizerWaveforms()

◆ saveTriggerWaveforms()

◆ TimeShiftAndSignalFluct()

◆ TriggerPath()

◆ WhichBandsPass()

◆ WhichBandsPassTrigger1()

◆ WhichBandsPassTrigger2()

Member Data Documentation

◆ bwslice_energy_pol0

◆ bwslice_energy_pol1

◆ bwslice_energy_pole

◆ bwslice_energy_polh

◆ bwslice_volts_pol0

◆ bwslice_volts_pol0_em

◆ bwslice_volts_pol1

◆ bwslice_volts_pol1_em

◆ bwslice_volts_pole

◆ bwslice_volts_polh

◆ cw_digPath

◆ e_component

◆ e_component_kvector

◆ flag_e

◆ flag_h

◆ h_component

◆ h_component_kvector

◆ HALFNFOUR

◆ hitangle_e

◆ hitangle_h

◆ integral_vmmhz

◆ justNoise_digPath

◆ justNoise_trigPath

◆ justSig_digPath

◆ justSig_trigPath

◆ n_component

◆ n_component_kvector

◆ NCHANNELS

◆ NFOUR

◆ noise_eachband

◆ NPOINTS

◆ NSURF

◆ NSURFMINUSONE

◆ NSURFPLUSONE

◆ passes_eachband

◆ signal_eachband

◆ thispowerthresh

◆ thisrate

◆ threshold_eachband