genie::flux::GPowerSpectrumAtmoFlux Class Reference

A driver for a power spectrum atmospheric neutrino flux. Elements extensively reused from GAtmoFlux. More...

#include "GPowerSpectrumAtmoFlux.h"

Inheritance diagram for genie::flux::GPowerSpectrumAtmoFlux:

Public Member Functions
	GPowerSpectrumAtmoFlux ()
	~GPowerSpectrumAtmoFlux ()
const PDGCodeList &	FluxParticles (void) override
	declare list of flux neutrinos that can be generated (for init. purposes) More...
double	MaxEnergy (void) override
	declare the max flux neutrino energy that can be generated (for init. purposes) More...
bool	GenerateNext (void) override
	generate the next flux neutrino (return false in err) More...
int	PdgCode (void) override
	returns the flux neutrino pdg code More...
double	Weight (void) override
	returns the flux neutrino weight (if any) More...
const TLorentzVector &	Momentum (void) override
	returns the flux neutrino 4-momentum More...
const TLorentzVector &	Position (void) override
	returns the flux neutrino 4-position (note: expect SI rather than physical units) More...
bool	End (void) override
	true if no more flux nu's can be thrown (eg reaching end of beam sim ntuples) More...
long int	Index (void) override
	returns corresponding index for current flux neutrino (e.g. for a flux ntuple returns the current entry number) More...
void	Clear (Option_t *opt) override
	reset state variables based on opt More...
void	GenerateWeighted (bool gen_weighted) override
	set whether to generate weighted or unweighted neutrinos More...
double	MinEnergy (void)
void	SetSpectralIndex (double index)
void	SetUserCoordSystem (TRotation &rotation)
void	SetRadii (double Rlongitudinal, double Rtransverse)
void	SetFlavors (std::vector< int > flavors)
void	SetMinEnergy (double Emin)
void	SetMaxEnergy (double Emax)
long int	NFluxNeutrinos (void) const
void	ResetNFluxNeutrinos (void)
void	AddFluxFile (int neutrino_pdg, string filename)
bool	LoadFluxData (void)
double	GetFlux (int flavour, double energy, double costh, double phi)
double	ComputeWeight (int flavour, double energy, double costh, double phi)
void	InitializeWeight ()

Private Member Functions
bool	FillFluxHisto (int nu_pdg, string filename)
void	AddAllFluxes (void)
TH3D *	CreateNormalisedFluxHisto (TH3D *hist)
void	ResetSelection (void)
void	Initialize (void)
void	CleanUp (void)

Private Attributes
PDGCodeList *	fPdgCList
	input list of neutrino pdg-codes More...
double	fMaxEvCut
	user-defined cut: maximum energy More...
double	fMinEvCut
	user-defined cut: minimum energy More...
int	fgPdgC
	current generated nu pdg-code More...
TLorentzVector	fgP4
	current generated nu 4-momentum More...
TLorentzVector	fgX4
	current generated nu 4-position More...
double	fWeight
	current generated nu weight More...
double	fSpectralIndex
	power law function More...
double	fRl
	defining flux neutrino generation surface: longitudinal radius More...
double	fRt
	defining flux neutrino generation surface: transverse radius More...
double	fGlobalGenWeight
	global generation weight to apply to all events More...
double	fAgen
	current generation area More...
TRotation	fRotTHz2User
	coord. system rotation: THZ -> Topocentric user-defined More...
long int	fNNeutrinos
	number of flux neutrinos thrown so far More...
vector< int >	fFluxFlavour
	input flux file for each neutrino species More...
vector< string >	fFluxFile
	input flux file for each neutrino species More...
map< int, TH3D * >	fRawFluxHistoMap
	flux = f(Ev,cos8,phi) for each neutrino species More...
map< int, TH2D * >	fRawFluxHistoMap2D
	flux = f(Ev,cos8) for each neutrino species More...

Detailed Description

A driver for a power spectrum atmospheric neutrino flux. Elements extensively reused from GAtmoFlux.

Author

Pierre Granger grang.nosp@m.er@a.nosp@m.pc.in.nosp@m.2p3..nosp@m.fr APC (CNRS)

December 16, 2022

Definition at line 28 of file GPowerSpectrumAtmoFlux.h.

Constructor & Destructor Documentation

GPowerSpectrumAtmoFlux::GPowerSpectrumAtmoFlux

(

)

Definition at line 27 of file GPowerSpectrumAtmoFlux.cxx.

References Initialize().

{
        LOG("Flux", pNOTICE)
        << "Instantiating the GENIE Power Spectrum atmospheric neutrino flux driver";
        this->Initialize();
}

GPowerSpectrumAtmoFlux::~GPowerSpectrumAtmoFlux

(

)

Definition at line 36 of file GPowerSpectrumAtmoFlux.cxx.

{

}

Member Function Documentation

void genie::flux::GPowerSpectrumAtmoFlux::AddAllFluxes ( void  )

private

void GPowerSpectrumAtmoFlux::AddFluxFile	(	int	neutrino_pdg,
		string	filename
	)

Definition at line 402 of file GPowerSpectrumAtmoFlux.cxx.

References f, fFluxFile, and fFluxFlavour.

Referenced by evgb::GENIEHelper::InitializeFluxDriver().

{
  // Check file exists
  std::ifstream f(filename.c_str());
  if (!f.good()) {
    LOG("Flux", pFATAL) << "Flux file does not exist "<<filename;
    exit(-1);
  }
  if ( pdg::IsNeutrino(nu_pdg) || pdg::IsAntiNeutrino(nu_pdg) ) {
    fFluxFlavour.push_back(nu_pdg);
    fFluxFile.push_back(filename);
  } else {
    LOG ("Flux", pWARN)
      << "Input particle code: " << nu_pdg << " not a neutrino!";
  }
}

void GPowerSpectrumAtmoFlux::CleanUp ( void  )

private

Definition at line 339 of file GPowerSpectrumAtmoFlux.cxx.

References fPdgCList.

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

        if (fPdgCList) delete fPdgCList;
}

void GPowerSpectrumAtmoFlux::Clear ( Option_t *  opt )

override

reset state variables based on opt

Definition at line 245 of file GPowerSpectrumAtmoFlux.cxx.

References ResetNFluxNeutrinos().

{
        LOG("Flux", pWARN) << "GPowerSpectrumAtmoFlux::Clear(" << opt << ") called";
        this->ResetNFluxNeutrinos();
}

double GPowerSpectrumAtmoFlux::ComputeWeight	(	int	flavour,
		double	energy,
		double	costh,
		double	phi
	)

Definition at line 492 of file GPowerSpectrumAtmoFlux.cxx.

References fAgen, fGlobalGenWeight, fSpectralIndex, and GetFlux().

Referenced by GenerateNext().

{
  double flux = this->GetFlux(flavour, energy, costh, phi);

  return flux*fAgen*fGlobalGenWeight*pow(energy, fSpectralIndex);
}

TH3D* genie::flux::GPowerSpectrumAtmoFlux::CreateNormalisedFluxHisto ( TH3D *  hist )

private

bool GPowerSpectrumAtmoFlux::End ( void  )

override

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

Definition at line 224 of file GPowerSpectrumAtmoFlux.cxx.

{
        return false;
}

bool GPowerSpectrumAtmoFlux::FillFluxHisto ( int  nu_pdg, string  filename  )

private

Definition at line 372 of file GPowerSpectrumAtmoFlux.cxx.

References f, fRawFluxHistoMap, fRawFluxHistoMap2D, and h2.

Referenced by LoadFluxData().

                                                                     {
        LOG("Flux", pNOTICE)
    << "Loading fine grained flux for neutrino: " << nu_pdg
    << " from file: " << filename;

        TH3D* histo = nullptr;

        TFile *f = new TFile(filename.c_str(), "READ");
        f->GetObject("flux", histo);

        if(!histo) {
                LOG("Flux", pERROR) << "Null flux histogram!";
                f->Close();
                return false;
        }

        TH3D* h = static_cast<TH3D*>(histo->Clone());
        f->Close();

        fRawFluxHistoMap.insert(std::make_pair(nu_pdg, h));

        TH2D* h2 = static_cast<TH2D*>(h->Project3D("yx"));

        fRawFluxHistoMap2D.insert(std::make_pair(nu_pdg, h2));

        return true;
}

const PDGCodeList & GPowerSpectrumAtmoFlux::FluxParticles ( void  )

override

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

Definition at line 44 of file GPowerSpectrumAtmoFlux.cxx.

References fPdgCList.

{
        return *fPdgCList;
}

bool GPowerSpectrumAtmoFlux::GenerateNext ( void  )

override

generate the next flux neutrino (return false in err)

Definition at line 65 of file GPowerSpectrumAtmoFlux.cxx.

References ComputeWeight(), fgP4, fgPdgC, fgX4, fNNeutrinos, fPdgCList, fRl, fRotTHz2User, fRt, fSpectralIndex, fWeight, MaxEnergy(), MinEnergy(), ResetSelection(), weight, x, y, and z.

{
        // Reset previously generated neutrino code / 4-p / 4-x
        this->ResetSelection();

        // Get a RandomGen instance
        RandomGen * rnd = RandomGen::Instance();

        // Generate (Ev, costheta, phi)
        double Ev       = 0.;
        double costheta = 0.;
        double phi      = 0;
        double weight   = 0;
        int    nu_pdg   = 0;

        // generate events according to a power law spectrum,
        // then weight events by inverse power law
        // (note: cannot use index alpha=1)
        double alpha = fSpectralIndex;

        double emin = TMath::Power(this->MinEnergy(),1.0-alpha);
        double emax = TMath::Power(this->MaxEnergy(),1.0-alpha);
        Ev          = TMath::Power(emin+(emax-emin)*rnd->RndFlux().Rndm(),1.0/(1.0-alpha));
        costheta    = -1+2*rnd->RndFlux().Rndm();
        phi         = 2.*kPi* rnd->RndFlux().Rndm();

        unsigned int nnu = fPdgCList->size();
        unsigned int inu = rnd->RndFlux().Integer(nnu);
        nu_pdg   = (*fPdgCList)[inu];
        weight = this->ComputeWeight(nu_pdg, Ev, costheta, phi);

        // Compute etc trigonometric numbers
        double sintheta  = TMath::Sqrt(1-costheta*costheta);
        double cosphi    = TMath::Cos(phi);
        double sinphi    = TMath::Sin(phi);

        // Set the neutrino pdg code
        fgPdgC = nu_pdg;

        // Set the neutrino weight
        fWeight = weight;

        // Compute the neutrino momentum
        // The `-1' means it is directed towards the detector.
        double pz = -1.* Ev * costheta;
        double py = -1.* Ev * sintheta * sinphi;
        double px = -1.* Ev * sintheta * cosphi;

        // Default vertex is at the origin
        double z = 0.0;
        double y = 0.0;
        double x = 0.0;

        // Shift the neutrino position onto the flux generation surface.
        // The position is computed at the surface of a sphere with R=fRl
        // at the topocentric horizontal (THZ) coordinate system.
        if( fRl>0.0 ){
        z += fRl * costheta;
        y += fRl * sintheta * sinphi;
        x += fRl * sintheta * cosphi;
        }

        // Apply user-defined rotation from THZ -> user-defined topocentric
        // coordinate system.
        if( !fRotTHz2User.IsIdentity() )
        {
        TVector3 tx3(x, y, z );
        TVector3 tp3(px,py,pz);

        tx3 = fRotTHz2User * tx3;
        tp3 = fRotTHz2User * tp3;

        x  = tx3.X();
        y  = tx3.Y();
        z  = tx3.Z();
        px = tp3.X();
        py = tp3.Y();
        pz = tp3.Z();
        }

        // If the position is left as is, then all generated neutrinos
        // would point towards the origin.
        // Displace the position randomly on the surface that is
        // perpendicular to the selected point P(xo,yo,zo) on the sphere
        if( fRt>0.0 ){
        TVector3 vec(x,y,z);               // vector towards selected point
        TVector3 dvec1 = vec.Orthogonal(); // orthogonal vector
        TVector3 dvec2 = dvec1;            // second orthogonal vector
        dvec2.Rotate(-kPi/2.0,vec);        // rotate second vector by 90deg,
                                           // now forming a new orthogonal cartesian coordinate system
        double psi = 2.*kPi* rnd->RndFlux().Rndm(); // rndm angle [0,2pi]
        double random = rnd->RndFlux().Rndm();      // rndm number  [0,1]
        dvec1.SetMag(TMath::Sqrt(random)*fRt*TMath::Cos(psi));
        dvec2.SetMag(TMath::Sqrt(random)*fRt*TMath::Sin(psi));
        x += dvec1.X() + dvec2.X();
        y += dvec1.Y() + dvec2.Y();
        z += dvec1.Z() + dvec2.Z();
        }

        // Set the neutrino momentum and position 4-vectors with values
        // calculated at previous steps.
        fgP4.SetPxPyPzE(px, py, pz, Ev);
        fgX4.SetXYZT   (x,  y,  z,  0.);

        // Increment flux neutrino counter used for sample normalization purposes.
        fNNeutrinos++;

        // Report and exit
        LOG("Flux", pINFO)
           << "Generated neutrino: "
           << "\n pdg-code: " << fgPdgC
           << "\n p4: " << utils::print::P4AsShortString(&fgP4)
           << "\n x4: " << utils::print::X4AsString(&fgX4);

        return true;
}

void GPowerSpectrumAtmoFlux::GenerateWeighted ( bool  gen_weighted )

override

set whether to generate weighted or unweighted neutrinos

Definition at line 253 of file GPowerSpectrumAtmoFlux.cxx.

{
// Dummy clear method needed to conform to GFluxI interface
//
  LOG("Flux", pERROR) <<
  "The neutrinos are always generated weighted with this implementation!! GenerateWeighted method not implemented.";
}

double GPowerSpectrumAtmoFlux::GetFlux	(	int	flavour,
		double	energy,
		double	costh,
		double	phi
	)

Definition at line 468 of file GPowerSpectrumAtmoFlux.cxx.

References fRawFluxHistoMap, fRawFluxHistoMap2D, and h2.

Referenced by ComputeWeight().

{
  TH3D* flux_hist = nullptr;
  std::map<int,TH3D*>::iterator it = fRawFluxHistoMap.find(flavour);
  if(it != fRawFluxHistoMap.end())
  {
    flux_hist = it->second;
  }

  LOG("Flux", pERROR) << "flux_hist: " << flux_hist;

  if(!flux_hist) return 0.0;

  if(flux_hist->GetZaxis()->GetNbins() == 1){ //no binning in phi, bilinear interpolation only so using the 2D hist
        TH2D* h2 = fRawFluxHistoMap2D[it->first];
        return h2->Interpolate(energy, costh);
  }
  else {
        return flux_hist->Interpolate(energy, costh, phi);
  }
}

long int GPowerSpectrumAtmoFlux::Index ( void  )

override

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

Definition at line 231 of file GPowerSpectrumAtmoFlux.cxx.

{
        return -1;
}

void GPowerSpectrumAtmoFlux::Initialize ( void  )

private

Definition at line 284 of file GPowerSpectrumAtmoFlux.cxx.

References fAgen, fMaxEvCut, fMinEvCut, fNNeutrinos, fPdgCList, fRl, fRotTHz2User, fRt, fSpectralIndex, InitializeWeight(), and ResetSelection().

{
        LOG("Flux", pNOTICE) << "Initializing power spectrum atmospheric flux driver";

        bool allow_dup = false;
        fPdgCList = new PDGCodeList(allow_dup);

        fSpectralIndex = 2.0;

        fMinEvCut = 0.01;
        fMaxEvCut = 9999999999;

        // Default radii
        fRl = 0.0;
        fRt = 0.0;
        fAgen = 0.0;

        // Default detector coord system: Topocentric Horizontal Coordinate system
        fRotTHz2User.SetToIdentity();

        // Reset `current' selected flux neutrino
        this->ResetSelection();

        // Reset number of neutrinos thrown so far
        fNNeutrinos = 0;

        // Init the global gen weight
        this->InitializeWeight();
}

void GPowerSpectrumAtmoFlux::InitializeWeight

(

)

Definition at line 501 of file GPowerSpectrumAtmoFlux.cxx.

References fGlobalGenWeight, fMaxEvCut, fMinEvCut, and fSpectralIndex.

Referenced by Initialize(), SetMaxEnergy(), SetMinEnergy(), and SetSpectralIndex().

{
        double IE = 1.;

        if(fSpectralIndex!=1){
                IE = (pow(fMaxEvCut,(1.-fSpectralIndex))-pow(fMinEvCut,(1.-fSpectralIndex)))/(1.-fSpectralIndex);
        }

        double ITheta = 4*kPi;

        fGlobalGenWeight=IE*ITheta;
}

bool GPowerSpectrumAtmoFlux::LoadFluxData

(

void

)

Definition at line 421 of file GPowerSpectrumAtmoFlux.cxx.

References fFluxFile, fFluxFlavour, FillFluxHisto(), fPdgCList, fRawFluxHistoMap, and n.

Referenced by evgb::GENIEHelper::InitializeFluxDriver().

{
  LOG("Flux", pNOTICE)
        << "Loading atmospheric neutrino flux simulation data";

  fPdgCList->clear();

  bool loading_status = true;

  for( unsigned int n=0; n<fFluxFlavour.size(); n++ ){
    int nu_pdg      = fFluxFlavour.at(n);
    string filename = fFluxFile.at(n);
    string pname = PDGLibrary::Instance()->Find(nu_pdg)->GetName();

    LOG("Flux", pNOTICE) << "Loading data for: " << pname;

    bool loaded = this->FillFluxHisto(nu_pdg, filename);

    loading_status = loading_status && loaded;

    if (!loaded) {
        LOG("Flux", pERROR)
          << "Error loading atmospheric neutrino flux simulation data from " << filename;
        break;
    }
  }

  if(loading_status) {
    map<int,TH3D*>::iterator hist_iter = fRawFluxHistoMap.begin();
    for ( ; hist_iter != fRawFluxHistoMap.end(); ++hist_iter) {
      int   nu_pdg = hist_iter->first;
      fPdgCList->push_back(nu_pdg);
    }

    LOG("Flux", pNOTICE)
          << "Atmospheric neutrino flux simulation data loaded!";
    return true;
  }

  LOG("Flux", pERROR)
    << "Error loading atmospheric neutrino flux simulation data";
  return false;
}

double GPowerSpectrumAtmoFlux::MaxEnergy ( void  )

override

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

Definition at line 51 of file GPowerSpectrumAtmoFlux.cxx.

References fMaxEvCut.

Referenced by GenerateNext().

{
  return fMaxEvCut;
}

double GPowerSpectrumAtmoFlux::MinEnergy

(

void

)

Definition at line 58 of file GPowerSpectrumAtmoFlux.cxx.

References fMinEvCut.

Referenced by GenerateNext().

{
  return fMinEvCut;
}

const TLorentzVector & GPowerSpectrumAtmoFlux::Momentum ( void  )

override

returns the flux neutrino 4-momentum

Definition at line 210 of file GPowerSpectrumAtmoFlux.cxx.

References fgP4.

{
        return fgP4;
}

long int GPowerSpectrumAtmoFlux::NFluxNeutrinos

(

void

)

const

Definition at line 365 of file GPowerSpectrumAtmoFlux.cxx.

References fNNeutrinos.

{
  return fNNeutrinos;
}

int GPowerSpectrumAtmoFlux::PdgCode ( void  )

override

returns the flux neutrino pdg code

Definition at line 196 of file GPowerSpectrumAtmoFlux.cxx.

References fgPdgC.

{
        return fgPdgC;
}

const TLorentzVector & GPowerSpectrumAtmoFlux::Position ( void  )

override

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

Definition at line 217 of file GPowerSpectrumAtmoFlux.cxx.

References fgX4.

{
        return fgX4;
}

void GPowerSpectrumAtmoFlux::ResetNFluxNeutrinos

(

void

)

Definition at line 238 of file GPowerSpectrumAtmoFlux.cxx.

References fNNeutrinos.

Referenced by Clear().

{
        fNNeutrinos = 0;
}

void GPowerSpectrumAtmoFlux::ResetSelection ( void  )

private

Definition at line 184 of file GPowerSpectrumAtmoFlux.cxx.

References fgP4, fgPdgC, fgX4, and fWeight.

Referenced by GenerateNext(), and Initialize().

{
// initializing running neutrino pdg-code, 4-position, 4-momentum

  fgPdgC = 0;
  fgP4.SetPxPyPzE (0.,0.,0.,0.);
  fgX4.SetXYZT    (0.,0.,0.,0.);
  fWeight = 0;
}

void GPowerSpectrumAtmoFlux::SetFlavors

(

std::vector< int >

flavors

)

Definition at line 329 of file GPowerSpectrumAtmoFlux.cxx.

References fPdgCList.

Referenced by evgb::GENIEHelper::InitializeFluxDriver().

{
        fPdgCList->clear();
        for(int flavor : flavors){
                fPdgCList->push_back(flavor);
        }
}

void GPowerSpectrumAtmoFlux::SetMaxEnergy

(

double

Emax

)

Definition at line 358 of file GPowerSpectrumAtmoFlux.cxx.

References fMaxEvCut, and InitializeWeight().

Referenced by evgb::GENIEHelper::InitializeFluxDriver().

                                                    {
        fMaxEvCut = Emax;
        this->InitializeWeight(); //Recompute the global weight
}

void GPowerSpectrumAtmoFlux::SetMinEnergy

(

double

Emin

)

Definition at line 348 of file GPowerSpectrumAtmoFlux.cxx.

References fMinEvCut, and InitializeWeight().

Referenced by evgb::GENIEHelper::InitializeFluxDriver().

{
        if(Emin > 0){
                fMinEvCut = Emin;
                this->InitializeWeight(); //Recompute the global weight
        }
}

void GPowerSpectrumAtmoFlux::SetRadii	(	double	Rlongitudinal,
		double	Rtransverse
	)

Definition at line 316 of file GPowerSpectrumAtmoFlux.cxx.

References fAgen, fRl, and fRt.

Referenced by evgb::GENIEHelper::InitializeFluxDriver().

{
  LOG ("Flux", pNOTICE) << "Setting R[longitudinal] = " << Rlongitudinal;
  LOG ("Flux", pNOTICE) << "Setting R[transverse]   = " << Rtransverse;

  fRl = Rlongitudinal;
  fRt = Rtransverse;

  fAgen = kPi*fRt*fRt;
}

void GPowerSpectrumAtmoFlux::SetSpectralIndex

(

double

index

)

Definition at line 263 of file GPowerSpectrumAtmoFlux.cxx.

References fSpectralIndex, and InitializeWeight().

Referenced by evgb::GENIEHelper::InitializeFluxDriver().

{
        if( index != 1.0 ){
                fSpectralIndex = index;
                this->InitializeWeight(); //Recompute the global weight
        }
        else {
                LOG("Flux", pWARN) << "Warning: cannot use a spectral index of unity";
        }
        LOG("Flux", pNOTICE) << "Using Spectral Index = " << index;
}

void GPowerSpectrumAtmoFlux::SetUserCoordSystem

(

TRotation &

rotation

)

Definition at line 277 of file GPowerSpectrumAtmoFlux.cxx.

References fRotTHz2User.

Referenced by evgb::GENIEHelper::InitializeFluxDriver().

{
  fRotTHz2User = rotation;
}

double GPowerSpectrumAtmoFlux::Weight ( void  )

override

returns the flux neutrino weight (if any)

Definition at line 203 of file GPowerSpectrumAtmoFlux.cxx.

References fWeight.

{
        return fWeight;
}

Member Data Documentation

double genie::flux::GPowerSpectrumAtmoFlux::fAgen

private

current generation area

Definition at line 74 of file GPowerSpectrumAtmoFlux.h.

Referenced by ComputeWeight(), Initialize(), and SetRadii().

vector<string> genie::flux::GPowerSpectrumAtmoFlux::fFluxFile

private

input flux file for each neutrino species

Definition at line 78 of file GPowerSpectrumAtmoFlux.h.

Referenced by AddFluxFile(), and LoadFluxData().

vector<int> genie::flux::GPowerSpectrumAtmoFlux::fFluxFlavour

private

input flux file for each neutrino species

Definition at line 77 of file GPowerSpectrumAtmoFlux.h.

Referenced by AddFluxFile(), and LoadFluxData().

double genie::flux::GPowerSpectrumAtmoFlux::fGlobalGenWeight

private

global generation weight to apply to all events

Definition at line 73 of file GPowerSpectrumAtmoFlux.h.

Referenced by ComputeWeight(), and InitializeWeight().

TLorentzVector genie::flux::GPowerSpectrumAtmoFlux::fgP4

private

current generated nu 4-momentum

Definition at line 67 of file GPowerSpectrumAtmoFlux.h.

Referenced by GenerateNext(), Momentum(), and ResetSelection().

int genie::flux::GPowerSpectrumAtmoFlux::fgPdgC

private

current generated nu pdg-code

Definition at line 66 of file GPowerSpectrumAtmoFlux.h.

Referenced by GenerateNext(), PdgCode(), and ResetSelection().

TLorentzVector genie::flux::GPowerSpectrumAtmoFlux::fgX4

private

current generated nu 4-position

Definition at line 68 of file GPowerSpectrumAtmoFlux.h.

Referenced by GenerateNext(), Position(), and ResetSelection().

double genie::flux::GPowerSpectrumAtmoFlux::fMaxEvCut

private

user-defined cut: maximum energy

Definition at line 64 of file GPowerSpectrumAtmoFlux.h.

Referenced by Initialize(), InitializeWeight(), MaxEnergy(), and SetMaxEnergy().

double genie::flux::GPowerSpectrumAtmoFlux::fMinEvCut

private

user-defined cut: minimum energy

Definition at line 65 of file GPowerSpectrumAtmoFlux.h.

Referenced by Initialize(), InitializeWeight(), MinEnergy(), and SetMinEnergy().

long int genie::flux::GPowerSpectrumAtmoFlux::fNNeutrinos

private

number of flux neutrinos thrown so far

Definition at line 76 of file GPowerSpectrumAtmoFlux.h.

Referenced by GenerateNext(), Initialize(), NFluxNeutrinos(), and ResetNFluxNeutrinos().

PDGCodeList* genie::flux::GPowerSpectrumAtmoFlux::fPdgCList

private

input list of neutrino pdg-codes

Definition at line 63 of file GPowerSpectrumAtmoFlux.h.

Referenced by CleanUp(), FluxParticles(), GenerateNext(), Initialize(), LoadFluxData(), and SetFlavors().

map<int, TH3D*> genie::flux::GPowerSpectrumAtmoFlux::fRawFluxHistoMap

private

flux = f(Ev,cos8,phi) for each neutrino species

Definition at line 79 of file GPowerSpectrumAtmoFlux.h.

Referenced by FillFluxHisto(), GetFlux(), and LoadFluxData().

map<int, TH2D*> genie::flux::GPowerSpectrumAtmoFlux::fRawFluxHistoMap2D

private

flux = f(Ev,cos8) for each neutrino species

Definition at line 80 of file GPowerSpectrumAtmoFlux.h.

Referenced by FillFluxHisto(), and GetFlux().

double genie::flux::GPowerSpectrumAtmoFlux::fRl

private

defining flux neutrino generation surface: longitudinal radius

Definition at line 71 of file GPowerSpectrumAtmoFlux.h.

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

TRotation genie::flux::GPowerSpectrumAtmoFlux::fRotTHz2User

private

coord. system rotation: THZ -> Topocentric user-defined

Definition at line 75 of file GPowerSpectrumAtmoFlux.h.

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

double genie::flux::GPowerSpectrumAtmoFlux::fRt

private

defining flux neutrino generation surface: transverse radius

Definition at line 72 of file GPowerSpectrumAtmoFlux.h.

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

double genie::flux::GPowerSpectrumAtmoFlux::fSpectralIndex

private

power law function

Definition at line 70 of file GPowerSpectrumAtmoFlux.h.

Referenced by ComputeWeight(), GenerateNext(), Initialize(), InitializeWeight(), and SetSpectralIndex().

double genie::flux::GPowerSpectrumAtmoFlux::fWeight

private

current generated nu weight

Definition at line 69 of file GPowerSpectrumAtmoFlux.h.

Referenced by GenerateNext(), ResetSelection(), and Weight().

---

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

• nugen/EventGeneratorBase/GENIE/GPowerSpectrumAtmoFlux.h
• nugen/EventGeneratorBase/GENIE/GPowerSpectrumAtmoFlux.cxx