HFShowerG4Hit.h

Go to the documentation of this file.

#ifndef SimG4CMS_ShowerLibraryProducer_HFShowerG4Hit_h
#define SimG4CMS_ShowerLibraryProducer_HFShowerG4Hit_h

#include "SimDataFormats/CaloHit/interface/HFShowerPhoton.h"
#include "DataFormats/Math/interface/Point3D.h"

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

#include <vector>

class HFShowerG4Hit : public G4VHit {
public:
  HFShowerG4Hit();
  HFShowerG4Hit(G4int hitId, G4int tkID, double edep, double time);
  ~HFShowerG4Hit() override;
  HFShowerG4Hit(const HFShowerG4Hit& right);
  const HFShowerG4Hit& operator=(const HFShowerG4Hit& right);
  G4int operator==(const HFShowerG4Hit& right) const;

  inline void* operator new(size_t);
  inline void operator delete(void* aHit);

private:
  G4int theHitId;
  G4int theTrackId;
  G4double theEdep;
  G4double theTime;
  G4ThreeVector localPos;
  G4ThreeVector globalPos;
  G4ThreeVector momDir;

public:
  inline void setHitId(G4int hitId) { theHitId = hitId; }
  inline void setTrackId(G4int trackId) { theTrackId = trackId; }
  inline void setEnergy(G4double edep) { theEdep = edep; }
  inline void updateEnergy(G4double edep) { theEdep += edep; }
  inline void setTime(G4double t) { theTime = t; }
  inline void setLocalPos(const G4ThreeVector& xyz) { localPos = xyz; }
  inline void setGlobalPos(const G4ThreeVector& xyz) { globalPos = xyz; }
  inline void setPrimMomDir(const G4ThreeVector& xyz) { momDir = xyz; }

  inline G4int hitId() const { return theHitId; }
  inline G4int trackId() const { return theTrackId; }
  inline G4double edep() const { return theEdep; };
  inline G4double time() const { return theTime; }
  inline G4ThreeVector localPosition() const { return localPos; }
  inline G4ThreeVector globalPosition() const { return globalPos; }
  inline G4ThreeVector primaryMomDir() const { return momDir; }
};

typedef G4THitsCollection<HFShowerG4Hit> HFShowerG4HitsCollection;

extern G4ThreadLocal G4Allocator<HFShowerG4Hit>* fHFShowerG4HitAllocator;

inline void* HFShowerG4Hit::operator new(size_t) {
  if (!fHFShowerG4HitAllocator)
    fHFShowerG4HitAllocator = new G4Allocator<HFShowerG4Hit>;
  return (void*)fHFShowerG4HitAllocator->MallocSingle();
}

inline void HFShowerG4Hit::operator delete(void* aHit) { fHFShowerG4HitAllocator->FreeSingle((HFShowerG4Hit*)aHit); }
#endif