TotemRPSD.h

Go to the documentation of this file.

#ifndef PPS_TotemRPSD_h
#define PPS_TotemRPSD_h

#include "SimG4CMS/PPS/interface/TotemRPG4Hit.h"
#include "SimG4CMS/PPS/interface/TotemRPG4HitCollection.h"
#include "SimG4CMS/PPS/interface/TotemRPVDetectorOrganization.h"

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

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/ParameterSet/interface/ParameterSetfwd.h"

#include <string>

class G4Step;
class G4HCofThisEvent;
class G4Track;
class G4StepPoint;
class TrackingSlaveSD;
class SimTrackManager;
class TotemTestHitHBNtuple;

class TotemRPSD : public SensitiveTkDetector, public Observer<const BeginOfEvent*>, public Observer<const EndOfEvent*> {
public:
  TotemRPSD(const std::string&, const SensitiveDetectorCatalog&, edm::ParameterSet const&, const SimTrackManager*);
  ~TotemRPSD() override;

  // Geant4 methods
  bool ProcessHits(G4Step* step, G4TouchableHistory* tHistory) override;

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

  // CMSSW methods
  uint32_t setDetUnitId(const G4Step* step) override;
  void fillHits(edm::PSimHitContainer&, const std::string&) override;
  void clearHits() override;

protected:
  void update(const BeginOfEvent*) override;
  void update(const ::EndOfEvent*) override;

private:
  G4ThreeVector setToLocal(const G4ThreeVector& globalPoint);
  void stepInfo(const G4Step* aStep);
  void createNewHit();
  void storeHit(TotemRPG4Hit*);
  void printHitInfo();

  std::unique_ptr<TrackingSlaveSD> slave_;
  std::unique_ptr<TotemRPVDetectorOrganization> numberingScheme_;

  TotemRPG4HitCollection* theHC_ = nullptr;
  TotemRPG4Hit* currentHit_ = nullptr;
  G4Track* theTrack_ = nullptr;
  G4VPhysicalVolume* currentPV_ = nullptr;
  G4int hcID_ = -1;
  G4int primaryID_ = 0;
  G4int parentID_ = 0;
  G4int tSliceID_ = 0;
  G4double tSlice_ = 0.0;

  G4StepPoint* preStepPoint_ = nullptr;
  G4StepPoint* postStepPoint_ = nullptr;
  G4ThreeVector hitPoint_;
  G4ThreeVector exitPoint_;
  G4ThreeVector theLocalEntryPoint_;
  G4ThreeVector theLocalExitPoint_;

  double incidentEnergy_ = 0.0;
  double pabs_ = 0.0;
  double thePx_ = 0.0;
  double thePy_ = 0.0;
  double thePz_ = 0.0;
  double tof_ = 0.0;
  double eloss_ = 0.0;

  double thetaAtEntry_ = 0.0;
  double phiAtEntry_ = 0.0;

  double vx_ = 0.0;
  double vy_ = 0.0;
  double vz_ = 0.0;

  unsigned int unitID_ = 0;
  int verbosity_;
  int eventno_ = 0;
  short particleType_ = 0;
};

#endif  // PPS_TotemRPSD_h