icemc
 All Classes Namespaces Files Functions Variables Friends Macros Pages
Functions | Variables
icemc.cc File Reference
#include <iostream>
#include <fstream>
#include <sstream>
#include <math.h>
#include <string>
#include <stdio.h>
#include <stdlib.h>
#include <vector>
#include <time.h>
#include "TTreeIndex.h"
#include "TChain.h"
#include "TH1.h"
#include "TF1.h"
#include "TF2.h"
#include "TFile.h"
#include "TRandom.h"
#include "TRandom2.h"
#include "TRandom3.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 "trigger.hh"
#include "Spectra.h"
#include "signal.hh"
#include "secondaries.hh"
#include "ray.hh"
#include "counting.hh"
#include "Primaries.h"
#include "Taumodel.hh"
#include <typeinfo>

Functions

 ClassImp (RX)
 
void GetSmearedIncidentAngle (Vector &specular, Vector &nsurf_rfexit, Vector &nrf_iceside, Vector &n_exit2bn, double SMEARINCIDENTANGLE, double theta_inc_smeared)
 
void GetScreenXY (int nrays, double screensize, Vector plusx_onscreen, Vector plusy_onscreen, Vector *position_on_screen, int &nx, Vector *coord_on_screen_vector)
 
double myfunction (double theta_inc, double theta_trans_perp, double theta_trans_par, double &norm)
 
double integratemyfunction (double theta_perp, double theta_par)
 
double GetAirDistance (double altitude_bn, double beta)
 
void SetupViewangles (Signal *sig1)
 
void ReadRoughnessData (Settings *settings1, double trans_angle[9][1000], double transmission[9][1000], double *minangle, double *maxangle, int *npoints, int &max_angles_backplane)
 
double GetTransmission (double inc_angle, double trans_angle_par, double trans_angle_perp, double trans_angle[9][1000], double transmission[9][1000], double *min_angle, double *max_angle, int *npoints, int max_angles_backplane)
 
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)
 
int GetDirection (Settings *settings1, Interaction *interaction1, const Vector &refr, double deltheta_em, double deltheta_had, double emfrac, double hadfrac, double pnu, double vmmhz1m_max, double r_fromballoon, double rough_sigma, Ray *ray1, Signal *sig1, Position posnu, Anita *anita1, Balloon *bn1, Vector &nnu, double &costhetanu, double &theta_threshold)
 
int GetDirectionRough (Settings *settings1, const Vector &refr, double deltheta_em, double deltheta_had, double emfrac, double hadfrac, double pnu, double vmmhz1m_max, double r_fromballoon, double rough_sigma, 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)
 
int Getibnposition ()
 
void IntegrateBands (Anita *anita1, int k, double *vmmhz, 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)
 
int main (int argc, char **argv)
 
void WriteNeutrinoInfo (Position &posnu, Vector &nnu, Position &r_bn, double altitude_int, double pnu, string nuflavor, string current, double elast_y, ofstream &nu_out)
 
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, ofstream &foutseckel)
 
void CloseTFile (TFile *hfile)
 
TStyle * RootStyle ()
 

Variables

Taumodel * TauPtr = NULL
 
const int NVIEWANGLE =100
 
int irays
 
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 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 logweight2 =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 posnu_down [3]
 
double gain_dipole =2.15
 
double changle_deg =0
 
double NMEDIUM_INTERACTION
 
int NNU
 
int whichray =0
 
double RANDOMISEPOL =0.
 
const double H_RX = 37.62 * CMINCH/100.
 
const double W_RX = 37.62 * CMINCH/100.
 
const double L_RX = 22.63 * CMINCH/100.
 
const double Z_RX = 50
 
double volume_thishorizon
 
int realtime_this
 
double longitude_this
 
double latitude_this
 
double altitude_this
 
double heading_this =0.
 
double pnu =pow(10.,20)
 
const bool SANITY = true
 
double trans_angle [9][1000]
 
double transmission [9][1000]
 
double minangle_roughness [9]
 
double maxangle_roughness [9]
 
int npoints_roughness [9]
 
int max_angles_backplane
 
double MAXMEASUREDANGLE =75.
 
Vector specular_vector
 
Vector reference_vector
 
Vector aimthescreen_vector
 
double SMEARINCIDENTANGLE =40.*RADDEG
 
double theta_inc_smeared
 
const int NRAYS =40000
 
int nrays =NRAYS
 
Vector diffx_on_screen_vector
 
Vector diffy_on_screen_vector
 
double posx_on_screen
 
double posy_on_screen
 
double posz_on_screen
 
double coordx_on_screen
 
double coordy_on_screen
 
double coordz_on_screen
 
double theta_inc_eachray [NRAYS]
 
double distance_eachray [NRAYS]
 
double time_eachray [NRAYS]
 
double time_eachray_fortree
 
double time_specular
 
double myfunction_eachray [NRAYS]
 
double myfunction_fortree
 
double theta_inc_specular
 
double theta_trans_specular
 
double theta_trans_eachray [NRAYS]
 
double distance_fortree
 
double theta_rfexit
 
double theta_specular
 
double theta_perptoplane
 
double theta_iceside
 
double duh
 
double theta_fromtir
 
Vector vector_temp
 
Vector plusz =Vector(0.,0.,1.)
 
Vector plusy =Vector(0.,1.,0.)
 
double phi_temp
 
Vector position_on_screen_vector [NRAYS]
 
Vector coord_on_screen_vector [NRAYS]
 
Vector rotation_axis_vector1 [NRAYS]
 
Vector rotation_axis_vector2 [NRAYS]
 
const double SCREENSIZE =100000.0
 
double screensize_fortree =SCREENSIZE
 
int nx
 
double rtemp
 
double Rtemp
 
double theta_temp
 
double sign_temp
 
double sintheta1
 
double costheta1
 
double sintheta2
 
double costheta2
 
double theta2
 
Vector n_exit2bn_eachray [NRAYS]
 
Vector rfexit_eachray [NRAYS]
 
double diffx_rfexit
 
double diffy_rfexit
 
double diffz_rfexit
 
double theta_trans_perp_fortree
 
double theta_trans_par_fortree
 
double dtheta_trans_perp_fortree
 
double dtheta_trans_par_fortree
 
double min_position =-SCREENSIZE/2.
 
double viewangle_eachray [NRAYS]
 
double viewangle_fortree
 
double theta_inc_fortree
 
double theta_trans_fortree
 
Vector nrf_iceside_eachray [NRAYS]
 
Vector component_perp_eachray [NRAYS]
 
Vector component_par_eachray [NRAYS]
 
Vector perp_to_plane_thisray
 
Vector perp_to_plane_specular
 
Vector plusx_onscreen
 
Vector plusy_onscreen
 
Vector temp_vector
 
Vector perp_to_nrf_iceside_specular_in_plane
 
Vector component_perp_to_nrf_iceside_specular_in_plane
 
double theta_trans_perp [NRAYS]
 
double theta_trans_par [NRAYS]
 
double cosine_temp
 
double dtheta_trans_perp [NRAYS]
 
double dtheta_trans_par [NRAYS]
 
double contribution_eachray [NRAYS]
 
double sumcontribution =0
 
double roughnessfactor =0
 
double contribution_total
 
double normalization
 
double contribution_max
 
double contribution_fortree
 
double integralmyfunction
 
double SIGMA_POSX =0.1
 
double SIGMA_NEGX =0.1
 
double SIGMA_POSY =0.1
 
double SIGMA_NEGY =0.1
 
double MEANX =0
 
double MEANY =0.
 
double SIGNALRADIUS =2.
 
double relativetaper_eachray [NRAYS]
 
double relativetaper_fortree
 
double relativetaper_specular
 
const double FREQ_LOW_DISCONES =120.E6
 
const double FREQ_HIGH_DISCONES =1000.E6
 
double bwslice_volts_db_pol0 [4]
 
double bwslice_volts_db_pol1 [4]
 
double bwslice_vnoise_thislayer [4]
 
double powerthresh_thischannel =0.
 
const double INTEGRATION_TIME =3.75E-9
 
int unmasked_thischannel =0
 
int passes_thisevent =0
 
int unmasked_thisevent =0
 
double sum_thisslice =0
 
int ibw_fortree =0
 
int isurf
 
int ichan
 
int lr =0
 
int channels_passing_db [4][16][2][4] = {{{{0}}}}
 
TF1 * f [4][4]
 
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
 
double fresnel1_eachboresight [Anita::NLAYERS_MAX][Anita::NPHI_MAX]
 
double fresnel2 =0
 
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_2tries =0
 
double diff_3tries =0
 
double diffnorm =0
 
double diffrefr =0
 
double costheta_inc =0
 
double costheta_refr =0
 
double costheta_exit =0
 
double theta_rf_atbn
 
double theta_rf_atbn_measured
 
double theta_inc =0
 
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 undetectable_counter =0
 
const int NTHETA_ROUGH =100
 
const int NPHI_ROUGH =100
 
double nrf_neighbors [10000]
 
int maxarray [NTHETA_ROUGH]
 
double SIGMA_HEIGHT =10.
 
double FEATURE_LENGTH =100.
 
double NEIGHBORS_RADIUS =1000.
 
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_db =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}
 
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 e_component_db =0
 
double h_component_db =0
 
double rflength_db =0
 
double ptau =0
 
int count_passestrigger_db =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
 
double costheta_inc_db =0
 
double costheta_refr_db =0
 
double costheta_exit_db =0
 
double diff_angle_db =0
 
Vector n_exit2bn_db [5]
 
Vector nrf_iceside_db [5]
 
double n_exit_phi
 
int count_nutaucc =0
 
int count_dbexitsice =0
 
int count_dboccurs =0
 
int count_dbinhorizon =0
 
int count_goodray_db =0
 
int count_firstbang =0
 
int count_interacts_in_earth =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 maxtaper_db =0
 
double dviewangle_deg =0
 
double dviewangle_db_deg =0
 
double downwardchord1 =0.
 
double downwardchord2 =0.
 
double downwardchord3 =0.
 
double downwardchord4 =0.
 
double diffdownwardchord1 =0
 
double diffdownwardchord2 =0
 
double diffdownwardchord3 =0
 
double diffdiffdownwardchord1 =0
 
double diffdiffdownwardchord2 =0
 
double forseckel [NVIEWANGLE][Anita::NFREQ]
 
double poyntingflux_1m =0.
 
double viewangles [NVIEWANGLE]
 
int err =0
 
int ilon_tmp =0
 
int ilat_tmp =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 total_kgm2_fordebug = 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 posnu2_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 secondary_tries =0
 
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
 
bool ABORT_EARLY = false
 
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 thresholds [NTHRESHOLDS]
 

Function Documentation

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 
)
ClassImp ( RX  )
void CloseTFile ( TFile *  hfile)
void GetAir ( double *  col1)
double GetAirDistance ( double  altitude_bn,
double  beta 
)
double GetAverageVoltageFromAntennasHit ( Settings *  settings1,
int *  nchannels_perrx_triggered,
double *  voltagearray,
double &  volts_rx_sum 
)
void GetCurrent ( string &  current)
int GetDirection ( Settings *  settings1,
Interaction *  interaction1,
const Vector &  refr,
double  deltheta_em,
double  deltheta_had,
double  emfrac,
double  hadfrac,
double  pnu,
double  vmmhz1m_max,
double  r_fromballoon,
double  rough_sigma,
Ray *  ray1,
Signal *  sig1,
Position  posnu,
Anita *  anita1,
Balloon *  bn1,
Vector &  nnu,
double &  costhetanu,
double &  theta_threshold 
)
int GetDirectionRough ( Settings *  settings1,
const Vector &  refr,
double  deltheta_em,
double  deltheta_had,
double  emfrac,
double  hadfrac,
double  pnu,
double  vmmhz1m_max,
double  r_fromballoon,
double  rough_sigma,
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 
)
int Getibnposition ( )
Vector GetPolarization ( const Vector &  nnu,
const Vector &  nrf2_iceside 
)
int GetRayIceSide ( const Vector &  n_exit2rx,
const Vector &  nsurf_rfexit,
double  nexit,
double  nenter,
Vector &  nrf2_iceside 
)
void GetScreenXY ( int  nrays,
double  screensize,
Vector  plusx_onscreen,
Vector  plusy_onscreen,
Vector *  position_on_screen,
int &  nx,
Vector *  coord_on_screen_vector 
)
void GetSmearedIncidentAngle ( Vector &  specular,
Vector &  nsurf_rfexit,
Vector &  nrf_iceside,
Vector &  n_exit2bn,
double  SMEARINCIDENTANGLE,
double  theta_inc_smeared 
)
double GetThisAirColumn ( Settings *  settings1,
Position  r_in,
Vector  nnu,
Position  posnu,
double *  col1,
double &  cosalpha,
double &  mytheta,
double &  cosbeta0,
double &  mybeta 
)
double GetTransmission ( double  inc_angle,
double  trans_angle_par,
double  trans_angle_perp,
double  trans_angle[9][1000],
double  transmission[9][1000],
double *  min_angle,
double *  max_angle,
int *  npoints,
int  max_angles_backplane 
)
double GetViewAngle ( const Vector &  nrf2_iceside,
const Vector &  nnu 
)
void Integrate ( Anita *  anita1,
int  j,
int  k,
double *  vmmhz,
double *  freq,
double  scalefactor,
double  sumsignal 
)
void IntegrateBands ( Anita *  anita1,
int  k,
double *  vmmhz,
double *  freq,
double  scalefactor,
double *  sumsignal 
)
double integratemyfunction ( double  theta_perp,
double  theta_par 
)
void interrupt_signal_handler ( int  sig)
void IsAbsorbed ( double  chord_kgm2,
double  len_int_kgm2,
double &  weight 
)
double IsItDoubleBang ( double  exitlength,
double  plepton 
)
int main ( int  argc,
char **  argv 
)

< tau trigger sets the chance to create tau particle

double myfunction ( double  theta_inc,
double  theta_trans_perp,
double  theta_trans_par,
double &  norm 
)
void ReadRoughnessData ( Settings *  settings1,
double  trans_angle[9][1000],
double  transmission[9][1000],
double *  minangle,
double *  maxangle,
int *  npoints,
int &  max_angles_backplane 
)
TStyle* RootStyle ( )
double ScaleVmMHz ( double  vmmhz1m_max,
const Position &  posnu,
const Position &  r_bn 
)
void SetupViewangles ( Signal *  sig1)
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,
ofstream &  foutseckel 
)
int TIR ( const Vector &  n_surf,
const Vector &  nrf2_iceside,
double  N_IN,
double  N_OUT 
)
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 WriteNeutrinoInfo ( Position &  posnu,
Vector &  nnu,
Position &  r_bn,
double  altitude_int,
double  pnu,
string  nuflavor,
string  current,
double  elast_y,
ofstream &  nu_out 
)

Variable Documentation

bool ABORT_EARLY = false
Vector aimthescreen_vector
int allcuts[2] ={0,0}
double allcuts_weighted[2] ={0,0}
double altitude_this
Vector ant_max_normal0
double ant_max_normal0_array[3]
Vector ant_max_normal1
double ant_max_normal1_array[3]
Vector ant_max_normal2
double ant_max_normal2_array[3]
double averaging_thetas1 =0
double averaging_thetas2 =0
double averaging_thetas3 =0
double BW_DISCONES =300.E6-120.E6
double bwslice_vnoise_thislayer[4]
double bwslice_volts_db_pol0[4]
double bwslice_volts_db_pol1[4]
int chanceinhell2 =0
double changle_deg =0
int channels_passing_db[4][16][2][4] = {{{{0}}}}
Vector component_par_eachray[NRAYS]
Vector component_perp_eachray[NRAYS]
Vector component_perp_to_nrf_iceside_specular_in_plane
double contribution_eachray[NRAYS]
double contribution_fortree
double contribution_max
double contribution_total
Vector coord_on_screen_vector[NRAYS]
double coordx_on_screen
double coordy_on_screen
double coordz_on_screen
int core_entered =0
double cosine_temp
double costheta1
double costheta2
double costheta_exit =0
double costheta_exit_db =0
double costheta_inc =0
double costheta_inc_db =0
double costheta_refr =0
double costheta_refr_db =0
double costhetanu =-1000
int count_asktrigger =0
int count_asktrigger_db =0
int count_asktrigger_nfb =0
int count_chanceinhell0 =0
double count_chanceinhell2_w =0
int count_chanceofsurviving =0
int count_chordgoodlength =0
int count_d2goodlength =0
int count_dbexitsice =0
int count_dbinhorizon =0
int count_dboccurs =0
int count_firstbang =0
int count_goodray_db =0
int count_interacts_in_earth =0
int count_inthisloop1 =0
int count_inthisloop2 =0
int count_inthisloop3 =0
int count_nutaucc =0
int count_pass =0
int count_passestrigger_db =0
int count_passestrigger_nfb =0
double count_passestrigger_w =0
int count_rx =0
int count_total =0
int crust_entered =0
double CUTONWEIGHTS =1.E-10
double denom_v_thresh[NTHRESHOLDS] ={0.}
double diff_2tries =0
double diff_3tries =0
double diff_angle_db =0
double diffdiffdownwardchord1 =0
double diffdiffdownwardchord2 =0
double diffdownwardchord1 =0
double diffdownwardchord2 =0
double diffdownwardchord3 =0
double diffexit =0
double diffnorm =0
double diffrefr =0
Vector diffx_on_screen_vector
double diffx_rfexit
Vector diffy_on_screen_vector
double diffy_rfexit
double diffz_rfexit
int discones_passing
double distance_eachray[NRAYS]
double distance_fortree
double djunk
double downwardchord1 =0.
double downwardchord2 =0.
double downwardchord3 =0.
double downwardchord4 =0.
double dtheta_trans_par[NRAYS]
double dtheta_trans_par_fortree
double dtheta_trans_perp[NRAYS]
double dtheta_trans_perp_fortree
double duh
double dviewangle_db_deg =0
double dviewangle_deg =0
double e_comp_max1 =0
double e_comp_max2 =0
double e_comp_max3 =0
double e_component_db =0
int err =0
double error_distance_minus[NBINS_DISTANCE] = {0.}
double error_distance_plus[NBINS_DISTANCE] = {0.}
double error_e_minus =0
double error_e_plus =0
double error_km3sr_nfb =0
double error_minus =0
double error_mu_minus =0
double error_mu_plus =0
double error_nfb =0
double error_percent_increase_nfb =0
double error_plus =0
double error_tau_minus =0
double error_tau_plus =0
double eventsfound =0
double eventsfound_beforetrigger =0.
double eventsfound_belowhorizon =0
double eventsfound_binned[NBINS]
double eventsfound_binned_distance[NBINS_DISTANCE] = {0.}
int eventsfound_binned_distance_forerror[NBINS_DISTANCE][NBINS] = {{0}}
double eventsfound_binned_e[NBINS]
double eventsfound_binned_mu[NBINS]
double eventsfound_binned_tau[NBINS]
double eventsfound_crust =0
double eventsfound_nfb_binned[NBINS]
double eventsfound_prob =0
double eventsfound_weightgt01 =0
TF1* f[4][4]
double FEATURE_LENGTH =100.
double forseckel[NVIEWANGLE][Anita::NFREQ]
int freq_bins = Anita::NFREQ
const double FREQ_HIGH_DISCONES =1000.E6
const double FREQ_LOW_DISCONES =120.E6
double fresnel1 =0
double fresnel1_eachboresight[Anita::NLAYERS_MAX][Anita::NPHI_MAX]
double fresnel2 =0
double gain_dipole =2.15
double h_comp_max1 =0
double h_comp_max2 =0
double h_comp_max3 =0
double h_component_db =0
const double H_RX = 37.62 * CMINCH/100.
double heading_this =0.
double heff_discone =0
double heff_max =0.62639
int ibw_fortree =0
int ichan
int ierr =0
int ilat_tmp =0
int ilon_tmp =0
int index_distance =0
int index_weights =0
double integralmyfunction
const double INTEGRATION_TIME =3.75E-9
int inu
int irays
int isurf
double km3sr = 0
double km3sr_db = 0
double km3sr_distance[NBINS_DISTANCE] = {0.}
double km3sr_e =0
double km3sr_mu =0
double km3sr_nfb =0
double km3sr_tau =0
const double L_RX = 22.63 * CMINCH/100.
double lat_in =0
double latitude_this
double len_int =0
double logweight =0.
double logweight2 =0
double longitude_this
int lr =0
double mag1 =0
double mag1_eachboresight[Anita::NLAYERS_MAX][Anita::NPHI_MAX]
double mag2 =0
int mantle_entered =0
int max_angles_backplane
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 MAX_ATTENLENGTH =1671
double MAX_LOGWEIGHT =-1
double maxangle_roughness[9]
int maxarray[NTHETA_ROUGH]
double MAXMEASUREDANGLE =75.
double maxtaper =0
double maxtaper_db =0
double MEANX =0
double MEANY =0.
double MIN_LOGWEIGHT =-3
double min_position =-SCREENSIZE/2.
double minangle_roughness[9]
double myfunction_eachray[NRAYS]
double myfunction_fortree
double n_bn_array[3]
double n_exit2bn_array[5][3]
Vector n_exit2bn_db[5]
double n_exit2bn_db_array[5][3]
Vector n_exit2bn_eachray[NRAYS]
double n_exit_phi
double n_pol_array[3]
const int NBINS =10
const int NBINS_DISTANCE =28
int NDISCONES =8
double nearthlayers =0
double NEIGHBORS_RADIUS =1000.
int neutrinos_passing_all_cuts =0
double NMEDIUM_INTERACTION
int NNU
double nnu_array[3]
double normalization
double npass_v_thresh[NTHRESHOLDS] ={0.}
const int NPHI_ROUGH =100
int npoints_roughness[9]
const int NRAYS =40000
int nrays =NRAYS
double nrf_iceside_array[5][3]
Vector nrf_iceside_db[5]
double nrf_iceside_db_array[5][3]
Vector nrf_iceside_eachray[NRAYS]
double nrf_neighbors[10000]
double nsurf_rfexit_array[3]
double nsurf_rfexit_db_array[3]
const int NTHETA_ROUGH =100
const int NTHRESHOLDS =20
const int NVIEWANGLE =100
int nx
int passes_thisevent =0
int passestrigger =0
double percent_increase_db =0
double percent_increase_nfb =0
double percent_increase_total =0
Vector perp_to_nrf_iceside_specular_in_plane
Vector perp_to_plane_specular
Vector perp_to_plane_thisray
double phi_temp
double pieceofkm2sr =0
Vector plusx_onscreen
Vector plusy =Vector(0.,1.,0.)
Vector plusy_onscreen
Vector plusz =Vector(0.,0.,1.)
double pnu =pow(10.,20)
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.}
double polarfactor_discone =0.
Vector position_on_screen_vector[NRAYS]
double posnu2_array[3]
double posnu_array[3]
double posnu_down[3]
double posx_on_screen
double posy_on_screen
double posz_on_screen
double powerthresh_thischannel =0.
double poyntingflux_1m =0.
double ptau =0
double r_bn_array[3]
double r_enterice_array[3]
double r_in_array[3]
double RANDOMISEPOL =0.
int realtime_this
Vector reference_vector
double relativetaper_eachray[NRAYS]
double relativetaper_fortree
double relativetaper_specular
double rfexit_array[5][3]
double rfexit_db_array[5][3]
Vector rfexit_eachray[NRAYS]
double rflength =0
double rflength_db =0
Vector rotation_axis_vector1[NRAYS]
Vector rotation_axis_vector2[NRAYS]
double roughnessfactor =0
double rtemp
double Rtemp
double rx0_noise_eachband[2][5]
int rx0_passes_eachband[2][5]
double rx0_signal_eachband[2][5]
double rx0_threshold_eachband[2][5]
const bool SANITY = true
double scalefactor_attenuation =0
double scalefactor_distance =0
const double SCREENSIZE =100000.0
double screensize_fortree =SCREENSIZE
int secondary_tries =0
double SIGMA_HEIGHT =10.
double SIGMA_NEGX =0.1
double SIGMA_NEGY =0.1
double SIGMA_POSX =0.1
double SIGMA_POSY =0.1
double sign_temp
double SIGNALRADIUS =2.
double sintheta1
double sintheta2
double SMEARINCIDENTANGLE =40.*RADDEG
Vector specular_vector
double sum[3]
double sum_frac[3]
double sum_frac_db[3]
double sum_thisslice =0
double sum_weights =0
double sumcontribution =0
double t_coeff_pokey
double t_coeff_slappy
Taumodel* TauPtr = NULL
Vector temp_vector
double theta2
double theta_fromtir
double theta_iceside
double theta_in =0
double theta_inc =0
double theta_inc_eachray[NRAYS]
double theta_inc_fortree
double theta_inc_smeared
double theta_inc_specular
double theta_perptoplane
double theta_rf_atbn
double theta_rf_atbn_measured
double theta_rfexit
double theta_specular
double theta_temp
double theta_trans_eachray[NRAYS]
double theta_trans_fortree
double theta_trans_par[NRAYS]
double theta_trans_par_fortree
double theta_trans_perp[NRAYS]
double theta_trans_perp_fortree
double theta_trans_specular
double thislambda =0
double threshold_end =-6.
double threshold_start =-1.
double threshold_step =(threshold_end-threshold_start)/(double)NTHRESHOLDS
double thresholds[NTHRESHOLDS]
double time_eachray[NRAYS]
double time_eachray_fortree
double time_specular
int times_core_entered_det =0
int times_crust_entered_det =0
int times_mantle_entered_det =0
double total_kgm2 = 0
double total_kgm2_fordebug = 0
double trans_angle[9][1000]
double transmission[9][1000]
int undetectable_counter =0
int unmasked_thischannel =0
int unmasked_thisevent =0
Vector vector_temp
double viewangle_eachray[NRAYS]
double viewangle_fortree
double viewangles[NVIEWANGLE]
double vmmhz1m_visible = 0
double vnoise_discone =0.
double voltagearray[Anita::NLAYERS_MAX *Anita::NPHI_MAX]
double volts_discone =0.
double volume_thishorizon
const double W_RX = 37.62 * CMINCH/100.
double weight =0.
double weight1 =0
double weight_prob =0.
int whichray =0
const double Z_RX = 50
Generated on Sat Nov 15 2014 22:18:21 for icemc by   doxygen 1.8.3.1