From 9b857100ccb342b6f3da3e3e3d237af44e58a5fc Mon Sep 17 00:00:00 2001
From: matta <matta@npt>
Date: Tue, 30 Mar 2010 08:26:53 +0000
Subject: [PATCH] * Fixing bug in MUST2 Geometry - No strip
inversion...
---
Inputs/EventGenerator/10He.reaction | 2 +-
NPLib/MUST2/TMust2Physics.cxx | 3 +--
NPSimulation/src/MUST2Array.cc | 4 ++--
3 files changed, 4 insertions(+), 5 deletions(-)
diff --git a/Inputs/EventGenerator/10He.reaction b/Inputs/EventGenerator/10He.reaction
index 138ab2f5c..2cbb5f098 100644
--- a/Inputs/EventGenerator/10He.reaction
+++ b/Inputs/EventGenerator/10He.reaction
@@ -19,7 +19,7 @@ TransfertToResonance
ResonanceDecayA= 8
CrossSectionPath= 11Li(d,3He)10He.txt
ShootLight= 1
- ShootHeavy= 0
+ ShootHeavy= 1
ShootDecayProduct= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
diff --git a/NPLib/MUST2/TMust2Physics.cxx b/NPLib/MUST2/TMust2Physics.cxx
index 4bf1c25c2..e64af68be 100644
--- a/NPLib/MUST2/TMust2Physics.cxx
+++ b/NPLib/MUST2/TMust2Physics.cxx
@@ -799,9 +799,8 @@ void TMust2Physics::AddTelescope( TVector3 C_X1_Y1 ,
for( int j = 0 ; j < 128 ; j++ )
{
- StripCenter = Strip_1_1 + StripPitch*( j*U + i*V ) ;
+ StripCenter = Strip_1_1 + StripPitch*( i*U + j*V ) ;
//StripCenter += -TargetPosition ;
-
lineX.push_back( StripCenter.X() ) ;
lineY.push_back( StripCenter.Y() ) ;
lineZ.push_back( StripCenter.Z() ) ;
diff --git a/NPSimulation/src/MUST2Array.cc b/NPSimulation/src/MUST2Array.cc
index 17825a63f..5f2a6ced9 100644
--- a/NPSimulation/src/MUST2Array.cc
+++ b/NPSimulation/src/MUST2Array.cc
@@ -879,10 +879,10 @@ void MUST2Array::ConstructDetector(G4LogicalVolume* world)
// (u,v) // to silicon plan
// w perpendicular to (u,v) plan and pointing CsI
MMu = m_X128_Y1[i] - m_X1_Y1[i] ;
- MMu = MMu.unit() ;
+ MMu = -MMu.unit() ;
MMv = m_X1_Y128[i] - m_X1_Y1[i] ;
- MMv = MMv.unit() ;
+ MMv = -MMv.unit() ;
MMw = MMv.cross(MMu) ;
// if (MMw.z() > 0)MMw = MMv.cross(MMu) ;
--
GitLab