source: trunk/source/event/include/G4AdjointPosOnPhysVolGenerator.hh@ 1272

Last change on this file since 1272 was 1228, checked in by garnier, 16 years ago

update geant4.9.3 tag
File size: 5.9 KB
RevLine 
[1197]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// $Id: G4AdjointPosOnPhysVolGenerator.hh,v 1.2 2009/11/18 17:57:59 gcosmo Exp $
[1228]27// GEANT4 tag $Name: geant4-09-03 $
[1197]28//
29/////////////////////////////////////////////////////////////////////////////////
30// Class Name: G4AdjointPosOnPhysVolGenerator
31// Author: L. Desorgher
32// Organisation: SpaceIT GmbH
33// Contract: ESA contract 21435/08/NL/AT
34// Customer: ESA/ESTEC
35/////////////////////////////////////////////////////////////////////////////////
36//
37// CHANGE HISTORY
38// --------------
39// ChangeHistory:
40// 1st June 2006 creation by L. Desorgher
41//
42//-------------------------------------------------------------
43// Documentation:
44// This class is responsible for the generation of primary adjoint particle on the external surface of a user selected volume.
45// The particle are generated uniformly on the surface with the angular distribution set to a cosine law relative to normal of the surface.
46// It is equivalent to the flux going in from the surface if an isotropic flux is considered outside.
47// It uses ray tracking technique and can be applied to all kind of convex volume. Uisng the ray tracking technique the area
48// of the external surface is also computed. The area is needed to fix the weight of the primary adjoint particle.
49// At the time of the development of this class, generation of particle on volume surface and computation of surface was limited in G4,
50// therfore the general ray tracking technique was adopted. It could be now (2009) that direct method of G4VSolid could be used instead. To be checked!
51//
52//
53//
54#ifndef G4AdjointPosOnPhysVolGenerator_h
55#define G4AdjointPosOnPhysVolGenerator_h 1
56
57#include "G4VPhysicalVolume.hh"
58#include "G4AffineTransform.hh"
59#include "G4ThreeVector.hh"
60
61class G4VSolid;
62
63class G4AdjointPosOnPhysVolGenerator
64///////////////////////
65{
66
67//--------
68 public: //without description
69//--------
70
71 static G4AdjointPosOnPhysVolGenerator* GetInstance();
72
73//--------
74 public: //public methods
75//--------
76 G4VPhysicalVolume* DefinePhysicalVolume(const G4String& aName);
77 void DefinePhysicalVolume1(const G4String& aName);
78 G4double ComputeAreaOfExtSurface();
79 G4double ComputeAreaOfExtSurface(G4int NStat);
80 G4double ComputeAreaOfExtSurface(G4double epsilon);
81 G4double ComputeAreaOfExtSurface(G4VSolid* aSolid);
82 G4double ComputeAreaOfExtSurface(G4VSolid* aSolid,G4int NStat);
83 G4double ComputeAreaOfExtSurface(G4VSolid* aSolid,G4double epsilon);
84
85 void GenerateAPositionOnTheExtSurfaceOfASolid(G4VSolid* aSolid,G4ThreeVector& p, G4ThreeVector& direction);
86 void GenerateAPositionOnTheExtSurfaceOfTheSolid(G4ThreeVector& p, G4ThreeVector& direction);
87 void GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume(G4ThreeVector& p, G4ThreeVector& direction);
88 void GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume(G4ThreeVector& p, G4ThreeVector& direction,
89 G4double& costh_to_normal);
90
91 //inline public methods
92
93 inline void SetSolid(G4VSolid* aSolid){theSolid=aSolid;}
94 inline G4double GetAreaOfExtSurfaceOfThePhysicalVolume(){return AreaOfExtSurfaceOfThePhysicalVolume;}
95 inline G4double GetCosThDirComparedToNormal(){return CosThDirComparedToNormal;}
96
97//---------
98 private: //private methods
99//---------
100 G4AdjointPosOnPhysVolGenerator();
101 ~G4AdjointPosOnPhysVolGenerator();
102 G4double ComputeAreaOfExtSurfaceStartingFromSphere(G4VSolid* aSolid,G4int NStat);
103 G4double ComputeAreaOfExtSurfaceStartingFromBox(G4VSolid* aSolid,G4int NStat);
104 void GenerateAPositionOnASolidBoundary(G4VSolid* aSolid,G4ThreeVector& p, G4ThreeVector& direction);
105 G4double GenerateAPositionOnASphereBoundary(G4VSolid* aSolid,G4ThreeVector& p, G4ThreeVector& direction);
106 G4double GenerateAPositionOnABoxBoundary(G4VSolid* aSolid,G4ThreeVector& p, G4ThreeVector& direction);
107 void ComputeTransformationFromPhysVolToWorld();
108
109//---------
110 private: //attributes
111//---------
112 static G4AdjointPosOnPhysVolGenerator* theInstance;
113 G4VSolid* theSolid;
114 G4VPhysicalVolume* thePhysicalVolume;
115 G4int NStat;
116 G4double epsilon;
117 G4bool UseSphere;
118 G4String ModelOfSurfaceSource;
119 G4double ExtSourceRadius;
120 G4double ExtSourceDx,ExtSourceDy,ExtSourceDz ;
121 G4AffineTransform theTransformationFromPhysVolToWorld;
122 G4double AreaOfExtSurfaceOfThePhysicalVolume;
123 G4double CosThDirComparedToNormal;
124};
125
126#endif
127
Note: See TracBrowser for help on using the repository browser.