| Geant4 Cross Reference |
>> 1 // This code implementation is the intellectual property of
>> 2 // the GEANT4 collaboration.
1 // 3 //
2 // ******************************************* << 4 // By copying, distributing or modifying the Program (or any work
3 // * License and Disclaimer << 5 // based on the Program) you indicate your acceptance of this statement,
4 // * << 6 // and all its terms.
5 // * The Geant4 software is copyright of th <<
6 // * the Geant4 Collaboration. It is provided <<
7 // * conditions of the Geant4 Software License <<
8 // * LICENSE and available at http://cern.ch/ <<
9 // * include a list of copyright holders. <<
10 // * <<
11 // * Neither the authors of this software syst <<
12 // * institutes,nor the agencies providing fin <<
13 // * work make any representation or warran <<
14 // * regarding this software system or assum <<
15 // * use. Please see the license in the file <<
16 // * for the full disclaimer and the limitatio <<
17 // * <<
18 // * This code implementation is the result <<
19 // * technical work of the GEANT4 collaboratio <<
20 // * By using, copying, modifying or distri <<
21 // * any work based on the software) you ag <<
22 // * use in resulting scientific publicati <<
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* <<
25 // <<
26 // 7 //
>> 8 // $Id: G4gsvolu.cc,v 1.7 1999/12/15 14:49:44 gunter Exp $
>> 9 // GEANT4 tag $Name: geant4-03-00 $
27 // 10 //
28 // by I.Hrivnacova, 13.10.99 11 // by I.Hrivnacova, 13.10.99
29 12
30 #include <iomanip> << 13 #include "g4std/iomanip"
31 #include "G3VolTable.hh" 14 #include "G3VolTable.hh"
32 #include "globals.hh" 15 #include "globals.hh"
33 #include "G3toG4.hh" 16 #include "G3toG4.hh"
34 #include "G3toG4MakeSolid.hh" 17 #include "G3toG4MakeSolid.hh"
35 18
36 void PG4gsvolu(G4String *tokens) { << 19 void PG4gsvolu(G4String tokens[]) {
37 // fill the parameter containers 20 // fill the parameter containers
38 G3fillParams(tokens,PTgsvolu); 21 G3fillParams(tokens,PTgsvolu);
39 22
40 // interpret the parameters 23 // interpret the parameters
41 G4String vname = Spar[0]; 24 G4String vname = Spar[0];
42 G4String shape = Spar[1]; 25 G4String shape = Spar[1];
43 G4int nmed = Ipar[0]; 26 G4int nmed = Ipar[0];
44 G4int npar = Ipar[1]; 27 G4int npar = Ipar[1];
45 G4double *pars = Rpar; 28 G4double *pars = Rpar;
46 29
47 G4gsvolu(vname, shape, nmed, pars, npar); 30 G4gsvolu(vname, shape, nmed, pars, npar);
48 } 31 }
49 32
50 G3VolTableEntry* G4CreateVTE(G4String vname, G 33 G3VolTableEntry* G4CreateVTE(G4String vname, G4String shape, G4int nmed,
51 G4double rpar[], << 34 G4double Rpar[], G4int npar)
52 { 35 {
53 // create the solid 36 // create the solid
54 G4bool hasNegPars; 37 G4bool hasNegPars;
55 G4bool deferred; 38 G4bool deferred;
56 G4bool okAxis[3]; 39 G4bool okAxis[3];
57 G4VSolid* solid 40 G4VSolid* solid
58 = G3toG4MakeSolid(vname, shape, rpar, npar << 41 = G3toG4MakeSolid(vname, shape, Rpar, npar, hasNegPars, deferred, okAxis);
59 42
60 // if solid has been deferred 43 // if solid has been deferred
61 // VTE is created with hasNegPars = true 44 // VTE is created with hasNegPars = true
62 if (deferred) hasNegPars = true; 45 if (deferred) hasNegPars = true;
63 46
64 // create VTE 47 // create VTE
65 G3VolTableEntry* vte 48 G3VolTableEntry* vte
66 = new G3VolTableEntry(vname, shape, rpar, << 49 = new G3VolTableEntry(vname, shape, Rpar, npar, nmed, solid, hasNegPars);
67 G3Vol.PutVTE(vte); 50 G3Vol.PutVTE(vte);
68 51
69 return vte; 52 return vte;
70 } 53 }
71 54
72 void G4gsvolu(G4String vname, G4String shape, << 55 void G4gsvolu(G4String vname, G4String shape, G4int nmed, G4double* Rpar,
73 G4int npar) 56 G4int npar)
74 { 57 {
75 /* 58 /*
76 G4cout << "Creating logical volume " << vnam 59 G4cout << "Creating logical volume " << vname << " shape " << shape
77 << " nmed " << nmed << " #pars "<< npar < 60 << " nmed " << nmed << " #pars "<< npar << " parameters (cm): ";
78 for (int ipar=0; ipar< npar; ipar++) G4cout << 61 for (int ipar=0; ipar< npar; ipar++) G4cout << G4std::setw(8) << Rpar[ipar];
79 G4cout << G4endl; 62 G4cout << G4endl;
80 */ 63 */
81 if (G3Vol.GetVTE(vname)) { 64 if (G3Vol.GetVTE(vname)) {
82 // abort if VTE with given name exists 65 // abort if VTE with given name exists
83 G4String text = "G4gsvolu: Attempt to crea << 66 G4Exception("G4gsvolu: Attempt to create volume " + vname + " twice.");
84 G4Exception("G4gsvolu()", "G3toG40024", Fa <<
85 return; <<
86 } 67 }
87 else { 68 else {
88 G4CreateVTE(vname, shape, nmed, rpar, npar << 69 G4CreateVTE(vname, shape, nmed, Rpar, npar);
89 } 70 }
90 } 71 }
91 72