GaspardTrackerTrapezoid.hh 7.98 KiB
/*****************************************************************************
* 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: *
* *
* *
*****************************************************************************/
#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.021 ;// = 52keV of Resolution // Unit is MeV/2.35
const G4double ResoSecondStage = 0.043 ;// = 130 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
// gaspard Trapezoid
// 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;
// gaspard Trapezoid
// 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;
// Second stage
const G4double SecondStageHeight = FirstStageHeight;
const G4double SecondStageBaseLarge = FirstStageBaseLarge;
const G4double SecondStageBaseSmall = FirstStageBaseSmall;
const G4double SecondStageThickness = 1.5*mm;
// Third stage
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;
}
#endif