GENIEGenerator
Loading...
Searching...
No Matches
Public Member Functions | Protected Member Functions | Protected Attributes | Private Attributes | List of all members
genie::hnl::SimpleHNL Class Reference

HNL object. More...

#include <SimpleHNL.h>

Public Member Functions

 SimpleHNL (std::string name, int index)
 SimpleHNL (std::string name, int index, const int PDG, const int parPDG, const double mass, const double Ue42, const double Umu42, const double Ut42, const bool IsMajorana)
 default c'tor
 ~SimpleHNL ()
 normal constructor
std::string GetName ()
int GetIndex ()
double GetMass ()
std::vector< double > GetCouplings ()
bool GetIsMajorana ()
double GetBeta ()
double GetGamma ()
double GetCoMLifetime ()
double GetLifetime ()
int GetPDG ()
int GetParentPDG ()
double GetDecayThrow ()
double GetSelectThrow ()
genie::hnl::HNLDecayMode_t GetDecayMode ()
std::vector< double > GetDecay4VX ()
std::vector< double > GetOrigin4VX ()
std::vector< double > Get4VP ()
std::vector< double > GetBetaVec ()
double GetMomentum ()
double GetPolarisationMag ()
std::vector< double > * GetPolarisationDir ()
std::map< genie::hnl::HNLDecayMode_t, double > GetValidChannels ()
std::map< genie::hnl::HNLDecayMode_t, double > GetInterestingChannels ()
std::vector< genie::hnl::HNLDecayMode_t > GetInterestingChannelsVec ()
int GetType ()
double GetAngularDeviation ()
std::vector< double > GetBeam2UserTranslation ()
std::vector< double > GetBeam2UserRotation ()
void SetName (const std::string name)
void SetIndex (const int idx)
void SetEnergy (const double E)
void SetBeta (const double bet)
void SetMomentumAngles (double theta, double phi)
void SetMomentumDirection (double ux, double uy, double uz)
void Set4Momentum (const std::vector< double > fourP)
void SetPolMag (const double pm)
void SetPolarisationDirection (const double plx, const double ply, const double plz)
void SetProdVtx (const std::vector< double > fourV)
void SetDecayVtx (const std::vector< double > fourV)
void SetInterestingChannels (const std::map< genie::hnl::HNLDecayMode_t, double > gammaMap)
void SetInterestingChannelsVec (const std::vector< genie::hnl::HNLDecayMode_t > decVec)
void SetDecayMode (const genie::hnl::HNLDecayMode_t decayMode)
void SetParentPDG (const int parPDG)
void SetPDG (const int PDG)
void SetType (const int type)
void SetAngularDeviation (const double adev)
void SetBeam2UserTranslation (const double tx, const double ty, const double tz)
void SetBeam2UserRotation (const double r1, const double r2, const double r3)

Protected Member Functions

double CalcBeta (const double E, const double P3)
double CalcGamma (const double bet)
double CalcLifetime (const double gam)

Protected Attributes

int defPDG = 2000020000
int defParPDG = genie::kPdgKP
double defMass = 0.250 * genie::units::GeV
double defUe42 = 1.0e-4
double defUmu42 = 1.0e-4
double defUt42 = 0.0

Private Attributes

std::string fName
int fIndex
int fPDG
int fParentPDG
double fMass
double fUe42
double fUmu42
double fUt42
bool fIsMajorana
int fType
std::map< genie::hnl::HNLDecayMode_t, double > fValidChannels
double fCoMLifetime
double fDecayThrow
double fSelectThrow
genie::hnl::HNLDecayMode_t fDecayMode
std::map< genie::hnl::HNLDecayMode_t, double > fInterestingChannels
std::vector< genie::hnl::HNLDecayMode_t > fInterestingChannelsVec
double fBeta
double fGamma
double fLifetime
double fT
double fX
double fY
double fZ
double fT0
double fX0
double fY0
double fZ0
double fE
double fPmag
double fPx
double fPy
double fPz
double fPol
double fPolUx
double fPolUy
double fPolUz
std::vector< double > * fPolDir
double fAngularDeviation
double fTx
double fTy
double fTz
double fR1
double fR2
double fR3
double fRM11
double fRM12
double fRM13
double fRM21
double fRM22
double fRM23
double fRM31
double fRM32
double fRM33

Detailed Description

HNL object.

Author
John Plows komni.nosp@m.nos-.nosp@m.john..nosp@m.plow.nosp@m.s@phy.nosp@m.sics.nosp@m..ox.a.nosp@m.c.uk
Created:\n December 10th, 2021
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 36 of file SimpleHNL.h.

Constructor & Destructor Documentation

◆ SimpleHNL() [1/2]

genie::hnl::SimpleHNL::SimpleHNL ( std::string name,
int index )
inline

Definition at line 39 of file SimpleHNL.h.

39 :
40 fName( name ), fIndex( index ),
41 fPDG( defPDG ), fParentPDG( defParPDG ), fMass( defMass ),
42 fUe42( defUe42 ), fUmu42( defUmu42 ), fUt42( defUt42 ),
43 fIsMajorana( false ),
44 fValidChannels(
45 genie::hnl::selector::GetValidChannelWidths( defMass,
46 defUe42, defUmu42, defUt42,
47 false ) ),
48 fCoMLifetime(
49 genie::hnl::selector::CalcCoMLifetime( defMass,
50 defUe42, defUmu42, defUt42,
51 false ) )
52 { } /// default c'tor
genie::hnl::SimpleHNL::fValidChannels
std::map< genie::hnl::HNLDecayMode_t, double > fValidChannels
Definition SimpleHNL.h:346
genie::hnl::selector::GetValidChannelWidths
std::map< genie::hnl::HNLDecayMode_t, double > GetValidChannelWidths(const double M, const double Ue42, const double Umu42, const double Ut42, const bool IsMajorana=false)
Definition HNLDecaySelector.cxx:18
genie::hnl::selector::CalcCoMLifetime
double CalcCoMLifetime(const double M, const double Ue42, const double Umu42, const double Ut42, const bool IsMajorana=false)
Definition HNLDecaySelector.cxx:194

References defMass, defParPDG, defPDG, defUe42, defUmu42, defUt42, fCoMLifetime, fIndex, fIsMajorana, fMass, fName, fParentPDG, fPDG, fUe42, fUmu42, fUt42, and fValidChannels.

◆ SimpleHNL() [2/2]

genie::hnl::SimpleHNL::SimpleHNL ( std::string name,
int index,
const int PDG,
const int parPDG,
const double mass,
const double Ue42,
const double Umu42,
const double Ut42,
const bool IsMajorana )
inline

default c'tor

Definition at line 54 of file SimpleHNL.h.

59 : fName( name ), fIndex( index ),
60 fPDG( PDG ), fParentPDG( parPDG ), fMass( mass ),
61 fUe42( Ue42 ), fUmu42( Umu42 ), fUt42( Ut42 ),
62 fIsMajorana( IsMajorana ),
63 fValidChannels(
64 genie::hnl::selector::GetValidChannelWidths( mass,
65 Ue42, Umu42, Ut42,
66 IsMajorana ) ),
67 fCoMLifetime(
68 genie::hnl::selector::CalcCoMLifetime( mass,
69 Ue42, Umu42, Ut42,
70 IsMajorana ) )
71 { LOG("HNL", pDEBUG) << "SimpleHNL built"; } /// normal constructor
pDEBUG
#define pDEBUG
Definition Messenger.h:63
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE,...
Definition Messenger.h:96

References fCoMLifetime, fIndex, fIsMajorana, fMass, fName, fParentPDG, fPDG, fUe42, fUmu42, fUt42, fValidChannels, LOG, and pDEBUG.

◆ ~SimpleHNL()

genie::hnl::SimpleHNL::~SimpleHNL ( )
inline

normal constructor

Definition at line 73 of file SimpleHNL.h.

73{ }

Member Function Documentation

◆ CalcBeta()

double genie::hnl::SimpleHNL::CalcBeta ( const double E,
const double P3 )
inlineprotected

Definition at line 327 of file SimpleHNL.h.

327 {
328 return P3 / E; }

Referenced by SetEnergy().

◆ CalcGamma()

double genie::hnl::SimpleHNL::CalcGamma ( const double bet)
inlineprotected

Definition at line 330 of file SimpleHNL.h.

330 {
331 return std::sqrt( 1.0 / ( 1.0 - bet*bet ) ); }

Referenced by SetBeta(), and SetEnergy().

◆ CalcLifetime()

double genie::hnl::SimpleHNL::CalcLifetime ( const double gam)
inlineprotected

Definition at line 333 of file SimpleHNL.h.

333 {
334 return fCoMLifetime * gam; } // rest-to-lab transf

References fCoMLifetime.

Referenced by GetLifetime().

◆ Get4VP()

std::vector< double > genie::hnl::SimpleHNL::Get4VP ( )
inline

Definition at line 127 of file SimpleHNL.h.

127 {
128 std::vector< double > momVec;
129 momVec.emplace_back(fE);
130 momVec.emplace_back(fPx);
131 momVec.emplace_back(fPy);
132 momVec.emplace_back(fPz);
133 return momVec; }

References fE, fPx, fPy, and fPz.

◆ GetAngularDeviation()

double genie::hnl::SimpleHNL::GetAngularDeviation ( )
inline

Definition at line 165 of file SimpleHNL.h.

165{ return fAngularDeviation; }

References fAngularDeviation.

◆ GetBeam2UserRotation()

std::vector< double > genie::hnl::SimpleHNL::GetBeam2UserRotation ( )
inline

Definition at line 172 of file SimpleHNL.h.

172 {
173 std::vector<double> rVec = { fR1, fR2, fR3 };
174 return rVec;
175 }

References fR1, fR2, and fR3.

◆ GetBeam2UserTranslation()

std::vector< double > genie::hnl::SimpleHNL::GetBeam2UserTranslation ( )
inline

Definition at line 167 of file SimpleHNL.h.

167 {
168 std::vector<double> tVec = { fTx, fTy, fTz };
169 return tVec;
170 }

References fTx, fTy, and fTz.

◆ GetBeta()

double genie::hnl::SimpleHNL::GetBeta ( )
inline

Definition at line 92 of file SimpleHNL.h.

92{ return fBeta; }

References fBeta.

◆ GetBetaVec()

std::vector< double > genie::hnl::SimpleHNL::GetBetaVec ( )
inline

Definition at line 135 of file SimpleHNL.h.

135 {
136 std::vector< double > betaVec;
137 const double mom = GetMomentum( );
138 const int pxMod = ( fPx < 0.0 ) ? -1 : 1;
139 const int pyMod = ( fPy < 0.0 ) ? -1 : 1;
140 const int pzMod = ( fPz < 0.0 ) ? -1 : 1;
141 betaVec.emplace_back( pxMod * fPx / fE * fPx / mom );
142 betaVec.emplace_back( pyMod * fPy / fE * fPy / mom );
143 betaVec.emplace_back( pzMod * fPz / fE * fPz / mom );
144 return betaVec; }

References fE, fPx, fPy, fPz, and GetMomentum().

◆ GetCoMLifetime()

double genie::hnl::SimpleHNL::GetCoMLifetime ( )
inline

Definition at line 96 of file SimpleHNL.h.

96{ return fCoMLifetime; }

References fCoMLifetime.

Referenced by genie::hnl::Decayer::LoadConfig(), main(), and SelectDecayMode().

◆ GetCouplings()

std::vector< double > genie::hnl::SimpleHNL::GetCouplings ( )
inline

Definition at line 83 of file SimpleHNL.h.

83 {
84 std::vector< double > coupVec;
85 coupVec.emplace_back( fUe42 );
86 coupVec.emplace_back( fUmu42 );
87 coupVec.emplace_back( fUt42 );
88 return coupVec; }

References fUe42, fUmu42, and fUt42.

◆ GetDecay4VX()

std::vector< double > genie::hnl::SimpleHNL::GetDecay4VX ( )
inline

Definition at line 111 of file SimpleHNL.h.

111 {
112 std::vector< double > decVec;
113 decVec.emplace_back(fT);
114 decVec.emplace_back(fX);
115 decVec.emplace_back(fY);
116 decVec.emplace_back(fZ);
117 return decVec; }

References fT, fX, fY, and fZ.

◆ GetDecayMode()

genie::hnl::HNLDecayMode_t genie::hnl::SimpleHNL::GetDecayMode ( )
inline

Definition at line 109 of file SimpleHNL.h.

109{ return fDecayMode; }
genie::hnl::SimpleHNL::fDecayMode
genie::hnl::HNLDecayMode_t fDecayMode
Definition SimpleHNL.h:351

References fDecayMode.

◆ GetDecayThrow()

double genie::hnl::SimpleHNL::GetDecayThrow ( )
inline

Definition at line 105 of file SimpleHNL.h.

105{ return fDecayThrow; }

References fDecayThrow.

◆ GetGamma()

double genie::hnl::SimpleHNL::GetGamma ( )
inline

Definition at line 94 of file SimpleHNL.h.

94{ return fGamma; }

References fGamma.

◆ GetIndex()

int genie::hnl::SimpleHNL::GetIndex ( )
inline

Definition at line 79 of file SimpleHNL.h.

79{ return fIndex; }

References fIndex.

◆ GetInterestingChannels()

std::map< genie::hnl::HNLDecayMode_t, double > genie::hnl::SimpleHNL::GetInterestingChannels ( )
inline

Definition at line 157 of file SimpleHNL.h.

157 {
158 return fInterestingChannels; }
genie::hnl::SimpleHNL::fInterestingChannels
std::map< genie::hnl::HNLDecayMode_t, double > fInterestingChannels
Definition SimpleHNL.h:353

References fInterestingChannels.

◆ GetInterestingChannelsVec()

std::vector< genie::hnl::HNLDecayMode_t > genie::hnl::SimpleHNL::GetInterestingChannelsVec ( )
inline

Definition at line 160 of file SimpleHNL.h.

160 {
161 return fInterestingChannelsVec; }
genie::hnl::SimpleHNL::fInterestingChannelsVec
std::vector< genie::hnl::HNLDecayMode_t > fInterestingChannelsVec
Definition SimpleHNL.h:354

References fInterestingChannelsVec.

◆ GetIsMajorana()

bool genie::hnl::SimpleHNL::GetIsMajorana ( )
inline

Definition at line 90 of file SimpleHNL.h.

90{ return fIsMajorana; }

References fIsMajorana.

◆ GetLifetime()

double genie::hnl::SimpleHNL::GetLifetime ( )
inline

Definition at line 98 of file SimpleHNL.h.

98 { /* return fLifetime; */
99 return CalcLifetime( fGamma ); }
genie::hnl::SimpleHNL::CalcLifetime
double CalcLifetime(const double gam)
Definition SimpleHNL.h:333

References CalcLifetime(), and fGamma.

◆ GetMass()

double genie::hnl::SimpleHNL::GetMass ( )
inline

Definition at line 81 of file SimpleHNL.h.

81{ return fMass; }

References fMass.

◆ GetMomentum()

double genie::hnl::SimpleHNL::GetMomentum ( )
inline

Definition at line 146 of file SimpleHNL.h.

146 {
147 return fPmag;
148 }

References fPmag.

Referenced by GetBetaVec(), and SetEnergy().

◆ GetName()

std::string genie::hnl::SimpleHNL::GetName ( )
inline

Definition at line 77 of file SimpleHNL.h.

77{ return fName; }

References fName.

◆ GetOrigin4VX()

std::vector< double > genie::hnl::SimpleHNL::GetOrigin4VX ( )
inline

Definition at line 119 of file SimpleHNL.h.

119 {
120 std::vector< double > oriVec;
121 oriVec.emplace_back(fT0);
122 oriVec.emplace_back(fX0);
123 oriVec.emplace_back(fY0);
124 oriVec.emplace_back(fZ0);
125 return oriVec; }

References fT0, fX0, fY0, and fZ0.

◆ GetParentPDG()

int genie::hnl::SimpleHNL::GetParentPDG ( )
inline

Definition at line 103 of file SimpleHNL.h.

103{ return fParentPDG; }

References fParentPDG.

◆ GetPDG()

int genie::hnl::SimpleHNL::GetPDG ( )
inline

Definition at line 101 of file SimpleHNL.h.

101{ return fPDG; }

References fPDG.

◆ GetPolarisationDir()

std::vector< double > * genie::hnl::SimpleHNL::GetPolarisationDir ( )
inline

Definition at line 151 of file SimpleHNL.h.

151 {
152 return fPolDir; }
genie::hnl::SimpleHNL::fPolDir
std::vector< double > * fPolDir
Definition SimpleHNL.h:363

References fPolDir.

◆ GetPolarisationMag()

double genie::hnl::SimpleHNL::GetPolarisationMag ( )
inline

Definition at line 150 of file SimpleHNL.h.

150{ return fPol; }

References fPol.

◆ GetSelectThrow()

double genie::hnl::SimpleHNL::GetSelectThrow ( )
inline

Definition at line 107 of file SimpleHNL.h.

107{ return fSelectThrow; }

References fSelectThrow.

◆ GetType()

int genie::hnl::SimpleHNL::GetType ( )
inline

Definition at line 163 of file SimpleHNL.h.

163{ return fType; }

References fType.

◆ GetValidChannels()

std::map< genie::hnl::HNLDecayMode_t, double > genie::hnl::SimpleHNL::GetValidChannels ( )
inline

Definition at line 154 of file SimpleHNL.h.

154 {
155 return fValidChannels; }

References fValidChannels.

Referenced by main(), SelectDecayMode(), and TestDecay().

◆ Set4Momentum()

void genie::hnl::SimpleHNL::Set4Momentum ( const std::vector< double > fourP)
inline

Definition at line 243 of file SimpleHNL.h.

243 {
244 SetEnergy( fourP.at(0) ); // also takes care of Pmag
245 SetMomentumDirection( fourP.at(1), fourP.at(2), fourP.at(3) ); }
genie::hnl::SimpleHNL::SetMomentumDirection
void SetMomentumDirection(double ux, double uy, double uz)
Definition SimpleHNL.h:231
genie::hnl::SimpleHNL::SetEnergy
void SetEnergy(const double E)
Definition SimpleHNL.h:183

References SetEnergy(), and SetMomentumDirection().

◆ SetAngularDeviation()

void genie::hnl::SimpleHNL::SetAngularDeviation ( const double adev)
inline

Definition at line 296 of file SimpleHNL.h.

296{ fAngularDeviation = adev; }

References fAngularDeviation.

Referenced by genie::hnl::Decayer::GetHNLInstance().

◆ SetBeam2UserRotation()

void genie::hnl::SimpleHNL::SetBeam2UserRotation ( const double r1,
const double r2,
const double r3 )
inline

Definition at line 302 of file SimpleHNL.h.

302 {
303 fR1 = r1; fR2 = r2; fR3 = r3;
304 // and the rotation matrix
305 fRM11 = std::cos( fR2 );
306 fRM12 = -std::cos( fR3 ) * std::sin( fR2 );
307 fRM13 = std::sin( fR2 ) * std::sin( fR3 );
308 fRM21 = std::cos( fR1 ) * std::sin( fR2 );
309 fRM22 = std::cos( fR1 ) * std::cos( fR2 ) * std::cos( fR3 ) - std::sin( fR1 ) * std::sin( fR3 );
310 fRM23 = -std::cos( fR3 ) * std::sin( fR1 ) - std::cos( fR1 ) * std::cos( fR2 ) * std::sin( fR3 );
311 fRM31 = std::sin( fR1 ) * std::sin( fR2 );
312 fRM32 = std::cos( fR1 ) * std::sin( fR3 ) + std::cos( fR2 ) * std::cos( fR3 ) * std::sin( fR1 );
313 fRM33 = std::cos( fR1 ) * std::cos( fR3 ) - std::cos( fR2 ) * std::sin( fR1 ) * std::sin( fR3 );
314 }

References fR1, fR2, fR3, fRM11, fRM12, fRM13, fRM21, fRM22, fRM23, fRM31, fRM32, and fRM33.

Referenced by genie::hnl::Decayer::GetHNLInstance().

◆ SetBeam2UserTranslation()

void genie::hnl::SimpleHNL::SetBeam2UserTranslation ( const double tx,
const double ty,
const double tz )
inline

Definition at line 298 of file SimpleHNL.h.

298 {
299 fTx = tx; fTy = ty; fTz = tz;
300 }

References fTx, fTy, and fTz.

Referenced by genie::hnl::Decayer::GetHNLInstance().

◆ SetBeta()

void genie::hnl::SimpleHNL::SetBeta ( const double bet)
inline

Definition at line 198 of file SimpleHNL.h.

198 {
199 fBeta = bet;
200 fGamma = CalcGamma( bet );
201 fE = fGamma * fMass;
202 double oldmom = fPmag;
203 fPmag = std::sqrt( fE*fE - fMass*fMass );
204
205 fPx *= fPmag / oldmom; fPy *= fPmag / oldmom; fPz *= fPmag / oldmom;
206
207 fLifetime = fCoMLifetime / ( fGamma * ( 1.0 + bet ) );
208 }
genie::hnl::SimpleHNL::CalcGamma
double CalcGamma(const double bet)
Definition SimpleHNL.h:330

References CalcGamma(), fBeta, fCoMLifetime, fE, fGamma, fLifetime, fMass, fPmag, fPx, fPy, and fPz.

◆ SetDecayMode()

void genie::hnl::SimpleHNL::SetDecayMode ( const genie::hnl::HNLDecayMode_t decayMode)
inline

Definition at line 285 of file SimpleHNL.h.

285 {
286 fDecayMode = decayMode; }

References fDecayMode.

◆ SetDecayVtx()

void genie::hnl::SimpleHNL::SetDecayVtx ( const std::vector< double > fourV)
inline

Definition at line 271 of file SimpleHNL.h.

271 {
272 fT = fourV.at(0); fX = fourV.at(1);
273 fY = fourV.at(2); fZ = fourV.at(3);
274 }

References fT, fX, fY, and fZ.

◆ SetEnergy()

void genie::hnl::SimpleHNL::SetEnergy ( const double E)
inline

Definition at line 183 of file SimpleHNL.h.

183 {
184 // updates beta, gamma, 4P, lifetime. Doesn't change angles.
185 if( E < fMass ) { LOG( "HNL", pERROR ) <<
186 "genie::hnl::SimpleHNL:: Set E too low." <<
187 "\nE = " << E << ", M = " << fMass; exit(3); }
188 double mom3 = std::sqrt( E*E - fMass*fMass );
189 __attribute__((unused)) double oldmom = GetMomentum( );
190 fPmag = mom3;
191 fPx = mom3 / std::sqrt(3.0); fPy = mom3 / std::sqrt(3.0); fPz = mom3 / std::sqrt(3.0);
192 fE = E;
193 fBeta = CalcBeta( E, mom3 );
194 fGamma = CalcGamma( fBeta );
195 fLifetime = fCoMLifetime / ( fGamma * ( 1.0 + fBeta ) );
196 }
pERROR
#define pERROR
Definition Messenger.h:59
genie::hnl::SimpleHNL::CalcBeta
double CalcBeta(const double E, const double P3)
Definition SimpleHNL.h:327

References CalcBeta(), CalcGamma(), fBeta, fCoMLifetime, fE, fGamma, fLifetime, fMass, fPmag, fPx, fPy, fPz, GetMomentum(), LOG, and pERROR.

Referenced by Set4Momentum(), and TestDecay().

◆ SetIndex()

void genie::hnl::SimpleHNL::SetIndex ( const int idx)
inline

Definition at line 181 of file SimpleHNL.h.

181{ fIndex = idx; }

References fIndex.

◆ SetInterestingChannels()

void genie::hnl::SimpleHNL::SetInterestingChannels ( const std::map< genie::hnl::HNLDecayMode_t, double > gammaMap)
inline

Definition at line 276 of file SimpleHNL.h.

277 {
278 fInterestingChannels = gammaMap; }

References fInterestingChannels.

Referenced by SelectDecayMode().

◆ SetInterestingChannelsVec()

void genie::hnl::SimpleHNL::SetInterestingChannelsVec ( const std::vector< genie::hnl::HNLDecayMode_t > decVec)
inline

Definition at line 281 of file SimpleHNL.h.

282 {
283 fInterestingChannelsVec = decVec; }

References fInterestingChannelsVec.

Referenced by genie::hnl::Decayer::GetHNLInstance().

◆ SetMomentumAngles()

void genie::hnl::SimpleHNL::SetMomentumAngles ( double theta,
double phi )
inline

does not change magnitude

Definition at line 210 of file SimpleHNL.h.

210 {
211 /// does not change magnitude
212 // bring angles into [0,\pi]*[0,2\pi)
213 if( std::abs( theta ) > genie::constants::kPi ) {
214 const int nabs = std::floor( std::abs( theta ) / genie::constants::kPi );
215 theta += ( theta < 0 ) ? nabs * genie::constants::kPi : -nabs * genie::constants::kPi; }
216 if( std::abs( phi ) > 2.0 * genie::constants::kPi ) {
217 const int nabs = std::floor( std::abs( phi ) / ( 2.0 * genie::constants::kPi ) );
218 phi += ( phi < 0 ) ? nabs * 2.0 * genie::constants::kPi : -nabs * 2.0 * genie::constants::kPi; }
219 phi += ( phi < 0 ) ? 2.0 * genie::constants::kPi : 0.0;
220
221 if( theta < 0 ){
222 double ap = ( phi < genie::constants::kPi ) ? genie::constants::kPi : -genie::constants::kPi;
223 theta *= -1.0; phi += ap;
224 }
225
226 fPx = fPmag * std::sin( theta ) * std::cos( phi );
227 fPy = fPmag * std::sin( theta ) * std::sin( phi );
228 fPz = fPmag * std::cos( theta );
229 }

References fPmag, fPx, fPy, fPz, and genie::constants::kPi.

◆ SetMomentumDirection()

void genie::hnl::SimpleHNL::SetMomentumDirection ( double ux,
double uy,
double uz )
inline

does not change magnitude

Definition at line 231 of file SimpleHNL.h.

231 {
232 /// does not change magnitude
233 if( ux == 0.0 && uy == 0.0 && uz == 0.0 ){
234 LOG( "HNL", pERROR ) <<
235 "genie::hnl::SimpleHNL::SetMomentumDirection:: " <<
236 "Zero vector entered. Exiting."; exit(3); }
237 const double umag = std::sqrt( ( ux*ux + uy*uy + uz*uz ) );
238 const double invu = 1.0 / umag;
239 ux *= invu; uy *= invu; uz *= invu;
240 fPx = fPmag * ux; fPy = fPmag * uy; fPz = fPmag * uz;
241 }

References fPmag, fPx, fPy, fPz, LOG, and pERROR.

Referenced by Set4Momentum(), and TestDecay().

◆ SetName()

void genie::hnl::SimpleHNL::SetName ( const std::string name)
inline

Definition at line 179 of file SimpleHNL.h.

179{ fName = name; }

References fName.

◆ SetParentPDG()

void genie::hnl::SimpleHNL::SetParentPDG ( const int parPDG)
inline

Definition at line 288 of file SimpleHNL.h.

288 {
289 fParentPDG = parPDG; }

References fParentPDG.

◆ SetPDG()

void genie::hnl::SimpleHNL::SetPDG ( const int PDG)
inline

Definition at line 291 of file SimpleHNL.h.

291 {
292 fPDG = PDG; }

References fPDG.

◆ SetPolarisationDirection()

void genie::hnl::SimpleHNL::SetPolarisationDirection ( const double plx,
const double ply,
const double plz )
inline

Definition at line 254 of file SimpleHNL.h.

255 {
256 if( plx == 0.0 && ply == 0.0 && plz == 0.0 ){
257 LOG( "HNL", pERROR ) <<
258 "genie::hnl::SimpleHNL::SetPolarisationDirection:: " <<
259 "Zero vector entered. Exiting."; exit(1); }
260 const double PM = std::sqrt( plx*plx + ply*ply * plz*plz );
261 fPolUx = plx / PM; fPolUy = ply / PM; fPolUz = plz / PM;
262 if( !fPolDir || fPolDir == 0 ) fPolDir = new std::vector< double >( );
263 else fPolDir->clear();
264 fPolDir->emplace_back( fPolUx ); fPolDir->emplace_back( fPolUy );
265 fPolDir->emplace_back( fPolUz ); }

References fPolDir, fPolUx, fPolUy, fPolUz, LOG, and pERROR.

◆ SetPolMag()

void genie::hnl::SimpleHNL::SetPolMag ( const double pm)
inline

Definition at line 247 of file SimpleHNL.h.

247 {
248 if( pm < -1.0 || pm > 1.0 ){
249 LOG( "HNL", pERROR ) <<
250 "genie::hnl::SimpleHNL::SetPolMag:: " <<
251 "Pol.vec. magnitude must be in [-1,1]. Exiting."; exit(3); }
252 fPol = pm; }

References fPol, LOG, and pERROR.

◆ SetProdVtx()

void genie::hnl::SimpleHNL::SetProdVtx ( const std::vector< double > fourV)
inline

Definition at line 267 of file SimpleHNL.h.

267 {
268 fT0 = fourV.at(0); fX0 = fourV.at(1);
269 fY0 = fourV.at(2); fZ0 = fourV.at(3); }

References fT0, fX0, fY0, and fZ0.

◆ SetType()

void genie::hnl::SimpleHNL::SetType ( const int type)
inline

Definition at line 294 of file SimpleHNL.h.

294{ fType = type; }

References fType.

Referenced by genie::hnl::Decayer::GetHNLInstance().

Member Data Documentation

◆ defMass

double genie::hnl::SimpleHNL::defMass = 0.250 * genie::units::GeV
protected

Definition at line 321 of file SimpleHNL.h.

Referenced by SimpleHNL().

◆ defParPDG

int genie::hnl::SimpleHNL::defParPDG = genie::kPdgKP
protected

Definition at line 320 of file SimpleHNL.h.

Referenced by SimpleHNL().

◆ defPDG

int genie::hnl::SimpleHNL::defPDG = 2000020000
protected

Definition at line 319 of file SimpleHNL.h.

Referenced by SimpleHNL().

◆ defUe42

double genie::hnl::SimpleHNL::defUe42 = 1.0e-4
protected

Definition at line 322 of file SimpleHNL.h.

Referenced by SimpleHNL().

◆ defUmu42

double genie::hnl::SimpleHNL::defUmu42 = 1.0e-4
protected

Definition at line 323 of file SimpleHNL.h.

Referenced by SimpleHNL().

◆ defUt42

double genie::hnl::SimpleHNL::defUt42 = 0.0
protected

Definition at line 324 of file SimpleHNL.h.

Referenced by SimpleHNL().

◆ fAngularDeviation

double genie::hnl::SimpleHNL::fAngularDeviation
mutableprivate

Definition at line 365 of file SimpleHNL.h.

Referenced by GetAngularDeviation(), and SetAngularDeviation().

◆ fBeta

double genie::hnl::SimpleHNL::fBeta
mutableprivate

Definition at line 356 of file SimpleHNL.h.

Referenced by GetBeta(), SetBeta(), and SetEnergy().

◆ fCoMLifetime

double genie::hnl::SimpleHNL::fCoMLifetime
mutableprivate

Definition at line 347 of file SimpleHNL.h.

Referenced by CalcLifetime(), GetCoMLifetime(), SetBeta(), SetEnergy(), SimpleHNL(), and SimpleHNL().

◆ fDecayMode

genie::hnl::HNLDecayMode_t genie::hnl::SimpleHNL::fDecayMode
mutableprivate

Definition at line 351 of file SimpleHNL.h.

Referenced by GetDecayMode(), and SetDecayMode().

◆ fDecayThrow

double genie::hnl::SimpleHNL::fDecayThrow
mutableprivate

Definition at line 349 of file SimpleHNL.h.

Referenced by GetDecayThrow().

◆ fE

double genie::hnl::SimpleHNL::fE
mutableprivate

Definition at line 360 of file SimpleHNL.h.

Referenced by Get4VP(), GetBetaVec(), SetBeta(), and SetEnergy().

◆ fGamma

double genie::hnl::SimpleHNL::fGamma
private

Definition at line 356 of file SimpleHNL.h.

Referenced by GetGamma(), GetLifetime(), SetBeta(), and SetEnergy().

◆ fIndex

int genie::hnl::SimpleHNL::fIndex
mutableprivate

Definition at line 339 of file SimpleHNL.h.

Referenced by GetIndex(), SetIndex(), SimpleHNL(), and SimpleHNL().

◆ fInterestingChannels

std::map< genie::hnl::HNLDecayMode_t, double > genie::hnl::SimpleHNL::fInterestingChannels
mutableprivate

Definition at line 353 of file SimpleHNL.h.

Referenced by GetInterestingChannels(), and SetInterestingChannels().

◆ fInterestingChannelsVec

std::vector< genie::hnl::HNLDecayMode_t > genie::hnl::SimpleHNL::fInterestingChannelsVec
mutableprivate

Definition at line 354 of file SimpleHNL.h.

Referenced by GetInterestingChannelsVec(), and SetInterestingChannelsVec().

◆ fIsMajorana

bool genie::hnl::SimpleHNL::fIsMajorana
mutableprivate

Definition at line 344 of file SimpleHNL.h.

Referenced by GetIsMajorana(), SimpleHNL(), and SimpleHNL().

◆ fLifetime

double genie::hnl::SimpleHNL::fLifetime
mutableprivate

Definition at line 357 of file SimpleHNL.h.

Referenced by SetBeta(), and SetEnergy().

◆ fMass

double genie::hnl::SimpleHNL::fMass
mutableprivate

Definition at line 342 of file SimpleHNL.h.

Referenced by GetMass(), SetBeta(), SetEnergy(), SimpleHNL(), and SimpleHNL().

◆ fName

std::string genie::hnl::SimpleHNL::fName
mutableprivate

Definition at line 338 of file SimpleHNL.h.

Referenced by GetName(), SetName(), SimpleHNL(), and SimpleHNL().

◆ fParentPDG

int genie::hnl::SimpleHNL::fParentPDG
mutableprivate

Definition at line 341 of file SimpleHNL.h.

Referenced by GetParentPDG(), SetParentPDG(), SimpleHNL(), and SimpleHNL().

◆ fPDG

int genie::hnl::SimpleHNL::fPDG
mutableprivate

Definition at line 340 of file SimpleHNL.h.

Referenced by GetPDG(), SetPDG(), SimpleHNL(), and SimpleHNL().

◆ fPmag

double genie::hnl::SimpleHNL::fPmag
private

Definition at line 360 of file SimpleHNL.h.

Referenced by GetMomentum(), SetBeta(), SetEnergy(), SetMomentumAngles(), and SetMomentumDirection().

◆ fPol

double genie::hnl::SimpleHNL::fPol
mutableprivate

Definition at line 361 of file SimpleHNL.h.

Referenced by GetPolarisationMag(), and SetPolMag().

◆ fPolDir

std::vector< double >* genie::hnl::SimpleHNL::fPolDir
mutableprivate

Definition at line 363 of file SimpleHNL.h.

Referenced by GetPolarisationDir(), and SetPolarisationDirection().

◆ fPolUx

double genie::hnl::SimpleHNL::fPolUx
mutableprivate

Definition at line 362 of file SimpleHNL.h.

Referenced by SetPolarisationDirection().

◆ fPolUy

double genie::hnl::SimpleHNL::fPolUy
private

Definition at line 362 of file SimpleHNL.h.

Referenced by SetPolarisationDirection().

◆ fPolUz

double genie::hnl::SimpleHNL::fPolUz
private

Definition at line 362 of file SimpleHNL.h.

Referenced by SetPolarisationDirection().

◆ fPx

double genie::hnl::SimpleHNL::fPx
private

Definition at line 360 of file SimpleHNL.h.

Referenced by Get4VP(), GetBetaVec(), SetBeta(), SetEnergy(), SetMomentumAngles(), and SetMomentumDirection().

◆ fPy

double genie::hnl::SimpleHNL::fPy
private

Definition at line 360 of file SimpleHNL.h.

Referenced by Get4VP(), GetBetaVec(), SetBeta(), SetEnergy(), SetMomentumAngles(), and SetMomentumDirection().

◆ fPz

double genie::hnl::SimpleHNL::fPz
private

Definition at line 360 of file SimpleHNL.h.

Referenced by Get4VP(), GetBetaVec(), SetBeta(), SetEnergy(), SetMomentumAngles(), and SetMomentumDirection().

◆ fR1

double genie::hnl::SimpleHNL::fR1
mutableprivate

Definition at line 367 of file SimpleHNL.h.

Referenced by GetBeam2UserRotation(), and SetBeam2UserRotation().

◆ fR2

double genie::hnl::SimpleHNL::fR2
private

Definition at line 367 of file SimpleHNL.h.

Referenced by GetBeam2UserRotation(), and SetBeam2UserRotation().

◆ fR3

double genie::hnl::SimpleHNL::fR3
private

Definition at line 367 of file SimpleHNL.h.

Referenced by GetBeam2UserRotation(), and SetBeam2UserRotation().

◆ fRM11

double genie::hnl::SimpleHNL::fRM11
mutableprivate

Definition at line 368 of file SimpleHNL.h.

Referenced by SetBeam2UserRotation().

◆ fRM12

double genie::hnl::SimpleHNL::fRM12
private

Definition at line 368 of file SimpleHNL.h.

Referenced by SetBeam2UserRotation().

◆ fRM13

double genie::hnl::SimpleHNL::fRM13
private

Definition at line 368 of file SimpleHNL.h.

Referenced by SetBeam2UserRotation().

◆ fRM21

double genie::hnl::SimpleHNL::fRM21
private

Definition at line 368 of file SimpleHNL.h.

Referenced by SetBeam2UserRotation().

◆ fRM22

double genie::hnl::SimpleHNL::fRM22
private

Definition at line 368 of file SimpleHNL.h.

Referenced by SetBeam2UserRotation().

◆ fRM23

double genie::hnl::SimpleHNL::fRM23
private

Definition at line 368 of file SimpleHNL.h.

Referenced by SetBeam2UserRotation().

◆ fRM31

double genie::hnl::SimpleHNL::fRM31
private

Definition at line 368 of file SimpleHNL.h.

Referenced by SetBeam2UserRotation().

◆ fRM32

double genie::hnl::SimpleHNL::fRM32
private

Definition at line 368 of file SimpleHNL.h.

Referenced by SetBeam2UserRotation().

◆ fRM33

double genie::hnl::SimpleHNL::fRM33
private

Definition at line 368 of file SimpleHNL.h.

Referenced by SetBeam2UserRotation().

◆ fSelectThrow

double genie::hnl::SimpleHNL::fSelectThrow
mutableprivate

Definition at line 350 of file SimpleHNL.h.

Referenced by GetSelectThrow().

◆ fT

double genie::hnl::SimpleHNL::fT
mutableprivate

Definition at line 358 of file SimpleHNL.h.

Referenced by GetDecay4VX(), and SetDecayVtx().

◆ fT0

double genie::hnl::SimpleHNL::fT0
mutableprivate

Definition at line 359 of file SimpleHNL.h.

Referenced by GetOrigin4VX(), and SetProdVtx().

◆ fTx

double genie::hnl::SimpleHNL::fTx
mutableprivate

Definition at line 366 of file SimpleHNL.h.

Referenced by GetBeam2UserTranslation(), and SetBeam2UserTranslation().

◆ fTy

double genie::hnl::SimpleHNL::fTy
private

Definition at line 366 of file SimpleHNL.h.

Referenced by GetBeam2UserTranslation(), and SetBeam2UserTranslation().

◆ fType

int genie::hnl::SimpleHNL::fType
mutableprivate

Definition at line 345 of file SimpleHNL.h.

Referenced by GetType(), and SetType().

◆ fTz

double genie::hnl::SimpleHNL::fTz
private

Definition at line 366 of file SimpleHNL.h.

Referenced by GetBeam2UserTranslation(), and SetBeam2UserTranslation().

◆ fUe42

double genie::hnl::SimpleHNL::fUe42
mutableprivate

Definition at line 343 of file SimpleHNL.h.

Referenced by GetCouplings(), SimpleHNL(), and SimpleHNL().

◆ fUmu42

double genie::hnl::SimpleHNL::fUmu42
private

Definition at line 343 of file SimpleHNL.h.

Referenced by GetCouplings(), SimpleHNL(), and SimpleHNL().

◆ fUt42

double genie::hnl::SimpleHNL::fUt42
private

Definition at line 343 of file SimpleHNL.h.

Referenced by GetCouplings(), SimpleHNL(), and SimpleHNL().

◆ fValidChannels

std::map< genie::hnl::HNLDecayMode_t, double > genie::hnl::SimpleHNL::fValidChannels
mutableprivate

Definition at line 346 of file SimpleHNL.h.

Referenced by GetValidChannels(), SimpleHNL(), and SimpleHNL().

◆ fX

double genie::hnl::SimpleHNL::fX
private

Definition at line 358 of file SimpleHNL.h.

Referenced by GetDecay4VX(), and SetDecayVtx().

◆ fX0

double genie::hnl::SimpleHNL::fX0
private

Definition at line 359 of file SimpleHNL.h.

Referenced by GetOrigin4VX(), and SetProdVtx().

◆ fY

double genie::hnl::SimpleHNL::fY
private

Definition at line 358 of file SimpleHNL.h.

Referenced by GetDecay4VX(), and SetDecayVtx().

◆ fY0

double genie::hnl::SimpleHNL::fY0
private

Definition at line 359 of file SimpleHNL.h.

Referenced by GetOrigin4VX(), and SetProdVtx().

◆ fZ

double genie::hnl::SimpleHNL::fZ
private

Definition at line 358 of file SimpleHNL.h.

Referenced by GetDecay4VX(), and SetDecayVtx().

◆ fZ0

double genie::hnl::SimpleHNL::fZ0
private

Definition at line 359 of file SimpleHNL.h.

Referenced by GetOrigin4VX(), and SetProdVtx().

The documentation for this class was generated from the following file:
Generated on for GENIEGenerator by doxygen 1.14.0