diff --git a/pandoraPFA/CLD_LAr/BeamInstrumentation_o3_v01_overlap.xml b/pandoraPFA/CLD_LAr/BeamInstrumentation_o3_v01_overlap.xml
new file mode 100644
index 0000000..634585f
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/BeamInstrumentation_o3_v01_overlap.xml
@@ -0,0 +1,26 @@
+
+
+
+ Beampipe Instrumentation
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/BeamInstrumentation_o3_v02_fitShield.xml b/pandoraPFA/CLD_LAr/BeamInstrumentation_o3_v02_fitShield.xml
new file mode 100644
index 0000000..dba09d3
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/BeamInstrumentation_o3_v02_fitShield.xml
@@ -0,0 +1,26 @@
+
+
+
+ Beampipe Instrumentation
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/Beampipe_o4_v04_noNotch_W_n02.xml b/pandoraPFA/CLD_LAr/Beampipe_o4_v04_noNotch_W_n02.xml
new file mode 100644
index 0000000..781ed74
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/Beampipe_o4_v04_noNotch_W_n02.xml
@@ -0,0 +1,140 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Part of beampipe made of Beryllium
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Golden foil in the inner part of the Be beampipe
+
+
+
+
+
+ Part of beampipe made of Copper
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+Full Cone Tungsten Shield
+
+
+
+ Before HOM space
+
+
+ After HOM space (1197.5*m - 1298.7*mm) +18 cm as solenoid is now closer to IP
+
+
+
+Asymmetric Tungsten Shield no Rotation
+
+
+
+
+
+
+ was 370
+
+
+ one degree less, to fit lumical window
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/ECalBarrel_o2_v01_03.xml b/pandoraPFA/CLD_LAr/ECalBarrel_o2_v01_03.xml
new file mode 100644
index 0000000..abe72e0
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/ECalBarrel_o2_v01_03.xml
@@ -0,0 +1,66 @@
+
+
+
+
+
+
+
+
+
+
+ system:5,side:2,module:8,stave:4,layer:9,submodule:4,x:32:-16,y:-16
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ EM Calorimeter Barrel
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/ECalEndcap_o2_v01_03.xml b/pandoraPFA/CLD_LAr/ECalEndcap_o2_v01_03.xml
new file mode 100644
index 0000000..dd84c50
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/ECalEndcap_o2_v01_03.xml
@@ -0,0 +1,57 @@
+
+
+
+
+
+
+
+
+ system:5,side:2,module:8,stave:4,layer:9,submodule:4,x:32:-16,y:-16
+
+
+
+
+
+
+
+
+ Electromagnetic Calorimeter Endcap
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/FCCee_o1_v05.xml b/pandoraPFA/CLD_LAr/FCCee_o1_v05.xml
new file mode 100644
index 0000000..9354561
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/FCCee_o1_v05.xml
@@ -0,0 +1,433 @@
+
+
+
+ The compact format for the FCCee Detector design
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ suggested naming convention:
+
+ main parameters:
+
+ DET_inner_radius : inner radius of tube like envelope ( inscribed cylinder )
+ DET_outer_radius : outer radius of tube like envelope ( circumscribed cylinder )
+ DET_half_length : half length along z axis
+ DET_min_z : smallest absolute value on z-axis
+ DET_max_z : largest absolute value on z-axis
+ DET_inner_symmetry : number of sides on the inside ( 0 for tube )
+ DET_outer_symmetry : number of sides on the inside ( 0 for tube )
+ DET_inner_phi0 : optional rotation of the inner polygon ( in r-phi plane )
+ DET_outer_phi0 : optional rotation of the outer polygon ( in r-phi plane )
+
+ additional parameters for cutting away volumes/shapes use one of the above with a number
+ appended and/or an extra specifiaction such as cone ( for a cut away cone )
+
+ DET_inner_radius_1
+ DET_outer_radius_2
+ DET_cone_min_z
+ DET_cone_max_z
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/HCalBarrel_o1_v01_01.xml b/pandoraPFA/CLD_LAr/HCalBarrel_o1_v01_01.xml
new file mode 100644
index 0000000..d0ade5f
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/HCalBarrel_o1_v01_01.xml
@@ -0,0 +1,58 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ system:5,side:2,module:8,stave:4,layer:9,submodule:4,x:32:-16,y:-16
+
+
+
+
+
+
+ Hadron Calorimeter Barrel
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/HCalEndcap_o1_v01_01.xml b/pandoraPFA/CLD_LAr/HCalEndcap_o1_v01_01.xml
new file mode 100644
index 0000000..ef70e41
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/HCalEndcap_o1_v01_01.xml
@@ -0,0 +1,131 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ HCalEndcap Assembly
+
+
+
+
+
+
+
+
+
+
+
+
+ system:5,side:2,module:8,stave:4,layer:9,submodule:4,x:32:-16,y:-16
+
+
+
+ system:5,side:2,module:8,stave:4,layer:9,submodule:4,x:32:-16,y:-16
+
+
+
+
+
+
+
+
+ Hadronic Calorimeter Endcap
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Hadronic Calorimeter Endcap
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/HOMAbsorber_v00.xml b/pandoraPFA/CLD_LAr/HOMAbsorber_v00.xml
new file mode 100644
index 0000000..f197ea2
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/HOMAbsorber_v00.xml
@@ -0,0 +1,53 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/InnerTrackerBarrelModuleDown.xml b/pandoraPFA/CLD_LAr/InnerTrackerBarrelModuleDown.xml
new file mode 100644
index 0000000..12d5d74
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/InnerTrackerBarrelModuleDown.xml
@@ -0,0 +1,25 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff --git a/pandoraPFA/CLD_LAr/InnerTrackerBarrelModuleUp.xml b/pandoraPFA/CLD_LAr/InnerTrackerBarrelModuleUp.xml
new file mode 100644
index 0000000..a814182
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/InnerTrackerBarrelModuleUp.xml
@@ -0,0 +1,25 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff --git a/pandoraPFA/CLD_LAr/InnerTracker_o2_v06_02.xml b/pandoraPFA/CLD_LAr/InnerTracker_o2_v06_02.xml
new file mode 100644
index 0000000..fd108c3
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/InnerTracker_o2_v06_02.xml
@@ -0,0 +1,596 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Tracking detectors
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Inner Tracker Assembly
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ ${GlobalTrackerReadoutID}
+
+
+ ${GlobalTrackerReadoutID}
+
+
+
+
+
+
+
+
+
+ Silicon Inner Tracker Barrel
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Silicon Inner Tracker Endcaps
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ 3.5959*cm = X0 for Carbon fibre
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ The next section is the cable for the vertex
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/LAr_ECalBarrel.xml b/pandoraPFA/CLD_LAr/LAr_ECalBarrel.xml
new file mode 100644
index 0000000..85ce8fe
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/LAr_ECalBarrel.xml
@@ -0,0 +1,144 @@
+
+
+
+
+
+ Settings for the inclined EM calorimeter.
+ The barrel is filled with liquid argon. Passive material includes lead in the middle and steal on the outside, glued together.
+ Passive plates are inclined by a certain angle from the radial direction.
+ In between of two passive plates there is a readout.
+ Space between the plate and readout is of trapezoidal shape and filled with liquid argon.
+ Definition of sizes, visualization settings, readout and longitudinal segmentation are specified.
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ system:4,cryo:1,type:3,subtype:3,layer:8,module:11,eta:9
+
+
+
+
+
+
+ system:4,cryo:1,type:3,subtype:3,layer:8,eta:9,phi:10
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/LAr_ECalEndcaps.xml b/pandoraPFA/CLD_LAr/LAr_ECalEndcaps.xml
new file mode 100644
index 0000000..e76463b
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/LAr_ECalEndcaps.xml
@@ -0,0 +1,167 @@
+
+
+
+
+
+ Liquid argon EM calorimeter endcap design.
+ Electromagnetic part (EMEC) includes lead+steel absorber.
+ Absorbers have shape of simple discs.
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ system:4,subsystem:1,type:3,subtype:3,layer:8,sublayer:8,eta:10,phi:10
+
+
+
+ system:4,subsystem:1,type:3,subtype:3,layer:8,eta:10,phi:10
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/LumiCal_o3_v02_03.xml b/pandoraPFA/CLD_LAr/LumiCal_o3_v02_03.xml
new file mode 100644
index 0000000..2ba7b40
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/LumiCal_o3_v02_03.xml
@@ -0,0 +1,180 @@
+
+
+
+
+
+
+
+
+
+ system:8,barrel:3,layer:8,slice:8,r:32:-16,phi:-16
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/OuterTrackerBarrelModuleDown.xml b/pandoraPFA/CLD_LAr/OuterTrackerBarrelModuleDown.xml
new file mode 100644
index 0000000..8c03038
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/OuterTrackerBarrelModuleDown.xml
@@ -0,0 +1,27 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff --git a/pandoraPFA/CLD_LAr/OuterTrackerBarrelModuleUp.xml b/pandoraPFA/CLD_LAr/OuterTrackerBarrelModuleUp.xml
new file mode 100644
index 0000000..12a1efe
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/OuterTrackerBarrelModuleUp.xml
@@ -0,0 +1,27 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff --git a/pandoraPFA/CLD_LAr/OuterTracker_o2_v06_02.xml b/pandoraPFA/CLD_LAr/OuterTracker_o2_v06_02.xml
new file mode 100644
index 0000000..272511d
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/OuterTracker_o2_v06_02.xml
@@ -0,0 +1,226 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Tracking detectors
+
+
+
+
+ Outer Tracker Assembly
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ ${GlobalTrackerReadoutID}
+
+
+ ${GlobalTrackerReadoutID}
+
+
+
+
+
+
+
+
+
+
+ Silicon Outer Tracker Barrel
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Silicon Outer Tracker Endcaps
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/Solenoid_o1_v01_02.xml b/pandoraPFA/CLD_LAr/Solenoid_o1_v01_02.xml
new file mode 100644
index 0000000..634f3ca
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/Solenoid_o1_v01_02.xml
@@ -0,0 +1,54 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Solenoid
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/TrackerDiskModuleIn.xml b/pandoraPFA/CLD_LAr/TrackerDiskModuleIn.xml
new file mode 100644
index 0000000..12a1efe
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/TrackerDiskModuleIn.xml
@@ -0,0 +1,27 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff --git a/pandoraPFA/CLD_LAr/TrackerDiskModuleOut.xml b/pandoraPFA/CLD_LAr/TrackerDiskModuleOut.xml
new file mode 100644
index 0000000..8c03038
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/TrackerDiskModuleOut.xml
@@ -0,0 +1,27 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff --git a/pandoraPFA/CLD_LAr/Vertex_o4_v05.xml b/pandoraPFA/CLD_LAr/Vertex_o4_v05.xml
new file mode 100644
index 0000000..7a67533
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/Vertex_o4_v05.xml
@@ -0,0 +1,223 @@
+
+
+
+
+ Tracking detectors
+
+
+
+
+
+
+
+
+
+
+
+ Vertex Assembly
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ ${GlobalTrackerReadoutID}
+
+
+ ${GlobalTrackerReadoutID}
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Vertex Detector Endcaps
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/YokeBarrel_o1_v01_02.xml b/pandoraPFA/CLD_LAr/YokeBarrel_o1_v01_02.xml
new file mode 100644
index 0000000..5d4e62d
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/YokeBarrel_o1_v01_02.xml
@@ -0,0 +1,68 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ system:5,side:2,layer:9,module:8,stave:4,submodule:4,x:32:-16,y:-16
+
+
+
+
+
+
+
+ Yoke Calorimeter Barrel
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/YokeEndcap_o1_v01_02.xml b/pandoraPFA/CLD_LAr/YokeEndcap_o1_v01_02.xml
new file mode 100644
index 0000000..e450384
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/YokeEndcap_o1_v01_02.xml
@@ -0,0 +1,67 @@
+
+
+
+
+
+
+ system:5,side:2,module:8,stave:4,layer:9,submodule:4,x:32:-16,y:-16
+
+
+
+
+
+
+
+
+
+
+
+
+
+ Encap Yoke
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/CLD_LAr/elements.xml b/pandoraPFA/CLD_LAr/elements.xml
new file mode 100644
index 0000000..8358bac
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/elements.xml
@@ -0,0 +1,884 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff --git a/pandoraPFA/CLD_LAr/materials.xml b/pandoraPFA/CLD_LAr/materials.xml
new file mode 100644
index 0000000..d1346b6
--- /dev/null
+++ b/pandoraPFA/CLD_LAr/materials.xml
@@ -0,0 +1,256 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/pandoraPFA/README.md b/pandoraPFA/README.md
new file mode 100644
index 0000000..80ad358
--- /dev/null
+++ b/pandoraPFA/README.md
@@ -0,0 +1,48 @@
+# Tools for running PandoraPFA on the noble liquid detector concept
+
+Particle Flow Algorithms (PFA) need the full detector information (or at least tracks and calorimeter hits).
+Since we currently do not have track in the drift chamber, one workaround to play with PFA is to take the CLD detector and replace the calorimeter by the Noble Liquid one.
+
+The `CLD_LAr` folder hosts a detector configuration corresponding to CLD from `k4geo` v00-18 where the SiW ECAL Barrel was replaced by the Noble Liquid Barrel one from `FCCDetectors` v0.1pre09 (it is not yet good for physics, just use it to start playing with PFA, once a full Noble Liquid based detector full sim is available we should move to this)
+
+## Recipe
+
+1) Make sure you have installed a local version of FCCDetectors and that you have set the environment variables properly:
+
+To be done only once:
+```
+git clone https://github.com/HEP-FCC/FCCDetectors
+cd FCCDetectors
+mkdir build install
+cd build
+cmake .. -DCMAKE_INSTALL_PREFIX=../install
+make install -j 8
+cd ../../
+```
+
+To be done each time (in the folder containing FCCDetectors)
+```
+export FCCDETECTORS=$PWD/FCCDetectors/
+PATH=$PWD/FCCDetectors/install/bin/:$PATH
+CMAKE_PREFIX_PATH=$PWD/FCCDetectors/install/:$CMAKE_PREFIX_PATH
+LD_LIBRARY_PATH=$PWD/FCCDetectors/install/lib:$LD_LIBRARY_PATH
+export PYTHONPATH=$PWD/FCCDetectors/install/python:$PYTHONPATH
+LD_LIBRARY_PATH=$PWD/FCCDetectors/install/lib64:$LD_LIBRARY_PATH
+```
+
+2) Produce events with Simhits from CLD_LAr:
+`./run_ddsim_cld_lar.sh`
+
+3) Run the reconstruction on it to get tracks and LAr ECAL deposits:
+`fccrun track_gaudi_produce_CLD_events_from_ddsim.py`
+
+4) TO BE IMPLEMENTED: run MarlinPandoraPFA on the output root file of the previous step
+
+## Ongoing attempt
+
+Currently trying to produce CLD_LAr events fully through Gaudi with
+
+`produce_CLD_LAr_events.py`
+
+but it is not working at the moment due to difference in the way Geant4 steps are saved in ddsim and in the FCC framework.
+
diff --git a/pandoraPFA/run_ddsim_cld_lar.sh b/pandoraPFA/run_ddsim_cld_lar.sh
new file mode 100755
index 0000000..7a42aab
--- /dev/null
+++ b/pandoraPFA/run_ddsim_cld_lar.sh
@@ -0,0 +1,7 @@
+ddsim --compactFile $PWD/CLD_LAr/FCCee_o1_v05.xml \
+ --enableGun \
+ --gun.distribution uniform \
+ --gun.energy "10*GeV" \
+ --gun.particle e- \
+ --numberOfEvents 10 \
+ --outputFile Step1_edm4hep.root
diff --git a/pandoraPFA/track_gaudi_produce_CLD_events_from_ddsim.py b/pandoraPFA/track_gaudi_produce_CLD_events_from_ddsim.py
new file mode 100644
index 0000000..e41a264
--- /dev/null
+++ b/pandoraPFA/track_gaudi_produce_CLD_events_from_ddsim.py
@@ -0,0 +1,767 @@
+import os
+
+from GaudiKernel.SystemOfUnits import MeV, GeV, tesla
+
+use_pythia = False
+addNoise = False
+
+# Input for simulations (momentum is expected in GeV!)
+# Parameters for the particle gun simulations, dummy if use_pythia is set to True
+# theta from 80 to 100 degrees corresponds to -0.17 < eta < 0.17
+momentum = 10 # in GeV
+#thetaMin = 90.25 # degrees
+#thetaMax = 90.25 # degrees
+thetaMin = 20 # degrees
+thetaMax = 130 # degrees
+pdgCode = 11 # 11 electron, 13 muon, 22 photon, 111 pi0, 211 pi+
+magneticField = False
+
+from Gaudi.Configuration import *
+
+from Configurables import FCCDataSvc
+podioevent = FCCDataSvc("EventDataSvc")
+
+################## Particle gun setup
+_pi = 3.14159
+
+from Configurables import GenAlg
+genAlg = GenAlg()
+if use_pythia:
+ from Configurables import PythiaInterface
+ pythia8gentool = PythiaInterface()
+ pythia8gentool.pythiacard = os.path.join(os.environ.get('PWD', ''), "MCGeneration/ee_Zgamma_inclusive.cmd")
+ #pythia8gentool.pythiacard = "MCGeneration/ee_Z_ee.cmd"
+ pythia8gentool.printPythiaStatistics = False
+ pythia8gentool.pythiaExtraSettings = [""]
+ genAlg.SignalProvider = pythia8gentool
+ # to smear the primary vertex position:
+ #from Configurables import GaussSmearVertex
+ #smeartool = GaussSmearVertex()
+ #smeartool.xVertexSigma = 0.5*units.mm
+ #smeartool.yVertexSigma = 0.5*units.mm
+ #smeartool.zVertexSigma = 40.0*units.mm
+ #smeartool.tVertexSigma = 180.0*units.picosecond
+ #genAlg.VertexSmearingTool = smeartool
+else:
+ from Configurables import MomentumRangeParticleGun
+ pgun = MomentumRangeParticleGun("ParticleGun_Electron")
+ pgun.PdgCodes = [pdgCode]
+ pgun.MomentumMin = momentum * GeV
+ pgun.MomentumMax = momentum * GeV
+ pgun.PhiMin = 0
+ #pgun.PhiMax = 0
+ pgun.PhiMax = 2 * _pi
+ pgun.ThetaMin = thetaMin * _pi / 180.
+ pgun.ThetaMax = thetaMax * _pi / 180.
+ genAlg.SignalProvider = pgun
+
+genAlg.hepmc.Path = "hepmc"
+
+from Configurables import HepMCToEDMConverter
+hepmc_converter = HepMCToEDMConverter()
+hepmc_converter.hepmc.Path="hepmc"
+genParticlesOutputName = "MCParticles"
+hepmc_converter.GenParticles.Path = genParticlesOutputName
+hepmc_converter.hepmcStatusList = []
+
+################## Simulation setup
+# Detector geometry
+from Configurables import GeoSvc
+geoservice = GeoSvc("GeoSvc")
+# if FCC_DETECTORS is empty, this should use relative path to working directory
+#path_to_detector = os.environ.get("FCCDETECTORS", "")
+path_to_detector = "./"
+print(path_to_detector)
+detectors_to_use=[
+ #'CLD_LAr/FCCee_o1_v05.xml',
+ #'../../FCCDetectors/Detector/DetFCCeeIDEA-LAr/compact/FCCee_DectMaster.xml',
+ '../../k4geo/FCCee/compact/FCCee_o1_v05/FCCee_o1_v05.xml',
+ ]
+# prefix all xmls with path_to_detector
+geoservice.detectors = [os.path.join(path_to_detector, _det) for _det in detectors_to_use]
+geoservice.OutputLevel = INFO
+
+# Geant4 service
+# Configures the Geant simulation: geometry, physics list and user actions
+
+from Configurables import SimG4FullSimActions, SimG4Alg, SimG4PrimariesFromEdmTool, SimG4SaveParticleHistory
+actions = SimG4FullSimActions()
+# Uncomment if history from Geant4 decays is needed (e.g. to get the photons from pi0) and set actions=actions in SimG4Svc + uncomment saveHistTool in SimG4Alg
+#actions.enableHistory=True
+#actions.energyCut = 0.2 * GeV
+#saveHistTool = SimG4SaveParticleHistory("saveHistory")
+
+from Configurables import SimG4Svc
+geantservice = SimG4Svc("SimG4Svc", detector='SimG4DD4hepDetector', physicslist="SimG4FtfpBert", actions=actions)
+
+# Fixed seed to have reproducible results, change it for each job if you split one production into several jobs
+# Mind that if you leave Gaudi handle random seed and some job start within the same second (very likely) you will have duplicates
+geantservice.randomNumbersFromGaudi = False
+geantservice.seedValue = 4242
+
+# Range cut
+geantservice.g4PreInitCommands += ["/run/setCut 0.1 mm"]
+
+# Magnetic field
+from Configurables import SimG4ConstantMagneticFieldTool
+if magneticField == 1:
+ field = SimG4ConstantMagneticFieldTool("SimG4ConstantMagneticFieldTool", FieldComponentZ=-2*tesla, FieldOn=True,IntegratorStepper="ClassicalRK4")
+else:
+ field = SimG4ConstantMagneticFieldTool("SimG4ConstantMagneticFieldTool",FieldOn=False)
+
+
+# input file from ddsim
+from Configurables import k4DataSvc, PodioInput
+evtsvc = k4DataSvc('EventDataSvc')
+evtsvc.input = 'Step1_edm4hep.root'
+inp = PodioInput('InputReader')
+inp.collections = [
+ 'EventHeader',
+ 'MCParticles',
+ 'VertexBarrelCollection',
+ 'VertexEndcapCollection',
+ 'InnerTrackerBarrelCollection',
+ 'OuterTrackerBarrelCollection',
+ 'ECalEndcapCollection',
+ 'ECalEndcapCollectionContributions',
+ 'ECalBarrelEta',
+ 'ECalBarrelEtaContributions',
+ #'ECalBarrelCollection',
+ #'ECalBarrelCollectionContributions',
+ 'HCalBarrelCollection',
+ 'HCalBarrelCollectionContributions',
+ 'InnerTrackerEndcapCollection',
+ 'OuterTrackerEndcapCollection',
+ 'HCalEndcapCollection',
+ 'HCalEndcapCollectionContributions',
+ 'HCalRingCollection',
+ 'HCalRingCollectionContributions',
+ 'YokeBarrelCollection',
+ 'YokeBarrelCollectionContributions',
+ 'YokeEndcapCollection',
+ 'YokeEndcapCollectionContributions',
+ 'LumiCalCollection',
+ 'LumiCalCollectionContributions',
+]
+#inp.OutputLevel = WARNING
+
+# Detector readouts
+from Configurables import SimG4SaveTrackerHits
+saveVTXBsimHitTool = SimG4SaveTrackerHits("saveVTXBsimHitTool", readoutNames=["VertexBarrelCollection"])
+saveVTXBsimHitTool.SimTrackHits.Path = "VertexBarrelCollection"
+
+saveVTXEsimHitTool = SimG4SaveTrackerHits("saveVTXEsimHitTool", readoutNames=["VertexEndcapCollection"])
+saveVTXEsimHitTool.SimTrackHits.Path = "VertexEndcapCollection"
+
+saveITKBsimHitTool = SimG4SaveTrackerHits("saveITKBsimHitTool", readoutNames=["InnerTrackerBarrelCollection"])
+saveITKBsimHitTool.SimTrackHits.Path = "InnerTrackerBarrelCollection"
+
+saveOTKBsimHitTool = SimG4SaveTrackerHits("saveOTKBsimHitTool", readoutNames=["OuterTrackerBarrelCollection"])
+saveOTKBsimHitTool.SimTrackHits.Path = "OuterTrackerBarrelCollection"
+
+saveITKEsimHitTool = SimG4SaveTrackerHits("saveITKEsimHitTool", readoutNames=["InnerTrackerEndcapCollection"])
+saveITKEsimHitTool.SimTrackHits.Path = "InnerTrackerEndcapCollection"
+
+saveOTKEsimHitTool = SimG4SaveTrackerHits("saveOTKEsimHitTool", readoutNames=["OuterTrackerEndcapCollection"])
+saveOTKEsimHitTool.SimTrackHits.Path = "OuterTrackerEndcapCollection"
+
+
+
+# ECAL
+ecalBarrelReadoutName = "ECalBarrelEta"
+ecalBarrelReadoutNamePhiEta = "ECalBarrelPhiEta"
+ecalEndcapReadoutName = "ECalEndcapPhiEta"
+
+# HCAL
+hcalBarrelReadoutName = "HCalBarrelReadout"
+hcalEndcapReadoutName = "HCalEndcapReadout"
+
+# Configure saving of calorimeter hits
+ecalBarrelHitsName = "ECalBarrelPositionedHits"
+from Configurables import SimG4SaveCalHits
+saveECalBarrelTool = SimG4SaveCalHits("saveECalBarrelHits", readoutNames = [ecalBarrelReadoutName])
+saveECalBarrelTool.CaloHits.Path = ecalBarrelHitsName
+
+hcalBarrelHitsName = "HCalBarrelPositionedHits"
+saveHCalTool = SimG4SaveCalHits("saveHCalBarrelHits", readoutNames = [hcalBarrelReadoutName])
+saveHCalTool.CaloHits.Path = hcalBarrelHitsName
+
+ecalEndcapHitsName = "ECalEndcapPositionedHits"
+saveECalEndcapTool = SimG4SaveCalHits("saveECalEndcapHits", readoutNames = [ecalEndcapReadoutName])
+saveECalEndcapTool.CaloHits.Path = ecalEndcapHitsName
+
+hcalEndcapHitsName = "HCalEndcapPositionedHits"
+saveHCalEndcapTool = SimG4SaveCalHits("saveHCalEndcapHits", readoutNames = [hcalEndcapReadoutName])
+saveHCalEndcapTool.CaloHits.Path = hcalEndcapHitsName
+
+# Geant4 algorithm
+# Translates EDM to G4Event, passes the event to G4, writes out outputs via tools
+# and a tool that saves the calorimeter hits
+# create the G4 algorithm, giving the list of names of tools ("XX/YY")
+from Configurables import SimG4PrimariesFromEdmTool
+particle_converter = SimG4PrimariesFromEdmTool("EdmConverter")
+particle_converter.GenParticles.Path = genParticlesOutputName
+
+from Configurables import SimG4Alg
+geantsim = SimG4Alg("SimG4Alg",
+ outputs= [#saveECalBarrelTool,
+ #saveHCalTool,
+ #saveECalEndcapTool,
+ #saveHCalEndcapTool,
+ #saveHistTool
+ saveVTXBsimHitTool, saveVTXEsimHitTool, saveITKBsimHitTool, saveOTKBsimHitTool, saveITKEsimHitTool, saveOTKEsimHitTool
+
+ ],
+ eventProvider=particle_converter,
+ OutputLevel=INFO)
+
+
+# EDM4hep to LCIO converter
+from Configurables import MarlinProcessorWrapper, Lcio2EDM4hepTool, EDM4hep2LcioTool
+MyAIDAProcessor = MarlinProcessorWrapper("MyAIDAProcessor")
+MyAIDAProcessor.OutputLevel = DEBUG
+MyAIDAProcessor.ProcessorType = "AIDAProcessor"
+MyAIDAProcessor.Parameters = {
+ "Compress": ["1"],
+ "FileName": ["histograms"],
+ "FileType": ["root"]
+ }
+edmConvTool = EDM4hep2LcioTool("EDM4hep2lcio")
+edmConvTool.convertAll = False
+edmConvTool.collNameMapping = {
+ 'MCParticles': 'MCParticles',
+ 'VertexBarrelCollection': 'VertexBarrelCollection',
+ 'VertexEndcapCollection': 'VertexEndcapCollection',
+ 'InnerTrackerBarrelCollection': 'InnerTrackerBarrelCollection',
+ 'OuterTrackerBarrelCollection': 'OuterTrackerBarrelCollection',
+ 'InnerTrackerEndcapCollection': 'InnerTrackerEndcapCollection',
+ 'OuterTrackerEndcapCollection': 'OuterTrackerEndcapCollection',
+ #'ECalEndcapCollection': 'ECalEndcapCollection',
+ 'ECalBarrelEta': 'ECalBarrelCollection',
+ #'ECalBarrelCollection': 'ECalBarrelCollection',
+ #'HCalBarrelCollection': 'HCalBarrelCollection',
+ #'HCalEndcapCollection': 'HCalEndcapCollection',
+ #'HCalRingCollection': 'HCalRingCollection',
+ #'YokeBarrelCollection': 'YokeBarrelCollection',
+ #'YokeEndcapCollection': 'YokeEndcapCollection',
+ #'LumiCalCollection': 'LumiCalCollection',
+}
+edmConvTool.OutputLevel = DEBUG
+MyAIDAProcessor.EDM4hep2LcioTool=edmConvTool
+
+
+InitDD4hep = MarlinProcessorWrapper("InitDD4hep")
+InitDD4hep.OutputLevel = WARNING
+InitDD4hep.ProcessorType = "InitializeDD4hep"
+InitDD4hep.Parameters = {
+ "DD4hepXMLFile": ["../../k4geo/FCCee/compact/FCCee_o1_v05/FCCee_o1_v05.xml"],
+ "EncodingStringParameter": ["GlobalTrackerReadoutID"]
+ }
+
+
+VXDBarrelDigitiser = MarlinProcessorWrapper("VXDBarrelDigitiser")
+VXDBarrelDigitiser.OutputLevel = DEBUG
+VXDBarrelDigitiser.ProcessorType = "DDPlanarDigiProcessor"
+VXDBarrelDigitiser.Parameters = {
+ "IsStrip": ["false"],
+ "ResolutionU": ["0.003", "0.003", "0.003", "0.003", "0.003", "0.003"],
+ "ResolutionV": ["0.003", "0.003", "0.003", "0.003", "0.003", "0.003"],
+ "SimTrackHitCollectionName": ["VertexBarrelCollection"],
+ "SimTrkHitRelCollection": ["VXDTrackerHitRelations"],
+ "SubDetectorName": ["Vertex"],
+ "TrackerHitCollectionName": ["VXDTrackerHits"]
+ }
+
+VXDEndcapDigitiser = MarlinProcessorWrapper("VXDEndcapDigitiser")
+VXDEndcapDigitiser.OutputLevel = DEBUG
+VXDEndcapDigitiser.ProcessorType = "DDPlanarDigiProcessor"
+VXDEndcapDigitiser.Parameters = {
+ "IsStrip": ["false"],
+ "ResolutionU": ["0.003", "0.003", "0.003", "0.003", "0.003", "0.003"],
+ "ResolutionV": ["0.003", "0.003", "0.003", "0.003", "0.003", "0.003"],
+ "SimTrackHitCollectionName": ["VertexEndcapCollection"],
+ "SimTrkHitRelCollection": ["VXDEndcapTrackerHitRelations"],
+ "SubDetectorName": ["Vertex"],
+ "TrackerHitCollectionName": ["VXDEndcapTrackerHits"]
+ }
+
+InnerPlanarDigiProcessor = MarlinProcessorWrapper("InnerPlanarDigiProcessor")
+InnerPlanarDigiProcessor.OutputLevel = WARNING
+InnerPlanarDigiProcessor.ProcessorType = "DDPlanarDigiProcessor"
+InnerPlanarDigiProcessor.Parameters = {
+ "IsStrip": ["false"],
+ "ResolutionU": ["0.007"],
+ "ResolutionV": ["0.09"],
+ "SimTrackHitCollectionName": ["InnerTrackerBarrelCollection"],
+ "SimTrkHitRelCollection": ["InnerTrackerBarrelHitsRelations"],
+ "SubDetectorName": ["InnerTrackers"],
+ "TrackerHitCollectionName": ["ITrackerHits"]
+ }
+
+InnerEndcapPlanarDigiProcessor = MarlinProcessorWrapper("InnerEndcapPlanarDigiProcessor")
+InnerEndcapPlanarDigiProcessor.OutputLevel = WARNING
+InnerEndcapPlanarDigiProcessor.ProcessorType = "DDPlanarDigiProcessor"
+InnerEndcapPlanarDigiProcessor.Parameters = {
+ "IsStrip": ["false"],
+ "ResolutionU": ["0.005", "0.007", "0.007", "0.007", "0.007", "0.007", "0.007"],
+ "ResolutionV": ["0.005", "0.09", "0.09", "0.09", "0.09", "0.09", "0.09"],
+ "SimTrackHitCollectionName": ["InnerTrackerEndcapCollection"],
+ "SimTrkHitRelCollection": ["InnerTrackerEndcapHitsRelations"],
+ "SubDetectorName": ["InnerTrackers"],
+ "TrackerHitCollectionName": ["ITrackerEndcapHits"]
+ }
+
+OuterPlanarDigiProcessor = MarlinProcessorWrapper("OuterPlanarDigiProcessor")
+OuterPlanarDigiProcessor.OutputLevel = WARNING
+OuterPlanarDigiProcessor.ProcessorType = "DDPlanarDigiProcessor"
+OuterPlanarDigiProcessor.Parameters = {
+ "IsStrip": ["false"],
+ "ResolutionU": ["0.007", "0.007", "0.007"],
+ "ResolutionV": ["0.09", "0.09", "0.09"],
+ "SimTrackHitCollectionName": ["OuterTrackerBarrelCollection"],
+ "SimTrkHitRelCollection": ["OuterTrackerBarrelHitsRelations"],
+ "SubDetectorName": ["OuterTrackers"],
+ "TrackerHitCollectionName": ["OTrackerHits"]
+ }
+
+OuterEndcapPlanarDigiProcessor = MarlinProcessorWrapper("OuterEndcapPlanarDigiProcessor")
+OuterEndcapPlanarDigiProcessor.OutputLevel = WARNING
+OuterEndcapPlanarDigiProcessor.ProcessorType = "DDPlanarDigiProcessor"
+OuterEndcapPlanarDigiProcessor.Parameters = {
+ "IsStrip": ["false"],
+ "ResolutionU": ["0.007", "0.007", "0.007", "0.007", "0.007"],
+ "ResolutionV": ["0.09", "0.09", "0.09", "0.09", "0.09"],
+ "SimTrackHitCollectionName": ["OuterTrackerEndcapCollection"],
+ "SimTrkHitRelCollection": ["OuterTrackerEndcapHitsRelations"],
+ "SubDetectorName": ["OuterTrackers"],
+ "TrackerHitCollectionName": ["OTrackerEndcapHits"]
+ }
+
+MyTruthTrackFinder = MarlinProcessorWrapper("MyTruthTrackFinder")
+MyTruthTrackFinder.OutputLevel = WARNING
+MyTruthTrackFinder.ProcessorType = "TruthTrackFinder"
+MyTruthTrackFinder.Parameters = {
+ "FitForward": ["true"],
+ "MCParticleCollectionName": ["MCParticle"],
+ "SiTrackCollectionName": ["SiTracks"],
+ "SiTrackRelationCollectionName": ["SiTrackRelations"],
+ "SimTrackerHitRelCollectionNames": ["VXDTrackerHitRelations", "InnerTrackerBarrelHitsRelations", "OuterTrackerBarrelHitsRelations", "VXDEndcapTrackerHitRelations", "InnerTrackerEndcapHitsRelations", "OuterTrackerEndcapHitsRelations"],
+ "TrackerHitCollectionNames": ["VXDTrackerHits", "ITrackerHits", "OTrackerHits", "VXDEndcapTrackerHits", "ITrackerEndcapHits", "OTrackerEndcapHits"],
+ "UseTruthInPrefit": ["false"]
+ }
+
+MyConformalTracking = MarlinProcessorWrapper("MyConformalTracking")
+MyConformalTracking.OutputLevel = DEBUG
+MyConformalTracking.ProcessorType = "ConformalTrackingV2"
+MyConformalTracking.Parameters = {
+ "DebugHits": ["DebugHits"],
+ "DebugPlots": ["false"],
+ "DebugTiming": ["false"],
+ "MCParticleCollectionName": ["MCParticle"],
+ "MaxHitInvertedFit": ["0"],
+ "MinClustersOnTrackAfterFit": ["3"],
+ "RelationsNames": ["VXDTrackerHitRelations", "VXDEndcapTrackerHitRelations", "InnerTrackerBarrelHitsRelations", "OuterTrackerBarrelHitsRelations", "InnerTrackerEndcapHitsRelations", "OuterTrackerEndcapHitsRelations"],
+ "RetryTooManyTracks": ["false"],
+ "SiTrackCollectionName": ["SiTracksCT"],
+ "SortTreeResults": ["true"],
+ "Steps": ["[VXDBarrel]", "@Collections", ":", "VXDTrackerHits", "@Parameters", ":", "MaxCellAngle", ":", "0.01;", "MaxCellAngleRZ", ":", "0.01;", "Chi2Cut", ":", "100;", "MinClustersOnTrack", ":", "4;", "MaxDistance", ":", "0.03;", "SlopeZRange:", "10.0;", "HighPTCut:", "10.0;", "@Flags", ":", "HighPTFit,", "VertexToTracker", "@Functions", ":", "CombineCollections,", "BuildNewTracks", "[VXDEncap]", "@Collections", ":", "VXDEndcapTrackerHits", "@Parameters", ":", "MaxCellAngle", ":", "0.01;", "MaxCellAngleRZ", ":", "0.01;", "Chi2Cut", ":", "100;", "MinClustersOnTrack", ":", "4;", "MaxDistance", ":", "0.03;", "SlopeZRange:", "10.0;", "HighPTCut:", "10.0;", "@Flags", ":", "HighPTFit,", "VertexToTracker", "@Functions", ":", "CombineCollections,", "ExtendTracks", "[LowerCellAngle1]", "@Collections", ":", "VXDTrackerHits,", "VXDEndcapTrackerHits", "@Parameters", ":", "MaxCellAngle", ":", "0.05;", "MaxCellAngleRZ", ":", "0.05;", "Chi2Cut", ":", "100;", "MinClustersOnTrack", ":", "4;", "MaxDistance", ":", "0.03;", "SlopeZRange:", "10.0;HighPTCut:", "10.0;", "@Flags", ":", "HighPTFit,", "VertexToTracker,", "RadialSearch", "@Functions", ":", "CombineCollections,", "BuildNewTracks", "[LowerCellAngle2]", "@Collections", ":", "@Parameters", ":", "MaxCellAngle", ":", "0.1;", "MaxCellAngleRZ", ":", "0.1;", "Chi2Cut", ":", "2000;", "MinClustersOnTrack", ":", "4;", "MaxDistance", ":", "0.03;", "SlopeZRange:", "10.0;", "HighPTCut:", "10.0;", "@Flags", ":", "HighPTFit,", "VertexToTracker,", "RadialSearch", "@Functions", ":", "BuildNewTracks,", "SortTracks", "[Tracker]", "@Collections", ":", "ITrackerHits,", "OTrackerHits,", "ITrackerEndcapHits,", "OTrackerEndcapHits", "@Parameters", ":", "MaxCellAngle", ":", "0.1;", "MaxCellAngleRZ", ":", "0.1;", "Chi2Cut", ":", "2000;", "MinClustersOnTrack", ":", "4;", "MaxDistance", ":", "0.03;", "SlopeZRange:", "10.0;", "HighPTCut:", "1.0;", "@Flags", ":", "HighPTFit,", "VertexToTracker,", "RadialSearch", "@Functions", ":", "CombineCollections,", "ExtendTracks", "[Displaced]", "@Collections", ":", "VXDTrackerHits,", "VXDEndcapTrackerHits,", "ITrackerHits,", "OTrackerHits,", "ITrackerEndcapHits,", "OTrackerEndcapHits", "@Parameters", ":", "MaxCellAngle", ":", "0.1;", "MaxCellAngleRZ", ":", "0.1;", "Chi2Cut", ":", "1000;", "MinClustersOnTrack", ":", "5;", "MaxDistance", ":", "0.015;", "SlopeZRange:", "10.0;", "HighPTCut:", "10.0;", "@Flags", ":", "OnlyZSchi2cut,", "RadialSearch", "@Functions", ":", "CombineCollections,", "BuildNewTracks"],
+ "ThetaRange": ["0.05"],
+ "TooManyTracks": ["100000"],
+ "TrackerHitCollectionNames": ["VXDTrackerHits", "VXDEndcapTrackerHits", "ITrackerHits", "OTrackerHits", "ITrackerEndcapHits", "OTrackerEndcapHits"],
+ "trackPurity": ["0.7"]
+ }
+
+ClonesAndSplitTracksFinder = MarlinProcessorWrapper("ClonesAndSplitTracksFinder")
+ClonesAndSplitTracksFinder.OutputLevel = WARNING
+ClonesAndSplitTracksFinder.ProcessorType = "ClonesAndSplitTracksFinder"
+ClonesAndSplitTracksFinder.Parameters = {
+ "EnergyLossOn": ["true"],
+ "InputTrackCollectionName": ["SiTracksCT"],
+ "MultipleScatteringOn": ["true"],
+ "OutputTrackCollectionName": ["SiTracks"],
+ "SmoothOn": ["false"],
+ "extrapolateForward": ["true"],
+ "maxSignificancePhi": ["3"],
+ "maxSignificancePt": ["2"],
+ "maxSignificanceTheta": ["3"],
+ "mergeSplitTracks": ["false"],
+ "minTrackPt": ["1"]
+ }
+
+Refit = MarlinProcessorWrapper("Refit")
+Refit.OutputLevel = WARNING
+Refit.ProcessorType = "RefitFinal"
+Refit.Parameters = {
+ "EnergyLossOn": ["true"],
+ "InputRelationCollectionName": ["SiTrackRelations"],
+ "InputTrackCollectionName": ["SiTracks"],
+ "Max_Chi2_Incr": ["1.79769e+30"],
+ "MinClustersOnTrackAfterFit": ["3"],
+ "MultipleScatteringOn": ["true"],
+ "OutputRelationCollectionName": ["SiTracks_Refitted_Relation"],
+ "OutputTrackCollectionName": ["SiTracks_Refitted"],
+ "ReferencePoint": ["-1"],
+ "SmoothOn": ["false"],
+ "extrapolateForward": ["true"]
+ }
+
+MyClicEfficiencyCalculator = MarlinProcessorWrapper("MyClicEfficiencyCalculator")
+MyClicEfficiencyCalculator.OutputLevel = WARNING
+MyClicEfficiencyCalculator.ProcessorType = "ClicEfficiencyCalculator"
+MyClicEfficiencyCalculator.Parameters = {
+ "MCParticleCollectionName": ["MCParticle"],
+ "MCParticleNotReco": ["MCParticleNotReco"],
+ "MCPhysicsParticleCollectionName": ["MCPhysicsParticles"],
+ "TrackCollectionName": ["SiTracks_Refitted"],
+ "TrackerHitCollectionNames": ["VXDTrackerHits", "VXDEndcapTrackerHits", "ITrackerHits", "OTrackerHits", "ITrackerEndcapHits", "OTrackerEndcapHits"],
+ "TrackerHitRelCollectionNames": ["VXDTrackerHitRelations", "VXDEndcapTrackerHitRelations", "InnerTrackerBarrelHitsRelations", "OuterTrackerBarrelHitsRelations", "InnerTrackerEndcapHitsRelations", "OuterTrackerEndcapHitsRelations"],
+ "efficiencyTreeName": ["trktree"],
+ "mcTreeName": ["mctree"],
+ "morePlots": ["false"],
+ "purityTreeName": ["puritytree"],
+ "reconstructableDefinition": ["ILDLike"],
+ "vertexBarrelID": ["1"]
+ }
+
+MyTrackChecker = MarlinProcessorWrapper("MyTrackChecker")
+MyTrackChecker.OutputLevel = WARNING
+MyTrackChecker.ProcessorType = "TrackChecker"
+MyTrackChecker.Parameters = {
+ "MCParticleCollectionName": ["MCParticle"],
+ "TrackCollectionName": ["SiTracks_Refitted"],
+ "TrackRelationCollectionName": ["SiTracksMCTruthLink"],
+ "TreeName": ["checktree"],
+ "UseOnlyTree": ["true"]
+ }
+
+EventNumber = MarlinProcessorWrapper("EventNumber")
+EventNumber.OutputLevel = WARNING
+EventNumber.ProcessorType = "Statusmonitor"
+EventNumber.Parameters = {
+ "HowOften": ["1"]
+ }
+
+Output_REC = MarlinProcessorWrapper("Output_REC")
+Output_REC.OutputLevel = WARNING
+Output_REC.ProcessorType = "LCIOOutputProcessor"
+Output_REC.Parameters = {
+ "DropCollectionNames": [],
+ "DropCollectionTypes": [],
+ "FullSubsetCollections": ["EfficientMCParticles", "InefficientMCParticles"],
+ "KeepCollectionNames": [],
+ "LCIOOutputFile": ["Output_REC.slcio"],
+ "LCIOWriteMode": ["WRITE_NEW"]
+ }
+
+# LCIO to EDM4hep converter
+lcioConvTool = Lcio2EDM4hepTool("lcio2EDM4hep")
+lcioConvTool.convertAll = True
+lcioConvTool.OutputLevel = DEBUG
+Output_REC.Lcio2EDM4hepTool=lcioConvTool
+
+Output_DST = MarlinProcessorWrapper("Output_DST")
+Output_DST.OutputLevel = WARNING
+Output_DST.ProcessorType = "LCIOOutputProcessor"
+Output_DST.Parameters = {
+ "DropCollectionNames": [],
+ "DropCollectionTypes": ["MCParticle", "LCRelation", "SimCalorimeterHit", "CalorimeterHit", "SimTrackerHit", "TrackerHit", "TrackerHitPlane", "Track", "ReconstructedParticle", "LCFloatVec"],
+ "FullSubsetCollections": ["EfficientMCParticles", "InefficientMCParticles", "MCPhysicsParticles"],
+ "KeepCollectionNames": ["MCParticlesSkimmed", "MCPhysicsParticles", "RecoMCTruthLink", "SiTracks", "SiTracks_Refitted", "PandoraClusters", "PandoraPFOs", "SelectedPandoraPFOs", "LooseSelectedPandoraPFOs", "TightSelectedPandoraPFOs", "RefinedVertexJets", "RefinedVertexJets_rel", "RefinedVertexJets_vtx", "RefinedVertexJets_vtx_RP", "BuildUpVertices", "BuildUpVertices_res", "BuildUpVertices_RP", "BuildUpVertices_res_RP", "BuildUpVertices_V0", "BuildUpVertices_V0_res", "BuildUpVertices_V0_RP", "BuildUpVertices_V0_res_RP", "PrimaryVertices", "PrimaryVertices_res", "PrimaryVertices_RP", "PrimaryVertices_res_RP", "RefinedVertices", "RefinedVertices_RP"],
+ "LCIOOutputFile": ["Output_DST.slcio"],
+ "LCIOWriteMode": ["WRITE_NEW"]
+ }
+
+
+############## Digitization (Merging hits into cells, EM scale calibration)
+# EM scale calibration (sampling fraction)
+from Configurables import CalibrateInLayersTool
+calibEcalBarrel = CalibrateInLayersTool("CalibrateECalBarrel",
+ samplingFraction = [0.3632447480841956] * 1 + [0.13187261040190248] * 1 + [0.14349714292943705] * 1 + [0.150266118277841] * 1 + [0.15502683375826457] * 1 + [0.15954408786354762] * 1 + [0.16375302347299436] * 1 + [0.16840384714588075] * 1 + [0.17219540619311383] * 1 + [0.1755068643940401] * 1 + [0.17816980262822366] * 1 + [0.18131266048670405] * 1,
+ readoutName = ecalBarrelReadoutName,
+ layerFieldName = "layer")
+
+from Configurables import CalibrateCaloHitsTool
+calibHcells = CalibrateCaloHitsTool("CalibrateHCal", invSamplingFraction="41.66")
+calibEcalEndcap = CalibrateCaloHitsTool("CalibrateECalEndcap", invSamplingFraction="4.27")
+calibHcalEndcap = CalibrateCaloHitsTool("CalibrateHCalEndcap", invSamplingFraction="31.7")
+
+# Create cells in ECal barrel
+# 1. step - merge hits into cells with Eta and module segmentation (phi module is a 'physical' cell i.e. lead + LAr + PCB + LAr +lead)
+# 2. step - rewrite the cellId using the Eta-Phi segmentation (merging several modules into one phi readout cell). Add noise at this step if you derived the noise already assuming merged cells
+from Configurables import CreateCaloCells
+createEcalBarrelCellsStep1 = CreateCaloCells("CreateECalBarrelCellsStep1",
+ doCellCalibration=True,
+ calibTool = calibEcalBarrel,
+ addCellNoise=False, filterCellNoise=False,
+ addPosition=True,
+ OutputLevel=INFO,
+ hits=ecalBarrelHitsName,
+ cells="ECalBarrelCellsStep1")
+
+## Use Phi-Theta segmentation in ECal barrel
+from Configurables import RedoSegmentation
+resegmentEcalBarrel = RedoSegmentation("ReSegmentationEcal",
+ # old bitfield (readout)
+ oldReadoutName = ecalBarrelReadoutName,
+ # specify which fields are going to be altered (deleted/rewritten)
+ oldSegmentationIds = ["module"],
+ # new bitfield (readout), with new segmentation
+ newReadoutName = ecalBarrelReadoutNamePhiEta,
+ OutputLevel = INFO,
+ inhits = "ECalBarrelCellsStep1",
+ outhits = "ECalBarrelCellsStep2")
+
+# cells without noise
+EcalBarrelCellsName = "ECalBarrelCells"
+createEcalBarrelCells = CreateCaloCells("CreateECalBarrelCells",
+ doCellCalibration=False,
+ addCellNoise=False, filterCellNoise=False,
+ OutputLevel=INFO,
+ hits="ECalBarrelCellsStep2",
+ cells=EcalBarrelCellsName)
+cell_creator_to_use = createEcalBarrelCells
+
+# generate noise for each cell
+if addNoise:
+ ecalBarrelNoisePath = os.environ['FCCBASEDIR']+"/LAr_scripts/data/elecNoise_ecalBarrelFCCee.root"
+ ecalBarrelNoiseHistName = "h_elecNoise_fcc_"
+ from Configurables import NoiseCaloCellsFromFileTool
+ noiseBarrel = NoiseCaloCellsFromFileTool("NoiseBarrel",
+ readoutName = ecalBarrelReadoutNamePhiEta,
+ noiseFileName = ecalBarrelNoisePath,
+ elecNoiseHistoName = ecalBarrelNoiseHistName,
+ activeFieldName = "layer",
+ addPileup = False,
+ filterNoiseThreshold = 2,
+ scaleFactor = 1/1000.,
+ numRadialLayers = 12)
+
+ from Configurables import TubeLayerPhiEtaCaloTool
+ barrelGeometry = TubeLayerPhiEtaCaloTool("EcalBarrelGeo",
+ readoutName = ecalBarrelReadoutNamePhiEta,
+ activeVolumeName = "LAr_sensitive",
+ activeFieldName = "layer",
+ fieldNames = ["system"],
+ fieldValues = [4])
+ #activeVolumesNumber = 12)
+ # cells with noise not filtered
+ createEcalBarrelCellsNoise = CreateCaloCells("CreateECalBarrelCellsNoise",
+ doCellCalibration=False,
+ addCellNoise=True, filterCellNoise=False,
+ OutputLevel=INFO,
+ hits="ECalBarrelCellsStep2",
+ noiseTool = noiseBarrel,
+ geometryTool = barrelGeometry,
+ cells=EcalBarrelCellsName)
+
+ # cells with noise filtered
+ #createEcalBarrelCellsNoise = CreateCaloCells("CreateECalBarrelCellsNoise_filtered",
+ # doCellCalibration=False,
+ # addCellNoise=True, filterCellNoise=True,
+ # OutputLevel=INFO,
+ # hits="ECalBarrelCellsStep2",
+ # noiseTool = noiseBarrel,
+ # geometryTool = barrelGeometry,
+ # cells=EcalBarrelCellsName)
+
+ cell_creator_to_use = createEcalBarrelCellsNoise
+
+
+# Ecal barrel cell positions (good for physics, all coordinates set properly)
+from Configurables import CellPositionsECalBarrelTool
+cellPositionEcalBarrelTool = CellPositionsECalBarrelTool("CellPositionsECalBarrel", readoutName = ecalBarrelReadoutNamePhiEta, OutputLevel = INFO)
+
+from Configurables import CreateCaloCellPositionsFCCee
+createEcalBarrelPositionedCells = CreateCaloCellPositionsFCCee("ECalBarrelPositionedCells", OutputLevel = INFO)
+createEcalBarrelPositionedCells.systemIdECalBarrel = 20
+createEcalBarrelPositionedCells.positionsECalBarrelTool = cellPositionEcalBarrelTool
+createEcalBarrelPositionedCells.hits.Path = EcalBarrelCellsName
+createEcalBarrelPositionedCells.positionedHits.Path = "ECalBarrelPositionedCells"
+
+
+# Create cells in HCal
+# 1. step - merge hits into cells with the default readout
+HcalBarrelCellsName = "HCalBarrelCells"
+createHcalBarrelCells = CreateCaloCells("CreateHCaloCells",
+ doCellCalibration=True,
+ calibTool=calibHcells,
+ addCellNoise = False,
+ filterCellNoise = False,
+ addPosition=True,
+ OutputLevel = INFO,
+ hits=hcalBarrelHitsName,
+ cells=HcalBarrelCellsName)
+
+HcalEndcapCellsName = "HCalEndcapCells"
+createHcalEndcapCells = CreateCaloCells("CreateHcalEndcapCaloCells",
+ doCellCalibration=True,
+ calibTool=calibHcalEndcap,
+ addCellNoise=False,
+ filterCellNoise=False,
+ OutputLevel=INFO,
+ hits=hcalEndcapHitsName,
+ cells=HcalEndcapCellsName)
+
+"""
+from Configurables import CellPositionsHCalBarrelTool
+cellPositionHcalBarrelTool = CellPositionsHCalBarrelTool("CellPositionsHCalBarrel", readoutName = hcalBarrelReadoutName, OutputLevel = INFO)
+createHcalBarrelPositionedCells = CreateCaloCellPositionsFCCee("HCalBarrelPositionedCells", OutputLevel = INFO)
+createHcalBarrelPositionedCells.positionsHCalBarrelTool = cellPositionHcalBarrelTool
+createHcalBarrelPositionedCells.hits.Path = HcalBarrelCellsName
+HCalBarrelPositionedCellsName = "HCalBarrelPositionedCellsName"
+createHcalBarrelPositionedCells.positionedHits.Path = "HCalBarrelPositionedCells"
+"""
+EcalEndcapCellsName = "ECalEndcapCells"
+createEcalEndcapCells = CreateCaloCells("CreateEcalEndcapCaloCells",
+ doCellCalibration=True,
+ calibTool=calibEcalEndcap,
+ addCellNoise=False,
+ filterCellNoise=False,
+ addPosition=True,
+ OutputLevel = INFO,
+ hits=ecalEndcapHitsName,
+ cells=EcalEndcapCellsName)
+
+#Empty cells for parts of calorimeter not implemented yet
+from Configurables import CreateEmptyCaloCellsCollection
+createemptycells = CreateEmptyCaloCellsCollection("CreateEmptyCaloCells")
+createemptycells.cells.Path = "emptyCaloCells"
+
+from Configurables import CaloTowerTool
+towers = CaloTowerTool("towers",
+ deltaEtaTower = 0.01, deltaPhiTower = 2*_pi/768.,
+ radiusForPosition = 2160 + 40 / 2.0,
+ ecalBarrelReadoutName = ecalBarrelReadoutNamePhiEta,
+ ecalEndcapReadoutName = ecalEndcapReadoutName,
+ ecalFwdReadoutName = "",
+ hcalBarrelReadoutName = hcalBarrelReadoutName,
+ hcalExtBarrelReadoutName = "",
+ hcalEndcapReadoutName = "",
+ hcalFwdReadoutName = "",
+ OutputLevel = INFO)
+towers.ecalBarrelCells.Path = EcalBarrelCellsName
+towers.ecalEndcapCells.Path = EcalEndcapCellsName
+towers.ecalFwdCells.Path = "emptyCaloCells"
+towers.hcalBarrelCells.Path = HcalBarrelCellsName
+towers.hcalExtBarrelCells.Path = "emptyCaloCells"
+towers.hcalEndcapCells.Path = "emptyCaloCells"
+towers.hcalFwdCells.Path = "emptyCaloCells"
+
+# Cluster variables
+windE = 9
+windP = 17
+posE = 5
+posP = 11
+dupE = 7
+dupP = 13
+finE = 9
+finP = 17
+# Minimal energy to create a cluster in GeV (FCC-ee detectors have to reconstruct low energy particles)
+threshold = 0.05
+
+from Configurables import CreateCaloClustersSlidingWindow
+createClusters = CreateCaloClustersSlidingWindow("CreateClusters",
+ towerTool = towers,
+ nEtaWindow = windE, nPhiWindow = windP,
+ nEtaPosition = posE, nPhiPosition = posP,
+ nEtaDuplicates = dupE, nPhiDuplicates = dupP,
+ nEtaFinal = finE, nPhiFinal = finP,
+ energyThreshold = threshold,
+ energySharingCorrection = True,
+ attachCells = True,
+ OutputLevel = INFO
+ )
+createClusters.clusters.Path = "CaloClusters"
+createClusters.clusterCells.Path = "CaloClusterCells"
+
+createEcalBarrelPositionedCaloClusterCells = CreateCaloCellPositionsFCCee("ECalBarrelPositionedCaloClusterCells", OutputLevel = INFO)
+createEcalBarrelPositionedCaloClusterCells.positionsECalBarrelTool = cellPositionEcalBarrelTool
+createEcalBarrelPositionedCaloClusterCells.hits.Path = "CaloClusterCells"
+createEcalBarrelPositionedCaloClusterCells.positionedHits.Path = "PositionedCaloClusterCells"
+
+from Configurables import CorrectCaloClusters
+correctCaloClusters = CorrectCaloClusters("correctCaloClusters",
+ inClusters = createClusters.clusters.Path,
+ outClusters = "Corrected"+createClusters.clusters.Path,
+ numLayers = [12],
+ firstLayerIDs = [0],
+ lastLayerIDs = [11],
+ readoutNames = [ecalBarrelReadoutNamePhiEta],
+ upstreamParameters = [[0.033955208567442975, -3.818122686176795, -146.59497297249345, 0.563447903447204, -3.7906629536351906, -8.569962044554627]],
+ upstreamFormulas = [['[0]+[1]/(x-[2])', '[0]+[1]/(x-[2])']],
+ downstreamParameters = [[-0.00357017357914002, 0.006624434345822984, 1.0065650241358008, -1.285181650875406, -0.0707783194915608, 12.907319280196257]],
+ downstreamFormulas = [['[0]+[1]*x', '[0]+[1]/sqrt(x)', '[0]+[1]/x']],
+ OutputLevel = INFO
+ )
+
+################ Output
+from Configurables import PodioOutput
+out = PodioOutput("out",
+ OutputLevel=INFO)
+#out.outputCommands = ["drop *", "keep ECalBarrelPositionedCells", "keep gen*"]
+#out.outputCommands = ["keep *", "drop ECalBarrelHits", "drop HCal*", "drop ECalBarrelCellsStep*", "drop ECalBarrelPositionedHits", "drop emptyCaloCells", "drop CaloClusterCells"]
+out.outputCommands = ["keep *", "drop ECalBarrelHits", "drop ECalBarrelCellsStep*", "drop emptyCaloCells"]
+#out.outputCommands = ["keep *", "drop ECalBarrelHits", "drop HCal*", "drop ECalBarrelCellsStep*", "drop ECalBarrelPositionedHits", "drop emptyCaloCells", "drop CaloClusterCells", "drop %s"%EcalBarrelCellsName, "drop %s"%createEcalBarrelPositionedCells.positionedHits.Path]
+
+import uuid
+out.filename = "track_from_ddsim_MagneticField_"+str(magneticField)+"_pMin_"+str(momentum*1000)+"_MeV"+"_ThetaMinMax_"+str(thetaMin)+"_"+str(thetaMax)+"_pdgId_"+str(pdgCode)+"_pythia"+str(use_pythia)+"_Noise"+str(addNoise)+".root"
+
+#CPU information
+from Configurables import AuditorSvc, ChronoAuditor
+chra = ChronoAuditor()
+audsvc = AuditorSvc()
+audsvc.Auditors = [chra]
+genAlg.AuditExecute = True
+hepmc_converter.AuditExecute = True
+geantsim.AuditExecute = True
+createEcalBarrelCellsStep1.AuditExecute = True
+resegmentEcalBarrel.AuditExecute = True
+cell_creator_to_use.AuditExecute = True
+#createHcalBarrelCells.AuditExecute = True
+out.AuditExecute = True
+
+from Configurables import EventCounter
+event_counter = EventCounter('event_counter')
+event_counter.Frequency = 10
+
+from Configurables import ApplicationMgr
+ApplicationMgr(
+ TopAlg = [
+ event_counter,
+ inp,
+ #genAlg,
+ #hepmc_converter,
+ #geantsim,
+ MyAIDAProcessor,
+ #InitDD4hep,
+ VXDBarrelDigitiser,
+ VXDEndcapDigitiser,
+ InnerPlanarDigiProcessor,
+ InnerEndcapPlanarDigiProcessor,
+ OuterPlanarDigiProcessor,
+ OuterEndcapPlanarDigiProcessor,
+ MyConformalTracking,
+ ClonesAndSplitTracksFinder,
+ Refit,
+ Output_REC,
+ Output_DST,
+ #createEcalBarrelCellsStep1,
+ #resegmentEcalBarrel,
+ ##createEcalBarrelCells,
+ #cell_creator_to_use,
+ #createEcalBarrelPositionedCells,
+ #createHcalBarrelCells,
+ ##createHcalBarrelPositionedCells,
+ #createEcalEndcapCells,
+ #createHcalEndcapCells,
+ #createemptycells,
+ #createClusters,
+ #createEcalBarrelPositionedCaloClusterCells,
+ #correctCaloClusters,
+ out
+ ],
+ EvtSel = 'NONE',
+ EvtMax = 10,
+ ExtSvc = [evtsvc, geoservice, podioevent, geantservice, audsvc],
+ StopOnSignal = True,
+ )