Skip to content
Snippets Groups Projects
Commit 977bda7b authored by Pierre Morfouace's avatar Pierre Morfouace
Browse files

Adding Hira classes for simulation

parent b8a71b32
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
Inputs/DetectorConfiguration/hira.detector 0 → 100755
+ 120
− 0
View file @ 977bda7b
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
GeneralTarget
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Thickness in micrometer
% Radius in mm
% Temperature in K, Pressure in bar
Target
THICKNESS= 0.2
ANGLE= 0
RADIUS= 5
MATERIAL= CD2
X= 0
Y= 0
Z= 0
%%%%%%%%%%%%%%%%%%%%%
HIRAArray
%%%%%%%%%%%%%%%%%%%%%
%%%%%%% Telescope 1 %%%%%%%
HiraTelescope
A= 173.37 146.75 262.67
B= 131.14 113.17 299.62
C= 151.53 79.58 301.68
D= 193.66 113.27 265.03
%%%%%%% Telescope 2 %%%%%%%
HiraTelescope
A= 185.54 74.46 284.21
B= 139.49 39.19 313.95
C= 155.91 3.85 309.52
D= 202.07 39.18 280.08
%%%%%%% Telescope 3 %%%%%%%
HiraTelescope
A= 189.04 -1.58 291.62
B= 139.10 -36.55 314.66
C= 151.29 -71.91 302.98
D= 201.44 -36.97 280.35
%%%%%%% Telescope 4 %%%%%%%
HiraTelescope
A= 183.90 -77.07 283.93
B= 130.43 -110.63 300.17
C= 139.05 -144.12 281.73
D= 192.75 -110.68 265.85
%%%%%%% Telescope 5 %%%%%%%
HiraTelescope
A= 89.03 146.90 300.27
B= 37.97 113.38 323.02
C= 56.43 79.76 330.94
D= 107.59 113.39 308.52
%%%%%%% Telescope 6 %%%%%%%
HiraTelescope
A= 94.01 74.21 324.51
B= 41.45 38.75 339.56
C= 58.48 34.39 339.75
D= 111.12 38.96 325.13
%%%%%%% Telescope 7 %%%%%%%
HiraTelescope
A= 95.57 -1.96 332.31
B= 41.29 -37.33 339.64
C= 56.63 -72.62 332.23
D= 110.97 -37.25 325.09
%%%%%%% Telescope 8 %%%%%%%
HiraTelescope
A= 93.18 -77.88 323.56
B= 37.45 -111.64 323.69
C= 51.20 -145.11 308.69
D= 107.05 -111.48 308.92
%%%%%%% Telescope 9 %%%%%%%
HiraTelescope
A= -22.26 91.79 332.36
B= -77.07 56.46 330.57
C= -61.34 21.44 338.41
D= -6.58 56.83 340.66
%%%%%%% Telescope 10 %%%%%%%
HiraTelescope
A= -23.10 -19.16 344.31
B= -76.45 -54.13 330.49
C= -59.35 -89.22 327.12
D= -5.97 -54.26 341.22
%%%%%%% Telescope 11 %%%%%%%
HiraTelescope
A= -96.45 147.32 297.42
B= -151.2 113.36 286.61
C= -140.2 79.81 303.70
D= -85.74 113.90 314.86
%%%%%%% Telescope 12 %%%%%%%
HiraTelescope
A= -105.62 74.17 320.23
B= -157.09 38.47 302.22
C= -143.26 3.23 312.38
D= -91.94 38.96 330.70
%%%%%%% Telescope 13 %%%%%%%
HiraTelescope
A= -108.60 -1.94 327.52
B= -156.99 -37.61 302.29
C= -140.04 -72.81 305.17
D= -91.79 -37.71 330.80
%%%%%%% Telescope 14 %%%%%%%
HiraTelescope
A= -104.76 -78.01 319.19
B= -150.28 -111.80 287.13
C= -130.29 -145.22 282.83
D= -84.95 -111.59 315.32
NPLib/HiRA/CMakeLists.txt
+ 5
− 4
View file @ 977bda7b
add_custom_command(OUTPUT THiraDataDict.cxx THiraDataDict_rdict.pcm THiraData.rootmap COMMAND ../scripts/build_dict.sh THiraData.h THiraDataDict.cxx THiraData.rootmap libNPHiRA.dylib DEPENDS THiraData.h)
add_library(NPHiRA SHARED THiraData.cxx THiraDataDict.cxx )
target_link_libraries(NPHiRA ${ROOT_LIBRARIES} NPCore)
install(FILES THiraData.h DESTINATION ${CMAKE_INCLUDE_OUTPUT_DIRECTORY})
add_custom_command(OUTPUT THiraPhysicsDict.cxx THiraPhysicsDict_rdict.pcm THiraPhysics.rootmap COMMAND ../scripts/build_dict.sh THiraPhysics.h THiraPhysicsDict.cxx THiraPhysics.rootmap libNPHira.dylib DEPENDS THiraPhysics.h)
add_custom_command(OUTPUT THiraDataDict.cxx THiraDataDict_rdict.pcm THiraData.rootmap COMMAND ../scripts/build_dict.sh THiraData.h THiraDataDict.cxx THiraData.rootmap libNPHira.dylib DEPENDS THiraData.h)
add_library(NPHira SHARED THiraData.cxx THiraPhysics.cxx THiraDataDict.cxx THiraPhysicsDict.cxx )
target_link_libraries(NPHira ${ROOT_LIBRARIES} NPCore)
install(FILES THiraData.h THiraPhysics.h DESTINATION ${CMAKE_INCLUDE_OUTPUT_DIRECTORY})
NPLib/HiRA/THiraPhysics.cxx 0 → 100644
+ 314
− 0
View file @ 977bda7b
/*****************************************************************************
* Copyright (C) 2009-2013 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: Adrien MATTA contact address: matta@ipno.in2p3.fr *
* *
* Creation Date : november 2009 *
* Last update : *
*---------------------------------------------------------------------------*
* Decription: *
* This class hold Hira Physics *
* *
*---------------------------------------------------------------------------*
* Comment: *
* *
* *
*****************************************************************************/
// NPL
#include "THiraPhysics.h"
#include "RootInput.h"
#include "RootOutput.h"
#include "NPDetectorFactory.h"
#include "NPCalibrationManager.h"
// STL
#include <iostream>
#include <sstream>
#include <fstream>
#include <limits>
#include <stdlib.h>
using namespace std;
// ROOT
#include "TChain.h"
// tranform an integer to a string
string itoa(int value)
{
char buffer [33];
sprintf(buffer,"%d",value);
return buffer;
}
ClassImp(THiraPhysics)
///////////////////////////////////////////////////////////////////////////
THiraPhysics::THiraPhysics()
{
NumberOfDetector = 0 ;
EventData = new THiraData ;
EventPhysics = this ;
}
///////////////////////////////////////////////////////////////////////////
THiraPhysics::~THiraPhysics()
{}
///////////////////////////////////////////////////////////////////////////
void THiraPhysics::Clear()
{
}
///////////////////////////////////////////////////////////////////////////
void THiraPhysics::ReadConfiguration(string Path)
{
ifstream ConfigFile ;
ConfigFile.open(Path.c_str()) ;
string LineBuffer ;
string DataBuffer ;
bool check_Theta = false ;
bool check_Phi = false ;
bool check_R = false ;
bool check_Thickness = false ;
bool check_Radius = false ;
bool check_LeadThickness = false ;
bool check_Scintillator = false ;
bool check_Height = false ;
bool check_Width = false ;
bool check_Shape = false ;
bool check_X = false ;
bool check_Y = false ;
bool check_Z = false ;
bool ReadingStatus = false ;
while (!ConfigFile.eof())
{
getline(ConfigFile, LineBuffer);
// If line is a Start Up Hira bloc, Reading toggle to true
if (LineBuffer.compare(0, 3, "Hira") == 0)
{
cout << "///" << endl ;
cout << "Platic found: " << endl ;
ReadingStatus = true ;
}
// Else don't toggle to Reading Block Status
else ReadingStatus = false ;
// Reading Block
while(ReadingStatus)
{
// Pickup Next Word
ConfigFile >> DataBuffer ;
// Comment Line
if (DataBuffer.compare(0, 1, "%") == 0) { ConfigFile.ignore ( std::numeric_limits<std::streamsize>::max(), '\n' );}
// Finding another telescope (safety), toggle out
else if (DataBuffer.compare(0, 6, "Hira") == 0) {
cout << "WARNING: Another Detector is find before standard sequence of Token, Error may occured in Telecope definition" << endl ;
ReadingStatus = false ;
}
//Angle method
else if (DataBuffer=="THETA=") {
check_Theta = true;
ConfigFile >> DataBuffer ;
cout << "Theta: " << atof(DataBuffer.c_str()) << "deg" << endl;
}
else if (DataBuffer=="PHI=") {
check_Phi = true;
ConfigFile >> DataBuffer ;
cout << "Phi: " << atof( DataBuffer.c_str() ) << "deg" << endl;
}
else if (DataBuffer=="R=") {
check_R = true;
ConfigFile >> DataBuffer ;
cout << "R: " << atof( DataBuffer.c_str() ) << "mm" << endl;
}
//Position method
else if (DataBuffer=="X=") {
check_X = true;
ConfigFile >> DataBuffer ;
cout << "X: " << atof( DataBuffer.c_str() ) << "mm" << endl;
}
else if (DataBuffer=="Y=") {
check_Y = true;
ConfigFile >> DataBuffer ;
cout << "Y: " << atof( DataBuffer.c_str() ) << "mm"<< endl;
}
else if (DataBuffer=="Z=") {
check_Z = true;
ConfigFile >> DataBuffer ;
cout << "Z: " << atof( DataBuffer.c_str() ) << "mm" << endl;
}
//General
else if (DataBuffer=="Shape=") {
check_Shape = true;
ConfigFile >> DataBuffer ;
cout << "Shape: " << DataBuffer << endl;
}
// Cylindrical shape
else if (DataBuffer== "Radius=") {
check_Radius = true;
ConfigFile >> DataBuffer ;
cout << "Hira Radius: " << atof( DataBuffer.c_str() ) << "mm" << endl;
}
// Squared shape
else if (DataBuffer=="Width=") {
check_Width = true;
ConfigFile >> DataBuffer ;
cout << "Hira Width: " <<atof( DataBuffer.c_str() ) << "mm" << endl;
}
else if (DataBuffer== "Height=") {
check_Height = true;
ConfigFile >> DataBuffer ;
cout << "Hira Height: " << atof( DataBuffer.c_str() ) << "mm" << endl;
}
// Common
else if (DataBuffer=="Thickness=") {
check_Thickness = true;
ConfigFile >> DataBuffer ;
cout << "Hira Thickness: " << atof( DataBuffer.c_str() ) << "mm" << endl;
}
else if (DataBuffer== "Scintillator=") {
check_Scintillator = true ;
ConfigFile >> DataBuffer ;
cout << "Hira Scintillator type: " << DataBuffer << endl;
}
else if (DataBuffer=="LeadThickness=") {
check_LeadThickness = true;
ConfigFile >> DataBuffer ;
cout << "Lead Thickness : " << atof( DataBuffer.c_str() ) << "mm" << endl;
}
///////////////////////////////////////////////////
// If no Detector Token and no comment, toggle out
else
{ReadingStatus = false; cout << "Wrong Token Sequence: Getting out " << DataBuffer << endl ;}
/////////////////////////////////////////////////
// If All necessary information there, toggle out
if ( check_Theta && check_Phi && check_R && check_Thickness && check_Radius && check_LeadThickness && check_Scintillator && check_Height && check_Width && check_Shape && check_X && check_Y && check_Z )
{
NumberOfDetector++;
// Reinitialisation of Check Boolean
check_Theta = false ;
check_Phi = false ;
check_R = false ;
check_Thickness = false ;
check_Radius = false ;
check_LeadThickness = false ;
check_Scintillator = false ;
check_Height = false ;
check_Width = false ;
check_Shape = false ;
check_X = false ;
check_Y = false ;
check_Z = false ;
ReadingStatus = false ;
cout << "///"<< endl ;
}
}
}
}
///////////////////////////////////////////////////////////////////////////
void THiraPhysics::AddParameterToCalibrationManager()
{
CalibrationManager* Cal = CalibrationManager::getInstance();
for(int i = 0 ; i < NumberOfDetector ; i++)
{
for( int j = 0 ; j < 16 ; j++)
{
Cal->AddParameter("Hira", "Detector"+itoa(i+1)+"_E","Hira_Detector"+itoa(i+1)+"_E") ;
Cal->AddParameter("Hira", "Detector"+itoa(i+1)+"_T","Hira_Detector"+itoa(i+1)+"_T") ;
}
}
}
///////////////////////////////////////////////////////////////////////////
void THiraPhysics::InitializeRootInputRaw()
{
TChain* inputChain = RootInput::getInstance()->GetChain() ;
inputChain->SetBranchStatus ( "Hira" , true ) ;
inputChain->SetBranchStatus ( "fHira_*" , true ) ;
inputChain->SetBranchAddress( "Hira" , &EventData ) ;
}
///////////////////////////////////////////////////////////////////////////
void THiraPhysics::InitializeRootInputPhysics()
{
/*TChain* inputChain = RootInput::getInstance()->GetChain();
inputChain->SetBranchStatus ( "Hira", true );
inputChain->SetBranchStatus ( "DetectorNumber", true );
inputChain->SetBranchStatus ( "Energy", true );
inputChain->SetBranchStatus ( "Time", true );
inputChain->SetBranchAddress( "Hira", &EventPhysics );*/
}
///////////////////////////////////////////////////////////////////////////
void THiraPhysics::InitializeRootOutput()
{
TTree* outputTree = RootOutput::getInstance()->GetTree() ;
outputTree->Branch( "Hira" , "THiraPhysics" , &EventPhysics ) ;
}
///////////////////////////////////////////////////////////////////////////
void THiraPhysics::BuildPhysicalEvent()
{
BuildSimplePhysicalEvent() ;
}
///////////////////////////////////////////////////////////////////////////
void THiraPhysics::BuildSimplePhysicalEvent()
{
}
////////////////////////////////////////////////////////////////////////////////
// Construct Method to be pass to the DetectorFactory //
////////////////////////////////////////////////////////////////////////////////
NPL::VDetector* THiraPhysics::Construct(){
return (NPL::VDetector*) new THiraPhysics();
}
////////////////////////////////////////////////////////////////////////////////
// Registering the construct method to the factory //
////////////////////////////////////////////////////////////////////////////////
extern "C"{
class proxy{
public:
proxy(){
NPL::DetectorFactory::getInstance()->AddToken("HIRAArray","Hira");
NPL::DetectorFactory::getInstance()->AddDetector("HIRAArray",THiraPhysics::Construct);
}
};
proxy p;
}
NPLib/HiRA/THiraPhysics.h 0 → 100644
+ 100
− 0
View file @ 977bda7b
#ifndef __HiraPhysics__
#define __HiraPhysics__
/*****************************************************************************
* Copyright (C) 2009-2014 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: Adrien MATTA contact address: a.matta@surrey.ac.uk *
* *
* Creation Date : November 2009 *
* Last update : *
*---------------------------------------------------------------------------*
* Decription: *
* This class hold theHira Detector Physics *
* *
*---------------------------------------------------------------------------*
* Comment: *
* *
* *
*****************************************************************************/
// STL
#include <vector>
using namespace std ;
// ROOT
#include "TObject.h"
// NPL
#include "THiraData.h"
#include "NPVDetector.h"
//#include "../include/CalibrationManager.h"
class THiraPhysics : public TObject, public NPL::VDetector
{
public: // Constructor and Destructor
THiraPhysics();
~THiraPhysics();
public:
void Clear();
void Clear(const Option_t*) {};
public: // Calibrated Data
vector<UShort_t> DetectorNumber ;
vector<Double_t> Energy;
vector<Double_t> Time;
public: // inherrited from VDetector
// Read stream at ConfigFile to pick-up parameters of detector (Position,...) using Token
void ReadConfiguration(string);
// Add Parameter to the CalibrationManger
void AddParameterToCalibrationManager();
// Activated associated Branches and link it to the private member DetectorData address
// In this method mother Branches (Detector) AND daughter leaf (fDetector_parameter) have to be activated
void InitializeRootInputRaw() ;
// Activated associated Branches and link it to the private member DetectorPhysics address
// In this method mother Branches (Detector) AND daughter leaf (parameter) have to be activated
void InitializeRootInputPhysics() ;
// Create associated branches and associated private member DetectorPhysics address
void InitializeRootOutput();
// This method is called at each event read from the Input Tree. Aime is to build treat Raw dat in order to extract physical parameter.
void BuildPhysicalEvent();
// Same as above, but only the simplest event and/or simple method are used (low multiplicity, faster algorythm but less efficient ...).
// This method aimed to be used for analysis performed during experiment, when speed is requiered.
// NB: This method can eventually be the same as BuildPhysicalEvent.
void BuildSimplePhysicalEvent();
// Same as above but for online analysis
void BuildOnlinePhysicalEvent() {BuildPhysicalEvent();};
// Give and external THiraData object to THiraPhysics. Needed for online analysis for example.
void SetRawDataPointer(THiraData* rawDataPointer) {EventData = rawDataPointer;}
// Those two method all to clear the Event Physics or Data
void ClearEventPhysics() {Clear();}
void ClearEventData() {EventData->Clear();}
private: // Data not writted in the tree
int NumberOfDetector ;//!
THiraData* EventData ;//!
THiraPhysics* EventPhysics ;//!
public:// Static constructor to be passed to the Detector Factory
static NPL::VDetector* Construct();
ClassDef(THiraPhysics,1) // HiraPhysics structure
};
#endif
NPSimulation/Hira/Hira.cc 0 → 100644
+ 729
− 0
View file @ 977bda7b
This diff is collapsed.
NPSimulation/Hira/Hira.hh 0 → 100644
+ 235
− 0
View file @ 977bda7b
#ifndef Sharc_h
#define Sharc_h 1
/*****************************************************************************
* Copyright (C) 2009-2013 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: Pierre MORFOUACE contact address: *
* *
* Creation Date : Feb 2015 *
* Last update : *
*---------------------------------------------------------------------------*
* Decription: *
* This class describe the Sharc Silicon detector *
* *
*---------------------------------------------------------------------------*
* Comment: *
* *
*****************************************************************************/
// C++ header
#include <string>
#include <vector>
// G4 header defining G4 types
#include "globals.hh"
// G4 headers
#include "G4ThreeVector.hh"
#include "G4RotationMatrix.hh"
#include "G4LogicalVolume.hh"
#include "G4VisAttributes.hh"
#include "G4MultiFunctionalDetector.hh"
// NPSimulation header
#include "NPSVDetector.hh"
// NPLib
#include "THiraData.h"
using namespace std;
using namespace CLHEP;
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
namespace HIRA
{
// Resolution
const G4double ResoTime = 0.212765957 ;// = 500ps // Unit is ns/2.35
const G4double ResoCsI = 0.08 ;// = 188 kev of resolution // Unit is MeV/2.35
const G4double ResoThickSi = 0.022 ;// = 52keV of Resolution // Unit is MeV/2.35
const G4double ResoThinSi = 0.064 ;// = 150keV of Resolution // Unit is MeV/2.35
const G4double EnergyThreshold = 100*keV;
const G4double FaceFront = 7.*cm ;
const G4double FaceBack = 7*cm ;
const G4double Length = 7.*cm ;
const G4int NumberOfStrip = 32;
const G4double SiliconFace = 64*mm ;
const G4double AluStripThickness = 0.4*micrometer ;
const G4double ThinSiThickness = 65*micrometer ;
const G4double ThickSiThickness = 1500*micrometer ;
const G4double DistBetweenSi = 1*cm;
const G4double VacBoxThickness = 1*cm ;
const G4double MylarCsIThickness = 3*micrometer;
const G4double CsIThickness = 4.*cm + 2*MylarCsIThickness ;
const G4double CsIFaceFront = 3.2*cm;
const G4double CsIFaceBack = 3.7*cm;
const G4double DistInterCsI = 0.2*mm;
const G4double ClusterFaceFront = 7*cm;
const G4double ClusterFaceBack = 9*cm;
// Starting at the front and going to CsI
const G4double AluStripFrontThinSi_PosZ = Length* -0.5 + 0.5*AluStripThickness;
const G4double ThinSi_PosZ = AluStripFrontThinSi_PosZ + 0.5*AluStripThickness + 0.5*ThinSiThickness;
const G4double AluStripBackThinSi_PosZ = ThinSi_PosZ + 0.5*ThinSiThickness + 0.5*AluStripThickness;
const G4double AluStripFrontThickSi_PosZ = AluStripBackThinSi_PosZ + DistBetweenSi + 0.5*AluStripThickness + 0.5*AluStripThickness;
const G4double ThickSi_PosZ = AluStripFrontThickSi_PosZ + 0.5*AluStripThickness + 0.5*ThickSiThickness;
const G4double AluStripBackThickSi_PosZ = ThickSi_PosZ + 0.5*ThickSiThickness + 0.5*AluStripThickness;
const G4double VacBox_PosZ = AluStripBackThickSi_PosZ + 0.5*AluStripThickness + 0.5* VacBoxThickness;
const G4double CsI_PosZ = VacBox_PosZ + 0.5*VacBoxThickness + 0.5*CsIThickness;
}
using namespace HIRA ;
class Hira : public NPS::VDetector
{
////////////////////////////////////////////////////
/////// Default Constructor and Destructor /////////
////////////////////////////////////////////////////
public:
Hira() ;
~Hira() ;
////////////////////////////////////////////////////
//////// Specific Function of this Class ///////////
////////////////////////////////////////////////////
public:
// By Position Method
void AddTelescope( G4ThreeVector Pos1 ,
G4ThreeVector Pos2 ,
G4ThreeVector Pos3 ,
G4ThreeVector Pos4 );
// Effectively construct Volume
// Avoid to have two time same code for Angle and Point definition
void VolumeMaker( G4int DetectorNumber ,
G4ThreeVector MMpos ,
G4RotationMatrix* MMrot ,
G4LogicalVolume* world );
private:
G4LogicalVolume* m_LogicThinSi;
G4LogicalVolume* m_LogicThickSi;
G4LogicalVolume* m_LogicCsICrystal;
G4LogicalVolume* m_LogicCluster;
G4LogicalVolume* m_logicMotherVolume;
////////////////////////////////////////////////////
//////// Specific Function of this Class ///////////
////////////////////////////////////////////////////
public:
// To add a box detector
// Effectively construct Volume
void ConstructCsICrystal();
void ConstructThinSi();
void ConstructThickSi();
////////////////////////////////////////////////////
///////// Inherite from VDetector 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() ;
// Read sensitive part and fill the Root tree.
// Called at in the EventAction::EndOfEventAvtion
void ReadSensitive(const G4Event* event) ;
////////////////////////////////////////////////////
///////////Event class to store Data////////////////
////////////////////////////////////////////////////
private:
THiraData* m_EventHira ;
////////////////////////////////////////////////////
///////////////// Scorer Related ///////////////////
////////////////////////////////////////////////////
private:
// Initialize all Scorer
void InitializeScorers() ;
// Silicon Associate Scorer
G4MultiFunctionalDetector* m_ThinSiStripScorer ;
G4MultiFunctionalDetector* m_ThickSiStripScorer ;
// CsI Associate Scorer
G4MultiFunctionalDetector* m_CsIScorer ;
private:
// Initialize material used in detector definition
void InitializeMaterial();
// List of material
// Si
G4Material* m_MaterialSilicon;
// Al
G4Material* m_MaterialAluminium;
// CsI
G4Material* m_MaterialCsI;
// Vacuum
G4Material* m_MaterialVacuum ;
// Mylar
G4Material* m_MaterialMylar;
// PCB
G4Material* m_MaterialPCB;
////////////////////////////////////////////////////
///////////////Private intern Data//////////////////
////////////////////////////////////////////////////
private:
// True if the detector is a Cluster, false if a single Phoswich
vector<bool> m_Type;
// Detector position
vector<G4ThreeVector> m_Pos;
vector<G4RotationMatrix*> m_Rot;
// Used for "By Point Definition"
//vector<G4ThreeVector> m_TL ; // Top Left Corner Position Vector
//vector<G4ThreeVector> m_BL ; // Bottom Left Corner Position Vector
//vector<G4ThreeVector> m_BR ; // Bottom Right Corner Position Vector
//vector<G4ThreeVector> m_TR ; // 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 ; // |
private:/// Visualisation Attribute:
G4VisAttributes* m_SiliconVisAtt ;
G4VisAttributes* m_CsIVisAtt;
public:
static NPS::VDetector* Construct();
};
#endif
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment