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