GENIEGenerator
Loading...
Searching...
No Matches
RSHelicityAmplModelEMn.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 Costas Andreopoulos <c.andreopoulos \at cern.ch>
7 University of Liverpool
8*/
9//____________________________________________________________________________
10
11#include "Framework/Algorithm/AlgConfigPool.h"
12#include "Framework/ParticleData/BaryonResUtils.h"
13#include "Framework/Conventions/Constants.h"
14#include "Framework/Messenger/Messenger.h"
15#include "Physics/Resonance/XSection/RSHelicityAmplModelEMn.h"
16#include "Physics/Resonance/XSection/RSHelicityAmpl.h"
17
18using namespace genie;
19using namespace genie::constants;
20
21//____________________________________________________________________________
22RSHelicityAmplModelEMn::RSHelicityAmplModelEMn() :
23RSHelicityAmplModelI("genie::RSHelicityAmplModelEMn")
24{
25
26}
27//____________________________________________________________________________
28RSHelicityAmplModelEMn::RSHelicityAmplModelEMn(string config) :
29RSHelicityAmplModelI("genie::RSHelicityAmplModelEMn", config)
30{
31
32}
33//____________________________________________________________________________
38//____________________________________________________________________________
39const RSHelicityAmpl &
40 RSHelicityAmplModelEMn::Compute(
41 Resonance_t res, const FKR & fkr) const
42{
43 switch(res) {
44
45 case (kP33_1232) :
46 {
47 fAmpl.fPlus1 = kSqrt2 * fkr.R;
48 fAmpl.fPlus3 = kSqrt6 * fkr.R;
49 fAmpl.fMinus1 = -1 * fAmpl.fPlus1;
50 fAmpl.fMinus3 = -1 * fAmpl.fPlus3;
51 fAmpl.f0Minus = 0.;
52 fAmpl.f0Plus = 0.;
53 break;
54 }
55 case (kS11_1535) :
56 {
57 fAmpl.fPlus1 = kSqrt3 * fkr.T + k1_Sqrt6 * fkr.Lamda * fkr.R;
58 fAmpl.f0Minus = kSqrt3_2 * fkr.Lamda * fkr.S;
59 fAmpl.fMinus1 = -1 * fAmpl.fPlus1;
60 fAmpl.f0Plus = -1 * fAmpl.f0Minus;
61 fAmpl.fMinus3 = 0.;
62 fAmpl.fPlus3 = 0.;
63
64 break;
65 }
66 case (kD13_1520) :
67 {
68 fAmpl.fMinus1 = -kSqrt3_2 * fkr.T + k1_Sqrt3 * fkr.Lamda * fkr.R;
69 fAmpl.fMinus3 = -k3_Sqrt2 * fkr.T;
70 fAmpl.f0Minus = kSqrt3 * fkr.Lamda * fkr.S;
71 fAmpl.fPlus1 = fAmpl.fMinus1;
72 fAmpl.fPlus3 = fAmpl.fMinus3;
73 fAmpl.f0Plus = fAmpl.f0Minus;
74 break;
75 }
76 case (kS11_1650) :
77 {
78 fAmpl.fPlus1 = k1_Sqrt6 * fkr.Lamda * fkr.R;
79 fAmpl.fMinus1 = -1 * fAmpl.fPlus1;
80 fAmpl.fMinus3 = 0.;
81 fAmpl.fPlus3 = 0.;
82 fAmpl.f0Minus = 0.;
83 fAmpl.f0Plus = 0.;
84 break;
85 }
86 case (kD13_1700) :
87 {
88 double LR = fkr.Lamda * fkr.R;
89
90 fAmpl.fMinus1 = -(1./kSqrt30) * LR;
91 fAmpl.fMinus3 = -(3./kSqrt10) * LR;
92 fAmpl.fPlus1 = fAmpl.fMinus1;
93 fAmpl.fPlus3 = fAmpl.fMinus3;
94 fAmpl.f0Minus = 0.;
95 fAmpl.f0Plus = 0.;
96 break;
97 }
98 case (kD15_1675) :
99 {
100 double LR = fkr.Lamda * fkr.R;
101
102 fAmpl.fMinus1 = kSqrt3_10 * LR;
103 fAmpl.fMinus3 = kSqrt3_5 * LR;
104 fAmpl.fPlus1 = -1 * fAmpl.fMinus1;
105 fAmpl.fPlus3 = -1 * fAmpl.fMinus3;
106 fAmpl.f0Minus = 0.;
107 fAmpl.f0Plus = 0.;
108 break;
109 }
110 case (kS31_1620) :
111 {
112 fAmpl.fMinus1 = kSqrt3 * fkr.T - k1_Sqrt6 * fkr.Lamda * fkr.R;
113 fAmpl.f0Minus = -kSqrt3_2 * fkr.Lamda * fkr.S;
114 fAmpl.fPlus1 = -1. * fAmpl.fMinus1;
115 fAmpl.f0Plus = -1. * fAmpl.f0Minus;
116 fAmpl.fMinus3 = 0.;
117 fAmpl.fPlus3 = 0.;
118 break;
119 }
120 case (kD33_1700) :
121 {
122 fAmpl.fMinus1 = kSqrt3_2 * fkr.T + k1_Sqrt3 * fkr.Lamda * fkr.R;
123 fAmpl.fMinus3 = k3_Sqrt2 * fkr.T;
124 fAmpl.f0Minus = -kSqrt3 * fkr.Lamda * fkr.S;
125 fAmpl.fPlus1 = fAmpl.fMinus1;
126 fAmpl.fPlus3 = fAmpl.fMinus3;
127 fAmpl.f0Plus = fAmpl.f0Minus;
128 break;
129 }
130 case (kP11_1440) :
131 {
132 fAmpl.fMinus1 = k1_Sqrt3 * TMath::Power(fkr.Lamda, 2) * fkr.R;
133 fAmpl.fPlus1 = fAmpl.fMinus1;
134 fAmpl.fMinus3 = 0.;
135 fAmpl.fPlus3 = 0.;
136 fAmpl.f0Minus = 0.;
137 fAmpl.f0Plus = 0.;
138 break;
139 }
140 case (kP33_1600) :
141 {
142 double L2R = TMath::Power(fkr.Lamda, 2) * fkr.R;
143
144 fAmpl.fMinus1 = k1_Sqrt6 * L2R;
145 fAmpl.fMinus3 = k1_Sqrt2 * L2R;
146 fAmpl.fPlus1 = -1. * fAmpl.fMinus1;
147 fAmpl.fPlus3 = -1. * fAmpl.fMinus3;
148 fAmpl.f0Minus = 0.;
149 fAmpl.f0Plus = 0.;
150 break;
151 }
152 case (kP13_1720) :
153 {
154 fAmpl.fMinus1 = k2_Sqrt15 * TMath::Power(fkr.Lamda, 2) * fkr.R;
155 fAmpl.fPlus1 = -1 * fAmpl.fMinus1;
156 fAmpl.fMinus3 = 0.;
157 fAmpl.fPlus3 = 0.;
158 fAmpl.f0Minus = 0.;
159 fAmpl.f0Plus = 0.;
160 break;
161 }
162 case (kF15_1680) :
163 {
164 fAmpl.fMinus1 = -kSqrt2_5 * TMath::Power(fkr.Lamda, 2) * fkr.R;
165 fAmpl.fPlus1 = fAmpl.fMinus1;
166 fAmpl.fMinus3 = 0.;
167 fAmpl.fPlus3 = 0.;
168 fAmpl.f0Minus = 0.;
169 fAmpl.f0Plus = 0.;
170 break;
171 }
172 case (kP31_1910) :
173 {
174 fAmpl.fMinus1 = -k1_Sqrt15 * TMath::Power(fkr.Lamda, 2) * fkr.R;
175 fAmpl.fPlus1 = fAmpl.fMinus1;
176 fAmpl.fMinus3 = 0.;
177 fAmpl.fPlus3 = 0.;
178 fAmpl.f0Minus = 0.;
179 fAmpl.f0Plus = 0.;
180 break;
181 }
182 case (kP33_1920) :
183 {
184 double L2R = TMath::Power(fkr.Lamda, 2) * fkr.R;
185
186 fAmpl.fMinus1 = k1_Sqrt15 * L2R;
187 fAmpl.fMinus3 = -k1_Sqrt5 * L2R;
188 fAmpl.fPlus1 = -1.* fAmpl.fMinus1;
189 fAmpl.fPlus3 = -1.* fAmpl.fMinus3;
190 fAmpl.f0Minus = 0.;
191 fAmpl.f0Plus = 0.;
192 break;
193 }
194 case (kF35_1905) :
195 {
196 double L2R = TMath::Power(fkr.Lamda, 2) * fkr.R;
197
198 fAmpl.fMinus1 = k1_Sqrt35 * L2R;
199 fAmpl.fMinus3 = kSqrt18_35 * L2R;
200 fAmpl.fPlus1 = fAmpl.fMinus1;
201 fAmpl.fPlus3 = fAmpl.fMinus3;
202 fAmpl.f0Minus = 0.;
203 fAmpl.f0Plus = 0.;
204 break;
205 }
206 case (kF37_1950) :
207 {
208 double L2R = TMath::Power(fkr.Lamda, 2) * fkr.R;
209
210 fAmpl.fMinus1 = -kSqrt6_35 * L2R;
211 fAmpl.fMinus3 = -kSqrt2_7 * L2R;
212 fAmpl.fPlus1 = -1. * fAmpl.fMinus1;
213 fAmpl.fPlus3 = -1. * fAmpl.fMinus3;
214 fAmpl.f0Minus = 0.;
215 fAmpl.f0Plus = 0.;
216 break;
217 }
218 case (kP11_1710) :
219 {
220 double L2 = TMath::Power(fkr.Lamda, 2);
221
222 fAmpl.fMinus1 = -k1_Sqrt24 * L2 * fkr.R;
223 fAmpl.f0Minus = -kSqrt3_8 * L2 * fkr.S;
224 fAmpl.fPlus1 = fAmpl.fMinus1;
225 fAmpl.f0Plus = fAmpl.f0Minus;
226 fAmpl.fMinus3 = 0.;
227 fAmpl.fPlus3 = 0.;
228
229 break;
230 }
231 case (kF17_1970) :
232 {
233 double L2R = TMath::Power(fkr.Lamda, 2) * fkr.R;
234
235 fAmpl.fMinus1 = kSqrt3_35 * L2R;
236 fAmpl.fPlus1 = -1 * fAmpl.fMinus1;
237 fAmpl.fMinus3 = k1_Sqrt7 * L2R;
238 fAmpl.fPlus3 = -1 * fAmpl.fMinus3;
239 fAmpl.f0Minus = 0.;
240 fAmpl.f0Plus = 0.;
241 break;
242 }
243 default:
244 {
245 LOG("RSHAmpl", pWARN) << "*** UNRECOGNIZED RESONANCE!";
246 fAmpl.fMinus1 = 0.;
247 fAmpl.fPlus1 = 0.;
248 fAmpl.fMinus3 = 0.;
249 fAmpl.fPlus3 = 0.;
250 fAmpl.f0Minus = 0.;
251 fAmpl.f0Plus = 0.;
252 break;
253 }
254
255 }//switch
256
257 return fAmpl;
258}
259//____________________________________________________________________________
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE,...
Definition Messenger.h:96
pWARN
#define pWARN
Definition Messenger.h:60
genie::FKR
Simple struct-like class holding the Feynmann-Kislinger-Ravndall (FKR) baryon excitation model parame...
Definition FKR.h:31
genie::FKR::T
double T
Definition FKR.h:46
genie::FKR::S
double S
Definition FKR.h:40
genie::FKR::R
double R
Definition FKR.h:45
genie::FKR::Lamda
double Lamda
Definition FKR.h:37
genie::RSHelicityAmplModelEMn::Compute
const RSHelicityAmpl & Compute(Resonance_t res, const FKR &fkr) const
Definition RSHelicityAmplModelEMn.cxx:40
genie::RSHelicityAmpl
A class holding the Rein-Sehgal's helicity amplitudes.
Definition RSHelicityAmpl.h:40
genie::constants
Basic constants.
genie
THE MAIN GENIE PROJECT NAMESPACE
Definition AlgCmp.h:25
genie::Resonance_t
enum genie::EResonance Resonance_t
Generated on for GENIEGenerator by doxygen 1.14.0