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