diff --git a/Documentation/Examples/cascade_example.cc b/Documentation/Examples/cascade_example.cc index 269ee60c372e79b471510601fd8d8459e06946f8..03feb8f6eb88b8d5f023bbcaec7ede771e9fadfd 100644 --- a/Documentation/Examples/cascade_example.cc +++ b/Documentation/Examples/cascade_example.cc @@ -207,30 +207,6 @@ public: HEPEnergyType GetEmEnergy() const { return fEmEnergy; } }; -class ObservationLevel : public process::ContinuousProcess<ObservationLevel> { - - LengthType fHeight; - -public: - ObservationLevel(const LengthType vHeight) - : fHeight(vHeight) {} - - template <typename Particle> - LengthType MaxStepLength(Particle&, setup::Trajectory&) const { - return 1_m * std::numeric_limits<double>::infinity(); - } - - template <typename TParticle, typename TTrack> - EProcessReturn DoContinuous(TParticle&, TTrack& vT) { - if ((vT.GetPosition(0).GetZ() <= fHeight && vT.GetPosition(1).GetZ() > fHeight) || - (vT.GetPosition(0).GetZ() > fHeight && vT.GetPosition(1).GetZ() <= fHeight)) { - cout << "OBSERVED " << endl; - return EProcessReturn::eParticleAbsorbed; - } - return EProcessReturn::eOk; - } - void Init() {} -}; // // The example main program for a particle cascade @@ -286,7 +262,6 @@ int main() { process::sibyll::NuclearInteraction sibyllNuc(sibyll, env); process::sibyll::Decay decay(trackedHadrons); ProcessCut cut(20_GeV); - ObservationLevel obsLevel(height_atmosphere - 2000_m); // 1400_m); process::TrackWriter::TrackWriter trackWriter("tracks.dat"); process::EnergyLoss::EnergyLoss eLoss; diff --git a/Documentation/Examples/vertical_EAS.cc b/Documentation/Examples/vertical_EAS.cc index efeb157b75b066746aa0e7e0c1ae2506739e866a..e2d381914668b83431cccac1c2a51caa117ed22f 100644 --- a/Documentation/Examples/vertical_EAS.cc +++ b/Documentation/Examples/vertical_EAS.cc @@ -210,30 +210,6 @@ public: HEPEnergyType GetEmEnergy() const { return fEmEnergy; } }; -class ObservationLevel : public process::ContinuousProcess<ObservationLevel> { - - LengthType fHeight; - -public: - ObservationLevel(const LengthType vHeight) - : fHeight(vHeight) {} - - template <typename Particle> - LengthType MaxStepLength(Particle&, setup::Trajectory&) const { - return 1_m * std::numeric_limits<double>::infinity(); - } - - template <typename TParticle, typename TTrack> - EProcessReturn DoContinuous(TParticle&, TTrack& vT) { - if ((vT.GetPosition(0).GetZ() <= fHeight && vT.GetPosition(1).GetZ() > fHeight) || - (vT.GetPosition(0).GetZ() > fHeight && vT.GetPosition(1).GetZ() <= fHeight)) { - cout << "OBSERVED " << endl; - return EProcessReturn::eParticleAbsorbed; - } - return EProcessReturn::eOk; - } - void Init() {} -}; // // The example main program for a particle cascade @@ -280,7 +256,6 @@ int main() { // random::RNGManager::GetInstance().RegisterRandomStream("pythia"); // process::pythia::Decay decay(trackedHadrons); ProcessCut cut(20_GeV); - ObservationLevel obsLevel(1400_m); // random::RNGManager::GetInstance().RegisterRandomStream("HadronicElasticModel"); // process::HadronicElasticModel::HadronicElasticInteraction @@ -291,7 +266,7 @@ int main() { // assemble all processes into an ordered process list // cut << trackWriter; - auto sequence = sibyll << sibyllNuc << decay << eLoss << cut << obsLevel; + auto sequence = sibyll << sibyllNuc << decay << eLoss << cut; // cout << "decltype(sequence)=" << type_id_with_cvr<decltype(sequence)>().pretty_name() // << "\n";