#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 36 of file MuNuclearSplittingProcess.cxx.

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