IAP GITLAB

Skip to content
Snippets Groups Projects
Commit b2acc431 authored by ralfulrich's avatar ralfulrich
Browse files

remove ObservatioLevel protoype process

parent 21217662
No related branches found
No related tags found
No related merge requests found
...@@ -207,30 +207,6 @@ public: ...@@ -207,30 +207,6 @@ public:
HEPEnergyType GetEmEnergy() const { return fEmEnergy; } 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 // The example main program for a particle cascade
...@@ -286,7 +262,6 @@ int main() { ...@@ -286,7 +262,6 @@ int main() {
process::sibyll::NuclearInteraction sibyllNuc(sibyll, env); process::sibyll::NuclearInteraction sibyllNuc(sibyll, env);
process::sibyll::Decay decay(trackedHadrons); process::sibyll::Decay decay(trackedHadrons);
ProcessCut cut(20_GeV); ProcessCut cut(20_GeV);
ObservationLevel obsLevel(height_atmosphere - 2000_m); // 1400_m);
process::TrackWriter::TrackWriter trackWriter("tracks.dat"); process::TrackWriter::TrackWriter trackWriter("tracks.dat");
process::EnergyLoss::EnergyLoss eLoss; process::EnergyLoss::EnergyLoss eLoss;
......
...@@ -210,30 +210,6 @@ public: ...@@ -210,30 +210,6 @@ public:
HEPEnergyType GetEmEnergy() const { return fEmEnergy; } 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 // The example main program for a particle cascade
...@@ -280,7 +256,6 @@ int main() { ...@@ -280,7 +256,6 @@ int main() {
// random::RNGManager::GetInstance().RegisterRandomStream("pythia"); // random::RNGManager::GetInstance().RegisterRandomStream("pythia");
// process::pythia::Decay decay(trackedHadrons); // process::pythia::Decay decay(trackedHadrons);
ProcessCut cut(20_GeV); ProcessCut cut(20_GeV);
ObservationLevel obsLevel(1400_m);
// random::RNGManager::GetInstance().RegisterRandomStream("HadronicElasticModel"); // random::RNGManager::GetInstance().RegisterRandomStream("HadronicElasticModel");
// process::HadronicElasticModel::HadronicElasticInteraction // process::HadronicElasticModel::HadronicElasticInteraction
...@@ -291,7 +266,7 @@ int main() { ...@@ -291,7 +266,7 @@ int main() {
// assemble all processes into an ordered process list // assemble all processes into an ordered process list
// cut << trackWriter; // 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() // cout << "decltype(sequence)=" << type_id_with_cvr<decltype(sequence)>().pretty_name()
// << "\n"; // << "\n";
......
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