#include <GiBUURESFormFactor.h>
|
| FormFactors () |
|
| ~FormFactors () |
|
double | C3V (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | C4V (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | C5V (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | C6V (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | C3A (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | C4A (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | C5A (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | C6A (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | F1V (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | F2V (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | FA (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | FP (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it) const |
|
double | Q2min (void) const |
|
double | Q2max (void) const |
|
|
static double | fMinQ2 = 0.0 |
| min Q2 for which resonance f/f data are given More...
|
|
static double | fMaxQ2 = 4.0 |
| max Q2 for which resonance f/f data are given More...
|
|
static const int | kNRes = 18 |
|
static const int | kNCurr = 3 |
|
static const int | kNHitNuc = 2 |
|
static const int | kNFFRes = 12 |
|
Definition at line 47 of file GiBUURESFormFactor.h.
GiBUURESFormFactor::FormFactors::FormFactors |
( |
void |
| ) |
|
GiBUURESFormFactor::FormFactors::~FormFactors |
( |
void |
| ) |
|
Definition at line 83 of file GiBUURESFormFactor.cxx.
86 cout <<
"GiBUURESFormFactor singleton dtor: Deleting all f/f splines" <<
endl;
91 for(
int i=0; i<
kNCurr; i++) {
QTextStream & endl(QTextStream &s)
Definition at line 257 of file GiBUURESFormFactor.cxx.
261 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,8);
double Q2(const Interaction *const i)
bool IsN(Resonance_t res)
is it an N resonance?
Definition at line 229 of file GiBUURESFormFactor.cxx.
233 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,4);
double Q2(const Interaction *const i)
bool IsN(Resonance_t res)
is it an N resonance?
Definition at line 264 of file GiBUURESFormFactor.cxx.
268 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,9);
double Q2(const Interaction *const i)
bool IsN(Resonance_t res)
is it an N resonance?
Definition at line 236 of file GiBUURESFormFactor.cxx.
240 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,5);
double Q2(const Interaction *const i)
bool IsN(Resonance_t res)
is it an N resonance?
Definition at line 271 of file GiBUURESFormFactor.cxx.
275 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,10);
double Q2(const Interaction *const i)
bool IsN(Resonance_t res)
is it an N resonance?
Definition at line 243 of file GiBUURESFormFactor.cxx.
247 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,6);
double Q2(const Interaction *const i)
bool IsN(Resonance_t res)
is it an N resonance?
Definition at line 278 of file GiBUURESFormFactor.cxx.
282 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,11);
double Q2(const Interaction *const i)
bool IsN(Resonance_t res)
is it an N resonance?
Definition at line 250 of file GiBUURESFormFactor.cxx.
254 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,7);
double Q2(const Interaction *const i)
bool IsN(Resonance_t res)
is it an N resonance?
Definition at line 285 of file GiBUURESFormFactor.cxx.
289 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,0);
bool IsDelta(Resonance_t res)
is it a Delta resonance?
double Q2(const Interaction *const i)
Definition at line 292 of file GiBUURESFormFactor.cxx.
296 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,1);
bool IsDelta(Resonance_t res)
is it a Delta resonance?
double Q2(const Interaction *const i)
Definition at line 299 of file GiBUURESFormFactor.cxx.
303 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,2);
bool IsDelta(Resonance_t res)
is it a Delta resonance?
double Q2(const Interaction *const i)
func to retrieve interpolated form factor values
Definition at line 313 of file GiBUURESFormFactor.cxx.
317 if(Q2 < fMinQ2 || Q2 >
fMaxQ2)
return 0.;
319 int r = -1, i = -1, j = -1;
321 if(ffresid<0 || ffresid >=
kNFFRes)
return 0.;
324 if(r<0 || r >=
kNRes)
return 0.;
328 else if (it ==
kIntEM) { i = 2; }
331 else if (hit_nucleon_pdg ==
kPdgProton ) { j = 1; }
double Q2(const Interaction *const i)
A numeric analysis tool class for interpolating 1-D functions.
double Evaluate(double x) const
Definition at line 306 of file GiBUURESFormFactor.cxx.
310 return this->
FFRes(
Q2,res,hit_nucleon_pdg,it,3);
bool IsDelta(Resonance_t res)
is it a Delta resonance?
double Q2(const Interaction *const i)
void GiBUURESFormFactor::FormFactors::LoadTables |
( |
void |
| ) |
|
|
private |
load all form factor data tables
Definition at line 104 of file GiBUURESFormFactor.cxx.
109 for(
int i=0; i<
kNCurr; i++) {
119 string(
"/data/evgen/gibuu");
127 for(
int i=0; i<
kNCurr; i++) {
134 ostringstream datafile;
135 datafile << data_dir <<
"/form_factors/";
138 case ( 0): datafile <<
"P33_1232";
break;
139 case ( 1): datafile <<
"S11_1535";
break;
140 case ( 2): datafile <<
"D13_1520";
break;
141 case ( 3): datafile <<
"S11_1650";
break;
142 case ( 5): datafile <<
"D15_1675";
break;
143 case ( 6): datafile <<
"S31_1620";
break;
144 case ( 7): datafile <<
"D33_1700";
break;
145 case ( 8): datafile <<
"P11_1440";
break;
146 case (10): datafile <<
"P13_1720";
break;
147 case (11): datafile <<
"F15_1680";
break;
148 case (12): datafile <<
"P31_1910";
break;
149 case (14): datafile <<
"F35_1905";
break;
150 case (15): datafile <<
"F37_1950";
break;
151 default : skip =
true;
155 case (0): datafile <<
"_CC";
break;
156 case (1): datafile <<
"_NC";
break;
157 case (2): datafile <<
"_EM";
break;
158 default : skip =
true;
161 case (0): datafile <<
"_neutron";
break;
162 case (1): datafile <<
"_proton";
break;
163 default : skip =
true;
165 datafile <<
"_FFres.dat";
170 assert( ! gSystem->AccessPathName(datafile.str().c_str()) );
180 data_ffres.ReadFile(datafile.str().c_str(),
181 "Q2/D:f1/D:f2/D:f3/D:f4/D:f5/D:f6/D:f7/D:f8/D");
184 <<
"Number of data rows: " << data_ffres.GetEntries();
216 for(
int i=0; i<
kNCurr; i++) {
226 <<
"Done loading all resonance form factor files...";
bool IsDelta(Resonance_t res)
is it a Delta resonance?
A numeric analysis tool class for interpolating 1-D functions.
enum genie::EResonance Resonance_t
bool IsN(Resonance_t res)
is it an N resonance?
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
void YCanBeNegative(bool tf)
double genie::GiBUURESFormFactor::FormFactors::Q2max |
( |
void |
| ) |
const |
|
inline |
double genie::GiBUURESFormFactor::FormFactors::Q2min |
( |
void |
| ) |
const |
|
inline |
double GiBUURESFormFactor::FormFactors::fMaxQ2 = 4.0 |
|
staticprivate |
double GiBUURESFormFactor::FormFactors::fMinQ2 = 0.0 |
|
staticprivate |
const int genie::GiBUURESFormFactor::FormFactors::kNCurr = 3 |
|
staticprivate |
const int genie::GiBUURESFormFactor::FormFactors::kNFFRes = 12 |
|
staticprivate |
const int genie::GiBUURESFormFactor::FormFactors::kNHitNuc = 2 |
|
staticprivate |
const int genie::GiBUURESFormFactor::FormFactors::kNRes = 18 |
|
staticprivate |
The documentation for this class was generated from the following files: