A Neutrino Interaction Target. Is a transparent encapsulation of quite different physical systems such as a nuclear target, a 'spectator' nuclear target with a Hit nucleon, a free nucleon or a free particle (eg a e- target in the inverse muon decay reaction) More...

#include <Target.h>

Inheritance diagram for genie::Target:

Public Member Functions
	Target ()

	Target (int pdgc)

	Target (int Z, int A)

	Target (int Z, int A, int hit_nucleon_pdgc)

	Target (const Target &tgt)

	Target (TRootIOCtor *)

	~Target ()

void	SetId (int pdgc)

void	SetId (int Z, int A)

void	SetHitNucPdg (int pdgc)

void	SetHitNucP4 (const TLorentzVector &p4)

void	SetHitNucPosition (double r)

void	SetHitQrkPdg (int pdgc)

void	SetHitSeaQrk (bool tf)

void	ForceHitNucOnMassShell (void)

int	Z (void) const

int	N (void) const

int	A (void) const

int	Pdg (void) const

double	Mass (void) const

double	Charge (void) const

bool	IsFreeNucleon (void) const

bool	IsProton (void) const

bool	IsNeutron (void) const

bool	IsNucleus (void) const

bool	IsParticle (void) const

bool	IsValidNucleus (void) const

bool	HitNucIsSet (void) const

bool	HitQrkIsSet (void) const

bool	HitSeaQrk (void) const

bool	IsEvenEven (void) const

bool	IsEvenOdd (void) const

bool	IsOddOdd (void) const

int	HitNucPdg (void) const

int	HitQrkPdg (void) const

double	HitNucMass (void) const

double	HitNucPosition (void) const

const TLorentzVector &	HitNucP4 (void) const

TLorentzVector *	HitNucP4Ptr (void) const

void	Reset (void)

void	Copy (const Target &t)

bool	Compare (const Target &t) const

string	AsString (void) const

void	Print (ostream &stream) const

bool	operator== (const Target &t) const
	equal? More...

Target &	operator= (const Target &t)
	copy More...

Private Member Functions
void	Init (void)

void	CleanUp (void)

void	ForceNucleusValidity (void)

bool	ForceHitNucValidity (void)

void	AutoSetHitNuc (void)

Private Attributes
int	fZ
	nuclear target Z More...

int	fA
	nuclear target A More...

int	fTgtPDG
	nuclear target PDG code More...

int	fHitNucPDG
	hit nucleon PDG code More...

int	fHitQrkPDG
	hit quark PDG code More...

bool	fHitSeaQrk
	hit quark from sea? More...

TLorentzVector *	fHitNucP4
	hit nucleon 4p More...

double	fHitNucRad
	hit nucleon position More...

Friends
ostream &	operator<< (ostream &stream, const Target &t)
	print More...

Detailed Description

A Neutrino Interaction Target. Is a transparent encapsulation of quite different physical systems such as a nuclear target, a 'spectator' nuclear target with a Hit nucleon, a free nucleon or a free particle (eg a e- target in the inverse muon decay reaction)

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

May 03, 2004

Definition at line 40 of file Target.h.

Constructor & Destructor Documentation

Target::Target ( )

Definition at line 41 of file Target.cxx.

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

Target::Target ( int pdgc )

Definition at line 47 of file Target.cxx.

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

Target::Target	(	int	Z,
		int	A
	)

Definition at line 54 of file Target.cxx.

                              :
 TObject()
 {
   this->Init();
   this->SetId(ZZ,AA);
 }

Target::Target	(	int	Z,
		int	A,
		int	hit_nucleon_pdgc
	)

Definition at line 61 of file Target.cxx.

                                                    :
 TObject()
 {
   this->Init();
   this->SetId(ZZ,AA);
   this->SetHitNucPdg(hit_nucleon_pdgc);
 }

Target::Target ( const Target & tgt )

Definition at line 69 of file Target.cxx.

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

Target::Target ( TRootIOCtor * )

Definition at line 76 of file Target.cxx.

                            :
 TObject(),
 fZ(0),
 fA(0),
 fTgtPDG(0),
 fHitNucPDG(0),
 fHitSeaQrk(false),
 fHitNucP4(0)
 {
 
 }

Target::~Target ( )

Definition at line 88 of file Target.cxx.

 {
   this->CleanUp();
 }

Member Function Documentation

int genie::Target::A ( void ) const

inline

Definition at line 70 of file Target.h.

70 { return fA; }

genie::Target::fA

int fA

nuclear target A

Definition: Target.h:118

string Target::AsString ( void ) const

Definition at line 383 of file Target.cxx.

 {
   ostringstream s;
 
   s << this->Pdg();
   if(this->HitNucIsSet())
      s << "[N=" << this->HitNucPdg() << "]";
   if(this->HitQrkIsSet()) {
      s << "[q=" << this->HitQrkPdg();
      s << (this->HitSeaQrk() ? "(s)" : "(v)");
      s << "]";
   }
 
   return s.str();
 }

void Target::AutoSetHitNuc ( void )

private

Definition at line 373 of file Target.cxx.

 {
 // for free nucleon targets -> (auto)set struck nucleon = target
 
   if( this->IsFreeNucleon() ) {
     if( this->IsProton() ) this->SetHitNucPdg(kPdgProton);
     else                   this->SetHitNucPdg(kPdgNeutron);
   }
 }

double Target::Charge ( void ) const

Definition at line 215 of file Target.cxx.

 {
 // Shortcut for commonly used code for extracting the nucleus charge from PDG
 //
   TParticlePDG * p = PDGLibrary::Instance()->Find(fTgtPDG);
   if(p) return p->Charge() / 3.; // in +e
   return 0;
 }

void Target::CleanUp ( void )

private

Definition at line 111 of file Target.cxx.

 {
   delete fHitNucP4;
 }

bool Target::Compare ( const Target & t ) const

Definition at line 422 of file Target.cxx.

 {
   int  tgt_pdg        = target.Pdg();
   int  struck_nuc_pdg = target.HitNucPdg();
   int  struck_qrk_pdg = target.HitQrkPdg();
   bool struck_sea_qrk = target.HitSeaQrk();
 
   bool equal = ( fTgtPDG    == tgt_pdg        ) &&
                ( fHitNucPDG == struck_nuc_pdg ) &&
                ( fHitQrkPDG == struck_qrk_pdg ) &&
                ( fHitSeaQrk == struck_sea_qrk );
   return equal;
 }

void Target::Copy ( const Target & t )

Definition at line 116 of file Target.cxx.

 {
   fTgtPDG = tgt.fTgtPDG;
 
   if( pdg::IsIon(fTgtPDG) ) {
 
      fZ         = tgt.fZ; // copy A,Z
      fA         = tgt.fA;
      fHitNucPDG = tgt.fHitNucPDG; // struck nucleon PDG
      fHitQrkPDG = tgt.fHitQrkPDG; // struck quark PDG
      fHitSeaQrk = tgt.fHitSeaQrk; // struck quark is from sea?
 
      //// valgrind warns about this ... try something else
      // (*fHitNucP4) = (*tgt.fHitNucP4);
      const TLorentzVector& p4 = *(tgt.fHitNucP4);
      //  *fHitNucP4 = p4; // nope
      //// this works for valgrind
      fHitNucP4->SetX(p4.X());
      fHitNucP4->SetY(p4.Y());
      fHitNucP4->SetZ(p4.Z());
      fHitNucP4->SetT(p4.T());
 
      fHitNucRad = tgt.fHitNucRad;
 
      // look-up the nucleus in the isotopes chart
      this->ForceNucleusValidity();
 
      // make sure the hit nucleus constituent object is either
      // a nucleon (p or n) or a di-nucleon cluster (p+p, p+n, n+n)
      this->ForceHitNucValidity();
   }
 }

void Target::ForceHitNucOnMassShell ( void )

Definition at line 200 of file Target.cxx.

 {
   if(this->HitNucIsSet()) {
      double m = this->HitNucMass();
      double p = this->HitNucP4Ptr()->P();
      double e = TMath::Sqrt(p*p+m*m);
      this->HitNucP4Ptr()->SetE(e);
   }
 }

bool Target::ForceHitNucValidity ( void )

private

Definition at line 350 of file Target.cxx.

 {
 // resets the struck nucleon pdg-code if it is found not to be a valid one
 
   bool valid =
       pdg::IsNucleon(fHitNucPDG)          ||
       pdg::Is2NucleonCluster (fHitNucPDG) ||
       (fHitNucPDG==0); /* not set */
 
   return valid;
 }

void Target::ForceNucleusValidity ( void )

private

Definition at line 362 of file Target.cxx.

 {
 // resets the target pdg-code if it is found not to be a valid one
 
   if( ! this->IsValidNucleus() ) {
     LOG("Target", pWARN) << "Invalid target -- Reseting to Z = 0, A = 0";
     fZ = 0;
     fA = 0;
   }
 }

bool Target::HitNucIsSet ( void ) const

Definition at line 283 of file Target.cxx.

 {
   bool ok =
      pdg::IsNucleon(fHitNucPDG)          ||
      pdg::Is2NucleonCluster (fHitNucPDG);
 
   return ok;
 }

double Target::HitNucMass ( void ) const

Definition at line 233 of file Target.cxx.

 {
   if(!fHitNucPDG) {
     LOG("Target", pWARN) << "Returning struck nucleon mass = 0";
     return 0;
   }
   return PDGLibrary::Instance()->Find(fHitNucPDG)->Mass();
 }

const TLorentzVector& genie::Target::HitNucP4 ( void ) const

inline

Definition at line 91 of file Target.h.

91 { return *this->HitNucP4Ptr(); }

genie::Target::HitNucP4Ptr

TLorentzVector * HitNucP4Ptr(void) const

Definition: Target.cxx:247

TLorentzVector * Target::HitNucP4Ptr ( void ) const

Definition at line 247 of file Target.cxx.

 {
   if(!fHitNucP4) {
     LOG("Target", pWARN) << "Returning NULL struck nucleon 4-momentum";
     return 0;
   }
 
   return fHitNucP4;
 }

int Target::HitNucPdg ( void ) const

Definition at line 304 of file Target.cxx.

 {
   return fHitNucPDG;
 }

double genie::Target::HitNucPosition ( void ) const

inline

Definition at line 89 of file Target.h.

89 { return fHitNucRad; }

genie::Target::fHitNucRad

double fHitNucRad

hit nucleon position

Definition: Target.h:124

bool Target::HitQrkIsSet ( void ) const

Definition at line 292 of file Target.cxx.

 {
   return (
      pdg::IsQuark(fHitQrkPDG) || pdg::IsAntiQuark(fHitQrkPDG)
   );
 }

int Target::HitQrkPdg ( void ) const

Definition at line 242 of file Target.cxx.

 {
   return fHitQrkPDG;
 }

bool Target::HitSeaQrk ( void ) const

Definition at line 299 of file Target.cxx.

 {
   return fHitSeaQrk;
 }

void Target::Init ( void )

private

Definition at line 99 of file Target.cxx.

 {
   fZ         = 0;
   fA         = 0;
   fTgtPDG    = 0;
   fHitNucPDG = 0;
   fHitQrkPDG = 0;
   fHitSeaQrk = false;
   fHitNucP4  = new TLorentzVector(0,0,0,kNucleonMass);
   fHitNucRad = 0.;
 }

bool Target::IsEvenEven ( void ) const

Definition at line 322 of file Target.cxx.

 {
   if( this->IsNucleus() ) {
     int NN = this->N();
     int ZZ = this->Z();
     if( NN % 2 == 0 && ZZ % 2 == 0 ) return true;
   }
   return false;
 }

bool Target::IsEvenOdd ( void ) const

Definition at line 332 of file Target.cxx.

 {
   if( this->IsNucleus() ) {
       if(! this->IsEvenEven() && ! this->IsOddOdd() ) return true;
   }
   return false;
 }

bool Target::IsFreeNucleon ( void ) const

Definition at line 257 of file Target.cxx.

 {
   return (fA == 1 && (fZ == 0 || fZ == 1));
 }

bool Target::IsNeutron ( void ) const

Definition at line 267 of file Target.cxx.

 {
   return (fA == 1 && fZ == 0);
 }

bool Target::IsNucleus ( void ) const

Definition at line 272 of file Target.cxx.

 {
   return (fA > 1); // IsValidNucleus() was ensured when A,Z were set
 }

bool Target::IsOddOdd ( void ) const

Definition at line 340 of file Target.cxx.

 {
   if( this->IsNucleus() ) {
     int NN = this->N();
     int ZZ = this->Z();
     if( NN % 2 == 1 && ZZ % 2 == 1 ) return true;
   }
   return false;
 }

bool Target::IsParticle ( void ) const

Definition at line 277 of file Target.cxx.

 {
   TParticlePDG * p = PDGLibrary::Instance()->Find(fTgtPDG);
   return (p && fA==0 && fZ==0);
 }

bool Target::IsProton ( void ) const

Definition at line 262 of file Target.cxx.

 {
   return (fA == 1 && fZ == 1);
 }

bool Target::IsValidNucleus ( void ) const

Definition at line 309 of file Target.cxx.

 {
   //-- it is valid if it is a free nucleon...
   if(this->IsFreeNucleon()) return true;
 
   //-- ... or a nucleus that can be found in the MINOS ion PDG extensions
   int pdg_code = pdg::IonPdgCode(fA, fZ);
   TParticlePDG * p = PDGLibrary::Instance()->Find(pdg_code);
   if(p) return true;
 
   return false;
 }

double Target::Mass ( void ) const

Definition at line 224 of file Target.cxx.

 {
 // Shortcut for commonly used code for extracting the nucleus mass from PDG
 //
   TParticlePDG * p = PDGLibrary::Instance()->Find(fTgtPDG);
   if(p) return p->Mass(); // in GeV
   return 0.;
 }

int genie::Target::N ( void ) const

inline

Definition at line 69 of file Target.h.

69 { return fA-fZ; }

genie::Target::fZ

int fZ

nuclear target Z

Definition: Target.h:117

genie::Target::fA

int fA

nuclear target A

Definition: Target.h:118

Target & Target::operator= ( const Target & t )

copy

Definition at line 441 of file Target.cxx.

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

bool Target::operator== ( const Target & t ) const

equal?

Definition at line 436 of file Target.cxx.

 {
   return this->Compare(target);
 }

int genie::Target::Pdg ( void ) const

inline

Definition at line 71 of file Target.h.

71 { return fTgtPDG; }

genie::Target::fTgtPDG

int fTgtPDG

nuclear target PDG code

Definition: Target.h:119

void Target::Print ( ostream & stream ) const

Definition at line 399 of file Target.cxx.

 {
   stream << " target PDG code = " << fTgtPDG << endl;
 
   if( this->IsNucleus() || this->IsFreeNucleon() ) {
       stream << " Z = " << fZ << ", A = " << fA << endl;
   }
 
   if( this->HitNucIsSet() ) {
     TParticlePDG * p = PDGLibrary::Instance()->Find(fHitNucPDG);
     stream << " struck nucleon = " << p->GetName()
               << ", P4 = " << utils::print::P4AsString(fHitNucP4) << endl;
   }
 
   if( this->HitQrkIsSet() ) {
     TParticlePDG * q = PDGLibrary::Instance()->Find(fHitQrkPDG);
     stream << " struck quark = " << q->GetName()
            << " (from sea: "
            << utils::print::BoolAsYNString(this->HitSeaQrk())
            << ")";
   }
 }

void Target::Reset ( void )

Definition at line 93 of file Target.cxx.

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

void Target::SetHitNucP4 ( const TLorentzVector & p4 )

Definition at line 189 of file Target.cxx.

 {
   if(fHitNucP4) delete fHitNucP4;
   fHitNucP4 = new TLorentzVector(p4);
 }

void Target::SetHitNucPdg ( int pdgc )

Definition at line 171 of file Target.cxx.

 {
   fHitNucPDG = nucl_pdgc;
   bool is_valid = this->ForceHitNucValidity();  // p, n or a di-nucleon
 
   // If it is a valid struck nucleon pdg code, initialize its 4P:
   // at-rest + on-mass-shell
   if(is_valid) {
     double M = PDGLibrary::Instance()->Find(nucl_pdgc)->Mass();
     fHitNucP4->SetPxPyPzE(0,0,0,M);
   }
 }

void Target::SetHitNucPosition ( double r )

Definition at line 210 of file Target.cxx.

 {
   fHitNucRad = r;
 }

void Target::SetHitQrkPdg ( int pdgc )

Definition at line 184 of file Target.cxx.

 {
   if(pdg::IsQuark(pdgc) || pdg::IsAntiQuark(pdgc)) fHitQrkPDG = pdgc;
 }

void Target::SetHitSeaQrk ( bool tf )

Definition at line 195 of file Target.cxx.

 {
   fHitSeaQrk = tf;
 }

void Target::SetId ( int pdgc )

Definition at line 149 of file Target.cxx.

 {
   fTgtPDG = pdgc;
   if( pdg::IsIon(pdgc) ) {
      fZ = pdg::IonPdgCodeToZ(pdgc);
      fA = pdg::IonPdgCodeToA(pdgc);
   }
 
   this->ForceNucleusValidity(); // search at the isotopes chart
   //this->AutoSetHitNuc();      // struck nuc := tgt for free nucleon tgt
 }

void Target::SetId	(	int	Z,
		int	A
	)

Definition at line 161 of file Target.cxx.

 {
   fTgtPDG = pdg::IonPdgCode(AA,ZZ);
   fZ = ZZ;
   fA = AA;
 
   this->ForceNucleusValidity(); // search at the isotopes chart
   //this->AutoSetHitNuc();      // struck nuc := tgt for free nucleon tgt
 }

int genie::Target::Z ( void ) const

inline

Definition at line 68 of file Target.h.

68 { return fZ; }

genie::Target::fZ

int fZ

nuclear target Z

Definition: Target.h:117

Friends And Related Function Documentation

ostream& operator<<	(	ostream &	stream,
		const Target &	t
	)

friend

print

Member Data Documentation

int genie::Target::fA

private

nuclear target A

Definition at line 118 of file Target.h.

TLorentzVector* genie::Target::fHitNucP4

private

hit nucleon 4p

Definition at line 123 of file Target.h.

int genie::Target::fHitNucPDG

private

hit nucleon PDG code

Definition at line 120 of file Target.h.

double genie::Target::fHitNucRad

private

hit nucleon position

Definition at line 124 of file Target.h.

int genie::Target::fHitQrkPDG

private

hit quark PDG code

Definition at line 121 of file Target.h.

bool genie::Target::fHitSeaQrk

private

hit quark from sea?

Definition at line 122 of file Target.h.

int genie::Target::fTgtPDG

private

nuclear target PDG code

Definition at line 119 of file Target.h.

int genie::Target::fZ

private

nuclear target Z

Definition at line 117 of file Target.h.

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

Generator/src/Framework/Interaction/Target.h
Generator/src/Framework/Interaction/Target.cxx

Public Member Functions

Private Member Functions

Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation