ParametricGeo.h

Go to the documentation of this file.

/*
 * ParametricGeo.h
 *
 *  Created on: 28 Mar 2018
 *      Author: Nicola Mori
 */

#ifndef PARAMETRICGEO_H_
#define PARAMETRICGEO_H_

// ParamGeo headers
#include "geometry/GGSVGeometryConstruction.h"
#include <simulation/geometry/parametric/DetStructs.h>

// Geant4 headers
#include "globals.hh"

class G4VPhysicalVolume;
class G4LogicalVolume;

// clang-format off
// clang-format on

class ParametricGeo : public GGSVGeometryConstruction {
public:
  ParametricGeo();

public:
  G4VPhysicalVolume *Construct();

  bool ExportParameters();

  bool IsInsideAcceptance(const G4ThreeVector &generationPosition, const G4ThreeVector &direction) const;

  int GetAcceptanceType(const G4ThreeVector &generationPosition, const G4ThreeVector &direction) const;

  G4VPhysicalVolume *GetVolume() { return _worldPhysical; }

  const std::string GetVersion() { return "0"; }

private:
  std::vector<double> InterceptX(double Xplane, const G4ThreeVector &generationPosition,
                                 const G4ThreeVector &direction) const;
  std::vector<double> InterceptY(double Yplane, const G4ThreeVector &generationPosition,
                                 const G4ThreeVector &direction) const;

  std::vector<double> InterceptZ(double Zplane, const G4ThreeVector &generationPosition,
                                 const G4ThreeVector &direction) const;

  double PntDistance(std::vector<double> v, const G4ThreeVector a) const;

  std::vector<bool> ParseListOfWLayers(std::string str);

  /* ***** Parameters of the geometry ***** */
  // General
  GeneralOptions _generalOpt;
  // CALO
  CaloOptions _caloOpt;
  // Acceptance
  AcceptanceOptions _acceptanceOpt;
  // PSD
  PSDOptions _psdOpt;
  // STK
  STKOptions _stkOpt;
  // FIT
  FITOptions _fitOpt;
  // SCD
  SCDOptions _scdOpt;
  // TRD
  TRDOptions _trdOpt;
  // Shield
  ShieldOptions _shieldOpt;
  // CSS
  CSSOptions _cssOpt;

  // Methods for setting geometric properties via datacard commands
  void _SetAnticoincidenceDetector(const std::string &ac);
  void _SetPsdType(const std::string &psdType);
  void _SetCaloCrystalDensity(double density);
  void _SetCaloSideFitDistance(double d);

  G4VPhysicalVolume *_worldPhysical;

  struct NormalPositions {
    std::vector<float> topLayerZPos;
    std::vector<float> sideLayerXPos;
    std::vector<float> sideLayerYPos;
  } stkNormalPositions, scdNormalPositions, psdNormalPositions, fitNormalPositions;
};

#endif /* PARAMETRICGEO_H_ */