PdfInfo for SingleTopDPDMaker

PdfInfo or PdfTreeMaker is a tool that inserts PdfInfo using two different methods for mc08 and mc09.

mc09: uses genEvt->pdf_info()

mc08: as pdf_info not available, fall-back method fills first two partons (usually barcode 3 & 4) from TruthParticleContainer

variables

For the mc09 method five variables will be saved in the PdfTree:

  int m_id1;        // flavour code of first parton
  int m_id2;        // flavour code of second parton 
  double m_x1;      // fraction of beam momentum carried by first parton ("beam side") 
  double m_x2;      // fraction of beam momentum carried by second parton ("target side") 
  double m_scalePDF;  //  Q-scale used in evaluation of PDF's   (in GeV) 

For the mc08 method the kinematics and pdgId of the incoming partons are saved. px, py and pt should be zero. The number of partons should be two.

  std::vector<double>* m_InitialParton_px;
  std::vector<double>* m_InitialParton_py;
  std::vector<double>* m_InitialParton_pz;
  std::vector<double>* m_InitialParton_pt;
  std::vector<double>* m_InitialParton_e;
  std::vector<double>* m_InitialParton_pdgId;
  std::vector<double>* m_InitialParton_barcode;
  int m_InitialParton_N;

implementation

Files:

• PdfTreeMaker.cxx

• PdfTreeMaker.h

The following additions have to be made to the python files in the share directory:

runSingleTopDPDMaker.py:

# Pdf Info dumper

ExecutePdfTreeMaker    = True #write pdf info
PdfMethod = "mc08" #"mc08", "mc09" or "mc08,mc09"
Mc08_barcode = [1,2,3,4]
Mc08_pdgId = [1,2,3,4,5,21]

setSingleTopDPDMaker.py:

# Pdf Info container

PdfTreeMaker_TreeName = "PdfTree"

if not "PdfMethod" in dir() : PdfMethod = ""
if not "Mc08_pdgId" in dir() : Mc08_pdgId = []
if not "Mc08_barcode" in dir() : Mc08_barcode = []

PdfTreeMaker_McEventContainerName = McEventContainerName
PdfTreeMaker_TruthParticleContainerName = TruthParticleContainerName
PdfTreeMaker_PdfMethod = PdfMethod
PdfTreeMaker_Mc08_pdgId = Mc08_pdgId
PdfTreeMaker_Mc08_barcode = Mc08_barcode

exeSingleTopDPDMaker:

# Pdf Info

if ExecutePdfTreeMaker:

  PdfTreeMaker = PdfTreeMaker("PdfTreeMaker")
  Sequence += PdfTreeMaker

  PdfTreeMaker.OutputLevel = INFO

  PdfTreeMaker.TreeName             = PdfTreeMaker_TreeName
  PdfTreeMaker.McEventContainerName = PdfTreeMaker_McEventContainerName
  PdfTreeMaker.TruthParticleContainerName = PdfTreeMaker_TruthParticleContainerName
  PdfTreeMaker.PdfMethod = PdfTreeMaker_PdfMethod
  PdfTreeMaker.Mc08_pdgId = PdfTreeMaker_Mc08_pdgId
  PdfTreeMaker.Mc08_barcode = PdfTreeMaker_Mc08_barcode

CategorySingleTopDPDMaker