From 5e56ec75c0c0dafd82c864491d6b232528311db0 Mon Sep 17 00:00:00 2001
From: adrien matta <matta@lpccaen.in2p3.fr>
Date: Mon, 27 Jul 2020 17:14:56 +0200
Subject: [PATCH] *adjusting CATANA resolution - assuming proton reso
is similar to MUST2, to be checked
---
NPSimulation/Detectors/Catana/Catana.cc | 2 +-
Projects/Strasse/Analysis.cxx | 4 +---
Projects/Strasse/reaction/C12_p2p.reaction | 10 +++++-----
3 files changed, 7 insertions(+), 9 deletions(-)
diff --git a/NPSimulation/Detectors/Catana/Catana.cc b/NPSimulation/Detectors/Catana/Catana.cc
index 6373846ec..bc634efbe 100644
--- a/NPSimulation/Detectors/Catana/Catana.cc
+++ b/NPSimulation/Detectors/Catana/Catana.cc
@@ -66,7 +66,7 @@ namespace Catana_NS{
// Energy and time Resolution
const double EnergyThreshold = 0.1*MeV;
const double ResoTime = 4.5*ns ;
- const double ResoEnergy = 1.0*MeV ;
+ const double ResoEnergy = 0.08*MeV ;
double Length = 600*mm ;
string Material = "CsI";
}
diff --git a/Projects/Strasse/Analysis.cxx b/Projects/Strasse/Analysis.cxx
index 92dc2e844..1f5acb933 100644
--- a/Projects/Strasse/Analysis.cxx
+++ b/Projects/Strasse/Analysis.cxx
@@ -109,7 +109,6 @@ void Analysis::TreatEvent(){
E1 = ReconstructProtonEnergy(Vertex,Proton1,Catana->Energy[i1]);
E2 = ReconstructProtonEnergy(Vertex,Proton2,Catana->Energy[i2]);
double TA = BeamTarget.Slow(InitialBeamEnergy,abs(VertexZ-75),0);
- //double TA = RC->GetBeamEnergy();
// 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
Proton1=E1*Proton1.Unit();
@@ -126,7 +125,6 @@ void Analysis::TreatEvent(){
LV_B = LV_A + LV_T - LV_p1 - LV_p2;
//LV_B = RC->GetParticleMomentum(2);
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
// Strasse Chamber
E = protonAl.EvaluateInitialEnergy(E,3*mm,Theta);
// Outer Barrel
- E = protonSi.EvaluateInitialEnergy(E,400*micrometer,Theta);
+ E = protonSi.EvaluateInitialEnergy(E,300*micrometer,Theta);
// Inner Barrel
E = protonSi.EvaluateInitialEnergy(E,200*micrometer,Theta);
// LH2 target
diff --git a/Projects/Strasse/reaction/C12_p2p.reaction b/Projects/Strasse/reaction/C12_p2p.reaction
index 97a681ccc..47d81151a 100755
--- a/Projects/Strasse/reaction/C12_p2p.reaction
+++ b/Projects/Strasse/reaction/C12_p2p.reaction
@@ -4,10 +4,10 @@ Beam
Particle= 12C
Energy= 4800 MeV
SigmaEnergy= 0 MeV
- SigmaThetaX= 0 deg
- SigmaPhiY= 0 deg
- SigmaX= 0 mm
- SigmaY= 0 mm
+ SigmaThetaX= 0.1 deg
+ SigmaPhiY= 0.1 deg
+ SigmaX= 5 mm
+ SigmaY= 5 mm
MeanThetaX= 0 deg
MeanPhiY= 0 deg
MeanX= 0 mm
@@ -22,7 +22,7 @@ QFSReaction
Heavy= 11B
ExcitationEnergyBeam= 0.0 MeV
ExcitationEnergyHeavy= 0.0 MeV
- MomentumSigma= 0.0
+ MomentumSigma= 50.0
ShootHeavy= 1
ShootLight= 1
--
GitLab