G4CascadeHistory Class Reference

#include <G4CascadeHistory.hh>

Data Structures
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 41 of file G4CascadeHistory.hh.

Constructor & Destructor Documentation

G4CascadeHistory::G4CascadeHistory ( G4int  verbose = 0 )

inline

Definition at line 43 of file G4CascadeHistory.hh.

: verboseLevel(verbose) {;}

G4CascadeHistory::~G4CascadeHistory ( )

inline

Definition at line 44 of file G4CascadeHistory.hh.

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

Member Function Documentation

G4int G4CascadeHistory::AddEntry

(

G4CascadParticle &

cpart

)

Definition at line 106 of file G4CascadeHistory.cc.

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

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

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

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

Definition at line 59 of file G4CascadeHistory.cc.

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

Referenced by G4IntraNucleiCascader::generateCascade().

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

void G4CascadeHistory::AssignHistoryID ( G4CascadParticle &  cpart )

protected

Definition at line 137 of file G4CascadeHistory.cc.

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

Referenced by AddEntry().

                                                              {
  if (cpart.getHistoryId() >= 0) return;        // ID already assigned

  if (verboseLevel>2) {
    G4cout << " >>> G4CascadeHistory::NewHistoryID assigning ID "
       << size() << G4endl;
  }

  cpart.setHistoryId(size());
}

void G4CascadeHistory::Clear

(

)

Definition at line 50 of file G4CascadeHistory.cc.

References G4cout, and G4endl.

Referenced by G4IntraNucleiCascader::newCascade().

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

void G4CascadeHistory::DropEntry

(

const G4CascadParticle &

cpart

)

Definition at line 127 of file G4CascadeHistory.cc.

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

Referenced by G4IntraNucleiCascader::processTrappedParticle().

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

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

protected

Definition at line 80 of file G4CascadeHistory.cc.

References AddEntry(), G4cout, and G4endl.

Referenced by AddVertex().

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

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

protected

Definition at line 205 of file G4CascadeHistory.cc.

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

Referenced by PrintEntry().

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

void G4CascadeHistory::Print

(

std::ostream &

os

)

const

Definition at line 156 of file G4CascadeHistory.cc.

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

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

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

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

protected

Definition at line 169 of file G4CascadeHistory.cc.

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

Referenced by Print().

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

G4bool G4CascadeHistory::PrintingDone ( G4int  iEntry ) const

inlineprotected

Definition at line 83 of file G4CascadeHistory.hh.

Referenced by Print(), and PrintEntry().

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

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

protected

void G4CascadeHistory::setVerboseLevel ( G4int  verbose = 0 )

inline

Definition at line 46 of file G4CascadeHistory.hh.

Referenced by G4IntraNucleiCascader::setVerboseLevel().

{ verboseLevel = verbose; }

G4int G4CascadeHistory::size ( ) const

inlineprotected

Definition at line 90 of file G4CascadeHistory.hh.

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

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

---

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

• G4CascadeHistory.hh
• G4CascadeHistory.cc