genie::GiBUUData::FormFactors Class Reference

#include <GiBUUData.h>

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	GiBUUData

Detailed Description

Definition at line 45 of file GiBUUData.h.

Constructor & Destructor Documentation

GiBUUData::FormFactors::FormFactors

(

void

)

Definition at line 89 of file GiBUUData.cxx.

{

}

GiBUUData::FormFactors::~FormFactors

(

void

)

Definition at line 94 of file GiBUUData.cxx.

{
  if(!gAbortingInErr) {
    cout << "GiBUUData 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 GiBUUData::FormFactors::C3A	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 268 of file GiBUUData.cxx.

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

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

Definition at line 240 of file GiBUUData.cxx.

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

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

Definition at line 275 of file GiBUUData.cxx.

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

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

Definition at line 247 of file GiBUUData.cxx.

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

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

Definition at line 282 of file GiBUUData.cxx.

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

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

Definition at line 254 of file GiBUUData.cxx.

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

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

Definition at line 289 of file GiBUUData.cxx.

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

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

Definition at line 261 of file GiBUUData.cxx.

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

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

Definition at line 296 of file GiBUUData.cxx.

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

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

Definition at line 303 of file GiBUUData.cxx.

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

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

Definition at line 310 of file GiBUUData.cxx.

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

double GiBUUData::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 324 of file GiBUUData.cxx.

{
  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 GiBUUData::FormFactors::FP	(	double	Q2,
		Resonance_t	res,
		int	nucleon_pdg,
		InteractionType_t	it
	)		const

Definition at line 317 of file GiBUUData.cxx.

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

void GiBUUData::FormFactors::LoadTables ( void  )

private

load all form factor data tables

Definition at line 115 of file GiBUUData.cxx.

{
// 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("GiBUUData", 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("GiBUUData", 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("GiBUUData", pINFO)  
     << "Done loading all resonance form factor files...";
}

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

inline

Definition at line 69 of file GiBUUData.h.

{ return fMaxQ2; }

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

inline

Definition at line 68 of file GiBUUData.h.

{ return fMinQ2; }

Friends And Related Function Documentation

friend class GiBUUData

friend

Definition at line 113 of file GiBUUData.h.

Member Data Documentation

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

private

actual form factor data = f(Q2)

Definition at line 105 of file GiBUUData.h.

double GiBUUData::FormFactors::fMaxQ2 = 4.0

staticprivate

max Q2 for which resonance f/f data are given

Definition at line 74 of file GiBUUData.h.

double GiBUUData::FormFactors::fMinQ2 = 0.0

staticprivate

min Q2 for which resonance f/f data are given

Definition at line 73 of file GiBUUData.h.

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

staticprivate

Definition at line 100 of file GiBUUData.h.

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

staticprivate

Definition at line 102 of file GiBUUData.h.

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

staticprivate

Definition at line 101 of file GiBUUData.h.

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

staticprivate

Definition at line 99 of file GiBUUData.h.

---

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

• GENIE/src/GiBUU/GiBUUData.h
• GENIE/src/GiBUU/GiBUUData.cxx