Dear GEANT4 community,
I need to simulate C17 to shoot C target, so that C17 is in an excited state, then inelastic scattering occurs and fragmentation occurs, and then relevant information of fragmented fragments is collected. Since the probability of inelastic scattering is very low, I think only inelastic scattering will occur during the whole run.
I am not a physics major, so it is painful to write. Below is PhysicsList which I try to write. There may be a lot of problems. I hope you can help me correct them.
```
PhysicsList::PhysicsList()
: G4VModularPhysicsList() {
G4cout<<"00000000000000000000"<<G4endl;
// default cut value (0.1 mm)
defaultCutValue = 0.1* mm;
}
//…oooOO0OOooo…oooOO0OOooo…oooOO0OOooo…oooOO0OOooo…
PhysicsList::~PhysicsList()
{
}
//…oooOO0OOooo…oooOO0OOooo…oooOO0OOooo…oooOO0OOooo…
void PhysicsList::ConstructParticle()
{
G4IonConstructor IonConstructor;
IonConstructor.ConstructParticle();
G4BosonConstructor pBosonConstructor;
pBosonConstructor.ConstructParticle();
G4LeptonConstructor pLeptonConstructor;
pLeptonConstructor.ConstructParticle();
G4MesonConstructor pMesonConstructor;
pMesonConstructor.ConstructParticle();
G4BaryonConstructor BaryonConstructor;
BaryonConstructor.ConstructParticle();
G4ShortLivedConstructor pShortLivedConstructor;
pShortLivedConstructor.ConstructParticle();
}
//…oooOO0OOooo…oooOO0OOooo…oooOO0OOooo…oooOO0OOooo…
void PhysicsList::ConstructProcess()//…oooOO0OOooo…oooOO0OOooo…oooOO0OOooo…oooOO0OOooo…()
{
AddTransportation();
ConstructHadronInelasticProcess();
}
void PhysicsList::ConstructHadronInelasticProcess()
{
// Get pointer to C17
G4IonTable* ionTable = G4IonTable::GetIonTable();
G4ParticleDefinition* particle = ionTable->GetIon(6, 17);
G4HadronInelasticProcess* ionInelasticProcess = new G4HadronInelasticProcess("ionInelastic",G4GenericIon::GenericIon());
G4QMDReaction* qmd = new G4QMDReaction();
qmd->SetMinEnergy(0.*MeV);
qmd->SetMaxEnergy(10.*GeV);
G4IonsShenCrossSection* shenXS = new G4IonsShenCrossSection;
ionInelasticProcess->RegisterMe(qmd);
ionInelasticProcess->AddDataSet(shenXS);
G4ProcessManager* pmanager = G4GenericIon::GenericIon()->GetProcessManager();
pmanager->AddDiscreteProcess(ionInelasticProcess);
}
I used G4HadronInelasticProcess to achieve inelastic scattering, but after running the program almost always collected C17 behind the target, so it should prove that most of the C17 scatteredelastic. But I can't find the reason. I hope someone in the Geant4 community can help me find out what's going on!
Thank you in advance for your help and insights.
Best wishes,
Huiqin