Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4DPMJET2_5Model.hh
Go to the documentation of this file.
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 // * *
21 // * Parts of this code which have been developed by QinetiQ Ltd *
22 // * under contract to the European Space Agency (ESA) are the *
23 // * intellectual property of ESA. Rights to use, copy, modify and *
24 // * redistribute this software for general public use are granted *
25 // * in compliance with any licensing, distribution and development *
26 // * policy adopted by the Geant4 Collaboration. This code has been *
27 // * written by QinetiQ Ltd for the European Space Agency, under ESA *
28 // * contract 19770/06/NL/JD (Technology Research Programme). *
29 // * *
30 // * By using, copying, modifying or distributing the software (or *
31 // * any work based on the software) you agree to acknowledge its *
32 // * use in resulting scientific publications, and indicate your *
33 // * acceptance of all terms of the Geant4 Software license. *
34 // ********************************************************************
35 //
38 //
39 #ifndef G4DPMJET2_5Model_h
40 #define G4DPMJET2_5Model_h
41 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
42 //
43 // MODULE: G4DPMJET2_5Model.hh
44 //
45 // Version: 0.B
46 // Date: 02/04/08
47 // Author: P R Truscott
48 // Organisation: QinetiQ Ltd, UK
49 // Customer: ESA/ESTEC, NOORDWIJK
50 // Contract: 19770/06/NL/JD
51 //
52 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
53 //
54 // Class Description
55 //
56 //
57 // Class Description - End
58 //
59 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
61 //
62 #include "globals.hh"
63 
64 #include "G4ParticleTable.hh"
65 
66 #include "G4HadronicInteraction.hh"
67 #include "G4WilsonAblationModel.hh"
68 #include "G4ExcitationHandler.hh"
69 #include "G4VPreCompoundModel.hh"
70 #include "G4HadFinalState.hh"
71 #include "G4Track.hh"
72 #include "G4Nucleus.hh"
73 #include "G4Fragment.hh"
74 #include "G4HadProjectile.hh"
75 #include "G4Material.hh"
76 #include "G4Element.hh"
77 
78 #include "G4GlaubAADataSetHandler.hh"
79 #include <sstream>
80 
81 class G4GlaubAADataSetHandler;
82 enum G4DPMJET2_5InitialisationType { DEFAULT=1, DPM2_5=2, DPM3=3, CORSIKA=4 };
83 
85 //
86 class G4DPMJET2_5Model : public G4HadronicInteraction
87 {
88  public:
89 //
90 //
91 // Standard constructor, destructor, copy etc declarations.
92 //
93  G4DPMJET2_5Model ();
94  G4DPMJET2_5Model (const G4DPMJET2_5InitialisationType);
95  G4DPMJET2_5Model (G4ExcitationHandler *,
96  const G4DPMJET2_5InitialisationType initType = DEFAULT);
97  G4DPMJET2_5Model (G4VPreCompoundModel *,
98  const G4DPMJET2_5InitialisationType initType = DEFAULT);
99  ~G4DPMJET2_5Model ();
100 
101  G4DPMJET2_5Model(const G4DPMJET2_5Model &right);
102 
103  const G4DPMJET2_5Model& operator=(G4DPMJET2_5Model &right);
104 
105  G4GlaubAADataSetHandler *GetGlauberDataSetHandler ();
106 
107  G4bool IsApplicable (const G4HadProjectile &theTrack, G4Nucleus &theTarget);
108 
109  virtual G4HadFinalState *ApplyYourself
110  (const G4HadProjectile &, G4Nucleus &);
111 
112  void SetVerboseLevel (G4int);
113 
114  void SetExcitationHandler (G4ExcitationHandler *);
115  void SetNoDeexcitation ();
116  void SetDefaultDeexcitation ();
117  G4ExcitationHandler *GetExcitationHandler () const;
118 
119  void SetPreCompoundModel(G4VPreCompoundModel* value);
120  void SetNoPreCompoundModel ();
121  void SetDefaultPreCompoundModel ();
122  G4VPreCompoundModel* GetPreCompoundModel() const;
123 
124  void SetDPMInitialRandomSeeds (const G4int seed1, const G4int seed2);
125  G4int GetDPMInitialRandomSeeds (const G4int i) const;
126 
127  G4double GetMinEnergy( const G4Material *aMaterial,
128  const G4Element *anElement ) const;
129  G4double GetMaxEnergy( const G4Material *aMaterial,
130  const G4Element *anElement ) const;
131 
132  G4bool SetVerboseFortranOutput (const G4String filename);
133  G4String GetVerboseFortranOutput () const;
134 
135  void SetDPMVariablesTAUFOR (const G4double TAUFOR_P,
136  const G4int KTAUGE_P,
137  const G4int ITAUVE_P);
138 /* void SetDPMVariablesXCUTS (const G4double CVQ_P,
139  const G4double CDQ_P,
140  const G4double CSEA_P,
141  const G4double SSMIMA_P);
142  void SetDPMVariablesCRONINPT (const G4int MKCRON_P,
143  const G4double CRONCO_P);
144  void SetDPMVariablesSEADISTR (const G4double
145  const G4double UNON_P,
146  const G4double UNOM_P,
147  const G4double UNOSEA_P);*/
148 
149  private:
150  void DumpVerboseInformation1 (const G4int n) const;
151  void DumpVerboseInformation2 (const G4String particleName,
152  const G4ThreeVector p, const G4double E, const G4double T,
153  const G4ThreeVector pinit) const;
154  void DumpVerboseInformation3 (const G4int i, const G4int A, const G4int Z,
155  const G4ThreeVector p, const G4double E, const G4double T,
156  const G4ThreeVector pinit) const;
157  void DumpVerboseInformation4 (const G4int i, const G4String particleName,
158  const G4ThreeVector p, const G4double E, const G4double T,
159  const G4ThreeVector pinit) const;
160  void PrintWelcomeMessage () const;
161  void Initialise ();
162 
163  private:
164  G4DPMJET2_5InitialisationType
165  theInitType;
166 
167  G4GlaubAADataSetHandler *theGlauberDataSetHandler;
168 
169  G4ParticleTable *theParticleTable;
170  G4IonTable *theIonTable;
171 
172  G4bool debug;
173  G4int debug_level;
174  G4int lunber;
175  G4double dpmver;
176  G4String defaultDirName;
177 
178  G4ExcitationHandler *theExcitationHandler;
179  G4VPreCompoundModel *thePreComp;
180 
181  G4double TAUFOR;
182  G4int KTAUGE;
183  G4int ITAUVE;
184  G4double UNON;
185  G4double UNOM;
186  G4double UNOSEA;
187  G4double CVQ;
188  G4double CDQ;
189  G4double CSEA;
190  G4double SSMIMA;
191  G4double VVMTHR;
192  G4double SEASQ;
193  G4int MKCRON;
194  G4double CRONCO;
195  G4int ISINGD;
196  G4int ISINGX;
197  G4int IDUBLD;
198  G4double SDFRAC;
199 
200 
201  ftnlogical LTRUE;
202  ftnlogical LFALSE;
203 
204  G4String verboseFortranFile;
205 };
207 //
208 inline G4GlaubAADataSetHandler *G4DPMJET2_5Model::GetGlauberDataSetHandler ()
209  {return theGlauberDataSetHandler;}
211 //
212 inline void G4DPMJET2_5Model::SetExcitationHandler
213  (G4ExcitationHandler *aExcitationHandler)
214  {theExcitationHandler = aExcitationHandler;}
216 //
217 inline G4ExcitationHandler *G4DPMJET2_5Model::GetExcitationHandler () const
218  {return theExcitationHandler;}
220 //
221 inline void G4DPMJET2_5Model::SetPreCompoundModel
222  (G4VPreCompoundModel *aPreCompoundModel)
223  {thePreComp = aPreCompoundModel;}
225 //
226 inline G4VPreCompoundModel *G4DPMJET2_5Model::GetPreCompoundModel () const
227  {return thePreComp;}
229 //
230 inline void G4DPMJET2_5Model::SetVerboseLevel (G4int verboseLevel1)
231  {verboseLevel = verboseLevel1;}
233 //
234 inline G4double G4DPMJET2_5Model::GetMinEnergy( const G4Material *,
235  const G4Element * ) const
236  {return theMinEnergy;}
238 //
239 inline G4double G4DPMJET2_5Model::GetMaxEnergy( const G4Material *,
240  const G4Element * ) const
241  {return theMaxEnergy;}
243 //
244 // SetVerboseFortranOutput
245 //
246 inline G4String G4DPMJET2_5Model::GetVerboseFortranOutput () const
247  {return verboseFortranFile;}
249 //
250 inline void G4DPMJET2_5Model::DumpVerboseInformation2
251  (const G4String particleName, const G4ThreeVector p,
252  const G4double E, const G4double T, const G4ThreeVector pinit) const
253 {
254  G4cout <<"Name = " <<particleName <<G4endl;
255  G4cout <<" Momentum = " <<p/MeV <<" MeV/c" <<G4endl;
256  G4cout <<" T. Energy = " <<E/MeV <<" MeV" <<G4endl;
257  G4cout <<" K. Energy = " <<T/MeV <<" MeV" <<G4endl;
258  if (verboseLevel >= 3)
259  {
260  G4ThreeVector axis = pinit.unit();
261  G4double pz = p.dot(axis);
262  G4cout <<" Transverse mass = " <<std::sqrt(E*E-pz*pz)/MeV <<" MeV"
263  <<G4endl;
264  G4cout <<" Rapidity = "
265  <<0.5*std::log((E+pz)/(E-pz)) <<G4endl;
266  }
267 }
269 //
270 inline void G4DPMJET2_5Model::DumpVerboseInformation3 (const G4int i,
271  const G4int A, const G4int Z, const G4ThreeVector p,
272  const G4double E, const G4double T, const G4ThreeVector pinit) const
273 {
274  G4cout <<"----------------------------------------"
275  <<"----------------------------------------" <<G4endl;
276  G4cout <<"The nuclear fragment #" <<i <<" before" <<G4endl;
277  G4cout <<"----------------------------------------"
278  <<"----------------------------------------" <<G4endl;
279 
280  std::ostringstream tmpStream;
281  tmpStream <<"(A = " <<A <<", Z = " <<Z <<")";
282  G4String AZ = tmpStream.str();
283 
284  DumpVerboseInformation2(AZ, p, E, T, pinit);
285 }
287 //
288 inline void G4DPMJET2_5Model::DumpVerboseInformation4 (const G4int i,
289  const G4String particleName, const G4ThreeVector p,
290  const G4double E, const G4double T, const G4ThreeVector pinit) const
291 {
292  G4cout <<"----------------------------------------"
293  <<"----------------------------------------" <<G4endl;
294  G4cout <<"The nuclear fragment #" <<i <<" after" <<G4endl;
295  G4cout <<"----------------------------------------"
296  <<"----------------------------------------" <<G4endl;
297 
298  DumpVerboseInformation2(particleName, p, E, T, pinit);
299 }
301 //
302 inline void G4DPMJET2_5Model::SetDPMInitialRandomSeeds (const G4int seed1,
303  const G4int seed2)
304 {
305  G4int s1 = seed1;
306  G4int s2 = seed2;
307  rd2in_ (&s1, &s2);
308 }
310 //
311 inline G4int G4DPMJET2_5Model::GetDPMInitialRandomSeeds (const G4int i)
312  const
313 {
314  G4int seed1, seed2;
315  rd2out_ (&seed1, &seed2);
316  if (i == 0) return seed1;
317  else if (i == 2) return seed2;
318  else return 0;
319 }
321 //
322 inline void G4DPMJET2_5Model::SetDPMVariablesTAUFOR (const G4double TAUFOR_P,
323  const G4int KTAUGE_P, const G4int ITAUVE_P)
324 {
325  TAUFOR = TAUFOR_P;
326  KTAUGE = KTAUGE_P;
327  ITAUVE = ITAUVE_P;
328  taufo_.taufor = TAUFOR;
329  taufo_.ktauge = KTAUGE;
330  taufo_.itauve = ITAUVE;
331 }
333 //
334 #endif
Generated on Sat May 25 2013 14:32:56 for Geant4 by   doxygen 1.8.4