Skip to content
Snippets Groups Projects
Commit 40c1d483 authored by Unknown's avatar Unknown
Browse files

No commit message 

No commit message
parent 3b09f025
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 1347 additions and 6 deletions
Inputs/DetectorConfiguration/plastic.detector 0 → 100644
+ 21
− 0
View file @ 40c1d483
%Fichier de configuration manip E225
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ScintillatorPlastic
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Plastic
THETA= 0
PHI= 0
R= 150
Thickness= 20
Radius= 50
LeadThickness= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Plastic
THETA= 0
PHI= 0
R= 175
Thickness= 20
Radius= 50
LeadThickness= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
NPLib/Makefile
+ 10
− 1
View file @ 40c1d483
......@@ -269,7 +269,7 @@ INCLUDE = -I$(CLHEP_BASE_DIR)/include
#------------------------------------------------------------------------------
SHARELIB = Vdetec InputOutputRoot InitCond InterCoord \
Must2All GaspardData AnnularS1Data \
Must2All GaspardData AnnularS1Data PlasticData\
Reaction EnergyLoss
all: $(SHARELIB)
......@@ -311,6 +311,15 @@ ifeq ($(ARCH),macosx)
cd lib; ln -sf libGaspardPyysics.so libGaspardPhysics.dylib
endif
## Gaspard ##
PlasticData:
make -C ./Plastic
cp ./Plastic/*.so ./lib ; cp ./Plastic/*.h ./include
ifeq ($(ARCH),macosx)
cd lib; ln -sf libPlasticData.so libPlasticData.dylib
cd lib; ln -sf libPlasticPyysics.so libPlasticPhysics.dylib
endif
############# Simulation ##############
......
NPLib/Plastic/Makefile 0 → 100644
+ 302
− 0
View file @ 40c1d483
# Makefile for the ROOT test programs.
# This Makefile shows nicely how to compile and link applications
# using the ROOT libraries on all supported platforms.
#
# Copyright (c) 2000 Rene Brun and Fons Rademakers
#
# Author: Fons Rademakers, 29/2/2000
ROOTCONFIG := root-config
ARCH := $(shell $(ROOTCONFIG) --arch)
PLATFORM := $(shell $(ROOTCONFIG) --platform)
ALTCC := $(shell $(ROOTCONFIG) --cc)
ALTCXX := $(shell $(ROOTCONFIG) --cxx)
ALTF77 := $(shell $(ROOTCONFIG) --f77)
ALTLD := $(shell $(ROOTCONFIG) --ld)
#CXX =
ObjSuf = o
SrcSuf = cxx
ExeSuf =
DllSuf = so
OutPutOpt = -o # keep whitespace after "-o"
ifeq (debug,$(findstring debug,$(ROOTBUILD)))
OPT = -g
OPT2 = -g
else
ifneq ($(findstring debug, $(strip $(shell $(ROOTCONFIG) --config))),)
OPT = -g
OPT2 = -g
else
OPT = -O
OPT2 = -O2
endif
endif
ROOTCFLAGS := $(shell $(ROOTCONFIG) --cflags)
ROOTLDFLAGS := $(shell $(ROOTCONFIG) --ldflags)
ROOTLIBS := $(shell $(ROOTCONFIG) --libs)
ROOTGLIBS := $(shell $(ROOTCONFIG) --glibs)
HASTHREAD := $(shell $(ROOTCONFIG) --has-thread)
ROOTDICTTYPE := $(shell $(ROOTCONFIG) --dicttype)
NOSTUBS := $(shell $(ROOTCONFIG) --nostubs)
ROOTCINT := rootcint
ifeq ($(ARCH),linux)
# Linux with egcs, gcc 2.9x, gcc 3.x
CXX = g++
CXXFLAGS = $(OPT2) -Wall -fPIC
LD = g++
LDFLAGS = $(OPT2)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxkcc)
# Linux with the KAI compiler
CXX = KCC --one_instantiation_per_object
CXXFLAGS = $(OPT) -fPIC +K0
LD = KCC
LDFLAGS = $(OPT) $(shell $(ROOTCONFIG) --cflags)
SOFLAGS =
endif
ifeq ($(ARCH),linuxicc)
# Linux with Intel icc compiler
ICC_MAJOR := $(shell icc -v 2>&1 | awk '{ if (NR==1) print $$2 }' | \
cut -d'.' -f1)
ICC_MINOR := $(shell icc -v 2>&1 | awk '{ if (NR==1) print $$2 }' | \
cut -d'.' -f2)
CXX = icc
CXXFLAGS = $(OPT) -fPIC -wd1476
LD = icpc
LDFLAGS = $(OPT)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxppcgcc)
# PPC Linux with gcc and glibc
CXX = g++
CXXFLAGS = $(OPT2) -Wall -fPIC
LD = g++
LDFLAGS = $(OPT2)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxia64gcc)
# Itanium Linux with gcc 2.9x
CXX = g++
CXXFLAGS = $(OPT2) -Wall -fPIC
LD = g++
LDFLAGS = $(OPT2)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxia64sgi)
# Itanium Linux with sgiCC
CXX = sgiCC
CXXFLAGS = $(OPT) -Wall -fPIC
LD = gsgiCC
LDFLAGS = $(OPT)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxia64ecc)
# Itanium Linux with Intel icc (was ecc)
ICC_MAJOR := $(shell icc -v 2>&1 | awk '{ if (NR==1) print $$2 }' | \
cut -d'.' -f1)
ICC_MINOR := $(shell icc -v 2>&1 | awk '{ if (NR==1) print $$2 }' | \
cut -d'.' -f2)
CXX = icc
CXXFLAGS = $(OPT) -fPIC -wd1476 -ftz
LD = icpc
LDFLAGS = $(OPT)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxx8664gcc)
# AMD Opteron and Intel EM64T (64 bit mode) Linux with gcc 3.x
CXX = g++
CXXFLAGS = $(OPT2) -Wall -fPIC
LD = g++
LDFLAGS = $(OPT2)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxppc64gcc)
# PPC64 Linux with gcc 3.x
CXX = g++
CXXFLAGS = $(OPT) -Wall -fPIC
LD = g++
LDFLAGS = $(OPT)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxx8664icc)
# AMD Opteron and Intel EM64T (64 bit mode) Linux with Intel icc compiler
CXX = icc
CXXFLAGS = $(OPT) -fPIC -wd1476 -wd1572
LD = icpc
LDFLAGS = $(OPT)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxalphagcc)
# Alpha Linux with gcc
CXX = g++
CXXFLAGS = $(OPT2) -Wall -fPIC
LD = g++
LDFLAGS = $(OPT2)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxmips)
# GNU/Linux on mips (BE/LE, O32/N32/N64) with gcc
CXX = g++
CXXFLAGS = $(OPT2) -Wall -fPIC
LD = g++
LDFLAGS = $(OPT2)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxhppa)
# GNU/Linux on hppa with gcc
CXX = g++
CXXFLAGS = $(OPT2) -Wall -fPIC
LD = g++
LDFLAGS = $(OPT2)
SOFLAGS = -shared
endif
ifeq ($(ARCH),linuxarm)
# ARM Linux with egcs
CXX = g++
CXXFLAGS = $(OPT) -Wall -fPIC
LD = g++
LDFLAGS = $(OPT)
SOFLAGS = -shared
endif
ifeq ($(ARCH),macosx)
# MacOS X with cc (GNU cc 2.95.2 and gcc 3.3)
MACOSX_MINOR := $(shell sw_vers | sed -n 's/ProductVersion://p' | cut -d . -f 2)
MACOSXTARGET := MACOSX_DEPLOYMENT_TARGET=10.$(MACOSX_MINOR)
CXX = g++
CXXFLAGS = $(OPT2) -pipe -Wall -W -Woverloaded-virtual
LD = $(MACOSXTARGET) g++
LDFLAGS = $(OPT2) -bind_at_load
# The SOFLAGS will be used to create the .dylib,
# the .so will be created separately
ifeq ($(subst $(MACOSX_MINOR),,1234),1234)
DllSuf = so
else
DllSuf = dylib
endif
UNDEFOPT = dynamic_lookup
ifneq ($(subst $(MACOSX_MINOR),,12),12)
UNDEFOPT = suppress
LD = g++
endif
SOFLAGS = -dynamiclib -single_module -undefined $(UNDEFOPT)
endif
ifeq ($(ARCH),macosxicc)
# MacOS X with Intel icc compiler
MACOSX_MINOR := $(shell sw_vers | sed -n 's/ProductVersion://p' | cut -d . -f 2)
MACOSXTARGET := MACOSX_DEPLOYMENT_TARGET=10.$(MACOSX_MINOR)
ifeq ($(MACOSX_MINOR),5)
MACOSX_MINOR := 4
endif
CXX = icc
CXXFLAGS = $(OPT) -fPIC -wd1476
LD = $(MACOSXTARGET) icpc
LDFLAGS = $(OPT)
# The SOFLAGS will be used to create the .dylib,
# the .so will be created separately
ifeq ($(subst $(MACOSX_MINOR),,1234),1234)
DllSuf = so
else
DllSuf = dylib
endif
SOFLAGS = -dynamiclib -single_module -undefined dynamic_lookup
endif
ifeq ($(ARCH),macosx64)
# MacOS X >= 10.4 with gcc 64 bit mode (GNU gcc 4.*)
# Only specific option (-m64) comes from root-config
MACOSX_MINOR := $(shell sw_vers | sed -n 's/ProductVersion://p' | cut -d . -f 2)
MACOSXTARGET := MACOSX_DEPLOYMENT_TARGET=10.$(MACOSX_MINOR)
CXX = g++
CXXFLAGS = $(OPT2) -pipe -Wall -W -Woverloaded-virtual
LD = $(MACOSXTARGET) g++ -m64
LDFLAGS = $(OPT2) -bind_at_load
# The SOFLAGS will be used to create the .dylib,
# the .so will be created separately
ifeq ($(subst $(MACOSX_MINOR),,1234),1234)
DllSuf = so
else
DllSuf = dylib
endif
SOFLAGS = -m64 -dynamiclib -single_module -undefined dynamic_lookup
endif
ifeq ($(ARCH),macosxxlc)
# MacOS X with IBM xlC compiler
MACOSX_MINOR := $(shell sw_vers | sed -n 's/ProductVersion://p' | cut -d . -f 2)
MACOSXTARGET := MACOSX_DEPLOYMENT_TARGET=10.$(MACOSX_MINOR)
CXX = xlC
CXXFLAGS = $(OPT)
LD = $(MACOSXTARGET) xlC
LDFLAGS = $(OPT) -Wl,-bind_at_load
# The SOFLAGS will be used to create the .dylib,
# the .so will be created separately
DllSuf = dylib
UNDEFOPT = dynamic_lookup
ifneq ($(subst $(MACOSX_MINOR),,12),12)
UNDEFOPT = suppress
LD = xlC
endif
SOFLAGS = -qmkshrobj -single_module -undefined $(UNDEFOPT)
endif
CXXFLAGS += $(ROOTCFLAGS)
LDFLAGS += $(ROOTLDFLAGS)
LIBS = $(ROOTLIBS) $(SYSLIBS)
GLIBS = $(ROOTGLIBS) $(SYSLIBS)
INCLUDE = -I$(CLHEP_BASE_DIR)/include
#------------------------------------------------------------------------------
SHARELIB = libPlasticData.so
all: $(SHARELIB)
#------------------------------------------------------------------------------
############### Detector ##############
## Plastic ##
libPlasticData.so: TPlasticData.o TPlasticDataDict.o
$(LD) $(SOFLAGS) $^ $(OutPutOpt) $@
TPlasticDataDict.cxx: TPlasticData.h
rootcint -f $@ -c $^
# dependances
TPlasticData.o: TPlasticData.cxx TPlasticData.h
#######################################
############# Clean and More ##########
clean:
@rm -f core *~ *.o *Dict*
distclean:
make clean; rm -f *.so
.SUFFIXES: .$(SrcSuf)
###
.$(SrcSuf).$(ObjSuf):
$(CXX) $(CXXFLAGS) $(INCLUDE) -c $<
NPLib/Plastic/TPlasticData.cxx 0 → 100644
+ 39
− 0
View file @ 40c1d483
#include <iostream>
using namespace std;
#include "TPlasticData.h"
ClassImp(TPlasticData)
TPlasticData::TPlasticData()
{
}
TPlasticData::~TPlasticData()
{
}
void TPlasticData::Clear()
{
fPlastic_Energy.clear();
fPlastic_Number.clear();
}
void TPlasticData::Dump()
{
cout << "XXXXXXXXXXXXXXXXXXXXXXXX New Event XXXXXXXXXXXXXXXXX" << endl;
for(short i = 0 ; i<fPlastic_Energy.size() ; i ++)
{
cout << "Plastic Number " << fPlastic_Number[i] << " Energy: " << fPlastic_Energy[i] << " Time: "<< fPlastic_Time[i] << endl;
}
}
NPLib/Plastic/TPlasticData.h 0 → 100644
+ 41
− 0
View file @ 40c1d483
#ifndef __PLASTICDATA__
#define __PLASTICDATA__
#include <vector>
#include "TObject.h"
class TPlasticData : public TObject {
private:
// ADC
vector<double> fPlastic_Energy ;
vector<double> fPlastic_Time ;
vector<short> fPlastic_Number ;
public:
TPlasticData();
virtual ~TPlasticData();
void Clear();
void Dump();
///////////////////// GETTERS ////////////////////////
// (E)
double GetEnergy(int i) {return fPlastic_Energy[i];}
// (T)
double GetTime(int i) {return fPlastic_Time[i];}
// (N)
int GetPlasticNumber(int i) {return fPlastic_Number[i];}
///////////////////// SETTERS ////////////////////////
// (E)
void SetEnergy(double E) {fPlastic_Energy.push_back(E);}
void SetTime(double T) {fPlastic_Time.push_back(T);}
void SetPlasticNumber(int N) {fPlastic_Number.push_back(N);}
//
ClassDef(TPlasticData,1) // PlasticData structure
};
#endif
NPSimulation/GNUmakefile
+ 1
− 1
View file @ 40c1d483
......@@ -19,7 +19,7 @@ G4WORKDIR = ./
CPPFLAGS += $(shell root-config --cflags)
CPPFLAGS += -I$(NPLIB)/include
EXTRALIBS = $(shell root-config --glibs)
EXTRALIBS += -L$(NPLIB)/lib -lMust2Data -lAnnularS1Data -lGaspardData \
EXTRALIBS += -L$(NPLIB)/lib -lMust2Data -lAnnularS1Data -lGaspardData -lPlasticData\
-lInitialConditions -lInteractionCoordinates \
-lReaction -lIORoot
.PHONY: all
......
NPSimulation/include/GeneralScorers.hh
+ 1
− 1
View file @ 40c1d483
......@@ -22,7 +22,7 @@
* deposit, time of flight or position,... particle by particle for each *
* event.Because standard scorer provide by G4 don't work this way but using*
* a global ID for each event you should not use those scorer with some G4 *
* provided ones or being very carefull doing so. * *
* provided ones or being very carefull doing so. *s
*****************************************************************************/
#include "G4VPrimitiveScorer.hh"
#include "G4THitsMap.hh"
......
NPSimulation/include/MUST2Array.hh
+ 0
− 1
View file @ 40c1d483
......@@ -35,7 +35,6 @@ namespace MUST2
{
// Resolution
const G4double ResoTimeMust = 0.212765957 ;// = 500ps // Unit is ns/2.35
const G4double ResoTimePPAC = 0.106382979 ;// = 250ps // Unit is ns/2.35
const G4double ResoSiLi = 0.055 ;// = 130 keV of resolution // Unit is MeV/2.35
const G4double ResoCsI = 0.08 ;// = 188 kev of resolution // Unit is MeV/2.35
const G4double ResoStrip = 0.022 ;// = 52keV of Resolution // Unit is MeV/2.35
......
NPSimulation/include/Plastic.hh 0 → 100644
+ 120
− 0
View file @ 40c1d483
#ifndef Plastic_h
#define Plastic_h 1
/*****************************************************************************
* Copyright (C) 2009 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 : september 2009 *
* Last update : *
*---------------------------------------------------------------------------*
* Decription: *
* This class describe a modular plastic array *
* *
*---------------------------------------------------------------------------*
* 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"
// NPTool header
#include "VDetector.hh"
#include "TPlasticData.h"
using namespace std;
class Plastic : public VDetector
{
////////////////////////////////////////////////////
/////// Default Constructor and Destructor /////////
////////////////////////////////////////////////////
public:
Plastic() ;
virtual ~Plastic() ;
////////////////////////////////////////////////////
//////// Specific Function of this Class ///////////
////////////////////////////////////////////////////
public:
// By Angle Method
void AddPlastic( G4double R ,
G4double Theta ,
G4double Phi ,
G4double PlasticThickness ,
G4double PlasticRadius ,
G4double LeadThickness );
void VolumeMaker(G4ThreeVector Det_pos, int DetNumber,G4LogicalVolume* world) ;
////////////////////////////////////////////////////
///////// 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) ;
public: // Material
void InitializeMaterial() ;
// Platic
G4Material* m_MaterialPlastic ;
// Lead
G4Material* m_MaterialLead ;
public: // Scorer
// Initialize all Scorer used by the MUST2Array
void InitializeScorers() ;
// Silicon Associate Scorer
G4MultiFunctionalDetector* m_PlasticScorer ;
////////////////////////////////////////////////////
///////////Event class to store Data////////////////
////////////////////////////////////////////////////
private:
TPlasticData* m_Event ;
////////////////////////////////////////////////////
///////////////Private intern Data//////////////////
////////////////////////////////////////////////////
private:
// if true a Lead plate is added in front or back of the detector
vector<double> m_LeadThickness ;
vector<double> m_PlasticThickness ;
vector<double> m_PlasticRadius ;
// Used for "By Angle Definition"
vector<G4double> m_R ; // |
vector<G4double> m_Theta ; // > Spherical coordinate plastic volume center
vector<G4double> m_Phi ; // |
};
#endif
NPSimulation/include/PlasticScorers.hh 0 → 100644
+ 94
− 0
View file @ 40c1d483
#ifndef GeneralScorers_h
#define GeneralScorers_h 1
/*****************************************************************************
* Copyright (C) 2009 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 : September 2009 *
* Last update : *
*---------------------------------------------------------------------------*
* Decription: *
* This class old some of the Platic Scorer, shared by different detector. *
*---------------------------------------------------------------------------*
* Comment: *
*****************************************************************************/
#include "G4VPrimitiveScorer.hh"
#include "G4THitsMap.hh"
namespace PLASTIC {
class PSEnergy : public G4VPrimitiveScorer
{
public: // with description
PSEnergy(G4String name, G4int depth = 0);
virtual ~PSEnergy();
protected: // with description
virtual G4bool ProcessHits(G4Step*, G4TouchableHistory*);
public:
virtual void Initialize(G4HCofThisEvent*);
virtual void EndOfEvent(G4HCofThisEvent*);
virtual void clear();
virtual void DrawAll();
virtual void PrintAll();
private:
G4int HCID;
G4THitsMap<G4double>* EvtMap;
};
class PSTOF : public G4VPrimitiveScorer
{
public: // with description
PSTOF(G4String name, G4int depth = 0);
virtual ~PSTOF();
protected: // with description
virtual G4bool ProcessHits(G4Step*, G4TouchableHistory*);
public:
virtual void Initialize(G4HCofThisEvent*);
virtual void EndOfEvent(G4HCofThisEvent*);
virtual void clear();
virtual void DrawAll();
virtual void PrintAll();
private:
G4int HCID;
G4THitsMap<G4double>* EvtMap;
};
class PSDetectorNumber : public G4VPrimitiveScorer
{
public: // with description
PSDetectorNumber(G4String name, G4int depth = 0 );
virtual ~PSDetectorNumber();
protected: // with description
virtual G4bool ProcessHits(G4Step*, G4TouchableHistory*);
public:
virtual void Initialize(G4HCofThisEvent*);
virtual void EndOfEvent(G4HCofThisEvent*);
virtual void clear();
virtual void DrawAll();
virtual void PrintAll();
private:
G4int HCID;
G4THitsMap<G4int>* EvtMap;
};
}
#endif
NPSimulation/src/GeneralScorers.cc
+ 2
− 1
View file @ 40c1d483
......@@ -20,7 +20,7 @@
* deposit, time of flight or position,... particle by particle for each *
* event.Because standard scorer provide by G4 don't work this way but using*
* a global ID for each event you should not use those scorer with some G4 *
* provided ones or being very carefull doing so. * *
* provided ones or being very carefull doing so. *
*****************************************************************************/
#include "GeneralScorers.hh"
#include "G4UnitsTable.hh"
......@@ -438,6 +438,7 @@ G4bool PSDetectorNumber::ProcessHits(G4Step* aStep, G4TouchableHistory*)
if (edep < 100*keV) return FALSE;
G4int index = aStep->GetTrack()->GetTrackID();
EvtMap->set(index, DetNbr);
return TRUE;
}
......
NPSimulation/src/MUST2Array.cc
+ 0
− 1
View file @ 40c1d483
......@@ -1066,7 +1066,6 @@ void MUST2Array::ReadSensitive(const G4Event* event)
G4double T = *(Time_itr->second) ;
if (TTrackID == ETrackID) {
T = RandGauss::shoot(T, ResoTimePPAC) ;
m_Event->SetMMStripXTTime(RandGauss::shoot(T, ResoTimeMust)) ;
m_Event->SetMMStripYTTime(RandGauss::shoot(T, ResoTimeMust)) ;
}
......
NPSimulation/src/Plastic.cc 0 → 100644
+ 476
− 0
View file @ 40c1d483
/*****************************************************************************
* Copyright (C) 2009 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 : January 2009 *
* Last update : *
*---------------------------------------------------------------------------*
* Decription: *
* This class describe a 20um Silicium detector *
* *
*---------------------------------------------------------------------------*
* Comment: *
* *
*****************************************************************************/
// C++ headers
#include <sstream>
#include <cmath>
//G4 Geometry object
#include "G4Tubs.hh"
//G4 sensitive
#include "G4SDManager.hh"
#include "G4MultiFunctionalDetector.hh"
//G4 various object
#include "G4MaterialTable.hh"
#include "G4Element.hh"
#include "G4ElementTable.hh"
#include "G4Transform3D.hh"
#include "G4PVPlacement.hh"
#include "G4VisAttributes.hh"
#include "G4Colour.hh"
// NPTool header
#include "Plastic.hh"
#include "PlasticScorers.hh"
#include "RootOutput.h"
// CLHEP header
#include "CLHEP/Random/RandGauss.h"
using namespace std;
using namespace CLHEP;
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
namespace PLASTIC
{
// Energy and time Resolution
const G4double ResoTime = 4.2 ;// = 10ns of Resolution // Unit is MeV/2.35
const G4double ResoEnergy = 0.42 ;// = 1MeV of Resolution // Unit is MeV/2.35
}
using namespace PLASTIC ;
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
// Plastic Specific Method
Plastic::Plastic()
{
InitializeMaterial();
m_Event = new TPlasticData() ;
}
Plastic::~Plastic()
{
delete m_MaterialPlastic ;
delete m_MaterialLead ;
delete m_PlasticScorer ;
}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
void Plastic::AddPlastic( G4double R ,
G4double Theta ,
G4double Phi ,
G4double PlasticThickness ,
G4double PlasticRadius ,
G4double LeadThickness )
{
m_R.push_back(R) ;
m_Theta.push_back(Theta) ;
m_Phi.push_back(Phi) ;
m_PlasticThickness.push_back(PlasticThickness) ;
m_PlasticRadius.push_back(PlasticRadius) ;
m_LeadThickness.push_back(LeadThickness);
}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
// Virtual Method of VDetector class
// Read stream at Configfile to pick-up parameters of detector (Position,...)
// Called in DetecorConstruction::ReadDetextorConfiguration Method
void Plastic::ReadConfiguration(string Path)
{
ifstream ConfigFile ;
ConfigFile.open(Path.c_str()) ;
string LineBuffer ;
string DataBuffer ;
G4double Theta = 0 , Phi = 0 , R = 0 , Thickness = 0 , Radius = 0 , LeadThickness = 0;
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 ReadingStatus = false ;
while (!ConfigFile.eof())
{
getline(ConfigFile, LineBuffer);
// If line is a Start Up MUST2 bloc, Reading toggle to true
if (LineBuffer.compare(0, 7, "Plastic") == 0)
{
G4cout << "///" << G4endl ;
G4cout << "Platic found: " << G4endl ;
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) {/*do nothing */;}
// Finding another telescope (safety), toggle out
else if (DataBuffer.compare(0, 6, "Plastic") == 0) {
cout << "WARNING: Another Telescope is find before standard sequence of Token, Error may occured in Telecope definition" << endl ;
ReadingStatus = false ;
}
//Angle method
else if (DataBuffer.compare(0, 6, "THETA=") == 0) {
check_Theta = true;
ConfigFile >> DataBuffer ;
Theta = atof(DataBuffer.c_str()) ;
Theta = Theta * deg;
cout << "Theta: " << Theta / deg << endl;
}
else if (DataBuffer.compare(0, 4, "PHI=") == 0) {
check_Phi = true;
ConfigFile >> DataBuffer ;
Phi = atof(DataBuffer.c_str()) ;
Phi = Phi * deg;
cout << "Phi: " << Phi / deg << endl;
}
else if (DataBuffer.compare(0, 2, "R=") == 0) {
check_R = true;
ConfigFile >> DataBuffer ;
R = atof(DataBuffer.c_str()) ;
R = R * mm;
cout << "R: " << R/mm << endl;
}
else if (DataBuffer.compare(0, 7, "Radius=") == 0) {
check_Radius = true;
ConfigFile >> DataBuffer ;
Radius = atof(DataBuffer.c_str()) ;
Radius = Radius * mm;
cout << "Plastic Radius: " << Radius/mm << endl;
}
else if (DataBuffer.compare(0, 10, "Thickness=") == 0) {
check_Thickness = true;
ConfigFile >> DataBuffer ;
Thickness = atof(DataBuffer.c_str()) ;
Thickness = Thickness * mm;
cout << "Plastic Thickness: " << Thickness/mm << endl;
}
else if (DataBuffer.compare(0, 14, "LeadThickness=") == 0) {
check_LeadThickness = true;
ConfigFile >> DataBuffer ;
LeadThickness = atof(DataBuffer.c_str()) ;
LeadThickness = LeadThickness * mm;
cout << "Lead Thickness : " << LeadThickness/mm << endl;
}
///////////////////////////////////////////////////
// If no Detector Token and no comment, toggle out
else
{ReadingStatus = false; G4cout << "Wrong Token Sequence: Getting out " << DataBuffer << G4endl ;}
/////////////////////////////////////////////////
// If All necessary information there, toggle out
if ( check_Theta && check_Phi && check_R && check_Thickness && check_Radius && check_LeadThickness)
{
AddPlastic( R ,
Theta ,
Phi ,
Thickness ,
Radius ,
LeadThickness);
// Reinitialisation of Check Boolean
check_Theta = false ;
check_Phi = false ;
check_R = false ;
check_Thickness = false ;
check_Radius = false ;
check_LeadThickness = false ;
ReadingStatus = false ;
cout << "///"<< endl ;
}
}
}
}
// Construct detector and inialise sensitive part.
// Called After DetecorConstruction::AddDetector Method
void Plastic::ConstructDetector(G4LogicalVolume* world)
{
G4ThreeVector Det_pos = G4ThreeVector(0, 0, 0) ;
for (ushort i = 0 ; i < m_R.size() ; i++)
{
G4double wX = m_R[i] * sin(m_Theta[i] ) * cos(m_Phi[i] ) ;
G4double wY = m_R[i] * sin(m_Theta[i] ) * sin(m_Phi[i] ) ;
G4double wZ = m_R[i] * cos(m_Theta[i] ) ;
Det_pos = G4ThreeVector(wX, wY, wZ) ;
G4LogicalVolume* logicPlastic = NULL ;
VolumeMaker(Det_pos , i+1, world) ;
}
}
void Plastic::VolumeMaker(G4ThreeVector Det_pos, int DetNumber, G4LogicalVolume* world)
{
////////////////////////////////////////////////////////////////
////////////// Starting Volume Definition //////////////////////
////////////////////////////////////////////////////////////////
G4PVPlacement* PVPBuffer ;
// Name of the module
std::ostringstream DetectorNumber ;
DetectorNumber << DetNumber ;
G4String Name = "Plastic" + DetectorNumber.str() ;
int i = DetNumber-1;
// Definition of the volume containing the sensitive detector
if(m_PlasticThickness[i]>0 && m_PlasticRadius[i]>0)
{
G4Tubs* solidPlastic = new G4Tubs( Name ,
0 ,
m_PlasticRadius[i] ,
m_PlasticThickness[i]/2 ,
0*deg ,
360*deg );
G4LogicalVolume* logicPlastic = new G4LogicalVolume(solidPlastic, m_MaterialPlastic, Name, 0, 0, 0);
logicPlastic->SetSensitiveDetector(m_PlasticScorer);
G4VisAttributes* PlastVisAtt = new G4VisAttributes(G4Colour(0.0, 0.0, 0.9)) ;
logicPlastic->SetVisAttributes(PlastVisAtt) ;
PVPBuffer = new G4PVPlacement( 0 ,
Det_pos ,
logicPlastic ,
Name ,
world ,
false ,
0 );
}
if(m_LeadThickness[i]>0&& m_PlasticRadius[i]>0)
{
G4Tubs* solidLead = new G4Tubs( Name ,
0 ,
m_PlasticRadius[i] ,
m_LeadThickness[i]/2 ,
0*deg ,
360*deg );
G4LogicalVolume* logicLead = new G4LogicalVolume(solidLead, m_MaterialLead, Name, 0, 0, 0);
G4VisAttributes* LeadVisAtt = new G4VisAttributes(G4Colour(0.1, 0.1, 0.1)) ;
logicLead->SetVisAttributes(LeadVisAtt) ;
PVPBuffer = new G4PVPlacement( 0 ,
Det_pos+(m_PlasticThickness[i]/2+m_LeadThickness[i]/2)*Det_pos.unit() ,
logicLead ,
Name ,
world ,
false ,
0 );
}
}
// Add Detector branch to the EventTree.
// Called After DetecorConstruction::AddDetector Method
void Plastic::InitializeRootOutput()
{
RootOutput *pAnalysis = RootOutput::getInstance();
TTree *pTree = pAnalysis->GetTree();
pTree->Branch("Plastic", "TPlasticData", &m_Event) ;
}
// Read sensitive part and fill the Root tree.
// Called at in the EventAction::EndOfEventAvtion
void Plastic::ReadSensitive(const G4Event* event)
{
G4String DetectorNumber ;
m_Event->Clear() ;
//////////////////////////////////////////////////////////////////////////////////////
//////////////////////// Used to Read Event Map of detector //////////////////////////
//////////////////////////////////////////////////////////////////////////////////////
std::map<G4int, G4int*>::iterator DetectorNumber_itr ;
std::map<G4int, G4double*>::iterator Energy_itr ;
std::map<G4int, G4double*>::iterator Time_itr ;
G4THitsMap<G4int>* DetectorNumberHitMap ;
G4THitsMap<G4double>* EnergyHitMap ;
G4THitsMap<G4double>* TimeHitMap ;
//////////////////////////////////////////////////////////////////////////////////////
//////////////////////////////////////////////////////////////////////////////////////
// Read the Scorer associate to the Silicon Strip
//Detector Number
G4int StripDetCollectionID = G4SDManager::GetSDMpointer()->GetCollectionID("PlasticScorer/PlasticNumber") ;
DetectorNumberHitMap = (G4THitsMap<G4int>*)(event->GetHCofThisEvent()->GetHC(StripDetCollectionID)) ;
DetectorNumber_itr = DetectorNumberHitMap->GetMap()->begin() ;
//Energy
G4int StripEnergyCollectionID = G4SDManager::GetSDMpointer()->GetCollectionID("PlasticScorer/Energy") ;
EnergyHitMap = (G4THitsMap<G4double>*)(event->GetHCofThisEvent()->GetHC(StripEnergyCollectionID)) ;
Energy_itr = EnergyHitMap->GetMap()->begin() ;
//Time of Flight
G4int StripTimeCollectionID = G4SDManager::GetSDMpointer()->GetCollectionID("PlasticScorer/Time") ;
TimeHitMap = (G4THitsMap<G4double>*)(event->GetHCofThisEvent()->GetHC(StripTimeCollectionID)) ;
Time_itr = TimeHitMap->GetMap()->begin() ;
G4int sizeN = DetectorNumberHitMap->entries() ;
G4int sizeE = EnergyHitMap->entries() ;
G4int sizeT = TimeHitMap->entries() ;
// Loop on Plastic Number
for (G4int l = 0 ; l < sizeN ; l++) {
G4int N = *(DetectorNumber_itr->second) ;
G4int NTrackID = DetectorNumber_itr->first - N ;
if (N > 0) {
m_Event->SetPlasticNumber(N) ;
// Energy
Energy_itr = EnergyHitMap->GetMap()->begin();
for (G4int h = 0 ; h < sizeE ; h++) {
G4int ETrackID = Energy_itr->first - N ;
G4double E = *(Energy_itr->second) ;
if (ETrackID == NTrackID) {
m_Event->SetEnergy(RandGauss::shoot(E, ResoEnergy)) ;
}
Energy_itr++;
}
// Time
Time_itr = TimeHitMap->GetMap()->begin();
for (G4int h = 0 ; h < sizeT ; h++) {
G4int TTrackID = Time_itr->first - N ;
G4double T = *(Time_itr->second) ;
if (TTrackID == NTrackID) {
m_Event->SetTime(RandGauss::shoot(T, ResoTime)) ;
}
Time_itr++;
}
}
DetectorNumber_itr++;
}
// clear map for next event
DetectorNumberHitMap ->clear();
EnergyHitMap ->clear() ;
}
////////////////////////////////////////////////////////////////
void Plastic::InitializeMaterial()
{
////////////////////////////////////////////////////////////////
/////////////////Element Definition ///////////////////////////
////////////////////////////////////////////////////////////////
G4String symbol ;
G4double density = 0. , a = 0, z = 0 ;
G4int ncomponents = 0, natoms = 0 ;
// for Plastic
G4Element* H = new G4Element("Hydrogen" , symbol = "H" , z = 1 , a = 1.01 * g / mole);
G4Element* C = new G4Element("Carbon" , symbol = "C" , z = 6 , a = 12.011 * g / mole);
////////////////////////////////////////////////////////////////
/////////////////Material Definition ///////////////////////////
////////////////////////////////////////////////////////////////
// Pb
a = 207.2 * g / mole;
density = 11.34 * g / cm3;
m_MaterialLead = new G4Material("Lead", z = 82 , a, density);
// Plastic
density = 1.243 * g / cm3;
m_MaterialPlastic = new G4Material("Plastic", density, ncomponents = 2);
m_MaterialPlastic->AddElement(H , natoms = 10);
m_MaterialPlastic->AddElement(C , natoms = 14);
}
////////////////////////////////////////////////////////////////
void Plastic::InitializeScorers()
{
m_PlasticScorer = new G4MultiFunctionalDetector("PlasticScorer") ;
G4SDManager::GetSDMpointer()->AddNewDetector(m_PlasticScorer);
G4VPrimitiveScorer* DetNbr = new PSDetectorNumber("PlasticNumber", 0) ;
G4VPrimitiveScorer* Energy = new PSEnergy("Energy", 0) ;
G4VPrimitiveScorer* Time = new PSTOF("Time", 0) ;
//and register it to the multifunctionnal detector
m_PlasticScorer->RegisterPrimitive(DetNbr) ;
m_PlasticScorer->RegisterPrimitive(Energy) ;
m_PlasticScorer->RegisterPrimitive(Time) ;
}
////////////////////////////////////////////////////////////////
NPSimulation/src/PlasticScorer.cc 0 → 100644
+ 240
− 0
View file @ 40c1d483
/*****************************************************************************
* Copyright (C) 2009 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 : September 2009 *
* Last update : *
*---------------------------------------------------------------------------*
* Decription: *
* This class old some of the Platic Scorer, shared by different detector. *
*---------------------------------------------------------------------------*
* Comment: *
*****************************************************************************/
#include "PlasticScorers.hh"
#include "G4UnitsTable.hh"
using namespace PLASTIC;
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
// Energy Scorer (deal with multiple particle hit)
PSEnergy::PSEnergy(G4String name, G4int depth)
: G4VPrimitiveScorer(name, depth), HCID(-1)
{
}
PSEnergy::~PSEnergy()
{
}
G4bool PSEnergy::ProcessHits(G4Step* aStep, G4TouchableHistory*)
{
std::string name = aStep->GetTrack()->GetVolume()->GetName();
std::string nbr ;
size_t found;
found=name.find("Plastic");
found = found + 7;
int numberOfCharacterInDetectorNumber = name.length() - (int)found ;
for( int i = found ; i < found + numberOfCharacterInDetectorNumber ; i++ )
nbr += name[i] ;
G4int DetNbr = atoi( nbr.c_str() ) ;
G4double edep = aStep->GetTotalEnergyDeposit();
if (edep < 100*keV) return FALSE;
G4int index = aStep->GetTrack()->GetTrackID();
EvtMap->add(DetNbr+index, edep);
return TRUE;
}
void PSEnergy::Initialize(G4HCofThisEvent* HCE)
{
EvtMap = new G4THitsMap<G4double>(GetMultiFunctionalDetector()->GetName(), GetName());
if (HCID < 0) {
HCID = GetCollectionID(0);
}
HCE->AddHitsCollection(HCID, (G4VHitsCollection*)EvtMap);
}
void PSEnergy::EndOfEvent(G4HCofThisEvent*)
{
;
}
void PSEnergy::clear()
{
EvtMap->clear();
}
void PSEnergy::DrawAll()
{
;
}
void PSEnergy::PrintAll()
{
G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl;
G4cout << " PrimitiveScorer " << GetName() << G4endl;
G4cout << " Number of entries " << EvtMap->entries() << G4endl;
std::map<G4int, G4double*>::iterator itr = EvtMap->GetMap()->begin();
for (; itr != EvtMap->GetMap()->end(); itr++) {
G4cout << " copy no.: " << itr->first
<< " energy deposit: " << G4BestUnit(*(itr->second), "Energy")
<< G4endl;
}
}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
//Detector Number Scorer
PSDetectorNumber::PSDetectorNumber(G4String name, G4int depth)
: G4VPrimitiveScorer(name, depth), HCID(-1)
{
}
PSDetectorNumber::~PSDetectorNumber()
{
;
}
G4bool PSDetectorNumber::ProcessHits(G4Step* aStep, G4TouchableHistory*)
{
std::string name = aStep->GetTrack()->GetVolume()->GetName();
std::string nbr ;
size_t found;
found=name.find("Plastic");
found = found + 7;
int numberOfCharacterInDetectorNumber = name.length() - (int)found ;
for( int i = found ; i < found + numberOfCharacterInDetectorNumber ; i++ )
nbr += name[i] ;
G4int DetNbr = atoi( nbr.c_str() ) ;
G4double edep = aStep->GetTotalEnergyDeposit();
if (edep < 100*keV) return FALSE;
G4int index = aStep->GetTrack()->GetTrackID();
EvtMap->set(DetNbr + index, DetNbr);
return TRUE;
}
void PSDetectorNumber::Initialize(G4HCofThisEvent* HCE)
{
EvtMap = new G4THitsMap<G4int>(GetMultiFunctionalDetector()->GetName(), GetName());
if (HCID < 0) {
HCID = GetCollectionID(0);
}
HCE->AddHitsCollection(HCID, (G4VHitsCollection*)EvtMap);
}
void PSDetectorNumber::EndOfEvent(G4HCofThisEvent*)
{
;
}
void PSDetectorNumber::clear()
{
EvtMap->clear();
}
void PSDetectorNumber::DrawAll()
{
;
}
void PSDetectorNumber::PrintAll()
{
G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl;
G4cout << " PrimitiveScorer " << GetName() << G4endl;
G4cout << " Number of entries " << EvtMap->entries() << G4endl;
std::map<G4int, G4int*>::iterator itr = EvtMap->GetMap()->begin();
for (; itr != EvtMap->GetMap()->end(); itr++) {
G4cout << " copy no.: " << itr->first
<< " energy deposit: " << G4BestUnit(*(itr->second), "Energy")
<< G4endl;
}
}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
//TOF Scorer//
PSTOF::PSTOF(G4String name, G4int depth)
: G4VPrimitiveScorer(name, depth), HCID(-1)
{
;
}
PSTOF::~PSTOF()
{
;
}
G4bool PSTOF::ProcessHits(G4Step* aStep, G4TouchableHistory*)
{
std::string name = aStep->GetTrack()->GetVolume()->GetName();
std::string nbr ;
size_t found;
found=name.find("Plastic");
found = found + 7;
int numberOfCharacterInDetectorNumber = name.length() - (int)found ;
for( int i = found ; i < found + numberOfCharacterInDetectorNumber ; i++ )
nbr += name[i] ;
G4int DetNbr = atoi( nbr.c_str() ) ;
G4double TOF = aStep->GetPreStepPoint()->GetGlobalTime();
G4double edep = aStep->GetTotalEnergyDeposit();
if (edep < 100*keV) return FALSE;
G4int index = aStep->GetTrack()->GetTrackID();
EvtMap->set(DetNbr + index, TOF);
return TRUE;
}
void PSTOF::Initialize(G4HCofThisEvent* HCE)
{
EvtMap = new G4THitsMap<G4double>(GetMultiFunctionalDetector()->GetName(),
GetName());
if (HCID < 0) {
HCID = GetCollectionID(0);
}
HCE->AddHitsCollection(HCID, (G4VHitsCollection*)EvtMap);
}
void PSTOF::EndOfEvent(G4HCofThisEvent*)
{
;
}
void PSTOF::clear()
{
EvtMap->clear();
}
void PSTOF::DrawAll()
{
;
}
void PSTOF::PrintAll()
{
G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl ;
G4cout << " PrimitiveScorer " << GetName() << G4endl ;
G4cout << " Number of entries " << EvtMap->entries() << G4endl ;
}
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