Skip to content
Snippets Groups Projects
Commit 68603db1 authored by LEMAIR Louis's avatar LEMAIR Louis
Browse files

* Adding files to Nebula Project

  - Modified detector geometry
  - 8Hp2p source changed for a simple neutron source
  - Physics and source geometry updated
parent 067234ee
No related branches found
No related tags found
No related merge requests found
...@@ -40,7 +40,7 @@ ...@@ -40,7 +40,7 @@
// NPTool header // NPTool header
#include "Nebula.hh" #include "Nebula.hh"
#include "CalorimeterScorers.hh" #include "PlasticBar.hh"
#include "InteractionScorers.hh" #include "InteractionScorers.hh"
#include "RootOutput.h" #include "RootOutput.h"
#include "MaterialManager.hh" #include "MaterialManager.hh"
...@@ -216,7 +216,7 @@ void Nebula::ReadSensitive(const G4Event* ){ ...@@ -216,7 +216,7 @@ void Nebula::ReadSensitive(const G4Event* ){
/////////// ///////////
// Module scorer // Module scorer
CalorimeterScorers::PS_Calorimeter* Scorer= (CalorimeterScorers::PS_Calorimeter*) m_ModuleScorer->GetPrimitive(0); PlasticBar::PS_PlasticBar* Scorer= (PlasticBar::PS_PlasticBar*) m_ModuleScorer->GetPrimitive(0);
unsigned int size = Scorer->GetMult(); unsigned int size = Scorer->GetMult();
for(unsigned int i = 0 ; i < size ; i++){ for(unsigned int i = 0 ; i < size ; i++){
...@@ -245,18 +245,18 @@ void Nebula::InitializeScorers() { ...@@ -245,18 +245,18 @@ void Nebula::InitializeScorers() {
// Otherwise the scorer is initialise // Otherwise the scorer is initialise
// Module // Module
vector<int> level; level.push_back(0); vector<int> level; level.push_back(0);
G4VPrimitiveScorer* ModuleCalorimeter= new CalorimeterScorers::PS_Calorimeter("ModuleCalorimeter",level, 0) ; G4VPrimitiveScorer* ModulePlasticBar= new PlasticBar::PS_PlasticBar("ModulePlasticBar",level, 0) ;
G4VPrimitiveScorer* ModuleInteraction= new InteractionScorers::PS_Interactions("ModuleInteraction",ms_InterCoord, 0) ; G4VPrimitiveScorer* ModuleInteraction= new InteractionScorers::PS_Interactions("ModuleInteraction",ms_InterCoord, 0) ;
//and register it to the multifunctionnal detector //and register it to the multifunctionnal detector
m_ModuleScorer->RegisterPrimitive(ModuleCalorimeter); m_ModuleScorer->RegisterPrimitive(ModulePlasticBar);
m_ModuleScorer->RegisterPrimitive(ModuleInteraction); m_ModuleScorer->RegisterPrimitive(ModuleInteraction);
G4SDManager::GetSDMpointer()->AddNewDetector(m_ModuleScorer) ; G4SDManager::GetSDMpointer()->AddNewDetector(m_ModuleScorer) ;
// Veto // Veto
G4VPrimitiveScorer* VetoCalorimeter= new CalorimeterScorers::PS_Calorimeter("VetoCalorimeter",level, 0) ; G4VPrimitiveScorer* VetoPlasticBar= new PlasticBar::PS_PlasticBar("VetoPlasticBar",level, 0) ;
G4VPrimitiveScorer* VetoInteraction= new InteractionScorers::PS_Interactions("VetoInteraction",ms_InterCoord, 0) ; G4VPrimitiveScorer* VetoInteraction= new InteractionScorers::PS_Interactions("VetoInteraction",ms_InterCoord, 0) ;
//and register it to the multifunctionnal detector //and register it to the multifunctionnal detector
m_VetoScorer->RegisterPrimitive(VetoCalorimeter); m_VetoScorer->RegisterPrimitive(VetoPlasticBar);
m_VetoScorer->RegisterPrimitive(VetoInteraction); m_VetoScorer->RegisterPrimitive(VetoInteraction);
G4SDManager::GetSDMpointer()->AddNewDetector(m_VetoScorer) ; G4SDManager::GetSDMpointer()->AddNewDetector(m_VetoScorer) ;
......
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Beam
Particle= 8He
Energy= 2000 MeV
SigmaEnergy= 100 MeV
SigmaThetaX= 1 deg
SigmaPhiY= 1 deg
SigmaX= 5 mm
SigmaY= 5 mm
MeanThetaX= 0 deg
MeanPhiY= 0 deg
MeanX= 0 mm
MeanY= 0 mm
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
QFSReaction
Beam= 8He
Target= 1H
Scattered= 1H
KnockedOut= 1H
Heavy= 7H
ExcitationEnergyBeam= 0.0 MeV
ExcitationEnergyHeavy= 3.6800 MeV
MomentumSigma= 50.0
ShootHeavy= 1
ShootLight= 1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Target Target
THICKNESS= 10 micrometer THICKNESS= 0.49 cm
RADIUS= 20 mm RADIUS= 28 mm
MATERIAL= CD2 MATERIAL= Vacuum
ANGLE= 0 deg ANGLE= 0 deg
X= 0 mm X= 3.5 cm
Y= 0 mm Y= 0 mm
Z= -4000 mm Z= 0 mm
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
NEBULA NEBULA
POS= 0 0 5 m POS= 0 0 12 m
NumberOfModule= 30 NumberOfModule= 15
Veto= 1 Veto= 1
Frame= 1 Frame= 1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
NEBULA NEBULA
POS= 0 0 5.15 m POS= 0 0 12.15 m
NumberOfModule= 30 NumberOfModule= 15
Veto= 0 Veto= 0
Frame= 0 Frame= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
NEBULA NEBULA
POS= 0 0 10 m POS= 0 0 16 m
NumberOfModule= 30 NumberOfModule= 30
Veto= 1 Veto= 1
Frame= 1 Frame= 1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
NEBULA NEBULA
POS= 0 0 10.15 m POS= 0 0 16.15 m
NumberOfModule= 30 NumberOfModule= 30
Veto= 0 Veto= 0
Frame= 0 Frame= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
NEBULA NEBULA
POS= 0 0 3 m POS= 0 0 18 m
NumberOfModule= 15 NumberOfModule= 30
Veto= 1 Veto= 1
Frame= 1 Frame= 1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
NEBULA NEBULA
POS= 0 0 3.15 m POS= 0 0 18.15 m
NumberOfModule= 15 NumberOfModule= 30
Veto= 0 Veto= 0
Frame= 0 Frame= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Samurai %Samurai
% POS= 0 0 0 mm % POS= 0 0 0 mm
...@@ -78,4 +78,4 @@ NEBULA ...@@ -78,4 +78,4 @@ NEBULA
% %InvertY= 1 % %InvertY= 1
% %InvertD= 1 % %InvertD= 1
%
EmPhysicsList Option4
DefaultCutOff 1
IonBinaryCascadePhysics 0
NPIonInelasticPhysics 0
EmExtraPhysics 0
HadronElasticPhysics 0
StoppingPhysics 0
OpticalPhysics 0
HadronPhysicsINCLXX 0
HadronPhysicsQGSP_BIC_HP 1
Decay 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Beam
Particle= n
Energy= 196 MeV
SigmaEnergy= 0.5 MeV
SigmaThetaX= 0.01 deg
SigmaPhiY= 0.01 deg
SigmaX= 0.5 mm
SigmaY= 0.5 mm
MeanThetaX= 0 deg
MeanPhiY= 0 deg
MeanX= 3.5 cm
MeanY= 0 mm
ZEmission= -1000 mm
0 1
90 1
180 1
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment