From 79db293172963e79d57b35eebab38aa5fe19797a Mon Sep 17 00:00:00 2001 From: Vladimir Date: Mon, 10 Jul 2023 15:39:01 +0200 Subject: [PATCH 1/3] SIM scripts and MC truth clean-up --- .../HcalTestBeam/python/TB2006Analysis_cfi.py | 2 +- .../test/python/TB2006Analysis_cfi.py | 74 +++++++++++++++++++ SimG4Core/Application/python/g4SimHits_cfi.py | 2 +- .../Notification/interface/TmpSimEvent.h | 2 +- .../Notification/interface/TmpSimVertex.h | 8 +- .../Notification/interface/TrackWithHistory.h | 53 +++++++------ SimG4Core/Notification/src/SimTrackManager.cc | 7 +- SimG4Core/Notification/src/TmpSimEvent.cc | 4 +- .../Notification/src/TrackWithHistory.cc | 24 ++---- 9 files changed, 115 insertions(+), 61 deletions(-) create mode 100644 SimG4CMS/HcalTestBeam/test/python/TB2006Analysis_cfi.py diff --git a/SimG4CMS/HcalTestBeam/python/TB2006Analysis_cfi.py b/SimG4CMS/HcalTestBeam/python/TB2006Analysis_cfi.py index 715ae1b90518c..1cabd740603e5 100644 --- a/SimG4CMS/HcalTestBeam/python/TB2006Analysis_cfi.py +++ b/SimG4CMS/HcalTestBeam/python/TB2006Analysis_cfi.py @@ -40,7 +40,7 @@ def testbeam2006(process): process.common_maximum_time.MaxTrackTime = cms.double(1000.0) process.common_maximum_time.MaxTimeNames = cms.vstring() - process.common_maximum_time.MaxTrackTimes = cms.vstring() + process.common_maximum_time.MaxTrackTimes = cms.vdouble() process.common_maximum_time.DeadRegions = cms.vstring() process.g4SimHits.NonBeamEvent = True diff --git a/SimG4CMS/HcalTestBeam/test/python/TB2006Analysis_cfi.py b/SimG4CMS/HcalTestBeam/test/python/TB2006Analysis_cfi.py new file mode 100644 index 0000000000000..524e8c56a9651 --- /dev/null +++ b/SimG4CMS/HcalTestBeam/test/python/TB2006Analysis_cfi.py @@ -0,0 +1,74 @@ +import FWCore.ParameterSet.Config as cms + +from IOMC.EventVertexGenerators.VtxSmearedParameters_cfi import * + +def testbeam2006(process): + + process.load('FWCore.MessageService.MessageLogger_cfi') + process.load('Configuration.StandardSequences.Services_cff') + process.load('SimGeneral.HepPDTESSource.pdt_cfi') + process.load('Configuration.EventContent.EventContent_cff') + process.load("Geometry.EcalCommonData.ecalSimulationParameters_cff") + process.load('Geometry.HcalTestBeamData.hcalDDDSimConstants_cff') + process.load('IOMC.EventVertexGenerators.VtxSmearedFlat_cfi') + process.load('GeneratorInterface.Core.generatorSmeared_cfi') + process.load('SimG4Core.Application.g4SimHits_cfi') + process.load('IOMC.RandomEngine.IOMC_cff') + + process.RandomNumberGeneratorService.generator.initialSeed = 456789 + process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876 + process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789 + + process.source = cms.Source("EmptySource", + firstRun = cms.untracked.uint32(1), + firstEvent = cms.untracked.uint32(1) + ) + + process.load('SimG4CMS.HcalTestBeam.TBDirectionParameters_cfi') + process.load('SimG4CMS.HcalTestBeam.TBVtxSmeared_cfi') + process.load('SimG4CMS.HcalTestBeam.TB06Analysis_cfi') + + process.generator = cms.EDProducer("FlatRandomEGunProducer", + PGunParameters = cms.PSet( + process.common_beam_direction_parameters + ), + Verbosity = cms.untracked.int32(0), + AddAntiParticle = cms.bool(False) + ) + + process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.g4SimHits*process.testbeam) + + process.common_maximum_time.MaxTrackTime = cms.double(1000.0) +# process.common_maximum_time.MaxTimeNames = cms.vstring() +# process.common_maximum_time.MaxTrackTimes = cms.vdouble() +# process.common_maximum_time.DeadRegions = cms.vstring() + + process.g4SimHits.NonBeamEvent = True + process.g4SimHits.UseMagneticField = False + + process.g4SimHits.Physics.type = 'SimG4Core/Physics/FTFP_BERT_EMM' + process.g4SimHits.Physics.Region = 'HcalRegion' + process.g4SimHits.Physics.MaxTrackTime = cms.double(1000.0) + + process.g4SimHits.Generator.ApplyEtaCuts = cms.bool(False) + + process.g4SimHits.StackingAction.MaxTrackTime = cms.double(1000.0) + process.g4SimHits.StackingAction.MaxTimeNames = cms.vstring() + process.g4SimHits.StackingAction.MaxTrackTimes = cms.vdouble() + process.g4SimHits.StackingAction.DeadRegions = cms.vstring() + + process.g4SimHits.SteppingAction.MaxTrackTime = cms.double(1000.0) + process.g4SimHits.SteppingAction.MaxTimeNames = cms.vstring() + process.g4SimHits.SteppingAction.MaxTrackTimes = cms.vdouble() + process.g4SimHits.SteppingAction.DeadRegions = cms.vstring() + + process.g4SimHits.CaloSD.EminHits = cms.vdouble(0.0,0.0,0.0,0.0) + process.g4SimHits.CaloSD.TmaxHits = cms.vdouble(1000.0,1000.0,1000.0,1000.0) + process.g4SimHits.CaloTrkProcessing.TestBeam = True + + process.g4SimHits.HCalSD.UseShowerLibrary = False + process.g4SimHits.HCalSD.UseHF = False + process.g4SimHits.HCalSD.ForTBHCAL = True + process.g4SimHits.HCalSD.ForTBH2 = False + + return(process) diff --git a/SimG4Core/Application/python/g4SimHits_cfi.py b/SimG4Core/Application/python/g4SimHits_cfi.py index 63079a0bff09d..3939a7a0ade21 100644 --- a/SimG4Core/Application/python/g4SimHits_cfi.py +++ b/SimG4Core/Application/python/g4SimHits_cfi.py @@ -706,7 +706,7 @@ from Configuration.Eras.Modifier_phase2_common_cff import phase2_common phase2_common.toModify(g4SimHits, - OnlySDs = ['ZdcSensitiveDetector', 'TotemT2ScintSensitiveDetector', 'TotemSensitiveDetector', 'RomanPotSensitiveDetector', 'PLTSensitiveDetector', 'MuonSensitiveDetector', 'MtdSensitiveDetector', 'BCM1FSensitiveDetector', 'EcalSensitiveDetector', 'CTPPSSensitiveDetector', 'HGCalSensitiveDetector', 'BSCSensitiveDetector', 'CTPPSDiamondSensitiveDetector', 'FP420SensitiveDetector', 'BHMSensitiveDetector', 'HFNoseSensitiveDetector', 'HGCScintillatorSensitiveDetector', 'CastorSensitiveDetector', 'CaloTrkProcessing', 'HcalSensitiveDetector', 'TkAccumulatingSensitiveDetector'], + OnlySDs = ['ZdcSensitiveDetector', 'TotemT2ScintSensitiveDetector', 'TotemSensitiveDetector', 'RomanPotSensitiveDetector', 'PLTSensitiveDetector', 'MuonSensitiveDetector', 'MtdSensitiveDetector', 'BCM1FSensitiveDetector', 'EcalSensitiveDetector', 'CTPPSSensitiveDetector', 'HGCalSensitiveDetector', 'CTPPSDiamondSensitiveDetector', 'FP420SensitiveDetector', 'BHMSensitiveDetector', 'HFNoseSensitiveDetector', 'HGCScintillatorSensitiveDetector', 'CastorSensitiveDetector', 'CaloTrkProcessing', 'HcalSensitiveDetector', 'TkAccumulatingSensitiveDetector'], LHCTransport = False, MuonSD = dict( HaveDemoChambers = False ) diff --git a/SimG4Core/Notification/interface/TmpSimEvent.h b/SimG4Core/Notification/interface/TmpSimEvent.h index 2bd883e38452a..28e3e35538344 100644 --- a/SimG4Core/Notification/interface/TmpSimEvent.h +++ b/SimG4Core/Notification/interface/TmpSimEvent.h @@ -14,7 +14,7 @@ class TmpSimEvent { public: TmpSimEvent(); - virtual ~TmpSimEvent(); + ~TmpSimEvent(); void load(edm::SimTrackContainer& c) const; void load(edm::SimVertexContainer& c) const; unsigned int nTracks() const { return g4tracks_.size(); } diff --git a/SimG4Core/Notification/interface/TmpSimVertex.h b/SimG4Core/Notification/interface/TmpSimVertex.h index 613b8e52e4b55..d66417f423f72 100644 --- a/SimG4Core/Notification/interface/TmpSimVertex.h +++ b/SimG4Core/Notification/interface/TmpSimVertex.h @@ -7,20 +7,20 @@ class TmpSimVertex { public: - TmpSimVertex(const math::XYZVectorD& ip, double it, int iv, unsigned int typ = 0) - : ilv_(ip), itime_(it), itrack_(iv), procType_(typ) {} + TmpSimVertex(const math::XYZVectorD& ip, double it, int iv, int typ=0) + : ilv_(ip), itime_(it), itrack_(iv), ptype_(typ) {} ~TmpSimVertex() = default; /// index of the parent (-1 if no parent) const math::XYZVectorD& vertexPosition() const { return ilv_; } double vertexGlobalTime() const { return itime_; } int parentIndex() const { return itrack_; } - unsigned int processType() const { return procType_; } + int processType() const { return ptype_; } private: math::XYZVectorD ilv_; double itime_; int itrack_; - unsigned int procType_; + int ptype_; }; #endif diff --git a/SimG4Core/Notification/interface/TrackWithHistory.h b/SimG4Core/Notification/interface/TrackWithHistory.h index 80de4cd633981..1eb24570c0c63 100644 --- a/SimG4Core/Notification/interface/TrackWithHistory.h +++ b/SimG4Core/Notification/interface/TrackWithHistory.h @@ -7,7 +7,6 @@ #include "G4Allocator.hh" -class G4VProcess; class G4PrimaryParticle; /** The part of the information about a SimTrack that we need from * a G4Track @@ -18,33 +17,36 @@ class TrackWithHistory { /** The constructor is called at time, * when some of the information may not available yet. */ - TrackWithHistory(const G4Track *g4track, int pID); - TrackWithHistory(const G4PrimaryParticle *, int trackID, const math::XYZVectorD &pos, double time); + TrackWithHistory(const G4Track* g4track, int pID); + TrackWithHistory(const G4PrimaryParticle*, int trackID, const math::XYZVectorD &pos, const double time); ~TrackWithHistory() = default; - inline void *operator new(size_t); + inline void *operator new(std::size_t); inline void operator delete(void *TrackWithHistory); - void setToBeSaved() { saved_ = true; } int trackID() const { return trackID_; } - int particleID() const { return particleID_; } + int particleID() const { return pdgID_; } int parentID() const { return parentID_; } int genParticleID() const { return genParticleID_; } int vertexID() const { return vertexID_; } - const math::XYZVectorD &momentum() const { return momentum_; } - double totalEnergy() const { return totalEnergy_; } - const math::XYZVectorD &vertexPosition() const { return vertexPosition_; } - double globalTime() const { return globalTime_; } - double localTime() const { return localTime_; } - double properTime() const { return properTime_; } - const G4VProcess *creatorProcess() const { return creatorProcess_; } - double weight() const { return weight_; } + int processType() const { return procType_; } + int getIDAtBoundary() const { return idAtBoundary_; } + void setTrackID(int i) { trackID_ = i; } void setParentID(int i) { parentID_ = i; } void setVertexID(int i) { vertexID_ = i; } void setGenParticleID(int i) { genParticleID_ = i; } + + double totalEnergy() const { return totalEnergy_; } + double time() const { return time_; } + double weight() const { return weight_; } + void setToBeSaved() { saved_ = true; } bool storeTrack() const { return storeTrack_; } bool saved() const { return saved_; } + bool crossedBoundary() const { return crossedBoundary_; } + + const math::XYZVectorD &momentum() const { return momentum_; } + const math::XYZVectorD &vertexPosition() const { return vertexPosition_; } // Boundary crossing variables void setCrossedBoundaryPosMom(int id, @@ -55,10 +57,8 @@ class TrackWithHistory { positionAtBoundary_ = position; momentumAtBoundary_ = momentum; } - bool crossedBoundary() const { return crossedBoundary_; } const math::XYZTLorentzVectorF &getPositionAtBoundary() const { return positionAtBoundary_; } const math::XYZTLorentzVectorF &getMomentumAtBoundary() const { return momentumAtBoundary_; } - int getIDAtBoundary() const { return idAtBoundary_; } // tracker surface const math::XYZVectorD &trackerSurfacePosition() const { return tkSurfacePosition_; } @@ -70,27 +70,24 @@ class TrackWithHistory { private: int trackID_; - int particleID_; + int pdgID_; int parentID_; int genParticleID_{-1}; int vertexID_{-1}; - math::XYZVectorD momentum_; + int idAtBoundary_{-1}; + int procType_{0}; double totalEnergy_; - math::XYZVectorD vertexPosition_; - double globalTime_; - double localTime_; - double properTime_; - const G4VProcess *creatorProcess_; + double time_; // lab system double weight_; - bool storeTrack_; - bool saved_{false}; - - bool crossedBoundary_{false}; - int idAtBoundary_{-1}; + math::XYZVectorD momentum_; + math::XYZVectorD vertexPosition_; math::XYZTLorentzVectorF positionAtBoundary_{math::XYZTLorentzVectorF(0.f, 0.f, 0.f, 0.f)}; math::XYZTLorentzVectorF momentumAtBoundary_{math::XYZTLorentzVectorF(0.f, 0.f, 0.f, 0.f)}; math::XYZVectorD tkSurfacePosition_{math::XYZVectorD(0., 0., 0.)}; math::XYZTLorentzVectorD tkSurfaceMomentum_{math::XYZTLorentzVectorD(0., 0., 0., 0.)}; + bool storeTrack_{false}; + bool saved_{false}; + bool crossedBoundary_{false}; }; extern G4ThreadLocal G4Allocator *fpTrackWithHistoryAllocator; diff --git a/SimG4Core/Notification/src/SimTrackManager.cc b/SimG4Core/Notification/src/SimTrackManager.cc index 2979ea088293a..ed5a6818f5eea 100644 --- a/SimG4Core/Notification/src/SimTrackManager.cc +++ b/SimG4Core/Notification/src/SimTrackManager.cc @@ -189,12 +189,7 @@ int SimTrackManager::getOrCreateVertex(TrackWithHistory* trkH, int iParentID) { } } - unsigned int ptype = 0; - const G4VProcess* pr = trkH->creatorProcess(); - if (nullptr != pr) { - ptype = pr->GetProcessSubType(); - } - m_simEvent->add(new TmpSimVertex(trkH->vertexPosition(), trkH->globalTime(), parent, ptype)); + m_simEvent->add(new TmpSimVertex(trkH->vertexPosition(), trkH->time(), parent, trkH->processType())); m_vertexMap[parent].push_back(VertexPosition(m_nVertices, trkH->vertexPosition())); ++m_nVertices; return (m_nVertices - 1); diff --git a/SimG4Core/Notification/src/TmpSimEvent.cc b/SimG4Core/Notification/src/TmpSimEvent.cc index 673407c96de55..f05dab1437945 100644 --- a/SimG4Core/Notification/src/TmpSimEvent.cc +++ b/SimG4Core/Notification/src/TmpSimEvent.cc @@ -56,13 +56,13 @@ void TmpSimEvent::load(edm::SimVertexContainer& c) const { TmpSimVertex* vtx = g4vertices_[i]; auto pos = vtx->vertexPosition(); math::XYZVectorD v3(pos.x() * invcm, pos.y() * invcm, pos.z() * invcm); - float t = vtx->vertexGlobalTime() / CLHEP::second; + float t = (float)(vtx->vertexGlobalTime() / CLHEP::second); int iv = vtx->parentIndex(); // v3 = position in cm // t = global time in second // iv = index of the parent in the SimEvent SimTrack container (-1 if no parent) SimVertex v = SimVertex(v3, t, iv, i); - v.setProcessType(vtx->processType()); + v.setProcessType((unsigned int)vtx->processType()); v.setEventId(EncodedEventId(0)); c.push_back(v); } diff --git a/SimG4Core/Notification/src/TrackWithHistory.cc b/SimG4Core/Notification/src/TrackWithHistory.cc index 9aa4c62de3062..ea6a843d19d51 100644 --- a/SimG4Core/Notification/src/TrackWithHistory.cc +++ b/SimG4Core/Notification/src/TrackWithHistory.cc @@ -14,21 +14,18 @@ G4ThreadLocal G4Allocator* fpTrackWithHistoryAllocator = nullptr; -//#define DEBUG - TrackWithHistory::TrackWithHistory(const G4Track* g4trk, int pID) { trackID_ = g4trk->GetTrackID(); - particleID_ = G4TrackToParticleID::particleID(g4trk); + pdgID_ = G4TrackToParticleID::particleID(g4trk); parentID_ = pID; auto mom = g4trk->GetMomentum(); momentum_ = math::XYZVectorD(mom.x(), mom.y(), mom.z()); totalEnergy_ = g4trk->GetTotalEnergy(); auto pos = g4trk->GetPosition(); vertexPosition_ = math::XYZVectorD(pos.x(), pos.y(), pos.z()); - globalTime_ = g4trk->GetGlobalTime(); - localTime_ = g4trk->GetLocalTime(); - properTime_ = g4trk->GetProperTime(); - creatorProcess_ = g4trk->GetCreatorProcess(); + time_ = g4trk->GetGlobalTime(); + auto p = g4trk->GetCreatorProcess(); + procType_ = (nullptr != p) ? p->GetProcessSubType() : 0; TrackInformation* trkinfo = static_cast(g4trk->GetUserInformation()); storeTrack_ = trkinfo->storeTrack(); auto vgprimary = g4trk->GetDynamicParticle()->GetPrimaryParticle(); @@ -41,30 +38,21 @@ TrackWithHistory::TrackWithHistory(const G4Track* g4trk, int pID) { // V.I. weight is computed in the same way as before // without usage of G4Track::GetWeight() weight_ = 10000 * genParticleID_; -#ifdef DEBUG - LogDebug("TrackInformation") << " TrackWithHistory : created history for " << trackID_ << " with mother " - << parentID_; -#endif } TrackWithHistory::TrackWithHistory(const G4PrimaryParticle* ptr, int trackID, const math::XYZVectorD& pos, double time) { trackID_ = trackID; - particleID_ = G4TrackToParticleID::particleID(ptr, trackID_); + pdgID_ = G4TrackToParticleID::particleID(ptr, trackID_); parentID_ = trackID; auto mom = ptr->GetMomentum(); momentum_ = math::XYZVectorD(mom.x(), mom.y(), mom.z()); totalEnergy_ = ptr->GetTotalEnergy(); vertexPosition_ = math::XYZVectorD(pos.x(), pos.y(), pos.z()); - globalTime_ = localTime_ = properTime_ = time; - creatorProcess_ = nullptr; + time_ = time; storeTrack_ = true; auto priminfo = static_cast(ptr->GetUserInformation()); if (nullptr != priminfo) { genParticleID_ = priminfo->id(); } weight_ = 10000. * genParticleID_; -#ifdef DEBUG - LogDebug("TrackInformation") << " TrackWithHistory : created history for " << trackID_ << " with mother " - << parentID_; -#endif } From 6077f4d6035e47edb253a1f2ed884d9ca81d4edd Mon Sep 17 00:00:00 2001 From: Vladimir Date: Mon, 10 Jul 2023 15:41:02 +0200 Subject: [PATCH 2/3] removed unused --- .../test/python/TB2006Analysis_cfi.py | 74 ------------------- 1 file changed, 74 deletions(-) delete mode 100644 SimG4CMS/HcalTestBeam/test/python/TB2006Analysis_cfi.py diff --git a/SimG4CMS/HcalTestBeam/test/python/TB2006Analysis_cfi.py b/SimG4CMS/HcalTestBeam/test/python/TB2006Analysis_cfi.py deleted file mode 100644 index 524e8c56a9651..0000000000000 --- a/SimG4CMS/HcalTestBeam/test/python/TB2006Analysis_cfi.py +++ /dev/null @@ -1,74 +0,0 @@ -import FWCore.ParameterSet.Config as cms - -from IOMC.EventVertexGenerators.VtxSmearedParameters_cfi import * - -def testbeam2006(process): - - process.load('FWCore.MessageService.MessageLogger_cfi') - process.load('Configuration.StandardSequences.Services_cff') - process.load('SimGeneral.HepPDTESSource.pdt_cfi') - process.load('Configuration.EventContent.EventContent_cff') - process.load("Geometry.EcalCommonData.ecalSimulationParameters_cff") - process.load('Geometry.HcalTestBeamData.hcalDDDSimConstants_cff') - process.load('IOMC.EventVertexGenerators.VtxSmearedFlat_cfi') - process.load('GeneratorInterface.Core.generatorSmeared_cfi') - process.load('SimG4Core.Application.g4SimHits_cfi') - process.load('IOMC.RandomEngine.IOMC_cff') - - process.RandomNumberGeneratorService.generator.initialSeed = 456789 - process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876 - process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789 - - process.source = cms.Source("EmptySource", - firstRun = cms.untracked.uint32(1), - firstEvent = cms.untracked.uint32(1) - ) - - process.load('SimG4CMS.HcalTestBeam.TBDirectionParameters_cfi') - process.load('SimG4CMS.HcalTestBeam.TBVtxSmeared_cfi') - process.load('SimG4CMS.HcalTestBeam.TB06Analysis_cfi') - - process.generator = cms.EDProducer("FlatRandomEGunProducer", - PGunParameters = cms.PSet( - process.common_beam_direction_parameters - ), - Verbosity = cms.untracked.int32(0), - AddAntiParticle = cms.bool(False) - ) - - process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.g4SimHits*process.testbeam) - - process.common_maximum_time.MaxTrackTime = cms.double(1000.0) -# process.common_maximum_time.MaxTimeNames = cms.vstring() -# process.common_maximum_time.MaxTrackTimes = cms.vdouble() -# process.common_maximum_time.DeadRegions = cms.vstring() - - process.g4SimHits.NonBeamEvent = True - process.g4SimHits.UseMagneticField = False - - process.g4SimHits.Physics.type = 'SimG4Core/Physics/FTFP_BERT_EMM' - process.g4SimHits.Physics.Region = 'HcalRegion' - process.g4SimHits.Physics.MaxTrackTime = cms.double(1000.0) - - process.g4SimHits.Generator.ApplyEtaCuts = cms.bool(False) - - process.g4SimHits.StackingAction.MaxTrackTime = cms.double(1000.0) - process.g4SimHits.StackingAction.MaxTimeNames = cms.vstring() - process.g4SimHits.StackingAction.MaxTrackTimes = cms.vdouble() - process.g4SimHits.StackingAction.DeadRegions = cms.vstring() - - process.g4SimHits.SteppingAction.MaxTrackTime = cms.double(1000.0) - process.g4SimHits.SteppingAction.MaxTimeNames = cms.vstring() - process.g4SimHits.SteppingAction.MaxTrackTimes = cms.vdouble() - process.g4SimHits.SteppingAction.DeadRegions = cms.vstring() - - process.g4SimHits.CaloSD.EminHits = cms.vdouble(0.0,0.0,0.0,0.0) - process.g4SimHits.CaloSD.TmaxHits = cms.vdouble(1000.0,1000.0,1000.0,1000.0) - process.g4SimHits.CaloTrkProcessing.TestBeam = True - - process.g4SimHits.HCalSD.UseShowerLibrary = False - process.g4SimHits.HCalSD.UseHF = False - process.g4SimHits.HCalSD.ForTBHCAL = True - process.g4SimHits.HCalSD.ForTBH2 = False - - return(process) From 591a0f64482980ca45ddf5897b604208873cbed8 Mon Sep 17 00:00:00 2001 From: Vladimir Date: Mon, 10 Jul 2023 15:43:26 +0200 Subject: [PATCH 3/3] code format --- SimG4Core/Notification/interface/TmpSimVertex.h | 2 +- SimG4Core/Notification/interface/TrackWithHistory.h | 6 +++--- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/SimG4Core/Notification/interface/TmpSimVertex.h b/SimG4Core/Notification/interface/TmpSimVertex.h index d66417f423f72..c95f909f80e1d 100644 --- a/SimG4Core/Notification/interface/TmpSimVertex.h +++ b/SimG4Core/Notification/interface/TmpSimVertex.h @@ -7,7 +7,7 @@ class TmpSimVertex { public: - TmpSimVertex(const math::XYZVectorD& ip, double it, int iv, int typ=0) + TmpSimVertex(const math::XYZVectorD& ip, double it, int iv, int typ = 0) : ilv_(ip), itime_(it), itrack_(iv), ptype_(typ) {} ~TmpSimVertex() = default; /// index of the parent (-1 if no parent) diff --git a/SimG4Core/Notification/interface/TrackWithHistory.h b/SimG4Core/Notification/interface/TrackWithHistory.h index 1eb24570c0c63..a4982ecc34943 100644 --- a/SimG4Core/Notification/interface/TrackWithHistory.h +++ b/SimG4Core/Notification/interface/TrackWithHistory.h @@ -17,8 +17,8 @@ class TrackWithHistory { /** The constructor is called at time, * when some of the information may not available yet. */ - TrackWithHistory(const G4Track* g4track, int pID); - TrackWithHistory(const G4PrimaryParticle*, int trackID, const math::XYZVectorD &pos, const double time); + TrackWithHistory(const G4Track *g4track, int pID); + TrackWithHistory(const G4PrimaryParticle *, int trackID, const math::XYZVectorD &pos, const double time); ~TrackWithHistory() = default; inline void *operator new(std::size_t); @@ -77,7 +77,7 @@ class TrackWithHistory { int idAtBoundary_{-1}; int procType_{0}; double totalEnergy_; - double time_; // lab system + double time_; // lab system double weight_; math::XYZVectorD momentum_; math::XYZVectorD vertexPosition_;