Commit 2ebb9a4b authored by Morfouace's avatar Morfouace
Browse files

* updating scone detector

parent e6cdd725
Pipeline #95666 passed with stages
in 13 minutes and 1 second
......@@ -568,7 +568,7 @@ void Scone::ReadSensitive(const G4Event* ){
m_Event->SetProtonTime(proton_time[i]);
}
Process_scorer->clear();
//Process_scorer->clear();
///////////
// Process scorer for Gd
......@@ -583,25 +583,29 @@ void Scone::ReadSensitive(const G4Event* ){
GdCaptureTime = GdProcess_scorer->GetProcessTime(i);
}
}
if(kPlasticCapture){
m_Event->SetCapture(1);
m_Event->SetCaptureTime(PlasticCaptureTime);
}
else if(kGdCapture){
if(kGdCapture){
m_Event->SetCapture(2);
m_Event->SetCaptureTime(GdCaptureTime);
}
else if(kPlasticCapture){
m_Event->SetCapture(1);
m_Event->SetCaptureTime(PlasticCaptureTime);
}
else m_Event->SetCapture(0);
GdProcess_scorer->clear();
///////////
// Process scorer for fission chamber
ProcessScorers::PS_Process* FCProcess_scorer = (ProcessScorers::PS_Process*) m_FCScorer->GetPrimitive(0);
vector<int> FC_process = FCProcess_scorer->GetFCProcess();
for(unsigned int i=0; i<FC_process.size(); i++){
m_Event->SetFCProcess(FC_process[i]);
if(m_BuildFissionChamber==1){
ProcessScorers::PS_Process* FCProcess_scorer = (ProcessScorers::PS_Process*) m_FCScorer->GetPrimitive(0);
vector<int> FC_process = FCProcess_scorer->GetFCProcess();
for(unsigned int i=0; i<FC_process.size(); i++){
m_Event->SetFCProcess(FC_process[i]);
}
FCProcess_scorer->clear();
}
FCProcess_scorer->clear();
}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
......
......@@ -39,12 +39,12 @@ G4bool PS_Process::ProcessHits(G4Step* aStep, G4TouchableHistory*){
G4String processname;
if(aStep->GetPostStepPoint()->GetProcessDefinedStep() != NULL){
processname = aStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName();
t_processname.push_back(processname);
t_processtime.push_back(aStep->GetPreStepPoint()->GetGlobalTime());
}
G4int trackID;
/*G4int trackID;
//G4int parentID;
G4String volumename, particlename;
//G4double step_length;
......@@ -76,9 +76,7 @@ G4bool PS_Process::ProcessHits(G4Step* aStep, G4TouchableHistory*){
}
}
}
}
}*/
return TRUE;
}
......
......@@ -76,7 +76,8 @@ void Analysis::TreatEvent(){
Det.push_back(m_ChiNu->DetectorNumber[i]);
T.push_back(m_ChiNu->Time[i]);
double T_stop = m_ChiNu->Time[i]*1e-9;
neutron->SetTimeOfFlight((T_stop-beam_TOF)/(Rdet*1e-3));
//neutron->SetTimeOfFlight((T_stop-beam_TOF)/(Rdet*1e-3));
neutron->SetTimeOfFlight((T_stop)/(Rdet*1e-3));
E.push_back(m_ChiNu->Energy[i]);
Elab.push_back(neutron->GetEnergy());
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment