Newer
Older
/*****************************************************************************
* Copyright (C) 2009-2010 this file is part of the NPTool Project *
* *
* For the licensing terms see $NPTOOL/Licence/NPTool_Licence *
* For the list of contributors see $NPTOOL/Licence/Contributors *
*****************************************************************************/
/*****************************************************************************
* Original Author: N. de Sereville contact address: deserevi@ipno.in2p3.fr *
* *
* Creation Date : 15/07/09 *
* Last update : *
*---------------------------------------------------------------------------*
* Decription: Define a module of trapezoidal shape for the Gaspard tracker *
* *
*---------------------------------------------------------------------------*
* Comment: *
* *
* *
*****************************************************************************/
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
#ifndef GaspardTrackerTrapezoid_h
#define GaspardTrackerTrapezoid_h 1
// C++ headers
#include <vector>
// NPTool header
#include "GaspardTrackerModule.hh"
#include "TInteractionCoordinates.h"
using namespace std;
class GaspardTrackerTrapezoid : public GaspardTrackerModule
{
////////////////////////////////////////////////////
/////// Default Constructor and Destructor /////////
////////////////////////////////////////////////////
public:
GaspardTrackerTrapezoid();
virtual ~GaspardTrackerTrapezoid();
////////////////////////////////////////////////////
//////// Specific Function of this Class ///////////
////////////////////////////////////////////////////
public:
// By Position Method
void AddModule(G4ThreeVector TL ,
G4ThreeVector BL ,
G4ThreeVector BR ,
G4ThreeVector CT ,
bool wFirstStage ,
bool wSecondStage ,
bool wThirdStage);
// By Angle Method
void AddModule(G4double R ,
G4double Theta ,
G4double Phi ,
G4double beta_u ,
G4double beta_v ,
G4double beta_w ,
bool wFirstStage ,
bool wSecondStage ,
bool wThirdStage);
// Effectively construct Volume
// Avoid to have two time same code for Angle and Point definition
void VolumeMaker(G4int TelescopeNumber ,
G4ThreeVector MMpos ,
G4RotationMatrix* MMrot ,
bool wFirstStage ,
bool wSecondStage ,
bool wThirdStage ,
G4LogicalVolume* world);
////////////////////////////////////////////////////
//// Inherite from GaspardTrackerModule class /////
////////////////////////////////////////////////////
public:
// Read stream at Configfile to pick-up parameters of detector (Position,...)
// Called in DetecorConstruction::ReadDetextorConfiguration Method
void ReadConfiguration(string Path);
// Construct detector and inialise sensitive part.
// Called After DetecorConstruction::AddDetector Method
void ConstructDetector(G4LogicalVolume* world);
// Add Detector branch to the EventTree.
// Called After DetecorConstruction::AddDetector Method
void InitializeRootOutput();
// Initialize all scorers necessary for the detector
void InitializeScorers();
// Read sensitive part and fill the Root tree.
// Called at in the EventAction::EndOfEventAvtion
void ReadSensitive(const G4Event* event);
// Give the static TInteractionCoordinates from VDetector to the classes
// deriving from GaspardTrackerModule
// This is mandatory since the GaspardTracker*** does not derive from VDetector
void SetInterCoordPointer(TInteractionCoordinates* interCoord);
TInteractionCoordinates* GetInterCoordPointer() {return ms_InterCoord;};
////////////////////////////////////////////////////
///////////////Private intern Data//////////////////
////////////////////////////////////////////////////
private:
// Interaction Coordinates coming from VDetector through the
// SetInteractionCoordinatesPointer method
TInteractionCoordinates* ms_InterCoord;
// True if Define by Position, False is Define by angle
vector<bool> m_DefinitionType ;
// Used for "By Point Definition"
vector<G4ThreeVector> m_X1_Y1 ; // Top Left Corner Position Vector
vector<G4ThreeVector> m_X1_Y128 ; // Bottom Left Corner Position Vector
vector<G4ThreeVector> m_X128_Y1 ; // Bottom Right Corner Position Vector
vector<G4ThreeVector> m_X128_Y128 ; // Center Corner Position Vector
// Used for "By Angle Definition"
vector<G4double> m_R ; // |
vector<G4double> m_Theta ; // > Spherical coordinate of Strips Silicium Plate
vector<G4double> m_Phi ; // |
vector<G4double> m_beta_u ; // |
vector<G4double> m_beta_v ; // > Tilt angle of the Telescope
vector<G4double> m_beta_w ; // |
};
namespace GPDTRAP
{
// Energy/Time resolutions for the different layers
const G4double ResoFirstStage = 0 ;// = 52keV of Resolution // Unit is MeV/2.35
// const G4double ResoFirstStage = 0.022 ;// = 52keV of Resolution // Unit is MeV/2.35
const G4double ResoSecondStage = 0.055 ;// = 130 keV of resolution // Unit is MeV/2.35
const G4double ResoThirdStage = 0 ;// = 100 keV of resolution // Unit is MeV/2.35
// const G4double ResoThirdStage = 0.043 ;// = 100 kev of resolution // Unit is MeV/2.35
const G4double ResoTimeGpd = 0.212765957 ;// = 500ps // Unit is ns/2.35
const G4double ResoTimePPAC = 0.106382979 ;// = 250ps // Unit is ns/2.35
// Geometry
// const G4double Height = 90*mm;
// const G4double BaseLarge = 85*mm;
// const G4double BaseSmall = 35*mm;
// const G4double Length = 20*mm;
// gaspHyde
const G4double Height = 115*mm;
const G4double BaseLarge = 98*mm;
const G4double BaseSmall = 30*mm;
const G4double Length = 20*mm;
const G4double InterStageDistance = 7*mm;
// First stage
const G4double FirstStageThickness = 300*micrometer;
const G4int NumberOfStripsX = 128;
const G4int NumberOfStripsY = 128;
// const G4double FirstStageHeight = 81*mm;
// const G4double FirstStageBaseLarge = 78*mm;
// const G4double FirstStageBaseSmall = 29*mm;
// gaspHyde
const G4double FirstStageHeight = 113.5*mm;
const G4double FirstStageBaseLarge = 97.5*mm;
const G4double FirstStageBaseSmall = 28.5*mm;
const G4double SecondStageHeight = FirstStageHeight;
const G4double SecondStageBaseLarge = FirstStageBaseLarge;
const G4double SecondStageBaseSmall = FirstStageBaseSmall;
const G4double SecondStageThickness = 1.5*mm;
const G4double ThirdStageHeight = FirstStageHeight;
const G4double ThirdStageBaseLarge = FirstStageBaseLarge;
const G4double ThirdStageBaseSmall = FirstStageBaseSmall;
const G4double ThirdStageThickness = 1.5*mm;
// Starting at the front of the first stage and pointing to the third stage
const G4double FirstStage_PosZ = Length* -0.5 + 0.5*FirstStageThickness;
const G4double SecondStage_PosZ = Length* -0.5 + 0.5*SecondStageThickness + 1*InterStageDistance;
const G4double ThirdStage_PosZ = Length* -0.5 + 0.5*ThirdStageThickness + 2*InterStageDistance;