G4HadLeadBias Class Reference

#include <G4HadLeadBias.hh>

Inheritance diagram for G4HadLeadBias:

Public Member Functions
virtual G4HadFinalState *	Bias (G4HadFinalState *result)
virtual	~G4HadLeadBias ()

Detailed Description

Definition at line 35 of file G4HadLeadBias.hh.

Constructor & Destructor Documentation

virtual G4HadLeadBias::~G4HadLeadBias ( )

inlinevirtual

Definition at line 39 of file G4HadLeadBias.hh.

{};

Member Function Documentation

G4HadFinalState * G4HadLeadBias::Bias ( G4HadFinalState *  result )

virtual

Implements G4VLeadingParticleBiasing.

Definition at line 34 of file G4HadLeadBias.cc.

References G4HadFinalState::AddSecondaries(), buffer, G4HadFinalState::ClearSecondaries(), G4UniformRand, G4Gamma::Gamma(), G4ParticleDefinition::GetBaryonNumber(), G4DynamicParticle::GetDefinition(), G4HadFinalState::GetEnergyChange(), G4DynamicParticle::GetKineticEnergy(), G4ParticleDefinition::GetLeptonNumber(), G4HadFinalState::GetNumberOfSecondaries(), G4HadSecondary::GetParticle(), G4HadFinalState::GetSecondary(), G4HadFinalState::GetStatusChange(), G4HadSecondary::GetWeight(), isAlive, and G4PionZero::PionZero().

  {
    // G4cerr << "bias enter"<<G4endl;
    G4int nMeson(0), nBaryon(0), npi0(0), ngamma(0), nLepton(0);
    G4int i(0);
    G4int maxE = -1;
    G4double emax = 0;
    if(result->GetStatusChange()==isAlive) 
    {
      emax = result->GetEnergyChange();
    }
    //G4cout << "max energy "<<G4endl;
    for(i=0;i<result->GetNumberOfSecondaries();i++)
    {
      if(result->GetSecondary(i)->GetParticle()->GetKineticEnergy()>emax)
      {
        maxE = i;
    emax = result->GetSecondary(i)->GetParticle()->GetKineticEnergy();
      }
    }
    //G4cout <<"loop1"<<G4endl;
    for(i=0; i<result->GetNumberOfSecondaries(); i++)
    {
      const G4DynamicParticle* aSecTrack = result->GetSecondary(i)->GetParticle();
      if(i==maxE)
      {
      }
      else if(aSecTrack->GetDefinition()->GetBaryonNumber()!=0) 
      {
        nBaryon++;
      }
      else if(aSecTrack->GetDefinition()->GetLeptonNumber()!=0) 
      {
        nLepton++;
      }
      else if(aSecTrack->GetDefinition()==G4Gamma::Gamma())
      {
        ngamma++;
      }
      else if(aSecTrack->GetDefinition()==G4PionZero::PionZero())
      {
        npi0++;
      }
      else
      {
        nMeson++;
      }
    }
     //G4cout << "BiasDebug 1 = "<<result->GetNumberOfSecondaries()<<" "
     //       <<nMeson<<" "<< nBaryon<<" "<< npi0<<" "<< ngamma<<" "<< nLepton<<G4endl;
    G4double mesonWeight = nMeson;
    G4double baryonWeight = nBaryon;
    G4double gammaWeight = ngamma;
    G4double npi0Weight = npi0;
    G4double leptonWeight = nLepton;
    G4int randomMeson = static_cast<G4int>((nMeson+1)*G4UniformRand());
    G4int randomBaryon = static_cast<G4int>((nBaryon+1)*G4UniformRand());
    G4int randomGamma = static_cast<G4int>((ngamma+1)*G4UniformRand());
    G4int randomPi0 = static_cast<G4int>((npi0+1)*G4UniformRand());
    G4int randomLepton = static_cast<G4int>((nLepton+1)*G4UniformRand());
    
    std::vector<G4HadSecondary> buffer;
    G4int cMeson(0), cBaryon(0), cpi0(0), cgamma(0), cLepton(0);
    for(i=0; i<result->GetNumberOfSecondaries(); i++)
    {
      G4bool aCatch = false;
      G4double weight = 1;
      const G4HadSecondary* aSecTrack = result->GetSecondary(i);
      G4ParticleDefinition* aSecDef = aSecTrack->GetParticle()->GetDefinition();
      if(i==maxE)
      {
        aCatch = true;
    weight = 1;
      }
      else if(aSecDef->GetBaryonNumber()!=0) 
      {
    if(++cBaryon==randomBaryon) 
    {
      aCatch = true;
      weight = baryonWeight;
    }
      }
      else if(aSecDef->GetLeptonNumber()!=0) 
      {
        if(++cLepton==randomLepton) 
    {
      aCatch = true;
      weight = leptonWeight;
    }
      }
      else if(aSecDef==G4Gamma::Gamma())
      {
        if(++cgamma==randomGamma) 
    {
      aCatch = true;
      weight = gammaWeight;
    }
      }
      else if(aSecDef==G4PionZero::PionZero())
      {
        if(++cpi0==randomPi0) 
    {
      aCatch = true;
      weight = npi0Weight;
    }
      }
      else
      {
        if(++cMeson==randomMeson) 
    {
      aCatch = true;
      weight = mesonWeight;
    }
      }
      if(aCatch)
      {
    buffer.push_back(*aSecTrack);
    buffer.back().SetWeight(aSecTrack->GetWeight()*weight);
      }
      else
      {
        delete aSecTrack;
      }
    }
    result->ClearSecondaries();
    // G4cerr << "pre"<<G4endl;
    result->AddSecondaries(buffer);
     // G4cerr << "bias exit"<<G4endl;
    
    return result;
  }

---

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

• G4HadLeadBias.hh
• G4HadLeadBias.cc