| Geant4 Cross Reference |
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 /// \file eventgenerator/HepMC/HepMCEx01/src/HepMCG4PythiaInterface.cc
27 /// \brief Implementation of the HepMCG4PythiaInterface class
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........oooOO0OOooo........oooOO0OOooo......
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........oooOO0OOooo........oooOO0OOooo......
66 HepMCG4PythiaInterface::HepMCG4PythiaInterface() : verbose(0), mpylist(0)
67 {
68 # ifdef NEED_INITPYDATA
69 initpydata();
70 // Some platforms may require the initialization of pythia PYDATA block
71 // data as external - if you get pythia initialization errors try
72 // commenting in/out the below call to initpydata().
73 # endif
74
75 messenger = new HepMCG4PythiaMessenger(this);
76 }
77
78 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
79 HepMCG4PythiaInterface::~HepMCG4PythiaInterface()
80 {
81 delete messenger;
82 }
83
84 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
85 void HepMCG4PythiaInterface::CallPygive(G4String par)
86 {
87 call_pygive(par);
88 }
89
90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
91 void HepMCG4PythiaInterface::CallPyinit(G4String frame, G4String beam, G4String target,
92 G4double win)
93 {
94 call_pyinit(frame.c_str(), beam.c_str(), target.c_str(), win);
95 }
96
97 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
98 void HepMCG4PythiaInterface::CallPystat(G4int istat)
99 {
100 call_pystat(istat);
101 }
102
103 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
104 void HepMCG4PythiaInterface::SetRandomSeed(G4int iseed)
105 {
106 pydatr.mrpy[1 - 1] = iseed;
107 }
108
109 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
110 void HepMCG4PythiaInterface::CallPyrget(G4int lun, G4int move)
111 {
112 call_pyrget(lun, move);
113 }
114
115 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
116 void HepMCG4PythiaInterface::CallPyrset(G4int lun, G4int move)
117 {
118 call_pyrset(lun, move);
119 }
120
121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
122 void HepMCG4PythiaInterface::PrintRandomStatus(std::ostream& ostr) const
123 {
124 ostr << "# Pythia random numbers status" << G4endl;
125 for (G4int j = 0; j < 6; j++) {
126 ostr << "pydatr.mrpy[" << j << "]= " << pydatr.mrpy[j] << G4endl;
127 }
128 for (G4int k = 0; k < 100; k++) {
129 ostr << "pydatr.rrpy[" << k << "]= " << pydatr.rrpy[k] << G4endl;
130 }
131 }
132
133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
134 void HepMCG4PythiaInterface::SetUserParameters()
135 {
136 G4cout << "set user parameters of PYTHIA common." << G4endl << "nothing to be done in default."
137 << G4endl;
138 }
139
140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
141 HepMC::GenEvent* HepMCG4PythiaInterface::GenerateHepMCEvent()
142 {
143 static G4int nevent = 0; // event counter
144
145 call_pyevnt(); // generate one event with Pythia
146 if (mpylist >= 1 && mpylist <= 3) call_pylist(mpylist);
147
148 call_pyhepc(1);
149
150 HepMC::GenEvent* evt = hepevtio.read_next_event();
151 evt->set_event_number(nevent++);
152 if (verbose > 0) evt->print();
153
154 return evt;
155 }
156
157 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
158 void HepMCG4PythiaInterface::Print() const
159 {
160 G4cout << "PythiaInterface::Print()..." << G4endl;
161 }
162
163 #endif
164