Difference between revisions of "Dilepton"

From Hall A Wiki
Jump to: navigation, search
(→‎code)
Line 28: Line 28:
  
 
* generator by Marie Boer and Michel GUIDAL
 
* generator by Marie Boer and Michel GUIDAL
−
** input file [https://userweb.jlab.org/~zwzhao/TCS/sim/TCS.input]
+
** the generator code is at /work/halla/solid/mboer/Generator
−
** root file [https://ipnshare.in2p3.fr/owncloud/public.php?service=files&t=2224e3a32537708e35f2082e15a1be8f] (700MB on cloud)
+
*** see "Readme" file there
−
** note from Marie
+
*** The code is still in construction, some values in the Tree may not be filled. Please, follow the instructions in the Readme or use only the kinematical variables + the 4-vectors + W_tot_unpol or cross_tot_unpol to weight the events.
−
<pre>
+
** output root file from the generator  /u/scratch/mboer/
−
It is a root file with 2000000 events. You have 4-vectors that correspond to the outgoing particle energies and momenta,
+
** analysis code to add acceptance https://jlabsvn.jlab.org/svnroot/solid/evgen/dilepton/
−
kinematical variables: Phi_CMV and Theta_CMV are the gamma* CM angles, Qp2 is Q'^2,  
+
*** see "readme" file there
−
Egamma is energy of the beam and tt is t. ttmin is t_min.
+
** output root file with acceptance from the analysis code /u/scratch/zwzhao
−
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.
+
−
</pre>
+
  
 
* genTCS
 
* genTCS

Revision as of 20:24, 5 May 2015

Introduction

meeting

Date Every Day 11:30AM
Location F227 in general, except 5/5 in F224 and 5/11 in L207
Conference ID 	392477255
Desktop sharing and video https://bluejeans.com/392477255
(You don't need to call in by phone if your PC has a mic and a speaker, 
but do REMEMBER to turn off mic and speaker to avoid echo if you use phoneline for audio)
Audio Dial: 18882402560 (US or Canada only) (see all numbers) and enter the Conference ID
meeting system "bluejeans" howto https://hallaweb.jlab.org/wiki/index.php/Bluejeans_howto

program

SoLID_JPsi

TCS

DDVCS

emaillist

ee@jlab.org

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

code

  • generator by Marie Boer and Michel GUIDAL
    • the generator code is at /work/halla/solid/mboer/Generator
      • see "Readme" file there
      • The code is still in construction, some values in the Tree may not be filled. Please, follow the instructions in the Readme or use only the kinematical variables + the 4-vectors + W_tot_unpol or cross_tot_unpol to weight the events.
    • output root file from the generator /u/scratch/mboer/
    • analysis code to add acceptance https://jlabsvn.jlab.org/svnroot/solid/evgen/dilepton/
      • see "readme" file there
    • output root file with acceptance from the analysis code /u/scratch/zwzhao
  • 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.
Retrieved from "https://hallaweb.jlab.org/wiki/index.php?title=Dilepton&oldid=25226"