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 // GEANT4 tag $Name: BesSim-00-01-24 $ 00027 // 00028 //--------------------------------------------------------------------------- 00029 // 00030 // ClassName: G4MiscLHEPBuilder_CHIPS 00031 // 00032 // Author: 2002 J.P. Wellisch 00033 // 00034 // Modified: 00035 // 16.11.2005 G.Folger: don't keep processes as data members, but new these 00036 // 13.06.2006 G.Folger: (re)move elastic scatterring 00037 // 00038 //---------------------------------------------------------------------------- 00039 // 00040 #include "G4MiscLHEPBuilder_CHIPS.hh" 00041 00042 #include "G4ParticleDefinition.hh" 00043 #include "G4ParticleTable.hh" 00044 #include "G4ProcessManager.hh" 00045 00046 G4MiscLHEPBuilder_CHIPS::G4MiscLHEPBuilder_CHIPS(): wasActivated(false) {} 00047 G4MiscLHEPBuilder_CHIPS::~G4MiscLHEPBuilder_CHIPS() 00048 { 00049 if(wasActivated) delete inelastic; 00050 } 00051 00052 void G4MiscLHEPBuilder_CHIPS::Build() 00053 { 00054 G4ProcessManager * aProcMan = 0; 00055 wasActivated = true; 00056 00057 // anti-Proton 00058 //theAntiProtonInelastic = new G4AntiProtonInelasticProcess(); 00059 aProcMan = G4AntiProton::AntiProton()->GetProcessManager(); 00060 //theLEAntiProtonModel = new G4LEAntiProtonInelastic(); 00061 //theHEAntiProtonModel = new G4HEAntiProtonInelastic(); 00062 //theHEAntiProtonModel->SetMaxEnergy(100*TeV); 00063 //theAntiProtonInelastic->RegisterMe(theLEAntiProtonModel); 00064 //theAntiProtonInelastic->RegisterMe(theHEAntiProtonModel); 00065 //aProcMan->AddDiscreteProcess(theAntiProtonInelastic); 00066 //CHIPS 00067 inelastic = new G4QInelastic(); 00068 aProcMan->AddDiscreteProcess(inelastic); 00069 00070 // AntiNeutron 00071 theAntiNeutronInelastic = new G4AntiNeutronInelasticProcess(); 00072 aProcMan = G4AntiNeutron::AntiNeutron()->GetProcessManager(); 00073 theLEAntiNeutronModel = new G4LEAntiNeutronInelastic(); 00074 theHEAntiNeutronModel = new G4HEAntiNeutronInelastic(); 00075 theHEAntiNeutronModel->SetMaxEnergy(100*TeV); 00076 theAntiNeutronInelastic->RegisterMe(theLEAntiNeutronModel); 00077 theAntiNeutronInelastic->RegisterMe(theHEAntiNeutronModel); 00078 aProcMan->AddDiscreteProcess(theAntiNeutronInelastic); 00079 00080 // Lambda 00081 theLambdaInelastic = new G4LambdaInelasticProcess(); 00082 aProcMan = G4Lambda::Lambda()->GetProcessManager(); 00083 theLELambdaModel = new G4LELambdaInelastic(); 00084 theHELambdaModel = new G4HELambdaInelastic(); 00085 theHELambdaModel->SetMaxEnergy(100*TeV); 00086 theLambdaInelastic->RegisterMe(theLELambdaModel); 00087 theLambdaInelastic->RegisterMe(theHELambdaModel); 00088 aProcMan->AddDiscreteProcess(theLambdaInelastic); 00089 00090 // AntiLambda 00091 theAntiLambdaInelastic = new G4AntiLambdaInelasticProcess(); 00092 aProcMan = G4AntiLambda::AntiLambda()->GetProcessManager(); 00093 theLEAntiLambdaModel = new G4LEAntiLambdaInelastic(); 00094 theHEAntiLambdaModel = new G4HEAntiLambdaInelastic(); 00095 theHEAntiLambdaModel->SetMaxEnergy(100*TeV); 00096 theAntiLambdaInelastic->RegisterMe(theLEAntiLambdaModel); 00097 theAntiLambdaInelastic->RegisterMe(theHEAntiLambdaModel); 00098 aProcMan->AddDiscreteProcess(theAntiLambdaInelastic); 00099 00100 // SigmaMinus 00101 theSigmaMinusInelastic = new G4SigmaMinusInelasticProcess(); 00102 aProcMan = G4SigmaMinus::SigmaMinus()->GetProcessManager(); 00103 theLESigmaMinusModel = new G4LESigmaMinusInelastic(); 00104 theHESigmaMinusModel = new G4HESigmaMinusInelastic(); 00105 theHESigmaMinusModel->SetMaxEnergy(100*TeV); 00106 theSigmaMinusInelastic->RegisterMe(theLESigmaMinusModel); 00107 theSigmaMinusInelastic->RegisterMe(theHESigmaMinusModel); 00108 aProcMan->AddDiscreteProcess(theSigmaMinusInelastic); 00109 00110 // anti-SigmaMinus 00111 theAntiSigmaMinusInelastic = new G4AntiSigmaMinusInelasticProcess(); 00112 aProcMan = G4AntiSigmaMinus::AntiSigmaMinus()->GetProcessManager(); 00113 theLEAntiSigmaMinusModel = new G4LEAntiSigmaMinusInelastic(); 00114 theHEAntiSigmaMinusModel = new G4HEAntiSigmaMinusInelastic(); 00115 theHEAntiSigmaMinusModel->SetMaxEnergy(100*TeV); 00116 theAntiSigmaMinusInelastic->RegisterMe(theLEAntiSigmaMinusModel); 00117 theAntiSigmaMinusInelastic->RegisterMe(theHEAntiSigmaMinusModel); 00118 aProcMan->AddDiscreteProcess(theAntiSigmaMinusInelastic); 00119 00120 // SigmaPlus 00121 theSigmaPlusInelastic = new G4SigmaPlusInelasticProcess(); 00122 aProcMan = G4SigmaPlus::SigmaPlus()->GetProcessManager(); 00123 theLESigmaPlusModel = new G4LESigmaPlusInelastic(); 00124 theHESigmaPlusModel = new G4HESigmaPlusInelastic(); 00125 theHESigmaPlusModel->SetMaxEnergy(100*TeV); 00126 theSigmaPlusInelastic->RegisterMe(theLESigmaPlusModel); 00127 theSigmaPlusInelastic->RegisterMe(theHESigmaPlusModel); 00128 aProcMan->AddDiscreteProcess(theSigmaPlusInelastic); 00129 00130 // anti-SigmaPlus 00131 theAntiSigmaPlusInelastic = new G4AntiSigmaPlusInelasticProcess(); 00132 aProcMan = G4AntiSigmaPlus::AntiSigmaPlus()->GetProcessManager(); 00133 theLEAntiSigmaPlusModel = new G4LEAntiSigmaPlusInelastic(); 00134 theHEAntiSigmaPlusModel = new G4HEAntiSigmaPlusInelastic(); 00135 theHEAntiSigmaPlusModel->SetMaxEnergy(100*TeV); 00136 theAntiSigmaPlusInelastic->RegisterMe(theLEAntiSigmaPlusModel); 00137 theAntiSigmaPlusInelastic->RegisterMe(theHEAntiSigmaPlusModel); 00138 aProcMan->AddDiscreteProcess(theAntiSigmaPlusInelastic); 00139 00140 // XiMinus 00141 theXiMinusInelastic = new G4XiMinusInelasticProcess(); 00142 aProcMan = G4XiMinus::XiMinus()->GetProcessManager(); 00143 theLEXiMinusModel = new G4LEXiMinusInelastic(); 00144 theHEXiMinusModel = new G4HEXiMinusInelastic(); 00145 theHEXiMinusModel->SetMaxEnergy(100*TeV); 00146 theXiMinusInelastic->RegisterMe(theLEXiMinusModel); 00147 theXiMinusInelastic->RegisterMe(theHEXiMinusModel); 00148 aProcMan->AddDiscreteProcess(theXiMinusInelastic); 00149 00150 // anti-XiMinus 00151 theAntiXiMinusInelastic = new G4AntiXiMinusInelasticProcess(); 00152 aProcMan = G4AntiXiMinus::AntiXiMinus()->GetProcessManager(); 00153 theLEAntiXiMinusModel = new G4LEAntiXiMinusInelastic(); 00154 theHEAntiXiMinusModel = new G4HEAntiXiMinusInelastic(); 00155 theHEAntiXiMinusModel->SetMaxEnergy(100*TeV); 00156 theAntiXiMinusInelastic->RegisterMe(theLEAntiXiMinusModel); 00157 theAntiXiMinusInelastic->RegisterMe(theHEAntiXiMinusModel); 00158 aProcMan->AddDiscreteProcess(theAntiXiMinusInelastic); 00159 00160 // XiZero 00161 theXiZeroInelastic = new G4XiZeroInelasticProcess(); 00162 aProcMan = G4XiZero::XiZero()->GetProcessManager(); 00163 theLEXiZeroModel = new G4LEXiZeroInelastic(); 00164 theHEXiZeroModel = new G4HEXiZeroInelastic(); 00165 theHEXiZeroModel->SetMaxEnergy(100*TeV); 00166 theXiZeroInelastic->RegisterMe(theLEXiZeroModel); 00167 theXiZeroInelastic->RegisterMe(theHEXiZeroModel); 00168 aProcMan->AddDiscreteProcess(theXiZeroInelastic); 00169 00170 // anti-XiZero 00171 theAntiXiZeroInelastic = new G4AntiXiZeroInelasticProcess(); 00172 aProcMan = G4AntiXiZero::AntiXiZero()->GetProcessManager(); 00173 theLEAntiXiZeroModel = new G4LEAntiXiZeroInelastic(); 00174 theHEAntiXiZeroModel = new G4HEAntiXiZeroInelastic(); 00175 theHEAntiXiZeroModel->SetMaxEnergy(100*TeV); 00176 theAntiXiZeroInelastic->RegisterMe(theLEAntiXiZeroModel); 00177 theAntiXiZeroInelastic->RegisterMe(theHEAntiXiZeroModel); 00178 aProcMan->AddDiscreteProcess(theAntiXiZeroInelastic); 00179 00180 // OmegaMinus 00181 theOmegaMinusInelastic = new G4OmegaMinusInelasticProcess(); 00182 aProcMan = G4OmegaMinus::OmegaMinus()->GetProcessManager(); 00183 theLEOmegaMinusModel = new G4LEOmegaMinusInelastic(); 00184 theHEOmegaMinusModel = new G4HEOmegaMinusInelastic(); 00185 theHEOmegaMinusModel->SetMaxEnergy(100*TeV); 00186 theOmegaMinusInelastic->RegisterMe(theLEOmegaMinusModel); 00187 theOmegaMinusInelastic->RegisterMe(theHEOmegaMinusModel); 00188 aProcMan->AddDiscreteProcess(theOmegaMinusInelastic); 00189 00190 // anti-OmegaMinus 00191 theAntiOmegaMinusInelastic = new G4AntiOmegaMinusInelasticProcess(); 00192 aProcMan = G4AntiOmegaMinus::AntiOmegaMinus()->GetProcessManager(); 00193 theLEAntiOmegaMinusModel = new G4LEAntiOmegaMinusInelastic(); 00194 theHEAntiOmegaMinusModel = new G4HEAntiOmegaMinusInelastic(); 00195 theHEAntiOmegaMinusModel->SetMaxEnergy(100*TeV); 00196 theAntiOmegaMinusInelastic->RegisterMe(theLEAntiOmegaMinusModel); 00197 theAntiOmegaMinusInelastic->RegisterMe(theHEAntiOmegaMinusModel); 00198 aProcMan->AddDiscreteProcess(theAntiOmegaMinusInelastic); 00199 } 00200 00201 // 2002 by J.P. Wellisch