#################################
####### GENERAL PARAMETERS ######
#################################
max_step = 2
amr.n_cell =  32 32
amr.max_grid_size = 16   # maximum size of each AMReX box, used to decompose the domain
amr.blocking_factor = 8 # minimum size of each AMReX box, used to decompose the domain
geometry.dims = 2
geometry.prob_lo     =  -0.5e-6   -0.5e-6  # physical domain
geometry.prob_hi     =  0.5e-6   0.5e-6
amr.max_level = 0 # Maximum level in hierarchy (1 might be unstable, >1 is not supported)

#################################
####### Boundary condition ######
#################################
boundary.field_lo = periodic periodic
boundary.field_hi = periodic periodic

#################################
############ NUMERICS ###########
#################################
algo.current_deposition = esirkepov
algo.charge_deposition = standard
algo.field_gathering = energy-conserving
algo.particle_pusher = boris
warpx.verbose = 1
warpx.do_dive_cleaning = 0
warpx.use_filter = 1
warpx.cfl = 1. # if 1., the time step is set to its CFL limit
warpx.serialize_initial_conditions = 1

# Order of particle shape factors
algo.particle_shape = 3

#################################
############ PLASMA #############
#################################
particles.species_names = p1 p2 p3 p4 ele1 ele2 ele3 ele4 pos1 pos2 pos3 pos4 dummy_phot
particles.photon_species = p1 p2 p3 p4 dummy_phot
#################################

p1.species_type = "photon"
p1.injection_style = "NUniformPerCell"
p1.profile = "constant"
p1.num_particles_per_cell_each_dim = 32 32
p1.density = 1e19
p1.momentum_distribution_type = "constant"
p1.ux = 2000.0
##########QED####################
p1.do_qed_breit_wheeler = 1
p1.qed_breit_wheeler_ele_product_species = ele1
p1.qed_breit_wheeler_pos_product_species = pos1
#################################

p2.species_type = "photon"
p2.injection_style = "NUniformPerCell"
p2.profile = "constant"
p2.num_particles_per_cell_each_dim = 32 32
p2.density = 1e19
p2.momentum_distribution_type = "constant"
p2.uy = 5000.0
##########QED####################
p2.do_qed_breit_wheeler = 1
p2.qed_breit_wheeler_ele_product_species = ele2
p2.qed_breit_wheeler_pos_product_species = pos2
#################################

p3.species_type = "photon"
p3.injection_style = "NUniformPerCell"
p3.profile = "constant"
p3.num_particles_per_cell_each_dim = 32 32
p3.density = 1e19
p3.momentum_distribution_type = "constant"
p3.uz = 10000.0
##########QED####################
p3.do_qed_breit_wheeler = 1
p3.qed_breit_wheeler_ele_product_species = ele3
p3.qed_breit_wheeler_pos_product_species = pos3
#################################

p4.species_type = "photon"
p4.injection_style = "NUniformPerCell"
p4.profile = "constant"
p4.num_particles_per_cell_each_dim = 32 32
p4.density = 1e19
p4.momentum_distribution_type = "constant"
p4.ux = 57735.02691896
p4.uy = 57735.02691896
p4.uz = 57735.02691896
##########QED####################
p4.do_qed_breit_wheeler = 1
p4.qed_breit_wheeler_ele_product_species = ele4
p4.qed_breit_wheeler_pos_product_species = pos4
#################################

### PRODUCT SPECIES ###
ele1.species_type = "electron"
ele1.injection_style = "none"
ele1.do_not_push = 1
ele1.do_qed_quantum_sync = 1
ele1.qed_quantum_sync_phot_product_species = dummy_phot

ele2.species_type = "electron"
ele2.injection_style = "none"
ele2.do_not_push = 1
ele2.do_qed_quantum_sync = 1
ele2.qed_quantum_sync_phot_product_species = dummy_phot

ele3.species_type = "electron"
ele3.injection_style = "none"
ele3.do_not_push = 1
ele3.do_qed_quantum_sync = 1
ele3.qed_quantum_sync_phot_product_species = dummy_phot

ele4.species_type = "electron"
ele4.injection_style = "none"
ele4.do_not_push = 1
ele4.do_qed_quantum_sync = 1
ele4.qed_quantum_sync_phot_product_species = dummy_phot

pos1.species_type = "positron"
pos1.injection_style = "none"
pos1.do_not_push = 1
pos1.do_qed_quantum_sync = 1
pos1.qed_quantum_sync_phot_product_species = dummy_phot

pos2.species_type = "positron"
pos2.injection_style = "none"
pos2.do_not_push = 1
pos2.do_qed_quantum_sync = 1
pos2.qed_quantum_sync_phot_product_species = dummy_phot

pos3.species_type = "positron"
pos3.injection_style = "none"
pos3.do_not_push = 1
pos3.do_qed_quantum_sync = 1
pos3.qed_quantum_sync_phot_product_species = dummy_phot

pos4.species_type = "positron"
pos4.injection_style = "none"
pos4.do_not_push = 1
pos4.do_qed_quantum_sync = 1
pos4.qed_quantum_sync_phot_product_species = dummy_phot

dummy_phot.species_type = "photon"
dummy_phot.injection_style = "none"

#################################

##########QED TABLES####################
qed_bw.chi_min = 0.001

qed_bw.lookup_table_mode = "builtin"

#qed_bw.lookup_table_mode = "generate"
#qed_bw.tab_dndt_chi_min = 0.01
#qed_bw.tab_dndt_chi_max = 1000.0
#qed_bw.tab_dndt_how_many = 256
#qed_bw.tab_pair_chi_min = 0.01
#qed_bw.tab_pair_chi_max = 1000.0
#qed_bw.tab_pair_chi_how_many = 256
#qed_bw.tab_pair_frac_how_many = 256
#qed_bw.save_table_in = "bw_table"

#qed_bw.lookup_table_mode = "load"
#qed_bw.load_table_from = "bw_table"

qed_qs.chi_min = 0.001

qed_qs.lookup_table_mode = "builtin"

qed_qs.photon_creation_energy_threshold = 2

#qed_qs.lookup_table_mode = "generate"
#qed_qs.tab_dndt_chi_min = 0.001
#qed_qs.tab_dndt_chi_max = 1000.0
#qed_qs.tab_dndt_how_many = 256
#qed_qs.tab_em_chi_min = 0.001
#qed_qs.tab_em_frac_min = 1.0e-12
#qed_qs.tab_em_chi_max = 1000.0
#qed_qs.tab_em_chi_how_many = 256
#qed_qs.tab_em_frac_how_many = 256
#qed_qs.save_table_in = "qs_table"

#qed_qs.lookup_table_mode = "load"
#qed_qs.load_table_from = "qs_table"
#################################

### EXTERNAL E FIELD ### (3e15 * [-0.811107105653813 0.324442842261525 0.486664263392288] )
particles.E_ext_particle_init_style = "constant"
particles.E_external_particle = -2433321316961438 973328526784575 1459992790176863
####

### EXTERNAL B FIELD ### (1e7 * [0.28571429 0.42857143 0.85714286] )
particles.B_ext_particle_init_style = "constant"
particles.B_external_particle = 2857142.85714286 4285714.28571428 8571428.57142857
####

# Diagnostics
diagnostics.diags_names = diag1
diag1.intervals = 2
diag1.diag_type = Full
diag1.fields_to_plot = Ex
diag1.p1.variables = ux uy uz w opticalDepthBW
diag1.p2.variables = ux uy uz w opticalDepthBW
diag1.p3.variables = ux uy uz w opticalDepthBW
diag1.p4.variables = ux uy uz w opticalDepthBW

diag1.ele1.variables = ux uy uz w opticalDepthQSR
diag1.ele2.variables = ux uy uz w opticalDepthQSR
diag1.ele3.variables = ux uy uz w opticalDepthQSR
diag1.ele4.variables = ux uy uz w opticalDepthQSR

diag1.pos1.variables = ux uy uz w opticalDepthQSR
diag1.pos2.variables = ux uy uz w opticalDepthQSR
diag1.pos3.variables = ux uy uz w opticalDepthQSR
diag1.pos4.variables = ux uy uz w opticalDepthQSR

diag1.format = plotfile