Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 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 // 27 // 28 // by I.Hrivnacova, V.Berejnoi, 29 Oct 99 29 30 #include "G3G4Interface.hh" 31 #include "G3Division.hh" 32 #include "G3VolTableEntry.hh" 33 #include "G3VolTable.hh" 34 #include "globals.hh" 35 #include "G3toG4.hh" 36 37 void PG4gsdvn(G4String *tokens) 38 { 39 // fill the parameter containers 40 G3fillParams(tokens,PTgsdvn); 41 42 // interpret the parameters 43 G4String vname = Spar[0]; 44 G4String vmoth = Spar[1]; 45 G4int ndiv = Ipar[0]; 46 G4int iaxis = Ipar[1]; 47 48 G4gsdvn(vname, vmoth, ndiv, iaxis); 49 } 50 51 void G4CreateCloneVTEWithDivision(G4String vna 52 G3DivType divType, G4int nofDivisions, 53 G4double c0, G4double step) 54 { 55 G3VolTableEntry* vte=0; 56 57 // loop over all mothers 58 for (G4int i=0; i<mvte->GetNoClones(); i++) 59 G3VolTableEntry* mvteClone = mvte->GetClo 60 G4String shape = mvteClone->GetShape(); 61 G4int nmed = mvteClone->GetNmed(); 62 G4String mvteName = mvteClone->GetName(); 63 64 G4String newName = vname; 65 if (i>0) { 66 newName += gSeparator; 67 newName = newName + std::to_string(i); 68 } 69 70 // create new VTE with 0 solid 71 // and let vol table know about it 72 G3VolTableEntry* vteClone 73 = new G3VolTableEntry(newName, shape, 0, 74 G3Vol.PutVTE(vteClone); 75 76 // set mother <-> daughter 77 // (mother/daughter are reset in case an e 78 // needs to be created in G3Division::Upd 79 mvteClone->AddDaughter(vteClone); 80 vteClone->AddMother(mvteClone); 81 82 // create new G3Division 83 G3Division* division 84 = new G3Division(divType, vteClone, mvte 85 nofDivisions, iaxis, nm 86 87 // set division to vte and update it 88 vteClone->SetDivision(division); 89 division->UpdateVTE(); 90 91 if (i == 0) { 92 // keep the first clone copy 93 vte = vteClone; 94 } 95 else { 96 // let vte know about this clone copy 97 vte->AddClone(vteClone); 98 } 99 } 100 } 101 102 void G4gsdvn(G4String vname, G4String vmoth, G 103 { 104 // find mother VTE 105 G3VolTableEntry* mvte = G3Vol.GetVTE(vmoth); 106 107 if (mvte == 0) { 108 G4String text = "G4gsdvn:'" + vmoth + "' h 109 G4Exception("G4gsdvn()", "G3toG40013", Fat 110 return; 111 } 112 else { 113 // a new vte clone copy with division is c 114 // for each mother (clone copy) 115 116 G4CreateCloneVTEWithDivision(vname, mvte, 117 } 118 } 119