[817] | 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 | * |
---|
| 27 | * $Id: g3tog4.F,v 1.5 2006/06/29 18:15:12 gunter Exp $ |
---|
[965] | 28 | * GEANT4 tag $Name: geant4-09-02-ref-02 $ |
---|
[817] | 29 | * |
---|
| 30 | * |
---|
| 31 | * G3toG4 |
---|
| 32 | * |
---|
| 33 | * Package to convert Geant3 Fortran geometry code to a call list |
---|
| 34 | * file to be interpreted by the Geant4 geometry call list |
---|
| 35 | * interpreter, or alternatively, directly to Geant4 code. |
---|
| 36 | * |
---|
| 37 | * This set of routines is to be linked in front of, so overriding, |
---|
| 38 | * the standard Geant library. |
---|
| 39 | * |
---|
| 40 | * It is possible to execute the Geant calls while at the same |
---|
| 41 | * time building the call list/Geant4 code. In order to do this, |
---|
| 42 | * these routines must occupy a different name space to that of |
---|
| 43 | * the real Geant routines. This is provided by the CALL_GEANT |
---|
| 44 | * cpp flag. If these routines are compiled with this flag, |
---|
| 45 | * the routine names begin with K rather than G. eg. GSVOLU |
---|
| 46 | * becomes KSVOLU. Routine names in your source code must be so |
---|
| 47 | * converted; a perl script is provided to do this. |
---|
| 48 | * $$$ provide the script |
---|
| 49 | * Under normal circumstances it should *not* be necessary to go |
---|
| 50 | * through this; it is only necessary if during the geometry |
---|
| 51 | * generation process your code extracts information from Geant |
---|
| 52 | * about material already generated. |
---|
| 53 | * |
---|
| 54 | * Torre Wenaus, LLNL 6/95 |
---|
| 55 | * |
---|
| 56 | * To Do |
---|
| 57 | * - option to divide generated Geant4 code into separate files/routines |
---|
| 58 | * based on context |
---|
| 59 | * |
---|
| 60 | ************************************************************************ |
---|
| 61 | * |
---|
| 62 | subroutine G3toG4(luni,lunc,chopt) |
---|
| 63 | ************************************************************************ |
---|
| 64 | * |
---|
| 65 | * G3toG4 |
---|
| 66 | * |
---|
| 67 | * Initialization/setup routine |
---|
| 68 | * |
---|
| 69 | * luni (call list), lunc (C++ code) logical unit numbers: |
---|
| 70 | * lun>0: Open output file on unit lun. Filenames used: |
---|
| 71 | * g3calls.dat Call list file |
---|
| 72 | * g4geom.cc Geant4 C++ geometry code |
---|
| 73 | * lun<0: File open has been done by the user. Just write to |lun| |
---|
| 74 | * lun=0: Don't generate this output. |
---|
| 75 | * ie. luni=0: Don't generate the call list |
---|
| 76 | * lunc=0: Don't generate the Geant4 code |
---|
| 77 | * |
---|
| 78 | * chopt options: |
---|
| 79 | * 'G' execute the actual Geant calls as well as building the |
---|
| 80 | * code/call list. In case users use info obtained from Geant |
---|
| 81 | * during the geometry building process. THIS IS THE DEFAULT |
---|
| 82 | * at present: |
---|
| 83 | #define CALL_GEANT |
---|
| 84 | * |
---|
| 85 | ************************************************************************ |
---|
| 86 | implicit none |
---|
| 87 | integer luni, lunc |
---|
| 88 | character chopt*(*) |
---|
| 89 | #include "G3toG4.inc" |
---|
| 90 | * |
---|
| 91 | print *,'Initializing Geant3 to Geant4 conversion' |
---|
| 92 | #ifdef CALL_GEANT |
---|
| 93 | c dogeom = index(chopt,'G') + index(chopt,'g') .ne. 0 |
---|
| 94 | dogeom = .true. |
---|
| 95 | #else |
---|
| 96 | dogeom = .false. |
---|
| 97 | #endif |
---|
| 98 | context = '----' |
---|
| 99 | if (luni.eq.0.and.lunc.eq.0) then |
---|
| 100 | print *,'G3TOG4: No output requested by user. No output'// |
---|
| 101 | + ' will be generated.' |
---|
| 102 | endif |
---|
| 103 | lunlist = abs(luni) |
---|
| 104 | luncode = abs(lunc) |
---|
| 105 | if (lunlist.ne.0) then |
---|
| 106 | doclist = .true. |
---|
| 107 | else |
---|
| 108 | doclist = .false. |
---|
| 109 | endif |
---|
| 110 | if (luncode.ne.0) then |
---|
| 111 | docode = .true. |
---|
| 112 | else |
---|
| 113 | docode = .false. |
---|
| 114 | endif |
---|
| 115 | *** If lun>0, open the file |
---|
| 116 | if (lunlist.gt.0) then |
---|
| 117 | open(unit=lunlist,file='g3calls.dat',status='unknown') |
---|
| 118 | endif |
---|
| 119 | if (luncode.gt.0) then |
---|
| 120 | nfile = 1 |
---|
| 121 | call g3source |
---|
| 122 | endif |
---|
| 123 | * |
---|
| 124 | end |
---|
| 125 | * |
---|
| 126 | subroutine g4init |
---|
| 127 | ************************************************************************ |
---|
| 128 | ************************************************************************ |
---|
| 129 | implicit none |
---|
| 130 | #include "G3toG4.inc" |
---|
| 131 | * |
---|
| 132 | if (luncode.ne.0) then |
---|
| 133 | write(luncode, |
---|
| 134 | + '(''//G4GeometryManager* GeoMgr = new G4GeometryManager();'')') |
---|
| 135 | * call ctocp('void G3G4init();') |
---|
| 136 | endif |
---|
| 137 | * |
---|
| 138 | end |
---|
| 139 | * |
---|
| 140 | subroutine g3header |
---|
| 141 | ************************************************************************ |
---|
| 142 | * |
---|
| 143 | ************************************************************************ |
---|
| 144 | implicit none |
---|
| 145 | call g4init |
---|
| 146 | end |
---|
| 147 | |
---|
| 148 | subroutine g3source |
---|
| 149 | ************************************************************************ |
---|
| 150 | * |
---|
| 151 | ************************************************************************ |
---|
| 152 | implicit none |
---|
| 153 | #include "G3toG4.inc" |
---|
| 154 | character fname*30 |
---|
| 155 | if (luncode.le.0) return |
---|
| 156 | if (nfile.gt.1) write(luncode,'(''}'')') |
---|
| 157 | close(luncode) |
---|
| 158 | write (fname,'(''G3toG4code_'',i2.2,''.cc'')') nfile |
---|
| 159 | open(unit=luncode,file=fname,status='unknown') |
---|
| 160 | write(luncode,'(''#include "G3toG4.hh"'')') |
---|
| 161 | if (nfile.eq.1) call g3header |
---|
| 162 | write(luncode,'(/''void G3toG4code_'',i2.2,''()'')') nfile |
---|
| 163 | write(luncode,'(''{'')') |
---|
| 164 | call ctocp('// init to 0 avoids "unused" warnings') |
---|
| 165 | call ctocp('G4int nd=0,nh=0,nv=0,imate=0,itmed=0,nmat=0,') |
---|
| 166 | call ctocp(' isvol=0,ifield=0,nwhi=0,nwdi=0,idtyp=0,ipart=0,') |
---|
| 167 | call ctocp(' itrtyp=0,nlmat=0,npar=0,ndvmx=0,numed=0,iaxis=0,') |
---|
| 168 | call ctocp( |
---|
| 169 | + ' ndiv=0,irot=0,ival=0,num=0,nmed=0,nbits[100],mode[6];') |
---|
| 170 | call ctocp('G4String chnam[100];') |
---|
| 171 | call ctocp('G4String name="",moth="",attr="",only="",shape="";') |
---|
| 172 | call ctocp('G4String chset="",chdet="",chali="",chpar="";') |
---|
| 173 | call ctocp('G4double amass=0.,charge=0.,tlife=0.,parval=0.;') |
---|
| 174 | call ctocp('G4double c0=0.,step=0.,a=0.,dens=0.,radl=0.,x=0.;') |
---|
| 175 | call ctocp('G4double y=0.,z=0.,theta1=0.,phi1=0.,theta2=0.;') |
---|
| 176 | call ctocp('G4double phi2=0.,theta3=0.,phi3=0.,fieldm=0.;') |
---|
| 177 | call ctocp('G4double tmaxfd=0.,stemax=0.,deemax=0.,epsil=0.;') |
---|
| 178 | call ctocp('G4double stmin=0.,par[100],fact[100],orig[100];') |
---|
| 179 | call ctocp('G4double bratio[6],aa[100],zz[100],wmat[100];') |
---|
| 180 | call ctocp('nbits[0]=mode[0]=0;chnam[0]="";par[0]=0.;') |
---|
| 181 | call ctocp('fact[0]=orig[0]=bratio[0]=aa[0]=zz[0]=wmat[0]=0.;') |
---|
| 182 | call ctocp(' ') |
---|
| 183 | if (nfile.eq.1) then |
---|
| 184 | * call ctocp('G3G4init();') |
---|
| 185 | call ctocp(' ') |
---|
| 186 | endif |
---|
| 187 | end |
---|
| 188 | |
---|
| 189 | subroutine g3main |
---|
| 190 | ************************************************************************ |
---|
| 191 | ************************************************************************ |
---|
| 192 | implicit none |
---|
| 193 | #include "G3toG4.inc" |
---|
| 194 | integer i |
---|
| 195 | * |
---|
| 196 | close(luncode) |
---|
| 197 | open(unit=luncode,file='G3toG4code.cc',status='unknown') |
---|
| 198 | do i=1,nfile |
---|
| 199 | write(luncode,'('' void G3toG4code_'',i2.2,''();'')') i |
---|
| 200 | enddo |
---|
| 201 | call ctocp('void G3toG4code()') |
---|
| 202 | call ctocp('{') |
---|
| 203 | do i=1,nfile |
---|
| 204 | write(luncode,'('' G3toG4code_'',i2.2,''();'')') i |
---|
| 205 | enddo |
---|
| 206 | call ctocp('}') |
---|
| 207 | close(luncode) |
---|
| 208 | end |
---|
| 209 | |
---|
| 210 | subroutine g3context(cntxt) |
---|
| 211 | ************************************************************************ |
---|
| 212 | * |
---|
| 213 | * g3context |
---|
| 214 | * |
---|
| 215 | * Set the current geometry code context. eg. context can be used |
---|
| 216 | * to distinguish code for different subdetectors. The Geant4 |
---|
| 217 | * call list interpreter can then execute the code selectively for |
---|
| 218 | * a particular context only, if desired. Spaces not allowed. |
---|
| 219 | * |
---|
| 220 | ************************************************************************ |
---|
| 221 | implicit none |
---|
| 222 | character*(*) cntxt |
---|
| 223 | #include "G3toG4.inc" |
---|
| 224 | context = cntxt |
---|
| 225 | end |
---|
| 226 | * |
---|
| 227 | subroutine ctocp(string) |
---|
| 228 | ************************************************************************ |
---|
| 229 | ************************************************************************ |
---|
| 230 | implicit none |
---|
| 231 | character*(*) string |
---|
| 232 | #include "G3toG4.inc" |
---|
| 233 | write (luncode,*) string |
---|
| 234 | end |
---|
| 235 | * |
---|
| 236 | subroutine rtocp(string,x) |
---|
| 237 | ************************************************************************ |
---|
| 238 | ************************************************************************ |
---|
| 239 | implicit none |
---|
| 240 | character*(*) string |
---|
| 241 | real x |
---|
| 242 | #include "G3toG4.inc" |
---|
| 243 | write(luncode,'(4x,a,'' = '',e14.8,'';'')') |
---|
| 244 | + string, x |
---|
| 245 | end |
---|
| 246 | * |
---|
| 247 | subroutine artocp(string,ax,n) |
---|
| 248 | ************************************************************************ |
---|
| 249 | ************************************************************************ |
---|
| 250 | implicit none |
---|
| 251 | character*(*) string |
---|
| 252 | real ax(*) |
---|
| 253 | integer n,i |
---|
| 254 | #include "G3toG4.inc" |
---|
| 255 | do i=1,n |
---|
| 256 | write(luncode,'('' '',a,''['',i3,''] = '',e14.8,'';'')') |
---|
| 257 | + string, i-1, ax(i) |
---|
| 258 | enddo |
---|
| 259 | end |
---|
| 260 | * |
---|
| 261 | subroutine aitocp(string,ai,n) |
---|
| 262 | ************************************************************************ |
---|
| 263 | ************************************************************************ |
---|
| 264 | implicit none |
---|
| 265 | character*(*) string |
---|
| 266 | integer ai(*) |
---|
| 267 | integer n,i |
---|
| 268 | #include "G3toG4.inc" |
---|
| 269 | do i=1,n |
---|
| 270 | write(luncode,'('' '',a,''['',i3,''] = '',i10,'';'')') |
---|
| 271 | + string, i-1, ai(i) |
---|
| 272 | enddo |
---|
| 273 | end |
---|
| 274 | * |
---|
| 275 | subroutine astocp(string,ac,n) |
---|
| 276 | ************************************************************************ |
---|
| 277 | ************************************************************************ |
---|
| 278 | implicit none |
---|
| 279 | character*(*) string, ac(*) |
---|
| 280 | integer n,i |
---|
| 281 | #include "G3toG4.inc" |
---|
| 282 | c write(luncode,'('' G4String '',a,''['',i3,''];'')') string, n |
---|
| 283 | do i=1,n |
---|
| 284 | write(luncode,'('' '',a,''['',i3,''] = "'',a,''";'')') |
---|
| 285 | + string, i-1, ac(i) |
---|
| 286 | enddo |
---|
| 287 | end |
---|
| 288 | * |
---|
| 289 | subroutine g3ldpar(par,npar) |
---|
| 290 | ************************************************************************ |
---|
| 291 | * |
---|
| 292 | * g3ldpar |
---|
| 293 | * |
---|
| 294 | ************************************************************************ |
---|
| 295 | implicit none |
---|
| 296 | * |
---|
| 297 | integer npar, i |
---|
| 298 | real par(*) |
---|
| 299 | #include "G3toG4.inc" |
---|
| 300 | * |
---|
| 301 | if (npar.gt.0) then |
---|
| 302 | write(luncode,'('' par['',i4,''] = '',e14.8,'';'')') |
---|
| 303 | + (i-1,par(i),i=1,npar) |
---|
| 304 | endif |
---|
| 305 | end |
---|
| 306 | * |
---|
| 307 | subroutine check_lines |
---|
| 308 | ************************************************************************ |
---|
| 309 | ************************************************************************ |
---|
| 310 | implicit none |
---|
| 311 | #include "G3toG4.inc" |
---|
| 312 | if (luncode.ne.0) then |
---|
| 313 | nlines = nlines +1 |
---|
| 314 | if (nlines.gt.maxlines) then |
---|
| 315 | nfile = nfile +1 |
---|
| 316 | call g3source |
---|
| 317 | nlines = 0 |
---|
| 318 | endif |
---|
| 319 | endif |
---|
| 320 | end |
---|