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

Updtaing Vendeta Inelastic project

parent 69a158dd
No related branches found
No related tags found
1 merge request!27Draft: [Epic] Preparation of the environement for the new GaseousDetectorScorers...
Pipeline #337383 passed
Stage: build-NPLib
Stage: build-NPSimulation
Stage: test
Hide whitespace changes
Inline Side-by-side
Projects/Inelastic/RunToTreat.txt
+ 1
− 2
View file @ e290603b
TTreeName
SimulatedTree
RootFileName
../../Outputs/Simulation/vendeta_el_*.root
../../Outputs/Simulation/vendeta_1st_inel_*.root
./root/simulation/neutron_56Fe_gs.root
Projects/Inelastic/Vendeta_inelastic.detector
+ 16
− 391
View file @ e290603b
......@@ -2,458 +2,83 @@
Target
THICKNESS= 10 micrometer
RADIUS= 20 mm
MATERIAL= 238U
MATERIAL= Fe
ANGLE= 0 deg
X= 0 mm
Y= 0 mm
Z= 0 mm
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 20 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 30 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 40 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 50 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 60 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 70 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 80 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 90 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 100 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 110 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 120 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 130 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 140 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 150 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
R= 2000 mm
THETA= 160 deg
PHI= 0 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 20 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 30 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 40 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 50 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 60 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 70 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 80 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 90 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 100 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 110 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 120 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 130 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 140 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 150 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 160 deg
PHI= 20 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 20 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 30 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 40 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 50 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 60 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 70 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 80 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 90 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 100 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 110 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 120 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 130 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 140 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 150 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 160 deg
PHI= 40 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 20 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 30 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 40 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 50 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 60 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 70 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 80 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 90 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 100 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 110 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 120 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 130 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 140 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 150 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 160 deg
PHI= 140 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 20 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 30 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 40 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 50 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 60 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 70 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 80 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 90 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 100 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 110 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 120 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 130 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 140 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 150 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 160 deg
PHI= 160 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 20 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 30 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 40 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 50 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 60 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 70 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 80 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 90 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 100 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 110 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 120 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 130 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 140 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 150 deg
PHI= 180 deg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Vendeta
R= 2500 mm
THETA= 160 deg
PHI= 180 deg
Projects/Inelastic/macro/CheckKinematics.C 0 → 100644
+ 41
− 0
View file @ e290603b
TRandom3 ra;
void CheckKinematics(){
TFile* ofile = new TFile("kine.root","recreate");
NPL::Reaction* r1 = new NPL::Reaction("7Li(p,1n)7Be@13.95");
double beam_energy_init = r1->GetBeamEnergy();
double eloss = 0.070;
double beam_energy;
double ThetaCM;
double E3, Theta3, E4, Theta4;
TTree* tree = new TTree("KineTree","KineTree");
tree->Branch("beam_energy",&beam_energy,"beam_energy/D");
tree->Branch("ThetaCM",&ThetaCM,"ThetaCM/D");
tree->Branch("E3",&E3,"E3/D");
tree->Branch("Theta3",&Theta3,"Theta3/D");
tree->Branch("E4",&E4,"E4/D");
tree->Branch("Theta4",&Theta4,"Theta4/D");
for(int i=0; i<1e6; i++){
beam_energy = ra.Uniform(beam_energy_init - eloss, beam_energy_init);
ThetaCM = ra.Uniform(0,180);
r1->SetThetaCM(ThetaCM*deg);
r1->SetBeamEnergy(beam_energy);
r1->KineRelativistic(Theta3,E3,Theta4,E4);
Theta3 = Theta3/deg;
Theta4 = Theta4/deg;
tree->Fill();
}
tree->Write();
ofile->Close();
}
Projects/Inelastic/macro/DirectKinematics.C 0 → 100644
+ 52
− 0
View file @ e290603b
TRandom3 ra;
double density_7Li = 0.534; // g/cm3
///////////////////////////////////////////////////
void DirectKinematics(double thickness_mg_cm2 = 0.3){
NPL::EnergyLoss* p7Li_loss;
p7Li_loss = new NPL::EnergyLoss("../EnergyLossTable/proton_7Li.G4table","G4Table",100);
TFile* ofile = new TFile("direct_kine.root","recreate");
NPL::Reaction* r1 = new NPL::Reaction("p(7Li,1n)7Be@2.8");
double thickness_micrometer;
thickness_micrometer = thickness_mg_cm2/(density_7Li*1e3)*1e4;
cout << thickness_mg_cm2 << " mg/cm2 -> " << thickness_micrometer << " um" << endl;
double beam_energy_init = r1->GetBeamEnergy();
double eloss;
eloss = p7Li_loss->EnergyLossCalculation(beam_energy_init,thickness_micrometer*micrometer,0);
cout << "Energy Loss = " << eloss << " MeV" << endl;
double beam_energy;
double ThetaCM;
double E3, Theta3, E4, Theta4;
TTree* tree = new TTree("KineTree","KineTree");
tree->Branch("beam_energy",&beam_energy,"beam_energy/D");
tree->Branch("ThetaCM",&ThetaCM,"ThetaCM/D");
tree->Branch("E3",&E3,"E3/D");
tree->Branch("Theta3",&Theta3,"Theta3/D");
tree->Branch("E4",&E4,"E4/D");
tree->Branch("Theta4",&Theta4,"Theta4/D");
for(int i=0; i<1e6; i++){
beam_energy = ra.Uniform(beam_energy_init - eloss, beam_energy_init);
ThetaCM = ra.Uniform(0,180);
r1->SetThetaCM(ThetaCM*deg);
r1->SetBeamEnergy(beam_energy);
r1->KineRelativistic(Theta3,E3,Theta4,E4);
Theta3 = Theta3/deg;
Theta4 = Theta4/deg;
tree->Fill();
}
tree->Write();
ofile->Close();
}
Projects/Inelastic/neutron_inelastic.reaction 0 → 100644
+ 28
− 0
View file @ e290603b
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%% Reaction file for 11Li(d,3He)10He reaction %%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Beam
Particle= 1n
Energy= 1.1 MeV
SigmaEnergy= 0.016 MeV
SigmaThetaX= 0. deg
SigmaPhiY= 0. deg
SigmaX= 0 mm
SigmaY= 0 mm
MeanThetaX= 0 deg
MeanPhiY= 0 deg
MeanX= 0 mm
MeanY= 0 mm
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
TwoBodyReaction
Beam= 1n
Target= 56Fe
Light= 1n
Heavy= 56Fe
ExcitationEnergy3= 0.0 MeV
ExcitationEnergy4= 0.0 MeV
CrossSectionPath= flat.txt CS10He
ShootLight= 1
ShootHeavy= 1
Projects/Inelastic/project.config 0 → 100644
+ 5
− 0
View file @ e290603b
Project e850
AnalysisOutput= ./root/analysis
SimulationOutput= ./root/simulation
EnergyLoss= ./EnergyLossTable
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