####################################################################################
#####                                2gNNLO                                   ######       
####################################################################################
#####      To select the PDFs go to ./src/pdf.f      (results in pb)          ######
#####      or go to                 ./src/pdf_lha.f  (results in pb)          ######
####################################################################################
# cm energy (in GeV)################################################################
####################################################################################
ecm 7000.0

####################################################################################
# kinematic cuts   #################################################################
####################################################################################
# photon pt (in GeV)################################################################
ph1_ptmin 25.0
ph2_ptmin 22.0
# photon rapidity###################################################################
# It is possible to remove the crack from user\_2fot.cc file                       #
# or simply put a negative number in crack_eta_min and crack_eta_min               #
####################################################################################
ph_eta_max 2.37 
crack_eta_min  -1.44
crack_eta_max  -1.57
# dR separation between the photons#################################################
ph_deltaR_min 0.4
####################################################################################
# invariant mass cuts  #############################################################
####################################################################################
# mgg_min has to be safely set in the perturbative region 
# mgg_max can be even larger than 1TeV but calculation time will be increased
####################################################################################
# Higgs
####################################################################################
# mgg_min 4.0
# mgg_max 800.0
####################################################################################
# SM diphoton
####################################################################################
mgg_min 1.0
mgg_max 800.0

####################################################################################
####################################################################################
# Write “Scale Mff” for Mgamgam as central scale, 
# Write “Scale MH” for Mgamgam=125GeV, as central scale, 
# or “Scale MT” for Mgamgam = sqrt(Mff^2+pTff^2), as central scale, 
####################################################################################
####################################################################################
Scale Mff
####################################################################################
muR 1.0
muF 1.0
####################################################################################

####################################################################################
# isolation cuts: R0 and E0 (in GeV) ###############################################
####################################################################################
# ph_iso_E0 is the fraction of the transverse momentum of the photon
# or the E_t_max parameter. This option is available in the user\_2fot.cc
# default: ph_iso_E0 is the E_t_max parameter
#isolType : Frixione or Hybrid
####################################################################################
ph_iso_R0 0.4 
ph_iso_E0 10.0
ph_iso_n 1.0
ph_iso_Ed 10.0
ph_iso_Rd 0.1
isolType Frixione

####################################################################################
# output option    on/off (1/0)      ###############################################
####################################################################################

# save or not the histograms in ROOT format
SaveRootHistos 1
# save or not the ntuple with Mgg and the weight
SaveRootNtuple 0
# save or not the histograms in TopDrawer format
SaveTopdHistos 1

####################################################################################
#histograms#########################################################################
#format: hist [id] [var] [nbins] [min] [max] [LIN/LOG]
# For more variables go to user\_2fot.cc
# pTgg : Transverse momentum of the diphoton pair
# Mgg  : Invariant mass of the diphoton pair
# dphi : Angular separation between the photons
# CosP : Cos_Theta_*
# CosA : Approximated Cos_Theta_* used in some Tevatron analyses
# pT1  : pT of the hardest photon
# pT2  : pT of the Softer photon
####################################################################################
hist 1 Mgg 100 0.0 200.0 LIN
hist 2 CosP 25 -1.0 1.0 LIN
hist 3 pTgg 100 0.0 400.0 LIN
hist 4 dphi 3 0. 1.5 LIN
hist 5 dphi 3 1.5 2.25 LIN
hist 6 dphi 4 2.25 2.65 LIN
hist 7 dphi 9 2.65 3.1 LIN
hist 8 dphi 1 3.1 3.15 LIN
       
####################################################################################
# 2gNNLO settings ##################################################################
####################################################################################
# contributions to switch on/off (1/0)
####################################################################################
## BORN ############################################################################
iifborn 0
## NLO  ############################################################################
iifnnlo 1
## NNLO ############################################################################
iinnlo 0
####################################################################################
# run option #######################################################################
####################################################################################
# 1: START
# 2: CONTINUE
# 3: RESTART WITH PREVIOUS GRIDS
####################################################################################
i0 1
ifnn 1
inn 1
####################################################################################
# Number of iterations  ############################################################
####################################################################################
n0 -200
fnn0 -300
nn0 -550
####################################################################################
# Number of calls per iteration  ###################################################
####################################################################################
ncl0 190000
fnncl 640000
nncl 2000000