genie::InitialState Class Reference

Initial State information. More...

#include <InitialState.h>

Inheritance diagram for genie::InitialState:

Public Member Functions
	InitialState ()
	InitialState (int tgt_pdgc, int probe_pdgc)
	InitialState (int Z, int A, int probe_pdgc)
	InitialState (const Target &tgt, int probe_pdgc)
	InitialState (const InitialState &initial_state)
	InitialState (TRootIOCtor *)
	~InitialState ()
TParticlePDG *	Probe (void) const
int	ProbePdg (void) const
int	TgtPdg (void) const
const Target &	Tgt (void) const
Target *	TgtPtr (void) const
TLorentzVector *	GetTgtP4 (RefFrame_t rf=kRfLab) const
TLorentzVector *	GetProbeP4 (RefFrame_t rf=kRfHitNucRest) const
double	ProbeE (RefFrame_t rf) const
double	CMEnergy () const
	centre-of-mass energy (sqrt s) More...
void	SetPdgs (int tgt_pdgc, int probe_pdgc)
void	SetProbePdg (int pdg_code)
void	SetTgtPdg (int pdg_code)
void	SetTgtP4 (const TLorentzVector &P4)
void	SetProbeP4 (const TLorentzVector &P4)
void	SetProbeE (double E)
bool	IsNuP (void) const
	is neutrino + proton? More...
bool	IsNuN (void) const
	is neutrino + neutron? More...
bool	IsNuBarP (void) const
	is anti-neutrino + proton? More...
bool	IsNuBarN (void) const
	is anti-neutrino + neutron? More...
bool	IsDMP (void) const
	is dark matter + proton? More...
bool	IsDMN (void) const
	is dark matter + neutron? More...
bool	IsDMBP (void) const
	is anti-dark matter + proton? More...
bool	IsDMBN (void) const
	is anti-dark matter + neutron? More...
void	Reset (void)
void	Copy (const InitialState &init_state)
bool	Compare (const InitialState &init_state) const
string	AsString (void) const
void	Print (ostream &stream) const
bool	operator== (const InitialState &i) const
	equal? More...
InitialState &	operator= (const InitialState &i)
	copy More...

Private Member Functions
void	Init (void)
void	Init (int target_pdgc, int probe_pdgc)
void	CleanUp (void)

Private Attributes
int	fProbePdg
	probe PDG code More...
Target *	fTgt
	nuclear target More...
TLorentzVector *	fProbeP4
	probe 4-momentum in LAB-frame More...
TLorentzVector *	fTgtP4
	nuclear target 4-momentum in LAB-frame More...

Friends
ostream &	operator<< (ostream &stream, const InitialState &i)
	print More...

Detailed Description

Initial State information.

Author

Costas Andreopoulos <constantinos.andreopoulos cern.ch> University of Liverpool & STFC Rutherford Appleton Laboratory

Changes required to implement the GENIE Boosted Dark Matter module were installed by Josh Berger (Univ. of Wisconsin)

Other minor changes / additions and fixes were installed by: Andy Furmanski (Univ. of Manchester) Joe Johnston (Univ of Pittsburgh)

May 02, 2004

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

Definition at line 48 of file InitialState.h.

Constructor & Destructor Documentation

InitialState::InitialState

(

)

Definition at line 46 of file InitialState.cxx.

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

InitialState::InitialState	(	int	tgt_pdgc,
		int	probe_pdgc
	)

Definition at line 52 of file InitialState.cxx.

                                                          :
TObject()
{
  this->Init(target_pdgc, probe_pdgc);
}

InitialState::InitialState	(	int	Z,
		int	A,
		int	probe_pdgc
	)

Definition at line 58 of file InitialState.cxx.

                                                       :
TObject()
{
  int target_pdgc = pdg::IonPdgCode(A,Z);
  this->Init(target_pdgc, probe_pdgc);
}

InitialState::InitialState	(	const Target &	tgt,
		int	probe_pdgc
	)

Definition at line 65 of file InitialState.cxx.

                                                             :
TObject()
{
  int target_pdgc = tgt.Pdg();
  this->Init(target_pdgc, probe_pdgc);
}

InitialState::InitialState

(

const InitialState &

initial_state

)

Definition at line 72 of file InitialState.cxx.

                                                          :
TObject()
{
  this->Init();
  this->Copy(init_state);
}

InitialState::InitialState

(

TRootIOCtor *

)

Definition at line 79 of file InitialState.cxx.

                                       :
TObject(),
fProbePdg(0),
fTgt(0),
fProbeP4(0),
fTgtP4(0)
{

}

InitialState::~InitialState

(

)

Definition at line 89 of file InitialState.cxx.

{
  this->CleanUp();
}

Member Function Documentation

string InitialState::AsString

(

void

)

const

Definition at line 408 of file InitialState.cxx.

{
// Code-ify the interaction in a string to be used as (part of a) keys
// Template:
//     nu_pdg:code;tgt-pdg:code;

  ostringstream init_state;

  if (this->Probe()->Mass() > 0) {
    init_state << "dm_mass:" << this->Probe()->Mass() << ";";
  }
  else {
    init_state << "nu-pdg:"  << this->ProbePdg()  << ";";
  }
  init_state << "tgt-pdg:" << this->Tgt().Pdg() << ";";

  return init_state.str();
}

void InitialState::CleanUp ( void  )

private

Definition at line 118 of file InitialState.cxx.

{
  delete fTgt;
  delete fProbeP4;
  delete fTgtP4;
}

double InitialState::CMEnergy

(

)

const

centre-of-mass energy (sqrt s)

Definition at line 394 of file InitialState.cxx.

{
  TLorentzVector * k4 = this->GetProbeP4(kRfLab);
  TLorentzVector * p4 = fTgt->HitNucP4Ptr();

  *k4 += *p4; // now k4 represents centre-of-mass 4-momentum
  double s = k4->Dot(*k4); // dot-product with itself
  double E = TMath::Sqrt(s);

  delete k4;

  return E;
}

bool InitialState::Compare

(

const InitialState &

init_state

)

const

Definition at line 497 of file InitialState.cxx.

{
  int            probe  = init_state.ProbePdg();
  const Target & target = init_state.Tgt();

  bool equal = (fProbePdg == probe) && (*fTgt == target);

  return equal;
}

void InitialState::Copy

(

const InitialState &

init_state

)

Definition at line 131 of file InitialState.cxx.

{
  fProbePdg = init_state.fProbePdg;

  fTgt->Copy(*init_state.fTgt);

  this -> SetProbeP4 ( *init_state.fProbeP4 );
  this -> SetTgtP4   ( *init_state.fTgtP4   );
}

TLorentzVector * InitialState::GetProbeP4

(

RefFrame_t

rf = kRfHitNucRest

)

const

Definition at line 325 of file InitialState.cxx.

{
// Return the probe 4-momentum in the specified reference frame
// Note: the caller adopts the TLorentzVector object

  switch (ref_frame) {

       //------------------ NUCLEAR TARGET REST FRAME:
       case (kRfTgtRest) :
       {
             // for now assure that the target nucleus is always at rest
             // and it is only the struck nucleons that can move:
             // so the [target rest frame] = [LAB frame]

             TLorentzVector * p4 = new TLorentzVector(*fProbeP4);
             return p4;
       }

       //------------------ STRUCK NUCLEON REST FRAME:
       case (kRfHitNucRest) :
       {
             // make sure that 'struck nucleon' properties were set in
             // the nuclear target object

             assert( fTgt->HitNucP4Ptr() != 0 );

             TLorentzVector * pnuc4 = fTgt->HitNucP4Ptr();

             // compute velocity vector (px/E, py/E, pz/E)

             double bx = pnuc4->Px() / pnuc4->Energy();
             double by = pnuc4->Py() / pnuc4->Energy();
             double bz = pnuc4->Pz() / pnuc4->Energy();

             // BOOST

             TLorentzVector * p4 = new TLorentzVector(*fProbeP4);

             p4->Boost(-bx,-by,-bz);

             return p4;

             break;
       }
       //------------------ LAB:
       case (kRfLab) :
       {
             TLorentzVector * p4 = new TLorentzVector(*fProbeP4);
             return p4;

             break;
       }
       default:

             LOG("Interaction", pERROR) << "Uknown reference frame";
  }
  return 0;
}

TLorentzVector * InitialState::GetTgtP4

(

RefFrame_t

rf = kRfLab

)

const

Definition at line 275 of file InitialState.cxx.

{
// Return the target 4-momentum in the specified reference frame
// Note: the caller adopts the TLorentzVector object

  switch (ref_frame) {

       //------------------ NUCLEAR TARGET REST FRAME:
       case (kRfTgtRest) :
       {
             // for now make sure that the target nucleus is always at
             // rest and it is only the struck nucleons that can move:
             // so the [target rest frame] = [LAB frame]

             return this->GetTgtP4(kRfLab);
       }

       //------------------ STRUCK NUCLEON REST FRAME:
       case (kRfHitNucRest) :
       {
             // make sure that 'struck nucleon' properties were set in
             // the nuclear target object
             assert(fTgt->HitNucIsSet());
             TLorentzVector * pnuc4 = fTgt->HitNucP4Ptr();

             // compute velocity vector (px/E, py/E, pz/E)
             double bx = pnuc4->Px() / pnuc4->Energy();
             double by = pnuc4->Py() / pnuc4->Energy();
             double bz = pnuc4->Pz() / pnuc4->Energy();

             // BOOST
             TLorentzVector * p4 = new TLorentzVector(*fTgtP4);
             p4->Boost(-bx,-by,-bz);

             return p4;
             break;
       }
       //------------------ LAB:
       case (kRfLab) :
       {
             TLorentzVector * p4 = new TLorentzVector(*fTgtP4);
             return p4;
             break;
       }
       default:
             LOG("Interaction", pERROR) << "Uknown reference frame";
  }
  return 0;
}

void InitialState::Init ( void  )

private

Definition at line 94 of file InitialState.cxx.

{
  fProbePdg  = 0;
  fTgt       = new Target();
  fProbeP4   = new TLorentzVector(0, 0, 0, 0);
  fTgtP4     = new TLorentzVector(0, 0, 0, 0);
}

void InitialState::Init ( int  target_pdgc, int  probe_pdgc  )

private

Definition at line 102 of file InitialState.cxx.

{
  TParticlePDG * t = PDGLibrary::Instance()->Find(target_pdgc);
  TParticlePDG * p = PDGLibrary::Instance()->Find(probe_pdgc );

  assert(t && p);

  double m = p->Mass();
  double M = t->Mass();

  fProbePdg  = probe_pdgc;
  fTgt       = new Target(target_pdgc);
  fProbeP4   = new TLorentzVector(0, 0, 0, m);
  fTgtP4     = new TLorentzVector(0, 0, 0, M);
}

bool InitialState::IsDMBN

(

void

)

const

is anti-dark matter + neutron?

Definition at line 265 of file InitialState.cxx.

{
// Check if DM - neutron interaction
  int  prob = fProbePdg;
  int  nucl = fTgt->HitNucPdg();
  bool isdn = pdg::IsAntiDarkMatter(prob) && pdg::IsNeutron(nucl);

  return isdn;
}

bool InitialState::IsDMBP

(

void

)

const

is anti-dark matter + proton?

Definition at line 255 of file InitialState.cxx.

{
// Check if DM - proton interaction
  int  prob = fProbePdg;
  int  nucl = fTgt->HitNucPdg();
  bool isdp = pdg::IsAntiDarkMatter(prob) && pdg::IsProton(nucl);

  return isdp;
}

bool InitialState::IsDMN

(

void

)

const

is dark matter + neutron?

Definition at line 245 of file InitialState.cxx.

{
// Check if DM - neutron interaction
  int  prob = fProbePdg;
  int  nucl = fTgt->HitNucPdg();
  bool isdn = pdg::IsDarkMatter(prob) && pdg::IsNeutron(nucl);

  return isdn;
}

bool InitialState::IsDMP

(

void

)

const

is dark matter + proton?

Definition at line 235 of file InitialState.cxx.

{
// Check if DM - proton interaction
  int  prob = fProbePdg;
  int  nucl = fTgt->HitNucPdg();
  bool isdp = pdg::IsDarkMatter(prob) && pdg::IsProton(nucl);

  return isdp;
}

bool InitialState::IsNuBarN

(

void

)

const

is anti-neutrino + neutron?

Definition at line 226 of file InitialState.cxx.

{
  int  prob  = fProbePdg;
  int  nucl  = fTgt->HitNucPdg();
  bool isvbn = pdg::IsAntiNeutrino(prob) && pdg::IsNeutron(nucl);

  return isvbn;
}

bool InitialState::IsNuBarP

(

void

)

const

is anti-neutrino + proton?

Definition at line 217 of file InitialState.cxx.

{
  int  prob  = fProbePdg;
  int  nucl  = fTgt->HitNucPdg();
  bool isvbp = pdg::IsAntiNeutrino(prob) && pdg::IsProton(nucl);

  return isvbp;
}

bool InitialState::IsNuN

(

void

)

const

is neutrino + neutron?

Definition at line 208 of file InitialState.cxx.

{
  int  prob = fProbePdg;
  int  nucl = fTgt->HitNucPdg();
  bool isvn = pdg::IsNeutrino(prob) && pdg::IsNeutron(nucl);

  return isvn;
}

bool InitialState::IsNuP

(

void

)

const

is neutrino + proton?

Definition at line 199 of file InitialState.cxx.

{
  int  prob = fProbePdg;
  int  nucl = fTgt->HitNucPdg();
  bool isvp = pdg::IsNeutrino(prob) && pdg::IsProton(nucl);

  return isvp;
}

InitialState & InitialState::operator=

(

const InitialState &

i

)

copy

Definition at line 512 of file InitialState.cxx.

{
  this->Copy(init_state);
  return (*this);
}

bool InitialState::operator==

(

const InitialState &

i

)

const

equal?

Definition at line 507 of file InitialState.cxx.

{
  return this->Compare(init_state);
}

void InitialState::Print

(

ostream &

stream

)

const

Definition at line 427 of file InitialState.cxx.

{
  stream << "[-] [Init-State] " << endl;

  stream << " |--> probe        : "
         << "PDG-code = " << fProbePdg
         << " (" << this->Probe()->GetName() << ")" << endl;

  stream << " |--> nucl. target : "
         << "Z = "          << fTgt->Z()
         << ", A = "        << fTgt->A()
         << ", PDG-Code = " << fTgt->Pdg();

  TParticlePDG * tgt = PDGLibrary::Instance()->Find( fTgt->Pdg() );
  if(tgt) {
    stream << " (" << tgt->GetName() << ")";
  }
  stream << endl;

  stream << " |--> hit nucleon  : ";
  int nuc_pdgc = fTgt->HitNucPdg();

  if ( pdg::IsNeutronOrProton(nuc_pdgc) ) {
    TParticlePDG * p = PDGLibrary::Instance()->Find(nuc_pdgc);
    stream << "PDC-Code = " << nuc_pdgc << " (" << p->GetName() << ")";
  } else {
    stream << "no set";
  }
  stream << endl;

  stream << " |--> hit quark    : ";
  int qrk_pdgc = fTgt->HitQrkPdg();

  if ( pdg::IsQuark(qrk_pdgc) || pdg::IsAntiQuark(qrk_pdgc)) {
    TParticlePDG * p = PDGLibrary::Instance()->Find(qrk_pdgc);
    stream << "PDC-Code = " << qrk_pdgc << " (" << p->GetName() << ") ";
    stream << (fTgt->HitSeaQrk() ? "[sea]" : "[valence]");
  } else {
    stream << "no set";
  }
  stream << endl;

  stream << " |--> probe 4P     : "
         << "(E = "   << setw(12) << setprecision(6) << fProbeP4->E()
         << ", Px = " << setw(12) << setprecision(6) << fProbeP4->Px()
         << ", Py = " << setw(12) << setprecision(6) << fProbeP4->Py()
         << ", Pz = " << setw(12) << setprecision(6) << fProbeP4->Pz()
         << ")"
         << endl;
  stream << " |--> target 4P    : "
         << "(E = "   << setw(12) << setprecision(6) << fTgtP4->E()
         << ", Px = " << setw(12) << setprecision(6) << fTgtP4->Px()
         << ", Py = " << setw(12) << setprecision(6) << fTgtP4->Py()
         << ", Pz = " << setw(12) << setprecision(6) << fTgtP4->Pz()
         << ")"
         << endl;

  if ( pdg::IsNeutronOrProton(nuc_pdgc) ) {

    TLorentzVector * nuc_p4 = fTgt->HitNucP4Ptr();

    stream << " |--> nucleon 4P   : "
           << "(E = "   << setw(12) << setprecision(6) << nuc_p4->E()
           << ", Px = " << setw(12) << setprecision(6) << nuc_p4->Px()
           << ", Py = " << setw(12) << setprecision(6) << nuc_p4->Py()
           << ", Pz = " << setw(12) << setprecision(6) << nuc_p4->Pz()
           << ")";
  }
}

TParticlePDG * InitialState::Probe

(

void

)

const

Definition at line 147 of file InitialState.cxx.

{
  TParticlePDG * p = PDGLibrary::Instance()->Find(fProbePdg);
  return p;
}

double InitialState::ProbeE

(

RefFrame_t

rf

)

const

Definition at line 384 of file InitialState.cxx.

{
  TLorentzVector * p4 = this->GetProbeP4(ref_frame);
  double E = p4->Energy();

  delete p4;
  return E;
}

int genie::InitialState::ProbePdg ( void  ) const

inline

Definition at line 64 of file InitialState.h.

{ return fProbePdg; }

void InitialState::Reset

(

void

)

Definition at line 125 of file InitialState.cxx.

{
  this->CleanUp();
  this->Init();
}

void InitialState::SetPdgs	(	int	tgt_pdgc,
		int	probe_pdgc
	)

Definition at line 153 of file InitialState.cxx.

{
  this->CleanUp();
  this->Init(tgt_pdgc, probe_pdgc);
}

void InitialState::SetProbeE

(

double

E

)

Definition at line 175 of file InitialState.cxx.

{
  fProbeP4 -> SetE  ( E );
  fProbeP4 -> SetPx ( 0.);
  fProbeP4 -> SetPy ( 0.);
  fProbeP4 -> SetPz ( E );
}

void InitialState::SetProbeP4

(

const TLorentzVector &

P4

)

Definition at line 183 of file InitialState.cxx.

{
  fProbeP4 -> SetE  ( P4.E()  );
  fProbeP4 -> SetPx ( P4.Px() );
  fProbeP4 -> SetPy ( P4.Py() );
  fProbeP4 -> SetPz ( P4.Pz() );
}

void InitialState::SetProbePdg

(

int

pdg_code

)

Definition at line 167 of file InitialState.cxx.

{
  int tgt_pdgc = this->TgtPdg();

  this->CleanUp();
  this->Init(tgt_pdgc, probe_pdgc);
}

void InitialState::SetTgtP4

(

const TLorentzVector &

P4

)

Definition at line 191 of file InitialState.cxx.

{
  fTgtP4 -> SetE  ( P4.E()  );
  fTgtP4 -> SetPx ( P4.Px() );
  fTgtP4 -> SetPy ( P4.Py() );
  fTgtP4 -> SetPz ( P4.Pz() );
}

void InitialState::SetTgtPdg

(

int

pdg_code

)

Definition at line 159 of file InitialState.cxx.

{
  int probe_pdgc = this->ProbePdg();

  this->CleanUp();
  this->Init(tgt_pdgc, probe_pdgc);
}

const Target& genie::InitialState::Tgt ( void  ) const

inline

Definition at line 66 of file InitialState.h.

{ return *fTgt; }

int InitialState::TgtPdg

(

void

)

const

Definition at line 141 of file InitialState.cxx.

{
  assert(fTgt);
  return fTgt->Pdg();
}

Target* genie::InitialState::TgtPtr ( void  ) const

inline

Definition at line 67 of file InitialState.h.

{ return  fTgt; }

Friends And Related Function Documentation

ostream& operator<< ( ostream &  stream, const InitialState &  i  )

friend

print

Member Data Documentation

TLorentzVector* genie::InitialState::fProbeP4

private

probe 4-momentum in LAB-frame

Definition at line 111 of file InitialState.h.

int genie::InitialState::fProbePdg

private

probe PDG code

Definition at line 109 of file InitialState.h.

Target* genie::InitialState::fTgt

private

nuclear target

Definition at line 110 of file InitialState.h.

TLorentzVector* genie::InitialState::fTgtP4

private

nuclear target 4-momentum in LAB-frame

Definition at line 112 of file InitialState.h.

---

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

• Generator/src/Framework/Interaction/InitialState.h
• Generator/src/Framework/Interaction/InitialState.cxx