| 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 //-------------------------------------------- 26 //---------------------------------------------------------------------------
27 // ClassName: G4HyperonQGSPBuilder 27 // ClassName: G4HyperonQGSPBuilder
28 // Author: Alberto Ribon 28 // Author: Alberto Ribon
29 // Date: May 2020 29 // Date: May 2020
30 // Modified: 30 // Modified:
31 //-------------------------------------------- 31 //---------------------------------------------------------------------------
32 32
33 #include "G4HyperonQGSPBuilder.hh" 33 #include "G4HyperonQGSPBuilder.hh"
34 #include "G4SystemOfUnits.hh" 34 #include "G4SystemOfUnits.hh"
35 #include "G4ParticleDefinition.hh" 35 #include "G4ParticleDefinition.hh"
36 #include "G4ParticleTable.hh" 36 #include "G4ParticleTable.hh"
37 #include "G4ProcessManager.hh" 37 #include "G4ProcessManager.hh"
38 #include "G4HadronInelasticProcess.hh" << 38 #include "G4LambdaInelasticProcess.hh"
>> 39 #include "G4AntiLambdaInelasticProcess.hh"
>> 40 #include "G4SigmaMinusInelasticProcess.hh"
>> 41 #include "G4AntiSigmaMinusInelasticProcess.hh"
>> 42 #include "G4SigmaPlusInelasticProcess.hh"
>> 43 #include "G4AntiSigmaPlusInelasticProcess.hh"
>> 44 #include "G4XiMinusInelasticProcess.hh"
>> 45 #include "G4AntiXiMinusInelasticProcess.hh"
>> 46 #include "G4XiZeroInelasticProcess.hh"
>> 47 #include "G4AntiXiZeroInelasticProcess.hh"
>> 48 #include "G4OmegaMinusInelasticProcess.hh"
>> 49 #include "G4AntiOmegaMinusInelasticProcess.hh"
39 #include "G4CrossSectionInelastic.hh" 50 #include "G4CrossSectionInelastic.hh"
40 #include "G4ComponentGGHadronNucleusXsc.hh" 51 #include "G4ComponentGGHadronNucleusXsc.hh"
41 #include "G4HadronicParameters.hh" 52 #include "G4HadronicParameters.hh"
42 #include "G4TheoFSGenerator.hh" 53 #include "G4TheoFSGenerator.hh"
43 #include "G4GeneratorPrecompoundInterface.hh" 54 #include "G4GeneratorPrecompoundInterface.hh"
44 #include "G4QGSModel.hh" 55 #include "G4QGSModel.hh"
45 #include "G4QGSParticipants.hh" 56 #include "G4QGSParticipants.hh"
46 #include "G4QGSMFragmentation.hh" 57 #include "G4QGSMFragmentation.hh"
47 #include "G4ExcitedStringDecay.hh" 58 #include "G4ExcitedStringDecay.hh"
48 #include "G4QuasiElasticChannel.hh" 59 #include "G4QuasiElasticChannel.hh"
49 60
50 61
51 G4HyperonQGSPBuilder::G4HyperonQGSPBuilder( G4 62 G4HyperonQGSPBuilder::G4HyperonQGSPBuilder( G4bool quasiElastic ) {
52 theHyperonQGSP = new G4TheoFSGenerator( "QGS 63 theHyperonQGSP = new G4TheoFSGenerator( "QGSP" );
53 G4QGSModel< G4QGSParticipants >* theStringMo 64 G4QGSModel< G4QGSParticipants >* theStringModel = new G4QGSModel< G4QGSParticipants >;
54 G4ExcitedStringDecay* theStringDecay = new G 65 G4ExcitedStringDecay* theStringDecay = new G4ExcitedStringDecay( new G4QGSMFragmentation );
55 theStringModel->SetFragmentationModel( theSt 66 theStringModel->SetFragmentationModel( theStringDecay );
56 G4GeneratorPrecompoundInterface* theCascade 67 G4GeneratorPrecompoundInterface* theCascade = new G4GeneratorPrecompoundInterface;
57 theHyperonQGSP->SetTransport( theCascade ); 68 theHyperonQGSP->SetTransport( theCascade );
58 theHyperonQGSP->SetHighEnergyGenerator( theS 69 theHyperonQGSP->SetHighEnergyGenerator( theStringModel );
59 if ( quasiElastic ) theHyperonQGSP->SetQuasi 70 if ( quasiElastic ) theHyperonQGSP->SetQuasiElasticChannel( new G4QuasiElasticChannel );
60 // It is recommended to use QGSP above >~ 12 71 // It is recommended to use QGSP above >~ 12 GeV .
61 theMin = G4HadronicParameters::Instance()->G 72 theMin = G4HadronicParameters::Instance()->GetMinEnergyTransitionQGS_FTF();
62 theMax = G4HadronicParameters::Instance()->G 73 theMax = G4HadronicParameters::Instance()->GetMaxEnergy();
63 theHyperonQGSP->SetMinEnergy( theMin ); 74 theHyperonQGSP->SetMinEnergy( theMin );
64 theHyperonQGSP->SetMaxEnergy( theMax ); 75 theHyperonQGSP->SetMaxEnergy( theMax );
65 theInelasticCrossSection = new G4CrossSectio 76 theInelasticCrossSection = new G4CrossSectionInelastic( new G4ComponentGGHadronNucleusXsc );
66 } 77 }
67 78
68 79
69 G4HyperonQGSPBuilder::~G4HyperonQGSPBuilder() 80 G4HyperonQGSPBuilder::~G4HyperonQGSPBuilder() {}
70 81
71 82
72 void G4HyperonQGSPBuilder::Build( G4HadronInel << 83 void G4HyperonQGSPBuilder::Build( G4LambdaInelasticProcess* aP ) {
>> 84 theHyperonQGSP->SetMinEnergy( theMin );
>> 85 theHyperonQGSP->SetMaxEnergy( theMax );
>> 86 aP->RegisterMe( theHyperonQGSP );
>> 87 aP->AddDataSet( theInelasticCrossSection );
>> 88 }
>> 89
>> 90 void G4HyperonQGSPBuilder::Build( G4AntiLambdaInelasticProcess* aP ) {
>> 91 theHyperonQGSP->SetMinEnergy( theMin );
>> 92 theHyperonQGSP->SetMaxEnergy( theMax );
>> 93 aP->RegisterMe( theHyperonQGSP );
>> 94 aP->AddDataSet( theInelasticCrossSection );
>> 95 }
>> 96
>> 97
>> 98 void G4HyperonQGSPBuilder::Build( G4SigmaMinusInelasticProcess* aP ) {
>> 99 theHyperonQGSP->SetMinEnergy( theMin );
>> 100 theHyperonQGSP->SetMaxEnergy( theMax );
>> 101 aP->RegisterMe( theHyperonQGSP );
>> 102 aP->AddDataSet( theInelasticCrossSection );
>> 103 }
>> 104
>> 105 void G4HyperonQGSPBuilder::Build( G4AntiSigmaMinusInelasticProcess* aP ) {
>> 106 theHyperonQGSP->SetMinEnergy( theMin );
>> 107 theHyperonQGSP->SetMaxEnergy( theMax );
>> 108 aP->RegisterMe( theHyperonQGSP );
>> 109 aP->AddDataSet( theInelasticCrossSection );
>> 110 }
>> 111
>> 112
>> 113 void G4HyperonQGSPBuilder::Build( G4SigmaPlusInelasticProcess* aP ) {
>> 114 theHyperonQGSP->SetMinEnergy( theMin );
>> 115 theHyperonQGSP->SetMaxEnergy( theMax );
>> 116 aP->RegisterMe( theHyperonQGSP );
>> 117 aP->AddDataSet( theInelasticCrossSection );
>> 118 }
>> 119
>> 120 void G4HyperonQGSPBuilder::Build( G4AntiSigmaPlusInelasticProcess* aP ) {
>> 121 theHyperonQGSP->SetMinEnergy( theMin );
>> 122 theHyperonQGSP->SetMaxEnergy( theMax );
>> 123 aP->RegisterMe( theHyperonQGSP );
>> 124 aP->AddDataSet( theInelasticCrossSection );
>> 125 }
>> 126
>> 127
>> 128 void G4HyperonQGSPBuilder::Build( G4XiMinusInelasticProcess* aP ) {
>> 129 theHyperonQGSP->SetMinEnergy( theMin );
>> 130 theHyperonQGSP->SetMaxEnergy( theMax );
>> 131 aP->RegisterMe( theHyperonQGSP );
>> 132 aP->AddDataSet( theInelasticCrossSection );
>> 133 }
>> 134
>> 135 void G4HyperonQGSPBuilder::Build( G4AntiXiMinusInelasticProcess* aP ) {
>> 136 theHyperonQGSP->SetMinEnergy( theMin );
>> 137 theHyperonQGSP->SetMaxEnergy( theMax );
>> 138 aP->RegisterMe( theHyperonQGSP );
>> 139 aP->AddDataSet( theInelasticCrossSection );
>> 140 }
>> 141
>> 142
>> 143 void G4HyperonQGSPBuilder::Build( G4XiZeroInelasticProcess* aP ) {
>> 144 theHyperonQGSP->SetMinEnergy( theMin );
>> 145 theHyperonQGSP->SetMaxEnergy( theMax );
>> 146 aP->RegisterMe( theHyperonQGSP );
>> 147 aP->AddDataSet( theInelasticCrossSection );
>> 148 }
>> 149
>> 150 void G4HyperonQGSPBuilder::Build( G4AntiXiZeroInelasticProcess* aP ) {
>> 151 theHyperonQGSP->SetMinEnergy( theMin );
>> 152 theHyperonQGSP->SetMaxEnergy( theMax );
>> 153 aP->RegisterMe( theHyperonQGSP );
>> 154 aP->AddDataSet( theInelasticCrossSection );
>> 155 }
>> 156
>> 157
>> 158 void G4HyperonQGSPBuilder::Build( G4OmegaMinusInelasticProcess* aP ) {
>> 159 theHyperonQGSP->SetMinEnergy( theMin );
>> 160 theHyperonQGSP->SetMaxEnergy( theMax );
>> 161 aP->RegisterMe( theHyperonQGSP );
>> 162 aP->AddDataSet( theInelasticCrossSection );
>> 163 }
>> 164
>> 165 void G4HyperonQGSPBuilder::Build( G4AntiOmegaMinusInelasticProcess* aP ) {
73 theHyperonQGSP->SetMinEnergy( theMin ); 166 theHyperonQGSP->SetMinEnergy( theMin );
74 theHyperonQGSP->SetMaxEnergy( theMax ); 167 theHyperonQGSP->SetMaxEnergy( theMax );
75 aP->RegisterMe( theHyperonQGSP ); 168 aP->RegisterMe( theHyperonQGSP );
76 aP->AddDataSet( theInelasticCrossSection ); 169 aP->AddDataSet( theInelasticCrossSection );
77 } 170 }
78 171