SuepDecay Class Reference

#include <SuepDecay.h>

Inheritance diagram for SuepDecay:

Public Member Functions
bool	canVetoProcessLevel () override
bool	doVetoProcessLevel (Pythia8::Event &event) override
bool	initAfterBeams () override
	SuepDecay (const edm::ParameterSet &iConfig)
	~SuepDecay () override

Protected Attributes
int	idDark_
int	idMediator_
float	mDark_
std::unique_ptr< SuepShower >	suep_shower_
float	temperature_

Detailed Description

Definition at line 12 of file SuepDecay.h.

Constructor & Destructor Documentation

SuepDecay()

SuepDecay::SuepDecay

(

const edm::ParameterSet &

iConfig

)

Definition at line 3 of file SuepDecay.cc.

    : idMediator_(iConfig.getParameter<int>("idMediator")),
      idDark_(iConfig.getParameter<int>("idDark")),
      temperature_(iConfig.getParameter<double>("temperature")) {}

~SuepDecay()

SuepDecay::~SuepDecay ( )

inlineoverride

Definition at line 15 of file SuepDecay.h.

{}

Member Function Documentation

canVetoProcessLevel()

bool SuepDecay::canVetoProcessLevel ( )

inlineoverride

Definition at line 19 of file SuepDecay.h.

{ return true; }

doVetoProcessLevel()

bool SuepDecay::doVetoProcessLevel ( Pythia8::Event &  event )

override

Definition at line 23 of file SuepDecay.cc.

References mps_fire::i, idDark_, idMediator_, mDark_, and suep_shower_.

                                                      {
  Pythia8::Vec4 pMediator, pDark;

  // Find the mediator in the event
  for (int i = 0; i < event.size(); ++i) {
    //mediator w/ distinct daughters = last copy (decayed)
    if (event[i].id() == idMediator_ && event[i].daughter1() != event[i].daughter2() && event[i].daughter1() > 0 &&
        event[i].daughter2() > 0) {
      pMediator = event[i].p();
      // undo mediator decay
      event[i].undoDecay();

      // Generate the shower, output are 4 vectors in the rest frame of the shower, adding energy here avoids issues if scalar is off-shell
      std::vector<Pythia8::Vec4> suep_shower_fourmomenta = suep_shower_->generateShower(pMediator.mCalc());
      // Loop over hidden sector mesons and append to the event
      int firstDaughter = event.size();
      for (auto& pDark : suep_shower_fourmomenta) {
        // Boost to the lab frame, i.e. apply the mediator boost
        pDark.bst(pMediator);
        // Append particle to the event w/ hidden meson pdg code. Magic number 91 means it is produced as a normal decay product
        event.append(idDark_, 91, i, 0, 0, 0, 0, 0, pDark.px(), pDark.py(), pDark.pz(), pDark.e(), mDark_);
      }

      // Change the status code of the mediator to reflect that it has decayed.
      event[i].statusNeg();

      //set daughters of the mediator: daughter1 < daughter2 > 0 -> the particle has a range of decay products from daughter1 to daughter2
      event[i].daughters(firstDaughter, event.size() - 1);
      break;
    }
  }

  //allow event to continue
  return false;
}

initAfterBeams()

bool SuepDecay::initAfterBeams ( )

override

Definition at line 8 of file SuepDecay.cc.

References idDark_, idMediator_, mDark_, suep_shower_, and temperature_.

                               {
  mDark_ = particleDataPtr->m0(idDark_);
  bool medDecay = particleDataPtr->mayDecay(idMediator_);
  if (!medDecay) {
    infoPtr->errorMsg("Error in SuepDecay::initAfterBeams: mediator decay should be enabled");
    return false;
  }

  //construct the shower helper
  suep_shower_ = std::make_unique<SuepShower>(mDark_, temperature_, rndmPtr);

  return true;
}

Member Data Documentation

idDark_

int SuepDecay::idDark_

protected

Definition at line 23 of file SuepDecay.h.

Referenced by doVetoProcessLevel(), and initAfterBeams().

idMediator_

int SuepDecay::idMediator_

protected

Definition at line 23 of file SuepDecay.h.

Referenced by doVetoProcessLevel(), and initAfterBeams().

mDark_

float SuepDecay::mDark_

protected

Definition at line 24 of file SuepDecay.h.

Referenced by doVetoProcessLevel(), and initAfterBeams().

suep_shower_

std::unique_ptr<SuepShower> SuepDecay::suep_shower_

protected

Definition at line 25 of file SuepDecay.h.

Referenced by doVetoProcessLevel(), and initAfterBeams().

temperature_

float SuepDecay::temperature_

protected

Definition at line 24 of file SuepDecay.h.

Referenced by initAfterBeams().