HllgamRepeatTimeShower Class Reference

Tool that will repeatedly perform time showering on the intial photons from a Higgs decay until at least one is offshell. More...

#include <HllgamRepeatTimeShower.h>

Inheritance diagram for HllgamRepeatTimeShower:

Collaboration diagram for HllgamRepeatTimeShower:

Public Member Functions
virtual	~HllgamRepeatTimeShower ()
	Destructor. More...
StatusCode	initialize () override
	AlgTool initialize method. More...
StatusCode	finalize () override
	AlgTool finalize method. More...
StatusCode	ModifyPythiaEvent (Pythia8::Pythia &pythia) const override
	Update the pythia event. More...
virtual double	CrossSectionScaleFactor () const override
	Return how much the cross section is modified. More...

Private Attributes
unsigned long	m_nPass {0}
unsigned long	m_nVetos {0}

Detailed Description

Tool that will repeatedly perform time showering on the intial photons from a Higgs decay until at least one is offshell.

Definition at line 16 of file HllgamRepeatTimeShower.h.

Constructor & Destructor Documentation

~HllgamRepeatTimeShower()

virtual HllgamRepeatTimeShower::~HllgamRepeatTimeShower ( )

inlinevirtual

Destructor.

Definition at line 24 of file HllgamRepeatTimeShower.h.

{};

Member Function Documentation

CrossSectionScaleFactor()

double HllgamRepeatTimeShower::CrossSectionScaleFactor ( ) const

overridevirtual

Return how much the cross section is modified.

Should only be called once all events have been processed

Definition at line 33 of file HllgamRepeatTimeShower.cxx.

                                                               {
    unsigned long sum = m_nVetos  +  m_nPass;
    if(sum == 0 || m_nPass==0 ){
      ATH_MSG_WARNING( "Are you sure that this code has been called correctly??  " );
      return 1.;
    }
    
    return double (m_nPass)/double ( sum );
  }

finalize()

StatusCode HllgamRepeatTimeShower::finalize ( )

override

AlgTool finalize method.

Definition at line 22 of file HllgamRepeatTimeShower.cxx.

                                             {
    
    ATH_MSG_INFO( "*************************************************************" );
    ATH_MSG_INFO( "N retries  " << m_nVetos );
    ATH_MSG_INFO( "N pass     " << m_nPass );
    ATH_MSG_INFO( "*************************************************************" );
    
    ATH_MSG_INFO( "Finalisation of " << name() << " was successful" );
    return StatusCode::SUCCESS;
  }

initialize()

StatusCode HllgamRepeatTimeShower::initialize ( )

override

AlgTool initialize method.

Definition at line 8 of file HllgamRepeatTimeShower.cxx.

                                                {
    ATH_MSG_INFO( "Initialisation of " << name() << " was successful" );
    
    ATH_MSG_INFO( "********************************************************************************" );
    ATH_MSG_INFO( "********************************************************************************" );
    ATH_MSG_INFO( "***               Will run with repeated time showers of H decay             ***" );
    ATH_MSG_INFO( "***            N.B. only make sense for gamma gamma decays                   ***" );
    ATH_MSG_INFO( "********************************************************************************" );
    ATH_MSG_INFO( "********************************************************************************" );
    
    return StatusCode::SUCCESS;
  }

ModifyPythiaEvent()

StatusCode HllgamRepeatTimeShower::ModifyPythiaEvent ( Pythia8::Pythia &  pythia ) const

override

Update the pythia event.

Definition at line 43 of file HllgamRepeatTimeShower.cxx.

                                                                                  {
    
    //Index of the photons for the Higgs decay
    std::vector<int> photons;
    
    int trials = 0;
    
    bool checkEvent = true;
    
    Pythia8::Event& event = pythia.event;
    
    while(checkEvent){
      
      int higgsIndex = -1;
      for(int ii=0; ii != event.size(); ++ii){
        // status 22 is intermediate hard process.
        // I think that is correct here because although the output
        // from Powheg is the Higgs, Pythia should still consider the
        // decay products of the Higgs to be part of the hard process
        if(event[ii].id() == 25 && abs(event[ii].status()) == 22){
          higgsIndex = ii;
          break;
        }
      }
      
      if(higgsIndex <0){
        // didn't find a Higgs - barf
        ATH_MSG_ERROR("No Higgs in event record");
        return StatusCode::FAILURE;
      }
      
      std::vector<int> daughters = event[higgsIndex].daughterList();
      while(daughters.size()==1){
        daughters=event[daughters[0]].daughterList();
      }
      
      std::vector<int> photons;
      
      for(const auto d: daughters){
        if(event[d].id() == 22){
          photons.push_back(d);
        }
      }
      
      if(photons.size() < 2){
        ATH_MSG_ERROR("Wrong number of photons from H decay. Expected >= 2, got " + std::to_string(photons.size()));
        return StatusCode::FAILURE;
      }
      
      for(const auto p: photons){
        daughters = event[p].daughterList();
        while(daughters.size() == 1){
          daughters = event[daughters[0]].daughterList();
        }
        
        if(daughters.size() < 2) continue;
        size_t nLeptons=0;
        
        for(const auto d: daughters){
          if(event[d].isLepton()) ++nLeptons;
        }
        
        if(nLeptons >= 2){
          checkEvent = false;
          break;
        }
      }
      
      if(checkEvent){
        pythia.forceTimeShower( photons[0], photons[1], 1000000);
        ++trials;
        ++m_nVetos;
        if(trials % 1000 == 0 ){
          ATH_MSG_WARNING("Has the decay been setup in Pythia correctly?  Possibly stuck in loop so far " <<  trials  << " attempts have been made");
        }
      }else{
       pythia.forceHadronLevel();
       ++m_nPass;
      }
    }
    
    return StatusCode::SUCCESS;
  }

Member Data Documentation

m_nPass

unsigned long HllgamRepeatTimeShower::m_nPass {0}

mutableprivate

Definition at line 40 of file HllgamRepeatTimeShower.h.

m_nVetos

unsigned long HllgamRepeatTimeShower::m_nVetos {0}

mutableprivate

Definition at line 41 of file HllgamRepeatTimeShower.h.

---

The documentation for this class was generated from the following files:

• HllgamRepeatTimeShower.h
• HllgamRepeatTimeShower.cxx