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 // 26 // >> 27 // $Id: G4gstmed.cc,v 1.7 2006-06-29 18:15:02 gunter Exp $ >> 28 // GEANT4 tag $Name: not supported by cvs2svn $ 27 // 29 // 28 // The last G4int argument of G4gstmed(..) is 30 // The last G4int argument of G4gstmed(..) is used for sending 29 // info whether the Geant3 tracking medium par 31 // info whether the Geant3 tracking medium parameters should 30 // be set (the max step, later: G3 default cut 32 // be set (the max step, later: G3 default cut values). 31 // 33 // 32 // by I.Hrivnacova, 27 Sep 99 34 // by I.Hrivnacova, 27 Sep 99 33 35 34 #include "G4SystemOfUnits.hh" << 35 #include "G4LogicalVolume.hh" 36 #include "G4LogicalVolume.hh" 36 #include "G3toG4.hh" 37 #include "G3toG4.hh" 37 #include "G3MatTable.hh" 38 #include "G3MatTable.hh" 38 #include "G3MedTable.hh" 39 #include "G3MedTable.hh" 39 #include "G4UserLimits.hh" 40 #include "G4UserLimits.hh" 40 #include "G4MagneticField.hh" 41 #include "G4MagneticField.hh" 41 #include "G4Material.hh" 42 #include "G4Material.hh" 42 43 43 void PG4gstmed(G4String *tokens) 44 void PG4gstmed(G4String *tokens) 44 { 45 { 45 // fill the parameter containers 46 // fill the parameter containers 46 G3fillParams(tokens,PTgstmed); 47 G3fillParams(tokens,PTgstmed); 47 48 48 // interpret the parameters 49 // interpret the parameters 49 G4String name = Spar[0]; 50 G4String name = Spar[0]; 50 G4int itmed = Ipar[0]; 51 G4int itmed = Ipar[0]; 51 G4int nmat = Ipar[1]; 52 G4int nmat = Ipar[1]; 52 G4int isvol = Ipar[2]; 53 G4int isvol = Ipar[2]; 53 G4int ifield = Ipar[3]; 54 G4int ifield = Ipar[3]; 54 G4int nwbuf = Ipar[4]; 55 G4int nwbuf = Ipar[4]; 55 G4double fieldm = Rpar[0]; 56 G4double fieldm = Rpar[0]; 56 G4double tmaxfd = Rpar[1]; 57 G4double tmaxfd = Rpar[1]; 57 G4double stemax = Rpar[2]; 58 G4double stemax = Rpar[2]; 58 G4double deemax = Rpar[3]; 59 G4double deemax = Rpar[3]; 59 G4double epsil = Rpar[4]; 60 G4double epsil = Rpar[4]; 60 G4double stmin = Rpar[5]; 61 G4double stmin = Rpar[5]; 61 G4double *ubuf = &Rpar[6]; 62 G4double *ubuf = &Rpar[6]; 62 63 63 G4gstmed(itmed,name,nmat,isvol,ifield,fiel 64 G4gstmed(itmed,name,nmat,isvol,ifield,fieldm,tmaxfd,stemax, 64 deemax,epsil,stmin,ubuf,nwbuf); 65 deemax,epsil,stmin,ubuf,nwbuf); 65 } 66 } 66 67 67 void G4gstmed(G4int itmed, G4String, G4int nma 68 void G4gstmed(G4int itmed, G4String, G4int nmat, G4int isvol, 68 G4int, G4double, G4double, 69 G4int, G4double, G4double, 69 G4double stemax, G4double, G4dou 70 G4double stemax, G4double, G4double, 70 G4double, G4double*, G4int useG3 71 G4double, G4double*, G4int useG3TMLimits) 71 { 72 { 72 // get the pointer to material nmat 73 // get the pointer to material nmat 73 G4Material* material = G3Mat.get(nmat); 74 G4Material* material = G3Mat.get(nmat); 74 75 75 // NB. there is the possibility for redund 76 // NB. there is the possibility for redundancy in the mag field 76 // and user limits objects. Who cares. 77 // and user limits objects. Who cares. 77 // Generate the mag field object 78 // Generate the mag field object 78 // $$$ G4MagneticField* field = new G4Magn 79 // $$$ G4MagneticField* field = new G4MagneticField(ifield, fieldm, tmaxfd); 79 G4MagneticField* field = 0; 80 G4MagneticField* field = 0; 80 81 81 // Generate the user limits object 82 // Generate the user limits object 82 // !!! only "stemax" has its equivalent in 83 // !!! only "stemax" has its equivalent in G4 83 84 84 G4UserLimits* limits = 0; 85 G4UserLimits* limits = 0; 85 if (useG3TMLimits) { 86 if (useG3TMLimits) { 86 limits = new G4UserLimits(); 87 limits = new G4UserLimits(); 87 limits->SetMaxAllowedStep(stemax*cm); 88 limits->SetMaxAllowedStep(stemax*cm); 88 // limits->SetG3DefaultCuts(); 89 // limits->SetG3DefaultCuts(); 89 // this is arranged globally by physi 90 // this is arranged globally by physics manager 90 } 91 } 91 92 92 // Store this medium in the G3Med structur 93 // Store this medium in the G3Med structure 93 G3Med.put(itmed, material, field, limits, 94 G3Med.put(itmed, material, field, limits, isvol); 94 } 95 } 95 96