Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 // << 27 //-------------------------------------------- << 28 // << 29 // ClassName: G4QGSPPiKBuilder << 30 // << 31 // Author: 2002 J.P. Wellisch << 32 // << 33 // Modified: << 34 // 17.11.2010 G.Folger, use G4CrossSectionPair << 35 // section at high energies. << 36 // 30.03.2009 V.Ivanchenko create cross sectio << 37 // << 38 //-------------------------------------------- << 39 // << 40 #include "G4QGSPPiKBuilder.hh" 26 #include "G4QGSPPiKBuilder.hh" 41 #include "G4SystemOfUnits.hh" << 42 #include "G4ParticleDefinition.hh" 27 #include "G4ParticleDefinition.hh" 43 #include "G4ParticleTable.hh" 28 #include "G4ParticleTable.hh" 44 #include "G4ProcessManager.hh" 29 #include "G4ProcessManager.hh" 45 #include "G4BGGPionInelasticXS.hh" << 46 #include "G4HadronicParameters.hh" << 47 << 48 30 49 G4QGSPPiKBuilder:: 31 G4QGSPPiKBuilder:: 50 G4QGSPPiKBuilder(G4bool quasiElastic) 32 G4QGSPPiKBuilder(G4bool quasiElastic) 51 { << 33 { 52 theMin = G4HadronicParameters::Instance()->G << 34 theMin = 12*GeV; 53 theModel = new G4TheoFSGenerator("QGSP"); << 35 theModel = new G4TheoFSGenerator; 54 << 36 55 G4QGSModel< G4QGSParticipants >* theStringMo << 37 theStringModel = new G4QGSModel< G4QGSParticipants >; 56 new G4QGSModel< G4QGSParticipants >; << 38 theStringDecay = new G4ExcitedStringDecay(new G4QGSMFragmentation); 57 G4ExcitedStringDecay* theStringDecay = << 58 new G4ExcitedStringDecay(new G4QGSMFragmen << 59 theStringModel->SetFragmentationModel(theStr 39 theStringModel->SetFragmentationModel(theStringDecay); >> 40 60 41 61 G4GeneratorPrecompoundInterface* theCascade << 42 theCascade = new G4GeneratorPrecompoundInterface; 62 new G4GeneratorPrecompoundInterface(); << 43 thePreEquilib = new G4PreCompoundModel(new G4ExcitationHandler); >> 44 theCascade->SetDeExcitation(thePreEquilib); 63 45 64 theModel->SetTransport(theCascade); << 65 theModel->SetHighEnergyGenerator(theStringMo 46 theModel->SetHighEnergyGenerator(theStringModel); 66 if (quasiElastic) 47 if (quasiElastic) 67 { << 48 { 68 theModel->SetQuasiElasticChannel(new G4Q << 49 theQuasiElastic=new G4QuasiElasticChannel; 69 } << 50 theModel->SetQuasiElasticChannel(theQuasiElastic); >> 51 } else >> 52 { theQuasiElastic=0;} >> 53 theModel->SetTransport(theCascade); 70 } 54 } 71 << 55 72 G4QGSPPiKBuilder:: 56 G4QGSPPiKBuilder:: 73 ~G4QGSPPiKBuilder() 57 ~G4QGSPPiKBuilder() 74 { 58 { >> 59 delete theCascade; >> 60 delete thePreEquilib; >> 61 if ( theQuasiElastic ) delete theQuasiElastic; >> 62 delete theStringDecay; >> 63 delete theStringModel; >> 64 delete theModel; >> 65 } >> 66 >> 67 void G4QGSPPiKBuilder:: >> 68 Build(G4HadronElasticProcess * ) {} >> 69 >> 70 void G4QGSPPiKBuilder:: >> 71 Build(G4PionPlusInelasticProcess * aP) >> 72 { >> 73 theModel->SetMinEnergy(theMin); >> 74 theModel->SetMaxEnergy(100*TeV); >> 75 aP->AddDataSet(&thePiData); >> 76 aP->RegisterMe(theModel); >> 77 } >> 78 >> 79 void G4QGSPPiKBuilder:: >> 80 Build(G4PionMinusInelasticProcess * aP) >> 81 { >> 82 theModel->SetMinEnergy(theMin); >> 83 theModel->SetMaxEnergy(100*TeV); >> 84 aP->AddDataSet(&thePiData); >> 85 aP->RegisterMe(theModel); 75 } 86 } 76 87 77 void G4QGSPPiKBuilder:: 88 void G4QGSPPiKBuilder:: 78 Build(G4HadronInelasticProcess * aP) << 89 Build(G4KaonPlusInelasticProcess * aP) 79 { 90 { 80 theModel->SetMinEnergy(theMin); 91 theModel->SetMinEnergy(theMin); 81 theModel->SetMaxEnergy( G4HadronicParameters << 92 theModel->SetMaxEnergy(100*TeV); 82 if ( aP->GetParticleDefinition() == G4PionPl << 83 aP->AddDataSet( new G4BGGPionInelasticXS( << 84 } else if ( aP->GetParticleDefinition() == G << 85 aP->AddDataSet( new G4BGGPionInelasticXS( << 86 } << 87 aP->RegisterMe(theModel); 93 aP->RegisterMe(theModel); 88 } 94 } >> 95 >> 96 void G4QGSPPiKBuilder:: >> 97 Build(G4KaonMinusInelasticProcess * aP) >> 98 { >> 99 theModel->SetMinEnergy(theMin); >> 100 theModel->SetMaxEnergy(100*TeV); >> 101 aP->RegisterMe(theModel); >> 102 } >> 103 >> 104 void G4QGSPPiKBuilder:: >> 105 Build(G4KaonZeroLInelasticProcess * aP) >> 106 { >> 107 theModel->SetMinEnergy(theMin); >> 108 theModel->SetMaxEnergy(100*TeV); >> 109 aP->RegisterMe(theModel); >> 110 } >> 111 >> 112 void G4QGSPPiKBuilder:: >> 113 Build(G4KaonZeroSInelasticProcess * aP) >> 114 { >> 115 theModel->SetMinEnergy(theMin); >> 116 theModel->SetMaxEnergy(100*TeV); >> 117 aP->RegisterMe(theModel); >> 118 } >> 119 >> 120 // 2002 by J.P. Wellisch 89 121