#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 20 of file MuNuclearSplittingProcess.h.

Constructor & Destructor Documentation

larg4::MuNuclearSplittingProcess::MuNuclearSplittingProcess ( )

inline

Definition at line 23 of file MuNuclearSplittingProcess.h.

23 {};

larg4::MuNuclearSplittingProcess::~MuNuclearSplittingProcess ( )

inline

Definition at line 24 of file MuNuclearSplittingProcess.h.

24 {};

Member Function Documentation

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

private

Definition at line 35 of file MuNuclearSplittingProcess.cxx.

References util::abs(), 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;
   }