| 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 //
29 // John Allison 12th September 2004
30
31 // Class Description:
32 //
33 // Calculates the mass of a geometry tree taking into account daughters
34 // up to the depth specified in the G4PhysicalVolumeModel. Culling is
35 // ignored so that all volumes are seen.
36 //
37 // Do not use this for a "parallel world" for which materials are not
38 // defined. Use only for the material world.
39 //
40 // The calculation is quite tricky, since it involves subtracting the
41 // mass of that part of the mother that is occupied by each daughter and
42 // then adding the mass of the daughter, and so on down the heirarchy.
43 //
44 // Usage for a given G4PhysicalVolumeModel* pvModel:
45 // G4PhysicalVolumeMassScene massScene(pvModel);
46 // pvModel->DescribeYourselfTo (massScene);
47 // G4double volume = massScene.GetVolume();
48 // G4double mass = massScene.GetMass();
49 // massScene.Reset();
50 // See, for example, G4ASCIITreeSceneHandler::EndModeling().
51
52 #ifndef G4PHYSICALVOLUMEMASSSCENE_HH
53 #define G4PHYSICALVOLUMEMASSSCENE_HH
54
55 #include "G4PseudoScene.hh"
56
57 #include <deque>
58
59 class G4VPhysicalVolume;
60 class G4LogicalVolume;
61 class G4PhysicalVolumeModel;
62 class G4Material;
63
64 class G4PhysicalVolumeMassScene: public G4PseudoScene
65 {
66
67 public:
68 G4PhysicalVolumeMassScene (G4PhysicalVolumeModel*);
69 virtual ~G4PhysicalVolumeMassScene ();
70
71 public: // With description
72
73 G4double GetVolume () const {return fVolume;}
74 // Overall volume.
75
76 G4double GetMass () const {return fMass;}
77 // Mass of whole tree, i.e., accounting for all daughters.
78
79 void Reset ();
80 // Reset for subsequent re-use.
81
82 private:
83 void ProcessVolume (const G4VSolid&);
84 G4PhysicalVolumeModel* fpPVModel;
85 G4double fVolume;
86 G4double fMass;
87 G4VPhysicalVolume* fpLastPV;
88 G4int fPVPCount;
89 G4int fLastDepth;
90 G4double fLastDensity;
91 std::deque<G4double> fDensityStack;
92 };
93
94 #endif
95