IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 5ddaccd3 authored by Maximilian Sackel's avatar Maximilian Sackel Committed by Ralf Ulrich
Browse files

try to modify cascade example so that proposal is used.

parent 4e190d56
No related branches found
No related tags found
1 merge request!245Include proposal process rebase
......@@ -137,6 +137,7 @@ int main() {
random::RNGManager::GetInstance().RegisterRandomStream("sibyll");
random::RNGManager::GetInstance().RegisterRandomStream("pythia");
random::RNGManager::GetInstance().RegisterRandomStream("proposal");
process::sibyll::Interaction sibyll;
process::sibyll::NuclearInteraction sibyllNuc(sibyll, env);
process::sibyll::Decay decay;
......@@ -148,21 +149,22 @@ int main() {
process::energy_loss::EnergyLoss eLoss{showerAxis};
// assemble all processes into an ordered process list
auto sequence = stackInspect << sibyll << sibyllNuc << decay << eLoss << cut
<< trackWriter;
auto sequence = stackInspect << sibyll << sibyllNuc << proposal << decay
/* << eLoss */
<< cut << trackWriter;
// define air shower object, run simulation
cascade::Cascade EAS(env, tracking, sequence, stack);
EAS.Run();
eLoss.PrintProfile(); // print longitudinal profile
/* eLoss.PrintProfile(); // print longitudinal profile */
cut.ShowResults();
const HEPEnergyType Efinal =
cut.GetCutEnergy() + cut.GetInvEnergy() + cut.GetEmEnergy();
cout << "total cut energy (GeV): " << Efinal / 1_GeV << endl
<< "relative difference (%): " << (Efinal / E0 - 1) * 100 << endl;
cout << "total dEdX energy (GeV): " << eLoss.GetTotal() / 1_GeV << endl
<< "relative difference (%): " << eLoss.GetTotal() / E0 * 100 << endl;
/* cout << "total dEdX energy (GeV): " << eLoss.GetTotal() / 1_GeV << endl */
/* << "relative difference (%): " << eLoss.GetTotal() / E0 * 100 << endl; */
}
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