| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // INCL++ intra-nuclear cascade model 26 // INCL++ intra-nuclear cascade model
27 // Alain Boudard, CEA-Saclay, France << 27 // Pekka Kaitaniemi, CEA and Helsinki Institute of Physics
28 // Joseph Cugnon, University of Liege, Belgium << 28 // Davide Mancusi, CEA
29 // Jean-Christophe David, CEA-Saclay, France << 29 // Alain Boudard, CEA
30 // Pekka Kaitaniemi, CEA-Saclay, France, and H << 30 // Sylvie Leray, CEA
31 // Sylvie Leray, CEA-Saclay, France << 31 // Joseph Cugnon, University of Liege
32 // Davide Mancusi, CEA-Saclay, France <<
33 // 32 //
34 #define INCLXX_IN_GEANT4_MODE 1 33 #define INCLXX_IN_GEANT4_MODE 1
35 34
36 #include "globals.hh" 35 #include "globals.hh"
37 36
38 #include "G4INCLRandom.hh" 37 #include "G4INCLRandom.hh"
39 #include "G4INCLPauliBlocking.hh" 38 #include "G4INCLPauliBlocking.hh"
40 #include "G4INCLGlobals.hh" 39 #include "G4INCLGlobals.hh"
41 40
42 #include "G4INCLPauliStrict.hh" <<
43 #include "G4INCLPauliStandard.hh" <<
44 #include "G4INCLPauliStrictStandard.hh" <<
45 #include "G4INCLPauliGlobal.hh" <<
46 #include "G4INCLCDPP.hh" <<
47 <<
48 namespace G4INCL { 41 namespace G4INCL {
49 42
50 namespace Pauli { << 43 G4INCL::IPauli const * Pauli::thePauliBlocker = 0;
51 << 44 G4INCL::IPauli const * Pauli::theCDPP = 0;
52 namespace { <<
53 G4ThreadLocal IPauli * thePauliBlocker = <<
54 G4ThreadLocal IPauli * theCDPP = NULL; <<
55 } <<
56 <<
57 IPauli * getBlocker() { return thePauliBlo <<
58 45
59 IPauli * getCDPP() { return theCDPP; } << 46 void Pauli::setBlocker(IPauli const * pauliBlocker) {
60 << 47 thePauliBlocker = pauliBlocker;
61 void setBlocker(IPauli * const pauliBlocke << 48 }
62 thePauliBlocker = pauliBlocker; <<
63 } <<
64 <<
65 void setCDPP(IPauli * const cdpp) { <<
66 theCDPP = cdpp; <<
67 } <<
68 49
69 G4bool isBlocked(ParticleList const &modif << 50 void Pauli::setCDPP(IPauli const * cdpp) {
70 G4bool isPauliBlocked = false; << 51 theCDPP = cdpp;
71 if(thePauliBlocker != 0) { << 52 }
72 isPauliBlocked = thePauliBlocker->isBl <<
73 } <<
74 53
75 return isPauliBlocked; << 54 G4bool Pauli::isBlocked(ParticleList const modifiedAndCreated, Nucleus const * const nucleus) {
>> 55 G4bool isPauliBlocked = false;
>> 56 if(thePauliBlocker != 0) {
>> 57 isPauliBlocked = thePauliBlocker->isBlocked(modifiedAndCreated, nucleus);
76 } 58 }
77 59
78 G4bool isCDPPBlocked(ParticleList const &c << 60 return isPauliBlocked;
79 G4bool isCDPPBlocked = false; << 61 }
80 if(theCDPP != 0) { <<
81 isCDPPBlocked = theCDPP->isBlocked(cre <<
82 } <<
83 <<
84 return isCDPPBlocked; <<
85 } <<
86 62
87 void deleteBlockers() { << 63 G4bool Pauli::isCDPPBlocked(ParticleList const created, Nucleus const * const nucleus) {
88 delete thePauliBlocker; << 64 G4bool isCDPPBlocked = false;
89 thePauliBlocker=NULL; << 65 if(theCDPP != 0) {
90 delete theCDPP; << 66 isCDPPBlocked = theCDPP->isBlocked(created, nucleus);
91 theCDPP=NULL; <<
92 } 67 }
93 68
94 void initialize(Config const * const aConf << 69 return isCDPPBlocked;
95 // Select the Pauli blocking algorithm: <<
96 PauliType pauli = aConfig->getPauliType( <<
97 if(pauli == StrictStatisticalPauli) <<
98 setBlocker(new PauliStrictStandard); <<
99 else if(pauli == StatisticalPauli) <<
100 setBlocker(new PauliStandard); <<
101 else if(pauli == StrictPauli) <<
102 setBlocker(new PauliStrict); <<
103 else if(pauli == GlobalPauli) <<
104 setBlocker(new PauliGlobal); <<
105 else if(pauli == NoPauli) <<
106 setBlocker(NULL); <<
107 <<
108 if(aConfig->getCDPP()) <<
109 setCDPP(new CDPP); <<
110 else <<
111 setCDPP(NULL); <<
112 <<
113 } <<
114 } 70 }
>> 71
115 } 72 }
116 73