| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 /// \file eventgenerator/HepMC/HepMCEx01/src/H
27 /// \brief Implementation of the HepMCG4Pythia
28 //
29 //
30
31 #ifdef G4LIB_USE_PYTHIA
32
33 # include "HepMCG4PythiaInterface.hh"
34
35 # include "HepMC/GenEvent.h"
36 # include "HepMC/PythiaWrapper6_4.h"
37 # include "HepMCG4PythiaMessenger.hh"
38
39 //....oooOO0OOooo........oooOO0OOooo........oo
40 // additional pythia calls
41 # define pygive pygive_
42 # define pyrget pyrget_
43 # define pyrset pyrset_
44
45 extern "C"
46 {
47 void pygive(const char*, int);
48 void pyrget(int*, int*);
49 void pyrset(int*, int*);
50 }
51
52 void call_pygive(G4String s)
53 {
54 pygive(s.c_str(), s.length());
55 }
56 void call_pyrget(int a, int b)
57 {
58 pyrget(&a, &b);
59 }
60 void call_pyrset(int a, int b)
61 {
62 pyrset(&a, &b);
63 }
64
65 //....oooOO0OOooo........oooOO0OOooo........oo
66 HepMCG4PythiaInterface::HepMCG4PythiaInterface
67 {
68 # ifdef NEED_INITPYDATA
69 initpydata();
70 // Some platforms may require the initializa
71 // data as external - if you get pythia init
72 // commenting in/out the below call to initp
73 # endif
74
75 messenger = new HepMCG4PythiaMessenger(this)
76 }
77
78 //....oooOO0OOooo........oooOO0OOooo........oo
79 HepMCG4PythiaInterface::~HepMCG4PythiaInterfac
80 {
81 delete messenger;
82 }
83
84 //....oooOO0OOooo........oooOO0OOooo........oo
85 void HepMCG4PythiaInterface::CallPygive(G4Stri
86 {
87 call_pygive(par);
88 }
89
90 //....oooOO0OOooo........oooOO0OOooo........oo
91 void HepMCG4PythiaInterface::CallPyinit(G4Stri
92 G4doub
93 {
94 call_pyinit(frame.c_str(), beam.c_str(), tar
95 }
96
97 //....oooOO0OOooo........oooOO0OOooo........oo
98 void HepMCG4PythiaInterface::CallPystat(G4int
99 {
100 call_pystat(istat);
101 }
102
103 //....oooOO0OOooo........oooOO0OOooo........oo
104 void HepMCG4PythiaInterface::SetRandomSeed(G4i
105 {
106 pydatr.mrpy[1 - 1] = iseed;
107 }
108
109 //....oooOO0OOooo........oooOO0OOooo........oo
110 void HepMCG4PythiaInterface::CallPyrget(G4int
111 {
112 call_pyrget(lun, move);
113 }
114
115 //....oooOO0OOooo........oooOO0OOooo........oo
116 void HepMCG4PythiaInterface::CallPyrset(G4int
117 {
118 call_pyrset(lun, move);
119 }
120
121 //....oooOO0OOooo........oooOO0OOooo........oo
122 void HepMCG4PythiaInterface::PrintRandomStatus
123 {
124 ostr << "# Pythia random numbers status" <<
125 for (G4int j = 0; j < 6; j++) {
126 ostr << "pydatr.mrpy[" << j << "]= " << py
127 }
128 for (G4int k = 0; k < 100; k++) {
129 ostr << "pydatr.rrpy[" << k << "]= " << py
130 }
131 }
132
133 //....oooOO0OOooo........oooOO0OOooo........oo
134 void HepMCG4PythiaInterface::SetUserParameters
135 {
136 G4cout << "set user parameters of PYTHIA com
137 << G4endl;
138 }
139
140 //....oooOO0OOooo........oooOO0OOooo........oo
141 HepMC::GenEvent* HepMCG4PythiaInterface::Gener
142 {
143 static G4int nevent = 0; // event counter
144
145 call_pyevnt(); // generate one event with P
146 if (mpylist >= 1 && mpylist <= 3) call_pylis
147
148 call_pyhepc(1);
149
150 HepMC::GenEvent* evt = hepevtio.read_next_ev
151 evt->set_event_number(nevent++);
152 if (verbose > 0) evt->print();
153
154 return evt;
155 }
156
157 //....oooOO0OOooo........oooOO0OOooo........oo
158 void HepMCG4PythiaInterface::Print() const
159 {
160 G4cout << "PythiaInterface::Print()..." << G
161 }
162
163 #endif
164