G4CascadeHistory Class Reference

#include <G4CascadeHistory.hh>

Classes
struct	HistoryEntry

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

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

Detailed Description

Definition at line 40 of file G4CascadeHistory.hh.

Constructor & Destructor Documentation

G4CascadeHistory()

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 ***

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

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

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

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

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
}

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

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
}

Referenced by PrintEntry().

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

Referenced by G4IntraNucleiCascader::collide(), operator<<(), 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;
}

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

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

Referenced by G4IntraNucleiCascader::setVerboseLevel().

size()

G4int G4CascadeHistory::size ( ) const

inlineprotected

Definition at line 89 of file G4CascadeHistory.hh.

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

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

---

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

• G4CascadeHistory.hh
• G4CascadeHistory.cc