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. More...

int	pdg
	HNL PDG code. More...

int	parPdg
	parent PDG code More...

int	lepPdg
	PDG code of lepton produced with HNL on parent decay. More...

int	nuPdg
	PDG code of SM neutrino that would have been produced. More...

int	prodChan
	Decay mode that produced HNL. More...

int	nuProdChan
	Decay mode that would have produced SM neutrino. More...

TVector3	startPoint
	parent decay vertex in NEAR coords [m] More...

TVector3	targetPoint
	point in detector HNL is forced towards in NEAR coords [m] More...

TVector3	startPointUser
	parent decay vertex in USER coords [m] More...

TVector3	targetPointUser
	point in detector HNL is forced towards in USER coords [m] More...

double	delay
	delay HNL would have wrt SMv [ns] More...

TVector3	polz
	HNL polarisation vector, in HNL rest frame, in NEAR coords. More...

TLorentzVector	p4
	HNL momentum in NEAR coords [GeV/c]. More...

TLorentzVector	parp4
	parent momentum at HNL production in NEAR coords [GeV/c] More...

TLorentzVector	p4User
	HNL momentum in USER coords [GeV/c]. More...

TLorentzVector	parp4User
	parent momentum at HNL production in USER coords [GeV/c] More...

double	Ecm
	Parent rest-frame energy of HNL [GeV]. More...

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

double	XYWgt
	geometric acceptance (angular size of detector in parent rest frame) More...

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

double	accCorr
	acceptance correction (collimation effect. SM v == 1) More...

double	zetaMinus
	minimum angular deviation from parent momentum to reach detector [deg] More...

double	zetaPlus
	maximum angular deviation from parent momentum to reach detector [deg] More...

double	acceptance
	full acceptance == XYWgt * boostCorr^2 * accCorr More...

double	nimpwt
	Weight of parent. More...

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:: Feb 16, 2023

License:: Copyright (c) 2003-2024, 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 ( )

Definition at line 61 of file HNLFluxContainer.cxx.

References ResetCopy().

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

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

inlinevirtual

Definition at line 66 of file HNLFluxContainer.h.

66 {};

Member Function Documentation

void FluxContainer::Print ( const Option_t * ) const

Definition at line 107 of file HNLFluxContainer.cxx.

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

void FluxContainer::ResetCopy ( ) const

Definition at line 68 of file HNLFluxContainer.cxx.

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().

 {
   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;
 }

Friends And Related Function Documentation

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;
   }

Member Data Documentation

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), genie::hnl::FluxCreator::ProcessEventRecord(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), genie::hnl::FluxCreator::ProcessEventRecord(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), genie::hnl::FluxCreator::ProcessEventRecord(), and ResetCopy().

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(), genie::hnl::operator<<(), genie::hnl::FluxCreator::ProcessEventRecord(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), genie::hnl::FluxCreator::ProcessEventRecord(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), genie::hnl::FluxCreator::ProcessEventRecord(), and ResetCopy().

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(), genie::hnl::operator<<(), genie::hnl::FluxCreator::ProcessEventRecord(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), genie::hnl::FluxCreator::ProcessEventRecord(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::operator<<(), and ResetCopy().

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(), genie::hnl::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

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation