No commit message

No commit message

Inputs/DetectorConfiguration/Riken_65mm.detector

@@ -75,10 +75,10 @@ VIS=	all
 			
 %%%%%%% Telescope 5 %%%%%%%			
 M2Telescope 			
-X1_Y1=		-50.23	-50.23	350
-X1_Y128=	50.23	-50.23	350
-X128_Y1=	-50.23	50.23	350
-X128_Y128=	50.23	50.23	350
+X1_Y1=		-50.23	-50.23	365
+X1_Y128=	50.23	-50.23	365
+X128_Y1=	-50.23	50.23	365
+X128_Y128=	50.23	50.23	365
 SI=	1		
 SILI=	0		
 CSI=	1		
@@ -150,7 +150,7 @@ ScintillatorPlastic
 Plastic
 	THETA= 0
 	PHI= 0
-	R= 303
+	R= 318
 	Thickness= 20
 	Radius= 20
 	Scintillator= BC452_10
@@ -159,7 +159,7 @@ Plastic
 Plastic
 	THETA= 0
 	PHI= 0
-	R= 328
+	R= 343
 	Thickness= 30
 	Radius= 20
 	Scintillator= BC452_10

Inputs/EventGenerator/10He.reaction

@@ -7,7 +7,7 @@ TransfertToResonance
 	Target= 2H
 	Light= 3He
 	Heavy= 10He
-	ExcitationEnergy= 2.0
+	ExcitationEnergy= 4.0
 	BeamEnergy= 550
 	BeamEnergySpread= 0
 	SigmaThetaX= 0.6921330164
@@ -16,7 +16,7 @@ TransfertToResonance
 	SigmaY= 9.069
 	ResonanceWidth= 0
 	ResonanceDecayZ= 2
-	ResonanceDecayA= 8
+	ResonanceDecayA= 6
 	CrossSectionPath= 11Li(d,3He)10He.txt	
 	ShootLight= 1
 	ShootHeavy= 1

NPAnalysis/10He_Riken/src/Analysis.cc

@@ -95,7 +95,7 @@ int main(int argc,char** argv)
 
 			// Minimum code
 			if( i%10000 == 0 && i!=0) 	{	
-											cout.precision(3);
+											cout.precision(5);
 											end=clock();										
 											double TimeElapsed = (end-begin)/CLOCKS_PER_SEC;
 											double percent = (double)i/N ;

NPSimulation/include/MUST2Array.hh

@@ -55,8 +55,9 @@ namespace MUST2
    const G4double MylarCsIThickness = 3*micrometer          ;
    const G4double CsIThickness      = 4.*cm + 2*MylarCsIThickness ;
    const G4double CsIFaceFront      = 12.2*cm                  ;
-   const G4double CsIFaceBack    = 16*cm                    ;
-
+   const G4double CsIFaceBack    	= 16*cm                    ;
+   const G4double DistInterCsI 		= 0.2 * mm;
+   
    // Starting at the front and going to CsI
    const G4double AluStripFront_PosZ   = Length* -0.5 + 0.5*AluStripThickness                       ;
    const G4double Silicon_PosZ      = AluStripFront_PosZ + 0.5*AluStripThickness + 0.5*SiliconThickness     ;

NPSimulation/src/MUST2Array.cc

@@ -382,9 +382,7 @@ void MUST2Array::VolumeMaker(G4int TelescopeNumber ,
 ////////////////////////////////////////////////////////////////
 
    if (wCsI) {
-      //G4double DistInterCsI= 1.*mm;//<-Set the distance between CsI cristal. Exact value yet unknown
-      G4double DistInterCsI = 0.2 * mm;
-
+   
       G4ThreeVector positionCsI = G4ThreeVector(0, 0, CsI_PosZ);
       G4Trd* solidCsI = new G4Trd("csI", 0.5*CsIFaceFront, 0.5*CsIFaceBack, 0.5*CsIFaceFront, 0.5*CsIFaceBack, 0.5*CsIThickness);