Skip to content
Snippets Groups Projects
Commit 5e56ec75 authored by Adrien Matta's avatar Adrien Matta :skull_crossbones:
Browse files

*adjusting CATANA resolution 

        - assuming proton reso is similar to MUST2, to be checked
parent fefe7880
No related branches found
No related tags found
No related merge requests found
Pipeline #77852 passed
Stage: build-NPLib
Stage: build-NPSimulation
Stage: test
Hide whitespace changes
Inline Side-by-side
NPSimulation/Detectors/Catana/Catana.cc
+ 1
− 1
View file @ 5e56ec75
...@@ -66,7 +66,7 @@ namespace Catana_NS{ ...@@ -66,7 +66,7 @@ namespace Catana_NS{
// Energy and time Resolution // Energy and time Resolution
const double EnergyThreshold = 0.1*MeV; const double EnergyThreshold = 0.1*MeV;
const double ResoTime = 4.5*ns ; const double ResoTime = 4.5*ns ;
const double ResoEnergy = 1.0*MeV ; const double ResoEnergy = 0.08*MeV ;
double Length = 600*mm ; double Length = 600*mm ;
string Material = "CsI"; string Material = "CsI";
} }
......
Projects/Strasse/Analysis.cxx
+ 1
− 3
View file @ 5e56ec75
...@@ -109,7 +109,6 @@ void Analysis::TreatEvent(){ ...@@ -109,7 +109,6 @@ void Analysis::TreatEvent(){
E1 = ReconstructProtonEnergy(Vertex,Proton1,Catana->Energy[i1]); E1 = ReconstructProtonEnergy(Vertex,Proton1,Catana->Energy[i1]);
E2 = ReconstructProtonEnergy(Vertex,Proton2,Catana->Energy[i2]); E2 = ReconstructProtonEnergy(Vertex,Proton2,Catana->Energy[i2]);
double TA = BeamTarget.Slow(InitialBeamEnergy,abs(VertexZ-75),0); double TA = BeamTarget.Slow(InitialBeamEnergy,abs(VertexZ-75),0);
//double TA = RC->GetBeamEnergy();
// setting up Lorentz Vector from measured trajectories and energies // setting up Lorentz Vector from measured trajectories and energies
TVector3 PA(0,0,sqrt(TA*(TA+2*m_QFS->GetNucleusA()->Mass()))); // for like there is no BDC TVector3 PA(0,0,sqrt(TA*(TA+2*m_QFS->GetNucleusA()->Mass()))); // for like there is no BDC
Proton1=E1*Proton1.Unit(); Proton1=E1*Proton1.Unit();
...@@ -126,7 +125,6 @@ void Analysis::TreatEvent(){ ...@@ -126,7 +125,6 @@ void Analysis::TreatEvent(){
LV_B = LV_A + LV_T - LV_p1 - LV_p2; LV_B = LV_A + LV_T - LV_p1 - LV_p2;
//LV_B = RC->GetParticleMomentum(2); //LV_B = RC->GetParticleMomentum(2);
Ex = LV_B.M() - m_QFS->GetNucleusB()->Mass(); Ex = LV_B.M() - m_QFS->GetNucleusB()->Mass();
cout << Ex << " " << m_QFS->GetNucleusB()->Mass() << endl;
} }
} }
...@@ -142,7 +140,7 @@ double Analysis::ReconstructProtonEnergy(const TVector3& x0, const TVector3& dir ...@@ -142,7 +140,7 @@ double Analysis::ReconstructProtonEnergy(const TVector3& x0, const TVector3& dir
// Strasse Chamber // Strasse Chamber
E = protonAl.EvaluateInitialEnergy(E,3*mm,Theta); E = protonAl.EvaluateInitialEnergy(E,3*mm,Theta);
// Outer Barrel // Outer Barrel
E = protonSi.EvaluateInitialEnergy(E,400*micrometer,Theta); E = protonSi.EvaluateInitialEnergy(E,300*micrometer,Theta);
// Inner Barrel // Inner Barrel
E = protonSi.EvaluateInitialEnergy(E,200*micrometer,Theta); E = protonSi.EvaluateInitialEnergy(E,200*micrometer,Theta);
// LH2 target // LH2 target
......
Projects/Strasse/reaction/C12_p2p.reaction
+ 5
− 5
View file @ 5e56ec75
...@@ -4,10 +4,10 @@ Beam ...@@ -4,10 +4,10 @@ Beam
Particle= 12C Particle= 12C
Energy= 4800 MeV Energy= 4800 MeV
SigmaEnergy= 0 MeV SigmaEnergy= 0 MeV
SigmaThetaX= 0 deg SigmaThetaX= 0.1 deg
SigmaPhiY= 0 deg SigmaPhiY= 0.1 deg
SigmaX= 0 mm SigmaX= 5 mm
SigmaY= 0 mm SigmaY= 5 mm
MeanThetaX= 0 deg MeanThetaX= 0 deg
MeanPhiY= 0 deg MeanPhiY= 0 deg
MeanX= 0 mm MeanX= 0 mm
...@@ -22,7 +22,7 @@ QFSReaction ...@@ -22,7 +22,7 @@ QFSReaction
Heavy= 11B Heavy= 11B
ExcitationEnergyBeam= 0.0 MeV ExcitationEnergyBeam= 0.0 MeV
ExcitationEnergyHeavy= 0.0 MeV ExcitationEnergyHeavy= 0.0 MeV
MomentumSigma= 0.0 MomentumSigma= 50.0
ShootHeavy= 1 ShootHeavy= 1
ShootLight= 1 ShootLight= 1
......
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