Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/human_phantom/src/G4MIRDRibCage.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/advanced/human_phantom/src/G4MIRDRibCage.cc (Version 11.3.0) and /examples/advanced/human_phantom/src/G4MIRDRibCage.cc (Version 8.0)


  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 // Previous authors: G. Guerrieri, S. Guatelli    
 28 // Authors (since 2007): S. Guatelli, Universi    
 29 //                                                
 30 //                                                
 31 #include "G4MIRDRibCage.hh"                       
 32                                                   
 33 #include "globals.hh"                             
 34 #include "G4SystemOfUnits.hh"                     
 35 #include "G4SDManager.hh"                         
 36 #include "G4VisAttributes.hh"                     
 37 #include "G4HumanPhantomMaterial.hh"              
 38 #include "G4SubtractionSolid.hh"                  
 39 #include "G4EllipticalTube.hh"                    
 40 #include "G4PVReplica.hh"                         
 41 #include "G4Box.hh"                               
 42 #include "G4PVPlacement.hh"                       
 43 #include "G4HumanPhantomColour.hh"                
 44 G4MIRDRibCage::G4MIRDRibCage():fPhysRib1(nullp    
 45               fPhysRib4(nullptr), fPhysRib5(nu    
 46               fPhysRib7(nullptr), fPhysRib8(nu    
 47               fPhysRib10(nullptr), fPhysRib11(    
 48 {;}                                               
 49                                                   
 50 G4VPhysicalVolume* G4MIRDRibCage::Construct(co    
 51               const G4String& colourName, G4bo    
 52                                                   
 53 {                                                 
 54   auto* material = new G4HumanPhantomMaterial(    
 55                                                   
 56   G4cout<<"Construct "<<volumeName<<" with mot    
 57                                                   
 58   auto* skeleton = material -> GetMaterial("sk    
 59   auto* soft = material -> GetMaterial("soft_t    
 60                                                   
 61   delete material;                                
 62                                                   
 63   G4double dx= 17. *cm; // a2                     
 64   G4double dy= 9.8 * cm; //b2                     
 65   G4double thickness= 32.4 * cm; // z2/2 of ca    
 66                                                   
 67   auto* outCage = new G4EllipticalTube("outCag    
 68                                                   
 69   dx = 16.4 * cm; // a1                           
 70   dy = 9.2 * cm; // b1                            
 71   G4double dz = 34. *cm; // z2/2                  
 72                                                   
 73   auto* inCage = new G4EllipticalTube("inCage"    
 74                                                   
 75   auto* cage = new G4SubtractionSolid("Cage",     
 76                 outCage,                          
 77                 inCage, nullptr, G4ThreeVector    
 78                                                   
 79                                                   
 80   auto* logicRibCage = new G4LogicalVolume(cag    
 81                                                   
 82   G4VPhysicalVolume* physRibCage = new G4PVPla    
 83                  // with respect to the trunk     
 84                  "physicalRibCage",               
 85                  logicRibCage,                    
 86                  mother,                          
 87                  false,                           
 88                  0, true);                        
 89                                                   
 90   G4double xx = 17.*cm;                           
 91   G4double yy = 9.8*cm;                           
 92   G4double ribThickness = 1.4*cm;                 
 93   auto* rib_out = new G4EllipticalTube("rib_ou    
 94                                                   
 95   xx = 16.5 *cm;                                  
 96   yy = 9.3 * cm;                                  
 97   G4double zz = 1.5 * cm;                         
 98   auto* rib_in = new G4EllipticalTube("rib_in"    
 99   auto* rib = new G4SubtractionSolid("rib",rib    
100                                                   
101   auto* logicRib= new G4LogicalVolume(rib, ske    
102                                                   
103   fPhysRib1 = new G4PVPlacement(nullptr,G4Thre    
104              // with respect to the trunk         
105              "physicalRib",                       
106              logicRib,                            
107              physRibCage,                         
108              false,                               
109              0, true);                            
110                                                   
111   fPhysRib2 = new G4PVPlacement(nullptr,G4Thre    
112              // with respect to the trunk         
113              "physicalRib",                       
114              logicRib,                            
115              physRibCage,                         
116              false,                               
117              0, true);                            
118                                                   
119   fPhysRib3 = new G4PVPlacement(nullptr,G4Thre    
120              // with respect to the trunk         
121              "physicalRib",                       
122              logicRib,                            
123              physRibCage,                         
124              false,                               
125              0, true);                            
126                                                   
127   fPhysRib4 = new G4PVPlacement(nullptr,G4Thre    
128              // with respect to the trunk         
129              "physicalRib",                       
130              logicRib,                            
131              physRibCage,                         
132              false,                               
133              0, true);                            
134                                                   
135   fPhysRib5 = new G4PVPlacement(nullptr,G4Thre    
136              // with respect to the trunk         
137              "physicalRib",                       
138              logicRib,                            
139              physRibCage,                         
140              false,                               
141              0, true);                            
142                                                   
143   fPhysRib6 = new G4PVPlacement(nullptr,G4Thre    
144              // with respect to the trunk         
145              "physicalRib",                       
146              logicRib,                            
147              physRibCage,                         
148              false,                               
149              0, true);                            
150                                                   
151   fPhysRib7 = new G4PVPlacement(nullptr,G4Thre    
152              // with respect to the trunk         
153              "physicalRib",                       
154              logicRib,                            
155              physRibCage,                         
156              false,                               
157              0, true);                            
158                                                   
159   fPhysRib8 = new G4PVPlacement(nullptr,G4Thre    
160              // with respect to the trunk         
161              "physicalRib",                       
162              logicRib,                            
163              physRibCage,                         
164              false,                               
165              0, true);                            
166                                                   
167   fPhysRib9 = new G4PVPlacement(nullptr,G4Thre    
168              // with respect to the trunk         
169              "physicalRib",                       
170              logicRib,                            
171              physRibCage,                         
172              false,                               
173              0, true);                            
174                                                   
175   fPhysRib10 = new G4PVPlacement(nullptr, G4Th    
176         // with respect to the trunk              
177         "physicalRib",                            
178         logicRib,                                 
179         physRibCage,                              
180         false,                                    
181         0, true);                                 
182                                                   
183   fPhysRib11 = new G4PVPlacement(nullptr,G4Thr    
184         // with respect to the trunk              
185         "physicalRib",                            
186         logicRib,                                 
187         physRibCage,                              
188         false,                                    
189         0, true);                                 
190                                                   
191   fPhysRib12 = new G4PVPlacement(nullptr,G4Thr    
192         // with respect to the trunk              
193         "physicalRib",                            
194         logicRib,                                 
195         physRibCage,                              
196         false,                                    
197         0, true);                                 
198                                                   
199   // Visualization Attributes                     
200   logicRibCage -> SetVisAttributes(G4VisAttrib    
201                                                   
202   //G4VisAttributes* RibCageVisAtt = new G4Vis    
203                                                   
204   auto* colourPointer = new G4HumanPhantomColo    
205   G4Colour colour = colourPointer -> GetColour    
206   auto* RibCageVisAtt = new G4VisAttributes(co    
207   RibCageVisAtt->SetForceSolid(wireFrame);        
208   logicRib->SetVisAttributes(RibCageVisAtt);      
209                                                   
210   G4cout << "RibCage created !!!!!!" << G4endl    
211   // Testing Pelvis Volume                        
212   G4double RibCageVol = logicRib->GetSolid()->    
213   G4cout << "Volume of RibCage = " << ((RibCag    
214                                                   
215   // Testing RibCage Material                     
216   G4String RibCageMat = logicRib->GetMaterial(    
217   G4cout << "Material of RibCage = " << RibCag    
218                                                   
219   // Testing Density                              
220   G4double RibCageDensity = logicRib->GetMater    
221   G4cout << "Density of Material = " << RibCag    
222                                                   
223   // Testing Mass                                 
224   G4double RibCageMass = (RibCageVol)* RibCage    
225   G4cout << "Mass of RibCage = " << (RibCageMa    
226                                                   
227   return physRibCage;                             
228 }                                                 
229