Geant4 11.1.1
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4CascadeMuMinusPChannel.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
28#include "G4CascadeMuMinusPChannel.hh"
29#include "G4InuclParticleNames.hh"
30using namespace G4InuclParticleNames;
31
32namespace {
33 // mu- p : Outgoing particle types of a given multiplicity
34 static const G4int mump2bfs[1][2] =
35 {{neu,mnu}};
36
37 static const G4int mump3bfs[1][3] =
38 {{neu,mnu,gam}};
39
40 static const G4int mump4bfs[1][4] =
41 {{neu,mnu,gam,gam}};
42
43 static const G4int mump5bfs[1][5] =
44 {{neu,mnu,gam,gam,gam}};
45
46 static const G4int mump6bfs[1][6] =
47 {{neu,mnu,gam,gam,gam,gam}};
48
49 static const G4int mump7bfs[1][7] =
50 {{neu,mnu,gam,gam,gam,gam,gam}};
51
52 static const G4int mump8bfs[1][8] =
53 {{neu,mnu,gam,gam,gam,gam,gam,gam}};
54
55 static const G4int mump9bfs[1][9] =
56 {{neu,mnu,gam,gam,gam,gam,gam,gam,gam}};
57}
58
59namespace {
60 // Total mu- p cross section as a function of kinetic energy
61 static const G4double mumpTotXSec[30] =
62 { 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01,
63 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01,
64 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01};
65/*
66 {100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0,
67 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0,
68 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0};
69*/
70 static const G4double mumpCrossSections[8][30] = {
71 //
72 // multiplicity 2 (1 channel)
73 //
74 // nu_mu n
75 { 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01,
76 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01,
77 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
78/*
79 {100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0,
80 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0,
81 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0},
82*/
83 //
84 // multiplicity 3 (1 channel)
85 //
86 // nu_mu n gam
87 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
88 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
89 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
90 //
91 // multiplicity 4 (1 channel)
92 //
93 // nu_mu n gam gam
94 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
95 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
96 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
97 //
98 // multiplicity 5 (1 channel)
99 //
100 // nu_mu n gam gam gam
101 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
102 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
103 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
104 //
105 // multiplicity 6 (1 channel)
106 //
107 // nu_mu n gam gam gam gam
108 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
109 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
110 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
111 //
112 // multiplicity 7 (1 channel)
113 //
114 // nu_mu n gam gam gam gam gam
115 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
116 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
117 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
118 //
119 // multiplicity 8 (1 channel)
120 //
121 // nu_mu n gam gam gam gam gam gam
122 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
123 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
124 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
125 //
126 // multiplicity 9 (1 channel)
127 //
128 // nu_mu n gam gam gam gam gam gam gam
129 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
130 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
131 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0}};
132}
133
134
135const G4CascadeMuMinusPChannelData::data_t
136G4CascadeMuMinusPChannelData::data(mump2bfs, mump3bfs, mump4bfs, mump5bfs,
137 mump6bfs, mump7bfs, mump8bfs, mump9bfs,
138 mumpCrossSections, mumpTotXSec, mum*pro,
139 "MuMinusProton");
140
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85