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 // ==================================================================== 27 /// \brief Implementation of the HepMCG4Pythia << 28 // 27 // >> 28 // HepMCG4PythiaInterface.cc >> 29 // $Id: HepMCG4PythiaInterface.cc,v 1.6 2006/07/05 12:04:13 gcosmo Exp $ 29 // 30 // >> 31 // ==================================================================== 30 32 31 #ifdef G4LIB_USE_PYTHIA 33 #ifdef G4LIB_USE_PYTHIA 32 34 33 # include "HepMCG4PythiaInterface.hh" << 35 #include "HepMCG4PythiaInterface.hh" >> 36 #include "HepMCG4PythiaMessenger.hh" 34 37 35 # include "HepMC/GenEvent.h" << 38 #include "HepMC/GenEvent.h" 36 # include "HepMC/PythiaWrapper6_4.h" << 39 #include "HepMC/PythiaWrapper6_152.h" 37 # include "HepMCG4PythiaMessenger.hh" << 38 40 39 //....oooOO0OOooo........oooOO0OOooo........oo << 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 } << 55 56 void call_pyrget(int a, int b) << 56 //////////////////////////////////////////////// 57 { << 57 HepMCG4PythiaInterface::HepMCG4PythiaInterface() 58 pyrget(&a, &b); << 58 : verbose(0), mpylist(0) 59 } << 59 //////////////////////////////////////////////// 60 void call_pyrset(int a, int b) << 61 { 60 { 62 pyrset(&a, &b); << 61 #ifdef NEED_INITPYDATA 63 } << 64 << 65 //....oooOO0OOooo........oooOO0OOooo........oo << 66 HepMCG4PythiaInterface::HepMCG4PythiaInterface << 67 { << 68 # ifdef NEED_INITPYDATA << 69 initpydata(); 62 initpydata(); 70 // Some platforms may require the initializa << 63 // Some platforms may require the initialization of pythia PYDATA block 71 // data as external - if you get pythia init << 64 // data as external - if you get pythia initialization errors try 72 // commenting in/out the below call to initp 65 // commenting in/out the below call to initpydata(). 73 # endif << 66 #endif 74 67 75 messenger = new HepMCG4PythiaMessenger(this) << 68 messenger= new HepMCG4PythiaMessenger(this); 76 } 69 } 77 70 78 //....oooOO0OOooo........oooOO0OOooo........oo << 71 ///////////////////////////////////////////////// 79 HepMCG4PythiaInterface::~HepMCG4PythiaInterfac 72 HepMCG4PythiaInterface::~HepMCG4PythiaInterface() >> 73 ///////////////////////////////////////////////// 80 { 74 { 81 delete messenger; 75 delete messenger; 82 } 76 } 83 77 84 //....oooOO0OOooo........oooOO0OOooo........oo << 78 ///////////////////////////////////////////////////// 85 void HepMCG4PythiaInterface::CallPygive(G4Stri 79 void HepMCG4PythiaInterface::CallPygive(G4String par) >> 80 ///////////////////////////////////////////////////// 86 { 81 { 87 call_pygive(par); 82 call_pygive(par); 88 } 83 } 89 84 90 //....oooOO0OOooo........oooOO0OOooo........oo << 85 ////////////////////////////////////////////////////////////////////// 91 void HepMCG4PythiaInterface::CallPyinit(G4Stri << 86 void HepMCG4PythiaInterface::CallPyinit(G4String frame, G4String beam, 92 G4doub << 87 G4String target, G4double win) >> 88 ////////////////////////////////////////////////////////////////////// 93 { 89 { 94 call_pyinit(frame.c_str(), beam.c_str(), tar 90 call_pyinit(frame.c_str(), beam.c_str(), target.c_str(), win); 95 } 91 } 96 92 97 //....oooOO0OOooo........oooOO0OOooo........oo << 93 //////////////////////////////////////////////////// 98 void HepMCG4PythiaInterface::CallPystat(G4int 94 void HepMCG4PythiaInterface::CallPystat(G4int istat) >> 95 //////////////////////////////////////////////////// 99 { 96 { 100 call_pystat(istat); 97 call_pystat(istat); 101 } 98 } 102 99 103 //....oooOO0OOooo........oooOO0OOooo........oo << 100 /////////////////////////////////////////////////////// 104 void HepMCG4PythiaInterface::SetRandomSeed(G4i 101 void HepMCG4PythiaInterface::SetRandomSeed(G4int iseed) >> 102 /////////////////////////////////////////////////////// 105 { 103 { 106 pydatr.mrpy[1 - 1] = iseed; << 104 pydatr.mrpy[1-1]= iseed; 107 } 105 } 108 106 109 //....oooOO0OOooo........oooOO0OOooo........oo << 107 ////////////////////////////////////////////////////////////// 110 void HepMCG4PythiaInterface::CallPyrget(G4int 108 void HepMCG4PythiaInterface::CallPyrget(G4int lun, G4int move) >> 109 ////////////////////////////////////////////////////////////// 111 { 110 { 112 call_pyrget(lun, move); 111 call_pyrget(lun, move); 113 } 112 } 114 113 115 //....oooOO0OOooo........oooOO0OOooo........oo << 114 ////////////////////////////////////////////////////////////// 116 void HepMCG4PythiaInterface::CallPyrset(G4int 115 void HepMCG4PythiaInterface::CallPyrset(G4int lun, G4int move) >> 116 ////////////////////////////////////////////////////////////// 117 { 117 { 118 call_pyrset(lun, move); 118 call_pyrset(lun, move); 119 } 119 } 120 120 121 //....oooOO0OOooo........oooOO0OOooo........oo << 121 ////////////////////////////////////////////////////////////////////////// 122 void HepMCG4PythiaInterface::PrintRandomStatus 122 void HepMCG4PythiaInterface::PrintRandomStatus(std::ostream& ostr) const >> 123 ////////////////////////////////////////////////////////////////////////// 123 { 124 { 124 ostr << "# Pythia random numbers status" << 125 ostr << "# Pythia random numbers status" << G4endl; 125 for (G4int j = 0; j < 6; j++) { << 126 for (G4int j=0; j<6; j++) { 126 ostr << "pydatr.mrpy[" << j << "]= " << py 127 ostr << "pydatr.mrpy[" << j << "]= " << pydatr.mrpy[j] << G4endl; 127 } 128 } 128 for (G4int k = 0; k < 100; k++) { << 129 for (G4int k=0; k<100; k++) { 129 ostr << "pydatr.rrpy[" << k << "]= " << py 130 ostr << "pydatr.rrpy[" << k << "]= " << pydatr.rrpy[k] << G4endl; 130 } 131 } 131 } 132 } 132 133 133 //....oooOO0OOooo........oooOO0OOooo........oo << 134 //////////////////////////////////////////////// 134 void HepMCG4PythiaInterface::SetUserParameters 135 void HepMCG4PythiaInterface::SetUserParameters() >> 136 //////////////////////////////////////////////// 135 { 137 { 136 G4cout << "set user parameters of PYTHIA com << 138 G4cout << "set user parameters of PYTHIA common." << G4endl 137 << G4endl; << 139 << "nothing to be done in default." >> 140 << G4endl; 138 } 141 } 139 142 140 //....oooOO0OOooo........oooOO0OOooo........oo << 143 ///////////////////////////////////////////////////////////// 141 HepMC::GenEvent* HepMCG4PythiaInterface::Gener 144 HepMC::GenEvent* HepMCG4PythiaInterface::GenerateHepMCEvent() >> 145 ///////////////////////////////////////////////////////////// 142 { 146 { 143 static G4int nevent = 0; // event counter << 147 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 >> 149 call_pyevnt(); // generate one event with Pythia >> 150 if(mpylist >=1 && mpylist<= 3) call_pylist(mpylist); >> 151 148 call_pyhepc(1); 152 call_pyhepc(1); 149 153 150 HepMC::GenEvent* evt = hepevtio.read_next_ev << 154 HepMC::GenEvent* evt= hepevtio.read_next_event(); 151 evt->set_event_number(nevent++); << 155 evt-> set_event_number(nevent++); 152 if (verbose > 0) evt->print(); << 156 if(verbose>0) evt-> print(); 153 157 154 return evt; 158 return evt; 155 } 159 } 156 160 157 //....oooOO0OOooo........oooOO0OOooo........oo << 161 ////////////////////////////////////////// 158 void HepMCG4PythiaInterface::Print() const 162 void HepMCG4PythiaInterface::Print() const >> 163 ////////////////////////////////////////// 159 { 164 { 160 G4cout << "PythiaInterface::Print()..." << G 165 G4cout << "PythiaInterface::Print()..." << G4endl; 161 } 166 } 162 167 163 #endif 168 #endif 164 169