G4QGSBinaryNeutronBuilder.cc

Go to the documentation of this file.
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$
00027 //
00028 //---------------------------------------------------------------------------
00029 //
00030 // ClassName:   G4QGSBinaryNeutronBuilder
00031 //
00032 // Author: 2006 G.Folger
00033 //
00034 // Modified:
00035 // 30.03.2009 V.Ivanchenko create cross section by new
00036 //
00037 //----------------------------------------------------------------------------
00038 //
00039 #include "G4QGSBinaryNeutronBuilder.hh"
00040 #include "G4SystemOfUnits.hh"
00041 #include "G4ParticleDefinition.hh"
00042 #include "G4ParticleTable.hh"
00043 #include "G4ProcessManager.hh"
00044 #include "G4NeutronInelasticCrossSection.hh"
00045 #include "G4BGGNucleonInelasticXS.hh"
00046 
00047 G4QGSBinaryNeutronBuilder::
00048 G4QGSBinaryNeutronBuilder(G4bool quasiElastic) 
00049 {
00050   theMin = 12*GeV;
00051   theModel = new G4TheoFSGenerator("QGSB");
00052 
00053   theStringModel = new G4QGSModel< G4QGSParticipants >;
00054   theStringDecay = new G4ExcitedStringDecay(new G4QGSMFragmentation);
00055   theStringModel->SetFragmentationModel(theStringDecay);
00056 
00057   theCascade = new G4BinaryCascade;
00058   thePreEquilib = new G4PreCompoundModel(new G4ExcitationHandler);
00059   theCascade->SetDeExcitation(thePreEquilib);  
00060 
00061   theModel->SetTransport(theCascade);
00062   theModel->SetHighEnergyGenerator(theStringModel);
00063   if (quasiElastic)
00064   {
00065      theQuasiElastic=new G4QuasiElasticChannel;
00066      theModel->SetQuasiElasticChannel(theQuasiElastic);
00067   } else 
00068   {  theQuasiElastic=0;}  
00069 }
00070 
00071 G4QGSBinaryNeutronBuilder::
00072 ~G4QGSBinaryNeutronBuilder() 
00073 {
00074   delete theStringDecay;
00075   delete theStringModel;
00076   delete thePreEquilib;
00077   delete theCascade;
00078   if ( theQuasiElastic ) delete theQuasiElastic;
00079   delete theModel;
00080 }
00081 
00082 void G4QGSBinaryNeutronBuilder::
00083 Build(G4HadronElasticProcess * )
00084 {
00085 }
00086 
00087 void G4QGSBinaryNeutronBuilder::
00088 Build(G4HadronFissionProcess * )
00089 {
00090 }
00091 
00092 void G4QGSBinaryNeutronBuilder::
00093 Build(G4HadronCaptureProcess * )
00094 {
00095 }
00096 
00097 void G4QGSBinaryNeutronBuilder::
00098 Build(G4NeutronInelasticProcess * aP)
00099 {
00100   theModel->SetMinEnergy(theMin);
00101   theModel->SetMaxEnergy(100*TeV);
00102   aP->RegisterMe(theModel);
00103     aP->AddDataSet(new G4BGGNucleonInelasticXS(G4Neutron::Neutron()));
00104 }
00105 

Generated on Mon May 27 17:49:33 2013 for Geant4 by  doxygen 1.4.7