G4CascadeHistory Class Reference

#include <G4CascadeHistory.hh>

Data Structures
struct	HistoryEntry

Public Member Functions
G4int	AddEntry (G4CascadParticle &cpart)
G4int	AddVertex (G4CascadParticle &cpart, std::vector< G4CascadParticle > &daug)
void	Clear ()
void	DropEntry (const G4CascadParticle &cpart)
	G4CascadeHistory (G4int verbose=0)
void	Print (std::ostream &os) const
void	setVerboseLevel (G4int verbose=0)
	~G4CascadeHistory ()

Protected Member Functions
void	AssignHistoryID (G4CascadParticle &cpart)
void	FillDaughters (G4int iEntry, std::vector< G4CascadParticle > &daug)
const char *	GuessTarget (const HistoryEntry &entry) const
void	PrintEntry (std::ostream &os, G4int iEntry) const
G4bool	PrintingDone (G4int iEntry) const
void	PrintParticle (std::ostream &os, const G4CascadParticle &cpart) const
G4int	size () const

Private Member Functions
	G4CascadeHistory (const G4CascadeHistory &rhs)
G4CascadeHistory &	operator= (const G4CascadeHistory &rhs)

Private Attributes
std::set< G4int >	entryPrinted
std::vector< HistoryEntry >	theHistory
G4int	verboseLevel

Detailed Description

Definition at line 40 of file G4CascadeHistory.hh.

Constructor & Destructor Documentation

G4CascadeHistory() [1/2]

G4CascadeHistory::G4CascadeHistory ( G4int  verbose = 0 )

inline

Definition at line 42 of file G4CascadeHistory.hh.

: verboseLevel(verbose) {;}

~G4CascadeHistory()

G4CascadeHistory::~G4CascadeHistory ( )

inline

Definition at line 43 of file G4CascadeHistory.hh.

{;}     // *** Do not want subclasses ***

G4CascadeHistory() [2/2]

G4CascadeHistory::G4CascadeHistory ( const G4CascadeHistory &  rhs )

private

Member Function Documentation

AddEntry()

G4int G4CascadeHistory::AddEntry

(

G4CascadParticle &

cpart

)

Definition at line 105 of file G4CascadeHistory.cc.

                                                        {
  AssignHistoryID(cpart);       // Make sure particle has index
 
  G4int id = cpart.getHistoryId();
  if (id < size()) {
    if (verboseLevel>2) 
      G4cout << " AddEntry updating " << id << " " << &theHistory[id] << G4endl;
    theHistory[id].cpart = cpart;       // Copies kinematics
  } else {
    theHistory.push_back(HistoryEntry(cpart));
    if (verboseLevel>2) 
      G4cout << " AddEntry creating " << id << " " << &theHistory.back() << G4endl;
  }
 
  if (verboseLevel>3) G4cout << theHistory[id].cpart << G4endl; // Sanity check
 
  return id;
}

References AssignHistoryID(), G4cout, G4endl, G4CascadParticle::getHistoryId(), size(), theHistory, and verboseLevel.

Referenced by AddVertex(), FillDaughters(), and G4IntraNucleiCascader::generateCascade().

AddVertex()

G4int G4CascadeHistory::AddVertex	(	G4CascadParticle &	cpart,
		std::vector< G4CascadParticle > &	daug
	)

Definition at line 58 of file G4CascadeHistory.cc.

                                                     {
  if (verboseLevel>1) G4cout << " >>> G4CascadeHistory::AddVertex" << G4endl;
 
  // Create new entry for vertex or update particle kinematics
  G4int id = AddEntry(cpart);
  FillDaughters(id, daug);
 
  if (verboseLevel>3) {
    G4cout << " entry " << id << " " << &theHistory[id] << " got "
       << theHistory[id].n << " daughters:";
    for (G4int i=0; i<theHistory[id].n; i++) {
      G4cout << " " << theHistory[id].dId[i];
    }
    G4cout << G4endl;
  }
 
  return id;
}

References AddEntry(), FillDaughters(), G4cout, G4endl, theHistory, and verboseLevel.

Referenced by G4IntraNucleiCascader::generateCascade().

AssignHistoryID()

void G4CascadeHistory::AssignHistoryID ( G4CascadParticle &  cpart )

protected

Definition at line 136 of file G4CascadeHistory.cc.

                                                              {
  if (cpart.getHistoryId() >= 0) return;        // ID already assigned
 
  if (verboseLevel>2) {
    G4cout << " >>> G4CascadeHistory::NewHistoryID assigning ID "
       << size() << G4endl;
  }
 
  cpart.setHistoryId(size());
}

References G4cout, G4endl, G4CascadParticle::getHistoryId(), G4CascadParticle::setHistoryId(), size(), and verboseLevel.

Referenced by AddEntry().

Clear()

void G4CascadeHistory::Clear

(

)

Definition at line 49 of file G4CascadeHistory.cc.

                             {
  if (verboseLevel>1) G4cout << " >>> G4CascadeHistory::Clear" << G4endl;
  theHistory.clear();
  entryPrinted.clear();
}

References entryPrinted, G4cout, G4endl, theHistory, and verboseLevel.

Referenced by G4IntraNucleiCascader::newCascade().

DropEntry()

void G4CascadeHistory::DropEntry

(

const G4CascadParticle &

cpart

)

Definition at line 126 of file G4CascadeHistory.cc.

                                                              {
  if (verboseLevel>1) G4cout << " >>> G4CascadeHistory::DropEntry" << G4endl;
 
  
  G4int id = cpart.getHistoryId();  // Particle must appear in history
  if (id>=0) theHistory[id].n = -1; // Special flag for absorbed particle
}

References G4cout, G4endl, G4CascadParticle::getHistoryId(), theHistory, and verboseLevel.

Referenced by G4IntraNucleiCascader::processTrappedParticle().

FillDaughters()

void G4CascadeHistory::FillDaughters ( G4int  iEntry, std::vector< G4CascadParticle > &  daug  )

protected

Definition at line 79 of file G4CascadeHistory.cc.

                                                                               {
  G4int nDaug = (G4int)daug.size();
 
  if (verboseLevel>1) 
    G4cout << " >>> G4CascadeHistory::FillDaughters " << iEntry << G4endl;
 
  // NOTE:  Cannot use reference to element, as push_back can invalidate refs!
  theHistory[iEntry].clear();
 
  theHistory[iEntry].n = nDaug;
  for (G4int i=0; i<nDaug; i++) {
    G4int id = AddEntry(daug[i]);
    theHistory[iEntry].dId[i] = id;
  }
 
  if (verboseLevel>3) {
    G4cout << " got " << theHistory[iEntry].n << " daughters:";
    for (G4int i=0; i<theHistory[iEntry].n; i++) {
      G4cout << " " << theHistory[iEntry].dId[i];
    }
    G4cout << G4endl;
  }
}

References AddEntry(), G4cout, G4endl, theHistory, and verboseLevel.

Referenced by AddVertex().

GuessTarget()

const char * G4CascadeHistory::GuessTarget ( const HistoryEntry &  entry ) const

protected

Definition at line 203 of file G4CascadeHistory.cc.

                                                             {
  if (verboseLevel>2) G4cout << " >>> G4CascadeHistory::GuessTarget" << G4endl;
 
  if (entry.n < 0) return "-O-";    // Exciton or trapped-decay
  if (entry.n == 0) return "***";   // Outgoing (final state) particle
 
  const G4CascadParticle& cpart = entry.cpart;  // For convenience
  if (verboseLevel>3) G4cout << "cpart: " << cpart;
 
  // Compute baryon number and charge from daughters minus projectile
  G4int targetB = -cpart.getParticle().baryon();
  G4int targetQ = (G4int)-cpart.getParticle().getCharge();
 
  for (G4int i=0; i<entry.n; i++) {
    const G4CascadParticle& cdaug = theHistory[entry.dId[i]].cpart;
    if (verboseLevel>3) 
      G4cout << "cdaug " << i << " ID " << entry.dId[i] << ": " << cdaug;
 
    targetB += cdaug.getParticle().baryon();
    targetQ += (G4int)cdaug.getParticle().getCharge();
  }
 
  // Target possibilities are proton, neutron or dibaryon (pp, nn, pn)
  if (targetB==1 && targetQ==0) return "n";
  if (targetB==1 && targetQ==1) return "p";
  if (targetB==2 && targetQ==0) return "nn";
  if (targetB==2 && targetQ==1) return "pn";
  if (targetB==2 && targetQ==2) return "pp";
 
  if (verboseLevel>2) {
    G4cout << " ERROR identifying target: deltaB " << targetB
       << " deltaQ " << targetQ << " from\n" << cpart << " to" << G4endl;
    for (G4int j=0; j<entry.n; j++) {
      G4cout << theHistory[entry.dId[j]].cpart;
    }
  }
 
  return "BAD TARGET";      // Should not get here if EPCollider worked right
}

References G4InuclElementaryParticle::baryon(), G4CascadeHistory::HistoryEntry::cpart, G4CascadeHistory::HistoryEntry::dId, G4cout, G4endl, G4InuclParticle::getCharge(), G4CascadParticle::getParticle(), G4CascadeHistory::HistoryEntry::n, theHistory, and verboseLevel.

Referenced by PrintEntry().

operator=()

G4CascadeHistory & G4CascadeHistory::operator= ( const G4CascadeHistory &  rhs )

private

Print()

void G4CascadeHistory::Print

(

std::ostream &

os

)

const

Definition at line 155 of file G4CascadeHistory.cc.

                                                 {
  if (verboseLevel) os << " >>> G4CascadeHistory::Print" << std::endl;
 
  os << " Cascade structure: vertices, (-O-) exciton, (***) outgoing"
     << std::endl;
 
  for (G4int i=0; i<size(); i++) {
    if (!PrintingDone(i)) PrintEntry(os, i);    
  }
}

References PrintEntry(), PrintingDone(), size(), and verboseLevel.

Referenced by G4IntraNucleiCascader::collide(), and G4IntraNucleiCascader::rescatter().

PrintEntry()

void G4CascadeHistory::PrintEntry ( std::ostream &  os, G4int  iEntry  ) const

protected

Definition at line 168 of file G4CascadeHistory.cc.

                                                                    {
  if (iEntry >= size()) return;         // Skip nonexistent entry
  if (PrintingDone(iEntry)) return;     // Skip entry already reported
 
  entryPrinted.insert(iEntry);
 
  const HistoryEntry& entry = theHistory[iEntry];   // For convenience
  const G4CascadParticle& cpart = entry.cpart;
 
  G4int indent = cpart.getGeneration()*2;
 
  // Index and indentation of cascade vertex
  std::ios::fmtflags osFlags = os.flags();
  os.setf(std::ios::left);  // Pushes all blanks to right end of output
  os << "#" << std::setw(3+indent) << iEntry;
  os.flags(osFlags);
 
  os << cpart.getParticle().getDefinition()->GetParticleName()
     << " p " << cpart.getMomentum() << " (cosTh "
     << cpart.getMomentum().vect().unit().z() << ")"
     << " @ " << cpart.getPosition()
     << " zone " << cpart.getCurrentZone();
 
  // Flag as final-state particle or report daughters iteratively
  os << " (" << GuessTarget(entry) << ")";
  if (entry.n > 0) {
    os << " -> N=" << entry.n << std::endl;
    for (G4int i=0; i<entry.n; i++) {
      PrintEntry(os, entry.dId[i]);
    }
  } else os << std::endl;
}

References G4CascadeHistory::HistoryEntry::cpart, G4CascadeHistory::HistoryEntry::dId, entryPrinted, G4CascadParticle::getCurrentZone(), G4InuclParticle::getDefinition(), G4CascadParticle::getGeneration(), G4CascadParticle::getMomentum(), G4CascadParticle::getParticle(), G4ParticleDefinition::GetParticleName(), G4CascadParticle::getPosition(), GuessTarget(), G4CascadeHistory::HistoryEntry::n, PrintEntry(), PrintingDone(), size(), theHistory, CLHEP::Hep3Vector::unit(), CLHEP::HepLorentzVector::vect(), and CLHEP::Hep3Vector::z().

Referenced by Print(), and PrintEntry().

PrintingDone()

G4bool G4CascadeHistory::PrintingDone ( G4int  iEntry ) const

inlineprotected

Definition at line 82 of file G4CascadeHistory.hh.

                                          {
    return (entryPrinted.find(iEntry) != entryPrinted.end());
  }

References entryPrinted.

Referenced by Print(), and PrintEntry().

PrintParticle()

void G4CascadeHistory::PrintParticle ( std::ostream &  os, const G4CascadParticle &  cpart  ) const

protected

setVerboseLevel()

void G4CascadeHistory::setVerboseLevel ( G4int  verbose = 0 )

inline

Definition at line 45 of file G4CascadeHistory.hh.

{ verboseLevel = verbose; }

References verboseLevel.

Referenced by G4IntraNucleiCascader::setVerboseLevel().

size()

G4int G4CascadeHistory::size ( ) const

inlineprotected

Definition at line 89 of file G4CascadeHistory.hh.

{ return (G4int)theHistory.size(); }

References theHistory.

Referenced by AddEntry(), AssignHistoryID(), Print(), and PrintEntry().

Field Documentation

entryPrinted

std::set<G4int> G4CascadeHistory::entryPrinted

mutableprivate

Definition at line 95 of file G4CascadeHistory.hh.

Referenced by Clear(), PrintEntry(), and PrintingDone().

theHistory

std::vector<HistoryEntry> G4CascadeHistory::theHistory

private

Definition at line 94 of file G4CascadeHistory.hh.

Referenced by AddEntry(), AddVertex(), Clear(), DropEntry(), FillDaughters(), GuessTarget(), PrintEntry(), and size().

verboseLevel

G4int G4CascadeHistory::verboseLevel

private

Definition at line 92 of file G4CascadeHistory.hh.

Referenced by AddEntry(), AddVertex(), AssignHistoryID(), Clear(), DropEntry(), FillDaughters(), GuessTarget(), Print(), and setVerboseLevel().

---

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

• source/processes/hadronic/models/cascade/cascade/include/G4CascadeHistory.hh
• source/processes/hadronic/models/cascade/cascade/src/G4CascadeHistory.cc