//SteppingAction
#include "Constants.hh"
using namespace ConstantSpace;

//Non funziona per via del Multi-threading
if (volume == fSD) {
	G4double ip,fractpart,jDET,iDET;
			
	G4double Lx = SCsize*0.5;
	G4double Ly = SCsize*0.5;
	
	G4double Dx = 2.*Lx/NPixelX;
	G4double Dy = 2.*Ly/NPixelY;
	
	G4ThreeVector pos = preStepPoint->GetPosition();
	
	ip = (pos.y()+Ly)/Dy;
 	fractpart = modf(ip, &iDET);
	if (fractpart>0.5) iDET = iDET+1;
	 
	ip = (pos.x()+Lx)/Dx;
	fractpart = modf(ip, &jDET);
 	if (fractpart>0.5) jDET = jDET+1;
 	
 	G4cout << int(iDET) << " " << int(jDET) << " " << preStepPoint->GetKineticEnergy() << G4endl;
   	
   	//fRunAction->UpdateRunScoringMatrix(int(iDET), int(jDET), preStepPoint->GetKineticEnergy());
   	fRunAction->UpdateRunScoringMatrix(0, 0, 1111);
   	
   	track->SetTrackStatus(fStopAndKill); 
}	


//RunAction
#include "Constants.hh"
using namespace ConstantSpace;

//initialization of the ScoringMatrix
for (unsigned int i=0; i<NPixelY; i++) {
	for (unsigned int j=0; j<NPixelX; j++) {
		RunScoringMatrix[i][j] = 0.;
	}
}	

//Write the scoring result of this Run
//if (IsMaster()) {
	std::ofstream fFileSM;
	fFileSM.open("ScoringMatrix.dat", std::ofstream::out);
	
	for (unsigned int i=0; i<NPixelY; i++) {
		for (unsigned int j=0; j<NPixelX; j++) {
			fFileSM << RunScoringMatrix[i][j] << " ";
		}
		fFileSM << G4endl;
	}	

	fFileSM.close();
//}


//costants.hh
#ifndef Constants_h
#define Constants_h 1

#include "G4Types.hh"
#include "G4SystemOfUnits.hh"
#include "G4UnitsTable.hh"

namespace ConstantSpace {

  //number of pixel of the Scoring Screen 
  const G4int NPixelX = 100;
  const G4int NPixelY = 100;

}

#endif