icemc
Functions | Variables
testThermalNoise.cc File Reference
#include <iostream>
#include <fstream>
#include <sstream>
#include <math.h>
#include <ctype.h>
#include <string>
#include <stdio.h>
#include <stdlib.h>
#include <vector>
#include <array>
#include <time.h>
#include "TTreeIndex.h"
#include "TChain.h"
#include "TH1.h"
#include "TF1.h"
#include "TF2.h"
#include "TFile.h"
#include "icemc_random.h"
#include "TTree.h"
#include "TLegend.h"
#include "TLine.h"
#include "TROOT.h"
#include "TPostScript.h"
#include "TCanvas.h"
#include "TH2F.h"
#include "TText.h"
#include "TProfile.h"
#include "TGraphErrors.h"
#include "TGraphAsymmErrors.h"
#include "TStyle.h"
#include "TMath.h"
#include <unistd.h>
#include "TVector3.h"
#include "TRotation.h"
#include "TSpline.h"
#include "Math/InterpolationTypes.h"
#include "Math/Interpolator.h"
#include "signal.h"
#include "Constants.h"
#include "Settings.h"
#include "position.hh"
#include "earthmodel.hh"
#include "Tools.h"
#include "vector.hh"
#include "roughness.hh"
#include "anita.hh"
#include "balloon.hh"
#include "icemodel.hh"
#include "Spectra.h"
#include "signal.hh"
#include "secondaries.hh"
#include "ray.hh"
#include "counting.hh"
#include "Primaries.h"
#include "Taumodel.hh"
#include "screen.hh"
#include "GlobalTrigger.h"
#include "ChanTrigger.h"
#include "SimulatedSignal.h"
#include "EnvironmentVariable.h"
#include <typeinfo>

Functions

 ClassImp (RX)
 
double GetAirDistance (double altitude_bn, double beta)
 
void SetupViewangles (Signal *sig1)
 
void GetAir (double *col1)
 
double GetThisAirColumn (Settings *, Position r_in, Vector nnu, Position posnu, double *col1, double &cosalpha, double &mytheta, double &cosbeta0, double &mybeta)
 
double ScaleVmMHz (double vmmhz1m_max, const Position &posnu, const Position &r_bn)
 
double IsItDoubleBang (double exitlength, double plepton)
 
int WhereIsSecondBang (const Position &posnu, const Vector &nnu, double nuexitlength, double pnu, IceModel *antarctica1, const Position &r_bn, Position &posnu2, Position &rfexit_db, Vector &n_exit2bn_db)
 
void GetCurrent (string &current)
 
double GetAverageVoltageFromAntennasHit (Settings *settings1, int *nchannels_perrx_triggered, double *voltagearray, double &volts_rx_sum)
 
Vector GetPolarization (const Vector &nnu, const Vector &nrf2_iceside)
 
void Attenuate (IceModel *antartica1, Settings *settings1, double &vmmhz_max, double rflength, const Position &posnu)
 
void Attenuate_down (IceModel *antarctica1, Settings *settings1, double &vmmhz_max, const Position &rfexit2, const Position &posnu, const Position &posnu_down)
 
void IsAbsorbed (double chord_kgm2, double len_int_kgm2, double &weight)
 
int GetRayIceSide (const Vector &n_exit2rx, const Vector &nsurf_rfexit, double nexit, double nenter, Vector &nrf2_iceside)
 
double GetViewAngle (const Vector &nrf2_iceside, const Vector &nnu)
 
int TIR (const Vector &n_surf, const Vector &nrf2_iceside, double N_IN, double N_OUT)
 
void IntegrateBands (Anita *anita1, int k, Screen *panel1, double *freq, double scalefactor, double *sumsignal)
 
void Integrate (Anita *anita1, int j, int k, double *vmmhz, double *freq, double scalefactor, double sumsignal)
 
void interrupt_signal_handler (int)
 
void WriteNeutrinoInfo (Position &, Vector &, Position &, double, string, string, double, ofstream &nu_out)
 
void CloseTFile (TFile *hfile)
 
int Getmine (double *, double *, double *, double *)
 
void Getearth (double *, double *, double *, double *)
 
int main (int argc, char **argv)
 

Variables

TaumodelTauPtr = NULL
 
const string ICEMC_SRC_DIR = EnvironmentVariable::ICEMC_SRC_DIR()
 
const int NVIEWANGLE =100
 
int inu
 
UInt_t eventNumber
 
double eventsfound_beforetrigger =0.
 
double eventsfound_crust =0
 
double eventsfound_weightgt01 =0
 
double eventsfound_belowhorizon =0
 
double eventsfound =0
 
double eventsfound_prob =0
 
double sum [3]
 
double poissonerror_minus [21] = {0.-0.00, 1.-0.37, 2.-0.74, 3.-1.10, 4.-2.34, 5.-2.75, 6.-3.82, 7.-4.25, 8.-5.30, 9.-6.33, 10.-6.78, 11.-7.81, 12.-8.83, 13.-9.28, 14.-10.30, 15.-11.32, 16.-12.33, 17.-12.79, 18.-13.81, 19.-14.82, 20.-15.83}
 
double poissonerror_plus [21] = {1.29-0., 2.75-1., 4.25-2., 5.30-3., 6.78-4., 7.81-5., 9.28-6., 10.30-7., 11.32-8., 12.79-9., 13.81-10., 14.82-11., 16.29-12., 17.30-13., 18.32-14., 19.32-15., 20.80-16., 21.81-17., 22.82-18., 23.82-19., 25.30-20.}
 
const int NBINS =10
 
double MIN_LOGWEIGHT =-3
 
double MAX_LOGWEIGHT =-1
 
int index_weights =0
 
double eventsfound_binned [NBINS]
 
double eventsfound_binned_e [NBINS]
 
double eventsfound_binned_mu [NBINS]
 
double eventsfound_binned_tau [NBINS]
 
double km3sr = 0
 
double km3sr_e =0
 
double km3sr_mu =0
 
double km3sr_tau =0
 
double error_plus =0
 
double error_e_plus =0
 
double error_mu_plus =0
 
double error_tau_plus =0
 
double error_minus =0
 
double error_e_minus =0
 
double error_mu_minus =0
 
double error_tau_minus =0
 
int ierr =0
 
double gain_dipole =2.15
 
double changle_deg =0
 
int NNU
 
int whichray =0
 
double RANDOMISEPOL =0.
 
double volume_thishorizon
 
int realtime_this
 
double longitude_this
 
double latitude_this
 
double altitude_this
 
double heading_this =0.
 
double gps_offset =0
 
double pnu =pow(10., 20)
 energy of neutrinos More...
 
double MEANX =0
 
double MEANY =0.
 
double SIGNALRADIUS =2.
 
double bwslice_vnoise_thislayer [4]
 
int passes_thisevent =0
 
int unmasked_thisevent =0
 
int discones_passing
 
int NDISCONES =8
 
double heff_discone =0
 
double polarfactor_discone =0.
 
double thislambda =0
 
double volts_discone =0.
 
double vnoise_discone =0.
 
double BW_DISCONES =300.E6-120.E6
 
double fresnel1 =0
 net fresnel factor on field at ice-firn interface More...
 
double fresnel1_eachboresight [Anita::NLAYERS_MAX][Anita::NPHI_MAX]
 
double fresnel2 =0
 net fresnel factor on field at firn-air More...
 
double mag1 =0
 
double mag1_eachboresight [Anita::NLAYERS_MAX][Anita::NPHI_MAX]
 
double mag2 =0
 
double t_coeff_pokey
 
double t_coeff_slappy
 
double rflength =0
 
double e_comp_max1 =0
 
double h_comp_max1 =0
 
double e_comp_max2 =0
 
double h_comp_max2 =0
 
double e_comp_max3 =0
 
double h_comp_max3 =0
 
double diffexit =0
 
double diff_3tries =0
 
double diffnorm =0
 
double diffrefr =0
 
double costheta_inc =0
 
double costheta_exit =0
 
double theta_rf_atbn
 
double theta_rf_atbn_measured
 
double costhetanu =-1000
 
double theta_in =0
 
double lat_in =0
 
double nearthlayers =0
 
double weight_prob =0.
 
double weight1 =0
 
double weight =0.
 
double logweight =0.
 
double len_int =0
 
double pieceofkm2sr =0
 
double CUTONWEIGHTS =1.E-10
 
int count_inthisloop1 =0
 
int count_inthisloop2 =0
 
int count_inthisloop3 =0
 
double averaging_thetas1 =0
 
double averaging_thetas2 =0
 
double averaging_thetas3 =0
 
int count_total =0
 
int count_asktrigger =0
 
int count_asktrigger_nfb =0
 
int count_pass =0
 
double count_passestrigger_w =0
 
int passestrigger =0
 
int allcuts [2] ={0, 0}
 
double allcuts_weighted [2] ={0, 0}
 
double allcuts_weighted_polarization [3] ={0, 0, 0}
 
int count_chanceofsurviving =0
 
int count_chanceinhell0 =0
 
double count_chanceinhell2_w =0
 
int chanceinhell2 =0
 
int count_chordgoodlength =0
 
int count_d2goodlength =0
 
int count_rx =0
 
double sum_frac [3]
 
double sum_frac_db [3]
 
const int NBINS_DISTANCE =28
 
double eventsfound_binned_distance [NBINS_DISTANCE] = {0.}
 
int index_distance =0
 
double km3sr_distance [NBINS_DISTANCE] = {0.}
 
double error_distance_plus [NBINS_DISTANCE] = {0.}
 
double error_distance_minus [NBINS_DISTANCE] = {0.}
 
int eventsfound_binned_distance_forerror [NBINS_DISTANCE][NBINS] = {{0}}
 
double km3sr_db = 0
 
double km3sr_nfb =0
 
double ptau =0
 
int count_passestrigger_nfb =0
 
double percent_increase_db =0
 
double percent_increase_nfb =0
 
double percent_increase_total =0
 
double error_nfb =0
 
double error_km3sr_nfb =0
 
double error_percent_increase_nfb =0
 
Vector n_exit2bn_db [5]
 
Vector nrf_iceside_db [5]
 
double n_exit_phi
 
int count_dbexitsice =0
 
double eventsfound_nfb_binned [NBINS]
 
double heff_max =0.62639
 
double scalefactor_distance =0
 
double scalefactor_attenuation =0
 
double MAX_ATTENLENGTH =1671
 
double maxtaper =0
 
double dviewangle_deg =0
 deviation from the cherenkov angle More...
 
double forseckel [NVIEWANGLE][Anita::NFREQ]
 
double viewangles [NVIEWANGLE]
 
int err =0
 
double djunk
 
int max_antenna0 = -1
 
int max_antenna1 = -1
 
int max_antenna2 = -1
 
double max_antenna_volts0 = 0
 
double max_antenna_volts0_em = 0
 
double max_antenna_volts1 = 0
 
double max_antenna_volts2 = 0
 
double rx0_signal_eachband [2][5]
 
double rx0_threshold_eachband [2][5]
 
double rx0_noise_eachband [2][5]
 
int rx0_passes_eachband [2][5]
 
double voltagearray [Anita::NLAYERS_MAX *Anita::NPHI_MAX]
 
Vector ant_max_normal0
 
Vector ant_max_normal1
 
Vector ant_max_normal2
 
double vmmhz1m_visible = 0
 
int freq_bins = Anita::NFREQ
 
double total_kgm2 = 0
 
double nnu_array [3]
 
double r_in_array [3]
 
double nsurf_rfexit_array [3]
 
double nsurf_rfexit_db_array [3]
 
double r_bn_array [3]
 
double n_bn_array [3]
 
double posnu_array [3]
 
double nrf_iceside_array [5][3]
 
double nrf_iceside_db_array [5][3]
 
double ant_max_normal0_array [3]
 
double ant_max_normal1_array [3]
 
double ant_max_normal2_array [3]
 
double n_pol_array [3]
 
double n_exit2bn_array [5][3]
 
double r_enterice_array [3]
 
double n_exit2bn_db_array [5][3]
 
double rfexit_array [5][3]
 
double rfexit_db_array [5][3]
 
int times_crust_entered_det =0
 
int times_mantle_entered_det =0
 
int times_core_entered_det =0
 
int crust_entered =0
 
int mantle_entered =0
 
int core_entered =0
 
int neutrinos_passing_all_cuts =0
 
double sum_weights =0
 
double justNoise_trig [2][48][512]
 
double justSignal_trig [2][48][512]
 
double justNoise_dig [2][48][512]
 
double justSignal_dig [2][48][512]
 
bool ABORT_EARLY = false
 
double thresholdsAnt [48][2][5]
 
double thresholdsAntPass [48][2][5]
 
double threshold_start =-1.
 
double threshold_end =-6.
 
const int NTHRESHOLDS =20
 
double threshold_step =(threshold_end-threshold_start)/(double)NTHRESHOLDS
 
double npass_v_thresh [NTHRESHOLDS] ={0.}
 
double denom_v_thresh [NTHRESHOLDS] ={0.}
 
double npass_h_thresh [NTHRESHOLDS] ={0.}
 
double denom_h_thresh [NTHRESHOLDS] ={0.}
 
double thresholds [NTHRESHOLDS]
 

Function Documentation

◆ Attenuate()

void Attenuate ( IceModel antartica1,
Settings settings1,
double &  vmmhz_max,
double  rflength,
const Position posnu 
)

◆ Attenuate_down()

void Attenuate_down ( IceModel antarctica1,
Settings settings1,
double &  vmmhz_max,
const Position rfexit2,
const Position posnu,
const Position posnu_down 
)

◆ ClassImp()

ClassImp ( RX  )

◆ CloseTFile()

void CloseTFile ( TFile *  hfile)

◆ GetAir()

void GetAir ( double *  col1)

◆ GetAirDistance()

double GetAirDistance ( double  altitude_bn,
double  beta 
)

◆ GetAverageVoltageFromAntennasHit()

double GetAverageVoltageFromAntennasHit ( Settings settings1,
int *  nchannels_perrx_triggered,
double *  voltagearray,
double &  volts_rx_sum 
)

◆ GetCurrent()

void GetCurrent ( string &  current)

◆ Getearth()

void Getearth ( double *  ,
double *  ,
double *  ,
double *   
)

◆ Getmine()

int Getmine ( double *  ,
double *  ,
double *  ,
double *   
)

◆ GetPolarization()

Vector GetPolarization ( const Vector nnu,
const Vector nrf2_iceside 
)

◆ GetRayIceSide()

int GetRayIceSide ( const Vector n_exit2rx,
const Vector nsurf_rfexit,
double  nexit,
double  nenter,
Vector nrf2_iceside 
)

◆ GetThisAirColumn()

double GetThisAirColumn ( Settings ,
Position  r_in,
Vector  nnu,
Position  posnu,
double *  col1,
double &  cosalpha,
double &  mytheta,
double &  cosbeta0,
double &  mybeta 
)

◆ GetViewAngle()

double GetViewAngle ( const Vector nrf2_iceside,
const Vector nnu 
)

◆ Integrate()

void Integrate ( Anita anita1,
int  j,
int  k,
double *  vmmhz,
double *  freq,
double  scalefactor,
double  sumsignal 
)

◆ IntegrateBands()

void IntegrateBands ( Anita anita1,
int  k,
Screen panel1,
double *  freq,
double  scalefactor,
double *  sumsignal 
)

◆ interrupt_signal_handler()

void interrupt_signal_handler ( int  sig)

◆ IsAbsorbed()

void IsAbsorbed ( double  chord_kgm2,
double  len_int_kgm2,
double &  weight 
)

◆ IsItDoubleBang()

double IsItDoubleBang ( double  exitlength,
double  plepton 
)

◆ main()

int main ( int  argc,
char **  argv 
)

◆ ScaleVmMHz()

double ScaleVmMHz ( double  vmmhz1m_max,
const Position posnu,
const Position r_bn 
)

◆ SetupViewangles()

void SetupViewangles ( Signal sig1)

◆ TIR()

int TIR ( const Vector n_surf,
const Vector nrf2_iceside,
double  N_IN,
double  N_OUT 
)

◆ WhereIsSecondBang()

int WhereIsSecondBang ( const Position posnu,
const Vector nnu,
double  nuexitlength,
double  pnu,
IceModel antarctica1,
const Position r_bn,
Position posnu2,
Position rfexit_db,
Vector n_exit2bn_db 
)

◆ WriteNeutrinoInfo()

void WriteNeutrinoInfo ( Position ,
Vector ,
Position ,
double  ,
string  ,
string  ,
double  ,
ofstream &  nu_out 
)

Variable Documentation

◆ ABORT_EARLY

bool ABORT_EARLY = false

◆ allcuts

int allcuts[2] ={0, 0}

◆ allcuts_weighted

double allcuts_weighted[2] ={0, 0}

◆ allcuts_weighted_polarization

double allcuts_weighted_polarization[3] ={0, 0, 0}

◆ altitude_this

double altitude_this

◆ ant_max_normal0

Vector ant_max_normal0

◆ ant_max_normal0_array

double ant_max_normal0_array[3]

◆ ant_max_normal1

Vector ant_max_normal1

◆ ant_max_normal1_array

double ant_max_normal1_array[3]

◆ ant_max_normal2

Vector ant_max_normal2

◆ ant_max_normal2_array

double ant_max_normal2_array[3]

◆ averaging_thetas1

double averaging_thetas1 =0

◆ averaging_thetas2

double averaging_thetas2 =0

◆ averaging_thetas3

double averaging_thetas3 =0

◆ BW_DISCONES

double BW_DISCONES =300.E6-120.E6

◆ bwslice_vnoise_thislayer

double bwslice_vnoise_thislayer[4]

◆ chanceinhell2

int chanceinhell2 =0

◆ changle_deg

double changle_deg =0

◆ core_entered

int core_entered =0

◆ costheta_exit

double costheta_exit =0

◆ costheta_inc

double costheta_inc =0

◆ costhetanu

double costhetanu =-1000

◆ count_asktrigger

int count_asktrigger =0

◆ count_asktrigger_nfb

int count_asktrigger_nfb =0

◆ count_chanceinhell0

int count_chanceinhell0 =0

◆ count_chanceinhell2_w

double count_chanceinhell2_w =0

◆ count_chanceofsurviving

int count_chanceofsurviving =0

◆ count_chordgoodlength

int count_chordgoodlength =0

◆ count_d2goodlength

int count_d2goodlength =0

◆ count_dbexitsice

int count_dbexitsice =0

◆ count_inthisloop1

int count_inthisloop1 =0

◆ count_inthisloop2

int count_inthisloop2 =0

◆ count_inthisloop3

int count_inthisloop3 =0

◆ count_pass

int count_pass =0

◆ count_passestrigger_nfb

int count_passestrigger_nfb =0

◆ count_passestrigger_w

double count_passestrigger_w =0

◆ count_rx

int count_rx =0

◆ count_total

int count_total =0

◆ crust_entered

int crust_entered =0

◆ CUTONWEIGHTS

double CUTONWEIGHTS =1.E-10

◆ denom_h_thresh

double denom_h_thresh[NTHRESHOLDS] ={0.}

◆ denom_v_thresh

double denom_v_thresh[NTHRESHOLDS] ={0.}

◆ diff_3tries

double diff_3tries =0

◆ diffexit

double diffexit =0

◆ diffnorm

double diffnorm =0

◆ diffrefr

double diffrefr =0

◆ discones_passing

int discones_passing

◆ djunk

double djunk

◆ dviewangle_deg

double dviewangle_deg =0

deviation from the cherenkov angle

◆ e_comp_max1

double e_comp_max1 =0

◆ e_comp_max2

double e_comp_max2 =0

◆ e_comp_max3

double e_comp_max3 =0

◆ err

int err =0

◆ error_distance_minus

double error_distance_minus[NBINS_DISTANCE] = {0.}

◆ error_distance_plus

double error_distance_plus[NBINS_DISTANCE] = {0.}

◆ error_e_minus

double error_e_minus =0

◆ error_e_plus

double error_e_plus =0

◆ error_km3sr_nfb

double error_km3sr_nfb =0

◆ error_minus

double error_minus =0

◆ error_mu_minus

double error_mu_minus =0

◆ error_mu_plus

double error_mu_plus =0

◆ error_nfb

double error_nfb =0

◆ error_percent_increase_nfb

double error_percent_increase_nfb =0

◆ error_plus

double error_plus =0

◆ error_tau_minus

double error_tau_minus =0

◆ error_tau_plus

double error_tau_plus =0

◆ eventNumber

UInt_t eventNumber

◆ eventsfound

double eventsfound =0

◆ eventsfound_beforetrigger

double eventsfound_beforetrigger =0.

◆ eventsfound_belowhorizon

double eventsfound_belowhorizon =0

◆ eventsfound_binned

double eventsfound_binned[NBINS]

◆ eventsfound_binned_distance

double eventsfound_binned_distance[NBINS_DISTANCE] = {0.}

◆ eventsfound_binned_distance_forerror

int eventsfound_binned_distance_forerror[NBINS_DISTANCE][NBINS] = {{0}}

◆ eventsfound_binned_e

double eventsfound_binned_e[NBINS]

◆ eventsfound_binned_mu

double eventsfound_binned_mu[NBINS]

◆ eventsfound_binned_tau

double eventsfound_binned_tau[NBINS]

◆ eventsfound_crust

double eventsfound_crust =0

◆ eventsfound_nfb_binned

double eventsfound_nfb_binned[NBINS]

◆ eventsfound_prob

double eventsfound_prob =0

◆ eventsfound_weightgt01

double eventsfound_weightgt01 =0

◆ forseckel

double forseckel[NVIEWANGLE][Anita::NFREQ]

◆ freq_bins

int freq_bins = Anita::NFREQ

◆ fresnel1

double fresnel1 =0

net fresnel factor on field at ice-firn interface

◆ fresnel1_eachboresight

double fresnel1_eachboresight[Anita::NLAYERS_MAX][Anita::NPHI_MAX]

◆ fresnel2

double fresnel2 =0

net fresnel factor on field at firn-air

◆ gain_dipole

double gain_dipole =2.15

◆ gps_offset

double gps_offset =0

◆ h_comp_max1

double h_comp_max1 =0

◆ h_comp_max2

double h_comp_max2 =0

◆ h_comp_max3

double h_comp_max3 =0

◆ heading_this

double heading_this =0.

◆ heff_discone

double heff_discone =0

◆ heff_max

double heff_max =0.62639

◆ ICEMC_SRC_DIR

const string ICEMC_SRC_DIR = EnvironmentVariable::ICEMC_SRC_DIR()

◆ ierr

int ierr =0

◆ index_distance

int index_distance =0

◆ index_weights

int index_weights =0

◆ inu

int inu

◆ justNoise_dig

double justNoise_dig[2][48][512]

◆ justNoise_trig

double justNoise_trig[2][48][512]

◆ justSignal_dig

double justSignal_dig[2][48][512]

◆ justSignal_trig

double justSignal_trig[2][48][512]

◆ km3sr

double km3sr = 0

◆ km3sr_db

double km3sr_db = 0

◆ km3sr_distance

double km3sr_distance[NBINS_DISTANCE] = {0.}

◆ km3sr_e

double km3sr_e =0

◆ km3sr_mu

double km3sr_mu =0

◆ km3sr_nfb

double km3sr_nfb =0

◆ km3sr_tau

double km3sr_tau =0

◆ lat_in

double lat_in =0

◆ latitude_this

double latitude_this

◆ len_int

double len_int =0

◆ logweight

double logweight =0.

◆ longitude_this

double longitude_this

◆ mag1

double mag1 =0

◆ mag1_eachboresight

double mag1_eachboresight[Anita::NLAYERS_MAX][Anita::NPHI_MAX]

◆ mag2

double mag2 =0

◆ mantle_entered

int mantle_entered =0

◆ max_antenna0

int max_antenna0 = -1

◆ max_antenna1

int max_antenna1 = -1

◆ max_antenna2

int max_antenna2 = -1

◆ max_antenna_volts0

double max_antenna_volts0 = 0

◆ max_antenna_volts0_em

double max_antenna_volts0_em = 0

◆ max_antenna_volts1

double max_antenna_volts1 = 0

◆ max_antenna_volts2

double max_antenna_volts2 = 0

◆ MAX_ATTENLENGTH

double MAX_ATTENLENGTH =1671

◆ MAX_LOGWEIGHT

double MAX_LOGWEIGHT =-1

◆ maxtaper

double maxtaper =0

◆ MEANX

double MEANX =0

◆ MEANY

double MEANY =0.

◆ MIN_LOGWEIGHT

double MIN_LOGWEIGHT =-3

◆ n_bn_array

double n_bn_array[3]

◆ n_exit2bn_array

double n_exit2bn_array[5][3]

◆ n_exit2bn_db

Vector n_exit2bn_db[5]

◆ n_exit2bn_db_array

double n_exit2bn_db_array[5][3]

◆ n_exit_phi

double n_exit_phi

◆ n_pol_array

double n_pol_array[3]

◆ NBINS

const int NBINS =10

◆ NBINS_DISTANCE

const int NBINS_DISTANCE =28

◆ NDISCONES

int NDISCONES =8

◆ nearthlayers

double nearthlayers =0

◆ neutrinos_passing_all_cuts

int neutrinos_passing_all_cuts =0

◆ NNU

int NNU

◆ nnu_array

double nnu_array[3]

◆ npass_h_thresh

double npass_h_thresh[NTHRESHOLDS] ={0.}

◆ npass_v_thresh

double npass_v_thresh[NTHRESHOLDS] ={0.}

◆ nrf_iceside_array

double nrf_iceside_array[5][3]

◆ nrf_iceside_db

Vector nrf_iceside_db[5]

◆ nrf_iceside_db_array

double nrf_iceside_db_array[5][3]

◆ nsurf_rfexit_array

double nsurf_rfexit_array[3]

◆ nsurf_rfexit_db_array

double nsurf_rfexit_db_array[3]

◆ NTHRESHOLDS

const int NTHRESHOLDS =20

◆ NVIEWANGLE

const int NVIEWANGLE =100

◆ passes_thisevent

int passes_thisevent =0

◆ passestrigger

int passestrigger =0

◆ percent_increase_db

double percent_increase_db =0

◆ percent_increase_nfb

double percent_increase_nfb =0

◆ percent_increase_total

double percent_increase_total =0

◆ pieceofkm2sr

double pieceofkm2sr =0

◆ pnu

double pnu =pow(10., 20)

energy of neutrinos

◆ poissonerror_minus

double poissonerror_minus[21] = {0.-0.00, 1.-0.37, 2.-0.74, 3.-1.10, 4.-2.34, 5.-2.75, 6.-3.82, 7.-4.25, 8.-5.30, 9.-6.33, 10.-6.78, 11.-7.81, 12.-8.83, 13.-9.28, 14.-10.30, 15.-11.32, 16.-12.33, 17.-12.79, 18.-13.81, 19.-14.82, 20.-15.83}

◆ poissonerror_plus

double poissonerror_plus[21] = {1.29-0., 2.75-1., 4.25-2., 5.30-3., 6.78-4., 7.81-5., 9.28-6., 10.30-7., 11.32-8., 12.79-9., 13.81-10., 14.82-11., 16.29-12., 17.30-13., 18.32-14., 19.32-15., 20.80-16., 21.81-17., 22.82-18., 23.82-19., 25.30-20.}

◆ polarfactor_discone

double polarfactor_discone =0.

◆ posnu_array

double posnu_array[3]

◆ ptau

double ptau =0

◆ r_bn_array

double r_bn_array[3]

◆ r_enterice_array

double r_enterice_array[3]

◆ r_in_array

double r_in_array[3]

◆ RANDOMISEPOL

double RANDOMISEPOL =0.

◆ realtime_this

int realtime_this

◆ rfexit_array

double rfexit_array[5][3]

◆ rfexit_db_array

double rfexit_db_array[5][3]

◆ rflength

double rflength =0

◆ rx0_noise_eachband

double rx0_noise_eachband[2][5]

◆ rx0_passes_eachband

int rx0_passes_eachband[2][5]

◆ rx0_signal_eachband

double rx0_signal_eachband[2][5]

◆ rx0_threshold_eachband

double rx0_threshold_eachband[2][5]

◆ scalefactor_attenuation

double scalefactor_attenuation =0

◆ scalefactor_distance

double scalefactor_distance =0

◆ SIGNALRADIUS

double SIGNALRADIUS =2.

◆ sum

double sum[3]

◆ sum_frac

double sum_frac[3]

◆ sum_frac_db

double sum_frac_db[3]

◆ sum_weights

double sum_weights =0

◆ t_coeff_pokey

double t_coeff_pokey

◆ t_coeff_slappy

double t_coeff_slappy

◆ TauPtr

Taumodel* TauPtr = NULL

◆ theta_in

double theta_in =0

◆ theta_rf_atbn

double theta_rf_atbn

◆ theta_rf_atbn_measured

double theta_rf_atbn_measured

◆ thislambda

double thislambda =0

◆ threshold_end

double threshold_end =-6.

◆ threshold_start

double threshold_start =-1.

◆ threshold_step

double threshold_step =(threshold_end-threshold_start)/(double)NTHRESHOLDS

◆ thresholds

double thresholds[NTHRESHOLDS]

◆ thresholdsAnt

double thresholdsAnt[48][2][5]

◆ thresholdsAntPass

double thresholdsAntPass[48][2][5]

◆ times_core_entered_det

int times_core_entered_det =0

◆ times_crust_entered_det

int times_crust_entered_det =0

◆ times_mantle_entered_det

int times_mantle_entered_det =0

◆ total_kgm2

double total_kgm2 = 0

◆ unmasked_thisevent

int unmasked_thisevent =0

◆ viewangles

double viewangles[NVIEWANGLE]

◆ vmmhz1m_visible

double vmmhz1m_visible = 0

◆ vnoise_discone

double vnoise_discone =0.

◆ voltagearray

double voltagearray[Anita::NLAYERS_MAX *Anita::NPHI_MAX]

◆ volts_discone

double volts_discone =0.

◆ volume_thishorizon

double volume_thishorizon

◆ weight

double weight =0.

◆ weight1

double weight1 =0

◆ weight_prob

double weight_prob =0.

◆ whichray

int whichray =0