| 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 /// \file Par01/src/Par01PiModel.cc 26 /// \file Par01/src/Par01PiModel.cc
27 /// \brief Implementation of the Par01PiModel 27 /// \brief Implementation of the Par01PiModel class
28 // 28 //
29 // 29 //
30 // 30 //
31 #include "Par01PiModel.hh" 31 #include "Par01PiModel.hh"
32 32
33 #include "G4Gamma.hh" <<
34 #include "G4PionMinus.hh" 33 #include "G4PionMinus.hh"
35 #include "G4PionPlus.hh" 34 #include "G4PionPlus.hh"
>> 35 #include "G4Gamma.hh"
36 36
37 //....oooOO0OOooo........oooOO0OOooo........oo 37 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
38 38
39 Par01PiModel::Par01PiModel(G4Region* anEnvelop << 39 Par01PiModel::Par01PiModel(G4Region *anEnvelope) :
40 : G4VFastSimulationModel("Par01PiModel", anE << 40 G4VFastSimulationModel("Par01PiModel",anEnvelope)
41 { << 41 {;}
42 ; <<
43 } <<
44 42
45 //....oooOO0OOooo........oooOO0OOooo........oo 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
46 44
47 Par01PiModel::~Par01PiModel() 45 Par01PiModel::~Par01PiModel()
48 { << 46 {;}
49 ; <<
50 } <<
51 47
52 //....oooOO0OOooo........oooOO0OOooo........oo 48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
53 49
54 G4bool Par01PiModel::IsApplicable(const G4Part 50 G4bool Par01PiModel::IsApplicable(const G4ParticleDefinition& particleType)
55 { 51 {
56 return &particleType == G4PionMinus::PionMin << 52 return
57 || &particleType == G4PionPlus::PionP << 53 &particleType == G4PionMinus::PionMinusDefinition() ||
>> 54 &particleType == G4PionPlus::PionPlusDefinition();
58 } 55 }
59 56
60 //....oooOO0OOooo........oooOO0OOooo........oo 57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
61 58
62 G4bool Par01PiModel::ModelTrigger(const G4Fast << 59 G4bool Par01PiModel::ModelTrigger(const G4FastTrack& fastTrack) {
63 { <<
64 //------------------------------------------ 60 //-------------------------------------------------------------
65 // UserTrigger() method: method which has to 61 // UserTrigger() method: method which has to decide if
66 // the parameterisation has to be applied. 62 // the parameterisation has to be applied.
67 // Here ModelTrigger() asks the user (ie you 63 // Here ModelTrigger() asks the user (ie you) a 0/1 answer.
68 // 64 //
69 // Note that quantities like the local/globa 65 // Note that quantities like the local/global position/direction etc..
70 // are available at this level via the fastT << 66 // are available at this level via the fastTrack parameter (allowing
71 // to check distance from boundaries, see be 67 // to check distance from boundaries, see below to allow the decision)
72 //------------------------------------------ 68 //--------------------------------------------------------------
73 G4cout << "\nPar01PiModel::ModelTrigger() ca 69 G4cout << "\nPar01PiModel::ModelTrigger() called:" << G4endl;
74 G4cout << "--------------------------------" << 70 G4cout << "--------------------------------" << G4endl;
75 G4cout << "(particle is a " << fastTrack.Get << 71 G4cout << "(particle is a " << fastTrack.GetPrimaryTrack()->
76 << " )\n" << 72 GetDefinition()->GetParticleName() << " )\n" << G4endl;
77 << G4endl; <<
78 73
79 // -- Examples of available informations: 74 // -- Examples of available informations:
80 75
81 // -- position: 76 // -- position:
82 G4cout << " Track position: " << fas << 77 G4cout << " Track position: " <<
83 << "(global coord.)" << fastTrack.Get << 78 fastTrack.GetPrimaryTrack()->GetPosition() << "(global coord.)" <<
84 << G4endl; << 79 fastTrack.GetPrimaryTrackLocalPosition() << "(in envelope coord.)"
>> 80 << G4endl;
85 81
86 // -- direction: 82 // -- direction:
87 G4cout << " Track direction:" << fas << 83 G4cout << " Track direction:" <<
88 << "(global coord.)" << fastTrack.Get << 84 fastTrack.GetPrimaryTrack()->GetMomentum().unit() <<
89 << G4endl; << 85 "(global coord.)" <<
>> 86 fastTrack.GetPrimaryTrackLocalDirection() << "(in envelope coord.)" <<
>> 87 G4endl;
90 88
91 return true; 89 return true;
92 } 90 }
93 91
94 //....oooOO0OOooo........oooOO0OOooo........oo 92 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
95 93
96 void Par01PiModel::DoIt(const G4FastTrack& fas << 94 void Par01PiModel::DoIt(const G4FastTrack& fastTrack,
97 //-------------------------------------------- << 95 G4FastStep& fastStep)
98 // << 96 //--------------------------------------------------------
99 // User method to code the parameterisation pr << 97 //
100 // << 98 // User method to code the parameterisation properly said.
101 //-------------------------------------------- << 99 //
>> 100 //--------------------------------------------------------
102 { 101 {
>> 102
103 //------------------------------------------ 103 //------------------------------------------------
104 // The primary track continues along its dir 104 // The primary track continues along its direction.
105 // One secondary (a photon) is added: 105 // One secondary (a photon) is added:
106 //------------------------------------------ 106 //------------------------------------------------
107 G4cout << " Pion `model' applied\n" << 107 G4cout << " Pion `model' applied\n" << G4endl;
108 108
109 //------------------------------ 109 //------------------------------
110 // Primary: 110 // Primary:
111 // idem as in "DefaultModel": 111 // idem as in "DefaultModel":
112 // 112 //
113 //------------------------------ 113 //------------------------------
114 G4ThreeVector position; 114 G4ThreeVector position;
115 G4double distance; 115 G4double distance;
116 distance = fastTrack.GetEnvelopeSolid()->Dis << 116 distance = fastTrack.GetEnvelopeSolid()->
117 << 117 DistanceToOut(fastTrack.GetPrimaryTrackLocalPosition(),
118 position = << 118 fastTrack.GetPrimaryTrackLocalDirection());
119 fastTrack.GetPrimaryTrackLocalPosition() + << 119 position = fastTrack.GetPrimaryTrackLocalPosition() +
120 << 120 distance*fastTrack.GetPrimaryTrackLocalDirection();
>> 121
121 // -- set final position: 122 // -- set final position:
122 fastStep.ProposePrimaryTrackFinalPosition(po 123 fastStep.ProposePrimaryTrackFinalPosition(position);
123 << 124
124 //--------------------------- 125 //---------------------------
125 // Secondary: 126 // Secondary:
126 // Adds one "secondary": 127 // Adds one "secondary":
127 //--------------------------- 128 //---------------------------
128 // -- First, user has to say how many second 129 // -- First, user has to say how many secondaries will be created:
129 fastStep.SetNumberOfSecondaryTracks(1); 130 fastStep.SetNumberOfSecondaryTracks(1);
130 << 131
131 //------------------------ 132 //------------------------
132 // -- Build the secondary: 133 // -- Build the secondary:
133 //------------------------ 134 //------------------------
134 // -- direction: 135 // -- direction:
135 G4ParticleMomentum direction(fastTrack.GetPr 136 G4ParticleMomentum direction(fastTrack.GetPrimaryTrackLocalDirection());
136 direction.setZ(direction.z() * 0.5); << 137 direction.setZ(direction.z()*0.5);
137 direction.setY(direction.y() + direction.z() << 138 direction.setY(direction.y()+direction.z()*0.1);
138 direction = direction.unit(); // necessary << 139 direction = direction.unit(); // necessary ?
139 << 140
140 // -- dynamics (Note that many constructors 141 // -- dynamics (Note that many constructors exists for G4DynamicParticle
141 // -- see prototype/particle+matter/particle 142 // -- see prototype/particle+matter/particles/management/include/G4DynamicParticle.hh)
142 G4DynamicParticle dynamique(G4Gamma::GammaDe << 143 G4DynamicParticle dynamique(G4Gamma::GammaDefinition(),
143 fastTrack.GetPri << 144 direction,
>> 145 fastTrack.GetPrimaryTrack()->
>> 146 GetKineticEnergy()/2.);
144 // -- position: 147 // -- position:
145 G4double Dist; 148 G4double Dist;
146 Dist = fastTrack.GetEnvelopeSolid()->Distanc << 149 Dist = fastTrack.GetEnvelopeSolid()->
147 << 150 DistanceToOut(fastTrack.GetPrimaryTrackLocalPosition(),
>> 151 direction);
148 G4ThreeVector posi; 152 G4ThreeVector posi;
149 posi = fastTrack.GetPrimaryTrackLocalPositio << 153 posi = fastTrack.GetPrimaryTrackLocalPosition() + Dist*direction;
150 << 154
151 //------------------------------------ 155 //------------------------------------
152 //-- Creation of the secondary Track: 156 //-- Creation of the secondary Track:
153 //------------------------------------ 157 //------------------------------------
154 fastStep.CreateSecondaryTrack(dynamique, pos << 158 fastStep.CreateSecondaryTrack(dynamique, posi,
>> 159 fastTrack.GetPrimaryTrack()->GetGlobalTime());
155 } 160 }
156 161