Taumodel.hh

#ifndef TAU_MODEL_H
#define TAU_MODEL_H
//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 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;/*!<Density of Standard Rock. g/cm^3> */
  
  //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);
  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?
 
  ClassDef(Taumodel,1);
}; //class Taumodel


#endif