CaptDriftRegionBuilder Class Reference

#include <CaptDriftRegionBuilder.hh>

Inheritance diagram for CaptDriftRegionBuilder:

Public Member Functions
	CaptDriftRegionBuilder (G4String name, EDepSim::Builder *parent)
virtual	~CaptDriftRegionBuilder ()
virtual G4LogicalVolume *	GetPiece (void)
double	GetHeight ()
	Get the total height of the drift region. More...
void	SetWirePlaneSpacing (double v)
	Set the spacing between the wire planes. More...
double	GetWirePlaneSpacing ()
G4ThreeVector	GetOffset ()
double	GetXPlaneOffset ()
	Get the offset of the X plane from the top of the drift region. More...
double	GetVPlaneOffset ()
	Get the offset of the X plane from the top of the drift region. More...
double	GetUPlaneOffset ()
	Get the offset of the X plane from the top of the drift region. More...
double	GetGridPlaneOffset ()
	Get the offset of the grid from the top of the drift region. More...
double	GetGroundPlaneOffset ()
void	SetApothem (double v)
double	GetApothem ()
void	SetRadius (double v)
double	GetRadius () const
void	SetDriftLength (double v)
double	GetDriftLength ()
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 Member Functions
void	Init (void)

Private Attributes
double	fApothem
	The radius of the circle that can fit inside. More...
double	fDriftLength
	The distance from the cathode to the first wire plane. More...
double	fWirePlaneSpacing
	The distance between the wire planes. 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 drift region of the TPC. The wire planes are located at the top of the drift region, with the first wire for each plane located at the most negative X coordinate. The U plane runs from negative Y to positive Y. The V plane runs from positive Y to negative Y. The planes are oriented from top to bottom X, U, V (in otherwords, the first plane seen by the electrons is the V plane). The X plane is the collection plane.

Definition at line 14 of file CaptDriftRegionBuilder.hh.

Constructor & Destructor Documentation

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

inline

Definition at line 16 of file CaptDriftRegionBuilder.hh.

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

CaptDriftRegionBuilder::~CaptDriftRegionBuilder ( )

virtual

Definition at line 97 of file CaptDriftRegionBuilder.cc.

{}

Member Function Documentation

double CaptDriftRegionBuilder::GetApothem ( )

inline

Definition at line 33 of file CaptDriftRegionBuilder.hh.

{return fApothem;}

double CaptDriftRegionBuilder::GetDriftLength ( )

inline

Definition at line 48 of file CaptDriftRegionBuilder.hh.

{return fDriftLength;}

double CaptDriftRegionBuilder::GetGridPlaneOffset

(

)

Get the offset of the grid from the top of the drift region.

Definition at line 130 of file CaptDriftRegionBuilder.cc.

                                                  {
    CaptWirePlaneBuilder& wires = Get<CaptWirePlaneBuilder>("XPlane");
    return 4.0*GetWirePlaneSpacing() + wires.GetHeight()/2;
}

double CaptDriftRegionBuilder::GetGroundPlaneOffset

(

)

Get the offset of the ground from the top of the drift region. The distance between the grid and the ground defines the drift distance for the electrons.

Definition at line 110 of file CaptDriftRegionBuilder.cc.

                                                    {
    CaptWirePlaneBuilder& wires = Get<CaptWirePlaneBuilder>("XPlane");
    return 0.0*GetWirePlaneSpacing() + wires.GetHeight()/2;
}

double CaptDriftRegionBuilder::GetHeight

(

)

Get the total height of the drift region.

Definition at line 99 of file CaptDriftRegionBuilder.cc.

                                         {
    CaptWirePlaneBuilder& wires = Get<CaptWirePlaneBuilder>("XPlane");
    return GetDriftLength() + GetGridPlaneOffset() + wires.GetHeight();
}

G4ThreeVector CaptDriftRegionBuilder::GetOffset

(

)

Return the offset of the intended origin of the volume relative to the center of the logical volume. To get the intended origin at the right location (say originPosition), the logical volume should be positioned at originPosition-GetOffset(). For the drift region, the intended origin is at the "center" of the bottom of the wire planes is at the global origin.

Definition at line 104 of file CaptDriftRegionBuilder.cc.

                                                {
    double zOffset = GetHeight()/2.0;
    zOffset -= GetGridPlaneOffset();
    return G4ThreeVector(0.0, 0.0, zOffset);
}

G4LogicalVolume * CaptDriftRegionBuilder::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 135 of file CaptDriftRegionBuilder.cc.

                                                      {

    double rInner[] = {0.0, 0.0};
    double rOuter[] = {GetApothem(), GetApothem()};
    double zPlane[] = {GetHeight()/2, -GetHeight()/2};

    EDepSimLog("Construct " << GetName()
            << " with " << GetHeight()/mm << " mm height"
            << " and " <<  GetDriftLength()/mm << " mm drift");

    G4LogicalVolume *logVolume
        = new G4LogicalVolume(new G4Polyhedra(GetName(),
                                              90*degree, 450*degree,
                                              6, 2,
                                              zPlane, rInner, rOuter),
                              FindMaterial("Argon_Liquid"),
                              GetName());
    logVolume->SetVisAttributes(GetColor(logVolume));
    G4FieldManager* manager = new G4FieldManager();
    G4ElectroMagneticField* field = new EDepSim::UniformField(
        G4ThreeVector(0,0,0), G4ThreeVector(0.0,0.0,500.0*volt/cm));
    EDepSim::EMFieldSetup* fieldSetup
        = new EDepSim::EMFieldSetup(field,manager);
    if (!fieldSetup) {
        EDepSimError("Field not properly created");
        throw std::runtime_error("Field not properly created");
    }
    logVolume->SetFieldManager(manager,true);
    if (GetSensitiveDetector()) {
        logVolume->SetSensitiveDetector(GetSensitiveDetector());
        logVolume->SetUserLimits(new G4UserLimits(1.0*mm));
    }

    // The wire planes are rotated so that the local Z axis points into the
    // drift volume, and the local X axis points along the original X axis
    // (this puts the local Y axis along the -Y global axis.
    G4RotationMatrix* xRotation = new G4RotationMatrix();
    xRotation->rotateX(180*degree);

    CaptWirePlaneBuilder& xPlane = Get<CaptWirePlaneBuilder>("XPlane");
    xPlane.SetApothem(GetApothem());
    G4LogicalVolume *logX = xPlane.GetPiece();
    new G4PVPlacement(xRotation,                // rotation.
                      G4ThreeVector(0,0,
                                    (GetHeight()/2 - GetXPlaneOffset())),
                      logX,                     // logical volume
                      logX->GetName(),          // name
                      logVolume,                // mother  volume
                      false,                    // (not used)
                      0,                        // Copy number (zero)
                      Check());                 // Check overlaps.


    G4RotationMatrix* vRotation = new G4RotationMatrix();
    vRotation->rotateX(180*degree);
    vRotation->rotateZ(-60*degree);

    CaptWirePlaneBuilder& vPlane = Get<CaptWirePlaneBuilder>("VPlane");
    vPlane.SetApothem(GetApothem());
    G4LogicalVolume *logV = vPlane.GetPiece();
    new G4PVPlacement(vRotation,                // rotation.
                      G4ThreeVector(0,0,
                                    (GetHeight()/2 - GetVPlaneOffset())),
                      logV,                     // logical volume
                      logV->GetName(),          // name
                      logVolume,                // mother  volume
                      false,                    // (not used)
                      0,                        // Copy number (zero)
                      Check());                 // Check overlaps.


    G4RotationMatrix* uRotation = new G4RotationMatrix();
    uRotation->rotateX(180*degree);
    uRotation->rotateZ(60*degree);

    CaptWirePlaneBuilder& uPlane = Get<CaptWirePlaneBuilder>("UPlane");
    uPlane.SetApothem(GetApothem());
    G4LogicalVolume *logU = uPlane.GetPiece();
    new G4PVPlacement(uRotation,                // rotation.
                      G4ThreeVector(0,0,
                                    (GetHeight()/2 - GetUPlaneOffset())),
                      logU,                     // logical volume
                      logU->GetName(),          // name
                      logVolume,                // mother  volume
                      false,                    // (not used)
                      0,                        // Copy number (zero)
                      Check());                 // Check overlaps.

    CaptWirePlaneBuilder& gridPlane = Get<CaptWirePlaneBuilder>("GridPlane");
    gridPlane.SetApothem(GetApothem());
    G4LogicalVolume *logGrid = gridPlane.GetPiece();
    new G4PVPlacement(xRotation,                // rotation.
                      G4ThreeVector(0,0,
                                    (GetHeight()/2 - GetGridPlaneOffset())),
                      logGrid,                     // logical volume
                      logGrid->GetName(),          // name
                      logVolume,                // mother  volume
                      false,                    // (not used)
                      0,                        // Copy number (zero)
                      Check());                 // Check overlaps.

    CaptWirePlaneBuilder& grndPlane = Get<CaptWirePlaneBuilder>("GroundPlane");
    grndPlane.SetApothem(GetApothem());
    G4LogicalVolume *logGrnd = grndPlane.GetPiece();
    new G4PVPlacement(xRotation,                // rotation.
                      G4ThreeVector(0,0,
                                    (GetHeight()/2 - GetGroundPlaneOffset())),
                      logGrnd,                     // logical volume
                      logGrnd->GetName(),          // name
                      logVolume,                // mother  volume
                      false,                    // (not used)
                      0,                        // Copy number (zero)
                      Check());                 // Check overlaps.

    return logVolume;
}

double CaptDriftRegionBuilder::GetRadius ( ) const

inline

Definition at line 41 of file CaptDriftRegionBuilder.hh.

{return fApothem/0.866025403784;}

double CaptDriftRegionBuilder::GetUPlaneOffset

(

)

Get the offset of the X plane from the top of the drift region.

Definition at line 125 of file CaptDriftRegionBuilder.cc.

                                               {
    CaptWirePlaneBuilder& wires = Get<CaptWirePlaneBuilder>("XPlane");
    return 3.0*GetWirePlaneSpacing() + wires.GetHeight()/2;
}

double CaptDriftRegionBuilder::GetVPlaneOffset

(

)

Get the offset of the X plane from the top of the drift region.

Definition at line 120 of file CaptDriftRegionBuilder.cc.

                                               {
    CaptWirePlaneBuilder& wires = Get<CaptWirePlaneBuilder>("XPlane");
    return 2.0*GetWirePlaneSpacing() + wires.GetHeight()/2;
}

double CaptDriftRegionBuilder::GetWirePlaneSpacing ( )

inline

Definition at line 53 of file CaptDriftRegionBuilder.hh.

{return fWirePlaneSpacing;}

double CaptDriftRegionBuilder::GetXPlaneOffset

(

)

Get the offset of the X plane from the top of the drift region.

Definition at line 115 of file CaptDriftRegionBuilder.cc.

                                               {
    CaptWirePlaneBuilder& wires = Get<CaptWirePlaneBuilder>("XPlane");
    return 1.0*GetWirePlaneSpacing() + wires.GetHeight()/2;
}

void CaptDriftRegionBuilder::Init ( void  )

private

Definition at line 81 of file CaptDriftRegionBuilder.cc.

                                      {
    SetMessenger(new CaptDriftRegionMessenger(this));
    SetApothem(1037.5*mm);
    SetDriftLength(1000*mm);
    SetWirePlaneSpacing(3.18*mm);
    SetSensitiveDetector("drift","segment");
    SetMaximumHitLength(1.0*mm);
    SetMaximumHitSagitta(0.25*mm);

    AddBuilder(new CaptWirePlaneBuilder("XPlane",this));
    AddBuilder(new CaptWirePlaneBuilder("VPlane",this));
    AddBuilder(new CaptWirePlaneBuilder("UPlane",this));
    AddBuilder(new CaptWirePlaneBuilder("GridPlane",this));
    AddBuilder(new CaptWirePlaneBuilder("GroundPlane",this));
}

void CaptDriftRegionBuilder::SetApothem ( double  v )

inline

Set the radius of the largest cylinder that can be contained in the drift region.

Definition at line 32 of file CaptDriftRegionBuilder.hh.

{fApothem = v;}

void CaptDriftRegionBuilder::SetDriftLength ( double  v )

inline

Set the length of the drift region from the cathode to the first wire plane. (i.e. the grid plane).

Definition at line 47 of file CaptDriftRegionBuilder.hh.

{fDriftLength = v;}

void CaptDriftRegionBuilder::SetRadius ( double  v )

inline

Set the radius of the smallest cylinder that contains the drift region. This is the maximum local Y dimension and is the apothem divided by the cosine of 30 degrees.

Definition at line 40 of file CaptDriftRegionBuilder.hh.

{fApothem = v*0.866025403784;}

void CaptDriftRegionBuilder::SetWirePlaneSpacing ( double  v )

inline

Set the spacing between the wire planes.

Definition at line 52 of file CaptDriftRegionBuilder.hh.

{fWirePlaneSpacing = v;}

Member Data Documentation

double CaptDriftRegionBuilder::fApothem

private

The radius of the circle that can fit inside.

Definition at line 85 of file CaptDriftRegionBuilder.hh.

double CaptDriftRegionBuilder::fDriftLength

private

The distance from the cathode to the first wire plane.

Definition at line 88 of file CaptDriftRegionBuilder.hh.

double CaptDriftRegionBuilder::fWirePlaneSpacing

private

The distance between the wire planes.

Definition at line 91 of file CaptDriftRegionBuilder.hh.

---

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

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