File indexing completed on 2024-04-06 12:30:13
0001 #ifndef SimG4Core_Phase2SteppingAction_H
0002 #define SimG4Core_Phase2SteppingAction_H
0003
0004 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0005 #include "SimG4Core/Notification/interface/SimActivityRegistry.h"
0006 #include "SimG4Core/Application/interface/SteppingTrackStatus.h"
0007
0008 #include "G4LogicalVolume.hh"
0009 #include "G4Region.hh"
0010 #include "G4UserSteppingAction.hh"
0011 #include "G4VPhysicalVolume.hh"
0012 #include "G4VTouchable.hh"
0013 #include "G4Track.hh"
0014
0015 #include <string>
0016 #include <vector>
0017
0018 class CMSSteppingVerbose;
0019
0020 class Phase2SteppingAction : public G4UserSteppingAction {
0021 public:
0022 explicit Phase2SteppingAction(const CMSSteppingVerbose*, const edm::ParameterSet&, bool, bool);
0023 ~Phase2SteppingAction() override = default;
0024
0025 void UserSteppingAction(const G4Step* aStep) final;
0026
0027 SimActivityRegistry::G4StepSignal m_g4StepSignal;
0028
0029 private:
0030 bool initPointer();
0031
0032 inline bool isInsideDeadRegion(const G4Region* reg) const;
0033 inline bool isOutOfTimeWindow(const G4Region* reg, const double& time) const;
0034
0035 bool isLowEnergy(const G4LogicalVolume*, const G4Track*) const;
0036 void PrintKilledTrack(const G4Track*, const TrackStatus&) const;
0037
0038 const G4VPhysicalVolume* cmse{nullptr};
0039 const G4VPhysicalVolume* tracker{nullptr};
0040 const G4VPhysicalVolume* calo{nullptr};
0041 const G4VPhysicalVolume* btl{nullptr};
0042 const CMSSteppingVerbose* steppingVerbose;
0043 double theCriticalEnergyForVacuum;
0044 double theCriticalDensity;
0045 double maxTrackTime;
0046 double maxTrackTimeForward;
0047 double maxZCentralCMS;
0048
0049 unsigned int numberTimes;
0050 unsigned int numberEkins;
0051 unsigned int numberPart;
0052 unsigned int ndeadRegions;
0053 unsigned int nWarnings{0};
0054 G4int maxNumberOfSteps;
0055
0056 bool initialized{false};
0057 bool killBeamPipe{false};
0058 bool hasWatcher;
0059 bool dd4hep_;
0060
0061 std::vector<double> maxTrackTimes, ekinMins;
0062 std::vector<std::string> maxTimeNames, ekinNames, ekinParticles;
0063 std::vector<std::string> deadRegionNames;
0064 std::vector<const G4Region*> maxTimeRegions;
0065 std::vector<const G4Region*> deadRegions;
0066 std::vector<G4LogicalVolume*> ekinVolumes;
0067 std::vector<int> ekinPDG;
0068 G4String cmseName_, trackerName_, caloName_, btlName_, cms2ZDCName_;
0069 };
0070
0071 inline bool Phase2SteppingAction::isInsideDeadRegion(const G4Region* reg) const {
0072 bool res = false;
0073 for (auto& region : deadRegions) {
0074 if (reg == region) {
0075 res = true;
0076 break;
0077 }
0078 }
0079 return res;
0080 }
0081
0082 inline bool Phase2SteppingAction::isOutOfTimeWindow(const G4Region* reg, const double& time) const {
0083 double tofM = maxTrackTime;
0084 for (unsigned int i = 0; i < numberTimes; ++i) {
0085 if (reg == maxTimeRegions[i]) {
0086 tofM = maxTrackTimes[i];
0087 break;
0088 }
0089 }
0090 return (time > tofM);
0091 }
0092
0093 #endif
