gen::JetMatchingAlpgen Class Reference

#include <JetMatchingAlpgen.h>

Inheritance diagram for gen::JetMatchingAlpgen:

Public Member Functions
	JetMatchingAlpgen (const edm::ParameterSet &params)
	~JetMatchingAlpgen () override
Public Member Functions inherited from gen::JetMatching
virtual void	beforeHadronisationExec ()
virtual const std::vector< int > *	getPartonList ()
virtual bool	initAfterBeams ()
bool	isMatchingDone ()
	JetMatching (const edm::ParameterSet &params)
void	resetMatchingStatus ()
virtual	~JetMatching ()

Private Member Functions
void	beforeHadronisation (const lhef::LHEEvent *event) override
std::set< std::string >	capabilities () const override
double	getJetEtaMax () const override
void	init (const lhef::LHERunInfo *runInfo) override
int	match (const lhef::LHEEvent *partonLevel, const std::vector< fastjet::PseudoJet > *jetInput) override

Private Attributes
bool	applyMatching
bool	eventInitialized
AlpgenHeader	header
bool	runInitialized

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

Detailed Description

Definition at line 10 of file JetMatchingAlpgen.h.

Constructor & Destructor Documentation

JetMatchingAlpgen()

gen::JetMatchingAlpgen::JetMatchingAlpgen

(

const edm::ParameterSet &

params

)

Definition at line 34 of file JetMatchingAlpgen.cc.

References ahopts_, AHOPTS::etclus, AHOPTS::iexc, submitPVValidationJobs::params, and AHOPTS::rclus.

      : JetMatching(params), applyMatching(params.getParameter<bool>("applyMatching")), runInitialized(false) {
    ahopts_.etclus = params.getParameter<double>("etMin");
    ahopts_.rclus = params.getParameter<double>("drMin");
    ahopts_.iexc = params.getParameter<bool>("exclusive");
  }

~JetMatchingAlpgen()

gen::JetMatchingAlpgen::~JetMatchingAlpgen ( )

override

Definition at line 42 of file JetMatchingAlpgen.cc.

{}

Member Function Documentation

beforeHadronisation()

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

overrideprivatevirtual

Reimplemented from gen::JetMatching.

Definition at line 113 of file JetMatchingAlpgen.cc.

References eventInitialized, Exception, and runInitialized.

                                                                       {
    // We can't continue if the run has not been initialized.
    if (!runInitialized)
      throw cms::Exception("Generator|PartonShowerVeto") << "Run not initialized in JetMatchingAlpgen" << std::endl;

    // We are called just after LHEInterface has filled in
    // the Fortran common block (and Pythia6 called UPEVNT).

    // Possibly not interesting for us.
    // (except perhaps for debugging?)
    //  pyupre_();
    //  dbpart_();
    eventInitialized = true;
  }

capabilities()

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

overrideprivatevirtual

Reimplemented from gen::JetMatching.

Definition at line 44 of file JetMatchingAlpgen.cc.

References mps_fire::result.

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

getJetEtaMax()

double gen::JetMatchingAlpgen::getJetEtaMax ( ) const

inlineoverrideprivatevirtual

Implements gen::JetMatching.

Definition at line 20 of file JetMatchingAlpgen.h.

{ return 0.; }

init()

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

overrideprivatevirtual

Reimplemented from gen::JetMatching.

Definition at line 54 of file JetMatchingAlpgen.cc.

References ahpars_, ahppara_, gen::alsetp_(), filterCSVwithJSON::copy, AlpgenHeader::ebeam, AHPPARA::ebeam, header, AHPPARA::ihrd, AlpgenHeader::ihrd, createfilelist::int, AlpgenHeader::MASS_MAX, AHPPARA::masses, AlpgenHeader::masses, AHPARS::nparam, AlpgenHeader::params, AlpgenHeader::parse(), AHPARS::parval, submitPVValidationJobs::runInfo, and runInitialized.

                                                            {
    // Read Alpgen run card stored in the LHERunInfo object.
    std::vector<std::string> headerLines = runInfo->findHeader("AlpgenUnwParFile");
    if (headerLines.empty())
      throw cms::Exception("Generator|PartonShowerVeto") << "In order to use Alpgen jet matching, "
                                                            "the input file has to contain the corresponding "
                                                            "Alpgen headers."
                                                         << std::endl;

    // Parse the header using its bultin function.
    header.parse(headerLines.begin(), headerLines.end());

    // I don't want to print this right now.
    //  std::cout << "Alpgen header" << std::endl;
    //  std::cout << "========================" << std::endl;
    //  std::cout << "\tihrd = " << header.ihrd << std::endl;
    //  std::cout << "\tmc = " << header.masses[AlpgenHeader::mc]
    //            << ", mb = " << header.masses[AlpgenHeader::mb]
    //            << ", mt = " << header.masses[AlpgenHeader::mt]
    //            << ", mw = " << header.masses[AlpgenHeader::mw]
    //            << ", mz = " << header.masses[AlpgenHeader::mz]
    //            << ", mh = " << header.masses[AlpgenHeader::mh]
    //            << std::endl;
    //  for(std::map<AlpgenHeader::Parameter, double>::const_iterator iter =
    //      header.params.begin(); iter != header.params.end(); ++iter)
    //      std::cout << "\t" << AlpgenHeader::parameterName(iter->first)
    //                << " = " << iter->second << std::endl;
    //  std::cout << "\txsec = " << header.xsec
    //            << " +-" << header.xsecErr << std::endl;
    //  std::cout << "========================" << std::endl;

    // Here we pass a few header variables to common block and
    // call Alpgen init routine to do the rest.
    // The variables passed first are the ones directly that
    // need to be set up "manually": IHRD and the masses.
    // (ebeam is set just to mimic the original code)
    // Afterwards, we pass the full spectrum of Alpgen
    // parameters directly into the AHPARS structure, to be
    // treated by AHSPAR which is called inside alsetp_().

    std::copy(header.masses, header.masses + AlpgenHeader::MASS_MAX, ahppara_.masses);
    ahppara_.ihrd = header.ihrd;
    ahppara_.ebeam = header.ebeam;

    for (std::map<AlpgenHeader::Parameter, double>::const_iterator iter = header.params.begin();
         iter != header.params.end();
         ++iter) {
      if (iter->first <= 0 || iter->first >= (int)AHPARS::nparam - 1)
        continue;
      ahpars_.parval[(int)iter->first - 1] = iter->second;
    }

    // Run the rest of the setup.
    alsetp_();

    // When we reach this point, the run is fully initialized.
    runInitialized = true;
  }

match()

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

overrideprivatevirtual

Implements gen::JetMatching.

Definition at line 133 of file JetMatchingAlpgen.cc.

References gen::alveto_(), applyMatching, eventInitialized, Exception, runInitialized, and relval_cleanedupgrade::veto.

                                                                                                           {
    /*
  if (!showeredFinalState)
    throw cms::Exception("Generator|PartonShowerVeto")
      << "Alpgen matching expected parton shower "
      "final state." << std::endl;
*/

    if (!runInitialized)
      throw cms::Exception("Generator|PartonShowerVeto") << "Run not initialized in JetMatchingAlpgen" << std::endl;

    if (!eventInitialized)
      throw cms::Exception("Generator|PartonShowerVeto") << "Event not initialized in JetMatchingAlpgen" << std::endl;

    // If matching not required (e.g., icckw = 0), don't run the
    // FORTRAN veto code.
    if (!applyMatching)
      return 0;

    // Call the Fortran veto code.
    int veto = 0;
    alveto_(&veto);

    eventInitialized = false;

    // If event was vetoed, the variable veto will contain the number 1.
    // In this case, we must return 1 - that will be used as the return value from UPVETO.
    // If event was accepted, the variable veto will contain the number 0.
    // In this case, we must return 0 - that will be used as the return value from UPVETO.
    return veto ? 1 : 0;
  }

Member Data Documentation

applyMatching

bool gen::JetMatchingAlpgen::applyMatching

private

Definition at line 28 of file JetMatchingAlpgen.h.

Referenced by match().

eventInitialized

bool gen::JetMatchingAlpgen::eventInitialized

private

Definition at line 30 of file JetMatchingAlpgen.h.

Referenced by beforeHadronisation(), and match().

header

AlpgenHeader gen::JetMatchingAlpgen::header

private

Definition at line 32 of file JetMatchingAlpgen.h.

Referenced by init().

runInitialized

bool gen::JetMatchingAlpgen::runInitialized

private

Definition at line 29 of file JetMatchingAlpgen.h.

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