DBremWatcher Class Reference

#include <DBremWatcher.h>

Inheritance diagram for DBremWatcher:

Public Member Functions
	DBremWatcher (edm::ParameterSet const &p)
void	produce (edm::Event &, const edm::EventSetup &) override
void	update (const BeginOfTrack *trk) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfEvent *event) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const EndOfEvent *event) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfRun *run) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const EndOfTrack *trk) override
	This routine will be called when the appropriate signal arrives. More...
	~DBremWatcher () override
Public Member Functions inherited from SimProducer
const SimProducer &	operator= (const SimProducer &)=delete
void	registerProducts (edm::ProducesCollector producesCollector)
	SimProducer ()
	SimProducer (const SimProducer &)=delete
Public Member Functions inherited from SimWatcher
virtual void	beginRun (edm::EventSetup const &)
bool	isMT () const
const SimWatcher &	operator= (const SimWatcher &)=delete
virtual void	registerConsumes (edm::ConsumesCollector)
	SimWatcher ()
	SimWatcher (const SimWatcher &)=delete
virtual	~SimWatcher ()
Public Member Functions inherited from Observer< const BeginOfTrack *>
	Observer ()
void	slotForUpdate (const BeginOfTrack * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfEvent *>
	Observer ()
void	slotForUpdate (const BeginOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfRun *>
	Observer ()
void	slotForUpdate (const BeginOfRun * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const EndOfEvent *>
	Observer ()
void	slotForUpdate (const EndOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const EndOfTrack *>
	Observer ()
void	slotForUpdate (const EndOfTrack * iT)
virtual	~Observer ()

Private Attributes
G4ThreeVector	aPrimeTraj
double	biasFactor
float	f_energy
G4ThreeVector	finaltraj
bool	foundbrem
float	m_weight
int	MotherId
std::vector< int >	pdgs_
G4ThreeVector	VertexPos

Additional Inherited Members
Protected Member Functions inherited from SimProducer
template<class T >
void	produces ()
template<class T >
void	produces (const std::string &instanceName)
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 17 of file DBremWatcher.h.

Constructor & Destructor Documentation

DBremWatcher()

DBremWatcher::DBremWatcher

(

edm::ParameterSet const &

p

)

Definition at line 28 of file DBremWatcher.cc.

References aPrimeTraj, biasFactor, f_energy, finaltraj, foundbrem, edm::ParameterSet::getUntrackedParameter(), dqmdumpme::k, m_weight, AlCaHLTBitMon_ParallelJobs::p, pdgs_, and VertexPos.

                                                   {
  edm::ParameterSet ps = p.getParameter<edm::ParameterSet>("DBremWatcher");
  biasFactor = ps.getUntrackedParameter<double>("DBremBiasFactor", 1);
  m_weight = 0;
  foundbrem = false;
  finaltraj = G4ThreeVector(0, 0, 0);
  aPrimeTraj = G4ThreeVector(0, 0, 0);
  VertexPos = G4ThreeVector(0, 0, 0);
  f_energy = 0;
  produces<float>("DBremEventWeight");
  produces<float>("DBremLocationX");
  produces<float>("DBremLocationY");
  produces<float>("DBremLocationZ");
  //produces<std::string>("DBremMaterial");
  produces<float>("DBremAngle");
  produces<float>("DBremInitialEnergy");
  produces<float>("DBremFinalEnergy");
  produces<float>("BiasFactor");
  pdgs_ = ps.getUntrackedParameter<std::vector<int>>("PDGCodes");
  edm::LogInfo("DBremWatcher") << "DBremWatcher:: Save Sim Track if PDG code "
                               << "is one from the list of " << pdgs_.size() << " items";
  for (unsigned int k = 0; k < pdgs_.size(); ++k)
    edm::LogInfo("DBremWatcher") << "[" << k << "] " << pdgs_[k];
}

~DBremWatcher()

DBremWatcher::~DBremWatcher ( )

override

Definition at line 53 of file DBremWatcher.cc.

{}

Member Function Documentation

produce()

void DBremWatcher::produce ( edm::Event &  fEvent, const edm::EventSetup &    )

overridevirtual

Implements SimProducer.

Definition at line 116 of file DBremWatcher.cc.

References aPrimeTraj, biasFactor, f_energy, hcaldqm::fEvent, finaltraj, foundbrem, m_weight, eostools::move(), funct::pow(), mathSSE::sqrt(), and VertexPos.

                                                                 {
  if (foundbrem) {
    std::unique_ptr<float> weight = std::make_unique<float>(m_weight);
    fEvent.put(std::move(weight), "DBremEventWeight");
    std::unique_ptr<float> vtxposx = std::make_unique<float>(VertexPos.x());
    std::unique_ptr<float> vtxposy = std::make_unique<float>(VertexPos.y());
    std::unique_ptr<float> vtxposz = std::make_unique<float>(VertexPos.z());
    fEvent.put(std::move(vtxposx), "DBremLocationX");
    fEvent.put(std::move(vtxposy), "DBremLocationY");
    fEvent.put(std::move(vtxposz), "DBremLocationZ");
    std::unique_ptr<float> finalE = std::make_unique<float>(f_energy / GeV);
    fEvent.put(std::move(finalE), "DBremFinalEnergy");
    float deflectionAngle = -1;
    float initialEnergy = sqrt(pow(finaltraj.x() + aPrimeTraj.x(), 2) + pow(finaltraj.y() + aPrimeTraj.y(), 2) +
                               pow(finaltraj.z() + aPrimeTraj.z(), 2) + pow(0.1056 * GeV, 2));
    G4ThreeVector mother(
        finaltraj.x() + aPrimeTraj.x(), finaltraj.y() + aPrimeTraj.y(), finaltraj.z() + aPrimeTraj.z());
    deflectionAngle = mother.angle(finaltraj);
    std::unique_ptr<float> dAngle = std::make_unique<float>(deflectionAngle);
    std::unique_ptr<float> initialE = std::make_unique<float>(initialEnergy / GeV);
    fEvent.put(std::move(dAngle), "DBremAngle");
    fEvent.put(std::move(initialE), "DBremInitialEnergy");
    std::unique_ptr<float> bias = std::make_unique<float>(biasFactor);
    fEvent.put(std::move(bias), "BiasFactor");
  } else {
    std::unique_ptr<float> weight = std::make_unique<float>(0.);
    fEvent.put(std::move(weight), "DBremEventWeight");
  }
}

update() [1/5]

void DBremWatcher::update ( const BeginOfTrack *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack *>.

Definition at line 55 of file DBremWatcher.cc.

References aPrimeTraj, f_energy, finaltraj, spr::find(), foundbrem, LogDebug, m_weight, dileptonTrigSettings_cff::pdg, pdgs_, TrackInformation::storeTrack(), and VertexPos.

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

                                                 {
  G4Track* theTrack = (G4Track*)((*trk)());
  TrackInformation* trkInfo = (TrackInformation*)(theTrack->GetUserInformation());
  if (trkInfo) {
    int pdg = theTrack->GetDefinition()->GetPDGEncoding();
    G4ThreeVector Vpos = theTrack->GetVertexPosition();
    const G4VProcess* TrPro = theTrack->GetCreatorProcess();
    if (TrPro != nullptr) {
      if ((theTrack->GetCreatorProcess()->GetProcessName()) == "muDBrem") {
        if (std::find(pdgs_.begin(), pdgs_.end(), pdg) == pdgs_.end()) {
          //Found the deflected muon track
          trkInfo->storeTrack(true);
          if (!theTrack->IsGoodForTracking()) {
            theTrack->SetGoodForTrackingFlag(true);
          }
          f_energy = theTrack->GetTotalEnergy();
          foundbrem = true;
          finaltraj = theTrack->GetMomentum();
        } else {
          m_weight = theTrack->GetWeight();
        }
      }
    }
    if (std::find(pdgs_.begin(), pdgs_.end(), pdg) != pdgs_.end()) {
      //Found an A'
      trkInfo->storeTrack(true);
      VertexPos = Vpos;
      aPrimeTraj = theTrack->GetMomentum();
      LogDebug("DBremWatcher") << "Save SimTrack the Track " << theTrack->GetTrackID() << " Type "
                               << theTrack->GetDefinition()->GetParticleName() << " Momentum "
                               << theTrack->GetMomentum() / MeV << " MeV/c";
    }
  }
}

update() [2/5]

void DBremWatcher::update ( const BeginOfEvent *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 92 of file DBremWatcher.cc.

References foundbrem, and unpackData-CaloStage2::pname.

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

                                                   {
  G4String pname = "muDBrem";
  G4ProcessTable* ptable = G4ProcessTable::GetProcessTable();
  G4bool state = true;
  ptable->SetProcessActivation(pname, state);
  foundbrem = false;
}

update() [3/5]

void DBremWatcher::update ( const EndOfEvent *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent *>.

Definition at line 100 of file DBremWatcher.cc.

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

{}

update() [4/5]

void DBremWatcher::update ( const BeginOfRun *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 90 of file DBremWatcher.cc.

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

{}

update() [5/5]

void DBremWatcher::update ( const EndOfTrack *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack *>.

Definition at line 102 of file DBremWatcher.cc.

References spr::find(), dileptonTrigSettings_cff::pdg, pdgs_, and TrackInformation::storeTrack().

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

                                               {
  G4Track* theTrack = (G4Track*)((*trk)());
  TrackInformation* trkInfo = (TrackInformation*)(theTrack->GetUserInformation());
  const G4VProcess* TrPro = theTrack->GetCreatorProcess();
  if (trkInfo && TrPro != nullptr) {
    int pdg = theTrack->GetDefinition()->GetPDGEncoding();

    if (std::find(pdgs_.begin(), pdgs_.end(), pdg) == pdgs_.end() &&
        (theTrack->GetCreatorProcess()->GetProcessName()) == "muDBrem") {
      trkInfo->storeTrack(true);
    }
  }
}

Member Data Documentation

aPrimeTraj

G4ThreeVector DBremWatcher::aPrimeTraj

private

Definition at line 39 of file DBremWatcher.h.

Referenced by DBremWatcher(), produce(), and update().

biasFactor

double DBremWatcher::biasFactor

private

Definition at line 37 of file DBremWatcher.h.

Referenced by DBremWatcher(), and produce().

f_energy

float DBremWatcher::f_energy

private

Definition at line 42 of file DBremWatcher.h.

Referenced by DBremWatcher(), produce(), and update().

finaltraj

G4ThreeVector DBremWatcher::finaltraj

private

Definition at line 40 of file DBremWatcher.h.

Referenced by DBremWatcher(), produce(), and update().

foundbrem

bool DBremWatcher::foundbrem

private

Definition at line 38 of file DBremWatcher.h.

Referenced by DBremWatcher(), produce(), and update().

m_weight

float DBremWatcher::m_weight

private

Definition at line 36 of file DBremWatcher.h.

Referenced by DBremWatcher(), produce(), and update().

MotherId

int DBremWatcher::MotherId

private

Definition at line 35 of file DBremWatcher.h.

pdgs_

std::vector<int> DBremWatcher::pdgs_

private

Definition at line 34 of file DBremWatcher.h.

Referenced by DBremWatcher(), and update().

VertexPos

G4ThreeVector DBremWatcher::VertexPos

private

Definition at line 41 of file DBremWatcher.h.

Referenced by DBremWatcher(), produce(), and update().