gwm17/Mask
1
0
Fork 0
mirror of https://github.com/gwm17/Mask.git synced 2024-11-22 10:18:50 -05:00
Code Issues Projects Releases Wiki Activity

Changed executable name

Browse Source
This commit is contained in:
Gordon McCann 2021-03-16 09:46:29 -04:00
parent 27712396c6
commit 5785c0535b
3 changed files with 1 additions and 151 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

46
input.txt.orig
View File

@ -1,46 +0,0 @@
----------Data Information----------
OutputFile: /data1/gwm17/test_dead.root
SaveTree: yes
SavePlots: yes
----------Reaction Information----------
ReactionType: 2
Z A (order is target, projectile, ejectile, break1, break3)
<<<<<<< HEAD
6 12
2 3
1 2
1 1
=======
5 10
2 3
2 4
1 2
>>>>>>> 4ccaabb534f35f2cf36d375b74ac4ebf99fe7bb7
----------Target Information----------
Name: test_targ
Layers: 1
~Layer1
<<<<<<< HEAD
Thickness(ug/cm^2): 40
=======
Thickness(ug/cm^2): 10
>>>>>>> 4ccaabb534f35f2cf36d375b74ac4ebf99fe7bb7
Z A Stoich
6 12 1
0
~
<<<<<<< HEAD
=======
~Layer2
Thickness(ug/cm^2): 80
Z A Stoich
5 10 1
0
~
>>>>>>> 4ccaabb534f35f2cf36d375b74ac4ebf99fe7bb7
----------Sampling Information----------
NumberOfSamples: 1000000
BeamMeanEnergy(MeV): 24 BeamEnergySigma(MeV): 0.001
EjectileThetaMin(deg): 3.0 EjectileThetaMax(deg): 3.0
ResidualExMean(MeV): 2.364 ResidualExSigma(MeV): 0.0317
--------------------------------------

2
makefile
View File

@ -17,7 +17,7 @@ DICTOBJ=$(OBJDIR)/kinematics_dict.o
DICTSRC=$(SRCDIR)/kinematics_dict.cxx
DICT_PAGES=$(INCLDIR)/Kinematics.h $(INCLDIR)/LinkDef_Kinematics.h
EXE=$(BINDIR)/kinematics
EXE=$(BINDIR)/mask
CLEANUP=$(EXE) $(OBJS) $(DICTOBJ) $(DICTSRC)

104
src/Target.cpp.orig
View File

@ -1,104 +0,0 @@
/*
Target.cpp
A basic target unit for use in the SPANCRedux environment. A target
is defined as a single compound with elements Z,A of a given stoichiometry
Holds an energy loss class
Based on code by D.W. Visser written at Yale for the original SPANC
Written by G.W. McCann Aug. 2020
*/
#include "Target.h"
/*Targets must be of known thickness*/
Target::Target(double thick) {
thickness = thick;
}
Target::~Target() {
}
/*Set target elements of given Z, A, S*/
void Target::SetElements(std::vector<int>& z, std::vector<int>& a, std::vector<int>& stoich) {
Z = z;
A = a;
Stoich = stoich;
eloss.SetTargetComponents(Z, A, Stoich);
}
/*Element verification*/
bool Target::ContainsElement(int z, int a) {
for(unsigned int i=0; i<Z.size(); i++) {
if( z == Z[i] && a == A[i]) return true;
}
return false;
}
/*Calculates energy loss for travelling all the way through the target*/
double Target::getEnergyLossTotal(int zp, int ap, double startEnergy, double theta) {
if(theta == PI/2.) return startEnergy;
<<<<<<< HEAD
else if(theta > PI/2.) theta = PI - theta;
=======
else if (theta > PI/2.) theta = PI - theta;
>>>>>>> 4ccaabb534f35f2cf36d375b74ac4ebf99fe7bb7
return eloss.GetEnergyLoss(zp, ap, startEnergy, thickness/fabs(cos(theta)));
}
/*Calculates energy loss for travelling halfway through the target*/
double Target::getEnergyLossHalf(int zp, int ap, double startEnergy, double theta) {
if(theta == PI/2.) return startEnergy;
<<<<<<< HEAD
else if(theta > PI/2.) theta = PI - theta;
=======
else if (theta > PI/2.) theta = PI - theta;
>>>>>>> 4ccaabb534f35f2cf36d375b74ac4ebf99fe7bb7
return eloss.GetEnergyLoss(zp, ap, startEnergy, thickness/(2.0*fabs(cos(theta))));
}
/*Calculates reverse energy loss for travelling all the way through the target*/
double Target::getReverseEnergyLossTotal(int zp, int ap, double finalEnergy, double theta) {
if(theta == PI/2.) return finalEnergy;
<<<<<<< HEAD
else if(theta > PI/2.) theta = PI - theta;
=======
else if (theta > PI/2.) theta = PI - theta;
>>>>>>> 4ccaabb534f35f2cf36d375b74ac4ebf99fe7bb7
return eloss.GetReverseEnergyLoss(zp, ap, finalEnergy, thickness/fabs(cos(theta)));
}
/*Calculates reverse energy loss for travelling half way through the target*/
double Target::getReverseEnergyLossHalf(int zp, int ap, double finalEnergy, double theta) {
if(theta == PI/2.) return finalEnergy;
<<<<<<< HEAD
else if(theta > PI/2.) theta = PI - theta;
=======
else if (theta > PI/2.) theta = PI - theta;
>>>>>>> 4ccaabb534f35f2cf36d375b74ac4ebf99fe7bb7
return eloss.GetReverseEnergyLoss(zp, ap, finalEnergy, thickness/(2.0*fabs(cos(theta))));
}
/*Getter functions*/
double& Target::GetThickness() {
return thickness;
}
int Target::GetNumberOfElements() {
return Z.size();
}
int Target::GetElementZ(int index) {
return Z[index];
}
int Target::GetElementA(int index) {
return A[index];
}
int Target::GetElementStoich(int index) {
return Stoich[index];
}