gen::JetMatchingAlpgen Class Reference

#include <JetMatchingAlpgen.h>

Inheritance diagram for gen::JetMatchingAlpgen:

Public Member Functions
	JetMatchingAlpgen (const edm::ParameterSet &params)
	~JetMatchingAlpgen ()
Public Member Functions inherited from gen::JetMatching
virtual void	beforeHadronisationExec ()
const std::vector < JetPartonMatch > &	getMatchSummary () const
bool	isMatchingDone ()
	JetMatching (const edm::ParameterSet &params)
void	resetMatchingStatus ()
virtual	~JetMatching ()

Private Member Functions
void	beforeHadronisation (const lhef::LHEEvent *event)
std::set< std::string >	capabilities () const
void	init (const lhef::LHERunInfo *runInfo)
int	match (const HepMC::GenEvent *partonLevel, const HepMC::GenEvent *finalState, bool showeredFinalState)

Private Attributes
bool	applyMatching
bool	eventInitialized
AlpgenHeader	header
bool	runInitialized

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

Detailed Description

Definition at line 12 of file JetMatchingAlpgen.h.

Constructor & Destructor Documentation

gen::JetMatchingAlpgen::JetMatchingAlpgen

(

const edm::ParameterSet &

params

)

Definition at line 34 of file JetMatchingAlpgen.cc.

References ahopts_, AHOPTS::etclus, edm::ParameterSet::getParameter(), AHOPTS::iexc, 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");
}

gen::JetMatchingAlpgen::~JetMatchingAlpgen

(

)

Definition at line 45 of file JetMatchingAlpgen.cc.

{
}

Member Function Documentation

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

privatevirtual

Reimplemented from gen::JetMatching.

Definition at line 120 of file JetMatchingAlpgen.cc.

References eventInitialized, edm::hlt::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;
}

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

privatevirtual

Reimplemented from gen::JetMatching.

Definition at line 49 of file JetMatchingAlpgen.cc.

References query::result.

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

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

privatevirtual

Reimplemented from gen::JetMatching.

Definition at line 60 of file JetMatchingAlpgen.cc.

References ahpars_, ahppara_, gen::alsetp_(), filterCSVwithJSON::copy, AlpgenHeader::ebeam, AHPPARA::ebeam, lhef::LHERunInfo::findHeader(), header, AHPPARA::ihrd, AlpgenHeader::ihrd, AlpgenHeader::MASS_MAX, AHPPARA::masses, AlpgenHeader::masses, AHPARS::nparam, AlpgenHeader::params, AlpgenHeader::parse(), AHPARS::parval, 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;
}

int gen::JetMatchingAlpgen::match ( const HepMC::GenEvent *  partonLevel, const HepMC::GenEvent *  finalState, bool  showeredFinalState  )

privatevirtual

Implements gen::JetMatching.

Definition at line 138 of file JetMatchingAlpgen.cc.

References gen::alveto_(), applyMatching, eventInitialized, edm::hlt::Exception, and runInitialized.

{
  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

bool gen::JetMatchingAlpgen::applyMatching

private

Definition at line 27 of file JetMatchingAlpgen.h.

Referenced by match().

bool gen::JetMatchingAlpgen::eventInitialized

private

Definition at line 29 of file JetMatchingAlpgen.h.

Referenced by beforeHadronisation(), and match().

AlpgenHeader gen::JetMatchingAlpgen::header

private

Definition at line 31 of file JetMatchingAlpgen.h.

Referenced by init().

bool gen::JetMatchingAlpgen::runInitialized

private

Definition at line 28 of file JetMatchingAlpgen.h.

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