Page MenuHomeHEPForge
Files F8723990

mcfmw_pdf.h
No OneTemporary
Actions

mcfmw_pdf.h
View Options

/** emacs: this is -*- c++ -*- **/
/**
** @file mcfmw_pdf.h
**
** @brief pdf transform functions for W production
**
** @author mark sutton
** @date Mon Dec 10 01:36:04 GMT 2007
**
** @copyright (C) 2002-2019 mark sutton (sutt @ cern.ch)
**
** $Id: mcfmw_pdf.h, v1.0 Mon Dec 10 01:36:04 GMT 2007 sutt $
**
**/
#ifndef __MCFMW_PDF_H
#define __MCFMW_PDF_H
#include "appl_grid/appl_pdf.h"
//
// MCFM W production
//
//
// W+
//
class mcfmwp_pdf : public appl::appl_pdf {
public:
mcfmwp_pdf(const std::string& s="mcfm-wp") : appl_pdf(s) {
m_Nproc = 6;
make_ckm( true );
}
~mcfmwp_pdf() { }
virtual void evaluate(const double* fA, const double* fB, double* H) const;
};
//
// W-
//
class mcfmwm_pdf : public appl::appl_pdf {
public:
mcfmwm_pdf(const std::string& s="mcfm-wm") : appl_pdf(s) {
m_Nproc = 6;
make_ckm( false );
}
~mcfmwm_pdf() { }
virtual void evaluate(const double* fA, const double* fB, double* H) const;
};
// actual funtion to evaluate the pdf combinations
// for the W+
inline void mcfmwp_pdf::evaluate(const double* fA, const double* fB, double* H) const {
const int nQuark = 6;
const int iQuark = 5;
// offset so we can use fA[-6]
// fA += 6;
// fB += 6;
double GA=fA[6];
double GB=fB[6];
double QA=0; double QB=0; double QbA=0; double QbB=0;
for(int i = 1; i <= iQuark; i++)
{
QA += fA[nQuark + i]*m_ckmsum[nQuark + i];
QB += fB[nQuark + i]*m_ckmsum[nQuark + i];
}
for(int i = -iQuark; i < 0; i++)
{
QbA += fA[nQuark + i]*m_ckmsum[nQuark + i];
QbB += fB[nQuark + i]*m_ckmsum[nQuark + i];
}
H[2]=QbA * GB;
H[3]= QA * GB;
H[4]= GA * QbB;
H[5]= GA * QB;
// have to zero H[0..1] first
H[0]=H[1]=0;
for (int i1 = 3; i1 <= 5; i1 += 2)
{
for(int i2 = 8; i2 <= 10; i2 += 2)
{
H[0] += fA[i1]*fB[i2]*m_ckm2[i1][i2];
H[1] += fA[i2]*fB[i1]*m_ckm2[i2][i1];
}
}
}
//
// actual funtion to evaluate the pdf combinations
// for the W-
//
inline void mcfmwm_pdf::evaluate(const double* fA, const double* fB, double* H) const {
const int nQuark = 6;
const int iQuark = 5;
// offset so we can use fA[-6]
// fA += 6;
// fB += 6;
double GA=fA[6];
double GB=fB[6];
double QA=0; double QB=0; double QbA=0; double QbB=0;
for(int i = 1; i <= iQuark; i++)
{
QA += fA[nQuark + i]*m_ckmsum[nQuark + i];
QB += fB[nQuark + i]*m_ckmsum[nQuark + i];
}
for(int i = -iQuark; i < 0; i++)
{
QbA += fA[nQuark + i]*m_ckmsum[nQuark + i];
QbB += fB[nQuark + i]*m_ckmsum[nQuark + i];
}
H[2]= QA * GB;
H[3]= QbA * GB;
H[4]= GA * QB;
H[5]= GA * QbB;
// have to zero H[0..1] first
H[0]=H[1]=0;
for (int i1 = 7; i1 <= 9; i1 += 2)
{
for(int i2 = 2; i2 <= 4; i2 += 2)
{
H[0] += fA[i1]*fB[i2]*m_ckm2[i1][i2];
H[1] += fA[i2]*fB[i1]*m_ckm2[i2][i1];
}
}
}
// fortran callable wrappers
extern "C" void fmcfmwp_pdf__(const double* fA, const double* fB, double* H);
extern "C" void fmcfmwm_pdf__(const double* fA, const double* fB, double* H);
#endif // __MCFMW_PDF_H

File Metadata

Mime Type
text/x-c++
Expires
Mon, Jan 20, 10:13 PM (1 d, 10 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4242594
Default Alt Text
mcfmw_pdf.h (3 KB)

Event Timeline