A GENIE flux container specific for HNL containers. Based on the dk2nu flux paradigm and genie::flux::GNuMIFluxPassThroughInfo. More...

#include <HNLFluxContainer.h>

Inheritance diagram for genie::hnl::FluxContainer:

[legend]

Collaboration diagram for genie::hnl::FluxContainer:

[legend]

Public Member Functions
	FluxContainer ()
virtual	~FluxContainer ()
void	ResetCopy () const
void	Print (const Option_t *) const

Public Attributes
int	evtno
	Event number.
int	pdg
	HNL PDG code.
int	parPdg
	parent PDG code
int	lepPdg
	PDG code of lepton produced with HNL on parent decay.
int	nuPdg
	PDG code of SM neutrino that would have been produced.
int	prodChan
	Decay mode that produced HNL.
int	nuProdChan
	Decay mode that would have produced SM neutrino.
TVector3	startPoint
	parent decay vertex in NEAR coords [m]
TVector3	targetPoint
	point in detector HNL is forced towards in NEAR coords [m]
TVector3	startPointUser
	parent decay vertex in USER coords [m]
TVector3	targetPointUser
	point in detector HNL is forced towards in USER coords [m]
double	delay
	delay HNL would have wrt SMv [ns]
TVector3	polz
	HNL polarisation vector, in HNL rest frame, in NEAR coords.
TLorentzVector	p4
	HNL momentum in NEAR coords [GeV/c].
TLorentzVector	parp4
	parent momentum at HNL production in NEAR coords [GeV/c]
TLorentzVector	p4User
	HNL momentum in USER coords [GeV/c].
TLorentzVector	parp4User
	parent momentum at HNL production in USER coords [GeV/c]
double	Ecm
	Parent rest-frame energy of HNL [GeV].
double	nuEcm
	Parent rest-frame energy of equivalent SM neutrino [GeV].
double	XYWgt
	geometric acceptance (angular size of detector in parent rest frame)
double	boostCorr
	boost correction wrt parent rest-frame (ELAB = ECM * boostCorr)
double	accCorr
	acceptance correction (collimation effect. SM v == 1)
double	zetaMinus
	minimum angular deviation from parent momentum to reach detector [deg]
double	zetaPlus
	maximum angular deviation from parent momentum to reach detector [deg]
double	acceptance
	full acceptance == XYWgt * boostCorr^2 * accCorr
double	nimpwt
	Weight of parent.

Friends
ostream &	operator<< (ostream &stream, const FluxContainer &gnmf)

Detailed Description

A GENIE flux container specific for HNL containers. Based on the dk2nu flux paradigm and genie::flux::GNuMIFluxPassThroughInfo.

FluxContainer:

A C-struct that is based on the flux::GNuMIFluxPassThroughInfo struct, but which is not GNuMI specific. Accepts enough information about the HNL fluxes and the base flux to pass all the necessary

information to the hnl::FluxCreator class for flux calculations.

Also see $GENIE/src/contrib/beamhnl/write_dk2nus.C for an example of expected flux-input structure

Author: John Plows

Created:\n Feb 16, 2023

License:\n Copyright (c) 2003-2025, The GENIE Collaboration: For the full text of the license visit http://copyright.genie-mc.org

Definition at line 62 of file HNLFluxContainer.h.

Constructor & Destructor Documentation

◆ FluxContainer()

FluxContainer::FluxContainer ( )

Definition at line 61 of file HNLFluxContainer.cxx.

  : TObject()
{
  this->ResetCopy();
}

References ResetCopy().

Referenced by operator<<.

◆ ~FluxContainer()

virtual genie::hnl::FluxContainer::~FluxContainer ( )

inlinevirtual

Definition at line 66 of file HNLFluxContainer.h.

66{};

Member Function Documentation

◆ Print()

void FluxContainer::Print ( const Option_t * ) const

Definition at line 107 of file HNLFluxContainer.cxx.

{
  std::cout << *this << std::endl;
}

◆ ResetCopy()

void FluxContainer::ResetCopy ( ) const

Definition at line 68 of file HNLFluxContainer.cxx.

{
  evtno = 0;
 
  pdg = 0;
  parPdg = 0;
  lepPdg = 0;
  nuPdg = 0;
 
  prodChan = 0;
  nuProdChan = 0;
 
  startPoint.SetXYZ(0.0, 0.0, 0.0);
  targetPoint.SetXYZ(0.0, 0.0, 0.0);
  startPointUser.SetXYZ(0.0, 0.0, 0.0);
  targetPointUser.SetXYZ(0.0, 0.0, 0.0);
  delay = 0.0;
 
  polz.SetXYZ(0.0, 0.0, 0.0);
 
  p4.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
  parp4.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
  p4User.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
  parp4User.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
 
  Ecm = 0.0;
  nuEcm = 0.0;
 
  XYWgt = 0.0;
  boostCorr = 0.0;
  accCorr = 0.0;
  zetaMinus = 0.0;
  zetaPlus = 0.0;
  acceptance = 0.0;
 
  nimpwt = 0.0;
}

References accCorr, acceptance, boostCorr, delay, Ecm, evtno, lepPdg, nimpwt, nuEcm, nuPdg, nuProdChan, p4, p4User, parp4, parp4User, parPdg, pdg, polz, prodChan, startPoint, startPointUser, targetPoint, targetPointUser, XYWgt, zetaMinus, and zetaPlus.

Referenced by FluxContainer().

◆ operator<<

ostream & operator<<	(	ostream &	stream,
		const FluxContainer &	gnmf )

friend

Definition at line 115 of file HNLFluxContainer.cxx.

    {
      // convert some stuff to string
      HNLProd_t hProdChan = static_cast<HNLProd_t>(info.prodChan);
 
      std::string sNuProdChan;
      int typeMod = (info.pdg > 0) ? 1 : -1;
      int switchChan = typeMod * info.nuProdChan;
      switch(switchChan){
      case 1:  sNuProdChan = std::string("K0L -> nue + pi- + e+"); break;
      case -1: sNuProdChan = std::string("K0L -> nuebar + pi+ + e-"); break;
      case 2:  sNuProdChan = std::string("K0L -> numu + pi- + mu"); break;
      case -2: sNuProdChan = std::string("K0L -> numubar + pi+ + mu-"); break;
      case 3:  sNuProdChan = std::string("K+ -> numu + mu+"); break;
      case -3: sNuProdChan = std::string("K- -> numubar + mu-"); break;
      case 4:  sNuProdChan = std::string("K+ -> nue + e+"); break;
      case -4: sNuProdChan = std::string("K- -> nuebar + e-"); break;
      case 5:  sNuProdChan = std::string("K+ -> numu + pi0 + mu+"); break;
      case -5: sNuProdChan = std::string("K- -> numubar + pi0 + mu-"); break;
      case 6:  sNuProdChan = std::string("K+ -> nue + pi0 + e+"); break;
      case -6: sNuProdChan = std::string("K- -> nuebar + pi0 + e-"); break;
      case 7:  sNuProdChan = std::string("pi+ -> numu + mu+"); break;
      case -7: sNuProdChan = std::string("pi- -> numubar + mu-"); break;
      case 8:  sNuProdChan = std::string("pi+ -> nue + e+"); break;
      case -8: sNuProdChan = std::string("pi- -> nuebar + e-"); break;
      case 9:  sNuProdChan = std::string("mu- -> numu + nuebar + e-"); break;
      case -9: sNuProdChan = std::string("mu+ -> numubar + nue + e+"); break;
      default: sNuProdChan = std::string("Unknown!"); break;
      }
      
      stream << "\nEvent number: " << info.evtno
             << "\nHNL    PDG code: " << info.pdg
             << "\nParent PDG code: " << info.parPdg
             << "\nCo-produced lepton PDG code: " << info.lepPdg
             << "\nParent weight: " << info.nimpwt
             << "\nHNL polarisation vector [HNL rest frame, NEAR coords]: " << utils::print::Vec3AsString(&info.polz)
             << "\nPDG code of equivalent SM neutrino: " << info.nuPdg
             << "\nProduction channel: " << utils::hnl::ProdAsString(hProdChan)
             << " ; code " << info.prodChan
             << "\nEquivalent neutrino production channel: " << sNuProdChan
             << " ; code " << info.nuProdChan
             << "\nHNL parent rest-frame energy [GeV]: " << info.Ecm
             << "\nEquivalent neutrino parent rest-frame energy [GeV]: " << info.nuEcm
             << "\nStart point [NEAR, m]: " << utils::print::Vec3AsString(&info.startPoint)
             << "\nStart point [USER, m]: " << utils::print::Vec3AsString(&info.startPointUser)
             << "\nFlux passes through point [NEAR, m]: " << utils::print::Vec3AsString(&info.targetPoint)
             << "\nFlux passes through point [USER, m]: " << utils::print::Vec3AsString(&info.targetPointUser)
             << "\nHNL momentum [NEAR, GeV]: " << utils::print::P4AsString(&info.p4)
             << "\nHNL momentum [USER, GeV]: " << utils::print::P4AsString(&info.p4User)
             << "\nHNL delay wrt SM neutrino [ns]: " << info.delay
             << "\nParent momentum [NEAR, GeV]: " << utils::print::P4AsString(&info.parp4)
             << "\nParent momentum [USER, GeV]: " << utils::print::P4AsString(&info.parp4User)
             << "\nDeviation angles zeta- = " << info.zetaMinus << ", zeta+ = " << info.zetaPlus
             << "\nGeometric acceptance: " << info.XYWgt
             << "\nBoost correction: " << info.boostCorr
             << "\nAcceptance correction: " << info.accCorr
             << "\nFull acceptance: " << info.acceptance;
 
     return stream;
  }

References accCorr, acceptance, boostCorr, delay, Ecm, evtno, FluxContainer(), lepPdg, nimpwt, nuEcm, nuPdg, nuProdChan, p4, genie::utils::print::P4AsString(), p4User, parp4, parp4User, parPdg, pdg, polz, genie::utils::hnl::ProdAsString(), prodChan, startPoint, startPointUser, targetPoint, targetPointUser, genie::utils::print::Vec3AsString(), XYWgt, zetaMinus, and zetaPlus.

Member Data Documentation

◆ accCorr

double genie::hnl::FluxContainer::accCorr

mutable

acceptance correction (collimation effect. SM v == 1)

Definition at line 104 of file HNLFluxContainer.h.

Referenced by ConvertToGRooTracker(), genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ acceptance

double genie::hnl::FluxContainer::acceptance

mutable

full acceptance == XYWgt * boostCorr^2 * accCorr

Definition at line 108 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ boostCorr

double genie::hnl::FluxContainer::boostCorr

mutable

boost correction wrt parent rest-frame (ELAB = ECM * boostCorr)

Definition at line 102 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ delay

double genie::hnl::FluxContainer::delay

mutable

delay HNL would have wrt SMv [ns]

Definition at line 89 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ Ecm

double genie::hnl::FluxContainer::Ecm

mutable

Parent rest-frame energy of HNL [GeV].

Definition at line 98 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ evtno

int genie::hnl::FluxContainer::evtno

mutable

Event number.

Definition at line 75 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ lepPdg

int genie::hnl::FluxContainer::lepPdg

mutable

PDG code of lepton produced with HNL on parent decay.

Definition at line 79 of file HNLFluxContainer.h.

Referenced by ConvertToGRooTracker(), genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ nimpwt

double genie::hnl::FluxContainer::nimpwt

mutable

Weight of parent.

Definition at line 110 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ nuEcm

double genie::hnl::FluxContainer::nuEcm

mutable

Parent rest-frame energy of equivalent SM neutrino [GeV].

Definition at line 99 of file HNLFluxContainer.h.

Referenced by ConvertToGRooTracker(), genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ nuPdg

int genie::hnl::FluxContainer::nuPdg

mutable

PDG code of SM neutrino that would have been produced.

Definition at line 80 of file HNLFluxContainer.h.

Referenced by ConvertToGRooTracker(), genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ nuProdChan

int genie::hnl::FluxContainer::nuProdChan

mutable

Decay mode that would have produced SM neutrino.

Definition at line 83 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ p4

TLorentzVector genie::hnl::FluxContainer::p4

mutable

HNL momentum in NEAR coords [GeV/c].

Definition at line 93 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ p4User

TLorentzVector genie::hnl::FluxContainer::p4User

mutable

HNL momentum in USER coords [GeV/c].

Definition at line 95 of file HNLFluxContainer.h.

Referenced by ConvertToGRooTracker(), genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ parp4

TLorentzVector genie::hnl::FluxContainer::parp4

mutable

parent momentum at HNL production in NEAR coords [GeV/c]

Definition at line 94 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ parp4User

TLorentzVector genie::hnl::FluxContainer::parp4User

mutable

parent momentum at HNL production in USER coords [GeV/c]

Definition at line 96 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ parPdg

int genie::hnl::FluxContainer::parPdg

mutable

parent PDG code

Definition at line 78 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ pdg

int genie::hnl::FluxContainer::pdg

mutable

HNL PDG code.

Definition at line 77 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ polz

TVector3 genie::hnl::FluxContainer::polz

mutable

HNL polarisation vector, in HNL rest frame, in NEAR coords.

Definition at line 91 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ prodChan

int genie::hnl::FluxContainer::prodChan

mutable

Decay mode that produced HNL.

Definition at line 82 of file HNLFluxContainer.h.

Referenced by ConvertToGRooTracker(), genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ startPoint

TVector3 genie::hnl::FluxContainer::startPoint

mutable

parent decay vertex in NEAR coords [m]

Definition at line 85 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ startPointUser

TVector3 genie::hnl::FluxContainer::startPointUser

mutable

parent decay vertex in USER coords [m]

Definition at line 87 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ targetPoint

TVector3 genie::hnl::FluxContainer::targetPoint

mutable

point in detector HNL is forced towards in NEAR coords [m]

Definition at line 86 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ targetPointUser

TVector3 genie::hnl::FluxContainer::targetPointUser

mutable

point in detector HNL is forced towards in USER coords [m]

Definition at line 88 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ XYWgt

double genie::hnl::FluxContainer::XYWgt

mutable

geometric acceptance (angular size of detector in parent rest frame)

Definition at line 101 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ zetaMinus

double genie::hnl::FluxContainer::zetaMinus

mutable

minimum angular deviation from parent momentum to reach detector [deg]

Definition at line 105 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

◆ zetaPlus

double genie::hnl::FluxContainer::zetaPlus

mutable

maximum angular deviation from parent momentum to reach detector [deg]

Definition at line 106 of file HNLFluxContainer.h.

Referenced by genie::hnl::FluxCreator::FillNonsense(), genie::hnl::FluxCreator::MakeTupleFluxEntry(), operator<<, and ResetCopy().

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

Public Member Functions

Public Attributes

Friends

Detailed Description

FluxContainer:

information to the hnl::FluxCreator class for flux calculations.

Constructor & Destructor Documentation

◆ FluxContainer()

◆ ~FluxContainer()

Member Function Documentation

◆ Print()

◆ ResetCopy()

Friends And Related Symbol Documentation

◆ operator<<

Member Data Documentation

◆ accCorr

◆ acceptance

◆ boostCorr

◆ delay

◆ Ecm

◆ evtno

◆ lepPdg

◆ nimpwt

◆ nuEcm

◆ nuPdg

◆ nuProdChan

◆ p4

◆ p4User

◆ parp4

◆ parp4User

◆ parPdg

◆ pdg

◆ polz

◆ prodChan

◆ startPoint

◆ startPointUser

◆ targetPoint

◆ targetPointUser

◆ XYWgt

◆ zetaMinus

◆ zetaPlus