// // ******************************************************************** // * License and Disclaimer * // * * // * The Geant4 software is copyright of the Copyright Holders of * // * the Geant4 Collaboration. It is provided under the terms and * // * conditions of the Geant4 Software License, included in the file * // * LICENSE and available at http://cern.ch/geant4/license . These * // * include a list of copyright holders. * // * * // * Neither the authors of this software system, nor their employing * // * institutes,nor the agencies providing financial support for this * // * work make any representation or warranty, express or implied, * // * regarding this software system or assume any liability for its * // * use. Please see the license in the file LICENSE and URL above * // * for the full disclaimer and the limitation of liability. * // * * // * This code implementation is the result of the scientific and * // * technical work of the GEANT4 collaboration. * // * By using, copying, modifying or distributing the software (or * // * any work based on the software) you agree to acknowledge its * // * use in resulting scientific publications, and indicate your * // * acceptance of all terms of the Geant4 Software license. * // ******************************************************************** // // $Id: G4EmStandardPhysics_option3.cc,v 1.21 2010/10/10 15:18:34 vnivanch Exp $ // GEANT4 tag $Name: phys-lists-V09-03-34 $ // //--------------------------------------------------------------------------- // // ClassName: G4EmStandardPhysics_option3 // // Author: V.Ivanchenko 13.03.2008 // // Modified: // 21.04.2008 V.Ivanchenko add long-lived D and B mesons; use spline // 28.05.2008 V.Ivanchenko linLossLimit=0.01 for ions 0.001 for others // //---------------------------------------------------------------------------- // #include "G4EmStandardPhysics_option3.hh" #include "G4ParticleDefinition.hh" #include "G4ProcessManager.hh" #include "G4LossTableManager.hh" #include "G4EmProcessOptions.hh" #include "G4ComptonScattering.hh" #include "G4GammaConversion.hh" #include "G4PhotoElectricEffect.hh" #include "G4RayleighScattering.hh" #include "G4eMultipleScattering.hh" #include "G4MuMultipleScattering.hh" #include "G4hMultipleScattering.hh" #include "G4MscStepLimitType.hh" #include "G4UrbanMscModel93.hh" #include "G4WentzelVIModel.hh" #include "G4CoulombScattering.hh" #include "G4eIonisation.hh" #include "G4eBremsstrahlung.hh" #include "G4eplusAnnihilation.hh" #include "G4MuIonisation.hh" #include "G4MuBremsstrahlung.hh" #include "G4MuPairProduction.hh" #include "G4hBremsstrahlung.hh" #include "G4hPairProduction.hh" #include "G4hIonisation.hh" #include "G4ionIonisation.hh" #include "G4IonParametrisedLossModel.hh" #include "G4NuclearStopping.hh" #include "G4Gamma.hh" #include "G4Electron.hh" #include "G4Positron.hh" #include "G4MuonPlus.hh" #include "G4MuonMinus.hh" #include "G4PionPlus.hh" #include "G4PionMinus.hh" #include "G4KaonPlus.hh" #include "G4KaonMinus.hh" #include "G4Proton.hh" #include "G4AntiProton.hh" #include "G4Deuteron.hh" #include "G4Triton.hh" #include "G4He3.hh" #include "G4Alpha.hh" #include "G4GenericIon.hh" //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... G4EmStandardPhysics_option3::G4EmStandardPhysics_option3(G4int ver) : G4VPhysicsConstructor("G4EmStandard_opt3"), verbose(ver) { G4LossTableManager::Instance(); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... G4EmStandardPhysics_option3::G4EmStandardPhysics_option3(G4int ver, const G4String&) : G4VPhysicsConstructor("G4EmStandard_opt3"), verbose(ver) { G4LossTableManager::Instance(); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... G4EmStandardPhysics_option3::~G4EmStandardPhysics_option3() {} //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... void G4EmStandardPhysics_option3::ConstructParticle() { // gamma G4Gamma::Gamma(); // leptons G4Electron::Electron(); G4Positron::Positron(); G4MuonPlus::MuonPlus(); G4MuonMinus::MuonMinus(); // mesons G4PionPlus::PionPlusDefinition(); G4PionMinus::PionMinusDefinition(); G4KaonPlus::KaonPlusDefinition(); G4KaonMinus::KaonMinusDefinition(); // barions G4Proton::Proton(); G4AntiProton::AntiProton(); // ions G4Deuteron::Deuteron(); G4Triton::Triton(); G4He3::He3(); G4Alpha::Alpha(); G4GenericIon::GenericIonDefinition(); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... void G4EmStandardPhysics_option3::ConstructProcess() { // Add standard EM Processes theParticleIterator->reset(); while( (*theParticleIterator)() ){ G4ParticleDefinition* particle = theParticleIterator->value(); G4ProcessManager* pmanager = particle->GetProcessManager(); G4String particleName = particle->GetParticleName(); if(verbose > 1) G4cout << "### " << GetPhysicsName() << " instantiates for " << particleName << G4endl; if (particleName == "gamma") { pmanager->AddDiscreteProcess(new G4PhotoElectricEffect); pmanager->AddDiscreteProcess(new G4ComptonScattering); pmanager->AddDiscreteProcess(new G4GammaConversion); pmanager->AddDiscreteProcess(new G4RayleighScattering); } else if (particleName == "e-") { G4eMultipleScattering* msc = new G4eMultipleScattering(); //msc->AddEmModel(0, new G4UrbanMscModel93()); msc->SetStepLimitType(fUseDistanceToBoundary); pmanager->AddProcess(msc, -1, 1, 1); G4eIonisation* eIoni = new G4eIonisation(); eIoni->SetStepFunction(0.2, 100*um); pmanager->AddProcess(eIoni, -1, 2, 2); pmanager->AddProcess(new G4eBremsstrahlung, -1,-3, 3); } else if (particleName == "e+") { G4eMultipleScattering* msc = new G4eMultipleScattering(); //msc->AddEmModel(0, new G4UrbanMscModel93()); msc->SetStepLimitType(fUseDistanceToBoundary); pmanager->AddProcess(msc, -1, 1, 1); G4eIonisation* eIoni = new G4eIonisation(); eIoni->SetStepFunction(0.2, 100*um); pmanager->AddProcess(eIoni, -1, 2, 2); pmanager->AddProcess(new G4eBremsstrahlung, -1,-3, 3); pmanager->AddProcess(new G4eplusAnnihilation,0,-1, 4); } else if (particleName == "mu+" || particleName == "mu-" ) { G4MuMultipleScattering* msc = new G4MuMultipleScattering(); msc->AddEmModel(0, new G4WentzelVIModel()); pmanager->AddProcess(msc, -1, 1, 1); G4MuIonisation* muIoni = new G4MuIonisation(); muIoni->SetStepFunction(0.2, 50*um); pmanager->AddProcess(muIoni, -1, 2, 2); pmanager->AddProcess(new G4MuBremsstrahlung, -1,-3, 3); pmanager->AddProcess(new G4MuPairProduction, -1,-4, 4); pmanager->AddDiscreteProcess(new G4CoulombScattering()); } else if (particleName == "alpha" || particleName == "He3") { pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1); G4ionIonisation* ionIoni = new G4ionIonisation(); ionIoni->SetStepFunction(0.1, 20*um); pmanager->AddProcess(ionIoni, -1, 2, 2); pmanager->AddProcess(new G4NuclearStopping(), -1, 3,-1); } else if (particleName == "GenericIon") { pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1); G4ionIonisation* ionIoni = new G4ionIonisation(); ionIoni->SetEmModel(new G4IonParametrisedLossModel()); ionIoni->SetStepFunction(0.1, 10*um); pmanager->AddProcess(ionIoni, -1, 2, 2); pmanager->AddProcess(new G4NuclearStopping(), -1, 3,-1); } else if (particleName == "pi+" || particleName == "pi-" || particleName == "kaon+" || particleName == "kaon-" || particleName == "proton" ) { pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1); G4hIonisation* hIoni = new G4hIonisation(); hIoni->SetStepFunction(0.2, 50*um); pmanager->AddProcess(hIoni, -1, 2, 2); pmanager->AddProcess(new G4hBremsstrahlung, -1,-3, 3); pmanager->AddProcess(new G4hPairProduction, -1,-4, 4); } else if (particleName == "B+" || particleName == "B-" || particleName == "D+" || particleName == "D-" || particleName == "Ds+" || particleName == "Ds-" || particleName == "anti_He3" || particleName == "anti_alpha" || particleName == "anti_deuteron" || particleName == "anti_lambda_c+" || particleName == "anti_omega-" || particleName == "anti_proton" || particleName == "anti_sigma_c+" || particleName == "anti_sigma_c++" || particleName == "anti_sigma+" || particleName == "anti_sigma-" || particleName == "anti_triton" || particleName == "anti_xi_c+" || particleName == "anti_xi-" || particleName == "deuteron" || particleName == "lambda_c+" || particleName == "omega-" || particleName == "sigma_c+" || particleName == "sigma_c++" || particleName == "sigma+" || particleName == "sigma-" || particleName == "tau+" || particleName == "tau-" || particleName == "triton" || particleName == "xi_c+" || particleName == "xi-" ) { pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1); pmanager->AddProcess(new G4hIonisation, -1, 2, 2); } } // Em options // G4EmProcessOptions opt; opt.SetVerbose(verbose); // Multiple Coulomb scattering // //opt.SetMscStepLimitation(fUseDistanceToBoundary); //opt.SetMscRangeFactor(0.02); // Physics tables // opt.SetMinEnergy(100*eV); opt.SetMaxEnergy(10*TeV); opt.SetDEDXBinning(220); opt.SetLambdaBinning(220); //opt.SetSplineFlag(true); opt.SetPolarAngleLimit(0.2); // Ionization // //opt.SetSubCutoff(true); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......