MaGeGeneratorPositionSampling.hh

Go to the documentation of this file.

//
// ********************************************************************
// * DISCLAIMER                                                       *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.                                                             *
// *                                                                  *
// ********************************************************************
//
// This class is supposed to sample the position of the primary 
// particle randomly in a volume or on a surface. The actual 
// algorithms are not implemented, but the interfaces are ok
//
// History:
// --------
// 28 Oct 2004   L.Pandola    First implementation (interfaces ok)
// 19 Nov 2004   Luciano      Implemented a modified 
//                            version of Dave Jordan algorithm

#ifndef MaGeGeneratorPositionSampling_h
#define MaGeGeneratorPositionSampling_h 1

#include "globals.hh"
#include "G4ThreeVector.hh"
#include "G4VSolid.hh"
#include "G4VPhysicalVolume.hh"
#include "G4RotationMatrix.hh"
#include "G4Point3D.hh"
#include <vector>
#include "G4Navigator.hh"
#include "G4TouchableHistory.hh"

//This class is supposed to sample one random point uniformly from 
//a given volume or from a surface

class MJGeneratorUtil;
class G4Navigator;
class G4VPhysicalVolume;
class MaGeGeneratorPositionSampling 
{
public:
  MaGeGeneratorPositionSampling(); 
  ~MaGeGeneratorPositionSampling();
  
  G4ThreeVector SampleUniformlyInVolume(G4String volName,G4int copynumber);
  G4ThreeVector SampleOnSurface(G4String volName,G4int copynumber);

  void SetGlobalTranslation(G4ThreeVector vect) {globalTranslation = vect;};
  void SetGlobalRotation(G4RotationMatrix mat) {globalRotation = mat;};
  G4ThreeVector GetGlobalTranslation() {return globalTranslation;};
  G4RotationMatrix GetGlobalRotation() {return globalRotation;};

private:
  G4String volumeName;
  MJGeneratorUtil* generatorUtil;
  G4bool InitializeSamplingVolume();
  G4int copyNumber;

  G4Navigator* myNavigator;
  G4VPhysicalVolume* myVolume;

  G4RotationMatrix globalRotation;
  G4ThreeVector globalTranslation;

  G4double radius;

  G4GeometryType fSolid_type;
  G4double fSolid_par_arr[6];

  G4VPhysicalVolume* FindTheDirectMother(G4VPhysicalVolume*);

};
#endif

---