/usr/include/HepMC/PdfInfo.h is in libhepmc-dev 2.06.09-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 | //--------------------------------------------------------------------------
#ifndef HEPMC_PDF_INFO_H
#define HEPMC_PDF_INFO_H
//////////////////////////////////////////////////////////////////////////
// garren@fnal.gov, July 2006
//
// Additional PDF information
//////////////////////////////////////////////////////////////////////////
namespace HepMC {
//! The PdfInfo class stores PDF information
///
/// \class PdfInfo
/// HepMC::PdfInfo stores additional PDF information for a GenEvent.
/// Creation and use of this information is optional.
///
/// - int id1; // flavour code of first parton
/// - int id2; // flavour code of second parton
/// - int pdf_id1; // LHAPDF set id of first parton (zero by default)
/// - int pdf_id2; // LHAPDF set id of second parton (zero by default)
/// - double x1; // fraction of beam momentum carried by first parton ("beam side")
/// - double x2; // fraction of beam momentum carried by second parton ("target side")
/// - double scalePDF; // Q-scale used in evaluation of PDF's (in GeV)
/// - double pdf1; // PDF (id1, x1, Q)
/// - double pdf2; // PDF (id2, x2, Q)
///
/// Input parton flavour codes id1 & id2 are expected to obey the
/// PDG code conventions, especially g = 21.
///
/// The contents of pdf1 and pdf2 are expected to be x*f(x).
/// The LHAPDF set ids are the entries in the first column of
/// http:///projects.hepforge.org/lhapdf/PDFsets.index
///
class PdfInfo {
public:
// --- birth/death:
//
/// default constructor
PdfInfo()
: m_id1(0),
m_id2(0),
m_pdf_id1(0),
m_pdf_id2(0),
m_x1(0),
m_x2(0),
m_scalePDF(0),
m_pdf1(0),
m_pdf2(0)
{}
/// all values EXCEPT pdf_id1 and pdf_id2 must be provided
PdfInfo( int i1, int i2, double x1, double x2,
double q, double p1, double p2,
int pdf_id1 = 0, int pdf_id2 = 0 );
~PdfInfo() {}
// --- copying:
//
PdfInfo( PdfInfo const & orig ); //!< copy constructor
PdfInfo & operator = ( PdfInfo const & rhs ); //!< make a copy
void swap( PdfInfo & other ); //!< swap two PdfInfo objects
// --- equivalence:
//
bool operator==( const PdfInfo& ) const; //!< check for equality
bool operator!=( const PdfInfo& ) const; //!< check for inequality
// --- accessors:
/// flavour code of first parton
int id1() const { return m_id1; }
/// flavour code of second parton
int id2() const { return m_id2; }
/// LHAPDF set id of first parton
int pdf_id1() const { return m_pdf_id1; }
/// LHAPDF set id of second parton
int pdf_id2() const { return m_pdf_id2; }
/// fraction of beam momentum carried by first parton ("beam side")
double x1() const { return m_x1; }
/// fraction of beam momentum carried by second parton ("target side")
double x2() const { return m_x2; }
/// Q-scale used in evaluation of PDF's (in GeV)
double scalePDF() const { return m_scalePDF; }
/// PDF (id1, x1, Q) - x*f(x)
double pdf1() const { return m_pdf1; }
/// PDF (id2, x2, Q) - x*f(x)
double pdf2() const { return m_pdf2; }
/// verify that the instance contains non-zero information
bool is_valid() const;
// --- mutators:
/// set flavour code of first parton
void set_id1(const int &i) { m_id1=i; }
/// set flavour code of second parton
void set_id2(const int &i) { m_id2=i; }
/// set LHAPDF set id of first parton
void set_pdf_id1(const int &i) { m_pdf_id1=i; }
/// set LHAPDF set id of second parton
void set_pdf_id2(const int &i) { m_pdf_id2=i; }
/// set fraction of beam momentum carried by first parton ("beam side")
void set_x1(const double &f) { m_x1=f; }
/// set fraction of beam momentum carried by second parton ("target side")
void set_x2(const double &f) { m_x2=f; }
/// set Q-scale used in evaluation of PDF's (in GeV)
void set_scalePDF(const double &f) { m_scalePDF=f; }
/// set x*f(x) of first parton
void set_pdf1(const double &f) { m_pdf1=f; }
/// set x*f(x) of second parton
void set_pdf2(const double &f) { m_pdf2=f; }
private: // data members
int m_id1;
int m_id2;
int m_pdf_id1;
int m_pdf_id2;
double m_x1;
double m_x2;
double m_scalePDF;
double m_pdf1;
double m_pdf2;
};
// Free Functions
// IO
std::ostream & operator << (std::ostream &, PdfInfo const *);
std::istream & operator >> (std::istream &, PdfInfo *);
// inline operators
inline PdfInfo::PdfInfo( int i1, int i2, double xin1, double xin2,
double q, double p1, double p2,
int pid1, int pid2 )
: m_id1(i1),
m_id2(i2),
m_pdf_id1(pid1),
m_pdf_id2(pid2),
m_x1(xin1),
m_x2(xin2),
m_scalePDF(q),
m_pdf1(p1),
m_pdf2(p2)
{}
inline PdfInfo::PdfInfo( PdfInfo const & orig )
: m_id1(orig.m_id1),
m_id2(orig.m_id2),
m_pdf_id1(orig.m_pdf_id1),
m_pdf_id2(orig.m_pdf_id2),
m_x1(orig.m_x1),
m_x2(orig.m_x2),
m_scalePDF(orig.m_scalePDF),
m_pdf1(orig.m_pdf1),
m_pdf2(orig.m_pdf2)
{}
inline PdfInfo & PdfInfo::operator = ( PdfInfo const & rhs )
{
PdfInfo temp( rhs );
swap( temp );
return *this;
}
inline void PdfInfo::swap( PdfInfo & other )
{
std::swap(m_id1, other.m_id1);
std::swap(m_id2, other.m_id2);
std::swap(m_pdf_id1, other.m_pdf_id1);
std::swap(m_pdf_id2, other.m_pdf_id2);
std::swap(m_x1, other.m_x1);
std::swap(m_x2, other.m_x2);
std::swap(m_scalePDF, other.m_scalePDF);
std::swap(m_pdf1, other.m_pdf1);
std::swap(m_pdf2, other.m_pdf2);
}
inline bool PdfInfo::operator==( const PdfInfo& a ) const
{
/// equality requires that each member match
return ( a.id1() == this->id1()
&& a.id2() == this->id2()
&& a.pdf_id1() == this->pdf_id1()
&& a.pdf_id2() == this->pdf_id2()
&& a.x1() == this->x1()
&& a.x2() == this->x2()
&& a.scalePDF() == this->scalePDF()
&& a.pdf1() == this->pdf1()
&& a.pdf2() == this->pdf2() );
}
inline bool PdfInfo::operator!=( const PdfInfo& a ) const
{
/// any nonmatching member generates inequality
return !( a == *this );
}
inline bool PdfInfo::is_valid() const
{
if( m_id1 != 0 ) return true;
if( m_id2 != 0 ) return true;
if( m_pdf_id1 != 0 ) return true;
if( m_pdf_id2 != 0 ) return true;
if( m_x1 != 0 ) return true;
if( m_x2 != 0 ) return true;
if( m_scalePDF != 0 ) return true;
if( m_pdf1 != 0 ) return true;
if( m_pdf2 != 0 ) return true;
return false;
}
} // HepMC
#endif // HEPMC_PDF_INFO_H
|