Jpp  18.4.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Classes | Functions | Variables
JSIRENE Namespace Reference

Detector simulations. More...

Classes

struct  JFoxWolframMoments
 Class for computing Fox-Wolfram moments. More...
 
struct  JSphericityTensor
 Class for sphericity tensor calculations. More...
 
struct  JHitInertiaTensor
 Class for hit inertia tensor calculations. More...
 
class  JPulse
 Auxiliary class for a time-over-threshold pulse from a PMT. More...
 
struct  JPythia
 Auxiliary class to determine EM-equivalent energy as a function of PDG particle code and energy. More...
 
class  JSeaWater
 Sea water composition. More...
 
struct  JParameters
 Detector simulation parameters. More...
 
struct  JHit_t
 Auxiliary class to set-up Hit. More...
 
struct  JHits_t
 Auxiliary data structure for list of hits with hit merging capability. More...
 
struct  JTrk_t
 Auxiliary class to set-up Trk. More...
 
struct  JPoint
 Point along muon trajectory. More...
 
struct  JTrack
 Muon trajectory. More...
 
struct  JVertex
 Vertex of energy loss of muon. More...
 

Functions

double getThrust (const std::vector< Trk >::const_iterator __begin, const std::vector< Trk >::const_iterator __end)
 Compute thrust for a given range of tracks. More...
 
double getThrust (const Evt &event)
 Compute thrust for a given event. More...
 
bool operator< (const JPulse &first, const JPulse &second)
 Compare Monte Carlo hit times. More...
 
bool operator< (const JPulse &hit, const double t0)
 Compare Monte Carlo hit times. More...
 
JHitType_t getHitType (const JPDFType_t pdf, const bool shower=false)
 Get hit type corresponding to given PDF type. More...
 
double getVisibleEnergy (const Trk &track, const JCylinder3D &can=getMaximumContainmentVolume())
 Get the visible energy of a track. More...
 
Vec getVisibleEnergyVector (const Trk &track, const JCylinder3D &can=getMaximumContainmentVolume())
 Get the visible energy vector of a track. More...
 
double getVisibleEnergy (std::vector< Trk >::const_iterator __begin, std::vector< Trk >::const_iterator __end, const JCylinder3D &can=getMaximumContainmentVolume())
 Get the visible energy of a given range of tracks. More...
 
Vec getVisibleEnergyVector (std::vector< Trk >::const_iterator __begin, std::vector< Trk >::const_iterator __end, const JCylinder3D &can=getMaximumContainmentVolume())
 Get the visible energy vector of a given range of tracks. More...
 
double getVisibleEnergy (const Evt &evt, const JCylinder3D &can=getMaximumContainmentVolume())
 Get the visible energy vector of an event. More...
 
Vec getVisibleEnergyVector (const Evt &evt, const JCylinder3D &can=getMaximumContainmentVolume())
 Get the visible energy vector of an event. More...
 
double opa_weight_high_e (const double ekin)
 
double ngamma_elec (const double ekin)
 
double weight_pion (const double ekin)
 
double weight_kaon (const double ekin)
 
double weight_kshort (const double ekin)
 
double weight_klong (const double ekin)
 
double weight_proton (const double ekin)
 
double weight_neutron (const double ekin)
 
double opa_efrac (const int ipart, const double ekin)
 
double pythia (const int type, const double E)
 Get equivalent EM-energy for given pion energy. More...
 

Variables

static const JPythia pythia
 Function object for relative light yield as a function of GEANT particle code. More...
 
struct {
   std::vector< size_t >   JSIRENE::ns
 
   std::vector< double >   JSIRENE::vs
 
} getNumberOfPhotoElectrons
 Auxiliary data structure for determination of number of photo-electrons. More...
 

Detailed Description

Detector simulations.

Author
mdejong

Function Documentation

double JSIRENE::getThrust ( const std::vector< Trk >::const_iterator  __begin,
const std::vector< Trk >::const_iterator  __end 
)
inline

Compute thrust for a given range of tracks.


The definition was taken from the description in section 15.2.2 of arXiv:hep-ph/0603175v2.

Parameters
__beginbeginning of track data
__endend of track data
Returns
thrust

< Thrust axis

< residual angle

Definition at line 42 of file JEventShapeVariables.hh.

44  {
45  using namespace std;
46  using namespace JPP;
47 
48  static const int MAX_STEPS = 20;
49  static const double epsilon = 1e-9;
50 
51 
52  Vec axis0(0.0, 0.0, 0.0); //!< Thrust axis
53 
54  // Set thrust axis start value (= summed final state particle momenta)
55 
56  for (vector<Trk>::const_iterator i = __begin; i != __end; ++i) {
57  if (is_finalstate(*i)) {
58  axis0 += getKineticEnergy(*i) * i->dir;
59  }
60  }
61 
62  if (axis0.len() < epsilon) { // If total final state momentum is zero, set start value equal to first final state track momentum
63 
64  vector<Trk>::const_iterator i = find_if(__begin, __end, &is_finalstate);
65 
66  if (i != __end) {
67  axis0 = getKineticEnergy(*i) * i->dir;
68  } else {
69  THROW(JValueOutOfRange, "getThrust(): No final state tracks given.");
70  }
71  }
72 
73  axis0.normalize();
74 
75  // Find axis which maximizes the thrust
76 
77  double res = numeric_limits<double>::max(); //!< residual angle
78 
79  for (int k = 0; k < MAX_STEPS && res > epsilon; ++k) {
80 
81  Vec axis1(0.0, 0.0, 0.0);
82 
83  for (vector<Trk>::const_iterator i = __begin; i != __end; ++i) {
84 
85  if (is_finalstate(*i)) {
86 
87  const Vec p = getKineticEnergy(*i) * i->dir;
88 
89  axis1 += (p.dot(axis0) > 0.0 ? p : -p);
90  }
91  }
92 
93  axis1.normalize();
94 
95  res = acos( axis0.dot(axis1) );
96 
97  axis0 = axis1;
98  }
99 
100  // Compute thrust
101 
102  double pSum = 0.0;
103  double norm = 0.0;
104 
105  for (vector<Trk>::const_iterator i = __begin; i != __end; ++i) {
106 
107  if (is_finalstate(*i)) {
108 
109  const Vec p = getKineticEnergy(*i) * i->dir;
110 
111  pSum += fabs(p.dot(axis0));
112  norm += p.len();
113  }
114  }
115 
116  return pSum / norm;
117  }
k
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
Definition: hydrophone-phi:fit.sh:35
THROW
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
Definition: JException.hh:712
Vec::len
double len() const
Get length.
Definition: Vec.hh:145
JAANET::is_finalstate
bool is_finalstate(const Trk &track)
Test whether given track corresponds to a final state particle.
Definition: JAAnetToolkit.hh:450
Vec::dot
double dot(const Vec &v) const
Get dot product.
Definition: Vec.hh:36
Vec
The Vec class is a straightforward 3-d vector, which also works in pyroot.
Definition: Vec.hh:12
JAANET::getKineticEnergy
double getKineticEnergy(const double E, const double m)
Get kinetic energy of particle with given mass.
Definition: JAAnetToolkit.hh:772
GAUSS_LEGENDRE::epsilon
const double epsilon
Definition: JQuadrature.cc:21
double JSIRENE::getThrust ( const Evt &  event)
inline

Compute thrust for a given event.


The definition was taken from the description in section 15.2.2 of arXiv:hep-ph/0603175v2.

Parameters
eventevent
Returns
thrust

Definition at line 127 of file JEventShapeVariables.hh.

128  {
129  return getThrust(event.mc_trks.begin(), event.mc_trks.end());
130  }
Evt::mc_trks
std::vector< Trk > mc_trks
MC: list of MC truth tracks.
Definition: Evt.hh:49
JSIRENE::getThrust
double getThrust(const std::vector< Trk >::const_iterator __begin, const std::vector< Trk >::const_iterator __end)
Compute thrust for a given range of tracks.
Definition: JEventShapeVariables.hh:42
bool JSIRENE::operator< ( const JPulse &  first,
const JPulse &  second 
)
inline

Compare Monte Carlo hit times.

Parameters
firstfirst hit
secondsecond hit
Returns
true if first hit earlier than second; else false

Definition at line 81 of file JPulse.hh.

82  {
83  return first.getLowerLimit() < second.getLowerLimit();
84  }
bool JSIRENE::operator< ( const JPulse &  hit,
const double  t0 
)
inline

Compare Monte Carlo hit times.

Parameters
hithit
t0time [ns]
Returns
true if hit earlier than given time; else false

Definition at line 94 of file JPulse.hh.

95  {
96  return hit.getLowerLimit() < t0;
97  }
JHitType_t JSIRENE::getHitType ( const JPDFType_t  pdf,
const bool  shower = false 
)
inline

Get hit type corresponding to given PDF type.

Parameters
pdfPDF type
showerforce origin from neutrino interaction shower
Returns
hit type

Definition at line 151 of file JSireneToolkit.hh.

152  {
153  using namespace JAANET;
154  using namespace JPHYSICS;
155 
156  switch (pdf) {
157 
158  case DIRECT_LIGHT_FROM_MUON:
159  return HIT_TYPE_MUON_DIRECT;
160 
161  case SCATTERED_LIGHT_FROM_MUON:
162  return HIT_TYPE_MUON_SCATTERED;
163 
164  case DIRECT_LIGHT_FROM_DELTARAYS:
165  return HIT_TYPE_DELTARAYS_DIRECT;
166 
167  case SCATTERED_LIGHT_FROM_DELTARAYS:
168  return HIT_TYPE_DELTARAYS_SCATTERED;
169 
170  case DIRECT_LIGHT_FROM_EMSHOWER:
171  return (shower ? HIT_TYPE_SHOWER_DIRECT : HIT_TYPE_BREMSSTRAHLUNG_DIRECT );
172 
173  case SCATTERED_LIGHT_FROM_EMSHOWER:
174  return (shower ? HIT_TYPE_SHOWER_SCATTERED : HIT_TYPE_BREMSSTRAHLUNG_SCATTERED);
175 
176  default:
177  return HIT_TYPE_UNKNOWN;
178  }
179  }
JAANET::HIT_TYPE_SHOWER_SCATTERED
Scattered light from primary shower.
Definition: JAAnetToolkit.hh:81
JAANET::HIT_TYPE_DELTARAYS_DIRECT
Direct light from delta-rays.
Definition: JAAnetToolkit.hh:76
JAANET::HIT_TYPE_MUON_SCATTERED
Scattered light from muon.
Definition: JAAnetToolkit.hh:75
JPHYSICS::SCATTERED_LIGHT_FROM_EMSHOWER
scattered light from EM shower
Definition: JPDFTypes.hh:38
JAANET::HIT_TYPE_BREMSSTRAHLUNG_SCATTERED
Scattered light from Bremsstrahlung.
Definition: JAAnetToolkit.hh:79
JAANET::HIT_TYPE_BREMSSTRAHLUNG_DIRECT
Direct light from Bremsstrahlung.
Definition: JAAnetToolkit.hh:78
JAANET::HIT_TYPE_SHOWER_DIRECT
Direct light from primary shower.
Definition: JAAnetToolkit.hh:80
JPHYSICS::DIRECT_LIGHT_FROM_MUON
direct light from muon
Definition: JPDFTypes.hh:26
JPHYSICS::SCATTERED_LIGHT_FROM_MUON
scattered light from muon
Definition: JPDFTypes.hh:27
JPHYSICS::SCATTERED_LIGHT_FROM_DELTARAYS
scattered light from delta-rays
Definition: JPDFTypes.hh:33
JPHYSICS::DIRECT_LIGHT_FROM_EMSHOWER
direct light from EM shower
Definition: JPDFTypes.hh:37
JPHYSICS::DIRECT_LIGHT_FROM_DELTARAYS
direct light from delta-rays
Definition: JPDFTypes.hh:32
JAANET::HIT_TYPE_MUON_DIRECT
Direct light from muon.
Definition: JAAnetToolkit.hh:74
JAANET::HIT_TYPE_DELTARAYS_SCATTERED
Scattered light from delta-rays.
Definition: JAAnetToolkit.hh:77
double JSIRENE::getVisibleEnergy ( const Trk &  track,
const JCylinder3D &  can = getMaximumContainmentVolume() 
)

Get the visible energy of a track.


This method accounts for muon radiative energy losses.

Parameters
tracktrack
candetector can
Returns
visible energy [GeV]

Definition at line 48 of file JVisibleEnergyToolkit.hh.

50  {
51  using namespace std;
52  using namespace JPP;
53 
54  double Evis = 0.0;
55 
56  if (is_finalstate(track)) {
57 
58  if (is_muon(track)) {
59 
60  // Determine muon pathlength inside detector [m]
61 
62  const JCylinder3D::intersection_type& intersection = can.getIntersection(getAxis(track));
63 
64  const double Lmuon = gWater.getX(track.E, MASS_MUON / getSinThetaC());
65  const double Leff = (intersection.first < 0.0 ?
66  min(Lmuon, intersection.second) :
67  min(Lmuon, intersection.second) - intersection.first);
68 
69  // Determine visible energy deposition [GeV]
70 
71  const double dEb = gWater.getEb(track.E, Leff);
72  const double dEc = Leff / geanc();
73 
74  Evis = dEb + dEc;
75 
76  } else if (!is_neutrino(track) && JPDB::getInstance().hasPDG(track.type)) {
77 
78  Evis = pythia(track.type, getKineticEnergy(track));
79  }
80  }
81 
82  return Evis;
83  }
JPHYSICS::geanc
double geanc()
Equivalent muon track length per unit shower energy.
Definition: JGeane.hh:28
JAANET::is_muon
bool is_muon(const Trk &track)
Test whether given track is a (anti-)muon.
Definition: JAAnetToolkit.hh:327
JPHYSICS::MASS_MUON
static const double MASS_MUON
muon mass [GeV]
Definition: JPhysics/JConstants.hh:69
JPHYSICS::gWater
static const JGeaneWater gWater
Function object for energy loss of muon in sea water.
Definition: JGeane.hh:381
JPHYSICS::JGeaneWater::getX
virtual double getX(const double E0, const double E1) const override
Get distance traveled by muon.
Definition: JGeane.hh:349
JSIRENE::pythia
static const JPythia pythia
Function object for relative light yield as a function of GEANT particle code.
Definition: JPythia.hh:96
Trk::E
double E
Energy [GeV] (either MC truth or reconstructed)
Definition: Trk.hh:20
JPHYSICS::JGeaneWater::getEb
double getEb(const double E, const double dx) const
Get energy loss due to pair production and bremsstrahlung.
Definition: JGeane.hh:334
first
then echo The file $DIR KM3NeT_00000001_00000000 root already please rename or remove it first
Definition: JDOMDAQDriver.sh:57
JAANET::is_neutrino
bool is_neutrino(const Trk &track)
Test whether given track is a neutrino.
Definition: JAAnetToolkit.hh:309
JAANET::is_finalstate
bool is_finalstate(const Trk &track)
Test whether given track corresponds to a final state particle.
Definition: JAAnetToolkit.hh:450
JAANET::getKineticEnergy
double getKineticEnergy(const double E, const double m)
Get kinetic energy of particle with given mass.
Definition: JAAnetToolkit.hh:772
JLANG::getInstance
T & getInstance(const T &object)
Get static instance from temporary object.
Definition: JObject.hh:75
JAANET::getAxis
JAxis3D getAxis(const Trk &track)
Get axis.
Definition: JAAnetToolkit.hh:247
Trk::type
int type
MC: particle type in PDG encoding.
Definition: Trk.hh:24
JPHYSICS::getSinThetaC
double getSinThetaC()
Get average sine of Cherenkov angle of water corresponding to group velocity.
Definition: JPhysics/JConstants.hh:178
Vec JSIRENE::getVisibleEnergyVector ( const Trk &  track,
const JCylinder3D &  can = getMaximumContainmentVolume() 
)
inline

Get the visible energy vector of a track.


This method accounts for muon radiative energy losses.

Parameters
tracktrack
candetector can
Returns
visible energy vector [GeV]

Definition at line 94 of file JVisibleEnergyToolkit.hh.

95  {
96  return getVisibleEnergy(track, can) * track.dir;
97  }
Trk::dir
Vec dir
track direction
Definition: Trk.hh:18
JSIRENE::getVisibleEnergy
double getVisibleEnergy(const Trk &track, const JCylinder3D &can=getMaximumContainmentVolume())
Get the visible energy of a track.
Definition: JVisibleEnergyToolkit.hh:48
double JSIRENE::getVisibleEnergy ( std::vector< Trk >::const_iterator  __begin,
std::vector< Trk >::const_iterator  __end,
const JCylinder3D &  can = getMaximumContainmentVolume() 
)
inline

Get the visible energy of a given range of tracks.


This method accounts for muon radiative energy losses.

Parameters
__beginstart of track data
__endend of track data
candetector can
Returns
visible energy [GeV]

Definition at line 109 of file JVisibleEnergyToolkit.hh.

112  {
113  using namespace std;
114 
115  double Evis = 0.0;
116 
117  for (vector<Trk>::const_iterator track = __begin; track != __end; ++track) {
118  Evis += getVisibleEnergy(*track, can);
119  }
120 
121  return Evis;
122  }
JSIRENE::getVisibleEnergy
double getVisibleEnergy(const Trk &track, const JCylinder3D &can=getMaximumContainmentVolume())
Get the visible energy of a track.
Definition: JVisibleEnergyToolkit.hh:48
Vec JSIRENE::getVisibleEnergyVector ( std::vector< Trk >::const_iterator  __begin,
std::vector< Trk >::const_iterator  __end,
const JCylinder3D &  can = getMaximumContainmentVolume() 
)
inline

Get the visible energy vector of a given range of tracks.


This method accounts for muon radiative energy losses.

Parameters
__beginstart of track data
__endend of track data
candetector can
Returns
visible energy vector [GeV]

Definition at line 134 of file JVisibleEnergyToolkit.hh.

137  {
138  using namespace std;
139 
140  Vec Evis(0.0, 0.0, 0.0);
141 
142  for (vector<Trk>::const_iterator track = __begin; track != __end; ++track) {
143  Evis += getVisibleEnergyVector(*track, can);
144  }
145 
146  return Evis;
147  }
JSIRENE::getVisibleEnergyVector
Vec getVisibleEnergyVector(const Trk &track, const JCylinder3D &can=getMaximumContainmentVolume())
Get the visible energy vector of a track.
Definition: JVisibleEnergyToolkit.hh:94
Vec
The Vec class is a straightforward 3-d vector, which also works in pyroot.
Definition: Vec.hh:12
double JSIRENE::getVisibleEnergy ( const Evt &  evt,
const JCylinder3D &  can = getMaximumContainmentVolume() 
)
inline

Get the visible energy vector of an event.


This method accounts for muon radiative energy losses.

Parameters
evtevent
candetector can
Returns
visible energy [GeV]

Definition at line 158 of file JVisibleEnergyToolkit.hh.

160  {
161  return getVisibleEnergy(evt.mc_trks.begin(), evt.mc_trks.end(), can);
162  }
JSIRENE::getVisibleEnergy
double getVisibleEnergy(const Trk &track, const JCylinder3D &can=getMaximumContainmentVolume())
Get the visible energy of a track.
Definition: JVisibleEnergyToolkit.hh:48
Evt::mc_trks
std::vector< Trk > mc_trks
MC: list of MC truth tracks.
Definition: Evt.hh:49
Vec JSIRENE::getVisibleEnergyVector ( const Evt &  evt,
const JCylinder3D &  can = getMaximumContainmentVolume() 
)
inline

Get the visible energy vector of an event.


This method accounts for muon radiative energy losses.

Parameters
evtevent
candetector can
Returns
visible energy vector [GeV]

Definition at line 173 of file JVisibleEnergyToolkit.hh.

175  {
176  return getVisibleEnergyVector(evt.mc_trks.begin(), evt.mc_trks.end(), can);
177  }
JSIRENE::getVisibleEnergyVector
Vec getVisibleEnergyVector(const Trk &track, const JCylinder3D &can=getMaximumContainmentVolume())
Get the visible energy vector of a track.
Definition: JVisibleEnergyToolkit.hh:94
Evt::mc_trks
std::vector< Trk > mc_trks
MC: list of MC truth tracks.
Definition: Evt.hh:49
double JSIRENE::opa_weight_high_e ( const double  ekin)
inline

Definition at line 20 of file pythia.hh.

21  {
22  static const double Ma = 72.425;
23  static const double Mb = -49.417;
24  static const double Mc = 5.858;
25  static const double Md = 207.252;
26  static const double Me = 132.784;
27  static const double Mf = -10.277;
28  static const double Mg = -19.441;
29  static const double Mh = 58.598;
30  static const double Mi = 53.161;
31  static const double Mkref = 2.698;
32  // leading_coef Mlc=(Ma-Mf)/Mkref
33  static const double Mlc = (Ma-Mf)/Mkref;
34 
35  double num, denom, lognp, Ee, logE;
36 
37  // implement the energy-fractions according to M. Dentler: ANTARES-SOFT-2012-009. E must be the kinetic energy in GeV. Here, 'everything else' is treated as a pion. Other than charged pions, mostly kaons of all varieties and protons fall into this category. These fits were made mostly at high energies - at low energies, other fits are used.
38 
39  if (ekin > 0.2)
40  logE=log10(ekin);
41  else
42  // nphotons constant below 0.2 GeV (lowest fit point)
43  return 0.292;
44 
45  num = Me + Md*logE + Mc*pow(logE,2) + Mb*pow(logE,3) + Ma *pow(logE,4) + Mlc*pow(logE,5);
46  denom = Mi + Mh*logE + Mg*pow(logE,2) + Mf*pow(logE,3) + Mlc*pow(logE,4);
47 
48  if (denom <= 0) {
49  return 0.;
50  }
51 
52  lognp = num/denom;
53  Ee = pow(10.0, lognp-Mkref);
54 
55  return Ee/ekin;
56  }
log10
set_variable E_E log10(E_{fit}/E_{#mu})"
JMATH::pow
T pow(const T &x, const double y)
Power .
Definition: JMath.hh:97
double JSIRENE::ngamma_elec ( const double  ekin)
inline

Definition at line 58 of file pythia.hh.

59  {
60  static const double eleca = 1.33356e5;
61  static const double elecb = 1.66113e2;
62  static const double elecc = 16.4949;
63  static const double elecd = 1.5385e5;
64  static const double elece = 6.04871e5;
65 
66  return (eleca*ekin+elecb) * exp(-ekin/elecc) + (elecd*ekin+elece)* (1.-exp(-ekin/elecc));
67  }
exp
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable STRING $argv[2] set_array QUANTILES set_variable FORMULA *[0] exp(-0.5 *(x-[1])*(x-[1])/([2]*[2]))" set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"` source JAcousticsToolkit.sh typeset -A TRIPODS get_tripods $WORKDIR/tripod.txt TRIPODS XMEAN
Definition: toashort_monitor.sh:46
double JSIRENE::weight_pion ( const double  ekin)
inline

Definition at line 69 of file pythia.hh.

70  {
71  static const double pia = 0.538346;
72  static const double pib = 1.32041;
73  static const double pic = 0.737415;
74  static const double pid = -0.813861;
75  static const double pie = -2.22444;
76  static const double pif = -2.15795;
77  static const double pig = -6.47242;
78  static const double pih = -2.7567;
79  static const double pix = 8.83498;
80 
81  if (ekin < 0.06)
82  return 1e4*pia/ngamma_elec(ekin);
83  else if (ekin < 0.15)
84  return pix*1e4*ekin/ngamma_elec(ekin);
85  else {
86  const double lek=log(ekin);
87  return (pib*1e5*ekin + (pow(ekin,(1.-1./pic))) * (pid*1e4 + 1e4*pie*lek + 1e4*pif*pow(lek,2) + 1e3*pig*pow(lek,3) + 1e2*pih*pow(lek,4)))/ngamma_elec(ekin);
88  }
89  }
log
then set_variable DIR else fatal Wrong number of arguments fi for INPUT_FILE in ls rt $DIR stage * log
Definition: plot-acoustics-fit.sh:30
JSIRENE::ngamma_elec
double ngamma_elec(const double ekin)
Definition: pythia.hh:58
JMATH::pow
T pow(const T &x, const double y)
Power .
Definition: JMath.hh:97
double JSIRENE::weight_kaon ( const double  ekin)
inline

Definition at line 92 of file pythia.hh.

93  {
94  static const double ka = 12.7537;
95  static const double kb = -1.052;
96  static const double kc = 3.49559;
97  static const double kd = 16.7914;
98  static const double ke = -0.355066;
99  static const double kf = 2.77116;
100 
101  if (ekin > 0.26)
102  return (1e4*ka*(ekin+kb)*(1.-exp(-ekin/kc)) + 1e4*(kd*ekin+ke)*exp(-ekin/kc))/ngamma_elec(ekin);
103  else
104  return kf*1e4/ngamma_elec(ekin);
105  }
JSIRENE::ngamma_elec
double ngamma_elec(const double ekin)
Definition: pythia.hh:58
exp
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable STRING $argv[2] set_array QUANTILES set_variable FORMULA *[0] exp(-0.5 *(x-[1])*(x-[1])/([2]*[2]))" set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"` source JAcousticsToolkit.sh typeset -A TRIPODS get_tripods $WORKDIR/tripod.txt TRIPODS XMEAN
Definition: toashort_monitor.sh:46
double JSIRENE::weight_kshort ( const double  ekin)
inline

Definition at line 108 of file pythia.hh.

109  {
110  static const double k0sa = 0.351242;
111  static const double k0sb = 0.613076;
112  static const double k0sc = 6.24682;
113  static const double k0sd = 2.8858;
114 
115  return (k0sa*1e5 + k0sb*1e5*ekin + ekin*k0sc*1e4*log(k0sd + 1./ekin))/ngamma_elec(ekin);
116  }
log
then set_variable DIR else fatal Wrong number of arguments fi for INPUT_FILE in ls rt $DIR stage * log
Definition: plot-acoustics-fit.sh:30
JSIRENE::ngamma_elec
double ngamma_elec(const double ekin)
Definition: pythia.hh:58
double JSIRENE::weight_klong ( const double  ekin)
inline

Definition at line 119 of file pythia.hh.

120  {
121  static const double k0la = 2.18152;
122  static const double k0lb = -0.632798;
123  static const double k0lc = 0.999514;
124  static const double k0ld = 1.36052;
125  static const double k0le = 4.22484;
126  static const double k0lf = -0.307859;
127 
128  if (ekin < 1.5)
129  return (1e4*k0la + ekin*1e5*k0lb*log(ekin) + 1e5*k0lc*pow(ekin,1.5))/ngamma_elec(ekin);
130  else
131  return (ekin*k0ld*1e5 + pow(ekin,1.-1./k0le)*k0lf*1e5)/ngamma_elec(ekin);
132  }
log
then set_variable DIR else fatal Wrong number of arguments fi for INPUT_FILE in ls rt $DIR stage * log
Definition: plot-acoustics-fit.sh:30
JSIRENE::ngamma_elec
double ngamma_elec(const double ekin)
Definition: pythia.hh:58
JMATH::pow
T pow(const T &x, const double y)
Power .
Definition: JMath.hh:97
double JSIRENE::weight_proton ( const double  ekin)
inline

Definition at line 135 of file pythia.hh.

136  {
137  static const double pa = 12.1281;
138  static const double pb = -17.1528;
139  static const double pc = 0.573158;
140  static const double pd = 34.1436;
141  static const double pe = -0.28944;
142  static const double pf = -13.2619;
143  static const double pg = 24.1357;
144 
145  return (1e4*(pa*ekin+pb)*(1.-exp(-ekin/pc)) + 1e4*(pd*ekin +pe+ pf*pow(ekin,2) + pg*pow(ekin,3))*exp(-ekin/pc)) / ngamma_elec(ekin);
146  }
JSIRENE::ngamma_elec
double ngamma_elec(const double ekin)
Definition: pythia.hh:58
JMATH::pow
T pow(const T &x, const double y)
Power .
Definition: JMath.hh:97
exp
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable STRING $argv[2] set_array QUANTILES set_variable FORMULA *[0] exp(-0.5 *(x-[1])*(x-[1])/([2]*[2]))" set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"` source JAcousticsToolkit.sh typeset -A TRIPODS get_tripods $WORKDIR/tripod.txt TRIPODS XMEAN
Definition: toashort_monitor.sh:46
double JSIRENE::weight_neutron ( const double  ekin)
inline

Definition at line 149 of file pythia.hh.

150  {
151  static const double na = 1.24605;
152  static const double nb = 0.63819;
153  static const double nc = -0.802822;
154  static const double nd = -0.935327;
155  static const double ne = -6.1126;
156  static const double nf = -1.96894;
157  static const double ng = 0.716954;
158  static const double nh = 2.68246;
159  static const double ni = -9.39464;
160  static const double nj = 15.0737;
161 
162  if (ekin > 0.5) {
163  const double lek=log(ekin);
164  return (na*1e5*ekin + 1e3*pow(ekin,1.-1./nb) * (100.*nc+100.*nd*lek + 10.*ne*pow(lek,2) + 11.*nf*pow(lek,3)))/ngamma_elec(ekin);
165  } else
166  return (1e3*ng + 1e4*nh*ekin + 1e4*ni*pow(ekin,2) + 1e4*nj*pow(ekin,3))/ngamma_elec(ekin);
167  }
log
then set_variable DIR else fatal Wrong number of arguments fi for INPUT_FILE in ls rt $DIR stage * log
Definition: plot-acoustics-fit.sh:30
JSIRENE::ngamma_elec
double ngamma_elec(const double ekin)
Definition: pythia.hh:58
JMATH::pow
T pow(const T &x, const double y)
Power .
Definition: JMath.hh:97
double JSIRENE::opa_efrac ( const int  ipart,
const double  ekin 
)
inline

Definition at line 173 of file pythia.hh.

174  {
175  using namespace JPP;
176 
177  double etemp, weight_part;
178 
179  // neutrino or muon
180  if (ipart == GEANT4_TYPE_NEUTRINO ||
181  ipart == GEANT4_TYPE_ANTIMUON ||
182  ipart == GEANT4_TYPE_MUON ||
183  ipart == GEANT4_TYPE_ANTITAU ||
184  ipart == GEANT4_TYPE_TAU) {
185  return 0.0;
186  }
187 
188  // 100% energy to EM showers if electron, positron, pi0, or gamma
189  if (ipart == GEANT4_TYPE_PHOTON ||
190  ipart == GEANT4_TYPE_ANTIELECTRON ||
191  ipart == GEANT4_TYPE_ELECTRON ||
192  ipart == GEANT4_TYPE_NEUTRAL_PION ||
193  ipart == GEANT4_TYPE_ETA) {
194  return 1.0;
195  }
196 
197  // low-energy fits went up to 40 GeV
198  if (ekin <= 40)
199  etemp=ekin;
200  else
201  etemp=40.;
202 
203  // otherwise, calculate fractions the long way. First we get the
204  // expected number of photons from the particle:
205  if (ipart == GEANT4_TYPE_PION_PLUS ||
206  ipart == GEANT4_TYPE_PION_MINUS)
207  weight_part=weight_pion(etemp);
208  else if (ipart == GEANT4_TYPE_KAON_LONG)
209  weight_part=weight_klong(etemp);
210  else if (ipart == GEANT4_TYPE_KAON_SHORT)
211  weight_part=weight_kshort(etemp);
212  else if (ipart == GEANT4_TYPE_KAON_PLUS ||
213  ipart == GEANT4_TYPE_KAON_MINUS)
214  weight_part=weight_kaon(etemp);
215  else if (ipart == GEANT4_TYPE_NEUTRON)
216  weight_part=weight_neutron(etemp);
217  else if (ipart == GEANT4_TYPE_PROTON ||
218  ipart == GEANT4_TYPE_ANTIPROTON)
219  weight_part=weight_proton(etemp);
220  else
221  // when in doubt, treat it as a proton
222  weight_part=weight_proton(etemp);
223 
224  if (ekin <= 40)
225  return weight_part;
226  else if (ekin >= 1e7)
227  return opa_weight_high_e(ekin);
228  else {
229  // smooth transition between low-energy weights and high-energy pion weight
230  const double wp40 = weight_part;
231  const double whe40 = opa_weight_high_e(40.);
232  const double whex = opa_weight_high_e(ekin);
233  return whex - (whe40-wp40)*(7.-log10(ekin))/5.398;
234  }
235  }
JSIRENE::weight_klong
double weight_klong(const double ekin)
Definition: pythia.hh:119
JSIRENE::opa_weight_high_e
double opa_weight_high_e(const double ekin)
Definition: pythia.hh:20
JSIRENE::weight_proton
double weight_proton(const double ekin)
Definition: pythia.hh:135
JSIRENE::weight_pion
double weight_pion(const double ekin)
Definition: pythia.hh:69
JSIRENE::weight_kaon
double weight_kaon(const double ekin)
Definition: pythia.hh:92
log10
set_variable E_E log10(E_{fit}/E_{#mu})"
JSIRENE::weight_kshort
double weight_kshort(const double ekin)
Definition: pythia.hh:108
JSIRENE::weight_neutron
double weight_neutron(const double ekin)
Definition: pythia.hh:149
double JSIRENE::pythia ( const int  type,
const double  E 
)
inline

Get equivalent EM-energy for given pion energy.

Parameters
typeparticle type [PDG]
Eparticle energy [GeV]
Returns
EM-equivalent energy [GeV]

Definition at line 245 of file pythia.hh.

246  {
247  int ipart = JPDB::getInstance().getPDG(type).geant;
248  float ekin = (float) E;
249 
250  return opa_efrac(ipart, ekin) * E;
251  }
E
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
Definition: JMuonPostfit.sh:40
JLANG::getInstance
T & getInstance(const T &object)
Get static instance from temporary object.
Definition: JObject.hh:75
JSIRENE::opa_efrac
double opa_efrac(const int ipart, const double ekin)
Definition: pythia.hh:173

Variable Documentation

const JPythia JSIRENE::pythia
static

Function object for relative light yield as a function of GEANT particle code.

Definition at line 96 of file JPythia.hh.

std::vector<size_t> JSIRENE::ns
mutable

Definition at line 138 of file JSireneToolkit.hh.

std::vector<double> JSIRENE::vs
mutable

Definition at line 139 of file JSireneToolkit.hh.

struct { ... } JSIRENE::getNumberOfPhotoElectrons

Auxiliary data structure for determination of number of photo-electrons.

Generated by   doxygen 1.8.5