From 54e6a3566a23fc242c404ce7e72f0f82b0f021d8 Mon Sep 17 00:00:00 2001
From: Shaojun Lu <shaojun.lu@desy.de>
Date: Mon, 29 Feb 2016 15:30:49 +0000
Subject: [PATCH]  Added CALICE testbeam 'TBModel2015' compact, material and
 steering files.

---
 .../CaloTB_EPT_AHCAL/TBModel2015_steering.py  |  38 ++++++
 .../CaloTB_EPT_AHCAL/compact/TBModel2015.xml  | 110 +++++++++++++++++
 CaloTB/CaloTB_EPT_AHCAL/compact/TBecal4d.xml  |  30 +++++
 CaloTB/CaloTB_EPT_AHCAL/compact/TBhcal4d.xml  | 111 ++++++++++++++++++
 .../compact/materials_TBModel2015.xml         | 100 ++++++++++++++++
 5 files changed, 389 insertions(+)
 create mode 100644 CaloTB/CaloTB_EPT_AHCAL/TBModel2015_steering.py
 create mode 100644 CaloTB/CaloTB_EPT_AHCAL/compact/TBModel2015.xml
 create mode 100644 CaloTB/CaloTB_EPT_AHCAL/compact/TBecal4d.xml
 create mode 100644 CaloTB/CaloTB_EPT_AHCAL/compact/TBhcal4d.xml
 create mode 100644 CaloTB/CaloTB_EPT_AHCAL/compact/materials_TBModel2015.xml

diff --git a/CaloTB/CaloTB_EPT_AHCAL/TBModel2015_steering.py b/CaloTB/CaloTB_EPT_AHCAL/TBModel2015_steering.py
new file mode 100644
index 000000000..b547c31e3
--- /dev/null
+++ b/CaloTB/CaloTB_EPT_AHCAL/TBModel2015_steering.py
@@ -0,0 +1,38 @@
+from DDSim.DD4hepSimulation import DD4hepSimulation
+from SystemOfUnits import mm, GeV, MeV
+
+SIM = DD4hepSimulation()
+
+SIM.compactFile = "../../DD4HEP/compact/TBModel2015.xml"
+SIM.runType = "batch"
+SIM.macroFile = "vis.mac"
+#SIM.inputFiles = "mcparticles.slcio"
+SIM.outputFile = "DD4hep_mu-_8GeV_QGSP_BERT_10k.slcio"
+
+SIM.numberOfEvents = 10000
+SIM.skipNEvents = 0
+SIM.physicsList = "QGSP_BERT"
+SIM.dumpSteeringFile = "TBModel2015_dump.xml"
+SIM.enableDetailedShowerMode=True
+
+SIM.random.seed = "0123456789"
+SIM.field.eps_min = 1*mm
+SIM.part.minimalKineticEnergy = 1*MeV
+
+SIM.action.calo = "Geant4ScintillatorCalorimeterAction"
+
+## set the particle.tbl file to add extra particles to DDsim (B-Baryons)
+## use the power of python to get the file from DD4hep wherever it is
+import os
+if os.path.exists( os.path.join( os.environ.get("DD4hepINSTALL"), "examples/DDG4/examples/particle.tbl") ):
+  SIM.physics.pdgfile = os.path.join( os.environ.get("DD4hepINSTALL"), "examples/DDG4/examples/particle.tbl")
+
+
+SIM.enableGun = True
+SIM.gun.particle = "mu-"
+SIM.gun.energy = 8*GeV
+SIM.gun.position = "0, 0, -1000"
+SIM.gun.direction = "0,0,1"
+
+#SIM.gun.isotrop
+SIM.gun.multiplicity = 1
diff --git a/CaloTB/CaloTB_EPT_AHCAL/compact/TBModel2015.xml b/CaloTB/CaloTB_EPT_AHCAL/compact/TBModel2015.xml
new file mode 100644
index 000000000..540c0adaf
--- /dev/null
+++ b/CaloTB/CaloTB_EPT_AHCAL/compact/TBModel2015.xml
@@ -0,0 +1,110 @@
+<lccdd xmlns:compact="http://www.lcsim.org/schemas/compact/1.0"
+       xmlns:xs="http://www.w3.org/2001/XMLSchema"
+       xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/compact/1.0/compact.xsd">
+  
+  <info name="TBhcal4d"
+        title="TBhcal4d"
+        author="Eldwan Brianne"
+        url="flc.desy.de"
+        status="development"
+        version="$Id$">
+    <comment>The compact File for the AHCAL TB 2015 Prototype</comment>
+  </info>
+  
+  <includes>
+    <gdmlFile  ref="${DD4hepINSTALL}/DDDetectors/compact/elements.xml"/>
+    <comment>Additional design specific material definitions</comment>
+    <gdmlFile  ref="materials_TBModel2015.xml"/>
+    <gdmlFile  ref="${DD4hepINSTALL}/DDDetectors/compact/materials.xml"/>
+  </includes>
+
+  <define>
+    <!-- world -->
+    <constant name="world_side" value="3000*mm"/>
+    <constant name="world_x" value="world_side"/>
+    <constant name="world_y" value="world_side"/>
+    <constant name="world_z" value="world_side"/>
+
+    <!-- Materials -->
+    <constant name="Hcal_radiator_thickness" value="20*mm"/>
+    <constant name="airgap_thickness" value="0.5*mm"/>
+    <constant name="cassette_thickness" value="0.5*mm"/>
+    <constant name="cable_mix_thickness" value="1.5*mm"/>
+    <constant name="pcb_thickness" value="0.7*mm"/>
+    <constant name="foil_thickness" value="0.115*mm"/>
+    <constant name="Hcal_scintillator_thickness" value="3*mm"/>
+
+    <constant name="Ecal_radiator_thickness" value="20*mm"/>
+    <constant name="Ecal_scintillator_thickness" value="2*mm"/>
+
+    <constant name="env_safety" value="0.01*mm"/>
+
+    <!-- ECAL -->
+    <constant name="ECAL_nlayers" value="2"/>
+    <constant name="Ecal_layer_thickness" value="Ecal_radiator_thickness+2*airgap_thickness+2*cassette_thickness+cable_mix_thickness+pcb_thickness+2*foil_thickness+Ecal_scintillator_thickness"/>
+    <constant name="EBU_dim_x" value="18*cm"/>
+    <constant name="EBU_dim_y" value="18*cm"/>
+    <constant name="ECAL_ncell_x" value="4"/>
+    <constant name="ECAL_ncell_y" value="18"/>
+    <constant name="Ecal_cell_size_x" value="45*mm"/>
+    <constant name="Ecal_cell_size_y" value="5*mm"/>
+
+    <!-- HCAL -->
+    <constant name="HCAL_SSF_nlayers" value="8"/>
+    <constant name="HCAL_BL_nlayers" value="22"/>
+    <constant name="Hcal_layer_thickness" value="Hcal_radiator_thickness+2*airgap_thickness+2*cassette_thickness+cable_mix_thickness+pcb_thickness+2*foil_thickness+Hcal_scintillator_thickness"/>
+    <constant name="HBU_SSF_dim_x" value="36*cm"/>
+    <constant name="HBU_SSF_dim_y" value="36*cm"/>
+    <constant name="HCAL_SSF_ncell_x" value="12"/>
+    <constant name="HCAL_SSF_ncell_y" value="12"/>
+    <constant name="HBU_BL_dim_x" value="72*cm"/>
+    <constant name="HBU_BL_dim_y" value="72*cm"/>
+    <constant name="HCAL_BL_ncell_x" value="24"/>
+    <constant name="HCAL_BL_ncell_y" value="24"/>
+    <constant name="Hcal_cell_size_x" value="30*mm"/>
+    <constant name="Hcal_cell_size_y" value="30*mm"/>
+  </define>
+  
+  
+  <limits>
+    <limitset name="cal_limits">
+      <limit name="step_length_max" particles="*" value="5.0" unit="mm" />
+    </limitset>
+  </limits>
+  
+  <comment>Common Generic visualization attributes</comment>
+  <display>
+    <vis name="WhiteVis"         alpha="0.0" r=".96" g=".96"  b=".96"   showDaughters="true"  visible="true"/>
+    <vis name="LightGrayVis"     alpha="0.0" r=".75" g=".75"  b=".75"   showDaughters="true"  visible="true"/>
+    <vis name="Invisible"        alpha="0.0" r="0" g="0"  b="0"         showDaughters="false"  visible="false"/>
+    <vis name="SeeThrough"       alpha="0.0" r="0" g="0"  b="0"         showDaughters="true"  visible="false"/>
+    <vis name="RedVis"           alpha="1.0" r="1.0" g="0.0"  b="0.0"   showDaughters="true"  visible="true"/>
+    <vis name="GreenVis"         alpha="1.0" r="0.0" g="1.0"  b="0.0"   showDaughters="true"  visible="true"/>
+    <vis name="BlueVis"          alpha="1.0" r="0.0" g="0.0"  b="1.0"   showDaughters="true"  visible="true"/>
+    <vis name="CyanVis"          alpha="1.0" r="0.0" g="1.0"  b="1.0"   showDaughters="true"  visible="true"/>
+    <vis name="MagentaVis"       alpha="1.0" r="1.0" g="0.0"  b="1.0"   showDaughters="true"  visible="true"/>
+    <vis name="YellowVis"        alpha="1.0" r="1.0" g="1.0"  b="0.0"   showDaughters="true"  visible="true"/>
+    <vis name="BlackVis"         alpha="1.0" r="0.0" g="0.0"  b="0.0"   showDaughters="true"  visible="true"/>
+    <vis name="GrayVis"          alpha="1.0" r="0.5" g="0.5"  b="0.5"   showDaughters="true"  visible="true"/>
+    <vis name="ECALVis"          alpha="0.5" r="0.0" g="0.0"  b="1.0" showDaughters="true"  visible="true"/>
+    <vis name="HCALVis"       alpha="0.5" r="1.0" g="0.0"  b="0.0" showDaughters="true"  visible="true"/>
+  </display>
+  
+  <detectors>
+    <comment>Calorimeters</comment>
+    <include ref="TBecal4d.xml"/>
+    <include ref="TBhcal4d.xml"/>
+  </detectors>
+
+  <readouts>
+    <readout name="ecalSD">
+      <segmentation type="TiledLayerGridXY" grid_size_x="45*mm" grid_size_y="5*mm" offset_x="-EBU_dim_x/2-4.5/2*cm" offset_y="-EBU_dim_y/2-0.5/2*cm" identifier_x="I" identifier_y="J" identifier_layer="K"/>
+      <id>system:0,K:8,J:8,I:8</id>
+    </readout>
+    <readout name="hcalSD">
+      <segmentation type="TiledLayerGridXY" grid_size_x="3*cm" grid_size_y="3*cm" offset_x="-HBU_BL_dim_x/2-3/2*cm" offset_y="-HBU_BL_dim_y/2-3/2*cm" identifier_x="I" identifier_y="J" identifier_layer="K"/>
+      <id>system:0,K:8,J:8,I:8</id>
+    </readout>
+  </readouts>
+  
+</lccdd>
diff --git a/CaloTB/CaloTB_EPT_AHCAL/compact/TBecal4d.xml b/CaloTB/CaloTB_EPT_AHCAL/compact/TBecal4d.xml
new file mode 100644
index 000000000..957913900
--- /dev/null
+++ b/CaloTB/CaloTB_EPT_AHCAL/compact/TBecal4d.xml
@@ -0,0 +1,30 @@
+<!-- 
+   Calorimeter
+  -->
+
+<detector id="1" name="TBecal4d" type="TBecal4d" vis="GreenVis" readout="ecalSD" insideTrackingVolume="false">
+
+  <comment>AHCAL TB 2015 Model ECAL</comment>
+
+  <envelope vis="ECALVis">
+    <shape type="Box" dx="EBU_dim_x*4/2 + env_safety" dy="EBU_dim_y*4/2 + env_safety"  dz="ECAL_nlayers*Ecal_layer_thickness/2 + env_safety" material="Air" />
+    <rotation x="0" y="0" z="0"/>
+    <position x="0" y="0" z="ECAL_nlayers*Ecal_layer_thickness/2"/>
+  </envelope>
+
+  <type_flags type="1" />
+
+  <layer repeat="ECAL_nlayers" vis="SeeThrough">
+    <slice material = "Steel235"       thickness = "Ecal_radiator_thickness"       vis="BlueVis"                          />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "CFM"            thickness = "cable_mix_thickness"           vis="Invisible"                        />
+    <slice material = "PCB"            thickness = "pcb_thickness"                 vis="GreenVis"                        />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "polystyrene"    thickness = "Ecal_scintillator_thickness"   vis="CyanVis"    sensitive = "yes"     />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+  </layer>
+
+</detector>
diff --git a/CaloTB/CaloTB_EPT_AHCAL/compact/TBhcal4d.xml b/CaloTB/CaloTB_EPT_AHCAL/compact/TBhcal4d.xml
new file mode 100644
index 000000000..b9db06cc8
--- /dev/null
+++ b/CaloTB/CaloTB_EPT_AHCAL/compact/TBhcal4d.xml
@@ -0,0 +1,111 @@
+<!-- 
+   Calorimeter
+  -->
+
+<detector id="2" name="TBhcal4d" type="TBhcal4d" vis="GreenVis" readout="hcalSD" insideTrackingVolume="false">
+
+  <comment>AHCAL TB 2015 Model SSF</comment>
+
+  <envelope vis="HCALVis">
+    <shape type="Box" dx="HBU_BL_dim_x/2 + env_safety" dy="HBU_BL_dim_y/2 + env_safety"  dz="(HCAL_SSF_nlayers*Hcal_layer_thickness + HCAL_BL_nlayers*Hcal_layer_thickness)/2 + env_safety" material="Air" />
+    <rotation x="0" y="0" z="0"/>
+    <position x="0" y="0" z="ECAL_nlayers*Ecal_layer_thickness + (HCAL_SSF_nlayers*Hcal_layer_thickness + HCAL_BL_nlayers*Hcal_layer_thickness)/2"/>
+  </envelope>
+
+  <type_flags type="1" />
+
+  <!-- slot 3 to 10 -->
+  <layer repeat="HCAL_SSF_nlayers" vis="SeeThrough">
+    <slice material = "Steel235"       thickness = "Hcal_radiator_thickness"       vis="BlueVis"                          />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "CFM"            thickness = "cable_mix_thickness"           vis="Invisible"                        />
+    <slice material = "PCB"            thickness = "pcb_thickness"                 vis="GreenVis"                        />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "polystyrene"    thickness = "Hcal_scintillator_thickness"   vis="CyanVis"    sensitive = "yes"     />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+  </layer>
+
+  <!-- slot 11 -->
+  <layer repeat="1" vis="SeeThrough">
+    <slice material = "Steel235"       thickness = "Hcal_radiator_thickness"       vis="BlueVis"                          />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "CFM"            thickness = "cable_mix_thickness"           vis="Invisible"                        />
+    <slice material = "PCB"            thickness = "pcb_thickness"                 vis="GreenVis"                        />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "polystyrene"    thickness = "Hcal_scintillator_thickness"   vis="CyanVis"    sensitive = "yes"     />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+  </layer>
+
+  <!-- slot 12 -->
+  <layer repeat="1" vis="SeeThrough">
+    <slice material = "Steel235"       thickness = "Hcal_radiator_thickness"       vis="BlueVis"                          />
+    <slice material = "Air"            thickness = "airgap_thickness*2+cassette_thickness*2+cable_mix_thickness+pcb_thickness+foil_thickness*2+Hcal_scintillator_thickness"              vis="Invisible"                        />
+  </layer>
+
+  <!-- slot 13 -->
+  <layer repeat="1" vis="SeeThrough">
+    <slice material = "Steel235"       thickness = "Hcal_radiator_thickness"       vis="BlueVis"                          />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "CFM"            thickness = "cable_mix_thickness"           vis="Invisible"                        />
+    <slice material = "PCB"            thickness = "pcb_thickness"                 vis="Invisible"                        />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "polystyrene"    thickness = "Hcal_scintillator_thickness"   vis="CyanVis"    sensitive = "yes"     />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+  </layer>
+
+  <!-- slot 14 to 20 -->
+  <layer repeat="7" vis="SeeThrough">
+    <slice material = "Steel235"       thickness = "Hcal_radiator_thickness"       vis="BlueVis"                          />
+    <slice material = "Air"            thickness = "airgap_thickness*2+cassette_thickness*2+cable_mix_thickness+pcb_thickness+foil_thickness*2+Hcal_scintillator_thickness"              vis="Invisible"                        />
+  </layer>
+
+  <!-- slot 21 -->
+  <layer repeat="1" vis="SeeThrough">
+    <slice material = "Steel235"       thickness = "Hcal_radiator_thickness"       vis="BlueVis"                          />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "CFM"            thickness = "cable_mix_thickness"           vis="Invisible"                        />
+    <slice material = "PCB"            thickness = "pcb_thickness"                 vis="Invisible"                        />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "polystyrene"    thickness = "Hcal_scintillator_thickness"   vis="CyanVis"    sensitive = "yes"     />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+  </layer>
+
+  <!-- slot 22 to 30 -->
+  <layer repeat="9" vis="SeeThrough">
+    <slice material = "Steel235"       thickness = "Hcal_radiator_thickness"       vis="BlueVis"                          />
+    <slice material = "Air"            thickness = "airgap_thickness*2+cassette_thickness*2+cable_mix_thickness+pcb_thickness+foil_thickness*2+Hcal_scintillator_thickness"              vis="Invisible"                        />
+  </layer>
+
+  <!-- slot 31 -->
+  <layer repeat="1" vis="SeeThrough">
+    <slice material = "Steel235"       thickness = "Hcal_radiator_thickness"       vis="BlueVis"                          />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "CFM"            thickness = "cable_mix_thickness"           vis="Invisible"                        />
+    <slice material = "PCB"            thickness = "pcb_thickness"                 vis="Invisible"                        />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "polystyrene"    thickness = "Hcal_scintillator_thickness"   vis="CyanVis"    sensitive = "yes"     />
+    <slice material = "Polystyrole"    thickness = "foil_thickness"                vis="Invisible"                        />
+    <slice material = "Steel235"       thickness = "cassette_thickness"            vis="Invisible"                        />
+    <slice material = "Air"            thickness = "airgap_thickness"              vis="Invisible"                        />
+  </layer>
+
+  <!-- slot 32 -->
+  <layer repeat="1" vis="SeeThrough">
+    <slice material = "Steel235"       thickness = "Hcal_radiator_thickness"       vis="BlueVis"                          />
+    <slice material = "Air"            thickness = "airgap_thickness*2+cassette_thickness*2+cable_mix_thickness+pcb_thickness+foil_thickness*2+Hcal_scintillator_thickness"              vis="Invisible"                        />
+  </layer>
+
+</detector>
diff --git a/CaloTB/CaloTB_EPT_AHCAL/compact/materials_TBModel2015.xml b/CaloTB/CaloTB_EPT_AHCAL/compact/materials_TBModel2015.xml
new file mode 100644
index 000000000..59e84fecc
--- /dev/null
+++ b/CaloTB/CaloTB_EPT_AHCAL/compact/materials_TBModel2015.xml
@@ -0,0 +1,100 @@
+<materials>
+
+  <!--The steel we are going to use in the Hcal: Material S235JR (old name St37)
+	  Numbers found under 
+	  http://n.ethz.ch/student/zwickers/ download/fs_pe_grundlagen_cyrill.pdf  -->
+  <material name="Steel235">
+    <D value="7.87" unit="g/cm3"/>
+    <fraction n="0.9843" ref="Fe"/>
+    <fraction n="0.0017" ref="C"/>
+    <fraction n="0.014" ref="Mn"/>
+  </material>
+
+  <!-- Tungsten is built with density 17.84*g/cm3 and it has 5.25% of nikel, 1.76% of copper and 92.99% of 19.3gr/cm3 tungsten -->
+  <material name="TungstenDens24">
+    <D value="17.84" unit="g/cm3"/>
+    <fraction n="0.9299" ref="W"/>
+    <fraction n="0.0525" ref="Ni"/>
+    <fraction n="0.0176" ref="Fe"/>
+  </material>
+
+  <!-- PCB (Printed Circuit Board) Material FR4
+       Composition and density found under 
+       http://pcba10.ba.infn.it/temp/ddd/ma/materials/list.html 
+    -->
+  <material name="PCB" state="solid">
+    <MEE unit="eV" value="88.255598548367"/>
+    <D unit="g/cm3" value="1.7"/>
+    <fraction n="0.180774" ref="Si"/>
+    <fraction n="0.405633" ref="O"/>
+    <fraction n="0.278042" ref="C"/>
+    <fraction n="0.0684428" ref="H"/>
+    <fraction n="0.0671091" ref="Br"/>
+  </material>
+
+  <!-- Scintillator -->
+  <!--material name="G4_POLYSTYRENE" state="solid">
+  <MEE unit="eV" value="68.7"/>
+  <D unit="g/cm3" value="1.06"/>
+  <fraction n="0.077418" ref="H"/>
+  <fraction n="0.922582" ref="C"/>
+  </material-->
+  <material name="polystyrene" state="solid">
+    <MEE unit="eV" value="68.7"/>
+    <D unit="g/cm3" value="1.032"/>
+    <fraction n="0.5246" ref="H"/>
+    <fraction n="0.4754" ref="C"/>
+  </material>
+
+  <!-- Numbers from http://www.elpac.de/Kunststoffkleinteile/Kleines_Kunststoff-Know-_How/PVC-P/pvc-p.html -->
+  <material name="PVC">
+    <D value="1.35" unit="g/cm3"/>
+    <fraction n="0.5" ref="H"/>
+    <fraction n="0.33" ref="C"/>
+    <fraction n="0.17" ref="Cl"/>
+  </material>
+
+  <!-- 3M foil / Numbers from http://de.wikipedia.org/wiki/Polystyrol
+       the structural formula for the Styrene Polymer is C6H5CH=CH2
+       The difference to Styropor definition in CGAGeometryManager
+       comes since we do not have the material in a foamed form -->
+  <material name="Polystyrole">
+    <D value="1.065" unit="g/cm3"/>
+    <fraction n="0.5" ref="H"/>
+    <fraction n="0.5" ref="C"/>
+  </material>
+
+  <!-- -->
+  <material name="Air">
+    <D value="0.00129" unit="g/cm3"/>
+    <fraction n="0.3" ref="O"/>
+    <fraction n="0.7" ref="N"/>
+  </material>
+
+  <!-- Now we define the material cf_mix
+       We assume the following:
+       a) a layer has a volume of V_total = 90x90x0.15 cm^3 = 1215 cm^3(last number is
+       longitudinal space reserved for cable fibre mix)
+       b) coax cable has diameter of 0.12 cm 
+       fibre has diameter of 0.5 cm
+       The cables are on average 45 cm long 
+       => V_coax = PI*(0.06 cm)^2*45 cm = 0.510 cm^3
+  => V_fibre = PI*(0.025 cm)^2*45 cm = 0.088 cm^3
+  ...The rest is occupied by air 
+  V_air = V_total - V_coax - V_fibre = 1214.402
+  There is one coax. cable and one fibre per tile
+  and we have on average 185 tiles per layer
+  => Total mass of coax cable (fibre), m = density*V
+  m_coax = (1.35*0.510)*185 = 127.37 g
+  m_fibre = (1.065*0.088)*185 = 17.33 g
+  ... and
+  m_air = 1214.402*1.29e-03 = 1.45 g
+  total density = (m_air + m_coax + m_fibre)/1215. = 0.120 g/cm^3 -->
+  <material name="CFM">
+    <D value="0.12" unit="g/cm3"/>
+    <fraction n="0.009" ref="Air"/>
+    <fraction n="0.872" ref="PVC"/>
+    <fraction n="0.119" ref="Polystyrole"/>
+  </material>
+  
+</materials>