(material_howto_core)= # Howto run the material mapping and validation

ACTS Material Mapping Tutorial

(This how to uses deprecated executables and will be uptaded soon)

When performing track reconstruction, the proper amount of material crossed by the particle needs to be accounted for. This material is originally available in the detector simulation with a lot of details, which would make it expensive to directly use. To circumvent this issue, the material is mapped onto different surfaces in the tracking geometry. This process will be performed in 3 steps:

  • first, a JSON geometry file is created, it will be used to configure which surface the material is mapped onto and with which binning.

  • second, a Geant4 simulation is used to collect the material inside the detector from the detailed geometry.

  • third, all the steps are projected onto the closest surfaces (or volume in case of volume mapping) and averaged out over many events to create a map.

This page will explain how to perform the material mapping with the ACTS Examples. For this example we will use the Open Data Detector (ODD) the last paragraph will explain what needs to be changed if you want to perform the material mapping with another detector.

Prerequisites

As a prerequisite you will need to build ACTS with the Examples, Geant4 and the JSON plugin (ACTS_BUILD_EXAMPLES, ACTS_BUILD_EXAMPLES_GEANT4 and ACTS_BUILD_PLUGIN_JSON) enabled, please refer to the general how-to ACTS guide. Depending on the type of detector you want to map you will need to use some additional packages, in our case ACTS_BUILD_EXAMPLES_DD4HEP and ACTS_BUILD_PLUGIN_TGEO are needed.

For this particular example the ODD will also be needed. To use it, don’t forget to get the corresponding submodule and then recompile the ACTS code if needed.

$ git submodule init
$ git submodule update

Once Acts has been built we can start the mapping. The mapping is divided in two aspects: the surface mapping in which the material is mapped onto the closest surfaces (following the propagation direction) and the volume mapping in which the material is mapped onto a 3D (or 2D) grid associated to a volume. The first step is to select which surfaces and volumes we will want to map material onto. This is done by association of an Acts::ProtoSurfaceMaterial (or an Acts::ProtoVolumeMaterial) with the surfaces (or volumes) of interest. In the case of the ODD and some other DD4hep detectors this is done at the building step. For other detectors, or if one wants to be able to control precisely which layer will be mapped on and with which binning, an additional step is required.

Mapping and configuration

First we need to extract the list of all the surfaces and volumes in our detector, to do so we will use the GeometryExample:

$ <build>/bin/ActsExampleGeometryDD4hep -n1 -j1 \
    --mat-output-file geometry-map \
    --dd4hep-input <source>/thirdparty/OpenDataDetector/xml/OpenDataDetector.xml \
    --output-json \
    --mat-output-allmaterial true \
    --mat-output-sensitives false

This algorithm is useful to obtain a visualisation of your detector using the different types of output available (output-obj gives .obj with a 3D representation of the different subdetectors, for example). Here, we use output-json to obtain a map of all the surfaces and volumes in the detector with a ProtoSurfaceMaterial (or a ProtoVolumeMaterial), mat-output-allmaterial ensure that a ProtoSurfaceMaterial (or a ProtoVolumeMaterial) is associated to all the surfaces (or volumes), enforcing all of them to be written. Four types of surfaces exist:

  • boundaries which represent the boundaries of the different volumes

  • approaches which correspond to the entrance and exist of the detector layer

  • representing which correspond to the surface representation of a layer (often positioned at the middle of the 2 approaches)

  • sensitives which correspond to the active part of the detector (sensors)

By default, all the surfaces will be written but one can turn a specific type off (for example the sensitive) by using the appropriate option: mat-output-XXX false

The JSON file can now be edited to select which surfaces and volumes you want to have material mapped on. The JSON file is comprise of two parts, the first one contain a list of surfaces and the second a list of volumes. Information of the surface and volumes such as their type, range, id and position are available. To add one surface to the material mapping, one simply needs to switch the mapMaterial variable to true. The binning can then be changed by changing the number associated to bins, the type of bin can also be changed. For the volume, the same method can be applied, except that up to 3 bins can be associated. As a rule of thumb volume material should only be used for large homogeneous detector (like calorimeters and gaseous detectors), for the material mapping a good first try would be to use the representing surfaces of the layers with sensors. The binning depends heavily on the geometry and could be of the order of 100 (the more bins are used the more events need to be simulated to populate the bins).

Warning

When mapping onto a surface, the material inside volumes with material (or ProtoMaterial) will be ignored, you should thus avoid mapping material onto surfaces within material volumes. When mapping onto a volume, only the material within that volume will be used. If you have a large gap between the last material surface and the volume you might then want to also map material onto the boundary of the material volume.

In addition to this, the mapping type can be changed for surface mapping by changing the mappingType variable. Four different types of mapping are available:

  • PreMapping: Only map material from before the surface.

  • Default: Map material from both before and after the surface

  • PostMapping: Only map material from after the surface.

  • Sensor: Only map the last material hits before the surface. Used to map only the sensor material onto the sensors.

In case two different sufaces would receive a material hit (Default followed by Default or PreMapping for example), the material hit is associated with the closest surface.

Warning

Due to the implementation, all the material hits need to be associated with a surface. If a PostMapping surface follows a PreMapping or a Sensor surface, all the material between the two surface will be mapped onto the PostMapping surface and a warning will be issued. The same goes if the first surface encountered is a PostMapping surface.

This configuration can be cumbersome to do by hand especially when trying to map on sensitives surfaces. To simplify this task two python scripts are available in Examples/scripts/MaterialMapping:

  • writeMapConfig.py

  • configureMap.py

The first one take as an input the surfaces map previously generated and will return a JSON config file. In this file you can see all the different types of surfaces associated to each volume. You can then change the binning associated to a surface type. When the second script is called, the resulting map will have the binning applied to all the surface of the corresponding type and "mapMaterial" will be changed to true. Attention: the configureMap.py will modify the surfaces map used in input, so you might want to make a copy of it beforehand.

$ python3 <source>/Examples/Scripts/MaterialMapping/writeMapConfig.py geometry-map.json config-map.json

Then edit the config-map.json file

$ python3 <source>/Examples/Scripts/MaterialMapping/configureMap.py geometry-map.json config-map.json

Geantino scan

The next step is to do a geantino scan of our detector. For this we will use the MaterialRecording application:

$ <build>/bin/ActsExampleMaterialRecordingDD4hep -n1000 -j1 \
    --dd4hep-input <source>/thirdparty/OpenDataDetector/xml/OpenDataDetector.xml \
    --output-root

The result of the geantino scan will be a root file containing material tracks. Those contain the direction and production vertex of the geantino, the total material accumulated and all the interaction points in the detector.

Material Mapping

With the surfaces map and the material track we can finally do the material mapping using the MaterialMapping application:

$ <build>/bin/ActsExampleMaterialMappingDD4hep -j1 \
    --input-root true \
    --input-files geant4_material_tracks.root \
    --mat-input-type file \
    --mat-input-file geometry-map.json \
    --output-root \
    --output-json \
    --output-cbor \
    --mat-output-file material-maps \
    --mat-mapping-surfaces true \
    --mat-mapping-volumes true \
    --mat-mapping-volume-stepsize 1 \
    --mat-mapping-read-surfaces false \
    --dd4hep-input <source>/thirdparty/OpenDataDetector/xml/OpenDataDetector.xml

Note that technically when using DD4hep (in particular for the ODD) using the option --mat-input-type is not strictly necessary as the DD4hep geometry can hold the information of which surface to map onto with which binning. We will ignore this option, since the goal of this guide is to explain how to make a material map regardless of the detector.

As an output you will obtain the material map as a root and JSON file and a new material track collection in a root file. This new collection adds to each material interaction the associated surface during the mapping. This can be used for the control plots. Depending on what you want to do there are three options you can change:

  • mat-mapping-surfaces: determine if material is mapped onto surfaces

  • mat-mapping-volumes: determine if material is mapped onto volumes

  • mat-mapping-volume-stepsize: determine the step size used in the sampling of the volume. This should be small compared to the bin size.

In addition to root and JSON output, one can also output the material map to a Cbor file (Concise Binary Object Representation). Doing so results in a file about 10 time smaller than the JSON one, but that file is no longer human-readable. This should be done once the map has been optimised and you want to export it.

Note

You can map onto surfaces and volumes separately (for example if you want to optimise first one then the other). In that case after mapping one of those you will need to use the resulting JSON material map as an input to the mat-input-file.

Note

You might need to run the material mapping multiple times in a row to optimise the binning for each surface, which can be time-consuming. To improve the speed of subsequent mapping one can turn the mat-mapping-read-surfaces option to true, with this option the intersection of the material track with the surfaces will be taken directly from the input root file (which would have been updated during the first mapping). This will save the time needed to compute them. This can save up to 50% of the mapping time. Please note that this can only be used if the mapping surfaces used didn’t change between this execution and the last.

Material Validation

Now that the map has been written, you may want to validate it. First you can use the MaterialValidation example. This will perform propagation throughout the detector once it has been decorated with the material map. It will then output material tracks with the same format as the one obtain with the Geantino.

By default, the Geantino scan is performed with no spread in \(z_0\) and \(d_0\), while the validation has a spread of 55 mm, to obtain meaningful results, use the same spread for both (in our example a spread of 0). Another difference between the scan and the validation is that the first uses a flat distribution in \(\theta\) while the second uses a flat distribution in \(\eta\), so some reweighing might be necessary when comparing some of the distributions.

$ <build>/bin/ActsExampleMaterialValidationDD4hep -n1000 \
    --mat-input-type file \
    --mat-input-file material-maps.json \
    --output-root \
    --mat-output-file val-mat-map \
    --dd4hep-input <source>/thirdparty/OpenDataDetector/xml/OpenDataDetector.xml \
    --prop-z0-sigma 0.0 \
    --prop-d0-sigma 0.0

To do the validation, five root macros are available in scripts/MaterialMapping:

  • Mat_map.C: general comparison at the track level and 2D map of the detector.

  • Mat_map_surface_plot.C: for each mapped surface show the position of the material.

  • Mat_map_surface_plot_ratio.C: material ratio between the truth and the validation for each surface.

  • Mat_map_surface_plot_dist.C: position of the Geantino interaction with respect to the surface they are mapped on.

  • Mat_map_surface_plot_1D.C: 1D distribution of the material in each surface.

mkdir Validation

root -l -b <source>/Examples/Scripts/MaterialMapping/Mat_map.C'("propagation-material.root","material-maps_tracks.root","Validation")'
.q

mkdir Surfaces
mkdir Surfaces/prop_plot
mkdir Surfaces/map_plot
mkdir Surfaces/ratio_plot
mkdir Surfaces/dist_plot
mkdir Surfaces/1D_plot

root -l -b <source>/Examples/Scripts/MaterialMapping/Mat_map_surface_plot_ratio.C'("propagation-material.root","material-maps_tracks.root",100000,"Surfaces/ratio_plot","Surfaces/prop_plot","Surfaces/map_plot")'
.q
root -l -b <source>/Examples/Scripts/MaterialMapping/Mat_map_surface_plot_dist.C'("material-maps_tracks.root",-1,"Surfaces/dist_plot")'
.q
root -l -b <source>/Examples/Scripts/MaterialMapping/Mat_map_surface_plot_1D.C'("material-maps_tracks.root",100000,"Surfaces/1D_plot")'
.q

Using the validation plots you can then adapt the binning and the mapped surface to improve the mapping.

On top of those plots:

root -l -b <source>/Examples/Scripts/MaterialMapping/Mat_map_detector_plot_ratio.C'("propagation-material.root","material-maps_tracks.root",{X,Y,Z},100000,"Det_ratio","Det_Acts","Det_G4")'
.q

Can be use with X,Y,Z is a list of volumes, this will plot the material ratio between the map and the Geantino scan for the given volumes.

Using a different detector

If you want to use a different type of detector, you will first need to ensure that the relevant packages were added during the compilation. After that, if your detector is a DD4hep detector you will just need to replace the path given to the --dd4hep-input option. In case it is another type of detector implementation, you can replace DD4hep in the name of the algorithm by what corresponds to your detector implementation. For more information on how to include your detector in that case you can refer to the documentation of the algorithm using the -h option.