GENIEGenerator
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
RSHelicityAmplModelEMn.cxx
Go to the documentation of this file.
1 //____________________________________________________________________________
2 /*
3  Copyright (c) 2003-2024, 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 
18 using namespace genie;
19 using namespace genie::constants;
20 
21 //____________________________________________________________________________
22 RSHelicityAmplModelEMn::RSHelicityAmplModelEMn() :
23 RSHelicityAmplModelI("genie::RSHelicityAmplModelEMn")
24 {
25 
26 }
27 //____________________________________________________________________________
28 RSHelicityAmplModelEMn::RSHelicityAmplModelEMn(string config) :
29 RSHelicityAmplModelI("genie::RSHelicityAmplModelEMn", config)
30 {
31 
32 }
33 //____________________________________________________________________________
34 RSHelicityAmplModelEMn::~RSHelicityAmplModelEMn()
35 {
36 
37 }
38 //____________________________________________________________________________
39 const 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 //____________________________________________________________________________
genie::FKR
Simple struct-like class holding the Feynmann-Kislinger-Ravndall (FKR) baryon excitation model parame...
Definition: FKR.h:31
genie::FKR::Lamda
double Lamda
Definition: FKR.h:37
genie::FKR::R
double R
Definition: FKR.h:45
genie::Resonance_t
enum genie::EResonance Resonance_t
genie::RSHelicityAmpl
A class holding the Rein-Sehgal&#39;s helicity amplitudes.
Definition: RSHelicityAmpl.h:40
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96
genie::FKR::T
double T
Definition: FKR.h:46
genie::RSHelicityAmplModelEMn::Compute
const RSHelicityAmpl & Compute(Resonance_t res, const FKR &fkr) const
Definition: RSHelicityAmplModelEMn.cxx:40
genie::RSHelicityAmplModelI
Pure abstract base class. Defines the RSHelicityAmplModelI interface.
Definition: RSHelicityAmplModelI.h:27
pWARN
#define pWARN
Definition: Messenger.h:60
genie::FKR::S
double S
Definition: FKR.h:40
Generated on Wed Jul 17 2024 14:46:45 for GENIEGenerator by   doxygen 1.8.5