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 "G4VHit.hh"
#include "G4THitsCollection.hh"
#include "G4Allocator.hh"
#include "G4ThreeVector.hh"
#include "G4VPhysicalVolume.hh"

#include "montecarlo/scoring/GGSPartHit.h"
#include "utils/GGSFactory.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 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 */