00001 // 00002 // ******************************************************************** 00003 // * License and Disclaimer * 00004 // * * 00005 // * The Geant4 software is copyright of the Copyright Holders of * 00006 // * the Geant4 Collaboration. It is provided under the terms and * 00007 // * conditions of the Geant4 Software License, included in the file * 00008 // * LICENSE and available at http://cern.ch/geant4/license . These * 00009 // * include a list of copyright holders. * 00010 // * * 00011 // * Neither the authors of this software system, nor their employing * 00012 // * institutes,nor the agencies providing financial support for this * 00013 // * work make any representation or warranty, express or implied, * 00014 // * regarding this software system or assume any liability for its * 00015 // * use. Please see the license in the file LICENSE and URL above * 00016 // * for the full disclaimer and the limitation of liability. * 00017 // * * 00018 // * This code implementation is the result of the scientific and * 00019 // * technical work of the GEANT4 collaboration. * 00020 // * By using, copying, modifying or distributing the software (or * 00021 // * any work based on the software) you agree to acknowledge its * 00022 // * use in resulting scientific publications, and indicate your * 00023 // * acceptance of all terms of the Geant4 Software license. * 00024 // ******************************************************************** 00025 // 00026 // $Id: HadronPhysicsQGSP_BERT_CHIPS.cc,v 1.3 2015/03/17 05:50:57 sunss Exp $ 00027 // GEANT4 tag $Name: BesSim-00-01-24 $ 00028 // 00029 //--------------------------------------------------------------------------- 00030 // 00031 // ClassName: HadronPhysicsQGSP_BERT_CHIPS 00032 // 00033 // Author: 2002 J.P. Wellisch 00034 // 00035 // Modified: 00036 // 23.11.2005 G.Folger: migration to non static particles 00037 // 08.06.2006 V.Ivanchenko: remove stopping 00038 // 20.06.2006 G.Folger: Bertini applies to Kaons, i.e. use SetMinEnergy instead of SetMinPionEnergy 00039 // 25.04.2007 G.Folger: Add code for quasielastic 00040 // 10.12.2007 G.Folger: Add projectilediffrative option for proton/neutron, off by default 00041 // 00042 //---------------------------------------------------------------------------- 00043 // 00044 #include "HadronPhysicsQGSP_BERT_CHIPS.hh" 00045 00046 #include "globals.hh" 00047 #include "G4ios.hh" 00048 #include <iomanip> 00049 #include "G4ParticleDefinition.hh" 00050 #include "G4ParticleTable.hh" 00051 00052 #include "G4MesonConstructor.hh" 00053 #include "G4BaryonConstructor.hh" 00054 #include "G4ShortLivedConstructor.hh" 00055 00056 HadronPhysicsQGSP_BERT_CHIPS::HadronPhysicsQGSP_BERT_CHIPS(const G4String& name, G4bool quasiElastic) 00057 : G4VPhysicsConstructor(name) , QuasiElastic(quasiElastic) 00058 { 00059 ProjectileDiffraction=false; 00060 } 00061 00062 void HadronPhysicsQGSP_BERT_CHIPS::CreateModels() 00063 { 00064 theNeutrons=new G4NeutronBuilder; 00065 theNeutrons->RegisterMe(theQGSPNeutron=new G4QGSPNeutronBuilder(QuasiElastic, ProjectileDiffraction)); 00066 theNeutrons->RegisterMe(theLEPNeutron=new G4LEPNeutronBuilder); 00067 theLEPNeutron->SetMinInelasticEnergy(9.5*GeV); 00068 theLEPNeutron->SetMaxInelasticEnergy(25*GeV); 00069 00070 theNeutrons->RegisterMe(theBertiniNeutron=new G4BertiniNeutronBuilder); 00071 theBertiniNeutron->SetMinEnergy(0.0*GeV); 00072 theBertiniNeutron->SetMaxEnergy(9.9*GeV); 00073 00074 thePro=new G4ProtonBuilder; 00075 thePro->RegisterMe(theQGSPPro=new G4QGSPProtonBuilder(QuasiElastic, ProjectileDiffraction)); 00076 thePro->RegisterMe(theLEPPro=new G4LEPProtonBuilder); 00077 theLEPPro->SetMinEnergy(9.5*GeV); 00078 theLEPPro->SetMaxEnergy(25*GeV); 00079 00080 thePro->RegisterMe(theBertiniPro=new G4BertiniProtonBuilder); 00081 theBertiniPro->SetMaxEnergy(9.9*GeV); 00082 00083 thePiK=new G4PiKBuilder; 00084 thePiK->RegisterMe(theQGSPPiK=new G4QGSPPiKBuilder(QuasiElastic)); 00085 thePiK->RegisterMe(theLEPPiK=new G4LEPPiKBuilder); 00086 theLEPPiK->SetMaxEnergy(25*GeV); 00087 theLEPPiK->SetMinEnergy(9.5*GeV); 00088 00089 thePiK->RegisterMe(theBertiniPiK=new G4BertiniPiKBuilder); 00090 theBertiniPiK->SetMaxEnergy(9.9*GeV); 00091 00092 theMiscLHEP=new G4MiscLHEPBuilder_CHIPS; 00093 } 00094 00095 HadronPhysicsQGSP_BERT_CHIPS::~HadronPhysicsQGSP_BERT_CHIPS() 00096 { 00097 delete theMiscLHEP; 00098 delete theQGSPNeutron; 00099 delete theLEPNeutron; 00100 delete theBertiniNeutron; 00101 delete theQGSPPro; 00102 delete theLEPPro; 00103 delete thePro; 00104 delete theBertiniPro; 00105 delete theQGSPPiK; 00106 delete theLEPPiK; 00107 delete theBertiniPiK; 00108 delete thePiK; 00109 } 00110 00111 void HadronPhysicsQGSP_BERT_CHIPS::ConstructParticle() 00112 { 00113 G4MesonConstructor pMesonConstructor; 00114 pMesonConstructor.ConstructParticle(); 00115 00116 G4BaryonConstructor pBaryonConstructor; 00117 pBaryonConstructor.ConstructParticle(); 00118 00119 G4ShortLivedConstructor pShortLivedConstructor; 00120 pShortLivedConstructor.ConstructParticle(); 00121 } 00122 00123 #include "G4ProcessManager.hh" 00124 void HadronPhysicsQGSP_BERT_CHIPS::ConstructProcess() 00125 { 00126 CreateModels(); 00127 theNeutrons->Build(); 00128 thePro->Build(); 00129 thePiK->Build(); 00130 theMiscLHEP->Build(); 00131 } 00132