gen::JetMatchingMadgraph Class Reference

#include <JetMatchingMadgraph.h>

Inheritance diagram for gen::JetMatchingMadgraph:

Public Member Functions
	JetMatchingMadgraph (const edm::ParameterSet &params)
template<>
std::string	parseParameter (const std::string &value)
template<>
bool	parseParameter (const std::string &value_)
	~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<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 112 of file JetMatchingMadgraph.cc.

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_, submitPVValidationJobs::params, gen::MEMAIN::qcut, gen::MEMAIN::rclmax, gen::MEMAIN::showerkt, soup, contentValuesCheck::ss, and AlCaHLTBitMon_QueryRunRegistry::string.

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

~JetMatchingMadgraph()

gen::JetMatchingMadgraph::~JetMatchingMadgraph ( )

override

Definition at line 153 of file JetMatchingMadgraph.cc.

{}

Member Function Documentation

beforeHadronisation()

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

overrideprotectedvirtual

Reimplemented from gen::JetMatching.

Definition at line 261 of file JetMatchingMadgraph.cc.

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

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

beforeHadronisationExec()

void gen::JetMatchingMadgraph::beforeHadronisationExec ( )

overrideprotectedvirtual

Reimplemented from gen::JetMatching.

Definition at line 295 of file JetMatchingMadgraph.cc.

References eventInitialized, and gen::mgevnt_().

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

capabilities()

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

overrideprotectedvirtual

Reimplemented from gen::JetMatching.

Definition at line 157 of file JetMatchingMadgraph.cc.

References mps_fire::result.

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

getJetEtaMax()

double gen::JetMatchingMadgraph::getJetEtaMax ( ) const

overrideprotectedvirtual

Implements gen::JetMatching.

Definition at line 155 of file JetMatchingMadgraph.cc.

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

{ return memain_.etaclmax; }

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 98 of file JetMatchingMadgraph.cc.

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

Referenced by getParameter(), and updateOrDie().

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

getParameter() [2/2]

template<typename T >

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

protected

Definition at line 108 of file JetMatchingMadgraph.cc.

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

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

init()

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

overrideprotectedvirtual

Reimplemented from gen::JetMatching.

Definition at line 208 of file JetMatchingMadgraph.cc.

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

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

match()

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

overrideprotectedvirtual

Implements gen::JetMatching.

Definition at line 306 of file JetMatchingMadgraph.cc.

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.

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

parseParameter() [1/3]

template<typename T >

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

staticprotected

Definition at line 73 of file JetMatchingMadgraph.cc.

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

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

parseParameter() [2/3]

template<>

std::string gen::JetMatchingMadgraph::parseParameter

(

const std::string &

value

)

Definition at line 81 of file JetMatchingMadgraph.cc.

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

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

parseParameter() [3/3]

template<>

bool gen::JetMatchingMadgraph::parseParameter

(

const std::string &

value_

)

Definition at line 91 of file JetMatchingMadgraph.cc.

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

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

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 191 of file JetMatchingMadgraph.cc.

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

Referenced by init().

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

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