BDSIM
BDSIM is a Geant4 extension toolkit for simulation of particle transport in accelerator beamlines.
Loading...
Searching...
No Matches
BDSProcessUserSpecialCutsPDG.cc
1/*
2Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway,
3University of London 2001 - 2023.
4
5This file is part of BDSIM.
6
7BDSIM is free software: you can redistribute it and/or modify
8it under the terms of the GNU General Public License as published
9by the Free Software Foundation version 3 of the License.
10
11BDSIM is distributed in the hope that it will be useful, but
12WITHOUT ANY WARRANTY; without even the implied warranty of
13MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14GNU General Public License for more details.
15
16You should have received a copy of the GNU General Public License
17along with BDSIM. If not, see <http://www.gnu.org/licenses/>.
18*/
19#include "BDSProcessUserSpecialCutsPDG.hh"
20
21#include "G4Types.hh"
22
23#include <limits>
24#include <set>
25
26BDSProcessUserSpecialCutsPDG::BDSProcessUserSpecialCutsPDG():
27 G4UserSpecialCuts("BDSUserSpecialCutsPDGID")
28{;}
29
30BDSProcessUserSpecialCutsPDG::BDSProcessUserSpecialCutsPDG(const std::set<G4int>& pdgIDsIn):
31 G4UserSpecialCuts("BDSUserSpecialCutsPDGID"),
32 pdgIDs(pdgIDsIn)
33{;}
34
35BDSProcessUserSpecialCutsPDG::~BDSProcessUserSpecialCutsPDG()
36{;}
37
38G4double BDSProcessUserSpecialCutsPDG::PostStepGetPhysicalInteractionLength(const G4Track& track,
39 G4double previousStepSize,
40 G4ForceCondition* condition)
41{
42 G4int pdgID = track.GetDynamicParticle()->GetPDGcode();
43 if (pdgIDs.count(pdgID) > 0)
44 {
45 *condition = NotForced;
46 return std::numeric_limits<double>::max();
47 }
48 else
49 {return G4UserSpecialCuts::PostStepGetPhysicalInteractionLength(track, previousStepSize, condition);}
50}