#include <GiBUURESFormFactor.h>

Collaboration diagram for genie::GiBUURESFormFactor::FormFactors:

Public Member Functions
	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

Private Member Functions
double	FFRes (double Q2, Resonance_t res, int nucleon_pdg, InteractionType_t it, int ffid) const
	func to retrieve interpolated form factor values More...

void	LoadTables (void)
	load all form factor data tables More...

Private Attributes
Spline *	fFFRes [kNRes][kNCurr][kNHitNuc][kNFFRes]
	actual form factor data = f(Q2) More...

Static Private Attributes
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

Friends
class	GiBUURESFormFactor

Detailed Description

Definition at line 47 of file GiBUURESFormFactor.h.

Constructor & Destructor Documentation

GiBUURESFormFactor::FormFactors::FormFactors ( void )

Definition at line 78 of file GiBUURESFormFactor.cxx.

 {
 
 }

GiBUURESFormFactor::FormFactors::~FormFactors ( void )

Definition at line 83 of file GiBUURESFormFactor.cxx.

References genie::gAbortingInErr.

 {
   if(!gAbortingInErr) {
     cout << "GiBUURESFormFactor singleton dtor: Deleting all f/f splines" << endl;
   }
 
   // resonance form factor splines
   for(int r=0; r<kNRes; r++) {
     for(int i=0; i<kNCurr; i++) {
       for(int j=0; j<kNHitNuc; j++) {
         for(int k=0; k<kNFFRes; k++) {
            if (fFFRes[r][i][j][k]) {
              delete fFFRes[r][i][j][k];
              fFFRes[r][i][j][k] = 0;
            }
         }//k
       }//j
     }//i
   }//r
 }

Member Function Documentation

double GiBUURESFormFactor::FormFactors::C3A	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 257 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsN().

Referenced by SaveToRootFile().

 {
   if(!res::IsN(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,8);
 }

double GiBUURESFormFactor::FormFactors::C3V	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 229 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsN().

Referenced by SaveToRootFile().

 {
   if(!res::IsN(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,4);
 }

double GiBUURESFormFactor::FormFactors::C4A	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 264 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsN().

Referenced by SaveToRootFile().

 {
   if(!res::IsN(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,9);
 }

double GiBUURESFormFactor::FormFactors::C4V	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 236 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsN().

Referenced by SaveToRootFile().

 {
   if(!res::IsN(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,5);
 }

double GiBUURESFormFactor::FormFactors::C5A	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 271 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsN().

Referenced by SaveToRootFile().

 {
   if(!res::IsN(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,10);
 }

double GiBUURESFormFactor::FormFactors::C5V	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 243 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsN().

Referenced by SaveToRootFile().

 {
   if(!res::IsN(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,6);
 }

double GiBUURESFormFactor::FormFactors::C6A	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 278 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsN().

Referenced by SaveToRootFile().

 {
   if(!res::IsN(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,11);
 }

double GiBUURESFormFactor::FormFactors::C6V	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 250 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsN().

Referenced by SaveToRootFile().

 {
   if(!res::IsN(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,7);
 }

double GiBUURESFormFactor::FormFactors::F1V	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 285 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsDelta().

Referenced by SaveToRootFile().

 {
   if(!res::IsDelta(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,0);
 }

double GiBUURESFormFactor::FormFactors::F2V	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 292 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsDelta().

Referenced by SaveToRootFile().

 {
   if(!res::IsDelta(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,1);
 }

double GiBUURESFormFactor::FormFactors::FA	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 299 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsDelta().

Referenced by SaveToRootFile().

 {
   if(!res::IsDelta(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,2);
 }

double GiBUURESFormFactor::FormFactors::FFRes	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it,
		int	ffid
	)		const

private

func to retrieve interpolated form factor values

Definition at line 313 of file GiBUURESFormFactor.cxx.

References genie::Spline::Evaluate(), genie::kIntEM, genie::kIntWeakCC, genie::kIntWeakNC, genie::kPdgNeutron, and genie::kPdgProton.

 {
   if(Q2 < fMinQ2 || Q2 > fMaxQ2) return 0.;
 
   int r = -1, i = -1, j = -1;
 
   if(ffresid<0 || ffresid >= kNFFRes) return 0.;
 
   r = (int)res;
   if(r<0 || r >= kNRes) return 0.;
 
   if      (it == kIntWeakCC) { i = 0; }
   else if (it == kIntWeakNC) { i = 1; }
   else if (it == kIntEM)     { i = 2; }
 
   if      (hit_nucleon_pdg == kPdgNeutron) { j = 0; }
   else if (hit_nucleon_pdg == kPdgProton ) { j = 1; }
 
   const Spline * spl = fFFRes[r][i][j][ffresid];
 
   if(!spl) return 0;
   else {
    return spl->Evaluate(Q2);
   }
 }

double GiBUURESFormFactor::FormFactors::FP	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 306 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsDelta().

Referenced by SaveToRootFile().

 {
   if(!res::IsDelta(res)) return 0.;
   return this->FFRes(Q2,res,hit_nucleon_pdg,it,3);
 }

void GiBUURESFormFactor::FormFactors::LoadTables ( void )

private

load all form factor data tables

Definition at line 104 of file GiBUURESFormFactor.cxx.

References genie::utils::res::IsDelta(), genie::utils::res::IsN(), LOG, pDEBUG, pINFO, and pNOTICE.

Referenced by genie::GiBUURESFormFactor::LoadTables().

 {
 // Loads hadronic x-section data
 
   for(int r=0; r<kNRes; r++) {
     for(int i=0; i<kNCurr; i++) {
       for(int j=0; j<kNHitNuc; j++) {
         for(int k=0; k<kNFFRes; k++) {
              fFFRes[r][i][j][k] = 0;
         }//k
       }//j
     }//i
   }//r
 
   string data_dir = string(gSystem->Getenv("GENIE")) +
                     string("/data/evgen/gibuu");
 
   LOG("GiBUURESFormFactor", pNOTICE) << "Loading GiBUU data from: " << data_dir;
 
   //
   // load resonance form factor data
   //
   for(int r=0; r<kNRes; r++) {
     for(int i=0; i<kNCurr; i++) {
       for(int j=0; j<kNHitNuc; j++) {
 
         bool skip = false;
 
         Resonance_t resonance = (Resonance_t)r;
 
         ostringstream datafile;
         datafile << data_dir << "/form_factors/";
 
         switch(r){
           case ( 0): datafile << "P33_1232"; break;
           case ( 1): datafile << "S11_1535"; break;
           case ( 2): datafile << "D13_1520"; break;
           case ( 3): datafile << "S11_1650"; break;
           case ( 5): datafile << "D15_1675"; break;
           case ( 6): datafile << "S31_1620"; break;
           case ( 7): datafile << "D33_1700"; break;
           case ( 8): datafile << "P11_1440"; break;
           case (10): datafile << "P13_1720"; break;
           case (11): datafile << "F15_1680"; break;
           case (12): datafile << "P31_1910"; break;
           case (14): datafile << "F35_1905"; break;
           case (15): datafile << "F37_1950"; break;
           default  : skip = true;
 
         }
         switch(i){
           case (0): datafile << "_CC"; break;
           case (1): datafile << "_NC"; break;
           case (2): datafile << "_EM"; break;
           default : skip = true;
         }
         switch(j){
           case (0): datafile << "_neutron"; break;
           case (1): datafile << "_proton";  break;
           default : skip = true;
         }
         datafile << "_FFres.dat";
 
         if(skip) continue;
 
         //-- Make sure that all data file exists
         assert( ! gSystem->AccessPathName(datafile.str().c_str()) );
 
         //-- Read the data and fill a tree
 
         //   The GiBUU files have the data organized in 9 columns:
         //              0    1      2      3      4      5      6      7      8
         //   I=3/2 res: Qs   C_3^V  C_4^V  C_5^V  C_6^V  C_3^A  C_4^A  C_5^A  C_6^A
         //   I=1/2 res: Qs   F_1^V  F_2^V  -----  -----  F_A    F_P    -----  ----
 
         TTree data_ffres;
         data_ffres.ReadFile(datafile.str().c_str(),
                             "Q2/D:f1/D:f2/D:f3/D:f4/D:f5/D:f6/D:f7/D:f8/D");
 
         LOG("GiBUURESFormFactor", pDEBUG)
            << "Number of data rows: " << data_ffres.GetEntries();
 
         //
         // I=3/2 resonances
         //
         if(res::IsDelta(resonance)) {
            fFFRes[r][i][j][0]  = new Spline(&data_ffres, "Q2:f1"); // F1V = f(Q2)
            fFFRes[r][i][j][1]  = new Spline(&data_ffres, "Q2:f2"); // F2V = f(Q2)
            fFFRes[r][i][j][2]  = new Spline(&data_ffres, "Q2:f5"); // FA  = f(Q2)
            fFFRes[r][i][j][3]  = new Spline(&data_ffres, "Q2:f6"); // FP  = f(Q2)
         } // Delta res
         else
         //
         // I=1/2 resonances
         //
         if(res::IsN(resonance)) {
            fFFRes[r][i][j][4]  = new Spline(&data_ffres, "Q2:f1"); // C3V = f(Q2)
            fFFRes[r][i][j][5]  = new Spline(&data_ffres, "Q2:f2"); // C4V = f(Q2)
            fFFRes[r][i][j][6]  = new Spline(&data_ffres, "Q2:f3"); // C5V = f(Q2)
            fFFRes[r][i][j][7]  = new Spline(&data_ffres, "Q2:f4"); // C6V = f(Q2)
            fFFRes[r][i][j][8]  = new Spline(&data_ffres, "Q2:f5"); // C3A = f(Q2)
            fFFRes[r][i][j][9]  = new Spline(&data_ffres, "Q2:f6"); // C4A = f(Q2)
            fFFRes[r][i][j][10] = new Spline(&data_ffres, "Q2:f7"); // C5A = f(Q2)
            fFFRes[r][i][j][11] = new Spline(&data_ffres, "Q2:f8"); // C6A = f(Q2)
         } //N res
 
       }//j
     }//i
   }//r
 
 
   for(int r=0; r<kNRes; r++) {
     for(int i=0; i<kNCurr; i++) {
       for(int j=0; j<kNHitNuc; j++) {
         for(int k=0; k<kNFFRes; k++) {
              if(fFFRes[r][i][j][k]) fFFRes[r][i][j][k]->YCanBeNegative(true);
         }//k
       }//j
     }//i
   }//r
 
   LOG("GiBUURESFormFactor", pINFO)
      << "Done loading all resonance form factor files...";
 }

double genie::GiBUURESFormFactor::FormFactors::Q2max ( void ) const

inline

Definition at line 71 of file GiBUURESFormFactor.h.

References fMaxQ2.

71 { return fMaxQ2; }

genie::GiBUURESFormFactor::FormFactors::fMaxQ2

static double fMaxQ2

max Q2 for which resonance f/f data are given

Definition: GiBUURESFormFactor.h:76

double genie::GiBUURESFormFactor::FormFactors::Q2min ( void ) const

inline

Definition at line 70 of file GiBUURESFormFactor.h.

References fMinQ2.

70 { return fMinQ2; }

genie::GiBUURESFormFactor::FormFactors::fMinQ2

static double fMinQ2

min Q2 for which resonance f/f data are given

Definition: GiBUURESFormFactor.h:75

Friends And Related Function Documentation

friend class GiBUURESFormFactor

friend

Definition at line 115 of file GiBUURESFormFactor.h.

Member Data Documentation

Spline* genie::GiBUURESFormFactor::FormFactors::fFFRes[kNRes][kNCurr][kNHitNuc][kNFFRes]

private

actual form factor data = f(Q2)

Definition at line 107 of file GiBUURESFormFactor.h.

double GiBUURESFormFactor::FormFactors::fMaxQ2 = 4.0

staticprivate

max Q2 for which resonance f/f data are given

Definition at line 76 of file GiBUURESFormFactor.h.

Referenced by Q2max().

double GiBUURESFormFactor::FormFactors::fMinQ2 = 0.0

staticprivate

min Q2 for which resonance f/f data are given

Definition at line 75 of file GiBUURESFormFactor.h.

Referenced by Q2min().

const int genie::GiBUURESFormFactor::FormFactors::kNCurr = 3

staticprivate

Definition at line 102 of file GiBUURESFormFactor.h.

const int genie::GiBUURESFormFactor::FormFactors::kNFFRes = 12

staticprivate

Definition at line 104 of file GiBUURESFormFactor.h.

const int genie::GiBUURESFormFactor::FormFactors::kNHitNuc = 2

staticprivate

Definition at line 103 of file GiBUURESFormFactor.h.

const int genie::GiBUURESFormFactor::FormFactors::kNRes = 18

staticprivate

Definition at line 101 of file GiBUURESFormFactor.h.

The documentation for this class was generated from the following files:

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation