| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 //---------------------------------------------------------------------------
28 //
29 // ClassName: G4NeutronTrackingCut
30 //
31 // Author: Nov 2006 G.Folger
32 //
33 //
34 //----------------------------------------------------------------------------
35 //
36
37 #include "G4NeutronTrackingCut.hh"
38
39 #include "G4SystemOfUnits.hh"
40 #include "G4ParticleDefinition.hh"
41 #include "G4ProcessManager.hh"
42
43 #include "G4Neutron.hh"
44 #include "G4NeutronKiller.hh"
45 #include "G4HadronicProcessStore.hh"
46 #include "G4PhysListUtil.hh"
47 #include "G4NeutronGeneralProcess.hh"
48
49 #include "G4BuilderType.hh"
50 #include "G4Threading.hh"
51
52 // factory
53 #include "G4PhysicsConstructorFactory.hh"
54 //
55 G4_DECLARE_PHYSCONSTR_FACTORY(G4NeutronTrackingCut);
56 //
57
58 G4NeutronTrackingCut::G4NeutronTrackingCut(G4int ver)
59 : G4VPhysicsConstructor("neutronTrackingCut"), verbose(ver)
60 {
61 timeLimit = 10*CLHEP::microsecond;
62 kineticEnergyLimit = 0.0;
63 SetPhysicsType(bUnknown);
64 }
65
66 G4NeutronTrackingCut::G4NeutronTrackingCut(const G4String&, G4int ver)
67 : G4NeutronTrackingCut(ver)
68 {}
69
70 G4NeutronTrackingCut::~G4NeutronTrackingCut()
71 {}
72
73 void G4NeutronTrackingCut::ConstructParticle()
74 {
75 G4Neutron::NeutronDefinition();
76 }
77
78 void G4NeutronTrackingCut::ConstructProcess()
79 {
80 auto particle = G4Neutron::Neutron();
81 auto proc = dynamic_cast<G4NeutronGeneralProcess*>(G4PhysListUtil::FindProcess(particle, fNeutronGeneral));
82 if(nullptr != proc) { return; }
83
84 G4NeutronKiller* pNeutronKiller = new G4NeutronKiller();
85 G4ProcessManager * pmanager = particle->GetProcessManager();
86
87 if(verbose > 0 && G4Threading::IsMasterThread()) {
88 G4String pn = particle->GetParticleName();//Avoid data-race when passing
89 //this string to G4MTcout
90 G4cout << "### Adding tracking cuts for " << pn
91 << " TimeCut(ns)= " << timeLimit/ns
92 << " KinEnergyCut(MeV)= " << kineticEnergyLimit/MeV
93 << G4endl;
94 }
95 pmanager->AddDiscreteProcess(pNeutronKiller);
96 pNeutronKiller->SetKinEnergyLimit(kineticEnergyLimit);
97 pNeutronKiller->SetTimeLimit(timeLimit);
98
99 G4HadronicProcessStore::Instance()->RegisterExtraProcess(pNeutronKiller);
100 G4HadronicProcessStore::Instance()->RegisterParticleForExtraProcess(pNeutronKiller, particle);
101 }
102
103
104