#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 ( 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 ( )

inlineoverride

Definition at line 15 of file SuepDecay.h.

References initAfterBeams().

15 {}

Member Function Documentation

bool SuepDecay::canVetoProcessLevel ( )

inlineoverride

Definition at line 19 of file SuepDecay.h.

References doVetoProcessLevel().

19 { return true; }

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_.

Referenced by canVetoProcessLevel().

                                                       {
   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;
 }

bool SuepDecay::initAfterBeams ( )

override

Definition at line 8 of file SuepDecay.cc.

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

Referenced by ~SuepDecay().

                                {
   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;
 }