MonopoleSteppingAction Class Reference

#include <MonopoleSteppingAction.h>

Inheritance diagram for MonopoleSteppingAction:

Public Member Functions
void	beginRun (edm::EventSetup const &) override
	MonopoleSteppingAction (edm::ParameterSet const &p)
void	registerConsumes (edm::ConsumesCollector) override
void	update (const BeginOfRun *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfTrack *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const G4Step *) override
	This routine will be called when the appropriate signal arrives. More...
	~MonopoleSteppingAction () override
Public Member Functions inherited from SimWatcher
bool	isMT () const
const SimWatcher &	operator= (const SimWatcher &)=delete
	SimWatcher ()
	SimWatcher (const SimWatcher &)=delete
virtual	~SimWatcher ()
Public Member Functions inherited from Observer< const BeginOfRun * >
	Observer ()
void	slotForUpdate (const BeginOfRun *iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfTrack * >
	Observer ()
void	slotForUpdate (const BeginOfTrack *iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const G4Step * >
	Observer ()
void	slotForUpdate (const G4Step *iT)
virtual	~Observer ()

Private Attributes
bool	actOnTrack
double	bZ
double	cInMByS
double	cMevToJ
double	cMeVToKgMByS
double	dirxStart
double	diryStart
double	dirzStart
double	eStart
double	magCharge
bool	mode
std::vector< int >	pdgCode
double	pxStart
double	pyStart
double	pzStart
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	tok_bFieldH_

Additional Inherited Members
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 16 of file MonopoleSteppingAction.h.

Constructor & Destructor Documentation

MonopoleSteppingAction::MonopoleSteppingAction

(

edm::ParameterSet const &

p

)

Definition at line 19 of file MonopoleSteppingAction.cc.

References edm::ParameterSet::getUntrackedParameter(), and mode.

                                                                       : actOnTrack(false), bZ(0) {
  mode = p.getUntrackedParameter<bool>("ChangeFromFirstStep", true);
  edm::LogVerbatim("SimG4CoreWatcher") << "MonopoleSeppingAction set mode for"
                                       << " start at first step to " << mode;
}

MonopoleSteppingAction::~MonopoleSteppingAction ( )

override

Definition at line 25 of file MonopoleSteppingAction.cc.

{}

Member Function Documentation

void MonopoleSteppingAction::beginRun ( edm::EventSetup const &  es )

overridevirtual

Reimplemented from SimWatcher.

Definition at line 32 of file MonopoleSteppingAction.cc.

References ecalTB2006H4_GenSimDigiReco_cfg::bField, bZ, edm::EventSetup::getData(), MagneticField::inTesla(), AlCaHLTBitMon_ParallelJobs::p, tok_bFieldH_, x, y, and z.

                                                             {
  const MagneticField *bField = &es.getData(tok_bFieldH_);
  const GlobalPoint p(0, 0, 0);
  bZ = (bField->inTesla(p)).z();
  edm::LogVerbatim("SimG4CoreWatcher") << "Magnetic Field (X): " << (bField->inTesla(p)).x()
                                       << " Y: " << (bField->inTesla(p)).y() << " Z: " << bZ;
}

void MonopoleSteppingAction::registerConsumes ( edm::ConsumesCollector  cc )

overridevirtual

Reimplemented from SimWatcher.

Definition at line 27 of file MonopoleSteppingAction.cc.

References edm::BeginRun, edm::ConsumesCollector::esConsumes(), dumpMFGeometry_cfg::MagneticField, and tok_bFieldH_.

                                                                     {
  tok_bFieldH_ = cc.esConsumes<MagneticField, IdealMagneticFieldRecord, edm::Transition::BeginRun>();
  edm::LogVerbatim("SimG4CoreWatcher") << "MonopoleSteppingAction::Initialize ESGetToken for MagneticField";
}

void MonopoleSteppingAction::update ( const BeginOfRun *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 40 of file MonopoleSteppingAction.cc.

References cInMByS, cMevToJ, cMeVToKgMByS, e_SI, cuy::ii, magCharge, pdgCode, edm::second(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                      {
  G4ParticleTable *partTable = G4ParticleTable::GetParticleTable();
  magCharge = CLHEP::e_SI / CLHEP::fine_structure_const * 0.5;
  for (int ii = 0; ii < partTable->size(); ii++) {
    G4ParticleDefinition *particle = partTable->GetParticle(ii);
    std::string particleName = (particle->GetParticleName()).substr(0, 8);
    if (strcmp(particleName.c_str(), "Monopole") == 0) {
      magCharge = CLHEP::e_SI * ((Monopole *)(particle))->MagneticCharge();
      pdgCode.push_back(particle->GetPDGEncoding());
    }
  }
  edm::LogVerbatim("SimG4CoreWatcher") << "MonopoleSeppingAction Finds " << pdgCode.size() << " candidates";
  for (unsigned int ii = 0; ii < pdgCode.size(); ++ii) {
    edm::LogVerbatim("SimG4CoreWatcher") << "PDG Code[" << ii << "] = " << pdgCode[ii];
  }
  cMevToJ = CLHEP::e_SI / CLHEP::eV;
  cMeVToKgMByS = CLHEP::e_SI * CLHEP::meter / (CLHEP::eV * CLHEP::c_light * CLHEP::second);
  cInMByS = CLHEP::c_light * CLHEP::second / CLHEP::meter;
  edm::LogVerbatim("SimG4CoreWatcher") << "MonopoleSeppingAction Constants:"
                                       << " MevToJoules  = " << cMevToJ << ", MevToKgm/s  = " << cMeVToKgMByS
                                       << ", c in m/s    = " << cInMByS << ", mag. charge = " << magCharge;
}

void MonopoleSteppingAction::update ( const BeginOfTrack *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 63 of file MonopoleSteppingAction.cc.

References actOnTrack, submitPVResolutionJobs::count, dirxStart, diryStart, dirzStart, eStart, GeV, LogDebug, pdgCode, pxStart, pyStart, and pzStart.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                           {
  actOnTrack = false;
  if (!pdgCode.empty()) {
    const G4Track *aTrack = (*trk)();
    int code = aTrack->GetDefinition()->GetPDGEncoding();
    if (std::count(pdgCode.begin(), pdgCode.end(), code) > 0) {
      actOnTrack = true;
      eStart = aTrack->GetTotalEnergy();
      pxStart = aTrack->GetMomentum().x();
      pyStart = aTrack->GetMomentum().y();
      pzStart = aTrack->GetMomentum().z();
      dirxStart = aTrack->GetMomentumDirection().x();
      diryStart = aTrack->GetMomentumDirection().y();
      dirzStart = aTrack->GetMomentumDirection().z();
      LogDebug("SimG4CoreWatcher") << "MonopoleSeppingAction Track " << code << " Flag " << actOnTrack
                                   << " (px,py,pz,E) = (" << pxStart / GeV << ", " << pyStart / GeV << ", "
                                   << pzStart / GeV << ", " << eStart / GeV << ")";
    }
  }
}

void MonopoleSteppingAction::update ( const G4Step *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 84 of file MonopoleSteppingAction.cc.

References actOnTrack, bZ, cInMByS, cMevToJ, cMeVToKgMByS, dirxStart, diryStart, eStart, isotrackTrainRegressor::fac0(), isotrackTrainRegressor::fac1(), isotrackTrainRegressor::fac2(), LogDebug, visualization-live-secondInstance_cfg::m, magCharge, mode, PVValHelper::pT, HLT_FULL_cff::pT1, pxStart, pyStart, pzStart, edm::second(), mathSSE::sqrt(), beampixel_dqm_sourceclient-live_cfg::xStep, beampixel_dqm_sourceclient-live_cfg::yStep, and beampixel_dqm_sourceclient-live_cfg::zStep.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                       {
  if (actOnTrack) {
    double eT, pT, pZ, tStep;
    G4Track *aTrack = aStep->GetTrack();
    G4ThreeVector initialPosition(0, 0, 0);
    if (mode) {
      tStep = aTrack->GetGlobalTime();
      eT = eStart * std::sqrt(dirxStart * dirxStart + diryStart * diryStart);
      pT = std::sqrt(pxStart * pxStart + pyStart * pyStart);
      pZ = pzStart;
    } else {
      const G4ThreeVector &dirStep = aTrack->GetMomentumDirection();
      double lStep = aTrack->GetStepLength();
      double xStep = aTrack->GetPosition().x() - lStep * dirStep.x();
      double yStep = aTrack->GetPosition().y() - lStep * dirStep.y();
      double zStep = aTrack->GetPosition().z() - lStep * dirStep.z();
      double vStep = aTrack->GetVelocity();
      initialPosition = G4ThreeVector(xStep, yStep, zStep);
      tStep = (vStep > 0. ? (lStep / vStep) : 0.);
      double eStep = aTrack->GetTotalEnergy();
      eT = eStep * dirStep.perp();
      pT = aTrack->GetMomentum().perp();
      pZ = aTrack->GetMomentum().z();
    }
    LogDebug("SimG4CoreWatcher") << "MonopoleSeppingAction: tStep " << tStep << " eT " << eT << " pT " << pT << " pZ "
                                 << pZ;
    double eT1 = eT * cMevToJ;
    double pT1 = pT * cMeVToKgMByS;
    double pZ1 = pZ * cMeVToKgMByS;
    double ts1 = tStep / CLHEP::second;
    double eT2 = (eT1 > 0. ? (1. / eT1) : 0.);
    double fac0 = magCharge * bZ * cInMByS;
    double fac2 = pZ1 * cInMByS * eT2;
    double fac1 = fac0 * cInMByS * ts1 * eT2 + fac2;
    double fact1 = (eT1 / fac0) * (std::sqrt(1 + fac1 * fac1) - std::sqrt(1 + fac2 * fac2));
    double fact2 = (pT1 / (magCharge * bZ) * (asinh(fac1) - asinh(fac2)));
    LogDebug("SimG4CoreWatcher") << "MonopoleSeppingAction: Factor eT = " << eT << " " << eT1 << " " << eT2
                                 << ", pT = " << pT << " " << pT1 << ", pZ = " << pZ << " " << pZ1
                                 << ", time = " << tStep << " " << ts1 << ", Factors " << fac0 << " " << fac1 << " "
                                 << fac2 << ", Final ones " << fact1 << " " << fact2;
    G4ThreeVector ez(0, 0, 1), et(std::sqrt(0.5), std::sqrt(0.5), 0);
    G4ThreeVector displacement = (fact1 * ez + fact2 * et) * CLHEP::m;

    LogDebug("SimG4CoreWatcher") << "MonopoleSeppingAction: Initial " << initialPosition << " Displacement "
                                 << displacement << " Final " << (initialPosition + displacement);
    aTrack->SetPosition(initialPosition + displacement);
  }
}

Member Data Documentation

bool MonopoleSteppingAction::actOnTrack

private

Definition at line 31 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::bZ

private

Definition at line 35 of file MonopoleSteppingAction.h.

Referenced by beginRun(), and update().

double MonopoleSteppingAction::cInMByS

private

Definition at line 35 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::cMevToJ

private

Definition at line 35 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::cMeVToKgMByS

private

Definition at line 35 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::dirxStart

private

Definition at line 34 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::diryStart

private

Definition at line 34 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::dirzStart

private

Definition at line 34 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::eStart

private

Definition at line 33 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::magCharge

private

Definition at line 35 of file MonopoleSteppingAction.h.

Referenced by update().

bool MonopoleSteppingAction::mode

private

Definition at line 31 of file MonopoleSteppingAction.h.

Referenced by MonopoleSteppingAction(), and update().

std::vector<int> MonopoleSteppingAction::pdgCode

private

Definition at line 32 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::pxStart

private

Definition at line 33 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::pyStart

private

Definition at line 33 of file MonopoleSteppingAction.h.

Referenced by update().

double MonopoleSteppingAction::pzStart

private

Definition at line 33 of file MonopoleSteppingAction.h.

Referenced by update().

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> MonopoleSteppingAction::tok_bFieldH_

private

Definition at line 30 of file MonopoleSteppingAction.h.

Referenced by beginRun(), and registerConsumes().