o2sim_geometry_doctor.cxx File Reference

Audits a placed geometry against the magnetic field it will be transported in. More...

#include "DetectorsBase/GeometryManager.h"
#include "Field/MagneticField.h"
#include <TFile.h>
#include <TGeoBBox.h>
#include <TGeoCone.h>
#include <TGeoEltu.h>
#include <TGeoManager.h>
#include <TGeoMatrix.h>
#include <TGeoMedium.h>
#include <TGeoNavigator.h>
#include <TGeoNode.h>
#include <TGeoPcon.h>
#include <TGeoPgon.h>
#include <TGeoShape.h>
#include <TGeoTube.h>
#include <TGeoVolume.h>
#include <TRandom3.h>
#include <TVectorD.h>
#include <boost/program_options.hpp>
#include <nlohmann/json.hpp>
#include <algorithm>
#include <atomic>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <ctime>
#include <fstream>
#include <iostream>
#include <map>
#include <set>
#include <string>
#include <thread>
#include <utility>
#include <vector>

Go to the source code of this file.

Typedefs
using	json = nlohmann::json

Functions
int	main (int argc, char **argv)

Detailed Description

Audits a placed geometry against the magnetic field it will be transported in.

The tool reads a geometry file and a magnetic field, and answers two questions that cannot be answered by looking at either one alone:

• where does the geometry fail to express the field's structure – a logical volume placed both inside and outside the field, or a mother whose own material straddles the boundary, so that no per-medium field flag can be right for all of its placements;
• where does it already express it wrongly – a medium with ifield == 0, which asks the transport engine for straight lines, sitting in real field.

It also answers one question about the geometry alone, and therefore runs that part without a field under –reachability-only: does every placement actually occupy the space it was built in? A daughter outside its mother, or one shadowed by an overlapping sibling, is never reached by the navigator, carries no material and produces no hits, and nothing in the construction code says so.

It detects and proposes. It never modifies a geometry.

The field enters in two ways. A support model (an outer bound on where |B| exceeds a threshold, maximised over phi) supplies the geometric argument: a placement whose extent misses every band of the model is field-free, and that is the only way this tool ever concludes "field-free". Sampling enters only in the opposite direction, to disprove – finding field inside a volume settles the question, finding none does not.

Reading the geometry needs no detector code: a medium's ifield flag and its sensitivity are GSTMED parameters 1 and 0 of the TGeoMedium, and are recovered from the file itself.

Definition in file o2sim_geometry_doctor.cxx.

Typedef Documentation

json

using json = nlohmann::json

Definition at line 83 of file o2sim_geometry_doctor.cxx.

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

Definition at line 1973 of file o2sim_geometry_doctor.cxx.