From f5a9c031aff1afe936735a6bfffa0def0204b897 Mon Sep 17 00:00:00 2001
From: Gitea <gitea@fake.local>
Date: Mon, 2 Apr 2018 23:16:42 +0200
Subject: [PATCH] Basic DetectorConstruction

---
 include/MedtechDetectorConstruction.hh | 35 ++++++++++++
 src/MedtechDetectorConstruction.cc     | 79 ++++++++++++++++++++++++++
 2 files changed, 114 insertions(+)
 create mode 100644 include/MedtechDetectorConstruction.hh
 create mode 100644 src/MedtechDetectorConstruction.cc

diff --git a/include/MedtechDetectorConstruction.hh b/include/MedtechDetectorConstruction.hh
new file mode 100644
index 0000000..f360cb3
--- /dev/null
+++ b/include/MedtechDetectorConstruction.hh
@@ -0,0 +1,35 @@
+//
+//  MedtechDetectorConstruction.hh
+//  medtech
+//
+//  Created by Baranyai David on 2018. 04. 01..
+//
+
+#ifndef MedtechDetectorConstruction_hh
+#define MedtechDetectorConstruction_hh
+
+#include <stdio.h>
+#include "G4VUserDetectorConstruction.hh"
+#include "G4Box.hh"
+#include "globals.hh"
+#include "G4SystemOfUnits.hh"
+#include "G4LogicalVolume.hh"
+#include "G4PVPlacement.hh"
+#include "G4ThreeVector.hh"
+#include "G4NistManager.hh"
+#include "G4Material.hh"
+
+class MedtechDetectorConstruction : public G4VUserDetectorConstruction
+{
+public:
+    MedtechDetectorConstruction();
+    virtual ~MedtechDetectorConstruction();
+    virtual G4VPhysicalVolume* Construct();
+    
+    G4LogicalVolume* GetScoringVolume() const { return fScoringVolume; }
+    
+protected:
+    G4LogicalVolume*  fScoringVolume;
+};
+
+#endif /* MedtechDetectorConstruction_hh */
diff --git a/src/MedtechDetectorConstruction.cc b/src/MedtechDetectorConstruction.cc
new file mode 100644
index 0000000..645aae3
--- /dev/null
+++ b/src/MedtechDetectorConstruction.cc
@@ -0,0 +1,79 @@
+//
+//  MedtechDetectorConstruction.cc
+//  medtech
+//
+//  Created by Baranyai David on 2018. 04. 01..
+//
+
+#include "MedtechDetectorConstruction.hh"
+
+MedtechDetectorConstruction::MedtechDetectorConstruction() : G4VUserDetectorConstruction()
+{
+    
+}
+
+MedtechDetectorConstruction::~MedtechDetectorConstruction()
+{
+    
+}
+
+G4VPhysicalVolume* MedtechDetectorConstruction::Construct()
+{
+    //Material Section
+    G4NistManager *manager = G4NistManager::Instance();
+    G4Material *air = manager -> FindOrBuildMaterial("G4_AIR");
+    G4Material *targetMaterial = manager -> FindOrBuildMaterial("G4_Pb");
+    G4Material *collMaterial = manager -> FindOrBuildMaterial("G4_W");
+    //End of Material Section
+    
+    //Physical world
+    G4int world_xyz = 1*m; // world is 1m x 1m x 2m for each
+    G4Box *solidWorld = new G4Box("World",0.5*world_xyz,0.5*world_xyz,4*m); //define a box for world
+    
+    G4LogicalVolume *logicWorld = new G4LogicalVolume(solidWorld, air, "World", 0, 0, 0); //define a logical volume for world
+    
+    G4PVPlacement *physicalWorld = new G4PVPlacement(0, //no rotation
+                                                     G4ThreeVector(), //at (0,0,0)
+                                                     logicWorld, //it's logical volume
+                                                     "World", //it's name
+                                                     0, //it's mother volume
+                                                     false, //no boolean operations
+                                                     0); //no magnetic field
+    //End of Physical world
+    
+    //Target box
+    G4Box *solidTarget = new G4Box("Target", 1*m, 1*m, 0.1*cm);
+    
+    G4LogicalVolume *logicTarget = new G4LogicalVolume(solidTarget, targetMaterial, "Collimator", 0, 0, 0);
+    
+    G4PVPlacement *physicalTarget = new G4PVPlacement(0, //no rotation
+                                                      G4ThreeVector(0, 0, 4*m), //at (0,0,0)
+                                                      logicTarget, //it's logical volume
+                                                      "Target", //it's name
+                                                      logicWorld, //it's mother volume
+                                                      false, //no boolean operations
+                                                      0); //no particular field
+    //End of Target box
+    
+    //Collimator
+    G4Box *solidColl = new G4Box("Collimator 1", 0.25*m, 0.25*m, 0.1*cm);
+    G4LogicalVolume *logicColl = new G4LogicalVolume(solidColl, collMaterial, "Collimator 1", 0, 0, 0);
+    G4PVPlacement *physicalColl = new G4PVPlacement(0, //no rotation
+                                                      G4ThreeVector(0.5*m, 0, 2*m), //at (0,0,0)
+                                                      logicColl, //it's logical volume
+                                                      "Collimator 1", //it's name
+                                                      logicWorld, //it's mother volume
+                                                      false, //no boolean operations
+                                                      0); //no particular field
+    
+    G4LogicalVolume *logicColl2 = new G4LogicalVolume(solidColl, collMaterial, "Collimator 2", 0, 0, 0);
+    G4PVPlacement *physicalColl2 = new G4PVPlacement(0, //no rotation
+                                                    G4ThreeVector(-0.5*m, 0, 2*m), //at (0,0,0)
+                                                    logicColl2, //it's logical volume
+                                                    "Collimator 2", //it's name
+                                                    logicWorld, //it's mother volume
+                                                    false, //no boolean operations
+                                                    0); //no particular field
+    
+    return physicalWorld;
+}