ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4XAqmTotal.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4XAqmTotal.cc
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 // -------------------------------------------------------------------
28 // GEANT4 Class file
29 //
30 // For information related to this code contact:
31 //
32 // File name: G4XAqmTotal
33 //
34 // Author:
35 //
36 // Creation date: 15 April 1999
37 //
38 // Modifications:
39 //
40 // Additive Quark Model total cross section
41 // (H.J. Lipkin and F. Scheck, Phys.Rev. 16 (1966) 71
42 //
43 // -------------------------------------------------------------------
44 
45 #include "globals.hh"
46 #include "G4ios.hh"
47 #include "G4Pow.hh"
48 #include "G4SystemOfUnits.hh"
49 #include "G4XAqmTotal.hh"
50 #include "G4KineticTrack.hh"
51 #include "G4ParticleDefinition.hh"
52 
53 
54 // Validity range of this cross-section
57 
59 {
60  // As a first approximation the model is assumed to be valid over
61  // the entire energy range
62 }
63 
64 
66 { }
67 
68 
70 {
71  return (this == (G4XAqmTotal *) &right);
72 }
73 
74 
76 {
77  return (this != (G4XAqmTotal *) &right);
78 }
79 
80 
82 {
83  G4double sigma = 0.;
84 
85  // Get strangeness content
86  const G4ParticleDefinition* def1 = trk1.GetDefinition();
87  G4int sTrk1 = def1->GetQuarkContent(3) + def1->GetAntiQuarkContent(3);
88  const G4ParticleDefinition* def2 = trk2.GetDefinition();
89  G4int sTrk2 = def2->GetQuarkContent(3) + def2->GetAntiQuarkContent(3);
90 
91  // Get non-strange quark content
92  G4int qTrk1 = def1->GetQuarkContent(1) + def1->GetAntiQuarkContent(1) +
93  def1->GetQuarkContent(2) + def1->GetAntiQuarkContent(2) +
94  def1->GetQuarkContent(4) + def1->GetAntiQuarkContent(4) +
95  def1->GetQuarkContent(5) + def1->GetAntiQuarkContent(5) +
96  def1->GetQuarkContent(6) + def1->GetAntiQuarkContent(6);
97 
98  G4int qTrk2 = def2->GetQuarkContent(1) + def2->GetAntiQuarkContent(1) +
99  def2->GetQuarkContent(2) + def2->GetAntiQuarkContent(2) +
100  def2->GetQuarkContent(4) + def2->GetAntiQuarkContent(4) +
101  def2->GetQuarkContent(5) + def2->GetAntiQuarkContent(5) +
102  def2->GetQuarkContent(6) + def2->GetAntiQuarkContent(6);
103 
104  G4double sRatio1 = 0.;
105  if (qTrk1 != 0) sRatio1 = sTrk1 / qTrk1;
106 
107  G4double sRatio2 = 0.;
108  if (qTrk2 != 0) sRatio2 = sTrk2 / qTrk2;
109 
110  // Calculate the number of colliding mesons
111  G4int nMesons = 0;
112  G4int nQ1 = sTrk1 + qTrk1;
113 
114  if (nQ1 == 2) nMesons++;
115  G4int nQ2 = sTrk2 + qTrk2;
116  if (nQ2 == 2) nMesons++;
117 
118  // Cross-section (units to be checked!)
119  sigma = 40. * G4Pow::GetInstance()->powN((2.0/3.0),nMesons) * (1. - 0.4 * sRatio1) * (1. - 0.4 * sRatio2) * millibarn;
120 
121  return sigma;
122 }
123 
124 
126 {
127  G4String name("AqmTotalCrossSection");
128  return name;
129 }
130 
131 
132 
134 {
135  G4bool answer = InLimits(e,_lowLimit,_highLimit);
136 
137  return answer;
138 }