doxygen_icemc/html/_taumodel_8hh_source.html

//class Taumodel:

#include <algorithm>

#include <numeric>

#include <iostream>

#include <string>

#include <sstream>

#include <fstream>

#include <iomanip>

#include <vector>

#ifndef __CINT__

#include "Primaries.h"

#include "Settings.h"

#endif

class TH1F;

class Vector;

class Settings;

class TRandom3;

class Primaries;

class Interaction;

class IceModel;


using std::vector;

using std::ifstream;

using std::string;


class Taumodel: public TObject

{


private: //stuff other programs arent allowed to touch


  double B0,B1,E0;


  double mT;

  double cT;

  double p;


  //p, the Density of Standard Rock. g/cm^3

  double A;

  double Mn;

  vector<double> mydensityvector;

  vector<double> myavgdensityvector;

  vector<double> myenergyvector;

  vector<double> myPsurvvector;

  vector<double> etaufarray;

  vector<double> PDFarray;

  void GetDensityVectors(IceModel *antarctica1,Interaction *interaction1, Vector nchord, double step, double Distance,int &totalnusteps,int &crust_entered, int &mantle_entered, int &core_entered);

  void GetEnergyVector(double Etau_final, double step,int totalnusteps, int &totalsteps, double &totaltaudistance, double pnu);


  void GetTauSurvVector(double step, int totalsteps);

public:

  Taumodel();

  double ptauf;

  double weight_tau_prob;

  int inu;

  double weight_nu_prob;


 double GetTauWeight(Primaries *primary1, Settings *settings1,IceModel*antarctica1,Interaction *interaction1,

             double pnu, int nu_nubar,

             double& ptauf,

             int& crust_entered, // 1 or 0

             int& mantle_entered, // 1 or 0

             int& core_entered);//include secondaries?


}; //class Taumodel