larg4::MuNuclearSplittingProcess Class Reference

#include "MuNuclearSplittingProcess.h"

Inheritance diagram for larg4::MuNuclearSplittingProcess:

Public Member Functions
	MuNuclearSplittingProcess ()
	~MuNuclearSplittingProcess ()
void	SetNSplit (G4int nTrx)
void	SetIsActive (G4bool doIt)

Private Member Functions
G4VParticleChange *	PostStepDoIt (const G4Track &track, const G4Step &step)

Private Attributes
G4int	fNSplit
G4bool	fActive

Detailed Description

Definition at line 29 of file MuNuclearSplittingProcess.h.

Constructor & Destructor Documentation

larg4::MuNuclearSplittingProcess::MuNuclearSplittingProcess ( )

inline

Definition at line 32 of file MuNuclearSplittingProcess.h.

{};

larg4::MuNuclearSplittingProcess::~MuNuclearSplittingProcess ( )

inline

Definition at line 33 of file MuNuclearSplittingProcess.h.

{};

Member Function Documentation

G4VParticleChange * larg4::MuNuclearSplittingProcess::PostStepDoIt ( const G4Track &  track, const G4Step &  step  )

private

Definition at line 39 of file MuNuclearSplittingProcess.cxx.

References fNSplit, and weight.

{

  G4VParticleChange* particleChange = new G4VParticleChange();
  G4double weight = track.GetWeight()/fNSplit;
  std::vector<G4Track*> secondaries; // Secondary store

  // Loop over PostStepDoIt method to generate multiple secondaries.
  // The point is for each value of ii up to Nsplit we want to re-toss
  // all secondaries for the muNucl process. Each toss gives back numSec
  // secondaries, which will be a different sample of species/momenta each time.
  for (unsigned int ii=0; ii<(unsigned int)fNSplit; ii++) {
    particleChange = pRegProcess->PostStepDoIt(track, step);
    if (!particleChange)
      throw std::runtime_error("MuNuclearSplittingProcess::PostStepDoIt(): no particle change");
    G4int j(0);
    G4int numSec(particleChange->GetNumberOfSecondaries());
    for (j=0; j<numSec; j++)
      {
  G4Track* newSec = new G4Track(*(particleChange->GetSecondary(j)));
  G4String pdgstr = newSec->GetParticleDefinition()->GetParticleName();
  G4double ke = newSec->GetKineticEnergy()/CLHEP::GeV;
  G4int pdg = newSec->GetParticleDefinition()->GetPDGEncoding();
  if (abs(pdg)==310 ||abs(pdg)==311 || abs(pdg)==3122 || abs(pdg)==2112)
    {
      //      std::cout << "MuNuclSplProc: Creating " << pdgstr << " of Kinetic Energy " << ke << " GeV. numSec is " << j << std::endl;

      if (pdg==G4KaonZeroShort::KaonZeroShort()->GetPDGEncoding()&&G4UniformRand()<0.50)
        {
    pdg = G4KaonZeroLong::KaonZeroLong()->GetPDGEncoding();
    pdgstr = G4KaonZeroLong::KaonZeroLong()->GetParticleName();
    G4LorentzVector pK0L(newSec->GetMomentum(),TMath::Sqrt(TMath::Power(G4KaonZeroLong::KaonZeroLong()->GetPDGMass(),2)+TMath::Power(newSec->GetMomentum().mag(),2)));
    G4DynamicParticle *newK0L = new G4DynamicParticle(G4KaonZeroLong::KaonZeroLong(),pK0L);

    G4Track* newSecK0L = new G4Track(newK0L,track.GetGlobalTime(),track.GetPosition());
    secondaries.push_back(newSecK0L);
        }
      else
        {
    secondaries.push_back(newSec);
        }

      if (abs(pdg)==130 ||abs(pdg)==310 ||abs(pdg)==311 || abs(pdg)==3122)
        {
    // WrappedmuNuclear always produces K0s if it produces a K0 at all.
    // Let's make half of these K0Ls.
    std::cout << "MuNuclSplProc: Creating " << pdgstr << " of Kinetic Energy " << ke << " GeV. numSec is " << j << std::endl;
        }

    }
      }
  }

  particleChange->SetNumberOfSecondaries(secondaries.size());
  particleChange->SetSecondaryWeightByProcess(true);
  //int numSecAdd = secondaries.size();
  std::vector<G4Track*>::iterator iter = secondaries.begin(); // Add all secondaries
  while (iter != secondaries.end()) {
    G4Track* myTrack = *iter;
    G4String pdgstr = myTrack->GetParticleDefinition()->GetParticleName();
    //G4double ke = myTrack->GetKineticEnergy()/GeV;
    G4int pdg = myTrack->GetParticleDefinition()->GetPDGEncoding();
    if (abs(pdg)==130 ||abs(pdg)==310 ||abs(pdg)==311 || abs(pdg)==3122 || abs(pdg)==2112)
      //      std::cout << "MuNuclSplProc: Adding " << pdgstr << " of Kinetic Energy " << ke << " GeV." << std::endl;
    // Will ask for PdgCode and only Add K0s.
    myTrack->SetWeight(weight);
    particleChange->AddSecondary(myTrack);
    iter++;
  }
  return particleChange;
}

void larg4::MuNuclearSplittingProcess::SetIsActive ( G4bool  doIt )

inline

Definition at line 36 of file MuNuclearSplittingProcess.h.

References fActive, and fNSplit.

{fActive = doIt;}; 

void larg4::MuNuclearSplittingProcess::SetNSplit ( G4int  nTrx )

inline

Definition at line 35 of file MuNuclearSplittingProcess.h.

References fNSplit.

{fNSplit = nTrx;};

Member Data Documentation

G4bool larg4::MuNuclearSplittingProcess::fActive

private

Definition at line 41 of file MuNuclearSplittingProcess.h.

Referenced by SetIsActive().

G4int larg4::MuNuclearSplittingProcess::fNSplit

private

Definition at line 36 of file MuNuclearSplittingProcess.h.

Referenced by PostStepDoIt(), SetIsActive(), and SetNSplit().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larsim/v07_09_00/source/larsim/LArG4/MuNuclearSplittingProcess.h
• /cvmfs/larsoft.opensciencegrid.org/products/larsim/v07_09_00/source/larsim/LArG4/MuNuclearSplittingProcess.cxx