diff --git a/NPSimulation/Detectors/Catana/Catana.cc b/NPSimulation/Detectors/Catana/Catana.cc index 6373846ec0aedd813086a87732889d29d6925cde..bc634efbeb90c3938d68840dd5fd4defa0cf2455 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 92dc2e8447b973b11b9fe89ede8cc9a75c0e3850..1f5acb9334ba72a1ba566747be1d52f91137d070 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 97a681cccf563e9ae57bda14d0fa2508f51b89ae..47d81151a1fb81d022acb2cef5a7c7b067aae6cf 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