Geant4
10.03
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
G4PSPassageCellFluxForCylinder3D.cc
Go to the documentation of this file.
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
// $Id: G4PSPassageCellFluxForCylinder3D.cc 67992 2013-03-13 10:59:57Z gcosmo $
28
// GEANT4 tag $Name: geant4-09-04 $
29
//
30
// G4PSPassageCellFluxForCylinder3D
31
#include "
G4PSPassageCellFluxForCylinder3D.hh
"
32
#include "
G4PhysicalConstants.hh
"
33
35
// (Description)
36
// This is a primitive scorer class for scoring cell flux.
37
// The Cell Flux is defined by a track length divided by a geometry
38
// volume, where only tracks passing through the geometry are taken
39
// into account. e.g. the unit of Cell Flux is mm/mm3.
40
//
41
// If you want to score all tracks in the geometry volume,
42
// please use G4PSCellFlux.
43
//
44
// Created: 2007-08-14 Tsukasa ASO
45
// 2011-03-24 Give Size and Segmentation for relicated volume in cylinder.
47
48
G4PSPassageCellFluxForCylinder3D::G4PSPassageCellFluxForCylinder3D
(
G4String
name
,
49
G4int
ni,
G4int
nj,
G4int
nk,
50
G4int
di,
G4int
dj,
G4int
dk)
51
:
G4PSPassageCellFlux3D
(name,ni,nj,nk,di,dj,dk)
52
{
53
nSegment
[0] =
nSegment
[1] =
nSegment
[2] = 0;
54
}
55
56
G4PSPassageCellFluxForCylinder3D::G4PSPassageCellFluxForCylinder3D
(
G4String
name
,
const
G4String
& unit,
57
G4int
ni,
G4int
nj,
G4int
nk,
58
G4int
di,
G4int
dj,
G4int
dk)
59
:
G4PSPassageCellFlux3D
(name,unit,ni,nj,nk,di,dj,dk)
60
{
61
nSegment
[0] =
nSegment
[1] =
nSegment
[2] = 0;
62
}
63
64
G4PSPassageCellFluxForCylinder3D::~G4PSPassageCellFluxForCylinder3D
()
65
{;}
66
67
void
G4PSPassageCellFluxForCylinder3D::SetCylinderSize
(
G4double
dr,
G4double
dz){
68
cylinderSize
.setX(dz);
//Z Phi R
69
cylinderSize
.setY(
twopi
);
//Z Phi R
70
cylinderSize
.setZ(dr);
//Z Phi R
71
}
72
void
G4PSPassageCellFluxForCylinder3D::SetNumberOfSegments
(
G4int
nSeg[3]){
73
nSegment
[0] = nSeg[0];
// Z
74
nSegment
[1] = nSeg[1];
// Phi
75
nSegment
[2] = nSeg[2];
// R
76
}
77
G4double
G4PSPassageCellFluxForCylinder3D::ComputeVolume
(
G4Step
*,
G4int
idx){
78
G4double
r0 = (
cylinderSize
[2]/
nSegment
[2])*(idx);
79
G4double
r1 = (
cylinderSize
[2]/
nSegment
[2])*(idx+1);
80
G4double
dRArea = (r1*r1-r0*r0)*
pi
;
81
82
// cylinderSize is given in Half Size
83
G4double
fullz =
cylinderSize
[0]/
nSegment
[0]*2.;
84
G4double
phiRatio = 1./
nSegment
[1];
85
G4double
v = dRArea*fullz*phiRatio;
86
87
return
v;
88
}
89
G4PSPassageCellFluxForCylinder3D::SetCylinderSize
void SetCylinderSize(G4double dr, G4double dz)
Definition:
G4PSPassageCellFluxForCylinder3D.cc:67
G4InuclParticleNames::name
const char * name(G4int ptype)
Definition:
G4InuclParticleNames.hh:77
G4PSPassageCellFluxForCylinder3D::ComputeVolume
virtual G4double ComputeVolume(G4Step *, G4int idx)
Definition:
G4PSPassageCellFluxForCylinder3D.cc:77
G4int
int G4int
Definition:
G4Types.hh:78
twopi
static constexpr double twopi
Definition:
G4SIunits.hh:76
G4PSPassageCellFluxForCylinder3D::cylinderSize
G4ThreeVector cylinderSize
Definition:
G4PSPassageCellFluxForCylinder3D.hh:70
G4PSPassageCellFluxForCylinder3D::~G4PSPassageCellFluxForCylinder3D
virtual ~G4PSPassageCellFluxForCylinder3D()
Definition:
G4PSPassageCellFluxForCylinder3D.cc:64
G4PSPassageCellFlux3D
Definition:
G4PSPassageCellFlux3D.hh:49
G4Step
Definition:
G4Step.hh:76
G4PhysicalConstants.hh
G4PSPassageCellFluxForCylinder3D::G4PSPassageCellFluxForCylinder3D
G4PSPassageCellFluxForCylinder3D(G4String name, G4int ni=1, G4int nj=1, G4int nk=1, G4int depi=2, G4int depj=1, G4int depk=0)
Definition:
G4PSPassageCellFluxForCylinder3D.cc:48
G4PSPassageCellFluxForCylinder3D::nSegment
G4int nSegment[3]
Definition:
G4PSPassageCellFluxForCylinder3D.hh:71
G4PSPassageCellFluxForCylinder3D.hh
pi
static constexpr double pi
Definition:
G4SIunits.hh:75
G4double
double G4double
Definition:
G4Types.hh:76
G4PSPassageCellFluxForCylinder3D::SetNumberOfSegments
void SetNumberOfSegments(G4int nSeg[3])
Definition:
G4PSPassageCellFluxForCylinder3D.cc:72
G4String
Definition:
G4String.hh:45
geant4.10.03
source
digits_hits
scorer
src
G4PSPassageCellFluxForCylinder3D.cc
Generated on Thu Feb 14 2002 02:28:03 for Geant4 by
1.8.8