PROCESS  GG
  collider = LHC
  sqrts = 13000.0
  photon_isolation = frix[R=0.1, ETmax_epsilon=0.15, n=1] + etsum[R=0.2, ETmax_epsilon=0.09]
  decay_type = 0
END_PROCESS

RUN  ATLAS_GG 
  PDF = NNPDF30_nnlo_as_0118[0]
  tcut = 1d-9
  reset_vegas_grid = .true.
  angular_average  = .false.
  multi_channel    = .false.
  lips_reduce      = .true.
  warmup = 10000[10]
  production = 10000[10]
END_RUN

PARAMETERS
  SCHEME[GV]  = CMS
  SCHEME[alpha] = Gmu 
  hard_photon_alpha0 = .true.
END_PARAMETERS

SELECTORS
    select dr_g1g2 min = 0.4 ! Angular separations for photons
    select ptg1 min = 40
    select ptg2 min = 30 ! Leading- and sub-leading-pt cuts for gammas.
    select photons_abs_y max = 2.37 
    reject photons_abs_y min = 1.37 max = 1.52 ! The upper bound for nontrivial cut
END_SELECTORS

HISTOGRAMS

	ptg1 > ph1_pt              [40, 43.6, 47.2, 51, 55, 59.2, 63.6, 68.2, 72.8, 77.6, 82.5, 88, 93.5, 99, 105, 111, 117.5, 124, 131, 138, 145.5, 153, 161, 169, 177.5, 186.5, 195.5, 205, 215, 225, 236, 247, 259, 270, 283, 296, 309, 323, 337, 352, 367, 383, 400, 417, 438, 466, 502, 550, 13000] output_type = 1 
	ptg2 > ph2_pt              [30, 33, 36.4, 39.8, 43.4, 47, 50.8, 54.8, 59, 63.4, 68, 72.6, 77.4, 82.5, 87.5, 93, 99, 104.5, 111, 117, 124, 130.5, 137.5, 145, 153, 160.5, 169, 177.5, 186, 195.5, 205, 214.5, 225, 235, 246, 260, 276, 300, 334, 13000] output_type = 1 
	Mg1g2 > yy_m             [13.5, 20, 22.6, 25, 27, 29.1, 31, 33.2, 35.4, 37.8, 40.4, 43, 46, 49, 52.4, 55.8, 59.6, 63.4, 67.6, 72.2, 77, 82, 87.5, 93, 99, 105, 112, 118.5, 126, 133.5, 141, 149.5, 158, 166.5, 175.5, 185, 194.5, 205, 215, 226, 238, 249, 262, 274, 287, 301, 315, 330, 345, 361, 377, 394, 411, 429, 448, 467, 487, 508, 530, 552, 574, 598, 622, 648, 674, 700, 732, 772, 825, 900, 1020, 13000] output_type = 1 
	abs_cos_thetaCS > yy_cosTS    nbins = 100     min = 0d0       max = 1d0 output_type = 1 

	pt_g1g2 > yy_pt           [1, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30, 33.2, 36.4, 39.8, 43.4, 47, 50.8, 55, 59.2, 63.4, 68, 72.8, 77.8, 83, 88.5, 94, 99.5, 106, 112, 118.5, 125.5, 132.5, 140, 147.5, 155.5, 164, 172.5, 181.5, 191, 201, 211, 222, 233, 244, 256, 268, 280, 293, 306, 320, 333, 348, 362, 379, 398, 424, 456, 500, 13000] output_type = 1 
	phi_star > yy_phiStar          [0.00138, 0.0028, 0.0043, 0.0059, 0.0074, 0.0089, 0.0105, 0.012, 0.0136, 0.0153, 0.0171, 0.0192, 0.0214, 0.024, 0.027, 0.0302, 0.034, 0.038, 0.0426, 0.048, 0.0535, 0.06, 0.067, 0.0755, 0.0845, 0.095, 0.106, 0.119, 0.133, 0.149, 0.167, 0.188, 0.21, 0.236, 0.264, 0.296, 0.332, 0.372, 0.418, 0.47, 0.525, 0.585, 0.66, 0.74, 0.825, 0.925, 1.04, 1.16, 1.31, 1.46, 1.64, 1.84, 2.06, 2.3, 2.58, 2.9, 3.26, 3.64, 4.08, 4.58, 5.15, 5.75, 6.45, 7.35, 8.6, 10.4, 13, 18, 30, 100000000] output_type = 1 
	pi_dphi_g1g2 > yy_phiMDphi      [0.0039, 0.0078, 0.0118, 0.0156, 0.0196, 0.0234, 0.0274, 0.0314, 0.0352, 0.0392, 0.043, 0.047, 0.051, 0.0548, 0.059, 0.0636, 0.0684, 0.0735, 0.0795, 0.0855, 0.092, 0.099, 0.1065, 0.1145, 0.1235, 0.1325, 0.143, 0.154, 0.165, 0.178, 0.192, 0.206, 0.222, 0.239, 0.257, 0.277, 0.298, 0.32, 0.346, 0.372, 0.4, 0.43, 0.464, 0.5, 0.538, 0.578, 0.622, 0.67, 0.72, 0.775, 0.835, 0.9, 0.97, 1.04, 1.12, 1.205, 1.3, 1.4, 1.51, 1.62, 1.74, 1.88, 2.02, 2.17, 2.34, 2.52, 2.71, 2.92, 3.1416] output_type = 1 
	pTt_g1g2 > yy_at           [0.19, 0.4, 0.62, 0.84, 1.08, 1.32, 1.56, 1.82, 2.1, 2.38, 2.66, 2.96, 3.28, 3.62, 3.96, 4.32, 4.7, 5.12, 5.54, 6, 6.48, 7, 7.5, 8.1, 8.7, 9.35, 10, 10.7, 11.5, 12.3, 13.15, 14.05, 15, 16, 17.1, 18.2, 19.5, 20.8, 22.1, 23.6, 25.1, 26.8, 28.6, 30.4, 32.4, 34.6, 36.8, 39.2, 41.8, 44.6, 47.4, 50.6, 54, 57.6, 61.4, 65.6, 70, 74.6, 79.4, 84.5, 89.5, 95, 100.5, 106.5, 112.5, 119, 125, 132, 138.5, 145.5, 153, 160.5, 168, 176, 184, 192, 200.5, 209.5, 219, 231, 248, 272, 13000] output_type = 1 


END_HISTOGRAMS

SCALES

	muf = 1.0 * Mg1g2     mur = 1.0 * Mg1g2
	muf = 0.5 * Mg1g2     mur = 0.5 * Mg1g2 
	muf = 2.0 * Mg1g2     mur = 2.0 * Mg1g2
	muf = 1.0 * Mg1g2     mur = 0.5 * Mg1g2
	muf = 1.0 * Mg1g2     mur = 2.0 * Mg1g2
	muf = 0.5 * Mg1g2     mur = 1.0 * Mg1g2
	muf = 2.0 * Mg1g2     mur = 1.0 * Mg1g2

END_SCALES


CHANNELS
  LO
END_CHANNELS