From fcd74c8147abf1cc8feb0a03d87d93f67890fd63 Mon Sep 17 00:00:00 2001
From: Morfouace <pierre.morfouace@gmail.com>
Date: Mon, 31 Jul 2023 13:47:59 +0200
Subject: [PATCH] Adding project for neutron inelastic
---
Projects/Inelastic/Analysis.cxx | 156 ++++++
Projects/Inelastic/Analysis.h | 67 +++
Projects/Inelastic/CMakeLists.txt | 5 +
.../Inelastic/Geometry/GenerateGeometry.C | 47 ++
Projects/Inelastic/PhysicsListOption.txt | 15 +
Projects/Inelastic/RunToTreat.txt | 4 +
Projects/Inelastic/Vendeta_inelastic.detector | 459 ++++++++++++++++++
Projects/Inelastic/run.mac | 1 +
Projects/Inelastic/sim.sh | 7 +
9 files changed, 761 insertions(+)
create mode 100644 Projects/Inelastic/Analysis.cxx
create mode 100644 Projects/Inelastic/Analysis.h
create mode 100644 Projects/Inelastic/CMakeLists.txt
create mode 100644 Projects/Inelastic/Geometry/GenerateGeometry.C
create mode 100644 Projects/Inelastic/PhysicsListOption.txt
create mode 100644 Projects/Inelastic/RunToTreat.txt
create mode 100644 Projects/Inelastic/Vendeta_inelastic.detector
create mode 100644 Projects/Inelastic/run.mac
create mode 100755 Projects/Inelastic/sim.sh
diff --git a/Projects/Inelastic/Analysis.cxx b/Projects/Inelastic/Analysis.cxx
new file mode 100644
index 000000000..452e686e3
--- /dev/null
+++ b/Projects/Inelastic/Analysis.cxx
@@ -0,0 +1,156 @@
+/*****************************************************************************
+ * 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 : march 2025 *
+ * Last update : *
+ *---------------------------------------------------------------------------*
+ * Decription: *
+ * Class describing the property of an Analysis object *
+ * *
+ *---------------------------------------------------------------------------*
+ * Comment: *
+ * *
+ * *
+ *****************************************************************************/
+#include<iostream>
+using namespace std;
+#include"Analysis.h"
+#include"NPAnalysisFactory.h"
+#include"NPDetectorManager.h"
+#include"NPOptionManager.h"
+#include"RootOutput.h"
+#include"RootInput.h"
+////////////////////////////////////////////////////////////////////////////////
+Analysis::Analysis(){
+}
+////////////////////////////////////////////////////////////////////////////////
+Analysis::~Analysis(){
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::Init(){
+ m_Vendeta= (TVendetaPhysics*) m_DetectorManager->GetDetector("Vendeta");
+ InitialConditions = new TInitialConditions();
+ InteractionCoordinates = new TInteractionCoordinates();
+ ReactionConditions = new TReactionConditions();
+
+ InitInputBranch();
+ InitOutputBranch();
+
+ my_Reaction = new NPL::Reaction("1n(238U,1n)238U@2.1");
+
+ neutron = new NPL::Nucleus("1n");
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::TreatEvent(){
+ ReInitValue();
+ Einit = InitialConditions->GetKineticEnergy(0);
+ double init_ThetaLab = ReactionConditions->GetTheta(0)*deg;
+ double init_BeamEnergy = ReactionConditions->GetBeamEnergy();
+ //my_Reaction->SetBeamEnergy(init_BeamEnergy);
+
+ neutron->SetKineticEnergy(Einit);
+ double beam_TOF = neutron->GetTimeOfFlight();
+
+ double Xtarget = InitialConditions->GetIncidentPositionX();
+ double Ytarget = InitialConditions->GetIncidentPositionY();
+ double Ztarget = 0;//InitialConditions->GetIncidentPositionZ();
+ TVector3 TargetPos = TVector3(Xtarget,Ytarget,Ztarget);
+
+ m_Vendeta->SetAnodeNumber(0);
+ m_Vendeta->BuildPhysicalEvent();
+ unsigned int size = m_Vendeta->LG_DetectorNumber.size();
+ for(int i=0; i<size; i++){
+ double Rdet, R;
+ Rdet = m_Vendeta->GetDistanceFromTarget(m_Vendeta->LG_DetectorNumber[i]);
+ TVector3 DetPos = m_Vendeta->GetVectorDetectorPosition(m_Vendeta->LG_DetectorNumber[i]);
+ //TVector3 HitPos = DetPos-TargetPos;
+
+ R= Rdet*mm;
+ Distance.push_back(R);
+ Det.push_back(m_Vendeta->LG_DetectorNumber[i]);
+ T.push_back(m_Vendeta->LG_Time[i]);
+ double T_stop = (m_Vendeta->LG_Time[i])*1e-9;
+ neutron->SetTimeOfFlight((T_stop-beam_TOF)/(Rdet*1e-3));
+ //neutron->SetTimeOfFlight((T_stop)/(Rdet*1e-3-8e-3));
+ Elab.push_back(neutron->GetEnergy());
+
+
+ double DeltaTheta = atan(63.5/Rdet);
+ double exp_ThetaLab = m_Vendeta->GetVectorDetectorPosition(m_Vendeta->LG_DetectorNumber[i]).Theta();
+ double random_ThetaLab = ra.Uniform(exp_ThetaLab-DeltaTheta, exp_ThetaLab+DeltaTheta);
+ double dEx = my_Reaction->ReconstructRelativistic(Elab[i], random_ThetaLab);
+
+ ThetaLab.push_back(random_ThetaLab/deg);
+ //ThetaLab.push_back(exp_ThetaLab/deg);
+ Ex.push_back(dEx);
+ }
+}
+
+///////////////////////////////////////////////////////////////////////////////
+void Analysis::InitOutputBranch() {
+ RootOutput::getInstance()->GetTree()->Branch("Einit",&Einit,"Einit/D");
+ RootOutput::getInstance()->GetTree()->Branch("ThetaLab",&ThetaLab);
+ RootOutput::getInstance()->GetTree()->Branch("Elab",&Elab);
+ RootOutput::getInstance()->GetTree()->Branch("Ex",&Ex);
+ RootOutput::getInstance()->GetTree()->Branch("T",&T);
+ RootOutput::getInstance()->GetTree()->Branch("Distance",&Distance);
+ RootOutput::getInstance()->GetTree()->Branch("Det",&Det);
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::InitInputBranch(){
+ RootInput:: getInstance()->GetChain()->SetBranchStatus("InitialConditions",true );
+ RootInput:: getInstance()->GetChain()->SetBranchStatus("fIC_*",true );
+ RootInput:: getInstance()->GetChain()->SetBranchAddress("InitialConditions",&InitialConditions);
+
+ RootInput:: getInstance()->GetChain()->SetBranchStatus("ReactionConditions",true );
+ RootInput:: getInstance()->GetChain()->SetBranchStatus("fRC_*",true );
+ RootInput:: getInstance()->GetChain()->SetBranchAddress("ReactionConditions",&ReactionConditions);
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::ReInitValue(){
+ Einit = -100;
+ Ex.clear();
+ ThetaLab.clear();
+ Elab.clear();
+ T.clear();
+ Distance.clear();
+ Det.clear();
+}
+
+////////////////////////////////////////////////////////////////////////////////
+void Analysis::End(){
+}
+
+
+////////////////////////////////////////////////////////////////////////////////
+// Construct Method to be pass to the DetectorFactory //
+////////////////////////////////////////////////////////////////////////////////
+NPL::VAnalysis* Analysis::Construct(){
+ return (NPL::VAnalysis*) new Analysis();
+}
+
+////////////////////////////////////////////////////////////////////////////////
+// Registering the construct method to the factory //
+////////////////////////////////////////////////////////////////////////////////
+extern "C"{
+ class proxy{
+ public:
+ proxy(){
+ NPL::AnalysisFactory::getInstance()->SetConstructor(Analysis::Construct);
+ }
+ };
+
+ proxy p;
+}
+
diff --git a/Projects/Inelastic/Analysis.h b/Projects/Inelastic/Analysis.h
new file mode 100644
index 000000000..0d321075e
--- /dev/null
+++ b/Projects/Inelastic/Analysis.h
@@ -0,0 +1,67 @@
+#ifndef Analysis_h
+#define Analysis_h
+/*****************************************************************************
+ * 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 : march 2025 *
+ * Last update : *
+ *---------------------------------------------------------------------------*
+ * Decription: *
+ * Class describing the property of an Analysis object *
+ * *
+ *---------------------------------------------------------------------------*
+ * Comment: *
+ * *
+ * *
+ *****************************************************************************/
+#include "NPVAnalysis.h"
+#include "TVendetaPhysics.h"
+#include "TInitialConditions.h"
+#include "TInteractionCoordinates.h"
+#include "TReactionConditions.h"
+#include "NPNucleus.h"
+#include "NPReaction.h"
+#include "TRandom3.h"
+
+class Analysis: public NPL::VAnalysis{
+ public:
+ Analysis();
+ ~Analysis();
+
+ public:
+ void Init();
+ void TreatEvent();
+ void End();
+ void InitOutputBranch();
+ void InitInputBranch();
+ void ReInitValue();
+ static NPL::VAnalysis* Construct();
+
+private:
+ double Einit;
+ vector<double> ThetaLab;
+ vector<double> Ex;
+ vector<double> Elab;
+ vector<double> T;
+ vector<double> Distance;
+ vector<int> Det;
+
+ private:
+ TVendetaPhysics* m_Vendeta;
+ TInitialConditions* InitialConditions;
+ TInteractionCoordinates* InteractionCoordinates;
+ TReactionConditions* ReactionConditions;
+
+ NPL::Reaction* my_Reaction;
+ NPL::Nucleus* neutron;
+
+ TRandom3 ra;
+};
+#endif
diff --git a/Projects/Inelastic/CMakeLists.txt b/Projects/Inelastic/CMakeLists.txt
new file mode 100644
index 000000000..22c74affd
--- /dev/null
+++ b/Projects/Inelastic/CMakeLists.txt
@@ -0,0 +1,5 @@
+cmake_minimum_required (VERSION 2.8)
+# Setting the policy to match Cmake version
+cmake_policy(VERSION ${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION})
+# include the default NPAnalysis cmake file
+include("../../NPLib/ressources/CMake/NPAnalysis.cmake")
diff --git a/Projects/Inelastic/Geometry/GenerateGeometry.C b/Projects/Inelastic/Geometry/GenerateGeometry.C
new file mode 100644
index 000000000..49a319934
--- /dev/null
+++ b/Projects/Inelastic/Geometry/GenerateGeometry.C
@@ -0,0 +1,47 @@
+double R = 2500; // mm
+
+void GenerateGeometry()
+{
+
+ ofstream ofile;
+ string ofilename = "Vendeta.detector";
+ ofile.open(ofilename.c_str());
+ ofile << "\%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" << endl;
+
+ ofile << "Target " << endl;
+ ofile << " THICKNESS= 10 micrometer" << endl;
+ ofile << " RADIUS= 20 mm" << endl;
+ ofile << " MATERIAL= 238U " << endl;
+ ofile << " ANGLE= 0 deg " << endl;
+ ofile << " X= 0 mm " << endl;
+ ofile << " Y= 0 mm " << endl;
+ ofile << " Z= 0 mm " << endl;
+
+ double Theta = 20;
+ double Phi[6] = {0, 20, 40, 140, 160, 180};
+ for(int p=0; p<6; p++){
+ if(p%2==0) Theta = 20;
+ else Theta = 15;
+ cout << p << " " << Theta << endl;
+ for(int i = 0; i <15; i++){
+ ofile << "\%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" << endl;
+ ofile << "Vendeta" << endl;
+ ofile << " R= " << R << " mm" << endl;
+ ofile << " THETA= " << Theta << " deg" << endl;
+ ofile << " PHI= " << Phi[p] << " deg" << endl;
+
+ //ofile << "\%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" << endl;
+ //ofile << "Vendeta" << endl;
+ //ofile << " R= " << R << " mm" << endl;
+ //ofile << " THETA= " << -Theta << " deg" << endl;
+ //ofile << " PHI= " << Phi << " deg" << endl;
+
+
+ Theta += 10;
+ }
+ }
+
+ ofile.close();
+
+
+}
diff --git a/Projects/Inelastic/PhysicsListOption.txt b/Projects/Inelastic/PhysicsListOption.txt
new file mode 100644
index 000000000..ec1a9581f
--- /dev/null
+++ b/Projects/Inelastic/PhysicsListOption.txt
@@ -0,0 +1,15 @@
+EmPhysicsList Option4
+DefaultCutOff 1
+DriftElectronPhysics 0
+IonBinaryCascadePhysics 0
+NPIonInelasticPhysics 0
+EmExtraPhysics 0
+HadronElasticPhysics 0
+StoppingPhysics 0
+OpticalPhysics 0
+HadronPhysicsINCLXX 0
+HadronPhysicsQGSP_BIC_HP 0
+HadronPhysicsQGSP_BERT_HP 0
+Decay 0
+Menate_R 0
+NeutronHP 1
diff --git a/Projects/Inelastic/RunToTreat.txt b/Projects/Inelastic/RunToTreat.txt
new file mode 100644
index 000000000..b1d637dcb
--- /dev/null
+++ b/Projects/Inelastic/RunToTreat.txt
@@ -0,0 +1,4 @@
+TTreeName
+ SimulatedTree
+RootFileName
+ ../../Outputs/Simulation/vendeta_sim_*.root
diff --git a/Projects/Inelastic/Vendeta_inelastic.detector b/Projects/Inelastic/Vendeta_inelastic.detector
new file mode 100644
index 000000000..777bad378
--- /dev/null
+++ b/Projects/Inelastic/Vendeta_inelastic.detector
@@ -0,0 +1,459 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Target
+ THICKNESS= 10 micrometer
+ RADIUS= 20 mm
+ MATERIAL= 238U
+ ANGLE= 0 deg
+ X= 0 mm
+ Y= 0 mm
+ Z= 0 mm
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 20 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 30 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 40 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 50 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 60 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 70 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 80 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 90 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 100 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 110 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 120 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 130 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 140 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 150 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 160 deg
+ PHI= 0 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 15 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 25 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 35 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 45 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 55 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 65 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 75 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 85 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 95 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 105 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 115 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 125 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 135 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 145 deg
+ PHI= 20 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 155 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= 15 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 25 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 35 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 45 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 55 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 65 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 75 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 85 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 95 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 105 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 115 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 125 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 135 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 145 deg
+ PHI= 140 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 155 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= 15 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 25 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 35 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 45 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 55 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 65 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 75 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 85 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 95 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 105 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 115 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 125 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 135 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 145 deg
+ PHI= 180 deg
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Vendeta
+ R= 2500 mm
+ THETA= 155 deg
+ PHI= 180 deg
diff --git a/Projects/Inelastic/run.mac b/Projects/Inelastic/run.mac
new file mode 100644
index 000000000..506556ed2
--- /dev/null
+++ b/Projects/Inelastic/run.mac
@@ -0,0 +1 @@
+/run/beamOn 100000
diff --git a/Projects/Inelastic/sim.sh b/Projects/Inelastic/sim.sh
new file mode 100755
index 000000000..9fc763e80
--- /dev/null
+++ b/Projects/Inelastic/sim.sh
@@ -0,0 +1,7 @@
+#!/bin/bash
+
+for i in {1..10}
+do
+ npsimulation -D Vendeta_inelastic.detector -E 238Uel.reaction --random-seed $i -O vendeta_el_$i -B run.mac &
+done
+
--
GitLab