icemc
Macros | Functions | Variables
icemc.cc File Reference
#include <iostream>
#include <fstream>
#include <sstream>
#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 "TH2.h"
#include "TF1.h"
#include "TF2.h"
#include "TFile.h"
#include "TObject.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 "TGraph.h"
#include "TGraphErrors.h"
#include "TGraphAsymmErrors.h"
#include "TStyle.h"
#include "TMath.h"
#include "TMarker.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 <cmath>
#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 "source.hh"
#include <typeinfo>
#include <fenv.h>

Macros

#define DO_SKIP   continue ;
 

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 &posnu1, const Position &r_bn, const Position &rfexit)
 
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)
 
void GetBalloonLocation (Interaction *interaction1, Ray *ray1, Balloon *bn1, IceModel *antarctica)
 
int GetRayIceSide (const Vector &n_exit2rx, const Vector &nsurf_rfexit, double nexit, double nenter, Vector &nrf2_iceside)
 
int GetDirection (Settings *settings1, Interaction *interaction1, const Vector &refr, double deltheta_em, double deltheta_had, double emfrac, double hadfrac, double vmmhz1m_max, double r_fromballoon, Ray *ray1, Signal *sig1, Position posnu, Anita *anita1, Balloon *bn1, Vector &nnu, double &costhetanu, double &theta_threshold)
 
void GetFresnel (Roughness *rough1, int ROUGHNESS_SETTING, const Vector &nsurf_rfexit, const Vector &n_exit2rx, Vector &n_pol, const Vector &nrf2_iceside, double efield, double emfrac, double hadfrac, double deltheta_em, double deltheta_had, double &t_coeff_pokey, double &t_coeff_slappy, double &fresnel, double &mag)
 
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 Summarize (Settings *settings1, Anita *anita1, Counting *count1, Spectra *spectra1, Signal *sig1, Primaries *primary1, double, double eventsfound, double, double, double, double *, double, double, double &, double &, double &, double &, ofstream &, ofstream &, TString, Balloon *, SourceModel *src_model)
 
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)
 
void GetSmearedIncidentAngle (Vector &specular, Vector &nrf_iceside, Vector &n_exit2bn, double SMEARINCIDENTANGLE)
 
TStyle * RootStyle ()
 

Variables

Taumodel * TauPtr = NULL
 
const string ICEMC_SRC_DIR = EnvironmentVariable::ICEMC_SRC_DIR()
 
const int NVIEWANGLE =100
 
int inu
 
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 sum_prob [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.
 
std::vector< int > NNU_per_source
 
std::vector< double > eventsfound_per_source
 
std::vector< double > eventsfound_prob_per_source
 
double volume_thishorizon
 
double 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.
 
const double FREQ_LOW_DISCONES =120.E6
 
const double FREQ_HIGH_DISCONES =1000.E6
 
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 theta_nnu_atbn
 
double theta_nnu_rf_diff_atbn
 
double phi_rf_atbn
 
double phi_nnu_atbn
 
double phi_nnu_rf_diff_atbn
 
Vector nnu_xy
 
Vector r_bn_xy
 
Vector rf_xy
 
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
 
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 = 1
 
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
 
int xsecParam_nutype = 0
 
int xsecParam_nuint = 1
 
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]
 

Macro Definition Documentation

◆ DO_SKIP

#define DO_SKIP   continue ;

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 
)

◆ GetBalloonLocation()

void GetBalloonLocation ( Interaction *  interaction1,
Ray *  ray1,
Balloon *  bn1,
IceModel *  antarctica 
)

◆ GetCurrent()

void GetCurrent ( string &  current)

◆ GetDirection()

int GetDirection ( Settings *  settings1,
Interaction *  interaction1,
const Vector &  refr,
double  deltheta_em,
double  deltheta_had,
double  emfrac,
double  hadfrac,
double  vmmhz1m_max,
double  r_fromballoon,
Ray *  ray1,
Signal *  sig1,
Position  posnu,
Anita *  anita1,
Balloon *  bn1,
Vector &  nnu,
double &  costhetanu,
double &  theta_threshold 
)

◆ Getearth()

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

◆ GetFresnel()

void GetFresnel ( Roughness *  rough1,
int  ROUGHNESS_SETTING,
const Vector &  nsurf_rfexit,
const Vector &  n_exit2rx,
Vector &  n_pol,
const Vector &  nrf2_iceside,
double  efield,
double  emfrac,
double  hadfrac,
double  deltheta_em,
double  deltheta_had,
double &  t_coeff_pokey,
double &  t_coeff_slappy,
double &  fresnel,
double &  mag 
)

◆ 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 
)

◆ GetSmearedIncidentAngle()

void GetSmearedIncidentAngle ( Vector &  specular,
Vector &  nrf_iceside,
Vector &  n_exit2bn,
double  SMEARINCIDENTANGLE 
)

◆ GetThisAirColumn()

double GetThisAirColumn ( Settings *  settings1,
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 
)

< alpha minus 90 degrees

< beta minus 90 degrees

< tau trigger sets the chance to create tau particle

◆ RootStyle()

TStyle* RootStyle ( )

◆ ScaleVmMHz()

double ScaleVmMHz ( double  vmmhz1m_max,
const Position &  posnu1,
const Position &  r_bn,
const Position &  rfexit 
)

◆ SetupViewangles()

void SetupViewangles ( Signal *  sig1)

◆ Summarize()

void Summarize ( Settings *  settings1,
Anita *  anita1,
Counting *  count1,
Spectra *  spectra1,
Signal *  sig1,
Primaries *  primary1,
double  pnu,
double  eventsfound,
double  eventsfound_db,
double  eventsfound_nfb,
double  sigma,
double *  sum,
double  volume,
double  ice_area,
double &  km3sr,
double &  km3sr_e,
double &  km3sr_mu,
double &  km3sr_tau,
ofstream &  foutput,
ofstream &  distanceout,
TString  outputdir,
Balloon *  bn1,
SourceModel *  src_model 
)

◆ 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 &  posnu,
Vector &  nnu,
Position &  r_bn,
double  altitude_int,
string  nuflavor,
string  current,
double  elast_y,
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

◆ 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 = 1

◆ 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

◆ 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_per_source

std::vector<double> eventsfound_per_source

◆ eventsfound_prob

double eventsfound_prob =0

◆ eventsfound_prob_per_source

std::vector<double> eventsfound_prob_per_source

◆ eventsfound_weightgt01

double eventsfound_weightgt01 =0

◆ forseckel

double forseckel[NVIEWANGLE][Anita::NFREQ]

◆ freq_bins

int freq_bins = Anita::NFREQ

◆ FREQ_HIGH_DISCONES

const double FREQ_HIGH_DISCONES =1000.E6

◆ FREQ_LOW_DISCONES

const double FREQ_LOW_DISCONES =120.E6

◆ 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]

◆ NNU_per_source

std::vector<int> NNU_per_source

◆ nnu_xy

Vector nnu_xy

◆ 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

◆ phi_nnu_atbn

double phi_nnu_atbn

◆ phi_nnu_rf_diff_atbn

double phi_nnu_rf_diff_atbn

◆ phi_rf_atbn

double phi_rf_atbn

◆ 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_bn_xy

Vector r_bn_xy

◆ r_enterice_array

double r_enterice_array[3]

◆ r_in_array

double r_in_array[3]

◆ RANDOMISEPOL

double RANDOMISEPOL =0.

◆ realtime_this

double realtime_this

◆ rf_xy

Vector rf_xy

◆ 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_prob

double sum_prob[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_nnu_atbn

double theta_nnu_atbn

◆ theta_nnu_rf_diff_atbn

double theta_nnu_rf_diff_atbn

◆ 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

◆ xsecParam_nuint

int xsecParam_nuint = 1

◆ xsecParam_nutype

int xsecParam_nutype = 0
Generated by   doxygen 1.8.14