genie::flux::GDk2NuFlux Class Reference

An implementation of the GENIE GFluxI interface ("flux driver") encapsulating reading/processing the "dk2nu" tree structure. More...

#include "/cvmfs/nova.opensciencegrid.org/externals/dk2nugenie/v01_08_00e/source/dk2nu/genie/GDk2NuFlux.h"

Inheritance diagram for genie::flux::GDk2NuFlux:

Public Member Functions
	GDk2NuFlux ()
	~GDk2NuFlux ()
const PDGCodeList &	FluxParticles (void)
	declare list of flux neutrinos that can be generated (for init. purposes) More...
double	MaxEnergy (void)
	declare the max flux neutrino energy that can be generated (for init. purposes) More...
bool	GenerateNext (void)
	generate the next flux neutrino (return false in err) More...
int	PdgCode (void)
	returns the flux neutrino pdg code More...
double	Weight (void)
	returns the flux neutrino weight (if any) More...
const TLorentzVector &	Momentum (void)
	returns the flux neutrino 4-momentum More...
const TLorentzVector &	Position (void)
	returns the flux neutrino 4-position (note: expect SI rather than physical units) More...
bool	End (void)
	true if no more flux nu's can be thrown (eg reaching end of beam sim ntuples) More...
long int	Index (void)
	returns corresponding index for current flux neutrino (e.g. for a flux ntuple returns the current entry number) More...
void	Clear (Option_t *opt)
	reset state variables based on opt More...
void	GenerateWeighted (bool gen_weighted)
	set whether to generate weighted or unweighted neutrinos More...
const bsim::NuChoice &	GetNuChoice (void)
const bsim::Dk2Nu &	GetDk2Nu (void)
const bsim::DkMeta &	GetDkMeta (void)
Long64_t	GetEntryNumber ()
	index in chain More...
double	GetDecayDist () const
	dist (user units) from dk to current pos More...
void	MoveToZ0 (double z0)
	move ray origin to user coord Z0 More...
virtual double	GetTotalExposure () const
double	POT_curr (void)
	current average POT (RWH?) More...
double	UsedPOTs (void) const
	of protons-on-target used More...
long int	NFluxNeutrinos (void) const
	number of flux neutrinos looped so far More...
double	SumWeight (void) const
	integrated weight for flux neutrinos looped so far More...
void	PrintCurrent (void)
	print current entry from leaves More...
void	PrintConfig ()
	print the current configuration More...
std::vector< std::string >	GetFileList ()
	list of files currently part of chain More...
virtual void	LoadBeamSimData (const std::vector< std::string > &filenames, const std::string &det_loc)
virtual void	GetBranchInfo (std::vector< std::string > &branchNames, std::vector< std::string > &branchClassNames, std::vector< void ** > &branchObjPointers)
virtual TTree *	GetMetaDataTree ()
void	SetTreeNames (std::string fname="dk2nuTree", std::string mname="dkmetaTree")
bool	LoadConfig (std::string cfg)
	load a named configuration More...
void	SetMaxEnergy (double Ev)
	specify max flx nu energy More...
void	SetGenWeighted (bool genwgt=false)
	toggle whether GenerateNext() returns weight=1 flux (initial default false) More...
void	SetEntryReuse (long int nuse=1)
	of times to use entry before moving to next More...
void	ScanForMaxWeight (void)
	scan for max flux weight (before generating unweighted flux neutrinos) More...
void	SetMaxWgtScan (double fudge=1.05, long int nentries=2500000)
void	SetMaxEFudge (double fudge=1.05)
void	SetMinMaxWeight (double minwgt)
void	SetMaxWeightFailModel (int i=0)
void	SetApplyWindowTiltWeight (bool apply=true)
void	SetLengthUnits (double user_units)
	Set units assumed by user. More...
double	LengthUnits (void) const
	Return user units. More...
void	SetBeamRotation (TRotation beamrot)
	< beam (0,0,0) relative to user frame, beam direction in user frame More...
void	SetBeamCenter (TVector3 beam0)
TRotation	GetBeamRotation () const
	rotation to apply from beam->user More...
TVector3	GetBeamCenter () const
	beam origin in user frame More...
bool	IsFluxSphere () const
	flat window or sphere More...
void	SetFluxWindow (TVector3 p1, TVector3 p2, TVector3 p3)
	3 points define a plane (by default in user coordinates) More...
void	GetFluxWindow (TVector3 &p1, TVector3 &p2, TVector3 &p3) const
	3 points define a plane in beam coordinate More...
void	SetFluxSphere (TVector3 center, double radius, bool inDetCoord=true)
	specification of a sphere More...
void	GetFluxSphere (TVector3 &center, double &radius, bool inDetCoord=true) const
	specification of a sphere More...
void	Beam2UserPos (const TLorentzVector &beamxyz, TLorentzVector &usrxyz) const
void	Beam2UserDir (const TLorentzVector &beamdir, TLorentzVector &usrdir) const
void	Beam2UserP4 (const TLorentzVector &beamp4, TLorentzVector &usrp4) const
void	User2BeamPos (const TLorentzVector &usrxyz, TLorentzVector &beamxyz) const
void	User2BeamDir (const TLorentzVector &usrdir, TLorentzVector &beamdir) const
void	User2BeamP4 (const TLorentzVector &usrp4, TLorentzVector &beamp4) const
TVector3	FluxWindowNormal () const
TVector3	SphereNormal () const
	GDk2NuFlux ()
	~GDk2NuFlux ()
const PDGCodeList &	FluxParticles (void)
	declare list of flux neutrinos that can be generated (for init. purposes) More...
double	MaxEnergy (void)
	declare the max flux neutrino energy that can be generated (for init. purposes) More...
bool	GenerateNext (void)
	generate the next flux neutrino (return false in err) More...
int	PdgCode (void)
	returns the flux neutrino pdg code More...
double	Weight (void)
	returns the flux neutrino weight (if any) More...
const TLorentzVector &	Momentum (void)
	returns the flux neutrino 4-momentum More...
const TLorentzVector &	Position (void)
	returns the flux neutrino 4-position (note: expect SI rather than physical units) More...
bool	End (void)
	true if no more flux nu's can be thrown (eg reaching end of beam sim ntuples) More...
long int	Index (void)
	returns corresponding index for current flux neutrino (e.g. for a flux ntuple returns the current entry number) More...
void	Clear (Option_t *opt)
	reset state variables based on opt More...
void	GenerateWeighted (bool gen_weighted)
	set whether to generate weighted or unweighted neutrinos More...
const bsim::NuChoice &	GetNuChoice (void)
const bsim::Dk2Nu &	GetDk2Nu (void)
const bsim::DkMeta &	GetDkMeta (void)
Long64_t	GetEntryNumber ()
	index in chain More...
double	GetDecayDist () const
	dist (user units) from dk to current pos More...
void	MoveToZ0 (double z0)
	move ray origin to user coord Z0 More...
virtual double	GetTotalExposure () const
double	POT_curr (void)
	current average POT (RWH?) More...
double	UsedPOTs (void) const
	of protons-on-target used More...
long int	NFluxNeutrinos (void) const
	number of flux neutrinos looped so far More...
double	SumWeight (void) const
	integrated weight for flux neutrinos looped so far More...
void	PrintCurrent (void)
	print current entry from leaves More...
void	PrintConfig ()
	print the current configuration More...
std::vector< std::string >	GetFileList ()
	list of files currently part of chain More...
virtual void	LoadBeamSimData (const std::vector< std::string > &filenames, const std::string &det_loc)
virtual void	GetBranchInfo (std::vector< std::string > &branchNames, std::vector< std::string > &branchClassNames, std::vector< void ** > &branchObjPointers)
virtual TTree *	GetMetaDataTree ()
void	SetTreeNames (std::string fname="dk2nuTree", std::string mname="dkmetaTree")
bool	LoadConfig (std::string cfg)
	load a named configuration More...
void	SetMaxEnergy (double Ev)
	specify max flx nu energy More...
void	SetGenWeighted (bool genwgt=false)
	toggle whether GenerateNext() returns weight=1 flux (initial default false) More...
void	SetEntryReuse (long int nuse=1)
	of times to use entry before moving to next More...
void	ScanForMaxWeight (void)
	scan for max flux weight (before generating unweighted flux neutrinos) More...
void	SetMaxWgtScan (double fudge=1.05, long int nentries=2500000)
void	SetMaxEFudge (double fudge=1.05)
void	SetMinMaxWeight (double minwgt)
void	SetMaxWeightFailModel (int i=0)
void	SetApplyWindowTiltWeight (bool apply=true)
void	SetLengthUnits (double user_units)
	Set units assumed by user. More...
double	LengthUnits (void) const
	Return user units. More...
void	SetBeamRotation (TRotation beamrot)
	< beam (0,0,0) relative to user frame, beam direction in user frame More...
void	SetBeamCenter (TVector3 beam0)
TRotation	GetBeamRotation () const
	rotation to apply from beam->user More...
TVector3	GetBeamCenter () const
	beam origin in user frame More...
bool	IsFluxSphere () const
	flat window or sphere More...
void	SetFluxWindow (TVector3 p1, TVector3 p2, TVector3 p3)
	3 points define a plane (by default in user coordinates) More...
void	GetFluxWindow (TVector3 &p1, TVector3 &p2, TVector3 &p3) const
	3 points define a plane in beam coordinate More...
void	SetFluxSphere (TVector3 center, double radius, bool inDetCoord=true)
	specification of a sphere More...
void	GetFluxSphere (TVector3 &center, double &radius, bool inDetCoord=true) const
	specification of a sphere More...
void	Beam2UserPos (const TLorentzVector &beamxyz, TLorentzVector &usrxyz) const
void	Beam2UserDir (const TLorentzVector &beamdir, TLorentzVector &usrdir) const
void	Beam2UserP4 (const TLorentzVector &beamp4, TLorentzVector &usrp4) const
void	User2BeamPos (const TLorentzVector &usrxyz, TLorentzVector &beamxyz) const
void	User2BeamDir (const TLorentzVector &usrdir, TLorentzVector &beamdir) const
void	User2BeamP4 (const TLorentzVector &usrp4, TLorentzVector &beamp4) const
TVector3	FluxWindowNormal () const
TVector3	SphereNormal () const
const char *	GetExposureUnits () const
	what units are returned by GetTotalExposure? More...
genie::flux::Exposure_t	GetExposureType () const
virtual void	LoadBeamSimData (const std::set< std::string > &filenames, const std::string &det_loc)
virtual void	LoadBeamSimData (const std::string &filename, const std::string &det_loc)
virtual void	SetXMLFileBase (std::string xmlbasename="")
virtual std::string	GetXMLFileBase () const
virtual void	SetFluxParticles (const PDGCodeList &particles)
	specify list of flux neutrino species More...
virtual void	SetUpstreamZ (double z0)
virtual void	SetNumOfCycles (long int ncycle)
	limit cycling through input files More...

Static Public Member Functions
static const char *	AsString (genie::flux::Exposure_t etype)
static genie::flux::Exposure_t	StringToEnum (const char *chars, int maxChar=0)

Protected Attributes
PDGCodeList *	fPdgCList
	list of neutrino pdg-codes to generate More...
PDGCodeList *	fPdgCListRej
	list of nu pdg-codes seen but rejected More...
std::string	fXMLbasename
	XML file that might hold config param_sets. More...
long int	fNCycles
	times to cycle through the ntuple(s) More...
long int	fICycle
double	fZ0

Private Member Functions
bool	GenerateNext_weighted (void)
void	Initialize (void)
void	SetDefaults (void)
void	CleanUp (void)
void	ResetCurrent (void)
void	AddFile (TTree *fluxtree, TTree *metatree, std::string fname)
void	CalcEffPOTsPerNu (void)
void	LoadDkMeta (void)
bool	GenerateNext_weighted (void)
void	Initialize (void)
void	SetDefaults (void)
void	CleanUp (void)
void	ResetCurrent (void)
void	AddFile (TTree *fluxtree, TTree *metatree, std::string fname)
void	CalcEffPOTsPerNu (void)
void	LoadDkMeta (void)

Private Attributes
double	fMaxEv
	maximum energy More...
bool	fEnd
	end condition reached More...
std::vector< std::string >	fNuFluxFilePatterns
	(potentially wildcarded) path(s) More...
std::string	fTreeNames [2]
	pair of names "dk2nuTree", "dkmetaTree" More...
TChain *	fNuFluxTree
	TTree // REF ONLY! More...
TChain *	fNuMetaTree
	TTree // REF ONLY! More...
bsim::Dk2Nu *	fCurDk2Nu
bsim::DkMeta *	fCurDkMeta
bsim::NuChoice *	fCurNuChoice
int	fNFiles
	number of files in chain More...
Long64_t	fNEntries
	number of flux ntuple entries More...
Long64_t	fIEntry
	current flux ntuple entry More...
Long64_t	fNuTot
	cummulative # of entries (=fNEntries) More...
Long64_t	fFilePOTs
	of protons-on-target represented by all files More...
std::map< int, int >	fJobToMetaIndex
	quick lookup from job# to meta chain More...
double	fWeight
	current neutrino weight, =1 if generating unweighted entries More...
double	fMaxWeight
	max flux neutrino weight in input file More...
double	fMinMaxWeight
	user set lower limit on estimate More...
double	fMaxWeightScan
	initial estimate from scan More...
double	fMaxWeightInit
	max of scan & minmaxweight More...
double	fMaxWeightMax
	if "frozen" this is what bump would given More...
double	fMaxWgtFudge
	fudge factor for estimating max wgt More...
long int	fMaxWgtEntries
	of entries in estimating max wgt More...
double	fMaxEFudge
	fudge factor for estmating max enu (0=> use fixed 120GeV) More...
long int	fMaxWgtExceeded
	track failures of estimate More...
int	fMaxWgtFailModel
	what to do ... 0=bump, 1=frozen, 2=abort More...
long int	fNUse
	how often to use same entry in a row More...
long int	fIUse
	current # of times an entry has been used More...
double	fSumWeight
	sum of weights for nus thrown so far More...
long int	fNNeutrinos
	number of flux neutrinos thrown so far More...
double	fEffPOTsPerNu
	what a entry is worth ... More...
double	fAccumPOTs
	POTs used so far. More...
bool	fGenWeighted
	does GenerateNext() give weights? More...
bool	fApplyTiltWeight
	wgt due to window normal not || beam More...
bool	fDetLocIsSet
	is a flux location (near/far) set? More...
double	fLengthUnits
	units for coord in user exchanges More...
double	fLengthScaleB2U
	scale factor beam (cm) –> user More...
double	fLengthScaleU2B
	scale factor beam user –> (cm) More...
TLorentzVector	fBeamZero
	beam origin in user coords More...
TLorentzRotation	fBeamRot
	rotation applied beam –> user coord More...
TLorentzRotation	fBeamRotInv
bool	fIsSphere
	doing this on a sphere rather than a flat window? More...
TVector3	fFluxWindowPtUser [3]
	user points of flux window More...
TLorentzVector	fFluxWindowBase
	base point for flux window - beam coord More...
TLorentzVector	fFluxWindowDir1
	extent for flux window (direction 1) More...
TLorentzVector	fFluxWindowDir2
	extent for flux window (direction 2) More...
double	fFluxWindowLen1
double	fFluxWindowLen2
TVector3	fFluxWindowNormal
	normal direction for flux window – beam coord More...
TVector3	fFluxSphereCenterUser
	center for flux sphere - user coords More...
TVector3	fFluxSphereCenterBeam
	center for flux sphere - beam coords More...
double	fFluxSphereRadius
	radius for flux sphere More...
TLorentzVector	fgX4dkvtx
	decay 4-position beam coord More...

Detailed Description

An implementation of the GENIE GFluxI interface ("flux driver") encapsulating reading/processing the "dk2nu" tree structure.

Author

Robert Hatcher <rhatcher fnal.gov> Fermi National Accelerator Laboratory

Nov 6, 2012

Copyright (c) 2012, GENIE Neutrino MC Generator Collaboration For the full text of the license visit http://copyright.genie-mc.org or see $GENIE/LICENSE

Definition at line 69 of file GDk2NuFlux.h.

Constructor & Destructor Documentation

GDk2NuFlux::GDk2NuFlux

(

)

Definition at line 130 of file GDk2NuFlux.cxx.

References Initialize().

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

  : GFluxExposureI(genie::flux::kPOTs)
// #endif
{
  this->Initialize();
}

GDk2NuFlux::~GDk2NuFlux

(

)

Definition at line 138 of file GDk2NuFlux.cxx.

References CleanUp().

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  this->CleanUp();
}

genie::flux::GDk2NuFlux::GDk2NuFlux

(

)

genie::flux::GDk2NuFlux::~GDk2NuFlux

(

)

Member Function Documentation

void GDk2NuFlux::AddFile ( TTree *  fluxtree, TTree *  metatree, std::string  fname  )

private

Definition at line 1192 of file GDk2NuFlux.cxx.

References dkmeta, fCurDk2Nu, fCurDkMeta, fCurNuChoice, fFilePOTs, plot_validation_datamc::fname, fNFiles, fNuFluxTree, fNuMetaTree, fNuTot, fTreeNames, LOG, nentries, pFATAL, pINFO, bsim::DkMeta::pots, and SLOG.

Referenced by LoadBeamSimData(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // Add a file to the chain

  ULong64_t nentries = ftree->GetEntries();

  // generally files will only have one meta data entry, but if they've
  // been combined (i.e. "hadd") there might be more than one.
  int nmeta = mtree->GetEntries();
  double potsum = 0;
  bsim::DkMeta* dkmeta = new bsim::DkMeta;
  mtree->SetBranchAddress("dkmeta",&dkmeta);
  for (int imeta = 0; imeta < nmeta; ++imeta ) {
    mtree->GetEntry(imeta);
    double potentry = dkmeta->pots;
    potsum += potentry;
  }
  delete dkmeta;

  // don't need these anymore
  delete ftree;
  delete mtree;

  // make sure the chains are defined and a branch object attached
  if ( ! fNuFluxTree ) {
    fNuFluxTree  = new TChain(fTreeNames[0].c_str());
    fNuMetaTree  = new TChain(fTreeNames[1].c_str());
    fCurDk2Nu    = new bsim::Dk2Nu;
    fCurDkMeta   = new bsim::DkMeta;
    fCurNuChoice = new bsim::NuChoice;
    fNuFluxTree->SetBranchAddress("dk2nu",&fCurDk2Nu);
    fNuMetaTree->SetBranchAddress("dkmeta",&fCurDkMeta);
  }

  // add the file to the chains
  int stat0 = fNuFluxTree->AddFile(fname.c_str());
  int stat1 = fNuMetaTree->AddFile(fname.c_str());

  LOG("Flux",pINFO)
    << "flux->AddFile() of " << nentries
    << " " << ((mtree)?"[+meta]":"[no-meta]")
    << " [status=" << stat0 << "," << stat1 << "]"
    << nentries << " (" << nmeta << ")"
    << " entries in file: " << fname;

  if ( stat0 != 1 || stat1 != 1 ) {
    SLOG("GDk2NuFlux", pFATAL) << "Add: \"" << fname << "\" failed";
  }

  fNuTot    += nentries;
  fFilePOTs += potsum;
  fNFiles++;

}

void genie::flux::GDk2NuFlux::AddFile ( TTree *  fluxtree, TTree *  metatree, std::string  fname  )

private

const char * genie::flux::GFluxExposureI::AsString ( genie::flux::Exposure_t  etype )

staticinherited

Definition at line 20 of file GFluxExposureI.cxx.

References genie::flux::kPOTs, genie::flux::kSeconds, and genie::flux::kUnknown.

Referenced by genie::flux::GFluxExposureI::GetExposureUnits().

  {
    switch (etype) {
    case kUnknown:  return "UnknownExposureUnits";   break;
    case kPOTs:     return "POT";                    break;
    case kSeconds:  return "Seconds";                break;
    default:        return "?UnknownExposureUnits?"; break;
    }
  }

void GDk2NuFlux::Beam2UserDir	(	const TLorentzVector &	beamdir,
		TLorentzVector &	usrdir
	)		const

Definition at line 1014 of file GDk2NuFlux.cxx.

References fBeamRot, and fLengthScaleB2U.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  usrdir = fLengthScaleB2U*(fBeamRot*beamdir);
}

void genie::flux::GDk2NuFlux::Beam2UserDir	(	const TLorentzVector &	beamdir,
		TLorentzVector &	usrdir
	)		const

void GDk2NuFlux::Beam2UserP4	(	const TLorentzVector &	beamp4,
		TLorentzVector &	usrp4
	)		const

Definition at line 1019 of file GDk2NuFlux.cxx.

References fBeamRot.

Referenced by GenerateNext_weighted(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  usrp4 = fBeamRot*beamp4;
}

void genie::flux::GDk2NuFlux::Beam2UserP4	(	const TLorentzVector &	beamp4,
		TLorentzVector &	usrp4
	)		const

void GDk2NuFlux::Beam2UserPos	(	const TLorentzVector &	beamxyz,
		TLorentzVector &	usrxyz
	)		const

Definition at line 1009 of file GDk2NuFlux.cxx.

References fBeamRot, fBeamZero, and fLengthScaleB2U.

Referenced by GenerateNext_weighted(), SetFluxSphere(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  usrxyz = fLengthScaleB2U*(fBeamRot*beamxyz) + fBeamZero;
}

void genie::flux::GDk2NuFlux::Beam2UserPos	(	const TLorentzVector &	beamxyz,
		TLorentzVector &	usrxyz
	)		const

void GDk2NuFlux::CalcEffPOTsPerNu ( void  )

private

Definition at line 481 of file GDk2NuFlux.cxx.

References e, fEffPOTsPerNu, fFilePOTs, fFluxWindowDir1, fFluxWindowDir2, fNEntries, fNuFluxTree, LOG, Mag(), pNOTICE, and pWARN.

Referenced by LoadBeamSimData(), SetFluxWindow(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // do this if flux window changes or # of files changes

  if (!fNuFluxTree) return;  // not yet fully configured

  // effpots = mc_pots * (wgtfunction-area) / window-area / wgt-max-est
  //   wgtfunction-area = pi * radius-det-element^2 = pi * (100.cm)^2

  // this should match what is used in the CalcEnuWgt()
  const double kRDET = 100.0;   // set to flux per 100 cm radius
  const double kRDET2 = kRDET * kRDET;
  double flux_area = fFluxWindowDir1.Vect().Cross(fFluxWindowDir2.Vect()).Mag();
  LOG("Flux",pNOTICE) << "in CalcEffPOTsPerNu, area = " << flux_area;

  if ( flux_area < 1.0e-30 ) {
    LOG("Flux", pWARN)
          << "CalcEffPOTsPerNu called with flux window area effectively zero";
    flux_area = 1;
  }
  double area_ratio = TMath::Pi() * kRDET2 / flux_area;
  fEffPOTsPerNu = area_ratio * ( (double)fFilePOTs / (double)fNEntries );
}

void genie::flux::GDk2NuFlux::CalcEffPOTsPerNu ( void  )

private

void genie::flux::GDk2NuFlux::CleanUp ( void  )

private

void GDk2NuFlux::CleanUp ( void  )

private

Definition at line 1178 of file GDk2NuFlux.cxx.

References fCurNuChoice, genie::flux::GFluxFileConfigI::fICycle, fIEntry, fIUse, genie::flux::GFluxFileConfigI::fNCycles, fNUse, genie::flux::GFluxFileConfigI::fPdgCList, genie::flux::GFluxFileConfigI::fPdgCListRej, LOG, and pNOTICE.

Referenced by ~GDk2NuFlux(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  LOG("Flux", pNOTICE) << "Cleaning up...";

  if ( fPdgCList )    delete fPdgCList;
  if ( fPdgCListRej ) delete fPdgCListRej;
  if ( fCurNuChoice ) delete fCurNuChoice;

  LOG("Flux", pNOTICE)
    << " flux file cycles: " << fICycle << " of " << fNCycles
    << ", entry " << fIEntry << " use: " << fIUse << " of " << fNUse;
}

void genie::flux::GDk2NuFlux::Clear ( Option_t *  opt )

virtual

reset state variables based on opt

Implements genie::GFluxI.

void GDk2NuFlux::Clear ( Option_t *  opt )

virtual

reset state variables based on opt

Implements genie::GFluxI.

Definition at line 1049 of file GDk2NuFlux.cxx.

References fAccumPOTs, genie::flux::GFluxFileConfigI::fICycle, fNNeutrinos, fSumWeight, LOG, and pWARN.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // Clear the driver state
  //
  LOG("Flux", pWARN) << "GSimpleNtpFlux::Clear(" << opt << ") called";
  // do it in all cases, but EVGDriver/GMCJDriver will pass "CycleHistory"

  fICycle     = 0;

  fSumWeight  = 0;
  fNNeutrinos = 0;
  fAccumPOTs  = 0;
}

bool genie::flux::GDk2NuFlux::End ( void  )

inlinevirtual

true if no more flux nu's can be thrown (eg reaching end of beam sim ntuples)

Implements genie::GFluxI.

Definition at line 89 of file GDk2NuFlux.h.

Referenced by GenerateNext().

{ return  fEnd;                 }

bool genie::flux::GDk2NuFlux::End ( void  )

inlinevirtual

true if no more flux nu's can be thrown (eg reaching end of beam sim ntuples)

Implements genie::GFluxI.

Definition at line 89 of file GDk2NuFlux.h.

{ return  fEnd;                 }

const PDGCodeList& genie::flux::GDk2NuFlux::FluxParticles ( void  )

inlinevirtual

declare list of flux neutrinos that can be generated (for init. purposes)

Implements genie::GFluxI.

Definition at line 82 of file GDk2NuFlux.h.

{ return *fPdgCList;            }

const PDGCodeList& genie::flux::GDk2NuFlux::FluxParticles ( void  )

inlinevirtual

declare list of flux neutrinos that can be generated (for init. purposes)

Implements genie::GFluxI.

Definition at line 82 of file GDk2NuFlux.h.

{ return *fPdgCList;            }

TVector3 genie::flux::GDk2NuFlux::FluxWindowNormal ( ) const

inline

Definition at line 237 of file GDk2NuFlux.h.

References make_static_page::AddFile(), plot_validation_datamc::fname, TMVAGlob::Initialize(), and string.

Referenced by GenerateNext_weighted().

{ return (fIsSphere) ? SphereNormal() : fFluxWindowNormal; }

TVector3 genie::flux::GDk2NuFlux::FluxWindowNormal ( ) const

inline

Definition at line 237 of file GDk2NuFlux.h.

References make_static_page::AddFile(), plot_validation_datamc::fname, TMVAGlob::Initialize(), and string.

{ return (fIsSphere) ? SphereNormal() : fFluxWindowNormal; }

bool genie::flux::GDk2NuFlux::GenerateNext ( void  )

virtual

generate the next flux neutrino (return false in err)

Implements genie::GFluxI.

bool GDk2NuFlux::GenerateNext ( void  )

virtual

generate the next flux neutrino (return false in err)

Implements genie::GFluxI.

Definition at line 166 of file GDk2NuFlux.cxx.

References bsim::NuChoice::AsString(), bsim::Dk2Nu::AsString(), om::cout, febshutoff_auto::end, End(), MakeMiniprodValidationCuts::f, fCurDk2Nu, fCurNuChoice, fGenWeighted, fMaxWeight, fMaxWeightMax, fMaxWgtExceeded, fMaxWgtFailModel, fMaxWgtFudge, fWeight, GenerateNext_weighted(), genie::RandomGen::Instance(), LOG, genie::utils::print::P4AsShortString(), bsim::NuChoice::p4NuBeam, bsim::NuChoice::p4NuUser, bsim::NuChoice::pdgNu, pERROR, pFATAL, pNOTICE, PrintConfig(), r(), generate_hists::rnd, genie::RandomGen::RndFlux(), Weight(), genie::utils::print::X4AsString(), bsim::NuChoice::x4NuBeam, and bsim::NuChoice::x4NuUser.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
// Get next (unweighted) flux ntuple entry on the specified detector location
//
  RandomGen* rnd = RandomGen::Instance();
  while ( true ) {
     // Check for end of flux ntuple
     bool end = this->End();
     if ( end ) {
       LOG("Flux", pNOTICE) << "GenerateNext signaled End() ";
       return false;
     }

     // Get next weighted flux ntuple entry
     bool nextok = this->GenerateNext_weighted();
     if ( fGenWeighted ) return nextok;
     if ( ! nextok ) continue;

     /* RWH - debug purposes
     if ( fNCycles == 0 ) {
       LOG("Flux", pNOTICE)
          << "Got flux entry: " << fIEntry
          << " - Cycle: "<< fICycle << "/ infinite";
     } else {
       LOG("Flux", pNOTICE)
          << "Got flux entry: "<< fIEntry
          << " - Cycle: "<< fICycle << "/"<< fNCycles;
     }
     */

     // Get fractional weight & decide whether to accept curr flux neutrino
     double f = this->Weight() / fMaxWeight;
     //LOG("Flux", pNOTICE)
     //   << "Curr flux neutrino fractional weight = " << f;
     if (f > 1.) {
       // Oh,dear me ... what the weight max _should_ have been
       double bumped = this->Weight() * fMaxWgtFudge;
       if ( bumped > fMaxWeightMax ) fMaxWeightMax = bumped;
       ++fMaxWgtExceeded;

       // what to do ...
       if ( fMaxWgtFailModel <= 0 ) {
         fMaxWeight = bumped;
         LOG("Flux", pERROR)
           << "** Fractional weight = " << f
           << " > 1 !! Bump fMaxWeight estimate to " << fMaxWeight << "\n"
           << fCurDk2Nu->AsString() << "\n" << fCurNuChoice->AsString();
         std::cout << std::flush;
       } else {
         LOG("Flux", pERROR)
           << "** Fractional weight = " << f
           << " > 1 !! Leave fMaxWeight frozen at " << fMaxWeight << " "
           << " not bumped to " << bumped << ", max so far: "
           << fMaxWeightMax << "\n"
           << fCurDk2Nu->AsString() << "\n" << fCurNuChoice->AsString();
         std::cout << std::flush;
         if ( fMaxWgtFailModel >= 2 ) {
           PrintConfig();
           std::cout << std::flush;
           LOG("Flux", pFATAL)
             << "User requested not to continue";
           abort();
         }
       }
     }
     double r = rnd->RndFlux().Rndm(); // (f < 1.) ? rnd->RndFlux().Rndm() : 0;
     bool accept = ( r < f );
     if ( accept ) {

#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
       LOG("Flux", pNOTICE)
         << "Generated beam neutrino: "
         << "\n pdg-code: " << fCurNuChoice->pdgNu
         << "\n p4: " << utils::print::P4AsShortString(&(fCurNuChoice->p4NuBeam))
         << "\n x4: " << utils::print::X4AsString(&(fCurNuChoice->x4NuBeam))
         << "\n p4: " << utils::print::P4AsShortString(&(fCurNuChoice->p4NuUser))
         << "\n x4: " << utils::print::X4AsString(&(fCurNuChoice->x4NuUser));
#endif

       fWeight = 1.;
       return true;
     }

     //LOG("Flux", pNOTICE)
     //  << "** Rejecting current flux neutrino based on the flux weight only";
  }
  return false;
}

bool GDk2NuFlux::GenerateNext_weighted ( void  )

private

user might modify list via SetFluxParticles() in order to reject certain flavors, even if they're found in the file. So don't make a big fuss. Spit out a single message and then stop reporting that flavor as problematic.

Definition at line 255 of file GDk2NuFlux.cxx.

References bsim::Dk2Nu::ancestor, ana::assert(), Beam2UserP4(), Beam2UserPos(), bsim::calcEnuWgt(), bsim::Dk2Nu::decay, genie::PDGCodeList::ExistsInPDGCodeList(), fAccumPOTs, fApplyTiltWeight, fCurDk2Nu, fCurNuChoice, fEffPOTsPerNu, fEnd, fFluxWindowBase, fFluxWindowDir1, fFluxWindowDir2, fgX4dkvtx, genie::flux::GFluxFileConfigI::fICycle, fIEntry, fIUse, FluxWindowNormal(), fMaxEv, fMaxWeight, genie::flux::GFluxFileConfigI::fNCycles, fNEntries, fNNeutrinos, fNuFluxTree, genie::flux::GFluxFileConfigI::fPdgCList, genie::flux::GFluxFileConfigI::fPdgCListRej, fSumWeight, fWeight, genie::flux::GFluxFileConfigI::fZ0, GetDecayDist(), bsim::NuChoice::impWgt, bsim::Dk2Nu::indxnu(), genie::RandomGen::Instance(), LOG, MoveToZ0(), bsim::Decay::nimpwt, bsim::Decay::ntype, bsim::Dk2Nu::overflow(), genie::utils::print::P4AsShortString(), bsim::NuChoice::p4NuBeam, bsim::NuChoice::p4NuUser, pDEBUG, bsim::NuChoice::pdgNu, pERROR, pFATAL, pINFO, pNOTICE, bsim::Dk2Nu::potnum, genie::PDGCodeList::push_back(), pWARN, ResetCurrent(), generate_hists::rnd, genie::RandomGen::RndFlux(), Unit(), bsim::Decay::vx, bsim::Decay::vy, bsim::Decay::vz, Weight(), genie::utils::print::X4AsString(), bsim::NuChoice::x4NuBeam, bsim::NuChoice::x4NuUser, and bsim::NuChoice::xyWgt.

Referenced by GenerateNext(), ScanForMaxWeight(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
// Get next (weighted) flux ntuple entry on the specified detector location
//

  // Check whether a flux ntuple has been loaded
  if ( ! fNuFluxTree ) {
     LOG("Flux", pERROR)
          << "The flux driver has not been properly configured";
     return false;
  }

  // Reuse an entry?
  //std::cout << " ***** iuse " << fIUse << " nuse " << fNUse
  //          << " ientry " << fIEntry << " nentry " << fNEntries
  //          << " icycle " << fICycle << " ncycle " << fNCycles << std::endl;
  if ( fIUse < fNUse && fIEntry >= 0 ) {
    // Reuse this entry
    fIUse++;
  } else {
    // Reset previously generated neutrino code / 4-p / 4-x
    this->ResetCurrent();
    // Move on, read next flux ntuple entry
    fIEntry++;
    if ( fIEntry >= fNEntries ) {
      // Ran out of entries @ the current cycle of this flux file
      // Check whether more (or infinite) number of cycles is requested
      if (fICycle < fNCycles || fNCycles == 0 ) {
        fICycle++;
        fIEntry=0;
      } else {
        LOG("Flux", pWARN)
          << "No more entries in input flux neutrino ntuple, cycle "
          << fICycle << " of " << fNCycles;
        fEnd = true;
        //assert(0);
        return false;
      }
    }

    fNuFluxTree->GetEntry(fIEntry);

#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("Flux",pDEBUG)
    << "got " << fNNeutrinos << " new fIEntry " << fIEntry
    << " pot# " << fCurDk2Nu->potnum;
#endif

    fIUse = 1;

    // here we might want to do flavor oscillations or simple mappings
    fCurNuChoice->pdgNu  = fCurDk2Nu->decay.ntype;
    fCurNuChoice->impWgt = fCurDk2Nu->decay.nimpwt;
  }

  // update the # POTs & number of neutrinos
  // Do this HERE (before rejecting flavors that users might be weeding out)
  // in order to keep the POT accounting correct.  This allows one to get
  // the right normalization for generating only events from the intrinsic
  // nu_e entries.
  fAccumPOTs += fEffPOTsPerNu / fMaxWeight;
  fNNeutrinos++;

  // Check neutrino pdg against declared list of neutrino species declared
  // by the current instance of the neutrino flux driver.
  // No undeclared neutrino species will be accepted at this point as GENIE
  // has already been configured to handle the specified list.
  // Make sure that the appropriate list of flux neutrino species was set at
  // initialization via GDk2NuFlux::SetFluxParticles(const PDGCodeList &)

  if ( ! fPdgCList->ExistsInPDGCodeList(fCurNuChoice->pdgNu) ) {
     /// user might modify list via SetFluxParticles() in order to reject certain
     /// flavors, even if they're found in the file.  So don't make a big fuss.
     /// Spit out a single message and then stop reporting that flavor as problematic.
     int badpdg = fCurNuChoice->pdgNu;
     if ( ! fPdgCListRej->ExistsInPDGCodeList(badpdg) ) {
       fPdgCListRej->push_back(badpdg);
       LOG("Flux", pWARN)
         << "Encountered neutrino specie (" << badpdg
         << " that wasn't in SetFluxParticles() list, "
         << "\nDeclared list of neutrino species: " << *fPdgCList;
     }
     return false;
  }

  // Update the curr neutrino weight and energy

  // Check current neutrino energy against the maximum flux neutrino energy
  // declared by the current instance of the neutrino flux driver.
  // No flux neutrino exceeding that maximum energy will be accepted at this
  // point as that maximum energy has already been used for normalizing the
  // interaction probabilities.
  // Make sure that the appropriate maximum flux neutrino energy was set at
  // initialization via GDk2NuFlux::SetMaxEnergy(double Ev)

  fCurNuChoice->x4NuBeam = fFluxWindowBase;

  double Ev = 0;
  double& wgt_xy = fCurNuChoice->xyWgt;
  // recalculate on x-y window
  RandomGen * rnd = RandomGen::Instance();
  fCurNuChoice->x4NuBeam += ( rnd->RndFlux().Rndm()*fFluxWindowDir1 +
                              rnd->RndFlux().Rndm()*fFluxWindowDir2   );
  bsim::calcEnuWgt(fCurDk2Nu->decay,fCurNuChoice->x4NuBeam.Vect(),Ev,wgt_xy);

  if (Ev > fMaxEv) {
     LOG("Flux", pWARN)
          << "Flux neutrino energy exceeds declared maximum neutrino energy"
          << "\nEv = " << Ev << "(> Ev{max} = " << fMaxEv << ")";
  }

  // Set the current flux neutrino 4-momentum
  // this is in *beam* coordinates
  fgX4dkvtx = TLorentzVector( fCurDk2Nu->decay.vx,
                              fCurDk2Nu->decay.vy,
                              fCurDk2Nu->decay.vz, 0.);
  // don't use TLorentzVector here for Mag() due to - on metric
  TVector3 dirNu = (fCurNuChoice->x4NuBeam.Vect() - fgX4dkvtx.Vect()).Unit();
  fCurNuChoice->p4NuBeam.SetPxPyPzE( Ev*dirNu.X(),
                                     Ev*dirNu.Y(),
                                     Ev*dirNu.Z(), Ev);

  // Set the current flux neutrino 4-position, direction in user coord
  Beam2UserP4(fCurNuChoice->p4NuBeam,fCurNuChoice->p4NuUser);
  Beam2UserPos(fCurNuChoice->x4NuBeam,fCurNuChoice->x4NuUser);

  fWeight = fCurNuChoice->impWgt * fCurNuChoice->xyWgt;  // full weight
  if ( fApplyTiltWeight ) {
    // additional weight due to window tilt relative to flux direction
    double tiltwgt = dirNu.Dot( FluxWindowNormal() );
    fWeight *= TMath::Abs( tiltwgt );
  }

  // set the time component of the Lorentz vectors
  double dist_dk2start = GetDecayDist();
  size_t inu = fCurDk2Nu->indxnu();
  double t_dk = 0;
  if ( ! fCurDk2Nu->overflow() && ! fCurDk2Nu->ancestor.empty() ) {
    t_dk = fCurDk2Nu->ancestor[inu].startt; // units?  seconds, hopefully
  }
  const double c_mbys  = 299792458;
  const double c_cmbys = c_mbys * 100.;
  double tstart = t_dk + ( dist_dk2start / c_cmbys );
  fCurNuChoice->x4NuBeam.SetT(tstart);
  fCurNuChoice->x4NuUser.SetT(tstart);
  if ( t_dk == 0 ) {
    // probably wasn't set
    static int nmsg = 5;
    if ( nmsg > 0 ) {
      --nmsg;
      LOG("Flux", pNOTICE)
        << "Setting time at flux window, \n"
        << "noticed that t_dk from ancestor list was 0, "
        << "this probably means that the calculated time is wrong";
    }
  }

  // if desired, move to user specified user coord z
  if ( TMath::Abs(fZ0) < 1.0e30 ) this->MoveToZ0(fZ0);

#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("Flux", pINFO)
    << "Generated neutrino: " << fIEntry << " " << fCurDk2Nu->potnum
    << "\n pdg-code: " << fCurNuChoice->pdgNu
    << "\n p4 beam: " << utils::print::P4AsShortString(&fCurNuChoice->p4NuBeam)
    << "\n x4 beam: " << utils::print::X4AsString(&fCurNuChoice->x4NuBeam)
    << "\n p4 user: " << utils::print::P4AsShortString(&(fCurNuChoice->p4NuUser))
    << "\n x4 user: " << utils::print::X4AsString(&(fCurNuChoice->x4NuUser));
#endif
  if ( Ev > fMaxEv ) {
    LOG("Flux", pFATAL)
      << "Generated neutrino had E_nu = " << Ev << " > " << fMaxEv
      << " maximum ";
    assert(0);
  }

  // update sum of weights
  fSumWeight += this->Weight();

  return true;
}

bool genie::flux::GDk2NuFlux::GenerateNext_weighted ( void  )

private

void genie::flux::GDk2NuFlux::GenerateWeighted ( bool  gen_weighted )

virtual

set whether to generate weighted or unweighted neutrinos

Implements genie::GFluxI.

void GDk2NuFlux::GenerateWeighted ( bool  gen_weighted )

virtual

set whether to generate weighted or unweighted neutrinos

Implements genie::GFluxI.

Definition at line 1063 of file GDk2NuFlux.cxx.

References fGenWeighted.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // Set whether to generate weighted rays
  //
  fGenWeighted = gen_weighted;
}

TVector3 GDk2NuFlux::GetBeamCenter

(

)

const

beam origin in user frame

Definition at line 971 of file GDk2NuFlux.cxx.

References fBeamZero.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  TVector3 beam0 = fBeamZero.Vect();
  return beam0;
}

TVector3 genie::flux::GDk2NuFlux::GetBeamCenter

(

)

const

beam origin in user frame

TRotation GDk2NuFlux::GetBeamRotation

(

)

const

rotation to apply from beam->user

Definition at line 958 of file GDk2NuFlux.cxx.

References fBeamRot.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // rotation is really only 3-d vector, but we'll be operating on LorentzV's
  // give people back the original TRotation ... not pretty
  // ... it think this is right
  TRotation rot3;
  const TLorentzRotation& rot4 = fBeamRot;
  TVector3 newX(rot4.XX(),rot4.XY(),rot4.XZ());
  TVector3 newY(rot4.YX(),rot4.YY(),rot4.YZ());
  TVector3 newZ(rot4.ZX(),rot4.ZY(),rot4.ZZ());
  rot3.RotateAxes(newX,newY,newZ);
  return rot3.Inverse();
}

TRotation genie::flux::GDk2NuFlux::GetBeamRotation

(

)

const

rotation to apply from beam->user

virtual void genie::flux::GDk2NuFlux::GetBranchInfo ( std::vector< std::string > &  branchNames, std::vector< std::string > &  branchClassNames, std::vector< void ** > &  branchObjPointers  )

virtual

allow caller to copy current status / ntuple entry info in the output file by providing copies of internal info

Assumes that branch object pointers will not change which may require either a copy be made or, if using the class directly for reading the branch, one must force ROOT to not autodelete: myns::MyClassType* fCurrMyClass = new myns::MyClassType; myTree->SetBranchAddress("bname",&fCurMyClass); //? TBranch* b = myTree->GetBranch("bname"); //? b->SetAutoDelete(false);

ensure vectors are sized sufficiently (or use .push_back()) branchNames[i] = "bname" branchClassNames[i] = "myns::MyClassType" branchObjPointers[i] = (void**)

Reimplemented from genie::flux::GFluxFileConfigI.

void GDk2NuFlux::GetBranchInfo ( std::vector< std::string > &  branchNames, std::vector< std::string > &  branchClassNames, std::vector< void ** > &  branchObjPointers  )

virtual

allow caller to copy current status / ntuple entry info in the output file by providing copies of internal info

Assumes that branch object pointers will not change which may require either a copy be made or, if using the class directly for reading the branch, one must force ROOT to not autodelete: myns::MyClassType* fCurrMyClass = new myns::MyClassType; myTree->SetBranchAddress("bname",&fCurMyClass); //? TBranch* b = myTree->GetBranch("bname"); //? b->SetAutoDelete(false);

ensure vectors are sized sufficiently (or use .push_back()) branchNames[i] = "bname" branchClassNames[i] = "myns::MyClassType" branchObjPointers[i] = (void**)

Reimplemented from genie::flux::GFluxFileConfigI.

Definition at line 740 of file GDk2NuFlux.cxx.

References fCurDk2Nu, and fCurNuChoice.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // allow flux driver to report back current status and/or ntuple entry
  // info for possible recording in the output file by supplying
  // the class name, and a pointer to the object that will be filled
  // as well as a suggested name for the branch.

  branchNames.push_back("dk2nu");
  branchClassNames.push_back("bsim::Dk2Nu");
  branchObjPointers.push_back((void**)&fCurDk2Nu);

  branchNames.push_back("nuchoice");
  branchClassNames.push_back("bsim::NuChoice");
  branchObjPointers.push_back((void**)&fCurNuChoice);

}

double GDk2NuFlux::GetDecayDist

(

)

const

dist (user units) from dk to current pos

Definition at line 437 of file GDk2NuFlux.cxx.

References fCurNuChoice, fgX4dkvtx, fLengthScaleB2U, and bsim::NuChoice::x4NuBeam.

Referenced by GenerateNext_weighted(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // return distance (user units) between dk point and start position
  // these are in beam units
  TVector3 x3diff = fCurNuChoice->x4NuBeam.Vect() - fgX4dkvtx.Vect();
  return x3diff.Mag() * fLengthScaleB2U;
}

double genie::flux::GDk2NuFlux::GetDecayDist

(

)

const

dist (user units) from dk to current pos

const bsim::Dk2Nu& genie::flux::GDk2NuFlux::GetDk2Nu ( void  )

inline

Definition at line 104 of file GDk2NuFlux.h.

{ return *fCurDk2Nu; };

const bsim::Dk2Nu& genie::flux::GDk2NuFlux::GetDk2Nu ( void  )

inline

Definition at line 104 of file GDk2NuFlux.h.

Referenced by evgen::GENIEGen::produce().

{ return *fCurDk2Nu; };

const bsim::DkMeta& genie::flux::GDk2NuFlux::GetDkMeta ( void  )

inline

Definition at line 105 of file GDk2NuFlux.h.

{ LoadDkMeta(); return *fCurDkMeta; };

const bsim::DkMeta& genie::flux::GDk2NuFlux::GetDkMeta ( void  )

inline

Definition at line 105 of file GDk2NuFlux.h.

Referenced by evgen::GENIEGen::produce().

{ LoadDkMeta(); return *fCurDkMeta; };

Long64_t genie::flux::GDk2NuFlux::GetEntryNumber ( )

inline

index in chain

Definition at line 107 of file GDk2NuFlux.h.

Long64_t genie::flux::GDk2NuFlux::GetEntryNumber ( )

inline

index in chain

Definition at line 107 of file GDk2NuFlux.h.

genie::flux::Exposure_t genie::flux::GFluxExposureI::GetExposureType ( ) const

inherited

Definition at line 56 of file GFluxExposureI.cxx.

References genie::flux::GFluxExposureI::fEType.

  { return fEType; }

const char * genie::flux::GFluxExposureI::GetExposureUnits ( ) const

inherited

what units are returned by GetTotalExposure?

Definition at line 59 of file GFluxExposureI.cxx.

References genie::flux::GFluxExposureI::AsString(), and genie::flux::GFluxExposureI::fEType.

Referenced by main().

  { return genie::flux::GFluxExposureI::AsString(fEType); }

std::vector<std::string> genie::flux::GDk2NuFlux::GetFileList

(

)

list of files currently part of chain

std::vector< std::string > GDk2NuFlux::GetFileList

(

)

list of files currently part of chain

Definition at line 1410 of file GDk2NuFlux.cxx.

References modify_metadata_with_upmu_sample_type::flist, fNuFluxTree, and next().

Referenced by PrintConfig(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  std::vector<std::string> flist;
  if ( ! fNuFluxTree ) return flist;
  TObjArray *fileElements=fNuFluxTree->GetListOfFiles();
  if ( ! fileElements ) return flist;
  TIter next(fileElements);
  TChainElement *chEl=0;
  while (( chEl=(TChainElement*)next() )) {
    flist.push_back(chEl->GetTitle());
  }
  return flist;
}

void GDk2NuFlux::GetFluxSphere	(	TVector3 &	center,
		double &	radius,
		bool	inDetCoord = true
	)		const

specification of a sphere

Definition at line 926 of file GDk2NuFlux.cxx.

References fFluxSphereCenterBeam, fFluxSphereCenterUser, and fFluxSphereRadius.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  center = ( inDetCoord ) ? fFluxSphereCenterUser : fFluxSphereCenterBeam;
  radius = fFluxSphereRadius;
}

void genie::flux::GDk2NuFlux::GetFluxSphere	(	TVector3 &	center,
		double &	radius,
		bool	inDetCoord = true
	)		const

specification of a sphere

void GDk2NuFlux::GetFluxWindow	(	TVector3 &	p1,
		TVector3 &	p2,
		TVector3 &	p3
	)		const

3 points define a plane in beam coordinate

Definition at line 934 of file GDk2NuFlux.cxx.

References fFluxWindowPtUser.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // return flux window points
  p0 = fFluxWindowPtUser[0];
  p1 = fFluxWindowPtUser[1];
  p2 = fFluxWindowPtUser[2];

}

void genie::flux::GDk2NuFlux::GetFluxWindow	(	TVector3 &	p1,
		TVector3 &	p2,
		TVector3 &	p3
	)		const

3 points define a plane in beam coordinate

virtual TTree* genie::flux::GDk2NuFlux::GetMetaDataTree ( )

virtual

Reimplemented from genie::flux::GFluxFileConfigI.

TTree * GDk2NuFlux::GetMetaDataTree ( )

virtual

Reimplemented from genie::flux::GFluxFileConfigI.

Definition at line 758 of file GDk2NuFlux.cxx.

References fNuMetaTree.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{ return fNuMetaTree; }

const bsim::NuChoice& genie::flux::GDk2NuFlux::GetNuChoice ( void  )

inline

Definition at line 103 of file GDk2NuFlux.h.

{ return *fCurNuChoice; };

const bsim::NuChoice& genie::flux::GDk2NuFlux::GetNuChoice ( void  )

inline

Definition at line 103 of file GDk2NuFlux.h.

Referenced by evgen::GENIEGen::produce().

{ return *fCurNuChoice; };

double GDk2NuFlux::GetTotalExposure ( ) const

virtual

Implements genie::flux::GFluxExposureI.

Definition at line 144 of file GDk2NuFlux.cxx.

References UsedPOTs().

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // complete the GFluxExposureI interface
  return UsedPOTs();
}

virtual double genie::flux::GDk2NuFlux::GetTotalExposure ( ) const

virtual

Implements genie::flux::GFluxExposureI.

virtual std::string genie::flux::GFluxFileConfigI::GetXMLFileBase ( ) const

inlinevirtualinherited

Definition at line 54 of file GFluxFileConfigI.h.

References genie::flux::GFluxFileConfigI::fXMLbasename, genie::flux::GFluxFileConfigI::GetBranchInfo(), genie::flux::GFluxFileConfigI::GetMetaDataTree(), genie::flux::GFluxFileConfigI::PrintConfig(), genie::flux::GFluxFileConfigI::SetFluxParticles(), genie::flux::GFluxFileConfigI::SetNumOfCycles(), and genie::flux::GFluxFileConfigI::SetUpstreamZ().

{ return fXMLbasename; }

long int genie::flux::GDk2NuFlux::Index ( void  )

inlinevirtual

returns corresponding index for current flux neutrino (e.g. for a flux ntuple returns the current entry number)

Implements genie::GFluxI.

Definition at line 90 of file GDk2NuFlux.h.

References MECModelEnuComparisons::opt.

{ return  fIEntry;              }

long int genie::flux::GDk2NuFlux::Index ( void  )

inlinevirtual

returns corresponding index for current flux neutrino (e.g. for a flux ntuple returns the current entry number)

Implements genie::GFluxI.

Definition at line 90 of file GDk2NuFlux.h.

References MECModelEnuComparisons::opt.

{ return  fIEntry;              }

void GDk2NuFlux::Initialize ( void  )

private

Definition at line 1070 of file GDk2NuFlux.cxx.

References fAccumPOTs, fApplyTiltWeight, fCurDk2Nu, fCurDkMeta, fCurNuChoice, fDetLocIsSet, fEffPOTsPerNu, fEnd, fFilePOTs, fGenWeighted, genie::flux::GFluxFileConfigI::fICycle, fIEntry, fIsSphere, fIUse, fMaxEFudge, fMaxEv, fMaxWeight, fMaxWgtEntries, fMaxWgtExceeded, fMaxWgtFailModel, fMaxWgtFudge, fMinMaxWeight, fNEntries, fNFiles, fNNeutrinos, fNuFluxTree, fNuMetaTree, fNUse, fNuTot, fSumWeight, fTreeNames, LOG, pNOTICE, ResetCurrent(), SetDefaults(), SetLengthUnits(), and genie::utils::units::UnitFromString().

Referenced by GDk2NuFlux(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  LOG("Flux", pNOTICE) << "Initializing GDk2NuFlux driver";

  fMaxEv           =  0;
  fEnd             =  false;

// #if __GENIE_RELEASE_CODE__ <= GRELCODE(2,9,0)
//   fPdgCList        = new PDGCodeList;
//   fPdgCListRej     = new PDGCodeList;
// #endif

  fTreeNames[0]    = "dk2nuTree";
  fTreeNames[1]    = "dkmetaTree";
  fNuFluxTree      =  0;
  fNuMetaTree      =  0;
  fCurDk2Nu        =  0;
  fCurDkMeta       =  0;
  fCurNuChoice     =  0;
  fNFiles          =  0;

  fNEntries        =  0;
  fIEntry          = -1;
  fICycle          =  0;
  fNUse            =  1;
  fIUse            =  999999;

  fNuTot           = 0;
  fFilePOTs        = 0;

  fMaxWeight       = -1;
  fMinMaxWeight    =  0;
  fMaxWgtFudge     =  1.05;
  fMaxWgtEntries   = 2500000;
  fMaxEFudge       =  0;

  fMaxWgtExceeded  =  0;
  fMaxWgtFailModel =  0;  // default: 0=bump, 1=leave frozen, 2=abort

  fSumWeight       =  0;
  fNNeutrinos      =  0;
  fEffPOTsPerNu    =  0;
  fAccumPOTs       =  0;

  fGenWeighted     = false;
  fApplyTiltWeight = true;
  fIsSphere        = false;
  fDetLocIsSet     = false;
  // by default assume user length is m
  SetLengthUnits(genie::utils::units::UnitFromString("m"));

// #if __GENIE_RELEASE_CODE__ < GRELCODE(2,9,0)
//   // migrated to GFluxFileConfigI
//   fXMLbasename     = "GNuMIFlux.xml";
//   fNCycles         =  0;
//   fZ0              =  -3.4e38;
// #endif

  this->SetDefaults();
  this->ResetCurrent();
}

void genie::flux::GDk2NuFlux::Initialize ( void  )

private

bool genie::flux::GDk2NuFlux::IsFluxSphere ( ) const

inline

flat window or sphere

Definition at line 208 of file GDk2NuFlux.h.

References plot_validation_datamc::p1, plot_validation_datamc::p2, make_associated_cosmic_defs::p3, and radius.

Referenced by PrintConfig().

bool genie::flux::GDk2NuFlux::IsFluxSphere ( ) const

inline

flat window or sphere

Definition at line 208 of file GDk2NuFlux.h.

References plot_validation_datamc::p1, plot_validation_datamc::p2, make_associated_cosmic_defs::p3, and radius.

double GDk2NuFlux::LengthUnits

(

void

)

const

Return user units.

Definition at line 1282 of file GDk2NuFlux.cxx.

References Munits::cm, fLengthScaleB2U, and genie::utils::units::UnitFromString().

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // Return the scale factor for lengths the user is getting
  double cm = genie::utils::units::UnitFromString("cm");
  return fLengthScaleB2U * cm ;
}

double genie::flux::GDk2NuFlux::LengthUnits

(

void

)

const

Return user units.

void genie::flux::GFluxFileConfigI::LoadBeamSimData ( const std::set< std::string > &  filenames, const std::string &  det_loc  )

virtualinherited

Definition at line 34 of file GFluxFileConfigI.cxx.

References genie::flux::GFluxFileConfigI::LoadBeamSimData().

  {
    // Loads a beam simulation root file into the GFluxFileConfig driver.
    // have a set<> want a vector<>
    std::vector<std::string> filevec;
    std::copy(fileset.begin(),fileset.end(),std::back_inserter(filevec));
    LoadBeamSimData(filevec,config); // call the one that takes a vector
  }

void genie::flux::GFluxFileConfigI::LoadBeamSimData ( const std::string &  filename, const std::string &  det_loc  )

virtualinherited

Definition at line 45 of file GFluxFileConfigI.cxx.

References genie::flux::GFluxFileConfigI::LoadBeamSimData().

  {
    // Loads a beam simulation root file into the GFluxFileConfig driver.
    std::vector<std::string> filevec;
    filevec.push_back(filename);
    LoadBeamSimData(filevec,config); // call the one that takes a vector
  }

virtual void genie::flux::GDk2NuFlux::LoadBeamSimData ( const std::vector< std::string > &  filenames, const std::string &  det_loc  )

virtual

first is primary method for loading root flux ntuple files and config others are alternatives that can be overloaded but have sensible defaults to fall back to calling the vector version

Implements genie::flux::GFluxFileConfigI.

void GDk2NuFlux::LoadBeamSimData ( const std::vector< std::string > &  filenames, const std::string &  det_loc  )

virtual

first is primary method for loading root flux ntuple files and config others are alternatives that can be overloaded but have sensible defaults to fall back to calling the vector version

Implements genie::flux::GFluxFileConfigI.

Definition at line 572 of file GDk2NuFlux.cxx.

References AddFile(), CalcEffPOTsPerNu(), exit(), fAccumPOTs, fDetLocIsSet, genie::flux::GFluxFileConfigI::fICycle, fIEntry, shutoffs::filename, fIUse, fMaxWeight, fNEntries, fNNeutrinos, fNuFluxFilePatterns, fNuFluxTree, fNuMetaTree, fSumWeight, fTreeNames, genie::RandomGen::Instance(), LoadConfig(), LOG, submit_syst::pattern, pERROR, pFATAL, pINFO, pNOTICE, generate_hists::rnd, genie::RandomGen::RndFlux(), ScanForMaxWeight(), and string.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
// Loads in a beam simulation root file into the GDk2NuFlux driver.

  bool found_cfg = this->LoadConfig(config);
  if ( ! found_cfg ) {
    LOG("Flux", pFATAL)
      << "LoadBeamSimData could not find XML config \"" << config << "\"\n";
    exit(1);
  }

  fNuFluxFilePatterns = patterns;
  std::vector<int> nfiles_from_pattern;

  // create a (sorted) set of file names
  // this also helps ensure that the same file isn't listed multiple times
  std::set<std::string> fnames;

  LOG("Flux", pINFO) << "LoadBeamSimData was passed " << patterns.size()
                     << " patterns";

  for (size_t ipatt = 0; ipatt < patterns.size(); ++ipatt ) {
    std::string pattern = patterns[ipatt];
    nfiles_from_pattern.push_back(0);

    LOG("Flux", pNOTICE)
      << "Loading dk2nu flux tree from ROOT file pattern ["
      << std::setw(3) << ipatt << "] \"" << pattern << "\"";

    // !WILDCARD only works for file name ... NOT directory
    std::string dirname = gSystem->UnixPathName(gSystem->WorkingDirectory());
    size_t slashpos = pattern.find_last_of("/");
    size_t fbegin;
    if ( slashpos != std::string::npos ) {
      dirname = pattern.substr(0,slashpos);
      LOG("Flux", pINFO) << "Look for flux using directory " << dirname;
      fbegin = slashpos + 1;
    } else { fbegin = 0; }

    void* dirp = gSystem->OpenDirectory(gSystem->ExpandPathName(dirname.c_str()));
    if ( dirp ) {
      std::string basename =
      pattern.substr(fbegin,pattern.size()-fbegin);
      TRegexp re(basename.c_str(),kTRUE);
      const char* onefile;
      while ( ( onefile = gSystem->GetDirEntry(dirp) ) ) {
        TString afile = onefile;
        if ( afile=="." || afile==".." ) continue;
        if ( basename!=afile && afile.Index(re) == kNPOS ) continue;
        std::string fullname = dirname + "/" + afile.Data();
        fnames.insert(fullname);
        nfiles_from_pattern[ipatt]++;
      }
      gSystem->FreeDirectory(dirp);
    } // legal directory
  } // loop over patterns

  // std::cout << "RWH start adding files" << std::endl << std::flush;

  size_t indx = 0;
  std::set<std::string>::const_iterator sitr = fnames.begin();
  for ( ; sitr != fnames.end(); ++sitr, ++indx ) {
    std::string filename = *sitr;
    // std::cout << "RWH  [" << std::setw(3) << indx << "]  \""
    //           << filename << "\"" << std::endl << std::flush;
    bool isok = true;
    // this next test only works for local files, so we can't do that
    // if we want to stream via xrootd
    // ! (gSystem->AccessPathName(filename.c_str()));
    if ( isok ) {
      TFile* tf = TFile::Open(filename.c_str(),"READ");
      TTree* ftree = (TTree*)tf->Get(fTreeNames[0].c_str());
      TTree* mtree = (TTree*)tf->Get(fTreeNames[1].c_str());
      if ( ftree && mtree ) {
        // found both trees
        this->AddFile(ftree,mtree,filename);
        // std::cout << "RWH  AddFile(" << ftree << ","
        //    << mtree << "," << filename << ")" << std::endl << std::flush;
      } else {
        // std::cout << "RWH  no AddFile" << std::endl << std::flush;
        LOG("Flux", pNOTICE) << "File " << filename << " lacked a tree: "
                             << " \"" << fTreeNames[0] << "\" " << ftree
                             << " \"" << fTreeNames[1] << "\" " << mtree
                             << " and was not added to the file list";
      }
      tf->Close();
      delete tf;
    } // loop over tree type
  } // loop over sorted file names

  if ( ! fNuFluxTree || ! fNuMetaTree ) {
    LOG("Flux", pERROR)
      << "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!";
    LOG("Flux", pERROR)
      << "Missing dk2nu (" << fNuFluxTree << ") or dkmeta ("
      << fNuMetaTree << ") tree";
    LOG("Flux", pERROR)
      << "please check that input files are valid";
    LOG("Flux", pERROR)
      << "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!";
  }

  // this will open all files and read header!!
  fNEntries = fNuFluxTree->GetEntries();

  if ( fNEntries == 0 ) {
    LOG("Flux", pERROR)
      << "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!";
    LOG("Flux", pERROR)
      << "Loaded flux tree contains " <<  fNEntries << " entries";
    LOG("Flux", pERROR)
      << "Was passed the file patterns: ";
    for (size_t ipatt = 0; ipatt < patterns.size(); ++ipatt ) {
      std::string pattern = patterns[ipatt];
      LOG("Flux", pERROR)
        << "  [" << std::setw(3) << ipatt <<"] " << pattern;
    }
    LOG("Flux", pERROR)
      << "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!";
  } else {
    LOG("Flux", pNOTICE)
      << "Loaded flux tree contains " <<  fNEntries << " entries"
      << " from " << fnames.size() << " unique files";
    for (size_t ipatt = 0; ipatt < patterns.size(); ++ipatt ) {
      std::string pattern = patterns[ipatt];
      LOG("Flux", pINFO)
        << " pattern: " << pattern << " yielded "
        << nfiles_from_pattern[ipatt] << " files";
    }
  }

  TBranch* bdk2nu = fNuFluxTree->GetBranch("dk2nu");
  if (bdk2nu) bdk2nu->SetAutoDelete(false);
  TBranch* bdkmeta = fNuMetaTree->GetBranch("dkmeta");
  if (bdkmeta) bdkmeta->SetAutoDelete(false);

  // we have a file we can work with
  if (!fDetLocIsSet) {
     LOG("Flux", pERROR)
       << "LoadBeamSimData left detector location unset";
  }
  if (fMaxWeight<=0) {
     LOG("Flux", pINFO)
       << "Run ScanForMaxWeight() as part of LoadBeamSimData";
     this->ScanForMaxWeight();
  }

  // current ntuple cycle # (flux ntuples may be recycled)
  fICycle =  0;
  // pick a starting entry index [0:fNEntries-1]
  // pretend we just used up the the previous one
  RandomGen* rnd = RandomGen::Instance();
  fIUse   =  9999999;
  fIEntry = rnd->RndFlux().Integer(fNEntries) - 1;

  // don't count things we used to estimate max weight
  fSumWeight  = 0;
  fNNeutrinos = 0;
  fAccumPOTs  = 0;

  LOG("Flux",pNOTICE) << "about to CalcEffPOTsPerNu";
  this->CalcEffPOTsPerNu();

}

bool genie::flux::GDk2NuFlux::LoadConfig

(

std::string

cfg

)

load a named configuration

bool GDk2NuFlux::LoadConfig

(

std::string

cfg

)

load a named configuration

Definition at line 1291 of file GDk2NuFlux.cxx.

References genie::flux::GDk2NuFluxXMLHelper::LoadConfig().

Referenced by LoadBeamSimData(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  genie::flux::GDk2NuFluxXMLHelper helper(this);
  return helper.LoadConfig(cfg);
}

void genie::flux::GDk2NuFlux::LoadDkMeta ( void  )

private

void GDk2NuFlux::LoadDkMeta ( void  )

private

Definition at line 506 of file GDk2NuFlux.cxx.

References ana::assert(), fCurDk2Nu, fCurDkMeta, fJobToMetaIndex, fNuMetaTree, bsim::DkMeta::job, site_stats_from_log::job, bsim::Dk2Nu::job, LOG, pERROR, and pFATAL.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // load the matching bsim::dkmeta entry that goes w/ bsim::dk2nu entry

  if ( fJobToMetaIndex.empty() ) {
    int nmeta = fNuMetaTree->GetEntries();
    for (int imeta =0; imeta < nmeta; ++imeta ) {
      fNuMetaTree->GetEntry(imeta);
      int mjob = fCurDkMeta->job;
      // there shouldn't already be an entry in the map
      // complain if there is
      std::map<int,int>::const_iterator mitr = fJobToMetaIndex.find(mjob);
      if ( mitr == fJobToMetaIndex.end() ) {
        fJobToMetaIndex[mjob] = imeta;  // make an entry
      } else {
        LOG("Flux", pERROR) << "LoadDkMeta already had an entry for job "
                            << mjob << " at " << mitr->second
                            << " which conflicts with new entry at "
                            << imeta;
      }
    }
  }

  int job = fCurDk2Nu->job;
  if ( fCurDkMeta->job == job ) return;  // already loaded

  std::map<int,int>::const_iterator lookat = fJobToMetaIndex.find(job);
  if ( lookat != fJobToMetaIndex.end() ) {
    int indx = lookat->second;
    fNuMetaTree->GetEntry(indx);
    if ( fCurDkMeta->job != job ) {
      LOG("Flux", pFATAL) << "Failed to get right metadata " << job
                          << " => " << fCurDkMeta->job
                          << " indx " << indx;
      assert(0);
    }
  } else {
    // wasn't already indexed
    LOG("Flux", pFATAL) << "Failed index metadata for job " << job
                        << " indx " << lookat->second;
    assert(0);
  }

}

double genie::flux::GDk2NuFlux::MaxEnergy ( void  )

inlinevirtual

declare the max flux neutrino energy that can be generated (for init. purposes)

Implements genie::GFluxI.

Definition at line 83 of file GDk2NuFlux.h.

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

{ return  fMaxEv;               }

double genie::flux::GDk2NuFlux::MaxEnergy ( void  )

inlinevirtual

declare the max flux neutrino energy that can be generated (for init. purposes)

Implements genie::GFluxI.

Definition at line 83 of file GDk2NuFlux.h.

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

{ return  fMaxEv;               }

const TLorentzVector& genie::flux::GDk2NuFlux::Momentum ( void  )

virtual

returns the flux neutrino 4-momentum

Implements genie::GFluxI.

const TLorentzVector & GDk2NuFlux::Momentum ( void  )

virtual

returns the flux neutrino 4-momentum

Implements genie::GFluxI.

Definition at line 156 of file GDk2NuFlux.cxx.

References fCurNuChoice, and bsim::NuChoice::p4NuUser.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  return fCurNuChoice->p4NuUser;
}

void genie::flux::GDk2NuFlux::MoveToZ0

(

double

z0

)

move ray origin to user coord Z0

void GDk2NuFlux::MoveToZ0

(

double

z0

)

move ray origin to user coord Z0

Definition at line 445 of file GDk2NuFlux.cxx.

References e, fCurNuChoice, fLengthScaleU2B, LOG, genie::utils::print::P4AsShortString(), bsim::NuChoice::p4NuBeam, bsim::NuChoice::p4NuUser, pDEBUG, pWARN, scale, genie::utils::print::X4AsString(), bsim::NuChoice::x4NuBeam, and bsim::NuChoice::x4NuUser.

Referenced by GenerateNext_weighted(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // move ray origin to specified user z0
  // move beam coord entry correspondingly

#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("Flux", pDEBUG)
    << "MoveToZ0 (z0usr=" << z0usr << ") before:"
    << "\n p4 user: " << utils::print::P4AsShortString(&(fCurNuChoice->p4NuUser))
    << "\n x4 user: " << utils::print::X4AsString(&(fCurNuChoice->x4NuUser));
#endif

  double pzusr    = fCurNuChoice->p4NuUser.Pz();
  if ( TMath::Abs(pzusr) < 1.0e-30 ) {
    // neutrino is moving almost entirely in x-y plane
    LOG("Flux", pWARN)
      << "MoveToZ0(" << z0usr << ") not possible due to pz_usr (" << pzusr << ")";
    return;
  }

  double scale = (z0usr - fCurNuChoice->x4NuUser.Z()) / pzusr;
  fCurNuChoice->x4NuUser += (scale*fCurNuChoice->p4NuUser);
  fCurNuChoice->x4NuBeam += ((fLengthScaleU2B*scale)*fCurNuChoice->p4NuBeam);
  // this scaling works for distances, but not the time component
  fCurNuChoice->x4NuBeam.SetT(0);
  fCurNuChoice->x4NuUser.SetT(0);

#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("Flux", pDEBUG)
    << "MoveToZ0 (" << z0usr << ") after:"
    << "\n x4 user: " << utils::print::X4AsString(&(fCurNuChoice->x4NuUser));
#endif

}

long int genie::flux::GDk2NuFlux::NFluxNeutrinos ( void  ) const

inlinevirtual

number of flux neutrinos looped so far

Implements genie::flux::GFluxExposureI.

Definition at line 126 of file GDk2NuFlux.h.

long int genie::flux::GDk2NuFlux::NFluxNeutrinos ( void  ) const

inlinevirtual

number of flux neutrinos looped so far

Implements genie::flux::GFluxExposureI.

Definition at line 126 of file GDk2NuFlux.h.

int GDk2NuFlux::PdgCode ( void  )

virtual

returns the flux neutrino pdg code

Implements genie::GFluxI.

Definition at line 151 of file GDk2NuFlux.cxx.

References fCurNuChoice, and bsim::NuChoice::pdgNu.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  return fCurNuChoice->pdgNu;
}

int genie::flux::GDk2NuFlux::PdgCode ( void  )

virtual

returns the flux neutrino pdg code

Implements genie::GFluxI.

const TLorentzVector& genie::flux::GDk2NuFlux::Position ( void  )

virtual

returns the flux neutrino 4-position (note: expect SI rather than physical units)

Implements genie::GFluxI.

const TLorentzVector & GDk2NuFlux::Position ( void  )

virtual

returns the flux neutrino 4-position (note: expect SI rather than physical units)

Implements genie::GFluxI.

Definition at line 161 of file GDk2NuFlux.cxx.

References fCurNuChoice, and bsim::NuChoice::x4NuUser.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  return fCurNuChoice->x4NuUser;
}

double GDk2NuFlux::POT_curr

(

void

)

current average POT (RWH?)

Definition at line 565 of file GDk2NuFlux.cxx.

References UsedPOTs().

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

                                {
  // RWH: Not sure what POT_curr is supposed to represent I'll guess for
  // now that that it means what I mean by UsedPOTs().
  return UsedPOTs();
}

double genie::flux::GDk2NuFlux::POT_curr

(

void

)

current average POT (RWH?)

void genie::flux::GDk2NuFlux::PrintConfig ( )

virtual

print the current configuration

Implements genie::flux::GFluxFileConfigI.

void GDk2NuFlux::PrintConfig ( void  )

virtual

print the current configuration

Implements genie::flux::GFluxFileConfigI.

Definition at line 1299 of file GDk2NuFlux.cxx.

References fAccumPOTs, fApplyTiltWeight, fBeamRot, fBeamRotInv, fBeamZero, fDetLocIsSet, fEffPOTsPerNu, fFilePOTs, fFluxSphereCenterBeam, fFluxSphereCenterUser, fFluxSphereRadius, fFluxWindowBase, fFluxWindowDir1, fFluxWindowDir2, fFluxWindowLen1, fFluxWindowLen2, fFluxWindowPtUser, fGenWeighted, genie::flux::GFluxFileConfigI::fICycle, fIEntry, fIUse, fLengthScaleB2U, fLengthScaleU2B, fLengthUnits, modify_metadata_with_upmu_sample_type::flist, fMaxEv, fMaxWeight, fMaxWeightInit, fMaxWeightMax, fMaxWeightScan, fMaxWgtEntries, fMaxWgtExceeded, fMaxWgtFailModel, fMaxWgtFudge, fMinMaxWeight, genie::flux::GFluxFileConfigI::fNCycles, fNEntries, fNFiles, fNNeutrinos, fNuFluxFilePatterns, fNUse, genie::flux::GFluxFileConfigI::fPdgCList, genie::flux::GFluxFileConfigI::fPdgCListRej, fSumWeight, genie::flux::GFluxFileConfigI::fZ0, GetFileList(), MECModelEnuComparisons::i, IsFluxSphere(), LOG, pNOTICE, User2BeamPos(), genie::utils::print::Vec3AsString(), w, and genie::utils::print::X4AsString().

Referenced by GenerateNext(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{

  std::ostringstream s;
  PDGCodeList::const_iterator itr = fPdgCList->begin();
  for ( ; itr != fPdgCList->end(); ++itr) s << (*itr) << " ";
  s << "[rejected: ";
  itr = fPdgCListRej->begin();
  for ( ; itr != fPdgCListRej->end(); ++itr) s << (*itr) << " ";
  s << " ] ";

  std::ostringstream fpattout;
  for (size_t i = 0; i < fNuFluxFilePatterns.size(); ++i)
    fpattout << "\n [" << std::setw(3) << i << "] " << fNuFluxFilePatterns[i];

  std::ostringstream flistout;
  std::vector<std::string> flist = GetFileList();
  for (size_t i = 0; i < flist.size(); ++i)
    flistout << "\n [" << std::setw(3) << i << "] " << flist[i];

  TLorentzVector usr0(0,0,0,0);
  TLorentzVector usr0asbeam;
  User2BeamPos(usr0,usr0asbeam);

  const int w=10, p=6;
  std::ostringstream beamrot_str, beamrotinv_str;
  beamrot_str
    << "fBeamRot: " << std::setprecision(p) << "\n"
    << "  [ "
    << std::setw(w) << fBeamRot.XX() << " "
    << std::setw(w) << fBeamRot.XY() << " "
    << std::setw(w) << fBeamRot.XZ() << " ]\n"
    << "  [ "
    << std::setw(w) << fBeamRot.YX() << " "
    << std::setw(w) << fBeamRot.YY() << " "
    << std::setw(w) << fBeamRot.YZ() << " ]\n"
    << "  [ "
    << std::setw(w) << fBeamRot.ZX() << " "
    << std::setw(w) << fBeamRot.ZY() << " "
    << std::setw(w) << fBeamRot.ZZ() << " ]";
  beamrotinv_str
    << "fBeamRotInv: " << std::setprecision(p) << "\n"
    << "  [ "
    << std::setw(w) << fBeamRotInv.XX() << " "
    << std::setw(w) << fBeamRotInv.XY() << " "
    << std::setw(w) << fBeamRotInv.XZ() << " ]\n"
    << "  [ "
    << std::setw(w) << fBeamRotInv.YX() << " "
    << std::setw(w) << fBeamRotInv.YY() << " "
    << std::setw(w) << fBeamRotInv.YZ() << " ]\n"
    << "  [ "
    << std::setw(w) << fBeamRotInv.ZX() << " "
    << std::setw(w) << fBeamRotInv.ZY() << " "
    << std::setw(w) << fBeamRotInv.ZZ() << " ]";

  std::ostringstream config_str;
  config_str
    << "GDk2NuFlux Config:"
    << "\n Enu_max " << fMaxEv
    << "\n pdg-codes: " << s.str() << "\n "
    << fNEntries << " entries"
    << " (FilePOTs " << fFilePOTs << ") "
    <<  "in " << fNFiles << " files: "
    << flistout.str()
    << "\n from file patterns:"
    << fpattout.str()
    << "\n wgt max=" << fMaxWeight
    << " min=" << fMinMaxWeight << " scan=" << fMaxWeightScan
    << " ( fudge=" << fMaxWgtFudge << " using scan of "
    << fMaxWgtEntries << " entries )"
    << "\n exceeded init=" << fMaxWeightInit << ", N=" << fMaxWgtExceeded
    << " times (fail model " << fMaxWgtFailModel << ") max = " << fMaxWeightMax
    << "\n Z0 pushback " << fZ0
    << "\n used entry " << fIEntry << " " << fIUse << "/" << fNUse
    << " times, in " << fICycle << "/" << fNCycles << " cycles"
    << "\n SumWeight " << fSumWeight << " for " << fNNeutrinos << " neutrino entries"
    << "\n EffPOTsPerNu " << fEffPOTsPerNu << " AccumPOTs " << fAccumPOTs
    << "\n GenWeighted: \"" << (fGenWeighted?"true":"false") << "\", "
    << "ApplyTiltWeight: \"" << (fApplyTiltWeight?"true":"false") << "\", "
    << "Detector location set: \"" << (fDetLocIsSet?"true":"false") << "\", "
    << "\n LengthUnits " << fLengthUnits << ", scale b2u " << fLengthScaleB2U
    << ", scale u2b " << fLengthScaleU2B;

  if ( IsFluxSphere() ) {
    config_str
      << "\n Flux Sphere radius: " << fFluxSphereRadius
      << "\n   User Center: " << utils::print::Vec3AsString(&fFluxSphereCenterUser)
      << "\n   Beam Center: " << utils::print::Vec3AsString(&fFluxSphereCenterBeam);
  } else {
    config_str
      << "\n User Flux Window: "
      << "\n       " << utils::print::Vec3AsString(&(fFluxWindowPtUser[0]))
      << "\n       " << utils::print::Vec3AsString(&(fFluxWindowPtUser[1]))
      << "\n       " << utils::print::Vec3AsString(&(fFluxWindowPtUser[2]))
      << "\n Flux Window (cm, beam coord): "
      << "\n  base " << utils::print::X4AsString(&fFluxWindowBase)
      << "\n  dir1 " << utils::print::X4AsString(&fFluxWindowDir1) << " len " << fFluxWindowLen1
      << "\n  dir2 " << utils::print::X4AsString(&fFluxWindowDir2) << " len " << fFluxWindowLen2;
  }
  config_str
      << "\n User Beam Origin: "
      << "\n  base " << utils::print::X4AsString(&fBeamZero)
      << "\n " << beamrot_str.str() << " "
      << "\n Detector Origin (beam coord): "
      << "\n  base " << utils::print::X4AsString(&usr0asbeam)
      << "\n " << beamrotinv_str.str();

  LOG("Flux", pNOTICE) << config_str.str();
}

void genie::flux::GDk2NuFlux::PrintCurrent

(

void

)

print current entry from leaves

void GDk2NuFlux::PrintCurrent

(

void

)

print current entry from leaves

Definition at line 1042 of file GDk2NuFlux.cxx.

References bsim::NuChoice::AsString(), bsim::Dk2Nu::AsString(), fCurDk2Nu, fCurNuChoice, LOG, and pNOTICE.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  LOG("Flux", pNOTICE) << "CurrentEntry:\n"
                       << fCurDk2Nu->AsString() << "\n"
                       << fCurNuChoice->AsString();
}

void genie::flux::GDk2NuFlux::ResetCurrent ( void  )

private

void GDk2NuFlux::ResetCurrent ( void  )

private

Definition at line 1169 of file GDk2NuFlux.cxx.

References bsim::NuChoice::clear(), bsim::Dk2Nu::clear(), fCurDk2Nu, and fCurNuChoice.

Referenced by GenerateNext_weighted(), Initialize(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
// reset running values of neutrino pdg-code, 4-position & 4-momentum
// and the input ntuple leaves

  if ( fCurDk2Nu )    fCurDk2Nu->clear();
  if ( fCurNuChoice ) fCurNuChoice->clear();
}

void genie::flux::GDk2NuFlux::ScanForMaxWeight

(

void

)

scan for max flux weight (before generating unweighted flux neutrinos)

void GDk2NuFlux::ScanForMaxWeight

(

void

)

scan for max flux weight (before generating unweighted flux neutrinos)

Definition at line 792 of file GDk2NuFlux.cxx.

References enu, fCurNuChoice, fDetLocIsSet, fMaxEFudge, fMaxEv, fMaxWeight, fMaxWeightInit, fMaxWeightMax, fMaxWeightScan, fMaxWgtEntries, fMaxWgtExceeded, fMaxWgtFudge, fMinMaxWeight, GenerateNext_weighted(), LOG, bsim::NuChoice::p4NuBeam, pERROR, pNOTICE, confusionMatrixTree::t, Weight(), and wgt.

Referenced by LoadBeamSimData(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  if (!fDetLocIsSet) {
     LOG("Flux", pERROR)
       << "Specify a detector location before scanning for max weight";
     return;
  }

  // scan for the maximum weight

  double wgtgenmx = 0, enumx = 0;
  TStopwatch t;
  t.Start();
  for (int itry=0; itry < fMaxWgtEntries; ++itry) {
    this->GenerateNext_weighted();
    double wgt = this->Weight();
    if ( wgt > wgtgenmx ) wgtgenmx = wgt;
    double enu = fCurNuChoice->p4NuBeam.Energy();
    if ( enu > enumx ) enumx = enu;
  }
  t.Stop();
  t.Print("u");
  LOG("Flux", pNOTICE) << "Maximum flux weight for spin = "
                       << wgtgenmx << ", energy = " << enumx
                       << " (" << fMaxWgtEntries << ")";

  if (wgtgenmx > fMaxWeight ) fMaxWeight = wgtgenmx;
  // apply a fudge factor to estimated weight
  fMaxWeight *= fMaxWgtFudge;
  if ( fMaxWgtFudge < 1 ) fMaxWgtFudge=1; // avoid run-away downward spiral

  fMaxWeightScan = fMaxWeight;
  // apply possible user supplied minimum (floor)
  if (fMinMaxWeight > fMaxWeight) {
    LOG("Flux", pNOTICE) << "Maximum flux weight set by user = "
                         << fMinMaxWeight << " exceeds estimate, use that";
    fMaxWeight = fMinMaxWeight;
  }

  fMaxWeightInit  = fMaxWeight; // what we started with before generation
  fMaxWeightMax   = fMaxWeight;
  fMaxWgtExceeded = 0;          // not yet exceeded

  // adjust max energy?
  if ( enumx*fMaxEFudge > fMaxEv ) {
    LOG("Flux", pNOTICE) << "Adjust max: was=" << fMaxEv
                         << " now " << enumx << "*" << fMaxEFudge
                         << " = " << enumx*fMaxEFudge;
    fMaxEv = enumx * fMaxEFudge;
  }

  LOG("Flux", pNOTICE) << "Maximum flux weight = " << fMaxWeight
                       << ", energy = " << fMaxEv;

}

void genie::flux::GDk2NuFlux::SetApplyWindowTiltWeight ( bool  apply = true )

inline

Parameters

apply

apply wgt due to tilt of flux window relative to beam

Definition at line 179 of file GDk2NuFlux.h.

References g4zmq::apply().

            { fApplyTiltWeight = apply; }

void genie::flux::GDk2NuFlux::SetApplyWindowTiltWeight ( bool  apply = true )

inline

Parameters

apply

apply wgt due to tilt of flux window relative to beam

Definition at line 179 of file GDk2NuFlux.h.

References g4zmq::apply().

            { fApplyTiltWeight = apply; }

void genie::flux::GDk2NuFlux::SetBeamCenter

(

TVector3

beam0

)

void GDk2NuFlux::SetBeamCenter

(

TVector3

beam0

)

Definition at line 950 of file GDk2NuFlux.cxx.

References fBeamZero.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // set coord transform between detector and beam
  // NOTE: internally these are in "cm", but user might have set a preference
  fBeamZero = TLorentzVector(beam0,0);  // no time shift
}

void GDk2NuFlux::SetBeamRotation

(

TRotation

beamrot

)

< beam (0,0,0) relative to user frame, beam direction in user frame

Definition at line 943 of file GDk2NuFlux.cxx.

References fBeamRot, and fBeamRotInv.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // rotation is really only 3-d vector, but we'll be operating on LorentzV's
  fBeamRot    = TLorentzRotation(beamrot);
  fBeamRotInv = fBeamRot.Inverse();
}

void genie::flux::GDk2NuFlux::SetBeamRotation

(

TRotation

beamrot

)

< beam (0,0,0) relative to user frame, beam direction in user frame

void genie::flux::GDk2NuFlux::SetDefaults ( void  )

private

void GDk2NuFlux::SetDefaults ( void  )

private

Definition at line 1132 of file GDk2NuFlux.cxx.

References genie::kPdgAntiNuE, genie::kPdgAntiNuMu, genie::kPdgNuE, genie::kPdgNuMu, LOG, pNOTICE, genie::PDGCodeList::push_back(), SetEntryReuse(), genie::flux::GFluxFileConfigI::SetFluxParticles(), SetMaxEnergy(), genie::flux::GFluxFileConfigI::SetNumOfCycles(), genie::flux::GFluxFileConfigI::SetUpstreamZ(), genie::flux::GFluxFileConfigI::SetXMLFileBase(), and string.

Referenced by Initialize(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
// - Set default neutrino species list (nue, nuebar, numu, numubar) and
//   maximum energy (120 GeV).
//   These defaults can be overwritten by user calls (at the driver init) to
//   GDk2NuFlux::SetMaxEnergy(double Ev) and
//   GDk2NuFlux::SetFluxParticles(const PDGCodeList & particles)
// - Set the default file normalization to 1E+21 POT
// - Set the default flux neutrino start z position at -5m (z=0 is the
//   detector centre).
// - Set number of cycles to 1

  LOG("Flux", pNOTICE) << "Setting default GDk2NuFlux driver options";

  PDGCodeList particles;
  particles.push_back(kPdgNuMu);
  particles.push_back(kPdgAntiNuMu);
  particles.push_back(kPdgNuE);
  particles.push_back(kPdgAntiNuE);

  this->SetFluxParticles (particles);
  this->SetMaxEnergy     (120./*GeV*/);  // was 200, but that would be wasteful
  this->SetUpstreamZ     (-3.4e38); // way upstream ==> use flux window
  this->SetNumOfCycles   (0);
  this->SetEntryReuse    (1);

  std::string xmlfile = "GNuMIFlux.xml";
  const char* altxml = gSystem->Getenv("GDK2NUFLUXXML");
  if ( altxml ) xmlfile = altxml;
// #if __GENIE_RELEASE_CODE__ >= GRELCODE(2,9,0)
  this->SetXMLFileBase(xmlfile);
// #else
//   this->SetXMLFile(xmlfile);
// #endif

}

void genie::flux::GDk2NuFlux::SetEntryReuse

(

long int

nuse = 1

)

of times to use entry before moving to next

void GDk2NuFlux::SetEntryReuse

(

long int

nuse = 1

)

of times to use entry before moving to next

Definition at line 856 of file GDk2NuFlux.cxx.

References fNUse, and CLHEP::L.

Referenced by SetDefaults(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
// With nuse > 1 then the same entry in the file is used "nuse" times
// before moving on to the next entry in the ntuple

  fNUse    = TMath::Max(1L, nuse);
}

void genie::flux::GFluxFileConfigI::SetFluxParticles ( const PDGCodeList &  particles )

virtualinherited

specify list of flux neutrino species

Definition at line 100 of file GFluxFileConfigI.cxx.

References genie::PDGCodeList::Copy(), genie::flux::GFluxFileConfigI::fPdgCList, LOG, and pINFO.

Referenced by genie::flux::GFluxFileConfigI::GetXMLFileBase(), main(), and SetDefaults().

  {
    fPdgCList->Copy(particles);
    
    LOG("Flux", pINFO)
      << "Declared list of neutrino species: " << *fPdgCList;
  }

void genie::flux::GDk2NuFlux::SetFluxSphere	(	TVector3	center,
		double	radius,
		bool	inDetCoord = true
	)

specification of a sphere

void GDk2NuFlux::SetFluxSphere	(	TVector3	center,
		double	radius,
		bool	inDetCoord = true
	)

specification of a sphere

Definition at line 902 of file GDk2NuFlux.cxx.

References ana::assert(), Beam2UserPos(), fFluxSphereCenterBeam, fFluxSphereCenterUser, fFluxSphereRadius, fIsSphere, LOG, pFATAL, radius, and User2BeamPos().

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  fIsSphere = true;
  fFluxSphereRadius = radius;
  if ( inDetCoord ) {
    fFluxSphereCenterUser = center;
    TLorentzVector v4beam;
    User2BeamPos(TLorentzVector(center,0),v4beam);
    fFluxSphereCenterBeam = v4beam.Vect();
  } else {
    fFluxSphereCenterBeam = center;
    TLorentzVector v4user;
    Beam2UserPos(TLorentzVector(center,0),v4user);
    fFluxSphereCenterBeam = v4user.Vect();
  }

  // not yet supported
  LOG("Flux", pFATAL) << "use of spherical flux window not yet fully supported";
  assert(0);

}

void genie::flux::GDk2NuFlux::SetFluxWindow	(	TVector3	p1,
		TVector3	p2,
		TVector3	p3
	)

3 points define a plane (by default in user coordinates)

void GDk2NuFlux::SetFluxWindow	(	TVector3	p1,
		TVector3	p2,
		TVector3	p3
	)

3 points define a plane (by default in user coordinates)

Definition at line 865 of file GDk2NuFlux.cxx.

References CalcEffPOTsPerNu(), stan::math::dot(), e, fDetLocIsSet, fFluxWindowBase, fFluxWindowDir1, fFluxWindowDir2, fFluxWindowLen1, fFluxWindowLen2, fFluxWindowNormal, fFluxWindowPtUser, fIsSphere, LOG, PandAna.Demos.pi0_spectra::p0, plot_validation_datamc::p1, plot_validation_datamc::p2, pNOTICE, pWARN, Unit(), and User2BeamPos().

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // set flux window
  // NOTE: internally these are in "cm", but user might have set a preference
  fDetLocIsSet         = true;
  fIsSphere            = false;

  fFluxWindowPtUser[0] = p0;
  fFluxWindowPtUser[1] = p1;
  fFluxWindowPtUser[2] = p2;

  // convert from user to beam coord and from 3 points to base + 2 directions
  // apply units conversion
  TLorentzVector ptbm0, ptbm1, ptbm2;
  User2BeamPos(TLorentzVector(fFluxWindowPtUser[0],0),ptbm0);
  User2BeamPos(TLorentzVector(fFluxWindowPtUser[1],0),ptbm1);
  User2BeamPos(TLorentzVector(fFluxWindowPtUser[2],0),ptbm2);

  fFluxWindowBase = ptbm0;
  fFluxWindowDir1 = ptbm1 - ptbm0;
  fFluxWindowDir2 = ptbm2 - ptbm0;

  fFluxWindowLen1   = fFluxWindowDir1.Mag();
  fFluxWindowLen2   = fFluxWindowDir2.Mag();
  fFluxWindowNormal = fFluxWindowDir1.Vect().Cross(fFluxWindowDir2.Vect()).Unit();

  double dot = fFluxWindowDir1.Dot(fFluxWindowDir2);
  if ( TMath::Abs(dot) > 1.0e-8 )
    LOG("Flux",pWARN) << "Dot product between window direction vectors was "
                      << dot << "; please check for orthoganality";

  LOG("Flux",pNOTICE) << "about to CalcEffPOTsPerNu";
  this->CalcEffPOTsPerNu();
}

void genie::flux::GDk2NuFlux::SetGenWeighted ( bool  genwgt = false )

inline

toggle whether GenerateNext() returns weight=1 flux (initial default false)

Definition at line 163 of file GDk2NuFlux.h.

void genie::flux::GDk2NuFlux::SetGenWeighted ( bool  genwgt = false )

inline

toggle whether GenerateNext() returns weight=1 flux (initial default false)

Definition at line 163 of file GDk2NuFlux.h.

void GDk2NuFlux::SetLengthUnits

(

double

user_units

)

Set units assumed by user.

Definition at line 1248 of file GDk2NuFlux.cxx.

References Munits::cm, fBeamZero, fCurNuChoice, fFluxWindowPtUser, fLengthScaleB2U, fLengthScaleU2B, fLengthUnits, genie::utils::units::UnitFromString(), and bsim::NuChoice::x4NuUser.

Referenced by Initialize(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // Set the scale factor for lengths going from beam (cm) to user coordinates

  // GDk2NuFlux uses "cm" as the length unit consistently internally (this is
  // the length units used by both the g3 and g4 ntuples).  User interactions
  // setting the beam-to-detector coordinate transform, flux window, and the
  // returned position might need to be in other units.  Use:
  //     double scale = genie::utils::units::UnitFromString("cm");
  // ( #include "Utils/UnitUtils.h for declaration )
  // to get the correct scale factor to pass in.

  double rescale = fLengthUnits / user_units;
  fLengthUnits = user_units;
  double cm = genie::utils::units::UnitFromString("cm");
  fLengthScaleB2U = cm / user_units;
  fLengthScaleU2B = user_units / cm;

  // in case we're changing units without resetting transform/window
  // not recommended, but should work
  if (fCurNuChoice) fCurNuChoice->x4NuUser  *= rescale;
  fBeamZero               *= rescale;
  fFluxWindowPtUser[0]    *= rescale;
  fFluxWindowPtUser[1]    *= rescale;
  fFluxWindowPtUser[2]    *= rescale;

  // case GDk2NuFlux::kmeter:  fLengthScaleB2U =   0.01  ; break;
  // case GDk2NuFlux::kcm:     fLengthScaleB2U =   1.    ; break;
  // case GDk2NuFlux::kmm:     fLengthScaleB2U =  10.    ; break;
  // case GDk2NuFlux::kfm:     fLengthScaleB2U =   1.e13 ; break;  // 10e-2m -> 10e-15m

}

void genie::flux::GDk2NuFlux::SetLengthUnits

(

double

user_units

)

Set units assumed by user.

void genie::flux::GDk2NuFlux::SetMaxEFudge ( double  fudge = 1.05 )

inline

Parameters

fudge

extra fudge factor in estimating maximum energy

Definition at line 170 of file GDk2NuFlux.h.

            { fMaxEFudge = fudge; }

void genie::flux::GDk2NuFlux::SetMaxEFudge ( double  fudge = 1.05 )

inline

Parameters

fudge

extra fudge factor in estimating maximum energy

Definition at line 170 of file GDk2NuFlux.h.

            { fMaxEFudge = fudge; }

void genie::flux::GDk2NuFlux::SetMaxEnergy

(

double

Ev

)

specify max flx nu energy

void GDk2NuFlux::SetMaxEnergy

(

double

Ev

)

specify max flx nu energy

Definition at line 848 of file GDk2NuFlux.cxx.

References fMaxEv, LOG, and pINFO.

Referenced by SetDefaults(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  fMaxEv = TMath::Max(0.,Ev);

  LOG("Flux", pINFO)
    << "Declared maximum flux neutrino energy: " << fMaxEv;
}

void genie::flux::GDk2NuFlux::SetMaxWeightFailModel ( int  i = 0 )

inline

Parameters

i	0=bump, 1=leave frozen, 2=abort

Definition at line 175 of file GDk2NuFlux.h.

References MECModelEnuComparisons::i.

            { fMaxWgtFailModel = i; }

void genie::flux::GDk2NuFlux::SetMaxWeightFailModel ( int  i = 0 )

inline

Parameters

i	0=bump, 1=leave frozen, 2=abort

Definition at line 175 of file GDk2NuFlux.h.

References MECModelEnuComparisons::i.

            { fMaxWgtFailModel = i; }

void genie::flux::GDk2NuFlux::SetMaxWgtScan ( double  fudge = 1.05, long int  nentries = 2500000  )

inline

Parameters

nentries

configuration when estimating max weight

Definition at line 168 of file GDk2NuFlux.h.

References nentries.

            { fMaxWgtFudge = fudge; fMaxWgtEntries = nentries; }

void genie::flux::GDk2NuFlux::SetMaxWgtScan ( double  fudge = 1.05, long int  nentries = 2500000  )

inline

Parameters

nentries

configuration when estimating max weight

Definition at line 168 of file GDk2NuFlux.h.

References nentries.

            { fMaxWgtFudge = fudge; fMaxWgtEntries = nentries; }

void genie::flux::GDk2NuFlux::SetMinMaxWeight ( double  minwgt )

inline

Parameters

minwgt

user floor on estimating maxweight (to avoid bumping)

Definition at line 173 of file GDk2NuFlux.h.

            { fMinMaxWeight = minwgt; }

void genie::flux::GDk2NuFlux::SetMinMaxWeight ( double  minwgt )

inline

Parameters

minwgt

user floor on estimating maxweight (to avoid bumping)

Definition at line 173 of file GDk2NuFlux.h.

            { fMinMaxWeight = minwgt; }

void genie::flux::GFluxFileConfigI::SetNumOfCycles ( long int  ncycle )

virtualinherited

limit cycling through input files

Definition at line 86 of file GFluxFileConfigI.cxx.

References genie::flux::GFluxFileConfigI::fNCycles, and CLHEP::L.

Referenced by genie::flux::GFluxFileConfigI::GetXMLFileBase(), main(), SetDefaults(), and genie::flux::GSimpleNtpFlux::SetDefaults().

  {
    // The flux ntuples can be recycled for a number of times to boost 
    // generated event statistics without requiring enormous beam simulation
    // statistics.
    // That option determines how many times the driver is going to cycle 
    // through the input flux ntuple.
    // With ncycle=0 the flux ntuple will be recycled an infinite amount of 
    // times so that the event generation loop can exit only on a POT or 
    // event num check.

    fNCycles = TMath::Max(0L, ncycle);
  }

void genie::flux::GDk2NuFlux::SetTreeNames ( std::string  fname = "dk2nuTree", std::string  mname = "dkmetaTree"  )

inline

Definition at line 154 of file GDk2NuFlux.h.

References plot_validation_datamc::fname, and string.

  { fTreeNames[0] = fname; fTreeNames[1] = mname; }

void genie::flux::GDk2NuFlux::SetTreeNames ( std::string  fname = "dk2nuTree", std::string  mname = "dkmetaTree"  )

inline

Definition at line 154 of file GDk2NuFlux.h.

References plot_validation_datamc::fname, and string.

  { fTreeNames[0] = fname; fTreeNames[1] = mname; }

void genie::flux::GFluxFileConfigI::SetUpstreamZ ( double  z0 )

virtualinherited

set flux neutrino initial z position (upstream of the detector) pushed back from the normal flux window

Definition at line 74 of file GFluxFileConfigI.cxx.

References genie::flux::GFluxFileConfigI::fZ0.

Referenced by fill_simple(), genie::flux::GFluxFileConfigI::GetXMLFileBase(), main(), SetDefaults(), and genie::flux::GSimpleNtpFlux::SetDefaults().

  {
    // The flux neutrino position (x,y) is given on the user specified 
    // flux window.  This method sets the preferred user coord starting z 
    // position upstream of detector face. Each flux neutrino will be 
    // backtracked from the initial flux window to the input z0.  
    // If the value is unreasonable (> 10^30) then the ray is left on 
    // the flux window.

    fZ0 = z0;
  }

void genie::flux::GFluxFileConfigI::SetXMLFileBase ( std::string  xmlbasename = "" )

virtualinherited

Definition at line 30 of file GFluxFileConfigI.cxx.

References genie::flux::GFluxFileConfigI::fXMLbasename.

Referenced by SetDefaults().

  { fXMLbasename = xmlbasename; }

TVector3 GDk2NuFlux::SphereNormal

(

)

const

Definition at line 978 of file GDk2NuFlux.cxx.

References dir, fCurNuChoice, fFluxSphereCenterBeam, and bsim::NuChoice::x4NuBeam.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  // return the normal for this current location on the sphere
  TVector3 dir = (fCurNuChoice->x4NuBeam.Vect() - fFluxSphereCenterBeam);
  return dir.Unit();
}

TVector3 genie::flux::GDk2NuFlux::SphereNormal

(

)

const

genie::flux::Exposure_t genie::flux::GFluxExposureI::StringToEnum ( const char *  chars, int  maxChar = 0  )

staticinherited

Definition at line 31 of file GFluxExposureI.cxx.

References MECModelEnuComparisons::i, genie::flux::kPOTs, genie::flux::kSeconds, and genie::flux::kUnknown.

  {
    int len = strlen(chars);
    if (maxChar == 0  ) maxChar = len;
    if (maxChar >  len) maxChar = len;
    if (maxChar <= 0  ) return genie::flux::kUnknown;
    
    char* lowchars = new char [maxChar];
    for (int i=0; i<maxChar; ++i) lowchars[i] = tolower(chars[i]);
    
    genie::flux::Exposure_t etype = genie::flux::kUnknown;
    if        ( 0 == strncmp(lowchars,"pot",TMath::Min(maxChar,3))) {
      etype = genie::flux::kPOTs;
    } else if ( 0 == strncmp(lowchars,"sec",TMath::Min(maxChar,3))) {
      etype = genie::flux::kSeconds;
    }
    delete [] lowchars;
    return etype;
  }

double genie::flux::GDk2NuFlux::SumWeight ( void  ) const

inline

integrated weight for flux neutrinos looped so far

Definition at line 127 of file GDk2NuFlux.h.

References site_stats_from_log::filenames, and string.

double genie::flux::GDk2NuFlux::SumWeight ( void  ) const

inline

integrated weight for flux neutrinos looped so far

Definition at line 127 of file GDk2NuFlux.h.

References site_stats_from_log::filenames, and string.

double genie::flux::GDk2NuFlux::UsedPOTs

(

void

)

const

of protons-on-target used

double GDk2NuFlux::UsedPOTs

(

void

)

const

of protons-on-target used

Definition at line 552 of file GDk2NuFlux.cxx.

References fAccumPOTs, fNuFluxTree, LOG, and pWARN.

Referenced by GetTotalExposure(), POT_curr(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
// Compute current number of flux POTs

  if (!fNuFluxTree) {
     LOG("Flux", pWARN)
          << "The flux driver has not been properly configured";
     return 0;
  }
  return fAccumPOTs;
}

void genie::flux::GDk2NuFlux::User2BeamDir	(	const TLorentzVector &	usrdir,
		TLorentzVector &	beamdir
	)		const

void GDk2NuFlux::User2BeamDir	(	const TLorentzVector &	usrdir,
		TLorentzVector &	beamdir
	)		const

Definition at line 1030 of file GDk2NuFlux.cxx.

References fBeamRotInv, and fLengthScaleU2B.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  beamdir = fLengthScaleU2B*(fBeamRotInv*usrdir);
}

void genie::flux::GDk2NuFlux::User2BeamP4	(	const TLorentzVector &	usrp4,
		TLorentzVector &	beamp4
	)		const

void GDk2NuFlux::User2BeamP4	(	const TLorentzVector &	usrp4,
		TLorentzVector &	beamp4
	)		const

Definition at line 1035 of file GDk2NuFlux.cxx.

References fBeamRotInv.

Referenced by genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  beamp4 = fBeamRotInv*usrp4;
}

void GDk2NuFlux::User2BeamPos	(	const TLorentzVector &	usrxyz,
		TLorentzVector &	beamxyz
	)		const

Definition at line 1025 of file GDk2NuFlux.cxx.

References fBeamRotInv, fBeamZero, and fLengthScaleU2B.

Referenced by PrintConfig(), SetFluxSphere(), SetFluxWindow(), and genie::flux::GDk2NuFluxXMLHelper::~GDk2NuFluxXMLHelper().

{
  beamxyz = fLengthScaleU2B*(fBeamRotInv*(usrxyz-fBeamZero));
}

void genie::flux::GDk2NuFlux::User2BeamPos	(	const TLorentzVector &	usrxyz,
		TLorentzVector &	beamxyz
	)		const

double genie::flux::GDk2NuFlux::Weight ( void  )

inlinevirtual

returns the flux neutrino weight (if any)

Implements genie::GFluxI.

Definition at line 86 of file GDk2NuFlux.h.

{ return  fWeight;              }

double genie::flux::GDk2NuFlux::Weight ( void  )

inlinevirtual

returns the flux neutrino weight (if any)

Implements genie::GFluxI.

Definition at line 86 of file GDk2NuFlux.h.

Referenced by GenerateNext(), GenerateNext_weighted(), and ScanForMaxWeight().

{ return  fWeight;              }

Member Data Documentation

double genie::flux::GDk2NuFlux::fAccumPOTs

private

POTs used so far.

Definition at line 304 of file GDk2NuFlux.h.

Referenced by Clear(), GenerateNext_weighted(), Initialize(), LoadBeamSimData(), PrintConfig(), and UsedPOTs().

bool genie::flux::GDk2NuFlux::fApplyTiltWeight

private

wgt due to window normal not || beam

Definition at line 307 of file GDk2NuFlux.h.

Referenced by GenerateNext_weighted(), Initialize(), and PrintConfig().

TLorentzRotation genie::flux::GDk2NuFlux::fBeamRot

private

rotation applied beam –> user coord

Definition at line 315 of file GDk2NuFlux.h.

Referenced by Beam2UserDir(), Beam2UserP4(), Beam2UserPos(), GetBeamRotation(), PrintConfig(), and SetBeamRotation().

TLorentzRotation genie::flux::GDk2NuFlux::fBeamRotInv

private

Definition at line 316 of file GDk2NuFlux.h.

Referenced by PrintConfig(), SetBeamRotation(), User2BeamDir(), User2BeamP4(), and User2BeamPos().

TLorentzVector genie::flux::GDk2NuFlux::fBeamZero

private

beam origin in user coords

Definition at line 314 of file GDk2NuFlux.h.

Referenced by Beam2UserPos(), GetBeamCenter(), PrintConfig(), SetBeamCenter(), SetLengthUnits(), and User2BeamPos().

bsim::Dk2Nu * genie::flux::GDk2NuFlux::fCurDk2Nu

private

Definition at line 275 of file GDk2NuFlux.h.

Referenced by AddFile(), GenerateNext(), GenerateNext_weighted(), GetBranchInfo(), Initialize(), LoadDkMeta(), PrintCurrent(), and ResetCurrent().

bsim::DkMeta * genie::flux::GDk2NuFlux::fCurDkMeta

private

Definition at line 276 of file GDk2NuFlux.h.

Referenced by AddFile(), Initialize(), and LoadDkMeta().

bsim::NuChoice * genie::flux::GDk2NuFlux::fCurNuChoice

private

Definition at line 277 of file GDk2NuFlux.h.

Referenced by AddFile(), CleanUp(), GenerateNext(), GenerateNext_weighted(), GetBranchInfo(), GetDecayDist(), Initialize(), Momentum(), MoveToZ0(), PdgCode(), Position(), PrintCurrent(), ResetCurrent(), ScanForMaxWeight(), SetLengthUnits(), and SphereNormal().

bool genie::flux::GDk2NuFlux::fDetLocIsSet

private

is a flux location (near/far) set?

Definition at line 308 of file GDk2NuFlux.h.

Referenced by Initialize(), LoadBeamSimData(), PrintConfig(), ScanForMaxWeight(), and SetFluxWindow().

double genie::flux::GDk2NuFlux::fEffPOTsPerNu

private

what a entry is worth ...

Definition at line 303 of file GDk2NuFlux.h.

Referenced by CalcEffPOTsPerNu(), GenerateNext_weighted(), Initialize(), and PrintConfig().

bool genie::flux::GDk2NuFlux::fEnd

private

end condition reached

Definition at line 257 of file GDk2NuFlux.h.

Referenced by GenerateNext_weighted(), and Initialize().

Long64_t genie::flux::GDk2NuFlux::fFilePOTs

private

of protons-on-target represented by all files

Definition at line 283 of file GDk2NuFlux.h.

Referenced by AddFile(), CalcEffPOTsPerNu(), Initialize(), and PrintConfig().

TVector3 genie::flux::GDk2NuFlux::fFluxSphereCenterBeam

private

center for flux sphere - beam coords

Definition at line 327 of file GDk2NuFlux.h.

Referenced by GetFluxSphere(), PrintConfig(), SetFluxSphere(), and SphereNormal().

TVector3 genie::flux::GDk2NuFlux::fFluxSphereCenterUser

private

center for flux sphere - user coords

Definition at line 326 of file GDk2NuFlux.h.

Referenced by GetFluxSphere(), PrintConfig(), and SetFluxSphere().

double genie::flux::GDk2NuFlux::fFluxSphereRadius

private

radius for flux sphere

Definition at line 328 of file GDk2NuFlux.h.

Referenced by GetFluxSphere(), PrintConfig(), and SetFluxSphere().

TLorentzVector genie::flux::GDk2NuFlux::fFluxWindowBase

private

base point for flux window - beam coord

Definition at line 320 of file GDk2NuFlux.h.

Referenced by GenerateNext_weighted(), PrintConfig(), and SetFluxWindow().

TLorentzVector genie::flux::GDk2NuFlux::fFluxWindowDir1

private

extent for flux window (direction 1)

Definition at line 321 of file GDk2NuFlux.h.

Referenced by CalcEffPOTsPerNu(), GenerateNext_weighted(), PrintConfig(), and SetFluxWindow().

TLorentzVector genie::flux::GDk2NuFlux::fFluxWindowDir2

private

extent for flux window (direction 2)

Definition at line 322 of file GDk2NuFlux.h.

Referenced by CalcEffPOTsPerNu(), GenerateNext_weighted(), PrintConfig(), and SetFluxWindow().

double genie::flux::GDk2NuFlux::fFluxWindowLen1

private

Definition at line 323 of file GDk2NuFlux.h.

Referenced by PrintConfig(), and SetFluxWindow().

double genie::flux::GDk2NuFlux::fFluxWindowLen2

private

Definition at line 324 of file GDk2NuFlux.h.

Referenced by PrintConfig(), and SetFluxWindow().

TVector3 genie::flux::GDk2NuFlux::fFluxWindowNormal

private

normal direction for flux window – beam coord

Definition at line 325 of file GDk2NuFlux.h.

Referenced by SetFluxWindow().

TVector3 genie::flux::GDk2NuFlux::fFluxWindowPtUser

private

user points of flux window

Definition at line 319 of file GDk2NuFlux.h.

Referenced by GetFluxWindow(), PrintConfig(), SetFluxWindow(), and SetLengthUnits().

bool genie::flux::GDk2NuFlux::fGenWeighted

private

does GenerateNext() give weights?

Definition at line 306 of file GDk2NuFlux.h.

Referenced by GenerateNext(), GenerateWeighted(), Initialize(), and PrintConfig().

TLorentzVector genie::flux::GDk2NuFlux::fgX4dkvtx

private

decay 4-position beam coord

Definition at line 330 of file GDk2NuFlux.h.

Referenced by GenerateNext_weighted(), and GetDecayDist().

long int genie::flux::GFluxFileConfigI::fICycle

protectedinherited

current file cycle default 0 = infinitely

Definition at line 98 of file GFluxFileConfigI.h.

Referenced by CleanUp(), genie::flux::GSimpleNtpFlux::CleanUp(), Clear(), genie::flux::GSimpleNtpFlux::Clear(), GenerateNext_weighted(), genie::flux::GSimpleNtpFlux::GenerateNext_weighted(), Initialize(), genie::flux::GSimpleNtpFlux::Initialize(), LoadBeamSimData(), genie::flux::GSimpleNtpFlux::LoadBeamSimData(), and PrintConfig().

Long64_t genie::flux::GDk2NuFlux::fIEntry

private

current flux ntuple entry

Definition at line 281 of file GDk2NuFlux.h.

Referenced by CleanUp(), GenerateNext_weighted(), Initialize(), LoadBeamSimData(), and PrintConfig().

bool genie::flux::GDk2NuFlux::fIsSphere

private

doing this on a sphere rather than a flat window?

Definition at line 318 of file GDk2NuFlux.h.

Referenced by Initialize(), SetFluxSphere(), and SetFluxWindow().

long int genie::flux::GDk2NuFlux::fIUse

private

current # of times an entry has been used

Definition at line 300 of file GDk2NuFlux.h.

Referenced by CleanUp(), GenerateNext_weighted(), Initialize(), LoadBeamSimData(), and PrintConfig().

std::map< int, int > genie::flux::GDk2NuFlux::fJobToMetaIndex

private

quick lookup from job# to meta chain

Definition at line 285 of file GDk2NuFlux.h.

Referenced by LoadDkMeta().

double genie::flux::GDk2NuFlux::fLengthScaleB2U

private

scale factor beam (cm) –> user

Definition at line 311 of file GDk2NuFlux.h.

Referenced by Beam2UserDir(), Beam2UserPos(), GetDecayDist(), LengthUnits(), PrintConfig(), and SetLengthUnits().

double genie::flux::GDk2NuFlux::fLengthScaleU2B

private

scale factor beam user –> (cm)

Definition at line 312 of file GDk2NuFlux.h.

Referenced by MoveToZ0(), PrintConfig(), SetLengthUnits(), User2BeamDir(), and User2BeamPos().

double genie::flux::GDk2NuFlux::fLengthUnits

private

units for coord in user exchanges

Definition at line 310 of file GDk2NuFlux.h.

Referenced by PrintConfig(), and SetLengthUnits().

double genie::flux::GDk2NuFlux::fMaxEFudge

private

fudge factor for estmating max enu (0=> use fixed 120GeV)

Definition at line 295 of file GDk2NuFlux.h.

Referenced by Initialize(), and ScanForMaxWeight().

double genie::flux::GDk2NuFlux::fMaxEv

private

maximum energy

Definition at line 255 of file GDk2NuFlux.h.

Referenced by GenerateNext_weighted(), Initialize(), PrintConfig(), ScanForMaxWeight(), and SetMaxEnergy().

double genie::flux::GDk2NuFlux::fMaxWeight

private

max flux neutrino weight in input file

Definition at line 288 of file GDk2NuFlux.h.

Referenced by GenerateNext(), GenerateNext_weighted(), Initialize(), LoadBeamSimData(), PrintConfig(), and ScanForMaxWeight().

double genie::flux::GDk2NuFlux::fMaxWeightInit

private

max of scan & minmaxweight

Definition at line 291 of file GDk2NuFlux.h.

Referenced by PrintConfig(), and ScanForMaxWeight().

double genie::flux::GDk2NuFlux::fMaxWeightMax

private

if "frozen" this is what bump would given

Definition at line 292 of file GDk2NuFlux.h.

Referenced by GenerateNext(), PrintConfig(), and ScanForMaxWeight().

double genie::flux::GDk2NuFlux::fMaxWeightScan

private

initial estimate from scan

Definition at line 290 of file GDk2NuFlux.h.

Referenced by PrintConfig(), and ScanForMaxWeight().

long int genie::flux::GDk2NuFlux::fMaxWgtEntries

private

of entries in estimating max wgt

Definition at line 294 of file GDk2NuFlux.h.

Referenced by Initialize(), PrintConfig(), and ScanForMaxWeight().

long int genie::flux::GDk2NuFlux::fMaxWgtExceeded

private

track failures of estimate

Definition at line 296 of file GDk2NuFlux.h.

Referenced by GenerateNext(), Initialize(), PrintConfig(), and ScanForMaxWeight().

int genie::flux::GDk2NuFlux::fMaxWgtFailModel

private

what to do ... 0=bump, 1=frozen, 2=abort

Definition at line 297 of file GDk2NuFlux.h.

Referenced by GenerateNext(), Initialize(), and PrintConfig().

double genie::flux::GDk2NuFlux::fMaxWgtFudge

private

fudge factor for estimating max wgt

Definition at line 293 of file GDk2NuFlux.h.

Referenced by GenerateNext(), Initialize(), PrintConfig(), and ScanForMaxWeight().

double genie::flux::GDk2NuFlux::fMinMaxWeight

private

user set lower limit on estimate

Definition at line 289 of file GDk2NuFlux.h.

Referenced by Initialize(), PrintConfig(), and ScanForMaxWeight().

long int genie::flux::GFluxFileConfigI::fNCycles

protectedinherited

times to cycle through the ntuple(s)

Definition at line 97 of file GFluxFileConfigI.h.

Referenced by CleanUp(), genie::flux::GSimpleNtpFlux::CleanUp(), GenerateNext_weighted(), genie::flux::GSimpleNtpFlux::GenerateNext_weighted(), PrintConfig(), and genie::flux::GFluxFileConfigI::SetNumOfCycles().

Long64_t genie::flux::GDk2NuFlux::fNEntries

private

number of flux ntuple entries

Definition at line 280 of file GDk2NuFlux.h.

Referenced by CalcEffPOTsPerNu(), GenerateNext_weighted(), Initialize(), LoadBeamSimData(), and PrintConfig().

int genie::flux::GDk2NuFlux::fNFiles

private

number of files in chain

Definition at line 279 of file GDk2NuFlux.h.

Referenced by AddFile(), Initialize(), and PrintConfig().

long int genie::flux::GDk2NuFlux::fNNeutrinos

private

number of flux neutrinos thrown so far

Definition at line 302 of file GDk2NuFlux.h.

Referenced by Clear(), GenerateNext_weighted(), Initialize(), LoadBeamSimData(), and PrintConfig().

std::vector< std::string > genie::flux::GDk2NuFlux::fNuFluxFilePatterns

private

(potentially wildcarded) path(s)

Definition at line 269 of file GDk2NuFlux.h.

Referenced by LoadBeamSimData(), and PrintConfig().

TChain * genie::flux::GDk2NuFlux::fNuFluxTree

private

TTree // REF ONLY!

Definition at line 272 of file GDk2NuFlux.h.

Referenced by AddFile(), CalcEffPOTsPerNu(), GenerateNext_weighted(), GetFileList(), Initialize(), LoadBeamSimData(), and UsedPOTs().

TChain * genie::flux::GDk2NuFlux::fNuMetaTree

private

TTree // REF ONLY!

Definition at line 273 of file GDk2NuFlux.h.

Referenced by AddFile(), GetMetaDataTree(), Initialize(), LoadBeamSimData(), and LoadDkMeta().

long int genie::flux::GDk2NuFlux::fNUse

private

how often to use same entry in a row

Definition at line 299 of file GDk2NuFlux.h.

Referenced by CleanUp(), Initialize(), PrintConfig(), and SetEntryReuse().

Long64_t genie::flux::GDk2NuFlux::fNuTot

private

cummulative # of entries (=fNEntries)

Definition at line 282 of file GDk2NuFlux.h.

Referenced by AddFile(), and Initialize().

PDGCodeList* genie::flux::GFluxFileConfigI::fPdgCList

protectedinherited

list of neutrino pdg-codes to generate

Definition at line 94 of file GFluxFileConfigI.h.

Referenced by CleanUp(), genie::flux::GSimpleNtpFlux::CleanUp(), GenerateNext_weighted(), genie::flux::GSimpleNtpFlux::GenerateNext_weighted(), PrintConfig(), genie::flux::GSimpleNtpFlux::ProcessMeta(), and genie::flux::GFluxFileConfigI::SetFluxParticles().

PDGCodeList* genie::flux::GFluxFileConfigI::fPdgCListRej

protectedinherited

list of nu pdg-codes seen but rejected

Definition at line 95 of file GFluxFileConfigI.h.

Referenced by CleanUp(), genie::flux::GSimpleNtpFlux::CleanUp(), GenerateNext_weighted(), genie::flux::GSimpleNtpFlux::GenerateNext_weighted(), and PrintConfig().

double genie::flux::GDk2NuFlux::fSumWeight

private

sum of weights for nus thrown so far

Definition at line 301 of file GDk2NuFlux.h.

Referenced by Clear(), GenerateNext_weighted(), Initialize(), LoadBeamSimData(), and PrintConfig().

std::string genie::flux::GDk2NuFlux::fTreeNames

private

pair of names "dk2nuTree", "dkmetaTree"

Definition at line 271 of file GDk2NuFlux.h.

Referenced by AddFile(), Initialize(), and LoadBeamSimData().

double genie::flux::GDk2NuFlux::fWeight

private

current neutrino weight, =1 if generating unweighted entries

Definition at line 287 of file GDk2NuFlux.h.

Referenced by GenerateNext(), and GenerateNext_weighted().

std::string genie::flux::GFluxFileConfigI::fXMLbasename

protectedinherited

XML file that might hold config param_sets.

Definition at line 96 of file GFluxFileConfigI.h.

Referenced by genie::flux::GFluxFileConfigI::GetXMLFileBase(), and genie::flux::GFluxFileConfigI::SetXMLFileBase().

double genie::flux::GFluxFileConfigI::fZ0

protectedinherited

configurable starting z position for each flux neutrino (in detector coord system)

Definition at line 100 of file GFluxFileConfigI.h.

Referenced by GenerateNext_weighted(), genie::flux::GSimpleNtpFlux::GenerateNext_weighted(), PrintConfig(), and genie::flux::GFluxFileConfigI::SetUpstreamZ().

---

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

• dk2nu/genie/GDk2NuFlux.h
• dk2nu/genie/GDk2NuFlux.cxx