Before tagging this version, added rest of files that potentially can be useful in the future.

Author: rafopar

AnaTwoPeg/AnaTwoPeg.cc

@@ -0,0 +1 @@
+g++ AnaTwoPeg.cc -o AnaTwoPeg.exe `root-config --cflags --libs`

AnaTwoPeg/compile_AnaTwoPeg.sh

@@ -0,0 +1,252 @@
+
+ Lepton type
+ For DDVCS: muon (2) or no mass (3) IN THE FINAL STATE
+
+ Choose one of the following menu items :
+ 1 : 5-fold DIFFERENTIAL CROSS SECTION for ep -> ep gamma as function of t
+ 2 : e+e- asymmetry for DVCS
+ 3 : electron single spin asymmetry (SSA)
+ 4 : 2-body DOUBLY POLARIZED cross sections for (D)DVCS
+     polarized electron, polarized target
+ 5 : DIFFERENTIAL CROSS SECTION for DVCS as function of Q2
+ 6 : DIFFERENTIAL CROSS SECTION for DVCS as function of t
+ 11 : parton distributions
+ 12 : GPDs
+ 13 : sum rules and positivity tests
+ 14 : form factors from GPDs and WACS
+ 15 : accuracy test of integrals for twist-2 amplitude
+ 16 : accuracy test of integrals for twist-3 amplitude
+ 17 : integration of BH over bin size
+ 18 : fit of kinematics
+ 19 : phase space for (TCS and DVCS) and DDVCS
+ 20 : 3-body DOUBLY POLARIZED cross sections for DDVCS
+     polarized electron, polarized target
+ 21 : 2-body DIFFERENTIAL CROSS SECTION for Rho as function of Q2 or theta_lab
+ 22 : Timelike Compton scattering
+ 30 : WACS
+ 40 : Fit of DVCS data
+ 50 : Transverse densities
+ 60 : Resolutions
+
+ Make a choice for the MECHANISM :
+ 1 : Bethe-Heitler contribution
+ 2 : DVCS contribution
+ 3 : Bethe-Heitler + DVCS contribution
+ 4 : rho contribution
+ 5 : DVCS + rho contribution
+ 6 : Bethe-Heitler + DVCS + rho contribution
+ 7 : Bethe-Heitler + rho contribution
+
+ Choice to calculate on proton or on neutron (default = 1)
+ 1 = proton
+ 2 = neutron
+
+ Choose model for GPD
+ 1 : xi INdependent parametrization with MRS S_o distribution
+ 2 : xi INdependent parametrization with MRS D_o distribution
+ 3 : xi INdependent parametrization with MRS S_o distribution
+ 4 : xi INdependent parametrization with MRS95G distribution
+ 5 : xi INdependent parametrization with MRST98 distribution
+ 6 : xi INdependent parametrization with MRST98 distribution at mu^2 = 1 GeV^2
+ 7 : xi INdependent parametrization with MRST98 distribution (upper gluon, k_perp_av = 0 GeV) at mu^2 = 1 GeV^2
+ 8 : xi INdependent parametrization with CTEQ4L distribution
+ 9 : xi INdependent parametrization with CTEQ4M distribution
+ 10 : xi INdependent parametrization with CTEQ5M distribution
+ 11 : xi INdependent parametrization with GRV98 NLO(MSbar) distribution
+ 12 : xi INdependent parametrization with GRV98 NLO(DIS) distribution
+ 13 : xi INdependent parametrization with GRV98 LO distribution
+ 14 : xi INdependent parametrization with MRST01 NLO distribution
+ 15 : xi INdependent parametrization with MRST02 NNLO distribution at mu^2 = 1 GeV^2
+ 16 : xi INdependent parametrization with MRST02 NNLO distribution
+ 17 : xi INdependent parametrization with CTEQ6M distribution
+ 21 : xi dependent parametrization with MRS S_o distribution 
+ 22 : xi dependent parametrization with MRS D_o distribution 
+ 24 : xi dependent parametrization with MRS95G distribution 
+ 25 : xi dependent parametrization with MRST98 distribution 
+ 26 : xi dependent parametrization with MRST98 distribution at mu^2 = 1 GeV^2
+ 27 : xi dependent parametrization with MRST98 distribution (upper gluon, k_perp_av = 0 GeV) at mu^2 = 1 GeV^2
+ 28 : xi dependent parametrization with CTEQ4L distribution 
+ 29 : xi dependent parametrization with CTEQ4M distribution 
+ 30 : xi dependent parametrization with CTEQ5M distribution 
+ 31 : xi dependent parametrization with GRV98 NLO(MSbar) distribution
+ 32 : xi dependent parametrization with GRV98 NLO(DIS) distribution
+ 33 : xi dependent parametrization with GRV98 LO distribution
+ 34 : xi dependent parametrization with MRST01 NLO distribution
+ 35 : xi dependent parametrization with MRST02 NNLO distribution at mu^2 = 1 GeV^2
+ 36 : xi dependent parametrization with MRST02 NNLO distribution
+ 37 : xi dependent parametrization with CTEQ6M distribution
+ 41 : k_perp dependence with MRS S_o distribution
+ 42 : k_perp dependence with MRS D_o distribution
+ 44 : k_perp dependence with MRS95G distribution
+ 45 : k_perp dependence with MRST98 distribution
+ 46 : k_perp dependence with MRST98 distribution at mu^2 = 1 GeV^2
+ 47 : k_perp dependence with MRST98 distribution (upper gluon, k_erp_av = 0 GeV) at mu^2 = 1 GeV^2
+ 48 : k_perp dependence with CTEQ4L distribution
+ 49 : k_perp dependence with CTEQ4M distribution
+ 50 : k_perp dependence with CTEQ5M distribution
+ 54 : k_perp dependence with MRST01 NLO distribution
+ 55 : k_perp dependence with MRST02 NNLO distribution at mu^2 = 1 GeV^2
+ 56 : k_perp dependence with MRST02 NNLO distribution
+ 57 : k_perp dependence with CTEQ6M distribution
+
+Give the value for the power b in the profile function for the valence contribution to H (e.g. 1.)
+
+Give the value for the power b in the profile function for the sea contribution to H (e.g. 1.)
+
+ Choose the model for the t-dependence of the GPD H
+1 = Factorized model for the t-dependence
+2 = Regge inspired ansatz for the t-dependence
+3 = Experimental exponential input by user
+4 = Experimental exponential fit (fit of all data)
+5 = Experimental exponential fit (fit of subrange -by Mick-)
+6 = Hybrid model (FF vor val and Regge for sea)
+7 = Diehl et al. model (t-dep in DDs)
+8 = R2 Regge ansatz model (t-dep in DDs)
+9 = Diehl et al. model (t-dep out of DDs)
+10 = R2 Regge ansatz model (t-dep out of DDs)
+
+Enter slope alphap (GeV-2)
+
+Do you want to evaluate the D-term contribution to the GPD H?
+1 = Yes
+2 = No
+
+Do you want to evaluate the GPD E?
+1 = only D-term contribution
+2 = double distribution contribution + D-term contribution
+3 = No, E = 0
+
+Do you want to evaluate the pi0 pole contribution (i.e. SPD Etilde)?
+1 = Yes
+2 = No
+
+ Do you want to include twist-3 corrections ?
+1: Do not include twist-3 corrections for L photon,
+   but include the minimal corrections to restore gauge invariance for T photon
+2: Include twist-3 corrections for L photon in Wandzura-Wilczek approximation
+
+ Give the polarization of the target proton
+ 1 : proton polarized along x-axis
+ 2 : proton polarized along y-axis (perpendicular to lepton plane)
+ 3 : proton polarized along z-axis (along the virtual photon direction)
+
+ Calculation for what LEPTON charge ?
+ 1 : negatively charged lepton (JLab)
+
+
+ Give the value of beam energy in GeV (e.g. 27.)
+
+
+ As a function of :
+ 1 : Q^2
+ 2 : x_B
+ 3 : scattering angle (or -t)
+ 4 : Phi_LH
+ 5 : Qprime^2
+ 6 : decay angles
+
+
+ Give the value of Q^2 in GeV^2 (e.g. 5.0)
+
+
+ Give the value of x_B (e.g. 0.3)
+
+
+ Give the value of Qprime^2 in GeV^2 (e.g. 2.0)
+
+
+ Decay angles
+ 0 : Fixed values of c.m. decay angles
+ 1 : Integration over the 2 angles
+
+
+ Theta_cm integration between 0 and PI (0) or between PI/4 and 3PI/4 (1) or user cut (3) ?
+
+
+  ANTI-symmetrisation or not ?
+ 1 : yes
+ 2 : no
+
+
+  Cross section in lab or invariant ?
+ 1 : in lab
+ 2 : invariant
+
+
+ Calculate only im.part (1), only real part (2) or both (3) ?
+
+
+ As a function of t or theta ?
+ 1 : t
+ 2 : theta
+
+
+ Give the first -t to calculate (e.g. .1)
+
+
+ Give the last -t to calculate (e.g. 10.)
+
+
+ Give the step in -t to calculate (e.g. 1.)
+
+
+ Give the out-of-plane angle (in deg) (e.g. 120.)
+
+ 0.000000e+00 	 1.000000e-02 	 nan  	 nan 
+
+ cos(theta) = 0.999790 , ga_en_out = 9.711111 GeV, ga_mom_out = 9.638759 GeV, theta_lab = 1.174746 deg
+ 2.050319e-02 	 1.100000e-01 	 1.891783e-06  	 5.596789e-17 
+
+ cos(theta) = 0.999383 , ga_en_out = 9.657822 GeV, ga_mom_out = 9.585068 GeV, theta_lab = 2.013469 deg
+ 3.514166e-02 	 2.100000e-01 	 8.921187e-07  	 1.780241e-16 
+
+ cos(theta) = 0.998971 , ga_en_out = 9.604532 GeV, ga_mom_out = 9.531371 GeV, theta_lab = 2.599129 deg
+ 4.536335e-02 	 3.100000e-01 	 4.248508e-07  	 3.115185e-17 
+
+ cos(theta) = 0.998556 , ga_en_out = 9.551243 GeV, ga_mom_out = 9.477671 GeV, theta_lab = 3.079723 deg
+ 5.375131e-02 	 4.100000e-01 	 2.033915e-07  	 -3.253550e-17 
+
+ cos(theta) = 0.998136 , ga_en_out = 9.497953 GeV, ga_mom_out = 9.423965 GeV, theta_lab = 3.498879 deg
+ 6.106696e-02 	 5.100000e-01 	 9.734270e-08  	 0.000000e+00 
+
+ cos(theta) = 0.997712 , ga_en_out = 9.444664 GeV, ga_mom_out = 9.370255 GeV, theta_lab = 3.876635 deg
+ 6.766004e-02 	 6.100000e-01 	 4.627972e-08  	 0.000000e+00 
+
+ cos(theta) = 0.997284 , ga_en_out = 9.391375 GeV, ga_mom_out = 9.316540 GeV, theta_lab = 4.224164 deg
+ 7.372558e-02 	 7.100000e-01 	 2.168795e-08  	 0.000000e+00 
+
+ cos(theta) = 0.996851 , ga_en_out = 9.338085 GeV, ga_mom_out = 9.262820 GeV, theta_lab = 4.548452 deg
+ 7.938546e-02 	 8.100000e-01 	 9.915066e-09  	 -4.171332e-17 
+
+ cos(theta) = 0.996413 , ga_en_out = 9.284796 GeV, ga_mom_out = 9.209095 GeV, theta_lab = 4.854206 deg
+ 8.472188e-02 	 9.100000e-01 	 4.356159e-09  	 -4.747190e-17 
+
+ cos(theta) = 0.995971 , ga_en_out = 9.231506 GeV, ga_mom_out = 9.155365 GeV, theta_lab = 5.144778 deg
+ 8.979332e-02 	 1.010000e+00 	 1.795639e-09  	 0.000000e+00 
+
+ cos(theta) = 0.995525 , ga_en_out = 9.178217 GeV, ga_mom_out = 9.101630 GeV, theta_lab = 5.422655 deg
+ 9.464318e-02 	 1.110000e+00 	 6.650238e-10  	 -7.773976e-17 
+
+ cos(theta) = 0.995073 , ga_en_out = 9.124927 GeV, ga_mom_out = 9.047889 GeV, theta_lab = 5.689740 deg
+ 9.930469e-02 	 1.210000e+00 	 2.016206e-10  	 -3.205202e-17 
+
+ cos(theta) = 0.994617 , ga_en_out = 9.071638 GeV, ga_mom_out = 8.994143 GeV, theta_lab = 5.947530 deg
+ 1.038040e-01 	 1.310000e+00 	 3.803139e-11  	 0.000000e+00 
+
+ cos(theta) = 0.994156 , ga_en_out = 9.018348 GeV, ga_mom_out = 8.940392 GeV, theta_lab = 6.197227 deg
+ 1.081620e-01 	 1.410000e+00 	 4.451147e-13  	 1.134249e-16 
+
+ cos(theta) = 0.993690 , ga_en_out = 8.965059 GeV, ga_mom_out = 8.886635 GeV, theta_lab = 6.439813 deg
+ 1.123959e-01 	 1.510000e+00 	 9.181520e-12  	 0.000000e+00 
+
+ cos(theta) = 0.993219 , ga_en_out = 8.911770 GeV, ga_mom_out = 8.832873 GeV, theta_lab = 6.676105 deg
+ 1.165200e-01 	 1.610000e+00 	 2.967123e-11  	 5.444961e-17 
+
+ cos(theta) = 0.992743 , ga_en_out = 8.858480 GeV, ga_mom_out = 8.779104 GeV, theta_lab = 6.906788 deg
+ 1.205462e-01 	 1.710000e+00 	 4.841967e-11  	 0.000000e+00 
+
+ cos(theta) = 0.992262 , ga_en_out = 8.805191 GeV, ga_mom_out = 8.725330 GeV, theta_lab = 7.132447 deg
+ 1.244847e-01 	 1.810000e+00 	 6.135448e-11  	 0.000000e+00 
+
+ cos(theta) = 0.991775 , ga_en_out = 8.751901 GeV, ga_mom_out = 8.671550 GeV, theta_lab = 7.353583 deg
+ 1.283442e-01 	 1.910000e+00 	 6.829013e-11  	 -4.731539e-17