GENIEGenerator
Loading...
Searching...
No Matches
NNBarOscUtils.cxx
Go to the documentation of this file.
1//____________________________________________________________________________
2/*
3 Copyright (c) 2003-2025, The GENIE Collaboration
4 For the full text of the license visit http://copyright.genie-mc.org
5
6 Jeremy Hewes, Georgia Karagiorgi
7 University of Manchester
8*/
9//____________________________________________________________________________
10
11#include "Framework/Messenger/Messenger.h"
12#include "Framework/GHEP/GHepParticle.h"
13#include "Framework/ParticleData/PDGCodes.h"
14#include "Framework/ParticleData/PDGUtils.h"
15#include "Framework/ParticleData/PDGLibrary.h"
16#include "Framework/Utils/PrintUtils.h"
17#include "Physics/NuclearState/NuclearUtils.h"
18#include "Physics/NNBarOscillation/NNBarOscUtils.h"
19
20using namespace genie;
21
22//____________________________________________________________________________
23string genie::utils::nnbar_osc::AsString(NNBarOscMode_t ndm)
24{
25 // this just maps the decay mode integers to string descriptors. replaced. -j
26 switch(ndm) {
27 case (kNORandom) : return "Random mode";
28 break;
29 case (kNOpto1pip1pi0) : return "p + nbar --> pi+ pi0";
30 break;
31 case (kNOpto1pip2pi0) : return "p + nbar --> pi+ 2pi0";
32 break;
33 case (kNOpto1pip3pi0) : return "p + nbar --> pi+ 3pi0";
34 break;
35 case (kNOpto2pip1pim1pi0) : return "p + nbar --> 2pi+ pi- pi0";
36 break;
37 case (kNOpto2pip1pim2pi0) : return "p + nbar --> 2pi+ pi- 2pi0";
38 break;
39 case (kNOpto2pip1pim2o) : return "p + nbar --> 2pi+ pi- 2omega";
40 break;
41 case (kNOpto3pip2pim1pi0) : return "p + nbar --> 3pi+ 2pi- pi0";
42 break;
43 case (kNOnto1pip1pim) : return "n + nbar --> pi+ pi-";
44 break;
45 case (kNOnto2pi0) : return "n + nbar --> 2pi0";
46 break;
47 case (kNOnto1pip1pim1pi0) : return "n + nbar --> pi+ pi- pi0";
48 break;
49 case (kNOnto1pip1pim2pi0) : return "n + nbar --> pi+ pi- 2pi0";
50 break;
51 case (kNOnto1pip1pim3pi0) : return "n + nbar --> pi+ pi- 3pi0";
52 break;
53 case (kNOnto2pip2pim) : return "n + nbar --> 2pi+ 2pi-";
54 break;
55 case (kNOnto2pip2pim1pi0) : return "n + nbar --> 2pi+ 2pi- pi0";
56 break;
57 case (kNOnto1pip1pim1o) : return "n + nbar --> pi+ pi- omega";
58 break;
59 case (kNOnto2pip2pim2pi0) : return "n + nbar --> 2pi+ 2pi- 2pi0";
60 break;
61 default : return "?";
62 break;
63 }
64 return "??";
65}
66//____________________________________________________________________________
67bool genie::utils::nnbar_osc::IsValidMode(NNBarOscMode_t ndm)
68{
69 // checks if a mode is valid. just straight replaced. -j
70 switch(ndm) {
71 case (kNORandom) :
72 case (kNOpto1pip1pi0) :
73 case (kNOpto1pip2pi0) :
74 case (kNOpto1pip3pi0) :
75 case (kNOpto2pip1pim1pi0) :
76 case (kNOpto2pip1pim2pi0) :
77 case (kNOpto2pip1pim2o) :
78 case (kNOpto3pip2pim1pi0) :
79 case (kNOnto1pip1pim) :
80 case (kNOnto2pi0) :
81 case (kNOnto1pip1pim1pi0) :
82 case (kNOnto1pip1pim2pi0) :
83 case (kNOnto1pip1pim3pi0) :
84 case (kNOnto2pip2pim) :
85 case (kNOnto2pip2pim1pi0) :
86 case (kNOnto1pip1pim1o) :
87 case (kNOnto2pip2pim2pi0) :
88 return true;
89 break;
90 default :
91 return false;
92 break;
93 }
94 return false;
95}
96//____________________________________________________________________________
97int genie::utils::nnbar_osc::AnnihilatingNucleonPdgCode(NNBarOscMode_t ndm)
98{
99 // name isn't really accurate any more. instead of decayed nucleon, function
100 // returns what particle the oscillated neutron annihilated with -j
101 switch(ndm) {
102 case (kNOpto1pip1pi0) : return kPdgProton; break;
103 case (kNOpto1pip2pi0) : return kPdgProton; break;
104 case (kNOpto1pip3pi0) : return kPdgProton; break;
105 case (kNOpto2pip1pim1pi0) : return kPdgProton; break;
106 case (kNOpto2pip1pim2pi0) : return kPdgProton; break;
107 case (kNOpto2pip1pim2o) : return kPdgProton; break;
108 case (kNOpto3pip2pim1pi0) : return kPdgProton; break;
109 case (kNOnto1pip1pim) : return kPdgNeutron; break;
110 case (kNOnto2pi0) : return kPdgNeutron; break;
111 case (kNOnto1pip1pim1pi0) : return kPdgNeutron; break;
112 case (kNOnto1pip1pim2pi0) : return kPdgNeutron; break;
113 case (kNOnto1pip1pim3pi0) : return kPdgNeutron; break;
114 case (kNOnto2pip2pim) : return kPdgNeutron; break;
115 case (kNOnto2pip2pim1pi0) : return kPdgNeutron; break;
116 case (kNOnto1pip1pim1o) : return kPdgNeutron; break;
117 case (kNOnto2pip2pim2pi0) : return kPdgNeutron; break;
118 default : return 0; break;
119 }
120 return 0;
121}
122//____________________________________________________________________________
123PDGCodeList genie::utils::nnbar_osc::DecayProductList(
124 NNBarOscMode_t ndm)
125{
126 // ok so i think this is the first function where a straight replacement
127 // isn't gonna cut it. all the nucleon decay modes are two-body, but that is
128 // painfully untrue for nnbar. i just threw aaaaaaallll of the final state
129 // particles into the vector, so let's just hope for the best -j
130
131 // need to implement a lorentz boost into rest frame of two nucleons -j
132
133 bool allow_duplicate = true;
134 PDGCodeList decay_products(allow_duplicate);
135
136 switch(ndm) {
137 case (kNOpto1pip1pi0) :
138 decay_products.push_back(kPdgPiP);
139 decay_products.push_back(kPdgPi0);
140 break;
141 case (kNOpto1pip2pi0) :
142 decay_products.push_back(kPdgPiP);
143 decay_products.push_back(kPdgPi0);
144 decay_products.push_back(kPdgPi0);
145 break;
146 case (kNOpto1pip3pi0) :
147 decay_products.push_back(kPdgPiP);
148 decay_products.push_back(kPdgPi0);
149 decay_products.push_back(kPdgPi0);
150 decay_products.push_back(kPdgPi0);
151 break;
152 case (kNOpto2pip1pim1pi0) :
153 decay_products.push_back(kPdgPiP);
154 decay_products.push_back(kPdgPiP);
155 decay_products.push_back(kPdgPiM);
156 decay_products.push_back(kPdgPi0);
157 break;
158 case (kNOpto2pip1pim2pi0) :
159 decay_products.push_back(kPdgPiP);
160 decay_products.push_back(kPdgPiP);
161 decay_products.push_back(kPdgPiM);
162 decay_products.push_back(kPdgPi0);
163 decay_products.push_back(kPdgPi0);
164 break;
165 case (kNOpto2pip1pim2o) :
166 decay_products.push_back(kPdgPiP);
167 decay_products.push_back(kPdgPiP);
168 decay_products.push_back(kPdgPiM);
169 decay_products.push_back(kPdgomega);
170 decay_products.push_back(kPdgomega);
171 break;
172 case (kNOpto3pip2pim1pi0) :
173 decay_products.push_back(kPdgPiP);
174 decay_products.push_back(kPdgPiP);
175 decay_products.push_back(kPdgPiP);
176 decay_products.push_back(kPdgPiM);
177 decay_products.push_back(kPdgPiM);
178 decay_products.push_back(kPdgPi0);
179 break;
180 case (kNOnto1pip1pim) :
181 decay_products.push_back(kPdgPiP);
182 decay_products.push_back(kPdgPiM);
183 break;
184 case (kNOnto2pi0) :
185 decay_products.push_back(kPdgPi0);
186 decay_products.push_back(kPdgPi0);
187 break;
188 case (kNOnto1pip1pim1pi0) :
189 decay_products.push_back(kPdgPiP);
190 decay_products.push_back(kPdgPiM);
191 decay_products.push_back(kPdgPi0);
192 break;
193 case (kNOnto1pip1pim2pi0) :
194 decay_products.push_back(kPdgPiP);
195 decay_products.push_back(kPdgPiM);
196 decay_products.push_back(kPdgPi0);
197 decay_products.push_back(kPdgPi0);
198 break;
199 case (kNOnto1pip1pim3pi0) :
200 decay_products.push_back(kPdgPiP);
201 decay_products.push_back(kPdgPiM);
202 decay_products.push_back(kPdgPi0);
203 decay_products.push_back(kPdgPi0);
204 decay_products.push_back(kPdgPi0);
205 break;
206 case (kNOnto2pip2pim) :
207 decay_products.push_back(kPdgPiP);
208 decay_products.push_back(kPdgPiP);
209 decay_products.push_back(kPdgPiM);
210 decay_products.push_back(kPdgPiM);
211 break;
212 case (kNOnto2pip2pim1pi0) :
213 decay_products.push_back(kPdgPiP);
214 decay_products.push_back(kPdgPiP);
215 decay_products.push_back(kPdgPiM);
216 decay_products.push_back(kPdgPiM);
217 decay_products.push_back(kPdgPi0);
218 break;
219 case (kNOnto1pip1pim1o) :
220 decay_products.push_back(kPdgPiP);
221 decay_products.push_back(kPdgPiM);
222 decay_products.push_back(kPdgomega);
223 break;
224 case (kNOnto2pip2pim2pi0) :
225 decay_products.push_back(kPdgPiP);
226 decay_products.push_back(kPdgPiP);
227 decay_products.push_back(kPdgPiM);
228 decay_products.push_back(kPdgPiM);
229 decay_products.push_back(kPdgPi0);
230 decay_products.push_back(kPdgPi0);
231 break;
232 default :
233 break;
234 }
235 return decay_products;
236}
237//____________________________________________________________________________
238GHepStatus_t genie::utils::nnbar_osc::DecayProductStatus(
239 bool in_nucleus, int pdgc)
240{
241 // took out all the irrelevant particles -j
242 if(in_nucleus) {
243 if( pdgc == kPdgPi0 ||
244 pdgc == kPdgPiM ||
245 pdgc == kPdgPiP ||
246 pdgc == kPdgomega)
247 {
248 return kIStHadronInTheNucleus;
249 }
250 }
251
252 return kIStStableFinalState;
253}
254//____________________________________________________________________________
PDGCodes.h
Most commonly used PDG codes. A set of utility functions to handle PDG codes is provided in PDGUtils.
genie::PDGCodeList
A list of PDG codes.
Definition PDGCodeList.h:32
genie::PDGCodeList::push_back
void push_back(int pdg_code)
Definition PDGCodeList.cxx:58
genie::utils::nnbar_osc::DecayProductStatus
GHepStatus_t DecayProductStatus(bool in_nucleus, int pdgc)
Definition NNBarOscUtils.cxx:238
genie::utils::nnbar_osc::IsValidMode
bool IsValidMode(NNBarOscMode_t ndm)
Definition NNBarOscUtils.cxx:67
genie::utils::nnbar_osc::AnnihilatingNucleonPdgCode
int AnnihilatingNucleonPdgCode(NNBarOscMode_t ndm)
Definition NNBarOscUtils.cxx:97
genie::utils::nnbar_osc::DecayProductList
PDGCodeList DecayProductList(NNBarOscMode_t ndm)
Definition NNBarOscUtils.cxx:123
genie::utils::nnbar_osc::AsString
string AsString(NNBarOscMode_t ndm)
Definition NNBarOscUtils.cxx:23
genie
THE MAIN GENIE PROJECT NAMESPACE
Definition AlgCmp.h:25
genie::kPdgPiM
const int kPdgPiM
Definition PDGCodes.h:159
genie::kIStHadronInTheNucleus
@ kIStHadronInTheNucleus
Definition GHepStatus.h:37
genie::kIStStableFinalState
@ kIStStableFinalState
Definition GHepStatus.h:30
genie::kNOnto2pip2pim2pi0
@ kNOnto2pip2pim2pi0
Definition NNBarOscMode.h:48
genie::kNOnto1pip1pim3pi0
@ kNOnto1pip1pim3pi0
Definition NNBarOscMode.h:44
genie::kNOpto1pip1pi0
@ kNOpto1pip1pi0
Definition NNBarOscMode.h:33
genie::kNOpto1pip3pi0
@ kNOpto1pip3pi0
Definition NNBarOscMode.h:35
genie::kNOnto2pip2pim
@ kNOnto2pip2pim
Definition NNBarOscMode.h:45
genie::kNOnto2pi0
@ kNOnto2pi0
Definition NNBarOscMode.h:41
genie::kNORandom
@ kNORandom
Definition NNBarOscMode.h:32
genie::kNOnto1pip1pim1pi0
@ kNOnto1pip1pim1pi0
Definition NNBarOscMode.h:42
genie::kNOpto3pip2pim1pi0
@ kNOpto3pip2pim1pi0
Definition NNBarOscMode.h:39
genie::kNOpto2pip1pim2o
@ kNOpto2pip1pim2o
Definition NNBarOscMode.h:38
genie::kNOnto1pip1pim
@ kNOnto1pip1pim
Definition NNBarOscMode.h:40
genie::kNOnto1pip1pim1o
@ kNOnto1pip1pim1o
Definition NNBarOscMode.h:47
genie::kNOpto2pip1pim2pi0
@ kNOpto2pip1pim2pi0
Definition NNBarOscMode.h:37
genie::kNOnto1pip1pim2pi0
@ kNOnto1pip1pim2pi0
Definition NNBarOscMode.h:43
genie::kNOpto2pip1pim1pi0
@ kNOpto2pip1pim1pi0
Definition NNBarOscMode.h:36
genie::kNOpto1pip2pi0
@ kNOpto1pip2pi0
Definition NNBarOscMode.h:34
genie::kNOnto2pip2pim1pi0
@ kNOnto2pip2pim1pi0
Definition NNBarOscMode.h:46
genie::NNBarOscMode_t
enum genie::ENNBarOscMode NNBarOscMode_t
genie::kPdgProton
const int kPdgProton
Definition PDGCodes.h:81
genie::kPdgPi0
const int kPdgPi0
Definition PDGCodes.h:160
genie::kPdgNeutron
const int kPdgNeutron
Definition PDGCodes.h:83
genie::GHepStatus_t
enum genie::EGHepStatus GHepStatus_t
genie::kPdgomega
const int kPdgomega
Definition PDGCodes.h:168
genie::kPdgPiP
const int kPdgPiP
Definition PDGCodes.h:158
Generated on for GENIEGenerator by doxygen 1.14.0