GGSIntHit.h

Go to the documentation of this file.

/*
 * GGSIntHit.h
 *
 *  Created on: 2010-10-17
 *      Authors: Franz Longo , Elena Vannuccini
 */

#ifndef GGSINTHIT_H
#define GGSINTHIT_H

#include "G4Allocator.hh"
#include "G4THitsCollection.hh"
#include "G4ThreeVector.hh"
#include "G4VHit.hh"
#include "G4VPhysicalVolume.hh"

#include "montecarlo/scoring/GGSPartHit.h"
#include "utils/GGSFactory.h"
#include "utils/GGSMiscUtils.h"

#include <numeric>

class GGSIntHit : public G4VHit {
public:
  GGSIntHit(const std::vector<G4double> &timeBins);

  ~GGSIntHit();

  GGSIntHit(const GGSIntHit &right);

  const GGSIntHit &operator=(const GGSIntHit &right);

  int operator==(const GGSIntHit &right) const;

  inline void *operator new(size_t size);

  inline void operator delete(void *aHit);

  virtual void UserInit(G4Step *aStep) { GGS_IGNORE_VAR(aStep); }

  virtual void AddStep(const G4Step &step);

  void SetVolume(const G4VPhysicalVolume *volume) { _volume = volume; }

  void SetPosHitsStorage(bool flag);

  void SetPartHitsStorage(bool flag);

  void SetAbsTranslation(const G4ThreeVector &pos) { _absTranslation = pos; }

  void SetID(G4int id) { _id = id; }

  G4double GetTotalEnergyDep() const { return accumulate(_eDep.begin(), _eDep.end(), 0.); }

  const std::vector<G4double> &GetEnergyDep() const { return _eDep; }

  G4double GetTime() const { return _time; }

  const std::vector<G4double> &GetTimeBins() const { return _timeBins; }

  G4double GetPathLength() const { return _pathLength; }

  const G4VPhysicalVolume *GetVolume() { return _volume; }

  GGSPartHitsCollection *GetPartHits() { return _partHits; }

  const G4ThreeVector &GetAbsTranslation() { return _absTranslation; }

  G4int GetID() { return _id; }

protected:
  std::vector<G4double> _eDep;     // Deposited energy
  std::vector<G4double> _timeBins; // Time bins
  G4double _time;                  // Time of first energy release
  G4double _pathLength;            // Path length of the Hit
  const G4VPhysicalVolume *_volume;
  G4ThreeVector _absTranslation;
  G4int _id;

  bool _isReset;
  GGSPartHitsCollection *_partHits;
  bool _storePosHits;
};

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
typedef G4THitsCollection<GGSIntHit> GGSIntHitsCollection;

extern G4Allocator<GGSIntHit> *GGSIntHitAllocator;

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

inline void *GGSIntHit::operator new(size_t) {
  void *aHit;
  aHit = (void *)GGSIntHitAllocator->MallocSingle();
  return aHit;
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

inline void GGSIntHit::operator delete(void *aHit) { GGSIntHitAllocator->FreeSingle((GGSIntHit *)aHit); }

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

#include "utils/GGSFactory.h"

using GGSIntHitFactory = GGSFactory<GGSIntHit, const std::vector<G4double> &>;

#define RegisterIntHit(className)                                                                                      \
  GGSIntHit *className##Builder(const std::vector<G4double> &timeBins) { return new className(timeBins); }             \
  class className##Proxy {                                                                                             \
  public:                                                                                                              \
    className##Proxy() { GGSIntHitFactory::GetInstance().RegisterBuilder(#className, className##Builder); }            \
  };                                                                                                                   \
  className##Proxy proxyFor##className
#endif /* GGSINTHIT_H */