DipoleAxialFormFactorModel.cxx
Go to the documentation of this file.
1 //____________________________________________________________________________
2 /*
3  Copyright (c) 2003-2019, The GENIE Collaboration
4  For the full text of the license visit http://copyright.genie-mc.org
5  or see $GENIE/LICENSE
6 
7 \author Aaron Meyer <asmeyer2012 \at uchicago.edu>
8  based off DipoleELFormFactorsModel by
9  Costas Andreopoulos <costas.andreopoulos \at stfc.ac.uk>
10  STFC, Rutherford Appleton Laboratory
11 
12  For the class documentation see the corresponding header file.
13 
14 */
15 //____________________________________________________________________________
16 
17 #include <TMath.h>
18 
23 
24 using namespace genie;
25 
26 //____________________________________________________________________________
28 AxialFormFactorModelI("genie::DipoleAxialFormFactorModel")
29 {
30 
31 }
32 //____________________________________________________________________________
34 AxialFormFactorModelI("genie::DipoleAxialFormFactorModel", config)
35 {
36 
37 }
38 //____________________________________________________________________________
40 {
41 
42 }
43 //____________________________________________________________________________
44 double DipoleAxialFormFactorModel::FA(const Interaction * interaction) const
45 {
46  // calculate and return FA
47  double q2 = interaction->Kine().q2();
48  double fa = fFA0 / TMath::Power(1-q2/fMa2, 2);
49 
50  LOG("AxialFormFactor", pDEBUG) << "FA(q^2 = " << q2 << ") = " << fa;
51  return fa;
52 }
53 //____________________________________________________________________________
55 {
56  Algorithm::Configure(config);
57  this->LoadConfig();
58 }
59 //____________________________________________________________________________
61 {
62  Algorithm::Configure(param_set);
63  this->LoadConfig();
64 }
65 //____________________________________________________________________________
67 {
68  GetParam( "QEL-Ma", fMa ) ;
69  fMa2 = TMath::Power(fMa,2);
70 
71  // FA(q2 = 0)
72  GetParam( "QEL-FA0", fFA0 ) ;
73 }
74 //____________________________________________________________________________
75 
THE MAIN GENIE PROJECT NAMESPACE
Definition: GeneratorBase.h:8
double FA(const Interaction *interaction) const
Compute the axial form factor.
Pure abstract base class. Defines the AxialFormFactorModelI interface to be implemented by LlewellynS...
Definition: config.py:1
Double_t q2[12][num]
Definition: f2_nu.C:137
Summary information for an interaction.
Definition: Interaction.h:56
double q2(bool selected=false) const
Definition: Kinematics.cxx:151
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:97
const Kinematics & Kine(void) const
Definition: Interaction.h:71
virtual void Configure(const Registry &config)
Definition: Algorithm.cxx:70
A registry. Provides the container for algorithm configuration parameters.
Definition: Registry.h:66
bool GetParam(const RgKey &name, T &p, bool is_top_call=true) const
#define pDEBUG
Definition: Messenger.h:64
TFile fa("Li7.root")