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