// // ******************************************************************** // * 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. * // ******************************************************************** // // // File name: RadmonPhysicsNuclear.cc // Creation date: Nov 2005 // Main author: Riccardo Capra // // Id: $Id: RadmonPhysicsNuclear.cc,v 1.3 2006/06/29 16:19:33 gunter Exp $ // Tag: $Name: $ // #include "RadmonPhysicsNuclear.hh" #include "G4Electron.hh" #include "G4Positron.hh" #include "G4Gamma.hh" #include "G4ProcessManager.hh" #include "G4ElectronNuclearProcess.hh" #include "G4PositronNuclearProcess.hh" #include "G4PhotoNuclearProcess.hh" // models #include "G4ElectroNuclearReaction.hh" #include "G4GammaNuclearReaction.hh" #include "G4TheoFSGenerator.hh" #include "G4GeneratorPrecompoundInterface.hh" #include "G4QGSMFragmentation.hh" #include "G4ExcitedStringDecay.hh" #include "G4GammaParticipants.hh" #include "G4QGSModel.hh" RadmonVSubPhysicsListWithLabel * RadmonPhysicsNuclear :: New(void) const { return new RadmonPhysicsNuclear; } void RadmonPhysicsNuclear :: ConstructParticle(void) { G4Positron::PositronDefinition(); G4Gamma::GammaDefinition(); G4Electron::ElectronDefinition(); } void RadmonPhysicsNuclear :: ConstructProcess(void) { G4ElectroNuclearReaction* electronReaction(new G4ElectroNuclearReaction); G4ProcessManager * manager(G4Electron::ElectronDefinition()->GetProcessManager()); G4ElectronNuclearProcess * electronNuclearProcess(new G4ElectronNuclearProcess); electronNuclearProcess->RegisterMe(electronReaction); manager->AddDiscreteProcess(electronNuclearProcess); manager=G4Positron::PositronDefinition()->GetProcessManager(); G4PositronNuclearProcess * positronNuclearProcess(new G4PositronNuclearProcess); positronNuclearProcess->RegisterMe(electronReaction); manager->AddDiscreteProcess(positronNuclearProcess); G4GammaNuclearReaction * lowEGammaModel(new G4GammaNuclearReaction); lowEGammaModel->SetMaxEnergy(3.5*GeV); G4TheoFSGenerator * highEGammaModel(new G4TheoFSGenerator); highEGammaModel->SetTransport(new G4GeneratorPrecompoundInterface); G4QGSModel * stringModel(new G4QGSModel); stringModel->SetFragmentationModel(new G4ExcitedStringDecay(new G4QGSMFragmentation)); highEGammaModel->SetHighEnergyGenerator(stringModel); highEGammaModel->SetMinEnergy(3.*GeV); highEGammaModel->SetMaxEnergy(100.*TeV); manager=G4Gamma::GammaDefinition()->GetProcessManager(); G4PhotoNuclearProcess * photoNuclearProcess(new G4PhotoNuclearProcess); photoNuclearProcess->RegisterMe(lowEGammaModel); photoNuclearProcess->RegisterMe(highEGammaModel); manager->AddDiscreteProcess(photoNuclearProcess); } void RadmonPhysicsNuclear :: SetCuts(void) { } const RadmonPhysicsInfoList & RadmonPhysicsNuclear :: Provides(void) const { if (infoList.GetNPhysicsInfos()==0) { RadmonPhysicsInfo info; info.SetProcessName("NuclearReaction"); info.SetParticleDefinition(G4Electron::ElectronDefinition()); info.SetMinEnergy(0.*eV); info.SetMaxEnergy(30.*TeV); infoList.InsertPhysicsInfo(info); info.SetParticleDefinition(G4Positron::PositronDefinition()); infoList.InsertPhysicsInfo(info); info.SetParticleDefinition(G4Gamma::GammaDefinition()); info.SetMaxEnergy(100.*TeV); infoList.InsertPhysicsInfo(info); } return infoList; }