source: trunk/source/processes/hadronic/models/high_energy/src/G4HESigmaZeroInelastic.cc @ 962

Last change on this file since 962 was 962, checked in by garnier, 15 years ago

update processes

File size: 3.9 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer                                           *
4// *                                                                  *
5// * The  Geant4 software  is  copyright of the Copyright Holders  of *
6// * the Geant4 Collaboration.  It is provided  under  the terms  and *
7// * conditions of the Geant4 Software License,  included in the file *
8// * LICENSE and available at  http://cern.ch/geant4/license .  These *
9// * include a list of copyright holders.                             *
10// *                                                                  *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work  make  any representation or  warranty, express or implied, *
14// * regarding  this  software system or assume any liability for its *
15// * use.  Please see the license in the file  LICENSE  and URL above *
16// * for the full disclaimer and the limitation of liability.         *
17// *                                                                  *
18// * This  code  implementation is the result of  the  scientific and *
19// * technical work of the GEANT4 collaboration.                      *
20// * By using,  copying,  modifying or  distributing the software (or *
21// * any work based  on the software)  you  agree  to acknowledge its *
22// * use  in  resulting  scientific  publications,  and indicate your *
23// * acceptance of all terms of the Geant4 Software license.          *
24// ********************************************************************
25//
26//
27// $Id: G4HESigmaZeroInelastic.cc,v 1.11 2006/06/29 20:30:44 gunter Exp $
28// GEANT4 tag $Name: geant4-09-02-ref-02 $
29//
30//
31
32#include "globals.hh"
33#include "G4ios.hh"
34
35//
36// G4 Process: Gheisha High Energy Collision model.
37// This includes the high energy cascading model, the two-body-resonance model
38// and the low energy two-body model. Not included are the low energy stuff like
39// nuclear reactions, nuclear fission without any cascading and all processes for
40// particles at rest. 
41// First work done by J.L.Chuma and F.W.Jones, TRIUMF, June 96. 
42// H. Fesefeldt, RWTH-Aachen, 23-October-1996
43// Last modified: 29-July-1998
44 
45#include "G4HESigmaZeroInelastic.hh"
46#include "G4Gamma.hh"
47
48G4HadFinalState *  G4HESigmaZeroInelastic::
49ApplyYourself( const G4HadProjectile &aTrack, G4Nucleus &targetNucleus )
50  {
51    G4HEVector * pv = new G4HEVector[MAXPART];
52    const G4HadProjectile *aParticle = &aTrack;
53//    G4DynamicParticle *originalTarget = targetNucleus.ReturnTargetParticle();
54    G4HEVector incidentParticle(aParticle);
55     
56    G4HELambdaInelastic theLambdaInelastic;
57    theLambdaInelastic.SetMaxNumberOfSecondaries(MAXPART);
58    theLambdaInelastic.SetVerboseLevel(verboseLevel);
59   
60    G4double incidentTotalMomentum = incidentParticle.getTotalMomentum();
61    G4double pgam = G4UniformRand()*incidentTotalMomentum*0.75;
62    G4HEVector incidentLambda; 
63    incidentLambda.SmulAndUpdate( incidentParticle, 
64                                 (incidentTotalMomentum - pgam)/incidentTotalMomentum);                   
65    G4DynamicParticle * aLambda = new G4DynamicParticle();
66    aLambda->SetDefinition(G4Lambda::Lambda());
67    aLambda->SetMomentum(incidentLambda.getMomentum());
68    G4HadProjectile aLambdaTrack(*aLambda);
69    G4HadFinalState * result = theLambdaInelastic.ApplyYourself(aLambdaTrack, targetNucleus);           
70    vecLength = theLambdaInelastic.GetNumberOfSecondaries();
71   
72    pv[vecLength] = Gamma;
73    pv[vecLength].setMomentum(incidentParticle.getMomentum());
74    pv[vecLength].SmulAndUpdate( pv[vecLength],pgam/incidentTotalMomentum);
75    G4DynamicParticle * aPhoton = new G4DynamicParticle();
76    aPhoton->SetDefinition(G4Gamma::Gamma());
77    aPhoton->SetMomentum(pv[vecLength].getMomentum());
78    result->AddSecondary(aPhoton);
79    delete [] pv;
80    return result;
81  } 
82       
83
84
85
86
87
88
89
90
Note: See TracBrowser for help on using the repository browser.