G4QGSPProtonBuilder.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:   G4QGSPProtonBuilder
00031 //
00032 // Author: 2002 J.P. Wellisch
00033 //
00034 // Modified:
00035 // 17.11.2010 G.Folger, use G4CrossSectionPairGG for relativistic rise of cross
00036 //             section at high energies.
00037 // 30.03.2009 V.Ivanchenko create cross section by new
00038 //
00039 //----------------------------------------------------------------------------
00040 //
00041 #include "G4QGSPProtonBuilder.hh"
00042 #include "G4SystemOfUnits.hh"
00043 #include "G4ParticleDefinition.hh"
00044 #include "G4ParticleTable.hh"
00045 #include "G4ProcessManager.hh"
00046 #include "G4ProtonInelasticCrossSection.hh"
00047 #include "G4BGGNucleonInelasticXS.hh"
00048 
00049 G4QGSPProtonBuilder::
00050 G4QGSPProtonBuilder(G4bool quasiElastic, G4bool projectileDiffraction) 
00051  {
00052    theMin = 12*GeV;
00053    theModel = new G4TheoFSGenerator("QGSP");
00054 
00055    theStringModel = new G4QGSModel< G4QGSParticipants >;
00056    theStringDecay = new G4ExcitedStringDecay(theQGSM = new G4QGSMFragmentation);
00057    theStringModel->SetFragmentationModel(theStringDecay);
00058 
00059    theCascade = new G4GeneratorPrecompoundInterface;
00060    thePreEquilib = new G4PreCompoundModel(theHandler = new G4ExcitationHandler);
00061    theCascade->SetDeExcitation(thePreEquilib);  
00062 
00063    theModel->SetTransport(theCascade);
00064    theModel->SetHighEnergyGenerator(theStringModel);
00065    if (quasiElastic)
00066    {
00067       theQuasiElastic=new G4QuasiElasticChannel;
00068       theModel->SetQuasiElasticChannel(theQuasiElastic);
00069    } else 
00070    {  theQuasiElastic=0;}  
00071    if ( projectileDiffraction )
00072    {
00073       theProjectileDiffraction=new G4ProjectileDiffractiveChannel;
00074       theModel->SetProjectileDiffraction(theProjectileDiffraction);
00075    } else 
00076    {  theProjectileDiffraction=0;}
00077  }
00078 
00079 void G4QGSPProtonBuilder::
00080 Build(G4ProtonInelasticProcess * aP)
00081  {
00082    aP->AddDataSet(new G4BGGNucleonInelasticXS(G4Proton::Proton()));
00083    theModel->SetMinEnergy(theMin);
00084    theModel->SetMaxEnergy(100*TeV);
00085    aP->RegisterMe(theModel);
00086  }
00087 
00088 void G4QGSPProtonBuilder::
00089 Build(G4HadronElasticProcess * )
00090  {
00091  }
00092 
00093 G4QGSPProtonBuilder::~G4QGSPProtonBuilder() 
00094  {
00095    delete thePreEquilib;
00096    delete theCascade;
00097    if ( theQuasiElastic ) delete theQuasiElastic;
00098    if ( theProjectileDiffraction ) delete theProjectileDiffraction;
00099    delete theStringDecay;
00100    delete theStringModel;
00101    delete theModel;
00102    delete theQGSM;
00103    //delete theHandler;
00104  }
00105 
00106  // 2002 by J.P. Wellisch

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