From 29a0170003789d2e127e16946c0cbec7f400a4a2 Mon Sep 17 00:00:00 2001
From: shuyaota <shuyan_timeslip@hotmail.com>
Date: Tue, 20 Feb 2018 10:03:06 -0600
Subject: [PATCH] Added Deadlayer for Hyball
---
NPSimulation/Detectors/Tiara/Tiara.cc | 18 ++++++++++++++++++
NPSimulation/Detectors/Tiara/Tiara.hh | 13 ++++++++++---
2 files changed, 28 insertions(+), 3 deletions(-)
diff --git a/NPSimulation/Detectors/Tiara/Tiara.cc b/NPSimulation/Detectors/Tiara/Tiara.cc
index f2d240f86..97b0b4128 100644
--- a/NPSimulation/Detectors/Tiara/Tiara.cc
+++ b/NPSimulation/Detectors/Tiara/Tiara.cc
@@ -719,6 +719,12 @@ void Tiara::ConstructHyball(G4LogicalVolume* world){
HYBALL_ActiveWafer_OuterRadius,0.6*HYBALL_ActiveWafer_Thickness,
-0.5*HYBALL_ActiveWafer_Angle,HYBALL_ActiveWafer_Angle);
+ //by Shuya 180219. Hyball DeadLayer
+ G4Tubs* Hyball_DeadLayer =
+ new G4Tubs("HyballDeadLayer",HYBALL_ActiveWafer_InnerRadius,
+ HYBALL_ActiveWafer_OuterRadius,0.5*HYBALL_DeadLayer_Thickness,
+ -0.5*HYBALL_ActiveWafer_Angle,HYBALL_ActiveWafer_Angle);
+
// Substract Active Wafer from Wafer
G4SubtractionSolid* InertWafer = new G4SubtractionSolid("Hyball_InertWafer", WaferFull, ActiveWaferShape,
@@ -751,6 +757,12 @@ void Tiara::ConstructHyball(G4LogicalVolume* world){
logicAW->SetVisAttributes(SiliconVisAtt);
logicAW->SetSensitiveDetector(m_HyballScorer);
+ //by Shuya 180219
+ // logic Hyball DeadLayer
+ G4LogicalVolume* logicDL =
+ new G4LogicalVolume(Hyball_DeadLayer,m_MaterialSilicon,"logicDL", 0, 0, 0);
+ logicDL->SetVisAttributes(SiliconVisAtt);
+
// Place all the Piece in the mother volume
new G4PVPlacement(new G4RotationMatrix(0,0,0),
G4ThreeVector(0,0,0),
@@ -767,6 +779,12 @@ void Tiara::ConstructHyball(G4LogicalVolume* world){
logicAW,"Hyball_ActiveWafer",
logicHyball,false,0);
+ //by Shuya 180219. Note beacause Hyball is placed in negative direction (such as -147 mm) in world coordinate, deadlayer (front side) positions must be added in positive direction.
+ new G4PVPlacement(new G4RotationMatrix(0,0,0),
+ G4ThreeVector(0,0,0.5*HYBALL_ActiveWafer_Thickness+0.5*HYBALL_DeadLayer_Thickness),
+ logicDL,"Hyball_DeadLayer",
+ logicHyball,false,0);
+
for(unsigned int i = 0 ; i < m_HyballZ.size() ; i++){
// Place mother volume
new G4PVPlacement(new G4RotationMatrix(0,0,m_HyballPhi[i]),
diff --git a/NPSimulation/Detectors/Tiara/Tiara.hh b/NPSimulation/Detectors/Tiara/Tiara.hh
index dd9544f40..e31169cd5 100644
--- a/NPSimulation/Detectors/Tiara/Tiara.hh
+++ b/NPSimulation/Detectors/Tiara/Tiara.hh
@@ -48,8 +48,10 @@ using namespace CLHEP;
namespace TIARA{
// Energy and time Resolution
const G4double ResoTime = 0 ;
- const G4double ResoEnergyInnerBarrel = 0.017*MeV ;// = 136keV FWHM
- const G4double ResoEnergyOuterBarrel = 0.017*MeV ;// = 136keV FWHM
+ //const G4double ResoEnergyInnerBarrel = 0.017*MeV ;// = 136keV FWHM
+ //const G4double ResoEnergyOuterBarrel = 0.017*MeV ;// = 136keV FWHM
+ const G4double ResoEnergyInnerBarrel = 0.050*MeV ;// = 136keV FWHM
+ const G4double ResoEnergyOuterBarrel = 0.050*MeV ;// = 136keV FWHM
const G4double ResoEnergyHyball = 0.017*MeV ;// = 70keV FWHM
const G4double EnergyThreshold = 200*keV;
@@ -86,7 +88,9 @@ namespace TIARA{
const G4double INNERBARREL_ActiveWafer_Thickness =400*um;
const G4double INNERBARREL_InertWafer_Length = 97.00*mm;
const G4double INNERBARREL_InertWafer_Width = 24.80*mm;
- const G4double INNERBARREL_ActiveWafer_DeadLayerThickness = 1*um;
+ //by Shuya 180219. Equivalent to thickness in Al (note material itself is Si).
+ //const G4double INNERBARREL_ActiveWafer_DeadLayerThickness = 1*um;
+ const G4double INNERBARREL_ActiveWafer_DeadLayerThickness = 0.3*um;
const G4double INNERBARREL_InertWafer_Thickness =
INNERBARREL_ActiveWafer_Thickness+ 2*INNERBARREL_ActiveWafer_DeadLayerThickness;
const G4int INNERBARREL_NumberOfStrip = 4;
@@ -117,6 +121,9 @@ namespace TIARA{
const G4double HYBALL_ActiveWafer_Angle = 54.8*deg;
const G4double HYBALL_NumberOfAnnularStrip = 16;
const G4double HYBALL_NumberOfRadialStrip = 8 ;
+ //by Shuya 180219. Equivalent to thickness in Al (note material itself is Si).
+ const G4double HYBALL_DeadLayer_Thickness = 0.61*um ;
+
}
--
GitLab