ECCE @ EIC Software
Reference for
ECCE @ EIC
simulation and reconstruction software on GitHub
Home page
Related Pages
Modules
Namespaces
Classes
Files
External Links
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
G4UCNMultiScattering.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file G4UCNMultiScattering.cc
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
//
28
// UCN Multiple Scattering Class Implementation
30
//
31
// File: G4UCNMultiScattering.cc
32
// Description: G4VDiscreteProcess -- MultiScattering of UCNs
33
// Version: 1.0
34
// Created: 2014-05-12
35
// Author: Peter Gumplinger
36
// adopted from Geant4UCN by Peter Fierlinger (7.9.04) and
37
// Marcin Kuzniak (21.4.06)
38
// Calculate "elastic scattering" inside materials
39
// Updated:
40
//
41
// mail: gum@triumf.ca
42
//
44
45
#include "
G4UCNProcessSubType.hh
"
46
47
#include "
G4UCNMultiScattering.hh
"
48
49
#include "
G4SystemOfUnits.hh
"
50
#include "
G4PhysicalConstants.hh
"
51
53
// Class Implementation
55
57
// Operators
59
60
// G4UCNMultiScattering::operator=(const G4UCNMultiScattering &right)
61
// {
62
// }
63
65
// Constructors
67
68
G4UCNMultiScattering::G4UCNMultiScattering
(
const
G4String
& processName,
69
G4ProcessType
type)
70
:
G4VDiscreteProcess
(processName, type)
71
{
72
if
(
verboseLevel
>0)
G4cout
<<
GetProcessName
() <<
" is created "
<<
G4endl
;
73
74
SetProcessSubType
(
fUCNMultiScattering
);
75
}
76
77
// G4UCNMultiScattering::G4UCNMultiScattering(const G4UCNMultiScattering &right)
78
// {
79
// }
80
82
// Destructors
84
85
G4UCNMultiScattering::~G4UCNMultiScattering
(){}
86
88
// Methods
90
91
// PostStepDoIt
92
// ------------
93
94
G4VParticleChange
*
95
G4UCNMultiScattering::PostStepDoIt
(
const
G4Track
& aTrack,
const
G4Step
& aStep)
96
{
97
aParticleChange
.
Initialize
(aTrack);
98
99
if
(
verboseLevel
> 0 )
G4cout
<<
"UCNMULTISCATTER at: "
100
<< aTrack.
GetProperTime
()/
s
<<
"s, "
101
<< aTrack.
GetGlobalTime
()/
s
<<
"s. "
102
<<
", after track length "
<< aTrack.
GetTrackLength
()/
cm
<<
"cm, "
103
<<
"in volume "
104
<< aStep.
GetPostStepPoint
()->
GetPhysicalVolume
()->
GetName
()
105
<<
G4endl
;
106
107
G4ThreeVector
scattered =
Scatter
();
108
109
aParticleChange
.
ProposeMomentumDirection
(-scattered);
110
111
return
G4VDiscreteProcess::PostStepDoIt
(aTrack, aStep);
112
}
113
114
// GetMeanFreePath
115
// ---------------
116
117
G4double
G4UCNMultiScattering::GetMeanFreePath
(
const
G4Track
& aTrack,
118
G4double
,
119
G4ForceCondition
*)
120
{
121
G4double
AttenuationLength =
DBL_MAX
;
122
123
const
G4Material
* aMaterial = aTrack.
GetMaterial
();
124
G4MaterialPropertiesTable
* aMaterialPropertiesTable =
125
aMaterial->
GetMaterialPropertiesTable
();
126
127
G4double
crossect = 0.0;
128
if
(aMaterialPropertiesTable) {
129
crossect = aMaterialPropertiesTable->
GetConstProperty
(
"SCATCS"
);
130
// if(crossect == 0.0)
131
// G4cout << "No UCN MultiScattering length specified" << G4endl;
132
}
133
// else G4cout << "No UCN MultiScattering length specified" << G4endl;
134
135
if
(crossect) {
136
137
// Calculate a UCN MultiScattering length for this cross section
138
139
G4double
density = aMaterial->
GetTotNbOfAtomsPerVolume
();
140
141
crossect *=
barn
;
142
143
// attenuation length in mm
144
AttenuationLength = 1./density/crossect;
145
}
146
147
return
AttenuationLength;
148
}
149
150
G4ThreeVector
G4UCNMultiScattering::Scatter
()
151
{
152
153
G4ThreeVector
final
(0.,0.,1.);
154
155
// Make a simple uniform distribution in 4 pi
156
// apply scattering, calculate angle phi, theta
157
158
G4double
theta
= std::acos(2*
G4UniformRand
()-1);
159
G4double
phi
=
G4UniformRand
() * 2 *
pi
;
160
161
final
.rotateY(theta);
162
final
.rotateZ(phi);
163
final
=
final
.unit();
164
165
return
final
;
166
}
geant4
tree
geant4-10.6-release
source
processes
hadronic
processes
src
G4UCNMultiScattering.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:51
using
1.8.2 with
ECCE GitHub integration