Dilepton

From Hall A Wiki
Revision as of 13:30, 17 April 2015 by Zwzhao (Talk | contribs) (code)

Jump to: navigation, search

Introduction

program

SoLID_JPsi

TCS

DDVCS

emaillist

ee@jlab.org

register at https://mailman.jlab.org/mailman/listinfo/ee

code

It is a root file with 2000000 events. You have 4-vectors that correspond to the outgoing particle energies and momenta, 
kinematical variables: Phi_CMV and Theta_CMV are the gamma* CM angles, Qp2 is Q'^2, 
Egamma is energy of the beam and tt is t. ttmin is t_min. 
You can select "VirtualFlag==0" correspond to event generated from a photon beam 
and "VirtualFlag==1" from quasi-real photon, with a cut maximal Q2=0.3 GeV2.

To weight events you should use "W_tot_unpol". I am not sure yet of the polarized weights in some of the bins, 
so please let me check it before using it. 
You can nevertheless also use "W_BH" and "W_TCS" that correspond to "BH only" or "TCS only" contribution to the unpolarized cross sections.

Polarization vectors are probably not good in this version, please don't use it. The code is still in construction... 

I am doing some checks with the weights of the event generator. It seems there are some shifts when attributing  weights in the non physical part of the phase space. It comes when I generate the grid, I removed all "nan" in calculations but I introduced this problem. I will solve it today but it means I have to run a new table for cross sections. It is around 3 days.
 
  • LPAIR
    • used by the experiments at DESY/HERA
  • Compton2.0
    • BH from photoproduction, used by ZEUS
(It may work on neutron also)
GiBUU does a good job with real photons on a proton (or heavier)
target.  It doesn't handle electroproduction as nicely (at least not the
last time I asked the Giessen group).
It is fairly straight forward to generate events, and there USED to be
good examples of input jobcards on the website, but they have been removed
recently while they upgrade.  I'd think the best route is to email them
directly to get the code and also ask for the example jobcards.