ClawCreator Class Reference

The creator for the CLAW geometry. More...

#include <ClawCreator.h>

Inheritance diagram for ClawCreator:

Collaboration diagram for ClawCreator:

Public Member Functions
	ClawCreator ()
	Constructor.
virtual	~ClawCreator ()
	Destructor.
virtual void	create (const GearDir &content, G4LogicalVolume &topVolume, geometry::GeometryTypes type)
	Creation of the detector geometry from Gearbox (XML). More...
	BELLE2_DEFINE_EXCEPTION (DBNotImplemented, "Cannot create geometry from Database.")
	Exception that will be thrown in createFromDB if member is not yet implemented by creator.
virtual void	createFromDB (const std::string &name, G4LogicalVolume &topVolume, GeometryTypes type)
	Function to create the geometry from the Database. More...
virtual void	createPayloads (const GearDir &content, const IntervalOfValidity &iov)
	Function to create the geometry database. More...

Protected Attributes
SensitiveDetector *	m_sensitive
	SensitiveDetector CLAW.

Detailed Description

The creator for the CLAW geometry.

Definition at line 27 of file ClawCreator.h.

Member Function Documentation

create()

void create ( const GearDir &  content, G4LogicalVolume &  topVolume, geometry::GeometryTypes  type  )

virtual

Creation of the detector geometry from Gearbox (XML).

Parameters

[in]	content	XML data directory.
[in]	topVolume	Geant world volume.
[in]	type	Geometry type.

Implements CreatorBase.

Definition at line 58 of file ClawCreator.cc.

    {
 
      m_sensitive = new SensitiveDetector();
 
      //lets get the stepsize parameter with a default value of 5 µm
      double stepSize = content.getLength("stepSize", 5 * CLHEP::um);
 
      G4VisAttributes* red = new G4VisAttributes(G4Colour(1, 0, 0));
      red->SetForceAuxEdgeVisible(true);
      G4VisAttributes* green = new G4VisAttributes(G4Colour(0, 1, 0));
      green->SetForceAuxEdgeVisible(true);
      G4VisAttributes* gray = new G4VisAttributes(G4Colour(.5, .5, .5));
      gray->SetForceAuxEdgeVisible(true);
 
      G4VisAttributes* coppercolor = new G4VisAttributes(G4Colour(218. / 255., 138. / 255., 103. / 255.));
      coppercolor->SetForceAuxEdgeVisible(true);
 
      int detID = 0;
      //Lets loop over all the Active nodes
      BOOST_FOREACH(const GearDir & activeParams, content.getNodes("Active")) {
        /*
              double x_pos[100];
              double y_pos[100];
              double x_off[100];
              double y_off[100];
              double z_pos[100];
              double phi[100];
              double thetaZ[100];
              double r[100];
              int dimx = 0;
              int dimy = 0;
              int dimx_offset = 0;
              int dimy_offset = 0;
              int dimz = 0;
              int dimr_dia = 0;
              for (double x_offset : activeParams.getArray("x_offset", {0})) {
                x_offset *= CLHEP::cm;
                x_off[dimx_offset] = x_offset;
                dimx_offset++;
              }
              for (double y_offset : activeParams.getArray("y_offset", {0})) {
                y_offset *= CLHEP::cm;
                y_off[dimy_offset] = y_offset;
                dimy_offset++;
              }
              for (double x : activeParams.getArray("x", {0})) {
                x *= CLHEP::cm;
                x_pos[dimx] = x + x_off[dimx];
                dimx++;
              }
              for (double y : activeParams.getArray("y", {0})) {
                y *= CLHEP::cm;
                y_pos[dimy] = y + y_off[dimy];
                dimy++;
              }
              for (double z : activeParams.getArray("z", {0})) {
                z *= CLHEP::cm;
                z_pos[dimz] = z;
                if (dimx != 0 && dimy != 0)
                  r[dimz] = sqrt(x_pos[dimz] * x_pos[dimz] + y_pos[dimz] * y_pos[dimz]);
                dimz++;
              }
              int dimPhi = 0;
              for (double Phi : activeParams.getArray("Phi", {0})) {
                phi[dimPhi] = Phi;
                dimPhi++;
              }
              int dimThetaZ = 0;
              for (double ThetaZ : activeParams.getArray("ThetaZ", {0})) {
                thetaZ[dimThetaZ] = ThetaZ;
                dimThetaZ++;
              }
              if (dimx == 0 && dimy == 0) {
                for (double r_dia : activeParams.getArray("r_dia", {0})) {
                  r_dia *= CLHEP::cm;
                  r[dimr_dia] = r_dia;
                  dimr_dia++;
                }
              }
        */
        int Nscint = activeParams.getInt("Nscint");
        G4double x = activeParams.getLength("x") * CLHEP::cm;
        G4double y = activeParams.getLength("y") * CLHEP::cm;
        G4double z = activeParams.getLength("z_claw") * CLHEP::cm;
        /*G4double r = 0;
        if (x == 0 && y == 0)
          r = activeParams.getLength("r_claw") * CLHEP::cm;
        else
        r = sqrt(x * x + y * y);*/
        //G4double phi = activeParams.getAngle("Phi");
        //G4double thetaZ = activeParams.getAngle("ThetaZ");
        G4double dx_board = activeParams.getLength("dx_board") / 2.*CLHEP::cm;
        G4double dy_board = activeParams.getLength("dy_board") / 2.*CLHEP::cm;
        G4double dz_board = activeParams.getLength("dz_board") / 2.*CLHEP::cm;
        G4double dz_Culayer =  activeParams.getLength("dz_Culayer") / 2.*CLHEP::cm;
        G4double dx_scint = activeParams.getLength("dx_scint") / 2.*CLHEP::cm;
        G4double dy_scint = activeParams.getLength("dy_scint") / 2.*CLHEP::cm;
        G4double dz_scint = activeParams.getLength("dz_scint") / 2.*CLHEP::cm;
        G4double Al_width = activeParams.getLength("Al_width") / 2.*CLHEP::cm;
        G4double dx_Al = dx_scint + Al_width;
        G4double dy_Al = dy_scint + Al_width;
        G4double dz_Al = dz_scint + Al_width;
        G4double dx_air = dx_Al;
        G4double dy_air = dy_Al * Nscint;
        G4double dz_air = dz_board + 2. * dz_Culayer + dz_Al;
        //cout << " dx sol 1 " << dx_air << " sol 2 " << dx_board << endl;
        //cout << " dy sol 1 " << dy_air << " sol 2 " << dy_board << endl;
 
        //create air volume engloging ladder + scintillator + Al foil
        G4Box* s_air = new G4Box("s_air", dx_air, dy_air, dz_air);
        G4LogicalVolume* l_air = new G4LogicalVolume(s_air, geometry::Materials::get("G4_AIR"), "l_air");
 
 
        //create claw G10 board
        G4Box* s_board = new G4Box("s_board", dx_board, dy_board, dz_board);
        G4LogicalVolume* l_board = new G4LogicalVolume(s_board, geometry::Materials::get("G10"), "l_board");
        l_board->SetVisAttributes(green);
 
        G4double r_board = dz_air - 2. * dz_Culayer - dz_board;
 
        //create copper layer on each side of the G10 board
        G4Box* s_Culayer = new G4Box("s_Culayer", dx_board, dy_board, dz_Culayer);
        G4LogicalVolume* l_Culayer = new G4LogicalVolume(s_Culayer, geometry::Materials::get("Cu"), "l_Culayer");
        l_Culayer->SetVisAttributes(coppercolor);
 
        G4double r_Culayer_bot = r_board - dz_board - dz_Culayer;
        G4double r_Culayer_top = r_board + dz_board + dz_Culayer;
        //cout << " cu bot sol 1 " << r_Culayer_bot << " sol 2 " << -dz_air + dz_Culayer<< endl;
 
        //create scintillator and Al foil around it
        G4VSolid* s_scint = new G4Box("s_scint", dx_scint, dy_scint, dz_scint);
        G4VSolid* s_Al = new G4Box("s_Al", dx_Al, dy_Al, dz_Al);
        s_Al = new G4SubtractionSolid("s_Al", s_Al, s_scint, 0, G4ThreeVector(0, 0, 0));
        G4LogicalVolume* l_Al = new G4LogicalVolume(s_Al, geometry::Materials::get("Al"), "l_Al");
        l_Al->SetVisAttributes(gray);
 
        G4LogicalVolume* l_scint = new G4LogicalVolume(s_scint, geometry::Materials::get("G4_POLYSTYRENE"), "l_scint", 0, m_sensitive);
        l_scint->SetVisAttributes(red);
        //Lets limit the Geant4 stepsize inside the volume
        l_scint->SetUserLimits(new G4UserLimits(stepSize));
 
        G4double r_Al = r_Culayer_bot - dz_Culayer - dz_Al;
        //cout << " dz_air " << dz_air << " x2 " << 2. * dz_air << " r_Al " << r_Al << " tot " << r_Al + dz_Al << " sol3 " << r_Culayer_top + dz_Culayer + dz_Al << endl;
        double z_0 = -dy_air + dy_Al;
 
        //for(int i = 0; i < dimz; i++) {
 
        //G4Transform3D transform = G4RotateZ3D(-M_PI / 2 + phi[i]) * G4Translate3D(0, r[i], z_pos[i]) * G4RotateX3D(-M_PI / 2 - thetaZ[i]);
        //G4Transform3D transform = G4RotateZ3D(-M_PI / 2 + phi) * G4Translate3D(0, r, z) * G4RotateX3D(-M_PI / 2 - thetaZ);
        G4Transform3D transform = G4Translate3D(x, y, z);
        new G4PVPlacement(transform, l_air, "p_air", &topVolume, false, 1);
        new G4PVPlacement(0, G4ThreeVector(0, 0, r_board), l_board, "p_board", l_air, false, 1);
        new G4PVPlacement(0, G4ThreeVector(0, 0, r_Culayer_bot), l_Culayer, "p_Culayer_bot", l_air, false, 1);
        new G4PVPlacement(0, G4ThreeVector(0, 0, r_Culayer_top), l_Culayer, "p_Culayer_top", l_air, false, 1);
 
        for (int j = 0; j < Nscint; j++) {
          double i_z = z_0 + j * 2. * dy_Al;
 
          new G4PVPlacement(0, G4ThreeVector(0, i_z, r_Al), l_Al, "p_Al", l_air, false, 1);
          new G4PVPlacement(0, G4ThreeVector(0, i_z, r_Al), l_scint, "p_scint", l_air, false, detID);
          B2INFO("Phase1-CLAWS-" << detID << " placed at: " << transform.getTranslation() << " mm");
          detID++;
          //}
        }
      }
    }

createFromDB()

void createFromDB ( const std::string &  name, G4LogicalVolume &  topVolume, GeometryTypes  type  )

virtualinherited

Function to create the geometry from the Database.

Parameters

name	name of the component in the database, could be used to disambiguate multiple components created with the same creator
topVolume	Top volume in which the geometry has to be placed
type	Type of geometry to be build

Reimplemented in GeoVXDServiceCreator, GeoHeavyMetalShieldCreator, GeoTOPCreator, GeoSVDCreator, GeoSTRCreator, GeoServiceMaterialCreator, GeoCOILCreator, GeoPXDCreator, GeoKLMCreator, GeoEKLMCreator, GeoBKLMCreator, GeoFarBeamLineCreator, GeoCryostatCreator, GeoBeamPipeCreator, MyDBCreator, GeoECLCreator, GeoCDCCreatorReducedCDC, GeoCDCCreator, BeamabortCreator, GeoARICHCreator, and GeoMagneticField.

Definition at line 17 of file CreatorBase.cc.

createPayloads()

void createPayloads ( const GearDir &  content, const IntervalOfValidity &  iov  )

virtualinherited

Function to create the geometry database.

This function should be implemented to convert Gearbox parameters to one ore more database payloads

Parameters

content	GearDir pointing to the parameters which should be used for construction
iov	interval of validity to use when generating payloads

Reimplemented in GeoVXDServiceCreator, GeoHeavyMetalShieldCreator, GeoTOPCreator, GeoSVDCreator, GeoSTRCreator, GeoServiceMaterialCreator, GeoCOILCreator, GeoPXDCreator, GeoKLMCreator, GeoEKLMCreator, GeoBKLMCreator, GeoFarBeamLineCreator, GeoCryostatCreator, GeoBeamPipeCreator, MyDBCreator, GeoMagneticField, GeoECLCreator, GeoCDCCreatorReducedCDC, GeoCDCCreator, BeamabortCreator, and GeoARICHCreator.

Definition at line 24 of file CreatorBase.cc.

---

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

• beast/claw/geometry/include/ClawCreator.h
• beast/claw/geometry/src/ClawCreator.cc