Hi everybody,
I’ve run into a weird issue during my most recent allpix simulations.
I’m trying to simulate the behaviour of a fixed set of charges deposited into a silicon strip sensor, under a TCAD simulated electric field and Ramo potential and a constant magnetic field. My configuration file for reference:
[Allpix]
log_level = “WARNING”
log_format = “DEFAULT”
model_paths = “/eos/user/d/damir/ALLPIX2_Remote/Prototype/Configurations”
detectors_file = “Configurations/detectors_file.conf”
number_of_events = 100
root_file = “histos_prototype_5strip_3DW_transprop_MIP”
random_seed = 0#The geometry is empty for now
[GeometryBuilderGeant4]#Magnetic Field Reader
[MagneticFieldReader]
model = “constant”
magnetic_field = 0T 2T 0T#Depositing the charges directly
[DepositionPointCharge]
source_type = “point”
model = “fixed”
position = 0um 0um 50um
number_of_charges = 1000
output_plots = 1#Electric Field Reader
[ElectricFieldReader]
model = “mesh”
file_name = “TCAD_Fields/fields_1strip_des_ElectricField_COARSE.init”
output_plots = 1#Ramo Potential Reader
[WeightingPotentialReader]
model = “mesh”
file_name = “TCAD_Fields/fields_1strip_des_Weighting_WeightingPotential_COARSE.init”
output_plots = 1#Using Transient Propagation since it uses the Ramo-Schokley Theorem
[TransientPropagation]
temperature = 253K
charge_per_step = 10
mobility_model = “canali”
trapping_model = “cmstracker”
fluence = 4e15/cm/cm
induction_matrix = 3 3
output_plots = 1#Pulse Transfer
[PulseTransfer]
max_depth_distance = 5um
output_plots = 1#Default Digitizer
[DefaultDigitizer]
output_plots = 1#Detector Histogrammer
[DetectorHistogrammer]#[VisualizationGeant4]
#mode = “gui”#Define the .root file that contains the tree, and what leaves the tree contains
[ROOTObjectWriter]
file_name = “trees_prototype_5strip_3DW_transprop_MIP”
include = “MCTrack”,“MCParticle”,“PixelCharge”,“PixelHit”,“PropagatedCharge”,“DepositedCharge”
As you can see at the end of the file, I explicitly ask allpix to save “PixelHit” to the output .root file. This tree is meant to contain the hit information I need to derive the charge collection efficiency.
However, when I look in the .root file, there is no such tree to be found. Even more bizarrely, one of the trees that come after it; “PropagatedCharge”, is present in the .root file, despite being dependent on “PixelHit”. I’ve attached a screenshot of the TBrowser I used to check:
Why is this happening? I’m tempted to say that this is some sort of bug, since I remember seeing someone posting on the forum about having a similar issue a couple of months ago.