Gabriele Cosmo
150 lines
4.7 KiB
C++
150 lines
4.7 KiB
C++
//
|
|
// ********************************************************************
|
|
// * License and Disclaimer *
|
|
// * *
|
|
// * The Geant4 software is copyright of the Copyright Holders of *
|
|
// * the Geant4 Collaboration. It is provided under the terms and *
|
|
// * conditions of the Geant4 Software License, included in the file *
|
|
// * LICENSE and available at http://cern.ch/geant4/license . These *
|
|
// * include a list of copyright holders. *
|
|
// * *
|
|
// * 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. Please see the license in the file LICENSE and URL above *
|
|
// * for the full disclaimer and the limitation of liability. *
|
|
// * *
|
|
// * This code implementation is the result of the scientific and *
|
|
// * technical work of the GEANT4 collaboration. *
|
|
// * By using, copying, modifying or distributing the software (or *
|
|
// * any work based on the software) you agree to acknowledge its *
|
|
// * use in resulting scientific publications, and indicate your *
|
|
// * acceptance of all terms of the Geant4 Software license. *
|
|
// ********************************************************************
|
|
//
|
|
// Gorad (Geant4 Open-source Radiation Analysis and Design)
|
|
//
|
|
// Author : Makoto Asai (SLAC National Accelerator Laboratory)
|
|
//
|
|
// Development of Gorad is funded by NASA Johnson Space Center (JSC)
|
|
// under the contract NNJ15HK11B.
|
|
//
|
|
// ********************************************************************
|
|
//
|
|
// GRDetectorConstruction.hh
|
|
// Header file of the detector construction. It reads a GDML file.
|
|
//
|
|
// History
|
|
// September 8th, 2020 : first implementation
|
|
//
|
|
// ********************************************************************
|
|
|
|
#ifndef GRDetectorConstruction_H
|
|
#define GRDetectorConstruction_H 1
|
|
|
|
#include "G4VUserDetectorConstruction.hh"
|
|
#include "globals.hh"
|
|
#include "G4ThreeVector.hh"
|
|
|
|
class G4GDMLParser;
|
|
class GRDetectorConstructionMessenger;
|
|
class GRGeomImpBiasWorld;
|
|
|
|
class GRDetectorConstruction : public G4VUserDetectorConstruction
|
|
{
|
|
friend class GRGeomImpBiasWorld;
|
|
|
|
public:
|
|
GRDetectorConstruction();
|
|
virtual ~GRDetectorConstruction();
|
|
virtual G4VPhysicalVolume* Construct();
|
|
virtual void ConstructSDAndField();
|
|
|
|
private:
|
|
GRDetectorConstructionMessenger* messenger;
|
|
G4GDMLParser* parser;
|
|
G4String gdmlFile;
|
|
G4VPhysicalVolume* fWorld;
|
|
G4bool initialized;
|
|
|
|
static G4double worldSize;
|
|
|
|
public:
|
|
G4bool SetGDMLFile(G4String&);
|
|
const G4String& GetGDMLFile()
|
|
{ return gdmlFile; }
|
|
|
|
static G4double GetWorldSize()
|
|
{ return worldSize; }
|
|
|
|
private:
|
|
void Read();
|
|
|
|
public:
|
|
void ListSolids(G4int);
|
|
void ListLogVols(G4int);
|
|
void ListPhysVols(G4int);
|
|
void ListRegions(G4int);
|
|
G4bool CreateRegion(G4String&,G4String&);
|
|
G4bool CheckOverlap(G4String&,G4int,G4int,G4double);
|
|
|
|
public:
|
|
void ListAllMaterial();
|
|
G4bool ListMaterial(G4String&);
|
|
void DumpNistMaterials();
|
|
G4bool CreateMaterial(G4String&);
|
|
G4bool GetMaterial(G4String&);
|
|
G4int SetMaterial(G4String&,G4String&);
|
|
|
|
private:
|
|
G4bool applyGeomImpBias = false;
|
|
GRGeomImpBiasWorld* geomImpBiasWorld = nullptr;
|
|
struct GeomImpParameters
|
|
{
|
|
G4double radius = -1.;
|
|
G4ThreeVector pos0;
|
|
G4int nLayer = 0;
|
|
G4double radiusT = -1.;
|
|
G4ThreeVector posT;
|
|
G4int factor = 2;
|
|
G4double prob = 1.;
|
|
} geoImpP;
|
|
|
|
public:
|
|
void GeomImp(G4int n,G4double r)
|
|
{
|
|
applyGeomImpBias = true;
|
|
geoImpP.nLayer = n;
|
|
geoImpP.radius = r;
|
|
}
|
|
G4bool ApplyGeomImpBias() const
|
|
{ return applyGeomImpBias; }
|
|
void GeomImpLocate(G4ThreeVector p0)
|
|
{ geoImpP.pos0 = p0; }
|
|
G4double GeomImpInnerRadius(G4double rt)
|
|
{
|
|
if(rt>geoImpP.radius)
|
|
{ return -rt; }
|
|
geoImpP.radiusT = rt;
|
|
return rt;
|
|
}
|
|
G4double GeomImpLocateTgt(G4ThreeVector pT)
|
|
{
|
|
G4ThreeVector dp = geoImpP.pos0 - pT;
|
|
G4double rt = geoImpP.radiusT;
|
|
if(rt<0.)
|
|
{ rt = geoImpP.radius / geoImpP.nLayer; }
|
|
if(dp.mag() >= (geoImpP.radius - rt))
|
|
{ return (geoImpP.radius - rt); }
|
|
geoImpP.posT = pT;
|
|
return 0.;
|
|
}
|
|
void GeomImpFactor(G4int f)
|
|
{ geoImpP.factor = f; }
|
|
void GeomImpProb(G4double p)
|
|
{ geoImpP.prob = p; }
|
|
};
|
|
|
|
#endif
|