CaptCryostatBuilder Class Reference

#include <CaptCryostatBuilder.hh>

Inheritance diagram for CaptCryostatBuilder:

Classes
class	Point

Public Member Functions
	CaptCryostatBuilder (G4String name, EDepSim::Builder *parent)
virtual	~CaptCryostatBuilder ()
virtual G4LogicalVolume *	GetPiece (void)
G4ThreeVector	GetOffset ()
G4ThreeVector	GetTPCOffset ()
double	GetArgonDepth () const
void	SetArgonDepth (double v)
double	GetTPCDepth () const
void	SetTPCDepth (double v)
std::string	GetVesselType () const
void	SetVesselType (std::string v)
Public Member Functions inherited from EDepSim::Builder
	Builder (G4String n, EDepSim::UserDetectorConstruction *c)
	Builder (G4String n, EDepSim::Builder *parent)
virtual	~Builder ()
G4String	GetName (void)
	Return the base name of the object that this builds. More...
G4String	GetLocalName (void)
	Return the base name of the object that this builds. More...
void	SetLocalName (const G4String &name)
	Set the base name of the logical volume that this builds. More...
void	SetOpacity (double v)
	Set the relative opacity of the constructed object. More...
double	GetOpacity (void) const
	Get the relative opacity of the constructed object. More...
void	SetDaughterOpacity (double v)
	Set the relative opacity of the object daughters. More...
EDepSim::UserDetectorConstruction *	GetConstruction (void)
	Return the detector construction that is building this piece. More...
virtual void	SetSensitiveDetector (G4VSensitiveDetector *s)
	Set the sensitive detector for this component. More...
virtual G4VSensitiveDetector *	GetSensitiveDetector (void)
	Get the sensitive detector for this component. More...
virtual void	SetSensitiveDetector (G4String name, G4String type)
	Set the sensitive detector for this component by name. More...
virtual void	SetMaximumHitSagitta (double sagitta)
virtual void	SetMaximumHitLength (double length)
	Set the maximum length of a single hit segment. More...
G4UImessenger *	GetMessenger (void)
	Return the messenger for this constructor. More...
void	SetMessenger (G4UImessenger *m)
	Set the messenger for this constructor. More...
void	AddBuilder (EDepSim::Builder *c)
template<class T >
T &	Get (G4String n)
template<class T >
T *	Find (G4String n)
bool	Check ()
void	SetCheck (bool v)
	Set the check value. More...

Private Types
typedef std::vector< Point >	Shape

Private Member Functions
void	Init (void)
void	DefineCAPTAINVessel ()
	Fill the vessel definition with values for CAPTAIN. More...
void	DefineMiniCAPTAINVessel ()
	Fill the vessel definition with values for miniCAPTAIN. More...

Private Attributes
G4LogicalVolume *	fUllageVolume
G4LogicalVolume *	fLiquidVolume
std::string	fVesselType
	The name of the vessel to be built (currently CAPTAIN or mCAPTAIN). More...
double	fArgonDepth
	The distance from the top of the flange to the liquid argon. More...
double	fTPCDepth
	The distance from the top of the flange to the TPC. More...
Shape	fInnerVessel
	The inner vessel polycone points. More...
Shape	fOuterVessel
	The outer vessel polycone points. More...
Shape	fVesselEnvelope
	The vessel envelop. More...

Additional Inherited Members
Protected Member Functions inherited from EDepSim::Builder
G4Material *	FindMaterial (G4String m)
G4VisAttributes	GetColor (G4LogicalVolume *volume, double opacity=0.0)
G4VisAttributes	GetColor (G4Material *volume, double opacity=0.0)

Detailed Description

Construct the cryostat. This builds the vessel, as well as the argon gas and liquid inside. This builds the vessel directly and fills it with argon gas. The liquid argon is build by a separate constructor which will contain the TPC.

Bug:

The current class is an extremely simplistic model. I think the cryostat is multi-walled, but I don't have the design documents right now. The invariants for this class are the inner diameter and inner heights which define the cold volume. The outer diameter and height should be calculated.

Definition at line 19 of file CaptCryostatBuilder.hh.

Member Typedef Documentation

typedef std::vector<Point> CaptCryostatBuilder::Shape

private

Definition at line 97 of file CaptCryostatBuilder.hh.

Constructor & Destructor Documentation

CaptCryostatBuilder::CaptCryostatBuilder ( G4String  name, EDepSim::Builder *  parent  )

inline

Definition at line 21 of file CaptCryostatBuilder.hh.

        : EDepSim::Builder(name,parent) {Init();};

CaptCryostatBuilder::~CaptCryostatBuilder ( )

virtual

Definition at line 94 of file CaptCryostatBuilder.cc.

{}

Member Function Documentation

void CaptCryostatBuilder::DefineCAPTAINVessel ( )

private

Fill the vessel definition with values for CAPTAIN.

Definition at line 104 of file CaptCryostatBuilder.cc.

                                              {
    SetArgonDepth(24*25.4*CLHEP::mm);  // The design spec for CAPTAIN (24 inch).
    SetTPCDepth(GetArgonDepth()+25*CLHEP::mm);   // I made this one up...

    fInnerVessel.clear();
#include "captainInnerVessel.hxx"
    
    fOuterVessel.clear();
#include "captainOuterVessel.hxx"

    for (Shape::reverse_iterator p = fOuterVessel.rbegin();
         p != fOuterVessel.rend(); ++p) {
        fVesselEnvelope.push_back(*p);
    }
    for (Shape::reverse_iterator p = fInnerVessel.rbegin();
         p != fInnerVessel.rend(); ++p) {
        if (p->fZ >= fVesselEnvelope.back().fZ) continue;
        fVesselEnvelope.push_back(*p);
    }
    std::reverse(fVesselEnvelope.begin(), fVesselEnvelope.end());
    
}

void CaptCryostatBuilder::DefineMiniCAPTAINVessel ( )

private

Fill the vessel definition with values for miniCAPTAIN.

Definition at line 127 of file CaptCryostatBuilder.cc.

                                                  {
    SetArgonDepth(9*25.4*CLHEP::mm);  // The design spec for CAPTAIN (24 inch).
    SetTPCDepth(GetArgonDepth()+25*CLHEP::mm);   // I made this one up...

    fInnerVessel.clear();
#include "miniCaptainInnerVessel.hxx"
    
    fOuterVessel.clear();
#include "miniCaptainOuterVessel.hxx"

    for (Shape::reverse_iterator p = fOuterVessel.rbegin();
         p != fOuterVessel.rend(); ++p) {
        fVesselEnvelope.push_back(*p);
    }
    for (Shape::reverse_iterator p = fInnerVessel.rbegin();
         p != fInnerVessel.rend(); ++p) {
        if (p->fZ >= fVesselEnvelope.back().fZ) continue;
        fVesselEnvelope.push_back(*p);
    }
    std::reverse(fVesselEnvelope.begin(), fVesselEnvelope.end());
    
}

double CaptCryostatBuilder::GetArgonDepth ( ) const

inline

Get (or set) the distance from the flange to the top of the liquid argon.

Definition at line 34 of file CaptCryostatBuilder.hh.

{return fArgonDepth;}

G4ThreeVector CaptCryostatBuilder::GetOffset

(

)

Return the offset of the intended origin of the volume relative to the center of the logical volume. To get the origin at the right location (say originPosition), the logical volume should be positioned at originPosition-GetOffset(). The origin of the cryostat is chosen so that it's at the center of flange (the top of the "bucket" flange" defines the Z origin of the cryostat). The bottom of the wire plane assembly is at the origin (a decision will be made in the future as to whether this is the bottom of the grid, or the bottom of the V plane. This means that the wires for the V plane are at a (very) small positive z coordinate.

Definition at line 96 of file CaptCryostatBuilder.cc.

                                             {
    return G4ThreeVector(0,0,0);
}

G4LogicalVolume * CaptCryostatBuilder::GetPiece ( void  )

virtual

Construct and return a G4 volume for the object. This is a pure virtual function, which means it must be implemented by the inheriting classes. This returns an unplaced logical volume which faces along the Z axis.

Implements EDepSim::Builder.

Definition at line 150 of file CaptCryostatBuilder.cc.

                                                   {

    if (fVesselType == "CAPTAIN") {
        DefineCAPTAINVessel();
    }
    else if (fVesselType == "mCAPTAIN") {
        DefineMiniCAPTAINVessel();
    }
    else {
        std::cout << "Undefine vessel type: " << fVesselType << std::endl;
        std::exit(0);
    }

    std::vector<double> conePlanes;
    std::vector<double> coneMax;
    std::vector<double> coneMin;

    ////////////////////////////////////////////////////////
    // Define the envelope to contain the vessel.
    ////////////////////////////////////////////////////////
    for (Shape::reverse_iterator p = fVesselEnvelope.rbegin();
         p != fVesselEnvelope.rend();
         ++p) {
        conePlanes.push_back(- p->fZ);
        coneMax.push_back(p->fOuter+10*CLHEP::cm);
        coneMin.push_back(0.0);
    }
    
    G4LogicalVolume* logVolume 
        = new G4LogicalVolume(
            new G4Polycone(GetName(),
                           0*CLHEP::degree, 360*CLHEP::degree, conePlanes.size(),
                           conePlanes.data(),
                           coneMin.data(), coneMax.data()),
            FindMaterial("Air"),
            GetName());
    logVolume->SetVisAttributes(G4VisAttributes::Invisible);
    
    ////////////////////////////////////////////////////////
    // Define the outer vessel.
    ////////////////////////////////////////////////////////
    conePlanes.clear();
    coneMax.clear();
    coneMin.clear();
    for (Shape::reverse_iterator p = fOuterVessel.rbegin();
         p != fOuterVessel.rend();
         ++p) {
        conePlanes.push_back(- p->fZ);
        coneMin.push_back(p->fInner);
        coneMax.push_back(p->fOuter);
    }

    G4LogicalVolume* logOuterVessel
        = new G4LogicalVolume(
            new G4Polycone(GetName()+"/OuterVessel",
                           0*CLHEP::degree, 360*CLHEP::degree, conePlanes.size(),
                           conePlanes.data(),
                           coneMin.data(), coneMax.data()),
            FindMaterial("SS_304"),
            GetName()+"/OuterVessel");
    logOuterVessel->SetVisAttributes(GetColor(logOuterVessel));
    
    new G4PVPlacement(NULL,                // rotation.
                      G4ThreeVector(0,0,0),
                      logOuterVessel,           // logical volume
                      logOuterVessel->GetName(), // name
                      logVolume,                // mother  volume
                      false,                    // (not used)
                      0,                        // Copy number (zero)
                      Check());                 // Check overlaps.

    
    ////////////////////////////////////////////////////////
    // Define the inner vessel.
    ////////////////////////////////////////////////////////
    conePlanes.clear();
    coneMax.clear();
    coneMin.clear();
    for (Shape::reverse_iterator p = fInnerVessel.rbegin();
         p != fInnerVessel.rend();
         ++p) {
        conePlanes.push_back(- p->fZ);
        coneMin.push_back(p->fInner);
        coneMax.push_back(p->fOuter);
    }

    G4LogicalVolume* logInnerVessel
        = new G4LogicalVolume(
            new G4Polycone(GetName()+"/InnerVessel",
                           0*CLHEP::degree, 360*CLHEP::degree, conePlanes.size(),
                           conePlanes.data(),
                           coneMin.data(), coneMax.data()),
            FindMaterial("SS_304"),
            GetName()+"/InnerVessel");
    logInnerVessel->SetVisAttributes(GetColor(logInnerVessel));
    
    new G4PVPlacement(NULL,                // rotation.
                      G4ThreeVector(0,0,0),
                      logInnerVessel,           // logical volume
                      logInnerVessel->GetName(), // name
                      logVolume,                // mother  volume
                      false,                    // (not used)
                      0,                        // Copy number (zero)
                      Check());                 // Check overlaps.

    ////////////////////////////////////////////////////////
    // Define the liquid volume.
    ////////////////////////////////////////////////////////
    conePlanes.clear();
    coneMax.clear();
    coneMin.clear();
    for (Shape::reverse_iterator p = fInnerVessel.rbegin();
         p != fInnerVessel.rend();
         ++p) {
        if (p->fZ < GetArgonDepth()) continue;
        conePlanes.push_back(- p->fZ);
        coneMin.push_back(0.0);
        coneMax.push_back(p->fInner);
    }
    if (conePlanes.back() < -GetArgonDepth()) {
        conePlanes.push_back(- GetArgonDepth());
        coneMin.push_back(0.0);
        coneMax.push_back(coneMax.back());
    }
        
    fLiquidVolume
        = new G4LogicalVolume(
            new G4Polycone(GetName()+"/Liquid",
                           0*CLHEP::degree, 360*CLHEP::degree, conePlanes.size(),
                           conePlanes.data(),
                           coneMin.data(), coneMax.data()),
            FindMaterial("Argon_Liquid"),
            GetName()+"/Liquid");
    fLiquidVolume->SetVisAttributes(GetColor(fLiquidVolume));
    
    new G4PVPlacement(NULL,                // rotation.
                      G4ThreeVector(0,0,0),
                      fLiquidVolume,           // logical volume
                      fLiquidVolume->GetName()+"/Liquid", // name
                      logVolume,                // mother  volume
                      false,                    // (not used)
                      0,                        // Copy number (zero)
                      Check());                 // Check overlaps.

    if (fVesselType == "CAPTAIN") {
        CaptImmersedBuilder& immersed = Get<CaptImmersedBuilder>("Immersed");
        G4LogicalVolume* logImmersed = immersed.GetPiece();
        G4ThreeVector p = GetTPCOffset() - immersed.GetOffset();
        new G4PVPlacement(NULL,                   // rotation.
                          p,                      // position
                          logImmersed,            // logical volume
                          logImmersed->GetName(), // name
                          fLiquidVolume,          // mother  volume
                          false,                  // (not used)
                          0,                      // Copy number (zero)
                          Check());               // Check overlaps.
    }
    else if (fVesselType == "mCAPTAIN") {
        MiniCaptImmersedBuilder& immersed
            = Get<MiniCaptImmersedBuilder>("mImmersed");
        G4LogicalVolume* logImmersed = immersed.GetPiece();
        G4ThreeVector p = GetTPCOffset() - immersed.GetOffset();
        new G4PVPlacement(NULL,                   // rotation.
                          p,                      // position
                          logImmersed,            // logical volume
                          logImmersed->GetName(), // name
                          fLiquidVolume,          // mother  volume
                          false,                  // (not used)
                          0,                      // Copy number (zero)
                          Check());               // Check overlaps.
    }
    else {
        std::cout << "Undefined immersed volume" << std::endl;
        std::exit(0);
    }

    
    ////////////////////////////////////////////////////////
    // Define the ullage volume.
    ////////////////////////////////////////////////////////
    conePlanes.clear();
    coneMax.clear();
    coneMin.clear();
    for (Shape::reverse_iterator p = fInnerVessel.rbegin();
         p != fInnerVessel.rend();
         ++p) {
        if (p->fZ >= GetArgonDepth()-1*CLHEP::mm) continue;
        if (conePlanes.empty() ) {
            conePlanes.push_back(-GetArgonDepth()+1*CLHEP::um);
            coneMin.push_back(0.0);
            coneMax.push_back(p->fInner);
        }
        conePlanes.push_back(- p->fZ);
        coneMin.push_back(0.0);
        coneMax.push_back(p->fInner);
    }
        
    fUllageVolume
        = new G4LogicalVolume(
            new G4Polycone(GetName()+"/Ullage",
                           0*CLHEP::degree, 360*CLHEP::degree, conePlanes.size(),
                           conePlanes.data(),
                           coneMin.data(), coneMax.data()),
            FindMaterial("Argon_Gas"),
            GetName()+"/Ullage");
    fUllageVolume->SetVisAttributes(GetColor(fUllageVolume));
    
    new G4PVPlacement(NULL,                // rotation.
                      G4ThreeVector(0,0,0),
                      fUllageVolume,           // logical volume
                      fUllageVolume->GetName(), // name
                      logVolume,                // mother  volume
                      false,                    // (not used)
                      0,                        // Copy number (zero)
                      Check());                 // Check overlaps.

    if (fVesselType == "CAPTAIN") {
        CaptExposedBuilder& exposed = Get<CaptExposedBuilder>("Exposed");
        G4LogicalVolume* logExposed = exposed.GetPiece();
        G4ThreeVector p(0.0,0.0,-GetArgonDepth());
        p -= exposed.GetOffset();
        new G4PVPlacement(NULL,                   // rotation.
                          p,                      // position
                          logExposed,            // logical volume
                          logExposed->GetName(), // name
                          fUllageVolume,          // mother  volume
                          false,                  // (not used)
                          0,                      // Copy number (zero)
                          Check());               // Check overlaps.
    }
    else if (fVesselType == "mCAPTAIN") {
        MiniCaptExposedBuilder& exposed
            = Get<MiniCaptExposedBuilder>("mExposed");
        G4LogicalVolume* logExposed = exposed.GetPiece();
        if (logExposed) {
            G4ThreeVector p(0.0,0.0,-GetArgonDepth());
            p -= exposed.GetOffset();
            new G4PVPlacement(NULL,                   // rotation.
                              p,                      // position
                              logExposed,            // logical volume
                              logExposed->GetName(), // name
                              fUllageVolume,          // mother  volume
                              false,                  // (not used)
                              0,                      // Copy number (zero)
                              Check());               // Check overlaps.
        }
    }
    else {
        std::cout << "Undefined exposed volume" << std::endl;
        std::exit(0);
    }

    return logVolume;
}

double CaptCryostatBuilder::GetTPCDepth ( ) const

inline

Get (or set) the distance from the flange to the TPC center.

Definition at line 40 of file CaptCryostatBuilder.hh.

{return fTPCDepth;}

G4ThreeVector CaptCryostatBuilder::GetTPCOffset

(

)

The position of the TPC origin (center of the grid plane defining the drift region) relative to the center of the Cryostat flange. The TPC origin should be arranged so that it will be at center of the global coordinate system.

Definition at line 100 of file CaptCryostatBuilder.cc.

                                                {
    return G4ThreeVector(0,0,-GetTPCDepth());
}

std::string CaptCryostatBuilder::GetVesselType ( ) const

inline

Get (or set) the type of vessel that should be built.

Definition at line 46 of file CaptCryostatBuilder.hh.

{return fVesselType;}

void CaptCryostatBuilder::Init ( void  )

private

Definition at line 81 of file CaptCryostatBuilder.cc.

                                   {
    SetMessenger(new CaptCryostatMessenger(this));

    SetVesselType("CAPTAIN");
    
    SetSensitiveDetector("cryo","segment");

    AddBuilder(new CaptImmersedBuilder("Immersed",this));
    AddBuilder(new CaptExposedBuilder("Exposed",this));
    AddBuilder(new MiniCaptImmersedBuilder("mImmersed",this));
    AddBuilder(new MiniCaptExposedBuilder("mExposed",this));
}

void CaptCryostatBuilder::SetArgonDepth ( double  v )

inline

Definition at line 35 of file CaptCryostatBuilder.hh.

{fArgonDepth = v;}

void CaptCryostatBuilder::SetTPCDepth ( double  v )

inline

Definition at line 41 of file CaptCryostatBuilder.hh.

{fTPCDepth = v;}

void CaptCryostatBuilder::SetVesselType ( std::string  v )

inline

Definition at line 47 of file CaptCryostatBuilder.hh.

{fVesselType=v;}

Member Data Documentation

double CaptCryostatBuilder::fArgonDepth

private

The distance from the top of the flange to the liquid argon.

Definition at line 85 of file CaptCryostatBuilder.hh.

Shape CaptCryostatBuilder::fInnerVessel

private

The inner vessel polycone points.

Definition at line 100 of file CaptCryostatBuilder.hh.

G4LogicalVolume* CaptCryostatBuilder::fLiquidVolume

private

The volume for the liquid argon. This is set when the cryostat vessel is constructed and fills the liquid part of the inner void. It can be passed to the daughter volumes to define the shape of the liquid.

Definition at line 79 of file CaptCryostatBuilder.hh.

Shape CaptCryostatBuilder::fOuterVessel

private

The outer vessel polycone points.

Definition at line 102 of file CaptCryostatBuilder.hh.

double CaptCryostatBuilder::fTPCDepth

private

The distance from the top of the flange to the TPC.

Definition at line 88 of file CaptCryostatBuilder.hh.

G4LogicalVolume* CaptCryostatBuilder::fUllageVolume

private

The volume for the ullage. This is set when the cryostat vessel is constructed and fills the gas part of the inner void. It can be passed to daughter volumes to define the shape of the ullage.

Definition at line 74 of file CaptCryostatBuilder.hh.

Shape CaptCryostatBuilder::fVesselEnvelope

private

The vessel envelop.

Definition at line 104 of file CaptCryostatBuilder.hh.

std::string CaptCryostatBuilder::fVesselType

private

The name of the vessel to be built (currently CAPTAIN or mCAPTAIN).

Definition at line 82 of file CaptCryostatBuilder.hh.

---

The documentation for this class was generated from the following files:

• edep-sim/src/captain/CaptCryostatBuilder.hh
• edep-sim/src/captain/CaptCryostatBuilder.cc