43 #include "Analysis.hh"
84 if (edep==0.)
return false;
108 <<
"-------->Hits Collection: in this event they are "
110 <<
" hits in the target volume " <<
G4endl;
128 G4int nbSampling = 1000;
133 for (
G4int k=0; k<nbSampling; k++ )
149 G4ThreeVector hitPos = (*fHitsCollection)[randHit]->GetPos();
168 for (
G4int i=0; i<nofHits; i++ )
172 (localPos.x()-randCenterPos.x()) * (localPos.x()-randCenterPos.x()) +
173 (localPos.y()-randCenterPos.y()) * (localPos.y()-randCenterPos.y()) +
174 (localPos.z()-randCenterPos.z()) * (localPos.z()-randCenterPos.z())
179 epsilon = epsilon + (*fHitsCollection)[i]->GetEdep() ;
211 analysisManager->FillNtupleDColumn(0, radius/
nm);
212 analysisManager->FillNtupleDColumn(1, (epsilon/
eV)/(chord/
nm));
213 analysisManager->AddNtupleRow();
G4String SensitiveDetectorName
void SetTrackID(G4int track)
virtual G4bool ProcessHits(G4Step *step, G4TouchableHistory *history)
G4THitsCollection< TrackerHit > TrackerHitsCollection
Definition of the TrackerSD class.
virtual void EndOfEvent(G4HCofThisEvent *hitCollection)
G4int GetCollectionID(G4String colName)
CLHEP::Hep3Vector G4ThreeVector
G4ThreeVector G4RandomDirection()
void SetPos(G4ThreeVector xyz)
virtual void Initialize(G4HCofThisEvent *hitCollection)
G4GLOB_DLL std::ostream G4cout
void insert(G4String str)
const G4ThreeVector & GetPosition() const
ExG4HbookAnalysisManager G4AnalysisManager
TrackerSD(const G4String &name, const G4String &hitsCollectionName)
G4double GetTotalEnergyDeposit() const
void AddHitsCollection(G4int HCID, G4VHitsCollection *aHC)
TrackerHitsCollection * fHitsCollection
G4StepPoint * GetPostStepPoint() const
static G4SDManager * GetSDMpointer()
void SetEdep(G4double de)
G4CollectionNameVector collectionName
G4Track * GetTrack() const
double epsilon(double density, double temperature)