root graphs being developed

parent 340c629a
...@@ -79,7 +79,7 @@ void B1SteppingAction::UserSteppingAction(const G4Step* step) ...@@ -79,7 +79,7 @@ void B1SteppingAction::UserSteppingAction(const G4Step* step)
//ATENÇÃO - Comments dentro dos else-if -> para histograma de particlas x z (BUG no Z=0 e resultado duvidoso -- MELHORAR) //ATENÇÃO - Comments dentro dos else-if -> para histograma de particlas x z (BUG no Z=0 e resultado duvidoso -- MELHORAR)
G4AnalysisManager* analysisManager = G4AnalysisManager::Instance(); G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
//collect particle for each step //collect particle for each step
G4String particle = step->GetTrack()->GetDynamicParticle()->GetDefinition() ->GetParticleName(); G4String particle = step->GetTrack()->GetDynamicParticle()->GetDefinition() ->GetParticleName();
G4double z =step->GetPreStepPoint()->GetPosition().z()/meter; G4double z =step->GetPreStepPoint()->GetPosition().z()/meter;
...@@ -112,7 +112,6 @@ G4AnalysisManager* analysisManager = G4AnalysisManager::Instance(); ...@@ -112,7 +112,6 @@ G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
} }
//Retire esse comment para histograma que leva em conta o numero de vezes que a particula foi identificado no evento. -> Por enquanto o mais confiável. //Retire esse comment para histograma que leva em conta o numero de vezes que a particula foi identificado no evento. -> Por enquanto o mais confiável.
/* /*
G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
analysisManager->FillNtupleDColumn(0, p); analysisManager->FillNtupleDColumn(0, p);
analysisManager->FillNtupleDColumn(1, e_p); analysisManager->FillNtupleDColumn(1, e_p);
analysisManager->FillNtupleDColumn(2, e_minus); analysisManager->FillNtupleDColumn(2, e_minus);
...@@ -124,7 +123,7 @@ G4AnalysisManager* analysisManager = G4AnalysisManager::Instance(); ...@@ -124,7 +123,7 @@ G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
analysisManager->FillNtupleDColumn(8, z); analysisManager->FillNtupleDColumn(8, z);
/**/ /**/
analysisManager->AddNtupleRow(); analysisManager->AddNtupleRow();
// collect energy deposited in this step // collect energy deposited in this step
G4double edepStep = step->GetTotalEnergyDeposit(); G4double edepStep = step->GetTotalEnergyDeposit();
fEventAction->AddEdep(edepStep); fEventAction->AddEdep(edepStep);
......
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