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