[Allpix]
log_level = "INFO"
log_format = "DEFAULT"
detectors_file = "planar_detector.conf" 
#detectors_file = "../detectors/normal_incidence_planar_detector.conf" 
number_of_events = 1000
#skip_events = 1000
root_file = "histos_250um"



[GeometryBuilderGeant4]
world_material = "air"

[DepositionGeant4]
physics_list = FTFP_BERT_LIV
particle_type = "Pi+"
number_of_particles = 1
source_energy = 120GeV
source_position = 0um 0um -200mm
source_type = "beam"
beam_size = 2mm
beam_divergence = 1mrad 1mrad
beam_direction = 0 0 1
max_step_length = 10.0um
output_plots=1



[ElectricFieldReader]
model = "custom"
field_function = "[0]*z*z*z*z*z+[1]*z*z*z*z+[2]*z*z*z+[3]*z*z+[4]*z+[5]"
#field_parameters = 878148V/mm/mm/mm/mm/mm, 73361.7V/mm/mm/mm/mm, -10198.9V/mm/mm/mm, 217.845V/mm/mm, -781.697V/mm, 603.915V
field_parameters = -865694.84V/mm/mm/mm/mm/mm/mm, -72321.346V/mm/mm/mm/mm/mm, 10054.227V/mm/mm/mm/mm, -214.75569V/mm/mm/mm, 770.61157V/mm/mm, -595.35017V/mm
output_plots=1


[WeightingPotentialReader]
model = "mesh"
#file_name = "/home/mbomben/work/Silvaco/Victory/230630_fei3_ramo/mirrored_shifted_completed_Ramo_Potential_RamoPotential.init"
file_name = "test_Blayer.init"
#file_name = "/home/mbomben/work/Digitizer/ID/Disks/Ramo/fei3_WP_weightingpotential.init"
output_plots = 1
field_mapping = PIXEL_FULL

#For TCAD Efield
[TransientPropagation] # only for TCAD efield
temperature = 266K 
charge_per_step = 10 
mobility_model = "canali"
trapping_model = "cmstracker"
fluence = 3e13/cm/cm
induction_matrix = 3 3  
output_plots = 1

[PulseTransfer]
max_depth_distance = 5um
output_plots = 1

[DefaultDigitizer]
output_plots = 1
threshold = 3500e

[DetectorHistogrammer]

#[Ignore]
[ROOTObjectWriter]
#file_name = "1kevents_ROOTObjectWriterData_TCADEfield_PulseTransf_4e15fluence_cmstrackerTrappingmodel.root"
file_name = "trees_250um.root"
include = "MCTrack","MCParticle","PixelCharge","PixelHit","PropagatedCharge","DepositedCharge"