LCDG4CalHit.hh

Go to the documentation of this file.

#ifndef LCDG4CalHit_h
#define LCDG4CalHit_h 1

#include "G4VHit.hh"
#include "G4THitsCollection.hh"
#include "G4Allocator.hh"
#include "G4ThreeVector.hh"
#include "G4LogicalVolume.hh"
#include <vector>

class LCDG4CalHit : public G4VHit {

public:
  // default constructor: creates an empty hit
  LCDG4CalHit();
  // destructor
  ~LCDG4CalHit();
  // copy constructor
  LCDG4CalHit(const LCDG4CalHit &right);
  // assignment operator
  const LCDG4CalHit& operator=(const LCDG4CalHit &right);
  // equality comparison operator
  bool operator==(const LCDG4CalHit &right) const;

  // customized new/delete operators
  inline void *operator new(size_t);
  inline void operator delete(void *aHit);

  void Draw();
  void Print();

public:
  // set hit layer index
  inline void SetLayer(const G4int layer) { m_layer = layer; }
  // return hit layer index
  inline G4int GetLayer() const { return m_layer; }
  // set hit SysNo
  inline void SetSystem(const G4int sysNo) { m_sysNo = sysNo; }
  // return hit SysNo
  inline G4int GetSystem() const { return m_sysNo; }
  // set hit barrel/endcap flag
  inline void SetBarEnd(const G4int barend) { m_barend = barend; }
  // return hit barrel/endcap flag
  inline G4int GetBarEnd() const { return m_barend; }
  // set hit north/south flag
  inline void SetNorth(const G4int north) { m_north = north; }
  // return hit north/south flag
  inline G4int GetNorth() const { return m_north; }
  // set hit phi index
  inline void SetPhi(const G4int phi) { m_phi = phi; }
  // return hit phi index
  inline G4int GetPhi() const { return m_phi; }
  // set hit theta index
  inline void SetTheta(const G4int theta) { m_theta = theta; }
  // return hit theta index
  inline G4int GetTheta() const { return m_theta; }
  // set total energy deposited into sensitive material for this hit
  inline void SetEdep(const G4double de) { m_edep = de; }
  // return total energy deposited into sensitive material for this hit
  inline G4double GetEdep() const { return m_edep; }
  // set total energy deposited into absorber for this hit
  inline void SetEdepAbs(const G4double de) { m_edepabs = de; }
  // return total energy deposited into absorber for this hit
  inline G4double GetEdepAbs() const { return m_edepabs; }

  // set ref position of this cell (center)
  inline void SetPos(const G4ThreeVector& xyz) { m_pos = xyz; }
  // return ref position of this cell
  inline G4ThreeVector GetPos() const { return m_pos; }
  // return number of particles contributing to total energy deposited
  inline G4int GetNMC() const {return m_nmc;}
  // return geant particle Id of i-th contributing particle
  inline G4int GetTrkID(G4int ith) const {return (*m_trkid)[ith];}
  // return energy deposited by idx-th contributing particle
  inline G4double GetEMC(G4int ith) const {return (*m_emc)[ith];}
  // return timing of energy deposition by idx-th contributing particle
  inline G4double GetTimeMC(G4int ith) const {return (*m_tmc)[ith];}

  // return cell towerID
  unsigned int GetCellID() const;
  unsigned int sioTagToLcioTag(unsigned int sioTag) const;

public:
  G4int AddLayerLog(G4LogicalVolume* g4vxdlog, G4int startlayerNo);
  // Add energy deposited into sensitive material by particle mcidx
  G4int AddMcPart(G4int mcidx, G4double emc, G4double tmc);
  // Hit on absorber: Add energy emc to energy deposited on absorber
  void  AddMcPartAbs(G4double emc) { m_edepabs += emc; }

private:
  // Layer index, = 1..Nlayers
  G4int    m_layer;
  // sysNo, whatever that means...
  G4int    m_sysNo;
  // barend = 0 for barrel, = 1 for endcap
  G4int    m_barend;
  // north = 0 for south (z>0), = 1 for north (z<0)
  G4int    m_north;
  // phi index, = 1..NphiCells
  G4int    m_phi;
  // (theta|z) index, = 1..N(th|z)Cells, (th|z) means either theta or z segmts
  G4int    m_theta;
  // total energy on sensitive material
  G4float  m_edep;
  // total energy deposited on absorber
  G4float  m_edepabs;
  G4ThreeVector m_pos;
  // number of geant particles contributing to this hit
  G4int    m_nmc;
  // vector of particles' Ids contributing to visible energy
  std::vector<G4int>    *m_trkid;
  // vector of energy contributions coming from each particle
  std::vector<G4double> *m_emc;
  // vector of timing of energy contributions coming from each particle
  std::vector<G4double> *m_tmc;
};

typedef G4THitsCollection<LCDG4CalHit> LCDG4CalHitsCollection;

extern G4Allocator<LCDG4CalHit> LCDG4CalHitAllocator;

inline void* LCDG4CalHit::operator new(size_t) {
  void *aHit;
  aHit = (void *) LCDG4CalHitAllocator.MallocSingle();
  return aHit;
}

inline void LCDG4CalHit::operator delete(void *aHit) {
  LCDG4CalHitAllocator.FreeSingle((LCDG4CalHit*) aHit);
}

#endif

---