FP420SD.h

Go to the documentation of this file.

//
#ifndef FP420_FP420SD_h
#define FP420_FP420SD_h
//

#include "SimG4Core/Notification/interface/Observer.h"
#include "SimG4Core/SensitiveDetector/interface/SensitiveTkDetector.h"

#include "SimG4Core/Notification/interface/BeginOfRun.h"
#include "SimG4Core/Notification/interface/BeginOfEvent.h"
#include "SimG4Core/Notification/interface/EndOfEvent.h"

#include "SimG4Core/Notification/interface/BeginOfTrack.h"
#include "SimG4Core/Notification/interface/BeginOfJob.h"

#include "DataFormats/GeometryVector/interface/LocalPoint.h"

#include "SimG4CMS/FP420/interface/FP420G4Hit.h"
#include "SimG4CMS/FP420/interface/FP420G4HitCollection.h"
#include "SimG4CMS/FP420/interface/FP420NumberingScheme.h"

#include "G4Step.hh"
#include "G4StepPoint.hh"
#include "G4Track.hh"
#include "G4VPhysicalVolume.hh"

#include <string>

class TrackingSlaveSD;
//AZ:
class FrameRotation;
class TrackInformation;
class SimTrackManager;
class TrackingSlaveSD;
class UpdatablePSimHit;
class G4ProcessTypeEnumerator;
class G4TrackToParticleID;

//-------------------------------------------------------------------

class FP420SD : public SensitiveTkDetector,
                public Observer<const BeginOfRun*>,
                public Observer<const BeginOfEvent*>,
                public Observer<const EndOfEvent*> {
public:
  FP420SD(const std::string&,
          const DDCompactView&,
          const SensitiveDetectorCatalog&,
          edm::ParameterSet const&,
          const SimTrackManager*);

  ~FP420SD() override;

  bool ProcessHits(G4Step*, G4TouchableHistory*) override;
  uint32_t setDetUnitId(const G4Step*) override;

  void Initialize(G4HCofThisEvent* HCE) override;
  void EndOfEvent(G4HCofThisEvent* eventHC) override;
  void clear() override;
  void DrawAll() override;
  void PrintAll() override;

  virtual double getEnergyDeposit(G4Step* step);
  void fillHits(edm::PSimHitContainer&, const std::string&) override;
  void clearHits() override;

private:
  void update(const BeginOfRun*) override;
  void update(const BeginOfEvent*) override;
  void update(const ::EndOfEvent*) override;

  G4ThreeVector SetToLocal(const G4ThreeVector& global);
  G4ThreeVector SetToLocalExit(const G4ThreeVector& globalPoint);
  void GetStepInfo(G4Step* aStep);
  G4bool HitExists();
  void CreateNewHit();
  void UpdateHit();
  void StoreHit(FP420G4Hit*);
  void ResetForNewPrimary();
  void Summarize();

private:
  //AZ:
  TrackingSlaveSD* slave;
  FP420NumberingScheme* numberingScheme;

  G4ThreeVector entrancePoint, exitPoint;
  G4ThreeVector theEntryPoint;
  G4ThreeVector theExitPoint;

  float incidentEnergy;

  G4int hcID;
  FP420G4HitCollection* theHC;
  const SimTrackManager* theManager;

  G4int tsID;
  FP420G4Hit* currentHit;
  G4Track* theTrack;
  G4VPhysicalVolume* currentPV;
  // unsigned int         unitID, previousUnitID;
  uint32_t unitID, previousUnitID;
  G4int tSliceID;
  unsigned int primaryID, primID;

  G4double tSlice;

  G4StepPoint* preStepPoint;
  G4StepPoint* postStepPoint;
  float edeposit;

  G4ThreeVector hitPoint;
  G4ThreeVector hitPointExit;
  G4ThreeVector hitPointLocal;
  G4ThreeVector hitPointLocalExit;
  float Pabs;
  float Tof;
  float Eloss;
  short ParticleType;

  float ThetaAtEntry;
  float PhiAtEntry;

  int ParentId;
  float Vx, Vy, Vz;
  float X, Y, Z;

  //
  // Hist
  //
  int eventno;

protected:
  float edepositEM, edepositHAD;
  G4int emPDG;
  G4int epPDG;
  G4int gammaPDG;
};

#endif  // FP420SD_h