CLAWSCreator.cc

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2010-2015 Belle II Collaboration                          *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Martin Ritter, Igal Jaegle                               *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#include <beast/claws/geometry/CLAWSCreator.h>
#include <beast/claws/simulation/SensitiveDetector.h>
 
#include <geometry/Materials.h>
#include <geometry/CreatorFactory.h>
#include <framework/gearbox/GearDir.h>
#include <framework/gearbox/Unit.h>
 
#include <G4LogicalVolume.hh>
#include <G4PVPlacement.hh>
#include <G4Box.hh>
#include <G4Tubs.hh>
#include <G4UserLimits.hh>
 
namespace Belle2 {
  namespace claws {
 
    // Register the creator
    geometry::CreatorFactory<CLAWSCreator> CLAWSFactory("CLAWSCreator");
 
    CLAWSCreator::CLAWSCreator(Simulation::SensitiveDetectorBase* sensitive): m_sensitive(sensitive)
    {
      if (!m_sensitive) {
        m_sensitive = new SensitiveDetector();
      }
    }
 
    CLAWSCreator::~CLAWSCreator()
    {
      delete m_sensitive;
      delete m_stepLength;
    }
 
    void CLAWSCreator::createShape(const std::string& prefix, const GearDir& params, G4LogicalVolume* parent, double roffset,
                                   bool check)
    {
      std::string name = params.getString("@name");
      if (!prefix.empty()) {
        name = prefix + "." + name;
      }
      const std::string type = params.getString("@type");
      const std::string material = params.getString("material", "");
      std::vector<double> ri = params.getArray("r", {0});
      int phase = 0;
      if (ri.size() > 1) phase = 1;
 
      //const double r = params.getLength("r", 0) / Unit::mm * CLHEP::mm;
      const double top = params.getLength("top", 0) / Unit::mm * CLHEP::mm;
      const double u = params.getLength("u", 0) / Unit::mm * CLHEP::mm;
      const bool active = params.getBool("active", false);
 
      G4Material* mat = m_topMaterial;
      if (!material.empty()) {
        mat = geometry::Materials::get(material);
      }
      G4VSolid* shape{nullptr};
      double height{0};
      if (type == "box") {
        const double length = params.getLength("length") / Unit::mm * CLHEP::mm;
        const double width = params.getLength("width") / Unit::mm * CLHEP::mm;
        height = params.getLength("height") / Unit::mm * CLHEP::mm;
        shape = new G4Box(name, width / 2, height / 2, length / 2);
      } else if (type == "tube") {
        const double length = params.getLength("length") / Unit::mm * CLHEP::mm;
        height = params.getLength("diameter") / Unit::mm * CLHEP::mm;
        shape = new G4Tubs(name, 0, height / 2, length / 2, 0, 2 * M_PI);
      }
      G4LogicalVolume* volume = new G4LogicalVolume(shape, mat, name);
      if (active) {
        volume->SetSensitiveDetector(m_sensitive);
        if (m_stepLength) {
          volume->SetUserLimits(m_stepLength);
        }
      }
      for (const GearDir& child : params.getNodes("shape")) {
        createShape(name, child, volume, height / 2, true);
      }
 
      int copyNo = 1;
 
      if (!phase) {
        double center = ri[0] / Unit::mm * CLHEP::mm  + roffset - height / 2 - top;
        //center *=CLHEP::mm / Unit::mm;
        for (double phi : params.getArray("phi", {M_PI / 2})) {
          for (double z : params.getArray("z", {0})) {
            z *= CLHEP::mm / Unit::mm;
            G4Transform3D transform = G4RotateZ3D(phi - M_PI / 2) * G4Translate3D(u, center, z);
            new G4PVPlacement(transform, volume, name, parent, false, copyNo++, check);
          }
        }
      } else {
        int i = 0;
        std::vector<double> alpha = params.getArray("alpha", {0});
        for (double z : params.getArray("z", {0})) {
          double center = ri[i] / Unit::mm * CLHEP::mm + roffset - height / 2 - top;
          //    center *= CLHEP::mm / Unit::mm;
          z *= CLHEP::mm / Unit::mm;
          for (double phi : params.getArray("phi", {M_PI / 2})) {
            G4Transform3D transform = G4RotateZ3D(phi - M_PI / 2) * G4Translate3D(u, center, z) * G4RotateY3D(alpha[i]);
            new G4PVPlacement(transform, volume, name, parent, false, copyNo++, check);
          }
          i++;
        }
 
      }
    }
 
 
 
    void CLAWSCreator::create(const GearDir& content, G4LogicalVolume& topVolume, geometry::GeometryTypes /* type */)
    {
      m_topMaterial = topVolume.GetMaterial();
      double stepLength = content.getLength("stepLength", 0) / Unit::mm * CLHEP::mm;
      if (stepLength > 0) {
        if (m_stepLength) delete m_stepLength;
        m_stepLength = new G4UserLimits(stepLength);
      }
      for (auto shape : content.getNodes("shape")) {
        createShape("", shape, &topVolume, 0, false);
      }
    }
  } // claws namespace
} // Belle2 namespace