LCDG4CalSD.hh

Go to the documentation of this file.

#ifndef LCDG4CalSD_h
#define LCDG4CalSD_h 1

#include "LCDG4CalHit.hh"
#include "G4VSensitiveDetector.hh"
#include "G4LogicalVolume.hh"
#include <map>

class G4Step;
class G4HCofThisEvent;
class G4TouchableHistory;

/* SensitiveDetector base class for calorimeter-like subdetectors */
// Some derived classes are LCDG4CalSDProj and LCDG4CalSDNonProj
class LCDG4CalSD : public G4VSensitiveDetector {

public:
  // constructor
  LCDG4CalSD(G4String name);
  // destructor
  ~LCDG4CalSD();

  // For each geant4 step, find which Cal cell was hit.
  // If that cell is not hit yet, create it.
  // Otherwise, just add up the deposited energy into that cell
  G4bool ProcessHits(G4Step* aStep,G4TouchableHistory* ROhist);

  // Begin of event processing
  void Initialize(G4HCofThisEvent* HCE);

  // end of event processing
  void EndOfEvent(G4HCofThisEvent* HCE);

  //*** Pure virtual methods
  // find calorimeter cell from a given position
  virtual unsigned int findCell(const G4ThreeVector& position) = 0;
  // get cell center
  virtual G4ThreeVector getCellCenter(unsigned int cellIndex) = 0;
  // get cell indices which depend on virtual cells
  virtual unsigned int getIndex1(unsigned int cellIndex) const = 0;
  virtual unsigned int getIndex2(unsigned int cellIndex) const = 0;
  virtual unsigned int getIndex3(unsigned int cellIndex) const = 0;

  // some convenience methods
  void clear();
  void DrawAll();
  void PrintAll();

public:
  G4int AddLogVolInfo(G4LogicalVolume* g4vxdlog, G4int startlayerNo,
                      G4int sysNo, G4int be, G4int ns, G4int fsens);

  // These are needed to calculate cell centers in projective geometry
  void SetLayerZeroCenterDepth(G4double cd) { layer0CenterDepth = cd; }
  void SetLayerThickness(G4double t) { layerThickness = t; }

protected:
  // hit collection ID associated to this sensitive detector
  G4int HCID;
  // pointer to calorimeter hits collection
  LCDG4CalHitsCollection* _hitCollection;
  // Cal layer map
  std::map<G4LogicalVolume*,G4int> _maplayer;
  // Cal sysNo map
  std::map<G4LogicalVolume*,G4int> _mapsys;
  // Cal barrel/endcap map
  std::map<G4LogicalVolume*,G4int> _mapbe;
  // Cal n/s map: -1,0,1 for north endcap, barrel, south endcap
  std::map<G4LogicalVolume*,G4int> _mapns;
  // Cal sensitive/absorber map
  std::map<G4LogicalVolume*,G4int> _mapsens;

  // Calorimeter hit collections handled by this SD
  std::map<G4int,LCDG4CalHit*> _cellhit;

  // These are needed to calculate radius for centers of sensitive layers
  // This is only used by projective SD...
  G4double layer0CenterDepth;
  G4double layerThickness;
  G4int north, barend;

  // debug flag
  int debugLevel;
  // convenience info: cell indices
  // index[0] is system number (2=LUM, 10=EM, 12=HAD, 7=MU)
  // index[1] is ilayer
  // index[2] is iphi
  // index[3] is itheta (proj) or iz (nonproj)
  unsigned int index[4];
};

#endif

---