IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 0a262e3d authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan: Committed by ralfulrich
Browse files

examples

parent e039c47a
No related branches found
No related tags found
1 merge request!280Refactory 2020
...@@ -46,7 +46,7 @@ struct MyBoundaryCrossingProcess ...@@ -46,7 +46,7 @@ struct MyBoundaryCrossingProcess
typename Particle::BaseNodeType const& to) { typename Particle::BaseNodeType const& to) {
std::cout << "boundary crossing! from: " << &from << "; to: " << &to << std::endl; std::cout << "boundary crossing! from: " << &from << "; to: " << &to << std::endl;
auto const& name = corsika::GetName(p.GetPID()); auto const& name = corsika::name(p.GetPID());
auto const start = p.GetPosition().GetCoordinates(); auto const start = p.GetPosition().GetCoordinates();
fFile << name << " " << start[0] / 1_m << ' ' << start[1] / 1_m << ' ' fFile << name << " " << start[0] / 1_m << ' ' << start[1] / 1_m << ' '
...@@ -116,7 +116,7 @@ int main() { ...@@ -116,7 +116,7 @@ int main() {
setup::Stack stack; setup::Stack stack;
stack.Clear(); stack.Clear();
const Code beamCode = Code::Proton; const Code beamCode = Code::Proton;
const HEPMassType mass = corsika::GetMass(Code::Proton); const HEPMassType mass = corsika::mass(Code::Proton);
const HEPEnergyType E0 = 50_TeV; const HEPEnergyType E0 = 50_TeV;
std::uniform_real_distribution distTheta(0., 180.); std::uniform_real_distribution distTheta(0., 180.);
......
...@@ -81,7 +81,7 @@ int main() { ...@@ -81,7 +81,7 @@ int main() {
const Code beamCode = Code::Nucleus; const Code beamCode = Code::Nucleus;
const int nuclA = 4; const int nuclA = 4;
const int nuclZ = int(nuclA / 2.15 + 0.7); const int nuclZ = int(nuclA / 2.15 + 0.7);
const HEPMassType mass = GetNucleusMass(nuclA, nuclZ); const HEPMassType mass = nucleus_mass(nuclA, nuclZ);
const HEPEnergyType E0 = nuclA * 1_TeV; const HEPEnergyType E0 = nuclA * 1_TeV;
double theta = 0.; double theta = 0.;
double phi = 0.; double phi = 0.;
......
...@@ -69,7 +69,7 @@ int main() { ...@@ -69,7 +69,7 @@ int main() {
setup::Stack stack; setup::Stack stack;
stack.Clear(); stack.Clear();
const Code beamCode = Code::Proton; const Code beamCode = Code::Proton;
const HEPMassType mass = corsika::Proton::GetMass(); const HEPMassType mass = corsika::mass(beamCode);
const HEPEnergyType E0 = 100_GeV; const HEPEnergyType E0 = 100_GeV;
double theta = 0.; double theta = 0.;
double phi = 0.; double phi = 0.;
......
...@@ -52,7 +52,7 @@ int main() { ...@@ -52,7 +52,7 @@ int main() {
for (HEPEnergyType E0 = 300_MeV; E0 < 1_PeV; E0 *= 1.05) { for (HEPEnergyType E0 = 300_MeV; E0 < 1_PeV; E0 *= 1.05) {
stack.Clear(); stack.Clear();
const Code beamCode = Code::MuPlus; const Code beamCode = Code::MuPlus;
const HEPMassType mass = GetMass(beamCode); const HEPMassType mass = corsika::mass(beamCode);
double theta = 0.; double theta = 0.;
double phi = 0.; double phi = 0.;
......
...@@ -77,7 +77,7 @@ int main() { ...@@ -77,7 +77,7 @@ int main() {
setup::Stack stack; setup::Stack stack;
stack.Clear(); stack.Clear();
const Code beamCode = Code::Proton; const Code beamCode = Code::Proton;
auto const mass = corsika::GetMass(beamCode); auto const mass = corsika::mass(beamCode);
const HEPEnergyType E0 = 0.1_PeV; const HEPEnergyType E0 = 0.1_PeV;
double theta = 0.; double theta = 0.;
double phi = 0.; double phi = 0.;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment