SuppressMPI.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#include "UserHooksUtils.h"
#include "UserSetting.h"
#include "Pythia8_i/UserHooksFactory.h"
#include <stdexcept>
#include <iostream>
 
namespace Pythia8{
  class SuppressMPI;
}
 
Pythia8_UserHooks::UserHooksFactory::Creator<Pythia8::SuppressMPI> SuppressMPICreator("SuppressMPI");
 
namespace Pythia8{
 
  class SuppressMPI : public UserHooks{
    
  public:
    
    SuppressMPI():
    m_pTCut("SuppressMPI:PTCut", 10.), m_nMPIVeto(3){
      
      std::cout<<"************************************************************"<<std::endl;
      std::cout<<"*                                                          *"<<std::endl;
      std::cout<<"*        Suppressing MPI emissions with UserHook!          *"<<std::endl;
      std::cout<<"*                                                          *"<<std::endl;
      std::cout<<"************************************************************"<<std::endl;
 
    }
    
    ~SuppressMPI(){}
    
    bool doVetoMPIStep(int nMPI, const Event &event){
            
      // MPI 1 is the hard process.  We do not veto that!
      if(nMPI < 2){
        return false;
      }
      
      // start at the end of the event record and work back
      // This is prior to showering, so there should be at most 2 new MPI emissions
      // event[0] is documentation, so stop before that.
      size_t nEmissions=0;
      for(int ii=event.size()-1; ii > 0 && nEmissions != 2; --ii){
        if(event[ii].status() != 33) continue;
        if(event[ii].pT() > m_pTCut(settingsPtr)){
          return true;
        }
        
        ++nEmissions;
      }
      
      return false;
    }
    
    bool canVetoMPIStep(){return true;}
    int numberVetoMPIStep(){return m_nMPIVeto;}
    bool canVetoISREmission(){return false;}
    bool canVetoFSREmission(){return false;}
    
  private:
    
    Pythia8_UserHooks::UserSetting<double> m_pTCut;
    
    int m_nMPIVeto;
    
    
  };
  
  
 
}