#include <JetMatchingMadgraph.h>

Inheritance diagram for gen::JetMatchingMadgraph:

Public Member Functions
	JetMatchingMadgraph (const edm::ParameterSet &params)

	~JetMatchingMadgraph () override

Public Member Functions inherited from gen::JetMatching
virtual const std::vector< int > *	getPartonList ()

virtual bool	initAfterBeams ()

bool	isMatchingDone ()

	JetMatching (const edm::ParameterSet &params)

void	resetMatchingStatus ()

virtual	~JetMatching ()

Protected Member Functions
void	beforeHadronisation (const lhef::LHEEvent *event) override

void	beforeHadronisationExec () override

std::set< std::string >	capabilities () const override

double	getJetEtaMax () const override

template<typename T >
T	getParameter (const std::string &var, const T &defValue=T()) const

void	init (const lhef::LHERunInfo *runInfo) override

int	match (const lhef::LHEEvent partonLevel, const std::vector< fastjet::PseudoJet > jetInput) override

Static Protected Member Functions
template<typename T >
static T	getParameter (const std::map< std::string, std::string > &params, const std::string &var, const T &defValue=T())

template<typename T >
static T	parseParameter (const std::string &value)

template<>
std::string	parseParameter (const std::string &value)

template<>
bool	parseParameter (const std::string &value_)

template<typename T >
static void	updateOrDie (const std::map< std::string, std::string > &params, T &param, const std::string &name)

Protected Attributes
bool	eventInitialized

bool	exclusive

std::map< std::string, std::string >	mgParams

bool	runInitialized

bool	soup

Protected Attributes inherited from gen::JetMatching
bool	fMatchingStatus

Additional Inherited Members
Static Public Member Functions inherited from gen::JetMatching
static std::unique_ptr< JetMatching >	create (const edm::ParameterSet &params)

Detailed Description

Definition at line 8 of file JetMatchingMadgraph.h.

Constructor & Destructor Documentation

◆ JetMatchingMadgraph()

gen::JetMatchingMadgraph::JetMatchingMadgraph ( const edm::ParameterSet & params )

Definition at line 110 of file JetMatchingMadgraph.cc.

       : JetMatching(params), runInitialized(false) {
     std::string mode = params.getParameter<std::string>("mode");
     if (mode == "inclusive") {
       soup = false;
       exclusive = false;
     } else if (mode == "exclusive") {
       soup = false;
       exclusive = true;
     } else if (mode == "auto")
       soup = true;
     else
       throw cms::Exception("Generator|LHEInterface") << "Madgraph jet matching scheme requires \"mode\" "
                                                         "parameter to be set to either \"inclusive\", "
                                                         "\"exclusive\" or \"auto\"."
                                                      << std::endl;
  
     memain_.etcjet = 0.;
     memain_.rclmax = 0.0;
     memain_.clfact = 0.0;
     memain_.ktsche = 0.0;
     memain_.etaclmax = params.getParameter<double>("MEMAIN_etaclmax");
     memain_.qcut = params.getParameter<double>("MEMAIN_qcut");
     memain_.minjets = params.getParameter<int>("MEMAIN_minjets");
     memain_.maxjets = params.getParameter<int>("MEMAIN_maxjets");
     memain_.showerkt = params.getParameter<double>("MEMAIN_showerkt");
     memain_.nqmatch = params.getParameter<int>("MEMAIN_nqmatch");
     outtree_.flag = params.getParameter<int>("outTree_flag");
     std::string list_excres = params.getParameter<std::string>("MEMAIN_excres");
     std::vector<std::string> elems;
     std::stringstream ss(list_excres);
     std::string item;
     int index = 0;
     while (std::getline(ss, item, ',')) {
       elems.push_back(item);
       memain_.excres[index] = std::atoi(item.c_str());
       index++;
     }
     memain_.nexcres = index;
   }

References gen::MEMAIN::clfact, gen::MEMAIN::etaclmax, gen::MEMAIN::etcjet, Exception, exclusive, gen::MEMAIN::excres, gen::OUTTREE::flag, B2GTnPMonitor_cfi::item, gen::MEMAIN::ktsche, gen::MEMAIN::maxjets, gen::memain_, gen::MEMAIN::minjets, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, gen::MEMAIN::nexcres, gen::MEMAIN::nqmatch, gen::outtree_, CalibrationSummaryClient_cfi::params, gen::MEMAIN::qcut, gen::MEMAIN::rclmax, gen::MEMAIN::showerkt, soup, contentValuesCheck::ss, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~JetMatchingMadgraph()

gen::JetMatchingMadgraph::~JetMatchingMadgraph ( )

override

Definition at line 151 of file JetMatchingMadgraph.cc.

151 {}

Member Function Documentation

◆ beforeHadronisation()

void gen::JetMatchingMadgraph::beforeHadronisation ( const lhef::LHEEvent * event )

overrideprotectedvirtual

Reimplemented from gen::JetMatching.

Definition at line 259 of file JetMatchingMadgraph.cc.

                                                                          {
     if (!runInitialized)
       throw cms::Exception("Generator|PartonShowerVeto") << "Run not initialized in JetMatchingMadgraph" << std::endl;
  
     if (uppriv_.ickkw) {
       std::vector<std::string> comments = event->getComments();
       if (comments.size() == 1) {
         std::istringstream ss(comments[0].substr(1));
         for (int i = 0; i < 1000; i++) {
           double pt;
           ss >> pt;
           if (!ss.good())
             break;
           pypart_.ptpart[i] = pt;
         }
       } else {
         edm::LogWarning("Generator|LHEInterface") << "Expected exactly one comment line per "
                                                      "event containing MadGraph parton scale "
                                                      "information."
                                                   << std::endl;
  
         const lhef::HEPEUP *hepeup = event->getHEPEUP();
         for (int i = 2; i < hepeup->NUP; i++) {
           double mt2 = hepeup->PUP[i][0] * hepeup->PUP[i][0] + hepeup->PUP[i][1] * hepeup->PUP[i][1] +
                        hepeup->PUP[i][4] * hepeup->PUP[i][4];
           pypart_.ptpart[i - 2] = std::sqrt(mt2);
         }
       }
     }
  
     // mgevnt_();
     eventInitialized = true;
   }

References indexGen::comments, eventInitialized, Exception, mps_fire::i, gen::UPPRIV::ickkw, MT2Analyzer::mt2, lhef::HEPEUP::NUP, DiDispStaMuonMonitor_cfi::pt, gen::PYPART::ptpart, lhef::HEPEUP::PUP, gen::pypart_, runInitialized, mathSSE::sqrt(), contentValuesCheck::ss, and gen::uppriv_.

◆ beforeHadronisationExec()

void gen::JetMatchingMadgraph::beforeHadronisationExec ( )

overrideprotectedvirtual

Reimplemented from gen::JetMatching.

Definition at line 293 of file JetMatchingMadgraph.cc.

                                                     {
     mgevnt_();
     eventInitialized = true;
     return;
   }

References eventInitialized, and gen::mgevnt_().

◆ capabilities()

std::set< std::string > gen::JetMatchingMadgraph::capabilities ( ) const

overrideprotectedvirtual

Reimplemented from gen::JetMatching.

Definition at line 155 of file JetMatchingMadgraph.cc.

                                                             {
     std::set<std::string> result;
     result.insert("psFinalState");
     result.insert("hepevt");
     result.insert("pythia6");
     return result;
   }

References mps_fire::result.

◆ getJetEtaMax()

double gen::JetMatchingMadgraph::getJetEtaMax ( ) const

overrideprotectedvirtual

Implements gen::JetMatching.

Definition at line 153 of file JetMatchingMadgraph.cc.

153 { return memain_.etaclmax; }

References gen::MEMAIN::etaclmax, and gen::memain_.

◆ getParameter() [1/2]

template<typename T >

T gen::JetMatchingMadgraph::getParameter	(	const std::map< std::string, std::string > &	params,
		const std::string &	var,
		const T &	defValue = `T()`
	)

staticprotected

Definition at line 96 of file JetMatchingMadgraph.cc.

                                                          {
     std::map<std::string, std::string>::const_iterator pos = params.find(var);
     if (pos == params.end())
       return defValue;
     return parseParameter<T>(pos->second);
   }

References CalibrationSummaryClient_cfi::params, and trigObjTnPSource_cfi::var.

Referenced by getParameter(), and updateOrDie().

◆ getParameter() [2/2]

template<typename T >

T gen::JetMatchingMadgraph::getParameter	(	const std::string &	var,
		const T &	defValue = `T()`
	)		const

protected

Definition at line 106 of file JetMatchingMadgraph.cc.

                                                                                    {
     return getParameter(mgParams, var, defValue);
   }

References getParameter(), mgParams, and trigObjTnPSource_cfi::var.

◆ init()

void gen::JetMatchingMadgraph::init ( const lhef::LHERunInfo * runInfo )

overrideprotectedvirtual

Reimplemented from gen::JetMatching.

Definition at line 206 of file JetMatchingMadgraph.cc.

                                                               {
     // read MadGraph run card
  
     std::map<std::string, std::string> parameters;
  
     std::vector<std::string> header = runInfo->findHeader("MGRunCard");
     if (header.empty())
       throw cms::Exception("Generator|PartonShowerVeto") << "In order to use MadGraph jet matching, "
                                                             "the input file has to contain the corresponding "
                                                             "MadGraph headers."
                                                          << std::endl;
  
     mgParams = parseHeader(header);
  
     // set variables in common block
  
     std::vector<Param> params;
     std::vector<Param> values;
     for (std::map<std::string, std::string>::const_iterator iter = mgParams.begin(); iter != mgParams.end(); ++iter) {
       params.push_back(" " + iter->first);
       values.push_back(iter->second);
     }
  
     // set MG matching parameters
  
     uppriv_.ickkw = getParameter<int>("ickkw", 0);
     memain_.mektsc = getParameter<int>("ktscheme", 0);
  
     header = runInfo->findHeader("MGParamCMS");
  
     std::map<std::string, std::string> mgInfoCMS = parseHeader(header);
  
     for (std::map<std::string, std::string>::const_iterator iter = mgInfoCMS.begin(); iter != mgInfoCMS.end(); ++iter) {
       std::cout << "mgInfoCMS: " << iter->first << " " << iter->second << std::endl;
     }
  
     updateOrDie(mgInfoCMS, memain_.etaclmax, "etaclmax");
     updateOrDie(mgInfoCMS, memain_.qcut, "qcut");
     updateOrDie(mgInfoCMS, memain_.minjets, "minjets");
     updateOrDie(mgInfoCMS, memain_.maxjets, "maxjets");
     updateOrDie(mgInfoCMS, memain_.showerkt, "showerkt");
     updateOrDie(mgInfoCMS, memain_.nqmatch, "nqmatch");
  
     // run Fortran initialization code
  
     int nparam = params.size();
     mginit_(&nparam, &params.front(), &values.front());
     runInitialized = true;
   }

References gather_cfg::cout, gen::MEMAIN::etaclmax, lhef::LHERunInfo::findHeader(), RecoTauValidation_cfi::header, gen::UPPRIV::ickkw, gen::MEMAIN::maxjets, gen::MEMAIN::mektsc, gen::memain_, gen::mginit_(), mgParams, gen::MEMAIN::minjets, gen::MEMAIN::nqmatch, CalibrationSummaryClient_cfi::params, gen::parseHeader(), gen::MEMAIN::qcut, runInitialized, gen::MEMAIN::showerkt, updateOrDie(), gen::uppriv_, and contentValuesCheck::values.

◆ match()

int gen::JetMatchingMadgraph::match	(	const lhef::LHEEvent *	partonLevel,
		const std::vector< fastjet::PseudoJet > *	jetInput
	)

overrideprotectedvirtual

Implements gen::JetMatching.

Definition at line 304 of file JetMatchingMadgraph.cc.

                                                                                                              {
     /*
     if (!showeredFinalState)
         throw cms::Exception("Generator|LHEInterface")
             << "MadGraph matching expected parton shower "
                "final state." << std::endl;
 */
  
     if (!runInitialized)
       throw cms::Exception("Generator|LHEInterface") << "Run not initialized in JetMatchingMadgraph" << std::endl;
  
     if (!eventInitialized)
       throw cms::Exception("Generator|LHEInterface") << "Event not initialized in JetMatchingMadgraph" << std::endl;
  
     if (soup)
       memaev_.iexc = (memaev_.nljets < memain_.maxjets);
     else
       memaev_.iexc = exclusive;
  
     int veto = 0;
     mgveto_(&veto);
     fMatchingStatus = true;
     eventInitialized = false;
  
     return veto;
   }

References eventInitialized, Exception, exclusive, gen::JetMatching::fMatchingStatus, gen::MEMAEV::iexc, gen::MEMAIN::maxjets, gen::memaev_, gen::memain_, gen::mgveto_(), gen::MEMAEV::nljets, runInitialized, soup, and PbPb_ZMuSkimMuonDPG_cff::veto.

◆ parseParameter() [1/3]

template<typename T >

T gen::JetMatchingMadgraph::parseParameter ( const std::string & value )

staticprotected

Definition at line 71 of file JetMatchingMadgraph.cc.

                                                               {
     std::istringstream ss(value);
     T result;
     ss >> result;
     return result;
   }

References mps_fire::result, and contentValuesCheck::ss.

◆ parseParameter() [2/3]

template<>

std::string gen::JetMatchingMadgraph::parseParameter ( const std::string & value )

staticprotected

Definition at line 79 of file JetMatchingMadgraph.cc.

                                                                       {
     std::string result;
     if (!result.empty() && result[0] == '\'')
       result = result.substr(1);
     if (!result.empty() && result[result.length() - 1] == '\'')
       result.resize(result.length() - 1);
     return result;
   }

References mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ parseParameter() [3/3]

template<>

bool gen::JetMatchingMadgraph::parseParameter ( const std::string & value_ )

staticprotected

Definition at line 89 of file JetMatchingMadgraph.cc.

                                                                   {
     std::string value(value_);
     std::transform(value.begin(), value.end(), value.begin(), (int (*)(int))std::toupper);
     return value == "T" || value == "Y" || value == "True" || value == "1" || value == ".TRUE.";
   }

References createfilelist::int, AlCaHLTBitMon_QueryRunRegistry::string, HcalDetIdTransform::transform(), and relativeConstraints::value.

◆ updateOrDie()

template<typename T >

void gen::JetMatchingMadgraph::updateOrDie	(	const std::map< std::string, std::string > &	params,
		T &	param,
		const std::string &	name
	)

staticprotected

Definition at line 189 of file JetMatchingMadgraph.cc.

                                                                {
     if (param < 0) {
       param = getParameter(params, name, param);
     }
     if (param < 0)
       throw cms::Exception("Generator|PartonShowerVeto") << "The MGParamCMS header does not specify the jet "
                                                             "matching parameter \""
                                                          << name
                                                          << "\", but it "
                                                             "is requested by the CMSSW configuration."
                                                          << std::endl;
   }

References Exception, getParameter(), Skims_PA_cff::name, and CalibrationSummaryClient_cfi::params.

Referenced by init().

Member Data Documentation

◆ eventInitialized

bool gen::JetMatchingMadgraph::eventInitialized

protected

Definition at line 44 of file JetMatchingMadgraph.h.

Referenced by beforeHadronisation(), beforeHadronisationExec(), and match().

◆ exclusive

bool gen::JetMatchingMadgraph::exclusive

protected

Definition at line 46 of file JetMatchingMadgraph.h.

Referenced by JetMatchingMadgraph(), and match().

◆ mgParams

std::map<std::string, std::string> gen::JetMatchingMadgraph::mgParams

protected

Definition at line 41 of file JetMatchingMadgraph.h.

Referenced by getParameter(), and init().

◆ runInitialized

bool gen::JetMatchingMadgraph::runInitialized

protected

Definition at line 43 of file JetMatchingMadgraph.h.

Referenced by beforeHadronisation(), init(), and match().

◆ soup

bool gen::JetMatchingMadgraph::soup

protected

Definition at line 45 of file JetMatchingMadgraph.h.

Referenced by JetMatchingMadgraph(), and match().

Public Member Functions

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ JetMatchingMadgraph()

◆ ~JetMatchingMadgraph()

Member Function Documentation

◆ beforeHadronisation()

◆ beforeHadronisationExec()

◆ capabilities()

◆ getJetEtaMax()

◆ getParameter() [1/2]

◆ getParameter() [2/2]

◆ init()

◆ match()

◆ parseParameter() [1/3]

◆ parseParameter() [2/3]

◆ parseParameter() [3/3]

◆ updateOrDie()

Member Data Documentation

◆ eventInitialized

◆ exclusive

◆ mgParams

◆ runInitialized

◆ soup