From 27ffe6e238ea1583a26439b83e05f4ea3945dac5 Mon Sep 17 00:00:00 2001
From: Markus Frank <Markus.Frank@cern.ch>
Date: Wed, 10 Aug 2022 11:09:22 +0200
Subject: [PATCH] Add channeling example with Geant4 specialization for
 G4ExtendedMaterial and G4LogicalCrystalVolume

---
 examples/DDG4/compact/Channeling.xml          |   87 +
 examples/DDG4/data/Si220pl/Si220pl_atd.txt    | 2050 +++++++++++++++++
 examples/DDG4/data/Si220pl/Si220pl_efx.txt    | 2050 +++++++++++++++++
 examples/DDG4/data/Si220pl/Si220pl_efy.txt    | 2050 +++++++++++++++++
 examples/DDG4/data/Si220pl/Si220pl_eld.txt    | 2050 +++++++++++++++++
 examples/DDG4/data/Si220pl/Si220pl_pot.txt    | 2050 +++++++++++++++++
 examples/DDG4/scripts/Channeling.py           |   81 +
 .../DDG4/src/ChannelingCrystalMaterial.cpp    |   81 +
 examples/DDG4/src/ChannelingCrystalVolume.cpp |   52 +
 9 files changed, 10551 insertions(+)
 create mode 100644 examples/DDG4/compact/Channeling.xml
 create mode 100644 examples/DDG4/data/Si220pl/Si220pl_atd.txt
 create mode 100644 examples/DDG4/data/Si220pl/Si220pl_efx.txt
 create mode 100644 examples/DDG4/data/Si220pl/Si220pl_efy.txt
 create mode 100644 examples/DDG4/data/Si220pl/Si220pl_eld.txt
 create mode 100644 examples/DDG4/data/Si220pl/Si220pl_pot.txt
 create mode 100644 examples/DDG4/scripts/Channeling.py
 create mode 100644 examples/DDG4/src/ChannelingCrystalMaterial.cpp
 create mode 100644 examples/DDG4/src/ChannelingCrystalVolume.cpp

diff --git a/examples/DDG4/compact/Channeling.xml b/examples/DDG4/compact/Channeling.xml
new file mode 100644
index 000000000..b3ece44da
--- /dev/null
+++ b/examples/DDG4/compact/Channeling.xml
@@ -0,0 +1,87 @@
+<lccdd>
+<!-- #==========================================================================
+     #  AIDA Detector description implementation 
+     #==========================================================================
+     # Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
+     # All rights reserved.
+     #
+     # For the licensing terms see $DD4hepINSTALL/LICENSE.
+     # For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
+     #
+     #==========================================================================
+-->
+
+  <info name="TestMaterialProperties"
+    title="Test reading of TGeo's Material Properties"
+    author="Markus Frank"
+    url="None"
+    status="development"
+    version="1.0">
+    <comment>Test reading of TGeo's Material Properties</comment>
+  </info>
+
+  <define>
+    <constant name="world_side" value="1*m"/>
+    <constant name="world_x"    value="world_side/2"/>
+    <constant name="world_y"    value="world_side/2"/>
+    <constant name="world_z"    value="world_side/2"/>
+  </define>
+  <debug>
+    <type name="surface" value="1"/>
+  </debug>
+
+  <includes>
+    <gdmlFile  ref="${DD4hepINSTALL}/DDDetectors/compact/elements.xml"/>
+  </includes>
+
+  <materials>
+    <material name="Air">
+      <D type="density" unit="g/cm3" value="0.0012"/>
+      <fraction n="0.754" ref="N"/>
+      <fraction n="0.234" ref="O"/>
+      <fraction n="0.012" ref="Ar"/>
+    </material>
+    <!-- We model vakuum just as very thin air -->
+    <material name="Vacuum">  
+      <D type="density" unit="g/cm3" value="0.0000000001" />
+      <fraction n="0.754" ref="N"/>
+      <fraction n="0.234" ref="O"/>
+      <fraction n="0.012" ref="Ar"/>
+    </material>
+    <material name="Ice">
+      <D type="density" value="1.0" unit="g/cm3"/>
+      <composite n="2" ref="H"/>
+      <composite n="1" ref="O"/>
+      <constant name="Geant4-plugin:material"            option="ChannelingCrystalMaterial"/>
+      <constant name="Geant4-ignore:crystal_data"        option="${DD4hepExamplesINSTALL}/examples/DDG4/data/Si220pl/Si220pl"/>
+      <constant name="Geant4-ignore:crystal_orientation" option="(1,0,0)"/>
+    </material>
+  </materials>
+
+  <display>
+    <vis name="EnvelopeVis"     alpha="0.1"   r="0"   g="0.0" b="1.0"  showDaughters="true"  visible="true"/>
+    <vis name="CrystalVis"  alpha="0.3"   r="1"   g="0.0" b="0.0"  showDaughters="true"  visible="true"/>
+  </display>
+
+  <detectors>
+    <detector id ="1" name="ChannelingDevice" type="DD4hep_SingleShape" vis="CrystalVis" sensitive="true" readout="ChannelingDeviceHits">
+      <box   x="1.3*m" y="1.3*m" z="1.3*m" vis="EnvelopeVis"/>
+      <shape type="Box" dx="50*cm" dy="50*cm" dz="20*cm">
+	<position x="0" y="0" z="-100*cm"/>
+	<position x="0" y="0" z="-50*cm"/>
+	<position x="0" y="0" z="0*cm"/>
+	<position x="0" y="0" z="50*cm"/>
+	<position x="0" y="0" z="100*cm"/>
+	<property name="Geant4-plugin" value="ChannelingCrystalVolume"/>
+      </shape>
+      <position x="0" y="0" z="0"/>
+      <material name="Ice"/>
+    </detector>
+  </detectors>
+
+  <readouts>
+    <readout name="ChannelingDeviceHits">
+      <id>system:8,module:8,x:32:-16,y:-16</id>
+    </readout>
+  </readouts>
+</lccdd>
diff --git a/examples/DDG4/data/Si220pl/Si220pl_atd.txt b/examples/DDG4/data/Si220pl/Si220pl_atd.txt
new file mode 100644
index 000000000..e1ade5a87
--- /dev/null
+++ b/examples/DDG4/data/Si220pl/Si220pl_atd.txt
@@ -0,0 +1,2050 @@
+2048 1 1
+1.91979e-10 2.21679e-10 3.13501e-10
+12.1802
+12.1788
+12.1748
+12.168
+12.1585
+12.1463
+12.1315
+12.114
+12.0938
+12.0709
+12.0454
+12.0173
+11.9866
+11.9534
+11.9175
+11.8791
+11.8382
+11.7949
+11.749
+11.7008
+11.6501
+11.5971
+11.5418
+11.4842
+11.4243
+11.3622
+11.2979
+11.2315
+11.163
+11.0925
+11.0199
+10.9454
+10.869
+10.7907
+10.7105
+10.6286
+10.545
+10.4598
+10.3729
+10.2844
+10.1944
+10.103
+10.0102
+9.91596
+9.82047
+9.72373
+9.6258
+9.52674
+9.42661
+9.32545
+9.22332
+9.12028
+9.01639
+8.9117
+8.80626
+8.70014
+8.59339
+8.48605
+8.37819
+8.26987
+8.16113
+8.05202
+7.94261
+7.83294
+7.72307
+7.61305
+7.50292
+7.39274
+7.28256
+7.17243
+7.06239
+6.95249
+6.84278
+6.7333
+6.6241
+6.51523
+6.40671
+6.2986
+6.19094
+6.08377
+5.97712
+5.87103
+5.76555
+5.6607
+5.55652
+5.45305
+5.35031
+5.24834
+5.14716
+5.04682
+4.94733
+4.84872
+4.75103
+4.65426
+4.55845
+4.46362
+4.36979
+4.27699
+4.18522
+4.09451
+4.00487
+3.91633
+3.82889
+3.74257
+3.65738
+3.57334
+3.49045
+3.40873
+3.32818
+3.24881
+3.17063
+3.09364
+3.01785
+2.94326
+2.86987
+2.79769
+2.72673
+2.65697
+2.58842
+2.52107
+2.45494
+2.39
+2.32627
+2.26373
+2.20238
+2.14222
+2.08324
+2.02543
+1.96879
+1.91331
+1.85898
+1.80578
+1.75372
+1.70279
+1.65296
+1.60423
+1.5566
+1.51004
+1.46455
+1.42012
+1.37672
+1.33436
+1.29301
+1.25266
+1.21331
+1.17492
+1.1375
+1.10103
+1.06549
+1.03086
+0.997144
+0.964312
+0.932353
+0.901254
+0.870998
+0.84157
+0.812956
+0.78514
+0.758107
+0.731842
+0.70633
+0.681556
+0.657504
+0.63416
+0.611508
+0.589535
+0.568225
+0.547563
+0.527535
+0.508127
+0.489324
+0.471112
+0.453476
+0.436404
+0.419881
+0.403894
+0.388429
+0.373474
+0.359014
+0.345037
+0.331531
+0.318482
+0.305879
+0.293709
+0.281961
+0.270623
+0.259682
+0.249129
+0.238951
+0.229138
+0.21968
+0.210564
+0.201782
+0.193324
+0.185179
+0.177337
+0.16979
+0.162528
+0.155542
+0.148823
+0.142362
+0.136152
+0.130184
+0.124449
+0.118941
+0.113651
+0.108572
+0.103698
+0.0990198
+0.0945319
+0.0902273
+0.0860995
+0.0821424
+0.0783496
+0.0747154
+0.0712339
+0.0678995
+0.0647068
+0.0616505
+0.0587255
+0.0559268
+0.0532497
+0.0506894
+0.0482415
+0.0459016
+0.0436655
+0.041529
+0.0394884
+0.0375396
+0.0356791
+0.0339032
+0.0322086
+0.0305918
+0.0290498
+0.0275793
+0.0261775
+0.0248414
+0.0235682
+0.0223553
+0.0212002
+0.0201002
+0.0190531
+0.0180565
+0.0171082
+0.0162061
+0.0153482
+0.0145324
+0.013757
+0.01302
+0.0123198
+0.0116546
+0.0110229
+0.0104231
+0.00985381
+0.0093135
+0.00880085
+0.00831458
+0.00785342
+0.00741619
+0.00700175
+0.006609
+0.00623689
+0.00588442
+0.00555064
+0.00523462
+0.0049355
+0.00465244
+0.00438463
+0.00413132
+0.00389178
+0.00366531
+0.00345126
+0.00324898
+0.00305788
+0.00287737
+0.00270692
+0.00254601
+0.00239412
+0.00225079
+0.00211558
+0.00198804
+0.00186778
+0.0017544
+0.00164754
+0.00154684
+0.00145198
+0.00136263
+0.00127849
+0.00119928
+0.00112473
+0.00105458
+0.000988582
+0.00092651
+0.000868142
+0.00081327
+0.000761697
+0.000713236
+0.000667709
+0.00062495
+0.000584798
+0.000547105
+0.000511727
+0.00047853
+0.000447388
+0.000418179
+0.00039079
+0.000365114
+0.000341049
+0.000318499
+0.000297374
+0.000277588
+0.000259061
+0.000241717
+0.000225484
+0.000210295
+0.000196085
+0.000182794
+0.000170366
+0.000158748
+0.00014789
+0.000137743
+0.000128264
+0.000119411
+0.000111144
+0.000103427
+9.62237e-05
+8.95023e-05
+8.3232e-05
+7.73837e-05
+7.19303e-05
+6.68464e-05
+6.2108e-05
+5.76926e-05
+5.35792e-05
+4.97481e-05
+4.61805e-05
+4.28593e-05
+3.97681e-05
+3.68917e-05
+3.42157e-05
+3.17267e-05
+2.94122e-05
+2.72606e-05
+2.52607e-05
+2.34023e-05
+2.16758e-05
+2.00722e-05
+1.85831e-05
+1.72007e-05
+1.59176e-05
+1.47269e-05
+1.36222e-05
+1.25976e-05
+1.16475e-05
+1.07666e-05
+9.95015e-06
+9.19356e-06
+8.49261e-06
+7.84336e-06
+7.24213e-06
+6.6855e-06
+6.17028e-06
+5.6935e-06
+5.25239e-06
+4.84438e-06
+4.46707e-06
+4.11823e-06
+3.79579e-06
+3.49782e-06
+3.22252e-06
+2.96822e-06
+2.73339e-06
+2.51658e-06
+2.31644e-06
+2.13175e-06
+1.96135e-06
+1.80417e-06
+1.65921e-06
+1.52557e-06
+1.40237e-06
+1.28884e-06
+1.18423e-06
+1.08788e-06
+9.99136e-07
+9.17431e-07
+8.42221e-07
+7.73004e-07
+7.09318e-07
+6.50734e-07
+5.96855e-07
+5.47316e-07
+5.01777e-07
+4.59925e-07
+4.2147e-07
+3.86144e-07
+3.537e-07
+3.2391e-07
+2.96564e-07
+2.71465e-07
+2.48436e-07
+2.27309e-07
+2.07933e-07
+1.90166e-07
+1.73879e-07
+1.58951e-07
+1.45273e-07
+1.32742e-07
+1.21265e-07
+1.10755e-07
+1.01134e-07
+9.23283e-08
+8.42704e-08
+7.68987e-08
+7.01562e-08
+6.39907e-08
+5.8354e-08
+5.3202e-08
+4.84941e-08
+4.41929e-08
+4.02643e-08
+3.66768e-08
+3.34014e-08
+3.04119e-08
+2.76837e-08
+2.51946e-08
+2.29243e-08
+2.08539e-08
+1.89663e-08
+1.72457e-08
+1.56777e-08
+1.42491e-08
+1.29477e-08
+1.17627e-08
+1.06837e-08
+9.7015e-09
+8.80766e-09
+7.99439e-09
+7.25461e-09
+6.58182e-09
+5.97009e-09
+5.41402e-09
+4.90865e-09
+4.44946e-09
+4.03233e-09
+3.65349e-09
+3.30951e-09
+2.99725e-09
+2.71384e-09
+2.45669e-09
+2.22341e-09
+2.01183e-09
+1.81999e-09
+1.64607e-09
+1.48844e-09
+1.3456e-09
+1.2162e-09
+1.099e-09
+9.92876e-10
+8.96799e-10
+8.09837e-10
+7.31147e-10
+6.59954e-10
+5.95561e-10
+5.37334e-10
+4.84691e-10
+4.37106e-10
+3.94107e-10
+3.55258e-10
+3.20167e-10
+2.8848e-10
+2.59869e-10
+2.34044e-10
+2.1074e-10
+1.89713e-10
+1.70745e-10
+1.53642e-10
+1.38219e-10
+1.24318e-10
+1.1179e-10
+1.00501e-10
+9.03321e-11
+8.11768e-11
+7.293e-11
+6.55067e-11
+5.88287e-11
+5.28168e-11
+4.74101e-11
+4.25481e-11
+3.81739e-11
+3.42432e-11
+3.0711e-11
+2.7534e-11
+2.46839e-11
+2.21222e-11
+1.98225e-11
+1.77564e-11
+1.59042e-11
+1.42388e-11
+1.27497e-11
+1.14107e-11
+1.02084e-11
+9.13514e-12
+8.17035e-12
+7.30604e-12
+6.53344e-12
+5.83977e-12
+5.21594e-12
+4.6626e-12
+4.165e-12
+3.71725e-12
+3.32212e-12
+2.96507e-12
+2.64366e-12
+2.36189e-12
+2.10665e-12
+1.87772e-12
+1.67588e-12
+1.49281e-12
+1.32949e-12
+1.18727e-12
+1.05582e-12
+9.40026e-13
+8.37885e-13
+7.45293e-13
+6.6358e-13
+5.91305e-13
+5.2558e-13
+4.67737e-13
+4.16778e-13
+3.68039e-13
+3.29625e-13
+2.92211e-13
+2.57905e-13
+2.31148e-13
+2.02727e-13
+1.80522e-13
+1.62315e-13
+1.42442e-13
+1.26343e-13
+1.15241e-13
+9.89209e-14
+8.90399e-14
+8.08242e-14
+6.90559e-14
+6.23945e-14
+5.68434e-14
+4.89608e-14
+4.30767e-14
+4.00791e-14
+3.41949e-14
+3.03091e-14
+2.67564e-14
+2.42029e-14
+2.13163e-14
+1.94289e-14
+1.70974e-14
+1.4766e-14
+1.35447e-14
+1.11022e-14
+1.07692e-14
+9.32587e-15
+7.10543e-15
+7.88258e-15
+6.55032e-15
+4.55191e-15
+5.77316e-15
+5.21805e-15
+3.44169e-15
+3.77476e-15
+4.10783e-15
+2.9976e-15
+3.55271e-15
+3.44169e-15
+2.77556e-15
+2.66454e-15
+2.66454e-15
+1.88738e-15
+3.66374e-15
+3.33067e-15
+1.44329e-15
+3.44169e-15
+2.44249e-15
+1.77636e-15
+3.55271e-15
+3.10862e-15
+1.55431e-15
+3.44169e-15
+2.66454e-15
+1.66533e-15
+-4.44089e-16
+-4.44089e-16
+1.33227e-15
+9.99201e-16
+1.66533e-15
+1.77636e-15
+-2.22045e-16
+-1.33227e-15
+1.77636e-15
+-2.22045e-16
+-4.44089e-16
+1.66533e-15
+7.77156e-16
+1.11022e-16
+1.11022e-15
+1.66533e-15
+-1.33227e-15
+2.22045e-15
+1.44329e-15
+-2.22045e-16
+1.55431e-15
+3.33067e-16
+-6.66134e-16
+1.77636e-15
+8.88178e-16
+0
+1.77636e-15
+7.77156e-16
+0
+1.11022e-16
+5.55112e-16
+-6.66134e-16
+1.11022e-15
+-1.11022e-15
+1.11022e-16
+-8.88178e-16
+6.66134e-16
+0
+-4.44089e-16
+-2.22045e-16
+1.11022e-16
+-1.33227e-15
+-8.88178e-16
+1.66533e-15
+-1.55431e-15
+-1.11022e-15
+1.11022e-16
+-1.55431e-15
+-2.44249e-15
+-1.55431e-15
+-6.66134e-16
+-1.9984e-15
+-1.11022e-15
+-2.66454e-15
+-1.11022e-15
+-6.66134e-16
+-1.33227e-15
+-1.55431e-15
+-6.66134e-16
+-3.10862e-15
+-1.11022e-15
+-2.88658e-15
+-2.22045e-16
+-1.9984e-15
+1.11022e-16
+-1.77636e-15
+-1.9984e-15
+-8.88178e-16
+-1.9984e-15
+-2.88658e-15
+-1.55431e-15
+-1.9984e-15
+-1.11022e-15
+-1.55431e-15
+-1.77636e-15
+-1.77636e-15
+-8.88178e-16
+-2.22045e-15
+-2.22045e-15
+-1.77636e-15
+-1.9984e-15
+-2.66454e-15
+-1.55431e-15
+-2.44249e-15
+-3.33067e-15
+-2.22045e-15
+-2.44249e-15
+-3.33067e-15
+-2.22045e-15
+-2.88658e-15
+-3.10862e-15
+-1.55431e-15
+-2.66454e-15
+-3.33067e-15
+-1.11022e-15
+-3.10862e-15
+-3.55271e-15
+-2.22045e-15
+-2.66454e-15
+-3.55271e-15
+-1.33227e-15
+-2.88658e-15
+-3.10862e-15
+-1.11022e-15
+-1.9984e-15
+-3.10862e-15
+-8.88178e-16
+-2.88658e-15
+-3.10862e-15
+-6.66134e-16
+-8.88178e-16
+-2.88658e-15
+-1.9984e-15
+-3.33067e-15
+-3.10862e-15
+-2.44249e-15
+-3.77476e-15
+-2.66454e-15
+-8.88178e-16
+-2.66454e-15
+-3.33067e-15
+-4.44089e-16
+-4.21885e-15
+-2.88658e-15
+-1.33227e-15
+-3.10862e-15
+-4.44089e-16
+-1.33227e-15
+-3.55271e-15
+9.99201e-16
+-1.33227e-15
+-3.55271e-15
+1.11022e-15
+-8.88178e-16
+-3.77476e-15
+-2.44249e-15
+-3.10862e-15
+-8.88178e-16
+9.99201e-16
+-2.88658e-15
+-1.11022e-15
+7.77156e-16
+-6.66134e-16
+-2.22045e-16
+1.33227e-15
+-6.66134e-16
+0
+1.22125e-15
+-1.9984e-15
+-2.22045e-16
+-6.66134e-16
+-1.9984e-15
+1.44329e-15
+-2.22045e-16
+-2.66454e-15
+3.33067e-16
+-1.9984e-15
+-1.55431e-15
+1.11022e-15
+2.55351e-15
+-8.88178e-16
+1.11022e-15
+-1.11022e-15
+-8.88178e-16
+1.55431e-15
+1.11022e-15
+-1.55431e-15
+1.22125e-15
+-1.33227e-15
+-8.88178e-16
+1.44329e-15
+5.55112e-16
+-4.44089e-16
+9.99201e-16
+-1.9984e-15
+-1.11022e-15
+1.88738e-15
+-2.44249e-15
+-2.22045e-16
+1.88738e-15
+-2.22045e-15
+-4.44089e-16
+1.9984e-15
+1.33227e-15
+5.55112e-16
+1.22125e-15
+-2.44249e-15
+2.9976e-15
+6.66134e-16
+-2.22045e-15
+4.44089e-16
+7.77156e-16
+-2.22045e-15
+1.11022e-16
+7.77156e-16
+-1.9984e-15
+5.55112e-16
+-4.44089e-16
+6.66134e-16
+3.21965e-15
+9.99201e-16
+1.44329e-15
+3.44169e-15
+7.77156e-16
+7.77156e-16
+3.10862e-15
+4.44089e-16
+-1.55431e-15
+2.77556e-15
+-4.44089e-16
+6.66134e-16
+1.11022e-16
+0
+4.44089e-16
+3.10862e-15
+0
+1.11022e-16
+3.10862e-15
+0
+5.55112e-16
+3.77476e-15
+9.99201e-16
+1.22125e-15
+-2.22045e-16
+0
+5.55112e-16
+2.10942e-15
+-4.44089e-16
+1.44329e-15
+-4.44089e-16
+2.22045e-16
+8.88178e-16
+8.88178e-16
+9.99201e-16
+1.9984e-15
+-1.33227e-15
+1.33227e-15
+7.77156e-16
+1.88738e-15
+1.9984e-15
+7.77156e-16
+-2.22045e-16
+1.55431e-15
+1.33227e-15
+1.77636e-15
+1.44329e-15
+1.22125e-15
+0
+1.11022e-15
+1.11022e-16
+-2.22045e-16
+1.44329e-15
+1.33227e-15
+-4.44089e-16
+1.22125e-15
+5.55112e-16
+1.11022e-16
+1.77636e-15
+9.99201e-16
+0
+1.33227e-15
+3.33067e-16
+0
+5.55112e-16
+4.44089e-16
+2.22045e-16
+1.22125e-15
+9.99201e-16
+5.55112e-16
+1.22125e-15
+3.33067e-16
+-4.44089e-16
+-2.22045e-16
+6.66134e-16
+-2.22045e-16
+2.22045e-16
+5.55112e-16
+0
+-2.22045e-16
+5.55112e-16
+6.66134e-16
+2.22045e-16
+0
+4.44089e-16
+1.11022e-16
+-2.22045e-16
+6.66134e-16
+1.11022e-16
+-2.22045e-16
+8.88178e-16
+1.11022e-15
+1.11022e-16
+4.44089e-16
+1.11022e-16
+8.88178e-16
+5.55112e-16
+5.55112e-16
+3.33067e-16
+8.88178e-16
+4.44089e-16
+1.44329e-15
+3.33067e-16
+1.22125e-15
+3.33067e-16
+5.55112e-16
+-2.22045e-16
+1.11022e-15
+8.88178e-16
+1.77636e-15
+1.11022e-15
+4.44089e-16
+2.10942e-15
+1.22125e-15
+5.55112e-16
+1.9984e-15
+1.11022e-15
+-2.22045e-16
+1.33227e-15
+6.66134e-16
+1.11022e-16
+1.9984e-15
+3.33067e-16
+9.99201e-16
+2.44249e-15
+-2.22045e-16
+6.66134e-16
+2.10942e-15
+8.88178e-16
+5.55112e-16
+2.10942e-15
+-4.44089e-16
+1.33227e-15
+-2.22045e-16
+-2.22045e-16
+0
+6.66134e-16
+-4.44089e-16
+8.88178e-16
+-4.44089e-16
+-8.88178e-16
+3.33067e-16
+-2.22045e-16
+-1.55431e-15
+0
+-4.44089e-16
+1.88738e-15
+1.11022e-15
+4.44089e-16
+3.33067e-16
+1.66533e-15
+0
+-6.66134e-16
+1.55431e-15
+-4.44089e-16
+0
+2.10942e-15
+1.11022e-16
+2.66454e-15
+2.33147e-15
+5.55112e-16
+8.88178e-16
+1.44329e-15
+1.22125e-15
+2.22045e-15
+2.55351e-15
+5.55112e-16
+2.22045e-16
+2.55351e-15
+5.55112e-16
+1.66533e-15
+2.44249e-15
+-4.44089e-16
+3.33067e-15
+2.9976e-15
+1.22125e-15
+1.77636e-15
+2.9976e-15
+2.22045e-16
+-4.44089e-16
+2.66454e-15
+1.22125e-15
+1.33227e-15
+-1.11022e-15
+1.44329e-15
+2.44249e-15
+2.33147e-15
+8.88178e-16
+1.66533e-15
+1.77636e-15
+6.66134e-16
+2.33147e-15
+2.22045e-16
+1.22125e-15
+-8.88178e-16
+-2.22045e-16
+1.22125e-15
+4.44089e-16
+9.99201e-16
+1.66533e-15
+1.77636e-15
+8.88178e-16
+1.66533e-15
+1.88738e-15
+-4.44089e-16
+7.77156e-16
+1.11022e-16
+1.77636e-15
+1.33227e-15
+6.66134e-16
+3.33067e-16
+2.22045e-15
+7.77156e-16
+-4.44089e-16
+1.22125e-15
+3.33067e-16
+5.55112e-16
+9.99201e-16
+9.99201e-16
+8.88178e-16
+1.22125e-15
+4.44089e-16
+6.66134e-16
+7.77156e-16
+-2.22045e-16
+4.44089e-16
+6.66134e-16
+5.55112e-16
+3.33067e-16
+1.11022e-15
+4.44089e-16
+0
+0
+9.99201e-16
+5.55112e-16
+3.33067e-16
+3.33067e-16
+7.77156e-16
+1.22125e-15
+5.55112e-16
+8.88178e-16
+-2.22045e-16
+-4.44089e-16
+1.11022e-16
+0
+4.44089e-16
+-4.44089e-16
+-4.44089e-16
+-6.66134e-16
+0
+-6.66134e-16
+3.33067e-16
+-2.22045e-16
+3.33067e-16
+-6.66134e-16
+5.55112e-16
+5.55112e-16
+-6.66134e-16
+-8.88178e-16
+3.33067e-16
+0
+0
+8.88178e-16
+1.11022e-16
+-1.11022e-15
+-2.22045e-16
+-4.44089e-16
+4.44089e-16
+3.33067e-16
+-1.77636e-15
+1.11022e-16
+-6.66134e-16
+3.33067e-16
+3.33067e-16
+1.11022e-15
+3.33067e-16
+4.44089e-16
+-2.22045e-16
+2.22045e-16
+-8.88178e-16
+-6.66134e-16
+-1.33227e-15
+1.11022e-16
+-6.66134e-16
+-8.88178e-16
+5.55112e-16
+-4.44089e-16
+1.11022e-16
+-2.22045e-16
+-2.22045e-16
+-1.33227e-15
+-1.11022e-15
+-6.66134e-16
+-2.44249e-15
+-8.88178e-16
+-1.33227e-15
+-2.44249e-15
+-1.33227e-15
+3.33067e-16
+-6.66134e-16
+-1.11022e-15
+-8.88178e-16
+-1.77636e-15
+-4.44089e-16
+-1.55431e-15
+-1.33227e-15
+-1.11022e-15
+-2.22045e-15
+-2.22045e-15
+-2.66454e-15
+-1.33227e-15
+-1.77636e-15
+-1.55431e-15
+-8.88178e-16
+5.55112e-16
+-2.22045e-16
+-1.11022e-15
+1.44329e-15
+2.22045e-16
+-2.22045e-15
+1.11022e-15
+-4.44089e-16
+-2.22045e-15
+-1.11022e-15
+-8.88178e-16
+-1.11022e-15
+1.55431e-15
+6.66134e-16
+8.88178e-16
+1.55431e-15
+7.77156e-16
+1.44329e-15
+1.55431e-15
+4.44089e-16
+3.33067e-16
+1.55431e-15
+0
+-8.88178e-16
+1.33227e-15
+1.33227e-15
+4.44089e-16
+6.66134e-16
+1.9984e-15
+8.88178e-16
+-1.9984e-15
+1.33227e-15
+-2.22045e-16
+1.11022e-15
+1.66533e-15
+5.55112e-16
+1.11022e-15
+-1.11022e-15
+-2.22045e-15
+1.9984e-15
+2.22045e-15
+-2.22045e-15
+-1.77636e-15
+2.10942e-15
+-2.22045e-15
+1.88738e-15
+2.22045e-15
+3.33067e-16
+-1.77636e-15
+1.11022e-16
+-1.77636e-15
+-1.33227e-15
+-2.22045e-16
+-2.88658e-15
+2.22045e-16
+1.11022e-16
+-1.9984e-15
+2.44249e-15
+-8.88178e-16
+-2.66454e-15
+-1.77636e-15
+1.11022e-16
+2.22045e-16
+-2.22045e-16
+-1.33227e-15
+-2.22045e-16
+1.11022e-16
+9.99201e-16
+1.11022e-16
+1.33227e-15
+-6.66134e-16
+1.44329e-15
+3.33067e-16
+0
+6.66134e-16
+-6.66134e-16
+-4.44089e-16
+-1.11022e-15
+6.66134e-16
+-8.88178e-16
+-2.22045e-16
+1.33227e-15
+4.44089e-16
+6.66134e-16
+-6.66134e-16
+-1.33227e-15
+1.11022e-16
+-8.88178e-16
+-1.55431e-15
+1.66533e-15
+5.55112e-16
+-1.77636e-15
+0
+-1.33227e-15
+-1.55431e-15
+9.99201e-16
+4.44089e-16
+-2.22045e-15
+1.11022e-16
+-1.11022e-15
+1.22125e-15
+7.77156e-16
+-2.44249e-15
+-2.44249e-15
+4.44089e-16
+-1.55431e-15
+1.44329e-15
+3.33067e-16
+-2.22045e-16
+-1.33227e-15
+3.33067e-16
+-6.66134e-16
+6.66134e-16
+-4.44089e-16
+-1.55431e-15
+-1.9984e-15
+9.99201e-16
+-2.22045e-16
+9.99201e-16
+-6.66134e-16
+-1.77636e-15
+-1.9984e-15
+1.11022e-15
+-6.66134e-16
+2.66454e-15
+2.22045e-16
+8.88178e-16
+7.77156e-16
+9.99201e-16
+-6.66134e-16
+2.77556e-15
+5.55112e-16
+-6.66134e-16
+8.88178e-16
+-3.9968e-15
+-8.88178e-16
+1.9984e-15
+-2.44249e-15
+0
+1.33227e-15
+7.77156e-16
+-2.22045e-16
+-1.9984e-15
+1.33227e-15
+-6.66134e-16
+-3.33067e-15
+-3.9968e-15
+-2.22045e-16
+1.77636e-15
+-2.66454e-15
+-1.11022e-15
+1.9984e-15
+-2.88658e-15
+-4.44089e-16
+1.66533e-15
+-2.66454e-15
+-1.11022e-15
+1.22125e-15
+-2.66454e-15
+1.77636e-15
+-3.55271e-15
+-1.55431e-15
+2.9976e-15
+-4.21885e-15
+7.77156e-16
+1.11022e-16
+-3.10862e-15
+-1.55431e-15
+-4.44089e-16
+-5.55112e-15
+0
+-8.88178e-16
+-8.88178e-16
+1.55431e-15
+-8.88178e-16
+4.44089e-16
+1.88738e-15
+-2.22045e-16
+-8.88178e-16
+1.88738e-15
+-8.88178e-16
+7.77156e-16
+1.88738e-15
+2.9976e-15
+8.88178e-16
+-1.55431e-15
+4.10783e-15
+7.77156e-16
+1.66533e-15
+2.9976e-15
+3.33067e-16
+-1.55431e-15
+4.88498e-15
+1.11022e-16
+2.33147e-15
+-8.88178e-16
+1.11022e-16
+1.66533e-15
+-2.44249e-15
+9.99201e-16
+3.21965e-15
+-4.44089e-16
+9.99201e-16
+-2.88658e-15
+-6.66134e-16
+-2.22045e-16
+-3.33067e-15
+2.44249e-15
+7.77156e-16
+-2.44249e-15
+-4.44089e-16
+2.22045e-15
+-3.33067e-15
+4.44089e-15
+1.33227e-15
+-2.66454e-15
+5.10703e-15
+0
+-2.88658e-15
+4.44089e-15
+1.55431e-15
+-1.55431e-15
+5.32907e-15
+2.33147e-15
+-3.10862e-15
+2.66454e-15
+3.9968e-15
+1.22125e-15
+-8.88178e-16
+4.55191e-15
+-2.22045e-15
+2.9976e-15
+4.66294e-15
+-1.11022e-15
+-4.44089e-16
+4.996e-15
+-1.33227e-15
+3.88578e-15
+4.66294e-15
+-1.33227e-15
+0
+4.32987e-15
+-2.22045e-16
+3.21965e-15
+4.77396e-15
+1.11022e-16
+3.33067e-16
+2.88658e-15
+-2.22045e-16
+3.10862e-15
+4.55191e-15
+1.33227e-15
+4.44089e-16
+3.66374e-15
+3.33067e-15
+6.66134e-16
+4.55191e-15
+2.33147e-15
+8.88178e-16
+2.44249e-15
+1.44329e-15
+3.88578e-15
+4.21885e-15
+2.66454e-15
+1.11022e-16
+2.44249e-15
+2.44249e-15
+4.10783e-15
+3.21965e-15
+1.9984e-15
+4.44089e-16
+2.77556e-15
+1.77636e-15
+3.88578e-15
+4.77396e-15
+3.33067e-15
+7.77156e-16
+-6.66134e-16
+1.77636e-15
+2.10942e-15
+3.66374e-15
+1.77636e-15
+1.33227e-15
+0
+2.88658e-15
+2.9976e-15
+4.21885e-15
+9.99201e-16
+8.88178e-16
+-4.44089e-16
+2.22045e-16
+2.44249e-15
+2.88658e-15
+8.88178e-16
+1.88738e-15
+3.33067e-16
+2.33147e-15
+2.77556e-15
+3.21965e-15
+4.44089e-16
+9.99201e-16
+-1.77636e-15
+-1.11022e-15
+1.66533e-15
+2.66454e-15
+2.22045e-16
+4.44089e-16
+-1.9984e-15
+-2.66454e-15
+-4.44089e-16
+2.10942e-15
+-4.44089e-16
+2.44249e-15
+-1.77636e-15
+-1.9984e-15
+2.77556e-15
+2.10942e-15
+-1.33227e-15
+1.55431e-15
+-2.44249e-15
+-2.66454e-15
+1.77636e-15
+8.88178e-16
+1.9984e-15
+1.11022e-16
+1.88738e-15
+-3.10862e-15
+1.22125e-15
+2.44249e-15
+-2.66454e-15
+-6.66134e-16
+2.33147e-15
+-3.33067e-15
+-2.66454e-15
+1.77636e-15
+-3.33067e-15
+-1.33227e-15
+1.22125e-15
+-3.77476e-15
+1.22125e-15
+2.22045e-15
+-3.10862e-15
+-8.88178e-16
+-1.9984e-15
+-3.33067e-15
+-2.22045e-15
+4.44089e-16
+-4.44089e-16
+-3.9968e-15
+2.22045e-15
+-2.22045e-15
+-3.9968e-15
+1.44329e-15
+3.33067e-16
+-3.9968e-15
+1.44329e-15
+-2.66454e-15
+-5.10703e-15
+2.10942e-15
+-4.44089e-16
+-8.88178e-16
+1.77636e-15
+-1.77636e-15
+2.22045e-16
+1.77636e-15
+-6.66134e-16
+-8.88178e-16
+2.44249e-15
+-4.44089e-16
+-4.44089e-16
+3.33067e-15
+-4.44089e-16
+-2.44249e-15
+9.99201e-16
+1.9984e-15
+-1.77636e-15
+1.44329e-15
+2.22045e-16
+-1.55431e-15
+2.33147e-15
+4.66294e-15
+1.55431e-15
+4.55191e-15
+4.55191e-15
+3.33067e-15
+7.77156e-15
+9.65894e-15
+7.54952e-15
+1.42109e-14
+1.18794e-14
+1.19904e-14
+2.0095e-14
+2.24265e-14
+2.15383e-14
+2.63123e-14
+2.65343e-14
+3.08642e-14
+3.76366e-14
+4.55191e-14
+5.05151e-14
+5.68434e-14
+6.15064e-14
+6.59472e-14
+7.79377e-14
+9.10383e-14
+1.02474e-13
+1.0858e-13
+1.26565e-13
+1.39555e-13
+1.61204e-13
+1.84075e-13
+2.03393e-13
+2.33369e-13
+2.58571e-13
+2.88214e-13
+3.28071e-13
+3.71814e-13
+4.15223e-13
+4.62408e-13
+5.24469e-13
+5.86198e-13
+6.60916e-13
+7.48401e-13
+8.41216e-13
+9.40026e-13
+1.05582e-12
+1.18039e-12
+1.32683e-12
+1.49503e-12
+1.67766e-12
+1.87683e-12
+2.10443e-12
+2.35545e-12
+2.64166e-12
+2.96663e-12
+3.32345e-12
+3.71791e-12
+4.16345e-12
+4.65705e-12
+5.21705e-12
+5.84088e-12
+6.5361e-12
+7.30138e-12
+8.16924e-12
+9.12959e-12
+1.02066e-11
+1.14102e-11
+1.27501e-11
+1.42391e-11
+1.58993e-11
+1.77564e-11
+1.98178e-11
+2.21206e-11
+2.46846e-11
+2.75338e-11
+3.07125e-11
+3.42381e-11
+3.81701e-11
+4.2545e-11
+4.74119e-11
+5.28152e-11
+5.88221e-11
+6.55012e-11
+7.29256e-11
+8.11731e-11
+9.03341e-11
+1.00499e-10
+1.1179e-10
+1.24312e-10
+1.38213e-10
+1.53636e-10
+1.70746e-10
+1.89714e-10
+2.10735e-10
+2.34043e-10
+2.59864e-10
+2.88477e-10
+3.20167e-10
+3.55258e-10
+3.94105e-10
+4.37103e-10
+4.84683e-10
+5.3733e-10
+5.95561e-10
+6.59954e-10
+7.31141e-10
+8.09834e-10
+8.96792e-10
+9.92872e-10
+1.099e-09
+1.2162e-09
+1.3456e-09
+1.48843e-09
+1.64606e-09
+1.81998e-09
+2.01183e-09
+2.22341e-09
+2.45669e-09
+2.71384e-09
+2.99724e-09
+3.3095e-09
+3.65349e-09
+4.03233e-09
+4.44945e-09
+4.90864e-09
+5.41401e-09
+5.97009e-09
+6.58182e-09
+7.25461e-09
+7.99439e-09
+8.80766e-09
+9.70149e-09
+1.06837e-08
+1.17627e-08
+1.29477e-08
+1.4249e-08
+1.56777e-08
+1.72457e-08
+1.89663e-08
+2.08539e-08
+2.29243e-08
+2.51946e-08
+2.76837e-08
+3.04119e-08
+3.34014e-08
+3.66768e-08
+4.02643e-08
+4.41929e-08
+4.84941e-08
+5.3202e-08
+5.8354e-08
+6.39907e-08
+7.01562e-08
+7.68987e-08
+8.42704e-08
+9.23283e-08
+1.01134e-07
+1.10755e-07
+1.21265e-07
+1.32742e-07
+1.45273e-07
+1.58951e-07
+1.73879e-07
+1.90166e-07
+2.07933e-07
+2.27309e-07
+2.48436e-07
+2.71465e-07
+2.96564e-07
+3.2391e-07
+3.537e-07
+3.86144e-07
+4.2147e-07
+4.59925e-07
+5.01777e-07
+5.47316e-07
+5.96855e-07
+6.50734e-07
+7.09318e-07
+7.73004e-07
+8.42221e-07
+9.17431e-07
+9.99136e-07
+1.08788e-06
+1.18423e-06
+1.28884e-06
+1.40237e-06
+1.52557e-06
+1.65921e-06
+1.80417e-06
+1.96135e-06
+2.13175e-06
+2.31644e-06
+2.51658e-06
+2.73339e-06
+2.96822e-06
+3.22252e-06
+3.49782e-06
+3.79579e-06
+4.11823e-06
+4.46707e-06
+4.84438e-06
+5.25239e-06
+5.6935e-06
+6.17028e-06
+6.6855e-06
+7.24213e-06
+7.84336e-06
+8.49261e-06
+9.19356e-06
+9.95015e-06
+1.07666e-05
+1.16475e-05
+1.25976e-05
+1.36222e-05
+1.47269e-05
+1.59176e-05
+1.72007e-05
+1.85831e-05
+2.00722e-05
+2.16758e-05
+2.34023e-05
+2.52607e-05
+2.72606e-05
+2.94122e-05
+3.17267e-05
+3.42157e-05
+3.68917e-05
+3.97681e-05
+4.28593e-05
+4.61805e-05
+4.97481e-05
+5.35792e-05
+5.76926e-05
+6.2108e-05
+6.68464e-05
+7.19303e-05
+7.73837e-05
+8.3232e-05
+8.95023e-05
+9.62237e-05
+0.000103427
+0.000111144
+0.000119411
+0.000128264
+0.000137743
+0.00014789
+0.000158748
+0.000170366
+0.000182794
+0.000196085
+0.000210295
+0.000225484
+0.000241717
+0.000259061
+0.000277588
+0.000297374
+0.000318499
+0.000341049
+0.000365114
+0.00039079
+0.000418179
+0.000447388
+0.00047853
+0.000511727
+0.000547105
+0.000584798
+0.00062495
+0.000667709
+0.000713236
+0.000761697
+0.00081327
+0.000868142
+0.00092651
+0.000988582
+0.00105458
+0.00112473
+0.00119928
+0.00127849
+0.00136263
+0.00145198
+0.00154684
+0.00164754
+0.0017544
+0.00186778
+0.00198804
+0.00211558
+0.00225079
+0.00239412
+0.00254601
+0.00270692
+0.00287737
+0.00305788
+0.00324898
+0.00345126
+0.00366531
+0.00389178
+0.00413132
+0.00438463
+0.00465244
+0.0049355
+0.00523462
+0.00555064
+0.00588442
+0.00623689
+0.006609
+0.00700175
+0.00741619
+0.00785342
+0.00831458
+0.00880085
+0.0093135
+0.00985381
+0.0104231
+0.0110229
+0.0116546
+0.0123198
+0.01302
+0.013757
+0.0145324
+0.0153482
+0.0162061
+0.0171082
+0.0180565
+0.0190531
+0.0201002
+0.0212002
+0.0223553
+0.0235682
+0.0248414
+0.0261775
+0.0275793
+0.0290498
+0.0305918
+0.0322086
+0.0339032
+0.0356791
+0.0375396
+0.0394884
+0.041529
+0.0436655
+0.0459016
+0.0482415
+0.0506894
+0.0532497
+0.0559268
+0.0587255
+0.0616505
+0.0647068
+0.0678995
+0.0712339
+0.0747154
+0.0783496
+0.0821424
+0.0860995
+0.0902273
+0.0945319
+0.0990198
+0.103698
+0.108572
+0.113651
+0.118941
+0.124449
+0.130184
+0.136152
+0.142362
+0.148823
+0.155542
+0.162528
+0.16979
+0.177337
+0.185179
+0.193324
+0.201782
+0.210564
+0.21968
+0.229138
+0.238951
+0.249129
+0.259682
+0.270623
+0.281961
+0.293709
+0.305879
+0.318482
+0.331531
+0.345037
+0.359014
+0.373474
+0.388429
+0.403894
+0.419881
+0.436404
+0.453476
+0.471112
+0.489324
+0.508127
+0.527535
+0.547563
+0.568225
+0.589535
+0.611508
+0.63416
+0.657504
+0.681556
+0.70633
+0.731842
+0.758107
+0.78514
+0.812956
+0.84157
+0.870998
+0.901254
+0.932353
+0.964312
+0.997144
+1.03086
+1.06549
+1.10103
+1.1375
+1.17492
+1.21331
+1.25266
+1.29301
+1.33436
+1.37672
+1.42012
+1.46455
+1.51004
+1.5566
+1.60423
+1.65296
+1.70279
+1.75372
+1.80578
+1.85898
+1.91331
+1.96879
+2.02543
+2.08324
+2.14222
+2.20238
+2.26373
+2.32627
+2.39
+2.45494
+2.52107
+2.58842
+2.65697
+2.72673
+2.79769
+2.86987
+2.94326
+3.01785
+3.09364
+3.17063
+3.24881
+3.32818
+3.40873
+3.49045
+3.57334
+3.65738
+3.74257
+3.82889
+3.91633
+4.00487
+4.09451
+4.18522
+4.27699
+4.36979
+4.46362
+4.55845
+4.65426
+4.75103
+4.84872
+4.94733
+5.04682
+5.14716
+5.24834
+5.35031
+5.45305
+5.55652
+5.6607
+5.76555
+5.87103
+5.97712
+6.08377
+6.19094
+6.2986
+6.40671
+6.51523
+6.6241
+6.7333
+6.84278
+6.95249
+7.06239
+7.17243
+7.28256
+7.39274
+7.50292
+7.61305
+7.72307
+7.83294
+7.94261
+8.05202
+8.16113
+8.26987
+8.37819
+8.48605
+8.59339
+8.70014
+8.80626
+8.9117
+9.01639
+9.12028
+9.22332
+9.32545
+9.42661
+9.52674
+9.6258
+9.72373
+9.82047
+9.91596
+10.0102
+10.103
+10.1944
+10.2844
+10.3729
+10.4598
+10.545
+10.6286
+10.7105
+10.7907
+10.869
+10.9454
+11.0199
+11.0925
+11.163
+11.2315
+11.2979
+11.3622
+11.4243
+11.4842
+11.5418
+11.5971
+11.6501
+11.7008
+11.749
+11.7949
+11.8382
+11.8791
+11.9175
+11.9534
+11.9866
+12.0173
+12.0454
+12.0709
+12.0938
+12.114
+12.1315
+12.1463
+12.1585
+12.168
+12.1748
+12.1788
+12.1802
diff --git a/examples/DDG4/data/Si220pl/Si220pl_efx.txt b/examples/DDG4/data/Si220pl/Si220pl_efx.txt
new file mode 100644
index 000000000..38e0ccf8f
--- /dev/null
+++ b/examples/DDG4/data/Si220pl/Si220pl_efx.txt
@@ -0,0 +1,2050 @@
+2048 1 1
+1.91979e-10 2.21679e-10 3.13501e-10
+-0.00123105
+9.5106e+09
+1.90184e+10
+2.85204e+10
+3.8014e+10
+4.74963e+10
+5.69643e+10
+6.64154e+10
+7.58468e+10
+8.52556e+10
+9.46391e+10
+1.03994e+11
+1.13319e+11
+1.2261e+11
+1.31865e+11
+1.41081e+11
+1.50256e+11
+1.59386e+11
+1.6847e+11
+1.77504e+11
+1.86487e+11
+1.95416e+11
+2.04288e+11
+2.13101e+11
+2.21852e+11
+2.3054e+11
+2.39162e+11
+2.47715e+11
+2.56199e+11
+2.64609e+11
+2.72945e+11
+2.81204e+11
+2.89385e+11
+2.97485e+11
+3.05502e+11
+3.13435e+11
+3.21282e+11
+3.29041e+11
+3.3671e+11
+3.44288e+11
+3.51774e+11
+3.59165e+11
+3.66461e+11
+3.73659e+11
+3.80759e+11
+3.8776e+11
+3.94659e+11
+4.01457e+11
+4.08152e+11
+4.14742e+11
+4.21227e+11
+4.27607e+11
+4.33879e+11
+4.40044e+11
+4.46101e+11
+4.52048e+11
+4.57886e+11
+4.63614e+11
+4.69232e+11
+4.74738e+11
+4.80132e+11
+4.85415e+11
+4.90587e+11
+4.95645e+11
+5.00592e+11
+5.05426e+11
+5.10148e+11
+5.14757e+11
+5.19254e+11
+5.23639e+11
+5.27912e+11
+5.32073e+11
+5.36123e+11
+5.40062e+11
+5.43889e+11
+5.47607e+11
+5.51215e+11
+5.54713e+11
+5.58103e+11
+5.61384e+11
+5.64558e+11
+5.67625e+11
+5.70586e+11
+5.73442e+11
+5.76194e+11
+5.78841e+11
+5.81386e+11
+5.83829e+11
+5.86172e+11
+5.88414e+11
+5.90557e+11
+5.92602e+11
+5.94551e+11
+5.96404e+11
+5.98162e+11
+5.99826e+11
+6.01398e+11
+6.02879e+11
+6.0427e+11
+6.05573e+11
+6.06787e+11
+6.07916e+11
+6.08959e+11
+6.09918e+11
+6.10796e+11
+6.11591e+11
+6.12308e+11
+6.12945e+11
+6.13505e+11
+6.1399e+11
+6.144e+11
+6.14736e+11
+6.15001e+11
+6.15196e+11
+6.15321e+11
+6.15379e+11
+6.1537e+11
+6.15296e+11
+6.15159e+11
+6.14959e+11
+6.14698e+11
+6.14377e+11
+6.13999e+11
+6.13563e+11
+6.13072e+11
+6.12526e+11
+6.11927e+11
+6.11277e+11
+6.10576e+11
+6.09826e+11
+6.09029e+11
+6.08185e+11
+6.07295e+11
+6.06362e+11
+6.05385e+11
+6.04367e+11
+6.03309e+11
+6.02211e+11
+6.01075e+11
+5.99903e+11
+5.98694e+11
+5.97451e+11
+5.96175e+11
+5.94866e+11
+5.93526e+11
+5.92155e+11
+5.90756e+11
+5.89328e+11
+5.87873e+11
+5.86392e+11
+5.84885e+11
+5.83355e+11
+5.81801e+11
+5.80225e+11
+5.78627e+11
+5.77009e+11
+5.75372e+11
+5.73715e+11
+5.72041e+11
+5.70349e+11
+5.68642e+11
+5.66918e+11
+5.6518e+11
+5.63428e+11
+5.61663e+11
+5.59885e+11
+5.58096e+11
+5.56295e+11
+5.54484e+11
+5.52663e+11
+5.50833e+11
+5.48994e+11
+5.47148e+11
+5.45294e+11
+5.43433e+11
+5.41567e+11
+5.39694e+11
+5.37817e+11
+5.35935e+11
+5.34048e+11
+5.32159e+11
+5.30266e+11
+5.2837e+11
+5.26472e+11
+5.24572e+11
+5.22671e+11
+5.20769e+11
+5.18866e+11
+5.16964e+11
+5.15061e+11
+5.13159e+11
+5.11257e+11
+5.09357e+11
+5.07458e+11
+5.05562e+11
+5.03667e+11
+5.01774e+11
+4.99885e+11
+4.97998e+11
+4.96114e+11
+4.94234e+11
+4.92358e+11
+4.90485e+11
+4.88617e+11
+4.86753e+11
+4.84894e+11
+4.83039e+11
+4.81189e+11
+4.79344e+11
+4.77505e+11
+4.75671e+11
+4.73843e+11
+4.7202e+11
+4.70203e+11
+4.68392e+11
+4.66588e+11
+4.6479e+11
+4.62998e+11
+4.61212e+11
+4.59433e+11
+4.57661e+11
+4.55896e+11
+4.54138e+11
+4.52386e+11
+4.50642e+11
+4.48905e+11
+4.47174e+11
+4.45452e+11
+4.43736e+11
+4.42028e+11
+4.40327e+11
+4.38634e+11
+4.36949e+11
+4.35271e+11
+4.336e+11
+4.31937e+11
+4.30282e+11
+4.28635e+11
+4.26995e+11
+4.25363e+11
+4.23739e+11
+4.22123e+11
+4.20514e+11
+4.18913e+11
+4.1732e+11
+4.15735e+11
+4.14158e+11
+4.12589e+11
+4.11027e+11
+4.09474e+11
+4.07928e+11
+4.0639e+11
+4.0486e+11
+4.03337e+11
+4.01823e+11
+4.00316e+11
+3.98818e+11
+3.97327e+11
+3.95843e+11
+3.94368e+11
+3.929e+11
+3.9144e+11
+3.89988e+11
+3.88543e+11
+3.87106e+11
+3.85677e+11
+3.84256e+11
+3.82842e+11
+3.81435e+11
+3.80036e+11
+3.78645e+11
+3.77261e+11
+3.75885e+11
+3.74516e+11
+3.73154e+11
+3.718e+11
+3.70453e+11
+3.69114e+11
+3.67782e+11
+3.66457e+11
+3.65139e+11
+3.63829e+11
+3.62526e+11
+3.6123e+11
+3.59941e+11
+3.58659e+11
+3.57384e+11
+3.56116e+11
+3.54855e+11
+3.53601e+11
+3.52354e+11
+3.51114e+11
+3.4988e+11
+3.48654e+11
+3.47434e+11
+3.46221e+11
+3.45015e+11
+3.43815e+11
+3.42622e+11
+3.41435e+11
+3.40255e+11
+3.39082e+11
+3.37915e+11
+3.36754e+11
+3.356e+11
+3.34452e+11
+3.33311e+11
+3.32176e+11
+3.31047e+11
+3.29924e+11
+3.28808e+11
+3.27697e+11
+3.26593e+11
+3.25495e+11
+3.24403e+11
+3.23316e+11
+3.22236e+11
+3.21162e+11
+3.20094e+11
+3.19031e+11
+3.17974e+11
+3.16923e+11
+3.15878e+11
+3.14838e+11
+3.13805e+11
+3.12776e+11
+3.11754e+11
+3.10737e+11
+3.09725e+11
+3.08719e+11
+3.07718e+11
+3.06723e+11
+3.05733e+11
+3.04748e+11
+3.03769e+11
+3.02795e+11
+3.01827e+11
+3.00863e+11
+2.99905e+11
+2.98951e+11
+2.98003e+11
+2.9706e+11
+2.96122e+11
+2.95189e+11
+2.94261e+11
+2.93338e+11
+2.92419e+11
+2.91506e+11
+2.90597e+11
+2.89693e+11
+2.88794e+11
+2.879e+11
+2.8701e+11
+2.86125e+11
+2.85244e+11
+2.84368e+11
+2.83497e+11
+2.8263e+11
+2.81768e+11
+2.8091e+11
+2.80056e+11
+2.79207e+11
+2.78363e+11
+2.77522e+11
+2.76686e+11
+2.75854e+11
+2.75026e+11
+2.74203e+11
+2.73384e+11
+2.72568e+11
+2.71757e+11
+2.7095e+11
+2.70147e+11
+2.69348e+11
+2.68553e+11
+2.67762e+11
+2.66975e+11
+2.66192e+11
+2.65412e+11
+2.64637e+11
+2.63865e+11
+2.63097e+11
+2.62333e+11
+2.61572e+11
+2.60815e+11
+2.60062e+11
+2.59312e+11
+2.58566e+11
+2.57824e+11
+2.57085e+11
+2.5635e+11
+2.55618e+11
+2.54889e+11
+2.54164e+11
+2.53442e+11
+2.52724e+11
+2.52009e+11
+2.51298e+11
+2.50589e+11
+2.49884e+11
+2.49183e+11
+2.48484e+11
+2.47789e+11
+2.47096e+11
+2.46407e+11
+2.45721e+11
+2.45038e+11
+2.44358e+11
+2.43682e+11
+2.43008e+11
+2.42337e+11
+2.41669e+11
+2.41004e+11
+2.40342e+11
+2.39683e+11
+2.39027e+11
+2.38373e+11
+2.37723e+11
+2.37075e+11
+2.3643e+11
+2.35788e+11
+2.35148e+11
+2.34511e+11
+2.33877e+11
+2.33246e+11
+2.32617e+11
+2.3199e+11
+2.31367e+11
+2.30745e+11
+2.30127e+11
+2.29511e+11
+2.28897e+11
+2.28286e+11
+2.27677e+11
+2.27071e+11
+2.26468e+11
+2.25866e+11
+2.25267e+11
+2.2467e+11
+2.24076e+11
+2.23484e+11
+2.22894e+11
+2.22307e+11
+2.21722e+11
+2.21139e+11
+2.20558e+11
+2.1998e+11
+2.19403e+11
+2.18829e+11
+2.18257e+11
+2.17687e+11
+2.17119e+11
+2.16554e+11
+2.1599e+11
+2.15428e+11
+2.14869e+11
+2.14311e+11
+2.13755e+11
+2.13202e+11
+2.1265e+11
+2.121e+11
+2.11553e+11
+2.11007e+11
+2.10463e+11
+2.09921e+11
+2.0938e+11
+2.08842e+11
+2.08305e+11
+2.0777e+11
+2.07237e+11
+2.06706e+11
+2.06177e+11
+2.05649e+11
+2.05123e+11
+2.04599e+11
+2.04076e+11
+2.03555e+11
+2.03036e+11
+2.02518e+11
+2.02002e+11
+2.01488e+11
+2.00975e+11
+2.00464e+11
+1.99954e+11
+1.99446e+11
+1.9894e+11
+1.98435e+11
+1.97932e+11
+1.9743e+11
+1.96929e+11
+1.9643e+11
+1.95933e+11
+1.95437e+11
+1.94942e+11
+1.94449e+11
+1.93957e+11
+1.93467e+11
+1.92978e+11
+1.92491e+11
+1.92004e+11
+1.9152e+11
+1.91036e+11
+1.90554e+11
+1.90073e+11
+1.89593e+11
+1.89115e+11
+1.88638e+11
+1.88162e+11
+1.87688e+11
+1.87215e+11
+1.86743e+11
+1.86272e+11
+1.85802e+11
+1.85334e+11
+1.84866e+11
+1.844e+11
+1.83936e+11
+1.83472e+11
+1.83009e+11
+1.82548e+11
+1.82087e+11
+1.81628e+11
+1.8117e+11
+1.80713e+11
+1.80257e+11
+1.79802e+11
+1.79348e+11
+1.78896e+11
+1.78444e+11
+1.77993e+11
+1.77544e+11
+1.77095e+11
+1.76647e+11
+1.76201e+11
+1.75755e+11
+1.7531e+11
+1.74866e+11
+1.74424e+11
+1.73982e+11
+1.73541e+11
+1.73101e+11
+1.72662e+11
+1.72224e+11
+1.71786e+11
+1.7135e+11
+1.70915e+11
+1.7048e+11
+1.70046e+11
+1.69613e+11
+1.69181e+11
+1.6875e+11
+1.6832e+11
+1.6789e+11
+1.67462e+11
+1.67034e+11
+1.66607e+11
+1.66181e+11
+1.65755e+11
+1.65331e+11
+1.64907e+11
+1.64484e+11
+1.64061e+11
+1.6364e+11
+1.63219e+11
+1.62799e+11
+1.62379e+11
+1.61961e+11
+1.61543e+11
+1.61126e+11
+1.60709e+11
+1.60293e+11
+1.59878e+11
+1.59464e+11
+1.5905e+11
+1.58637e+11
+1.58225e+11
+1.57813e+11
+1.57402e+11
+1.56992e+11
+1.56582e+11
+1.56173e+11
+1.55764e+11
+1.55357e+11
+1.54949e+11
+1.54543e+11
+1.54137e+11
+1.53731e+11
+1.53327e+11
+1.52922e+11
+1.52519e+11
+1.52116e+11
+1.51713e+11
+1.51311e+11
+1.5091e+11
+1.50509e+11
+1.50109e+11
+1.4971e+11
+1.4931e+11
+1.48912e+11
+1.48514e+11
+1.48116e+11
+1.47719e+11
+1.47323e+11
+1.46927e+11
+1.46532e+11
+1.46137e+11
+1.45743e+11
+1.45349e+11
+1.44955e+11
+1.44562e+11
+1.4417e+11
+1.43778e+11
+1.43387e+11
+1.42996e+11
+1.42605e+11
+1.42215e+11
+1.41826e+11
+1.41436e+11
+1.41048e+11
+1.4066e+11
+1.40272e+11
+1.39884e+11
+1.39498e+11
+1.39111e+11
+1.38725e+11
+1.3834e+11
+1.37954e+11
+1.3757e+11
+1.37185e+11
+1.36801e+11
+1.36418e+11
+1.36035e+11
+1.35652e+11
+1.3527e+11
+1.34888e+11
+1.34506e+11
+1.34125e+11
+1.33744e+11
+1.33364e+11
+1.32984e+11
+1.32604e+11
+1.32225e+11
+1.31846e+11
+1.31468e+11
+1.31089e+11
+1.30712e+11
+1.30334e+11
+1.29957e+11
+1.2958e+11
+1.29204e+11
+1.28828e+11
+1.28452e+11
+1.28077e+11
+1.27702e+11
+1.27327e+11
+1.26952e+11
+1.26578e+11
+1.26205e+11
+1.25831e+11
+1.25458e+11
+1.25085e+11
+1.24713e+11
+1.24341e+11
+1.23969e+11
+1.23597e+11
+1.23226e+11
+1.22855e+11
+1.22485e+11
+1.22114e+11
+1.21744e+11
+1.21375e+11
+1.21005e+11
+1.20636e+11
+1.20267e+11
+1.19898e+11
+1.1953e+11
+1.19162e+11
+1.18794e+11
+1.18427e+11
+1.1806e+11
+1.17693e+11
+1.17326e+11
+1.1696e+11
+1.16593e+11
+1.16227e+11
+1.15862e+11
+1.15496e+11
+1.15131e+11
+1.14766e+11
+1.14402e+11
+1.14038e+11
+1.13673e+11
+1.1331e+11
+1.12946e+11
+1.12583e+11
+1.12219e+11
+1.11857e+11
+1.11494e+11
+1.11131e+11
+1.10769e+11
+1.10407e+11
+1.10046e+11
+1.09684e+11
+1.09323e+11
+1.08962e+11
+1.08601e+11
+1.0824e+11
+1.0788e+11
+1.0752e+11
+1.0716e+11
+1.068e+11
+1.06441e+11
+1.06081e+11
+1.05722e+11
+1.05363e+11
+1.05005e+11
+1.04646e+11
+1.04288e+11
+1.0393e+11
+1.03572e+11
+1.03215e+11
+1.02857e+11
+1.025e+11
+1.02143e+11
+1.01786e+11
+1.01429e+11
+1.01073e+11
+1.00717e+11
+1.00361e+11
+1.00005e+11
+9.9649e+10
+9.92935e+10
+9.89381e+10
+9.8583e+10
+9.82281e+10
+9.78733e+10
+9.75187e+10
+9.71643e+10
+9.68101e+10
+9.6456e+10
+9.61022e+10
+9.57485e+10
+9.53949e+10
+9.50416e+10
+9.46884e+10
+9.43354e+10
+9.39826e+10
+9.36299e+10
+9.32774e+10
+9.29251e+10
+9.25729e+10
+9.22209e+10
+9.18691e+10
+9.15174e+10
+9.11659e+10
+9.08145e+10
+9.04633e+10
+9.01123e+10
+8.97614e+10
+8.94107e+10
+8.90601e+10
+8.87097e+10
+8.83595e+10
+8.80094e+10
+8.76594e+10
+8.73096e+10
+8.696e+10
+8.66105e+10
+8.62611e+10
+8.59119e+10
+8.55629e+10
+8.5214e+10
+8.48652e+10
+8.45166e+10
+8.41681e+10
+8.38198e+10
+8.34716e+10
+8.31235e+10
+8.27756e+10
+8.24278e+10
+8.20802e+10
+8.17327e+10
+8.13853e+10
+8.10381e+10
+8.0691e+10
+8.0344e+10
+7.99972e+10
+7.96505e+10
+7.93039e+10
+7.89575e+10
+7.86112e+10
+7.8265e+10
+7.7919e+10
+7.75731e+10
+7.72273e+10
+7.68816e+10
+7.65361e+10
+7.61907e+10
+7.58454e+10
+7.55002e+10
+7.51552e+10
+7.48103e+10
+7.44655e+10
+7.41208e+10
+7.37762e+10
+7.34318e+10
+7.30875e+10
+7.27433e+10
+7.23992e+10
+7.20552e+10
+7.17114e+10
+7.13676e+10
+7.1024e+10
+7.06805e+10
+7.03371e+10
+6.99938e+10
+6.96506e+10
+6.93076e+10
+6.89646e+10
+6.86218e+10
+6.82791e+10
+6.79364e+10
+6.75939e+10
+6.72515e+10
+6.69092e+10
+6.6567e+10
+6.62249e+10
+6.58829e+10
+6.55411e+10
+6.51993e+10
+6.48576e+10
+6.4516e+10
+6.41746e+10
+6.38332e+10
+6.34919e+10
+6.31507e+10
+6.28097e+10
+6.24687e+10
+6.21278e+10
+6.1787e+10
+6.14464e+10
+6.11058e+10
+6.07653e+10
+6.04249e+10
+6.00846e+10
+5.97444e+10
+5.94043e+10
+5.90642e+10
+5.87243e+10
+5.83845e+10
+5.80447e+10
+5.77051e+10
+5.73655e+10
+5.7026e+10
+5.66866e+10
+5.63473e+10
+5.60081e+10
+5.5669e+10
+5.533e+10
+5.4991e+10
+5.46521e+10
+5.43134e+10
+5.39747e+10
+5.3636e+10
+5.32975e+10
+5.2959e+10
+5.26207e+10
+5.22824e+10
+5.19442e+10
+5.16061e+10
+5.1268e+10
+5.093e+10
+5.05921e+10
+5.02543e+10
+4.99166e+10
+4.95789e+10
+4.92414e+10
+4.89039e+10
+4.85664e+10
+4.82291e+10
+4.78918e+10
+4.75546e+10
+4.72174e+10
+4.68804e+10
+4.65434e+10
+4.62065e+10
+4.58696e+10
+4.55329e+10
+4.51961e+10
+4.48595e+10
+4.45229e+10
+4.41864e+10
+4.385e+10
+4.35136e+10
+4.31773e+10
+4.28411e+10
+4.2505e+10
+4.21688e+10
+4.18328e+10
+4.14968e+10
+4.11609e+10
+4.08251e+10
+4.04893e+10
+4.01536e+10
+3.98179e+10
+3.94823e+10
+3.91468e+10
+3.88113e+10
+3.84759e+10
+3.81405e+10
+3.78052e+10
+3.747e+10
+3.71348e+10
+3.67997e+10
+3.64646e+10
+3.61296e+10
+3.57946e+10
+3.54597e+10
+3.51249e+10
+3.47901e+10
+3.44553e+10
+3.41206e+10
+3.3786e+10
+3.34514e+10
+3.31169e+10
+3.27824e+10
+3.24479e+10
+3.21136e+10
+3.17792e+10
+3.14449e+10
+3.11107e+10
+3.07765e+10
+3.04424e+10
+3.01083e+10
+2.97742e+10
+2.94402e+10
+2.91062e+10
+2.87723e+10
+2.84384e+10
+2.81046e+10
+2.77708e+10
+2.74371e+10
+2.71034e+10
+2.67697e+10
+2.64361e+10
+2.61025e+10
+2.5769e+10
+2.54355e+10
+2.5102e+10
+2.47686e+10
+2.44352e+10
+2.41018e+10
+2.37685e+10
+2.34353e+10
+2.3102e+10
+2.27688e+10
+2.24356e+10
+2.21025e+10
+2.17694e+10
+2.14363e+10
+2.11033e+10
+2.07703e+10
+2.04373e+10
+2.01044e+10
+1.97715e+10
+1.94386e+10
+1.91058e+10
+1.87729e+10
+1.84401e+10
+1.81074e+10
+1.77746e+10
+1.74419e+10
+1.71093e+10
+1.67766e+10
+1.6444e+10
+1.61114e+10
+1.57788e+10
+1.54462e+10
+1.51137e+10
+1.47812e+10
+1.44487e+10
+1.41162e+10
+1.37838e+10
+1.34514e+10
+1.3119e+10
+1.27866e+10
+1.24542e+10
+1.21219e+10
+1.17896e+10
+1.14573e+10
+1.1125e+10
+1.07927e+10
+1.04605e+10
+1.01282e+10
+9.79599e+09
+9.46379e+09
+9.13159e+09
+8.79941e+09
+8.46725e+09
+8.1351e+09
+7.80296e+09
+7.47083e+09
+7.13871e+09
+6.8066e+09
+6.47451e+09
+6.14242e+09
+5.81035e+09
+5.47828e+09
+5.14622e+09
+4.81417e+09
+4.48213e+09
+4.15009e+09
+3.81806e+09
+3.48603e+09
+3.15401e+09
+2.822e+09
+2.48999e+09
+2.15798e+09
+1.82598e+09
+1.49398e+09
+1.16198e+09
+8.29986e+08
+4.97991e+08
+1.65997e+08
+-1.65997e+08
+-4.97991e+08
+-8.29986e+08
+-1.16198e+09
+-1.49398e+09
+-1.82598e+09
+-2.15798e+09
+-2.48999e+09
+-2.822e+09
+-3.15401e+09
+-3.48603e+09
+-3.81806e+09
+-4.15009e+09
+-4.48213e+09
+-4.81417e+09
+-5.14622e+09
+-5.47828e+09
+-5.81035e+09
+-6.14242e+09
+-6.47451e+09
+-6.8066e+09
+-7.13871e+09
+-7.47083e+09
+-7.80296e+09
+-8.1351e+09
+-8.46725e+09
+-8.79941e+09
+-9.13159e+09
+-9.46379e+09
+-9.79599e+09
+-1.01282e+10
+-1.04605e+10
+-1.07927e+10
+-1.1125e+10
+-1.14573e+10
+-1.17896e+10
+-1.21219e+10
+-1.24542e+10
+-1.27866e+10
+-1.3119e+10
+-1.34514e+10
+-1.37838e+10
+-1.41162e+10
+-1.44487e+10
+-1.47812e+10
+-1.51137e+10
+-1.54462e+10
+-1.57788e+10
+-1.61114e+10
+-1.6444e+10
+-1.67766e+10
+-1.71093e+10
+-1.74419e+10
+-1.77746e+10
+-1.81074e+10
+-1.84401e+10
+-1.87729e+10
+-1.91058e+10
+-1.94386e+10
+-1.97715e+10
+-2.01044e+10
+-2.04373e+10
+-2.07703e+10
+-2.11033e+10
+-2.14363e+10
+-2.17694e+10
+-2.21025e+10
+-2.24356e+10
+-2.27688e+10
+-2.3102e+10
+-2.34353e+10
+-2.37685e+10
+-2.41018e+10
+-2.44352e+10
+-2.47686e+10
+-2.5102e+10
+-2.54355e+10
+-2.5769e+10
+-2.61025e+10
+-2.64361e+10
+-2.67697e+10
+-2.71034e+10
+-2.74371e+10
+-2.77708e+10
+-2.81046e+10
+-2.84384e+10
+-2.87723e+10
+-2.91062e+10
+-2.94402e+10
+-2.97742e+10
+-3.01083e+10
+-3.04424e+10
+-3.07765e+10
+-3.11107e+10
+-3.14449e+10
+-3.17792e+10
+-3.21136e+10
+-3.24479e+10
+-3.27824e+10
+-3.31169e+10
+-3.34514e+10
+-3.3786e+10
+-3.41206e+10
+-3.44553e+10
+-3.47901e+10
+-3.51249e+10
+-3.54597e+10
+-3.57946e+10
+-3.61296e+10
+-3.64646e+10
+-3.67997e+10
+-3.71348e+10
+-3.747e+10
+-3.78052e+10
+-3.81405e+10
+-3.84759e+10
+-3.88113e+10
+-3.91468e+10
+-3.94823e+10
+-3.98179e+10
+-4.01536e+10
+-4.04893e+10
+-4.08251e+10
+-4.11609e+10
+-4.14968e+10
+-4.18328e+10
+-4.21688e+10
+-4.2505e+10
+-4.28411e+10
+-4.31773e+10
+-4.35136e+10
+-4.385e+10
+-4.41864e+10
+-4.45229e+10
+-4.48595e+10
+-4.51961e+10
+-4.55329e+10
+-4.58696e+10
+-4.62065e+10
+-4.65434e+10
+-4.68804e+10
+-4.72174e+10
+-4.75546e+10
+-4.78918e+10
+-4.82291e+10
+-4.85664e+10
+-4.89039e+10
+-4.92414e+10
+-4.95789e+10
+-4.99166e+10
+-5.02543e+10
+-5.05921e+10
+-5.093e+10
+-5.1268e+10
+-5.16061e+10
+-5.19442e+10
+-5.22824e+10
+-5.26207e+10
+-5.2959e+10
+-5.32975e+10
+-5.3636e+10
+-5.39747e+10
+-5.43134e+10
+-5.46521e+10
+-5.4991e+10
+-5.533e+10
+-5.5669e+10
+-5.60081e+10
+-5.63473e+10
+-5.66866e+10
+-5.7026e+10
+-5.73655e+10
+-5.77051e+10
+-5.80447e+10
+-5.83845e+10
+-5.87243e+10
+-5.90642e+10
+-5.94043e+10
+-5.97444e+10
+-6.00846e+10
+-6.04249e+10
+-6.07653e+10
+-6.11058e+10
+-6.14464e+10
+-6.1787e+10
+-6.21278e+10
+-6.24687e+10
+-6.28097e+10
+-6.31507e+10
+-6.34919e+10
+-6.38332e+10
+-6.41746e+10
+-6.4516e+10
+-6.48576e+10
+-6.51993e+10
+-6.55411e+10
+-6.58829e+10
+-6.62249e+10
+-6.6567e+10
+-6.69092e+10
+-6.72515e+10
+-6.75939e+10
+-6.79364e+10
+-6.82791e+10
+-6.86218e+10
+-6.89646e+10
+-6.93076e+10
+-6.96506e+10
+-6.99938e+10
+-7.03371e+10
+-7.06805e+10
+-7.1024e+10
+-7.13676e+10
+-7.17114e+10
+-7.20552e+10
+-7.23992e+10
+-7.27433e+10
+-7.30875e+10
+-7.34318e+10
+-7.37762e+10
+-7.41208e+10
+-7.44655e+10
+-7.48103e+10
+-7.51552e+10
+-7.55002e+10
+-7.58454e+10
+-7.61907e+10
+-7.65361e+10
+-7.68816e+10
+-7.72273e+10
+-7.75731e+10
+-7.7919e+10
+-7.8265e+10
+-7.86112e+10
+-7.89575e+10
+-7.93039e+10
+-7.96505e+10
+-7.99972e+10
+-8.0344e+10
+-8.0691e+10
+-8.10381e+10
+-8.13853e+10
+-8.17327e+10
+-8.20802e+10
+-8.24278e+10
+-8.27756e+10
+-8.31235e+10
+-8.34716e+10
+-8.38198e+10
+-8.41681e+10
+-8.45166e+10
+-8.48652e+10
+-8.5214e+10
+-8.55629e+10
+-8.59119e+10
+-8.62611e+10
+-8.66105e+10
+-8.696e+10
+-8.73096e+10
+-8.76594e+10
+-8.80094e+10
+-8.83595e+10
+-8.87097e+10
+-8.90601e+10
+-8.94107e+10
+-8.97614e+10
+-9.01123e+10
+-9.04633e+10
+-9.08145e+10
+-9.11659e+10
+-9.15174e+10
+-9.18691e+10
+-9.22209e+10
+-9.25729e+10
+-9.29251e+10
+-9.32774e+10
+-9.36299e+10
+-9.39826e+10
+-9.43354e+10
+-9.46884e+10
+-9.50416e+10
+-9.53949e+10
+-9.57485e+10
+-9.61022e+10
+-9.6456e+10
+-9.68101e+10
+-9.71643e+10
+-9.75187e+10
+-9.78733e+10
+-9.82281e+10
+-9.8583e+10
+-9.89381e+10
+-9.92935e+10
+-9.9649e+10
+-1.00005e+11
+-1.00361e+11
+-1.00717e+11
+-1.01073e+11
+-1.01429e+11
+-1.01786e+11
+-1.02143e+11
+-1.025e+11
+-1.02857e+11
+-1.03215e+11
+-1.03572e+11
+-1.0393e+11
+-1.04288e+11
+-1.04646e+11
+-1.05005e+11
+-1.05363e+11
+-1.05722e+11
+-1.06081e+11
+-1.06441e+11
+-1.068e+11
+-1.0716e+11
+-1.0752e+11
+-1.0788e+11
+-1.0824e+11
+-1.08601e+11
+-1.08962e+11
+-1.09323e+11
+-1.09684e+11
+-1.10046e+11
+-1.10407e+11
+-1.10769e+11
+-1.11131e+11
+-1.11494e+11
+-1.11857e+11
+-1.12219e+11
+-1.12583e+11
+-1.12946e+11
+-1.1331e+11
+-1.13673e+11
+-1.14038e+11
+-1.14402e+11
+-1.14766e+11
+-1.15131e+11
+-1.15496e+11
+-1.15862e+11
+-1.16227e+11
+-1.16593e+11
+-1.1696e+11
+-1.17326e+11
+-1.17693e+11
+-1.1806e+11
+-1.18427e+11
+-1.18794e+11
+-1.19162e+11
+-1.1953e+11
+-1.19898e+11
+-1.20267e+11
+-1.20636e+11
+-1.21005e+11
+-1.21375e+11
+-1.21744e+11
+-1.22114e+11
+-1.22485e+11
+-1.22855e+11
+-1.23226e+11
+-1.23597e+11
+-1.23969e+11
+-1.24341e+11
+-1.24713e+11
+-1.25085e+11
+-1.25458e+11
+-1.25831e+11
+-1.26205e+11
+-1.26578e+11
+-1.26952e+11
+-1.27327e+11
+-1.27702e+11
+-1.28077e+11
+-1.28452e+11
+-1.28828e+11
+-1.29204e+11
+-1.2958e+11
+-1.29957e+11
+-1.30334e+11
+-1.30712e+11
+-1.31089e+11
+-1.31468e+11
+-1.31846e+11
+-1.32225e+11
+-1.32604e+11
+-1.32984e+11
+-1.33364e+11
+-1.33744e+11
+-1.34125e+11
+-1.34506e+11
+-1.34888e+11
+-1.3527e+11
+-1.35652e+11
+-1.36035e+11
+-1.36418e+11
+-1.36801e+11
+-1.37185e+11
+-1.3757e+11
+-1.37954e+11
+-1.3834e+11
+-1.38725e+11
+-1.39111e+11
+-1.39498e+11
+-1.39884e+11
+-1.40272e+11
+-1.4066e+11
+-1.41048e+11
+-1.41436e+11
+-1.41826e+11
+-1.42215e+11
+-1.42605e+11
+-1.42996e+11
+-1.43387e+11
+-1.43778e+11
+-1.4417e+11
+-1.44562e+11
+-1.44955e+11
+-1.45349e+11
+-1.45743e+11
+-1.46137e+11
+-1.46532e+11
+-1.46927e+11
+-1.47323e+11
+-1.47719e+11
+-1.48116e+11
+-1.48514e+11
+-1.48912e+11
+-1.4931e+11
+-1.4971e+11
+-1.50109e+11
+-1.50509e+11
+-1.5091e+11
+-1.51311e+11
+-1.51713e+11
+-1.52116e+11
+-1.52519e+11
+-1.52922e+11
+-1.53327e+11
+-1.53731e+11
+-1.54137e+11
+-1.54543e+11
+-1.54949e+11
+-1.55357e+11
+-1.55764e+11
+-1.56173e+11
+-1.56582e+11
+-1.56992e+11
+-1.57402e+11
+-1.57813e+11
+-1.58225e+11
+-1.58637e+11
+-1.5905e+11
+-1.59464e+11
+-1.59878e+11
+-1.60293e+11
+-1.60709e+11
+-1.61126e+11
+-1.61543e+11
+-1.61961e+11
+-1.62379e+11
+-1.62799e+11
+-1.63219e+11
+-1.6364e+11
+-1.64061e+11
+-1.64484e+11
+-1.64907e+11
+-1.65331e+11
+-1.65755e+11
+-1.66181e+11
+-1.66607e+11
+-1.67034e+11
+-1.67462e+11
+-1.6789e+11
+-1.6832e+11
+-1.6875e+11
+-1.69181e+11
+-1.69613e+11
+-1.70046e+11
+-1.7048e+11
+-1.70915e+11
+-1.7135e+11
+-1.71786e+11
+-1.72224e+11
+-1.72662e+11
+-1.73101e+11
+-1.73541e+11
+-1.73982e+11
+-1.74424e+11
+-1.74866e+11
+-1.7531e+11
+-1.75755e+11
+-1.76201e+11
+-1.76647e+11
+-1.77095e+11
+-1.77544e+11
+-1.77993e+11
+-1.78444e+11
+-1.78896e+11
+-1.79348e+11
+-1.79802e+11
+-1.80257e+11
+-1.80713e+11
+-1.8117e+11
+-1.81628e+11
+-1.82087e+11
+-1.82548e+11
+-1.83009e+11
+-1.83472e+11
+-1.83936e+11
+-1.844e+11
+-1.84866e+11
+-1.85334e+11
+-1.85802e+11
+-1.86272e+11
+-1.86743e+11
+-1.87215e+11
+-1.87688e+11
+-1.88162e+11
+-1.88638e+11
+-1.89115e+11
+-1.89593e+11
+-1.90073e+11
+-1.90554e+11
+-1.91036e+11
+-1.9152e+11
+-1.92004e+11
+-1.92491e+11
+-1.92978e+11
+-1.93467e+11
+-1.93957e+11
+-1.94449e+11
+-1.94942e+11
+-1.95437e+11
+-1.95933e+11
+-1.9643e+11
+-1.96929e+11
+-1.9743e+11
+-1.97932e+11
+-1.98435e+11
+-1.9894e+11
+-1.99446e+11
+-1.99954e+11
+-2.00464e+11
+-2.00975e+11
+-2.01488e+11
+-2.02002e+11
+-2.02518e+11
+-2.03036e+11
+-2.03555e+11
+-2.04076e+11
+-2.04599e+11
+-2.05123e+11
+-2.05649e+11
+-2.06177e+11
+-2.06706e+11
+-2.07237e+11
+-2.0777e+11
+-2.08305e+11
+-2.08842e+11
+-2.0938e+11
+-2.09921e+11
+-2.10463e+11
+-2.11007e+11
+-2.11553e+11
+-2.121e+11
+-2.1265e+11
+-2.13202e+11
+-2.13755e+11
+-2.14311e+11
+-2.14869e+11
+-2.15428e+11
+-2.1599e+11
+-2.16554e+11
+-2.17119e+11
+-2.17687e+11
+-2.18257e+11
+-2.18829e+11
+-2.19403e+11
+-2.1998e+11
+-2.20558e+11
+-2.21139e+11
+-2.21722e+11
+-2.22307e+11
+-2.22894e+11
+-2.23484e+11
+-2.24076e+11
+-2.2467e+11
+-2.25267e+11
+-2.25866e+11
+-2.26468e+11
+-2.27071e+11
+-2.27677e+11
+-2.28286e+11
+-2.28897e+11
+-2.29511e+11
+-2.30127e+11
+-2.30745e+11
+-2.31367e+11
+-2.3199e+11
+-2.32617e+11
+-2.33246e+11
+-2.33877e+11
+-2.34511e+11
+-2.35148e+11
+-2.35788e+11
+-2.3643e+11
+-2.37075e+11
+-2.37723e+11
+-2.38373e+11
+-2.39027e+11
+-2.39683e+11
+-2.40342e+11
+-2.41004e+11
+-2.41669e+11
+-2.42337e+11
+-2.43008e+11
+-2.43682e+11
+-2.44358e+11
+-2.45038e+11
+-2.45721e+11
+-2.46407e+11
+-2.47096e+11
+-2.47789e+11
+-2.48484e+11
+-2.49183e+11
+-2.49884e+11
+-2.50589e+11
+-2.51298e+11
+-2.52009e+11
+-2.52724e+11
+-2.53442e+11
+-2.54164e+11
+-2.54889e+11
+-2.55618e+11
+-2.5635e+11
+-2.57085e+11
+-2.57824e+11
+-2.58566e+11
+-2.59312e+11
+-2.60062e+11
+-2.60815e+11
+-2.61572e+11
+-2.62333e+11
+-2.63097e+11
+-2.63865e+11
+-2.64637e+11
+-2.65412e+11
+-2.66192e+11
+-2.66975e+11
+-2.67762e+11
+-2.68553e+11
+-2.69348e+11
+-2.70147e+11
+-2.7095e+11
+-2.71757e+11
+-2.72568e+11
+-2.73384e+11
+-2.74203e+11
+-2.75026e+11
+-2.75854e+11
+-2.76686e+11
+-2.77522e+11
+-2.78363e+11
+-2.79207e+11
+-2.80056e+11
+-2.8091e+11
+-2.81768e+11
+-2.8263e+11
+-2.83497e+11
+-2.84368e+11
+-2.85244e+11
+-2.86125e+11
+-2.8701e+11
+-2.879e+11
+-2.88794e+11
+-2.89693e+11
+-2.90597e+11
+-2.91506e+11
+-2.92419e+11
+-2.93338e+11
+-2.94261e+11
+-2.95189e+11
+-2.96122e+11
+-2.9706e+11
+-2.98003e+11
+-2.98951e+11
+-2.99905e+11
+-3.00863e+11
+-3.01827e+11
+-3.02795e+11
+-3.03769e+11
+-3.04748e+11
+-3.05733e+11
+-3.06723e+11
+-3.07718e+11
+-3.08719e+11
+-3.09725e+11
+-3.10737e+11
+-3.11754e+11
+-3.12776e+11
+-3.13805e+11
+-3.14838e+11
+-3.15878e+11
+-3.16923e+11
+-3.17974e+11
+-3.19031e+11
+-3.20094e+11
+-3.21162e+11
+-3.22236e+11
+-3.23316e+11
+-3.24403e+11
+-3.25495e+11
+-3.26593e+11
+-3.27697e+11
+-3.28808e+11
+-3.29924e+11
+-3.31047e+11
+-3.32176e+11
+-3.33311e+11
+-3.34452e+11
+-3.356e+11
+-3.36754e+11
+-3.37915e+11
+-3.39082e+11
+-3.40255e+11
+-3.41435e+11
+-3.42622e+11
+-3.43815e+11
+-3.45015e+11
+-3.46221e+11
+-3.47434e+11
+-3.48654e+11
+-3.4988e+11
+-3.51114e+11
+-3.52354e+11
+-3.53601e+11
+-3.54855e+11
+-3.56116e+11
+-3.57384e+11
+-3.58659e+11
+-3.59941e+11
+-3.6123e+11
+-3.62526e+11
+-3.63829e+11
+-3.65139e+11
+-3.66457e+11
+-3.67782e+11
+-3.69114e+11
+-3.70453e+11
+-3.718e+11
+-3.73154e+11
+-3.74516e+11
+-3.75885e+11
+-3.77261e+11
+-3.78645e+11
+-3.80036e+11
+-3.81435e+11
+-3.82842e+11
+-3.84256e+11
+-3.85677e+11
+-3.87106e+11
+-3.88543e+11
+-3.89988e+11
+-3.9144e+11
+-3.929e+11
+-3.94368e+11
+-3.95843e+11
+-3.97327e+11
+-3.98818e+11
+-4.00316e+11
+-4.01823e+11
+-4.03337e+11
+-4.0486e+11
+-4.0639e+11
+-4.07928e+11
+-4.09474e+11
+-4.11027e+11
+-4.12589e+11
+-4.14158e+11
+-4.15735e+11
+-4.1732e+11
+-4.18913e+11
+-4.20514e+11
+-4.22123e+11
+-4.23739e+11
+-4.25363e+11
+-4.26995e+11
+-4.28635e+11
+-4.30282e+11
+-4.31937e+11
+-4.336e+11
+-4.35271e+11
+-4.36949e+11
+-4.38634e+11
+-4.40327e+11
+-4.42028e+11
+-4.43736e+11
+-4.45452e+11
+-4.47174e+11
+-4.48905e+11
+-4.50642e+11
+-4.52386e+11
+-4.54138e+11
+-4.55896e+11
+-4.57661e+11
+-4.59433e+11
+-4.61212e+11
+-4.62998e+11
+-4.6479e+11
+-4.66588e+11
+-4.68392e+11
+-4.70203e+11
+-4.7202e+11
+-4.73843e+11
+-4.75671e+11
+-4.77505e+11
+-4.79344e+11
+-4.81189e+11
+-4.83039e+11
+-4.84894e+11
+-4.86753e+11
+-4.88617e+11
+-4.90485e+11
+-4.92358e+11
+-4.94234e+11
+-4.96114e+11
+-4.97998e+11
+-4.99885e+11
+-5.01774e+11
+-5.03667e+11
+-5.05562e+11
+-5.07458e+11
+-5.09357e+11
+-5.11257e+11
+-5.13159e+11
+-5.15061e+11
+-5.16964e+11
+-5.18866e+11
+-5.20769e+11
+-5.22671e+11
+-5.24572e+11
+-5.26472e+11
+-5.2837e+11
+-5.30266e+11
+-5.32159e+11
+-5.34048e+11
+-5.35935e+11
+-5.37817e+11
+-5.39694e+11
+-5.41567e+11
+-5.43433e+11
+-5.45294e+11
+-5.47148e+11
+-5.48994e+11
+-5.50833e+11
+-5.52663e+11
+-5.54484e+11
+-5.56295e+11
+-5.58096e+11
+-5.59885e+11
+-5.61663e+11
+-5.63428e+11
+-5.6518e+11
+-5.66918e+11
+-5.68642e+11
+-5.70349e+11
+-5.72041e+11
+-5.73715e+11
+-5.75372e+11
+-5.77009e+11
+-5.78627e+11
+-5.80225e+11
+-5.81801e+11
+-5.83355e+11
+-5.84885e+11
+-5.86392e+11
+-5.87873e+11
+-5.89328e+11
+-5.90756e+11
+-5.92155e+11
+-5.93526e+11
+-5.94866e+11
+-5.96175e+11
+-5.97451e+11
+-5.98694e+11
+-5.99903e+11
+-6.01075e+11
+-6.02211e+11
+-6.03309e+11
+-6.04367e+11
+-6.05385e+11
+-6.06362e+11
+-6.07295e+11
+-6.08185e+11
+-6.09029e+11
+-6.09826e+11
+-6.10576e+11
+-6.11277e+11
+-6.11927e+11
+-6.12526e+11
+-6.13072e+11
+-6.13563e+11
+-6.13999e+11
+-6.14377e+11
+-6.14698e+11
+-6.14959e+11
+-6.15159e+11
+-6.15296e+11
+-6.1537e+11
+-6.15379e+11
+-6.15321e+11
+-6.15196e+11
+-6.15001e+11
+-6.14736e+11
+-6.144e+11
+-6.1399e+11
+-6.13505e+11
+-6.12945e+11
+-6.12308e+11
+-6.11591e+11
+-6.10796e+11
+-6.09918e+11
+-6.08959e+11
+-6.07916e+11
+-6.06787e+11
+-6.05573e+11
+-6.0427e+11
+-6.02879e+11
+-6.01398e+11
+-5.99826e+11
+-5.98162e+11
+-5.96404e+11
+-5.94551e+11
+-5.92602e+11
+-5.90557e+11
+-5.88414e+11
+-5.86172e+11
+-5.83829e+11
+-5.81386e+11
+-5.78841e+11
+-5.76194e+11
+-5.73442e+11
+-5.70586e+11
+-5.67625e+11
+-5.64558e+11
+-5.61384e+11
+-5.58103e+11
+-5.54713e+11
+-5.51215e+11
+-5.47607e+11
+-5.43889e+11
+-5.40062e+11
+-5.36123e+11
+-5.32073e+11
+-5.27912e+11
+-5.23639e+11
+-5.19254e+11
+-5.14757e+11
+-5.10148e+11
+-5.05426e+11
+-5.00592e+11
+-4.95645e+11
+-4.90587e+11
+-4.85415e+11
+-4.80132e+11
+-4.74738e+11
+-4.69232e+11
+-4.63614e+11
+-4.57886e+11
+-4.52048e+11
+-4.46101e+11
+-4.40044e+11
+-4.33879e+11
+-4.27607e+11
+-4.21227e+11
+-4.14742e+11
+-4.08152e+11
+-4.01457e+11
+-3.94659e+11
+-3.8776e+11
+-3.80759e+11
+-3.73659e+11
+-3.66461e+11
+-3.59165e+11
+-3.51774e+11
+-3.44288e+11
+-3.3671e+11
+-3.29041e+11
+-3.21282e+11
+-3.13435e+11
+-3.05502e+11
+-2.97485e+11
+-2.89385e+11
+-2.81204e+11
+-2.72945e+11
+-2.64609e+11
+-2.56199e+11
+-2.47715e+11
+-2.39162e+11
+-2.3054e+11
+-2.21852e+11
+-2.13101e+11
+-2.04288e+11
+-1.95416e+11
+-1.86487e+11
+-1.77504e+11
+-1.6847e+11
+-1.59386e+11
+-1.50256e+11
+-1.41081e+11
+-1.31865e+11
+-1.2261e+11
+-1.13319e+11
+-1.03994e+11
+-9.46391e+10
+-8.52556e+10
+-7.58468e+10
+-6.64154e+10
+-5.69643e+10
+-4.74963e+10
+-3.8014e+10
+-2.85204e+10
+-1.90184e+10
+-9.5106e+09
+-0.00204791
diff --git a/examples/DDG4/data/Si220pl/Si220pl_efy.txt b/examples/DDG4/data/Si220pl/Si220pl_efy.txt
new file mode 100644
index 000000000..f3f304294
--- /dev/null
+++ b/examples/DDG4/data/Si220pl/Si220pl_efy.txt
@@ -0,0 +1,2050 @@
+2048 1 1
+1.91979e-10 2.21679e-10 3.13501e-10
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
+0
diff --git a/examples/DDG4/data/Si220pl/Si220pl_eld.txt b/examples/DDG4/data/Si220pl/Si220pl_eld.txt
new file mode 100644
index 000000000..918ad68be
--- /dev/null
+++ b/examples/DDG4/data/Si220pl/Si220pl_eld.txt
@@ -0,0 +1,2050 @@
+2048 1 1
+1.91979e-10 2.21679e-10 3.13501e-10
+4.1687
+4.16854
+4.16807
+4.16727
+4.16616
+4.16474
+4.16299
+4.16094
+4.15856
+4.15588
+4.15288
+4.14958
+4.14596
+4.14204
+4.13781
+4.13327
+4.12843
+4.1233
+4.11786
+4.11213
+4.10611
+4.09979
+4.09319
+4.0863
+4.07913
+4.07168
+4.06395
+4.05595
+4.04767
+4.03913
+4.03033
+4.02127
+4.01195
+4.00238
+3.99256
+3.9825
+3.97219
+3.96165
+3.95088
+3.93987
+3.92865
+3.9172
+3.90554
+3.89366
+3.88158
+3.8693
+3.85682
+3.84415
+3.83129
+3.81824
+3.80502
+3.79162
+3.77805
+3.76432
+3.75042
+3.73638
+3.72218
+3.70783
+3.69335
+3.67873
+3.66397
+3.64909
+3.63409
+3.61898
+3.60375
+3.58841
+3.57297
+3.55743
+3.54181
+3.52609
+3.51028
+3.4944
+3.47845
+3.46242
+3.44633
+3.43018
+3.41397
+3.3977
+3.38139
+3.36504
+3.34864
+3.33221
+3.31575
+3.29926
+3.28274
+3.26621
+3.24966
+3.2331
+3.21653
+3.19995
+3.18338
+3.1668
+3.15023
+3.13367
+3.11712
+3.10059
+3.08407
+3.06758
+3.05111
+3.03467
+3.01826
+3.00188
+2.98553
+2.96923
+2.95296
+2.93674
+2.92056
+2.90443
+2.88835
+2.87232
+2.85634
+2.84043
+2.82456
+2.80876
+2.79302
+2.77734
+2.76173
+2.74618
+2.7307
+2.71529
+2.69995
+2.68468
+2.66948
+2.65436
+2.63931
+2.62434
+2.60944
+2.59462
+2.57988
+2.56522
+2.55064
+2.53614
+2.52172
+2.50738
+2.49312
+2.47894
+2.46485
+2.45084
+2.43691
+2.42307
+2.40931
+2.39563
+2.38204
+2.36853
+2.3551
+2.34176
+2.3285
+2.31532
+2.30223
+2.28922
+2.27629
+2.26344
+2.25068
+2.238
+2.2254
+2.21288
+2.20044
+2.18809
+2.17581
+2.16361
+2.15149
+2.13945
+2.12748
+2.1156
+2.10379
+2.09206
+2.0804
+2.06882
+2.05731
+2.04587
+2.03451
+2.02323
+2.01201
+2.00086
+1.98979
+1.97879
+1.96785
+1.95698
+1.94619
+1.93545
+1.92479
+1.91419
+1.90366
+1.89319
+1.88279
+1.87245
+1.86217
+1.85195
+1.8418
+1.8317
+1.82167
+1.8117
+1.80178
+1.79192
+1.78212
+1.77238
+1.7627
+1.75307
+1.74349
+1.73397
+1.72451
+1.7151
+1.70574
+1.69643
+1.68718
+1.67797
+1.66882
+1.65972
+1.65067
+1.64166
+1.63271
+1.6238
+1.61495
+1.60614
+1.59738
+1.58866
+1.57999
+1.57137
+1.56279
+1.55426
+1.54577
+1.53733
+1.52893
+1.52057
+1.51226
+1.50399
+1.49576
+1.48758
+1.47943
+1.47133
+1.46327
+1.45526
+1.44728
+1.43934
+1.43144
+1.42359
+1.41577
+1.408
+1.40026
+1.39256
+1.3849
+1.37728
+1.3697
+1.36216
+1.35465
+1.34719
+1.33976
+1.33237
+1.32501
+1.3177
+1.31042
+1.30318
+1.29597
+1.28881
+1.28167
+1.27458
+1.26752
+1.2605
+1.25352
+1.24657
+1.23965
+1.23278
+1.22594
+1.21913
+1.21236
+1.20563
+1.19893
+1.19226
+1.18563
+1.17904
+1.17248
+1.16596
+1.15947
+1.15302
+1.1466
+1.14022
+1.13387
+1.12755
+1.12127
+1.11503
+1.10882
+1.10264
+1.0965
+1.09039
+1.08432
+1.07828
+1.07227
+1.0663
+1.06036
+1.05446
+1.04859
+1.04275
+1.03695
+1.03118
+1.02544
+1.01974
+1.01407
+1.00843
+1.00283
+0.997261
+0.991725
+0.986222
+0.980751
+0.975314
+0.969909
+0.964538
+0.959199
+0.953892
+0.948618
+0.943377
+0.938168
+0.932991
+0.927847
+0.922735
+0.917655
+0.912607
+0.90759
+0.902606
+0.897654
+0.892733
+0.887844
+0.882986
+0.87816
+0.873365
+0.868601
+0.863869
+0.859168
+0.854497
+0.849858
+0.845249
+0.840671
+0.836124
+0.831607
+0.82712
+0.822664
+0.818238
+0.813842
+0.809476
+0.80514
+0.800834
+0.796557
+0.792311
+0.788093
+0.783905
+0.779746
+0.775616
+0.771516
+0.767444
+0.763401
+0.759386
+0.755401
+0.751443
+0.747514
+0.743614
+0.739741
+0.735896
+0.732079
+0.72829
+0.724529
+0.720795
+0.717088
+0.713409
+0.709757
+0.706131
+0.702533
+0.698961
+0.695416
+0.691898
+0.688406
+0.68494
+0.6815
+0.678086
+0.674698
+0.671336
+0.667999
+0.664688
+0.661402
+0.658141
+0.654905
+0.651695
+0.648509
+0.645348
+0.642211
+0.639099
+0.63601
+0.632947
+0.629907
+0.626891
+0.623898
+0.62093
+0.617984
+0.615062
+0.612164
+0.609288
+0.606435
+0.603605
+0.600798
+0.598013
+0.595251
+0.592511
+0.589793
+0.587097
+0.584423
+0.581771
+0.57914
+0.57653
+0.573942
+0.571376
+0.56883
+0.566305
+0.563801
+0.561317
+0.558854
+0.556412
+0.55399
+0.551587
+0.549205
+0.546843
+0.5445
+0.542177
+0.539874
+0.53759
+0.535325
+0.533079
+0.530852
+0.528644
+0.526454
+0.524283
+0.522131
+0.519997
+0.517881
+0.515783
+0.513703
+0.511641
+0.509596
+0.507569
+0.50556
+0.503567
+0.501592
+0.499634
+0.497693
+0.495769
+0.493861
+0.49197
+0.490096
+0.488237
+0.486395
+0.484569
+0.482759
+0.480965
+0.479186
+0.477423
+0.475676
+0.473944
+0.472227
+0.470525
+0.468839
+0.467167
+0.46551
+0.463868
+0.46224
+0.460627
+0.459028
+0.457443
+0.455873
+0.454316
+0.452773
+0.451245
+0.449729
+0.448228
+0.44674
+0.445265
+0.443803
+0.442355
+0.440919
+0.439497
+0.438087
+0.43669
+0.435306
+0.433934
+0.432575
+0.431228
+0.429893
+0.428571
+0.42726
+0.425961
+0.424674
+0.423399
+0.422135
+0.420883
+0.419643
+0.418413
+0.417195
+0.415988
+0.414792
+0.413607
+0.412433
+0.41127
+0.410117
+0.408975
+0.407844
+0.406723
+0.405612
+0.404511
+0.403421
+0.40234
+0.40127
+0.40021
+0.399159
+0.398118
+0.397086
+0.396064
+0.395052
+0.394049
+0.393055
+0.392071
+0.391095
+0.390129
+0.389172
+0.388223
+0.387283
+0.386352
+0.38543
+0.384516
+0.383611
+0.382714
+0.381826
+0.380945
+0.380073
+0.379209
+0.378354
+0.377506
+0.376666
+0.375833
+0.375009
+0.374192
+0.373383
+0.372582
+0.371787
+0.371001
+0.370221
+0.369449
+0.368684
+0.367927
+0.367176
+0.366432
+0.365695
+0.364965
+0.364242
+0.363526
+0.362816
+0.362113
+0.361417
+0.360727
+0.360043
+0.359366
+0.358695
+0.35803
+0.357372
+0.356719
+0.356073
+0.355433
+0.354799
+0.35417
+0.353547
+0.352931
+0.35232
+0.351714
+0.351114
+0.35052
+0.349931
+0.349348
+0.34877
+0.348198
+0.34763
+0.347068
+0.346512
+0.34596
+0.345413
+0.344872
+0.344335
+0.343804
+0.343277
+0.342755
+0.342238
+0.341725
+0.341218
+0.340715
+0.340216
+0.339723
+0.339233
+0.338748
+0.338268
+0.337792
+0.33732
+0.336853
+0.33639
+0.335931
+0.335476
+0.335025
+0.334579
+0.334136
+0.333697
+0.333263
+0.332832
+0.332405
+0.331982
+0.331563
+0.331148
+0.330736
+0.330328
+0.329923
+0.329523
+0.329125
+0.328732
+0.328341
+0.327955
+0.327571
+0.327191
+0.326815
+0.326441
+0.326071
+0.325705
+0.325341
+0.324981
+0.324624
+0.324269
+0.323918
+0.32357
+0.323226
+0.322884
+0.322545
+0.322209
+0.321875
+0.321545
+0.321218
+0.320893
+0.320571
+0.320252
+0.319935
+0.319622
+0.31931
+0.319002
+0.318696
+0.318393
+0.318092
+0.317794
+0.317498
+0.317205
+0.316914
+0.316626
+0.31634
+0.316056
+0.315775
+0.315496
+0.315219
+0.314945
+0.314672
+0.314402
+0.314135
+0.313869
+0.313606
+0.313344
+0.313085
+0.312828
+0.312573
+0.31232
+0.312069
+0.31182
+0.311573
+0.311328
+0.311084
+0.310843
+0.310604
+0.310366
+0.310131
+0.309897
+0.309665
+0.309435
+0.309206
+0.30898
+0.308755
+0.308532
+0.30831
+0.30809
+0.307872
+0.307656
+0.307441
+0.307228
+0.307016
+0.306806
+0.306597
+0.30639
+0.306185
+0.305981
+0.305779
+0.305578
+0.305378
+0.30518
+0.304984
+0.304789
+0.304595
+0.304403
+0.304212
+0.304022
+0.303834
+0.303647
+0.303462
+0.303277
+0.303095
+0.302913
+0.302733
+0.302553
+0.302376
+0.302199
+0.302024
+0.301849
+0.301676
+0.301505
+0.301334
+0.301165
+0.300996
+0.300829
+0.300663
+0.300498
+0.300334
+0.300172
+0.30001
+0.29985
+0.29969
+0.299532
+0.299374
+0.299218
+0.299063
+0.298908
+0.298755
+0.298603
+0.298451
+0.298301
+0.298152
+0.298003
+0.297856
+0.29771
+0.297564
+0.297419
+0.297276
+0.297133
+0.296991
+0.29685
+0.29671
+0.296571
+0.296432
+0.296295
+0.296158
+0.296022
+0.295888
+0.295753
+0.29562
+0.295488
+0.295356
+0.295225
+0.295095
+0.294966
+0.294838
+0.29471
+0.294583
+0.294457
+0.294331
+0.294207
+0.294083
+0.29396
+0.293838
+0.293716
+0.293595
+0.293475
+0.293355
+0.293237
+0.293119
+0.293001
+0.292885
+0.292769
+0.292653
+0.292539
+0.292425
+0.292312
+0.292199
+0.292087
+0.291976
+0.291865
+0.291755
+0.291646
+0.291537
+0.291429
+0.291322
+0.291215
+0.291109
+0.291003
+0.290899
+0.290794
+0.290691
+0.290587
+0.290485
+0.290383
+0.290282
+0.290181
+0.290081
+0.289982
+0.289883
+0.289784
+0.289687
+0.289589
+0.289493
+0.289397
+0.289301
+0.289206
+0.289112
+0.289018
+0.288925
+0.288832
+0.28874
+0.288648
+0.288557
+0.288466
+0.288376
+0.288287
+0.288198
+0.28811
+0.288022
+0.287934
+0.287847
+0.287761
+0.287675
+0.28759
+0.287505
+0.287421
+0.287337
+0.287254
+0.287171
+0.287089
+0.287007
+0.286926
+0.286845
+0.286765
+0.286685
+0.286606
+0.286527
+0.286449
+0.286371
+0.286294
+0.286217
+0.286141
+0.286065
+0.28599
+0.285915
+0.285841
+0.285767
+0.285693
+0.28562
+0.285548
+0.285476
+0.285404
+0.285333
+0.285263
+0.285193
+0.285123
+0.285054
+0.284985
+0.284917
+0.284849
+0.284782
+0.284715
+0.284648
+0.284582
+0.284517
+0.284452
+0.284387
+0.284323
+0.28426
+0.284196
+0.284134
+0.284071
+0.284009
+0.283948
+0.283887
+0.283827
+0.283767
+0.283707
+0.283648
+0.283589
+0.283531
+0.283473
+0.283416
+0.283359
+0.283302
+0.283246
+0.283191
+0.283136
+0.283081
+0.283027
+0.282973
+0.282919
+0.282866
+0.282814
+0.282762
+0.28271
+0.282659
+0.282608
+0.282558
+0.282508
+0.282459
+0.28241
+0.282361
+0.282313
+0.282265
+0.282218
+0.282171
+0.282125
+0.282079
+0.282034
+0.281988
+0.281944
+0.2819
+0.281856
+0.281812
+0.28177
+0.281727
+0.281685
+0.281643
+0.281602
+0.281562
+0.281521
+0.281481
+0.281442
+0.281403
+0.281364
+0.281326
+0.281288
+0.281251
+0.281214
+0.281178
+0.281142
+0.281106
+0.281071
+0.281036
+0.281002
+0.280968
+0.280935
+0.280902
+0.280869
+0.280837
+0.280806
+0.280774
+0.280743
+0.280713
+0.280683
+0.280653
+0.280624
+0.280596
+0.280567
+0.28054
+0.280512
+0.280485
+0.280459
+0.280432
+0.280407
+0.280381
+0.280357
+0.280332
+0.280308
+0.280285
+0.280262
+0.280239
+0.280217
+0.280195
+0.280173
+0.280152
+0.280132
+0.280112
+0.280092
+0.280073
+0.280054
+0.280035
+0.280017
+0.28
+0.279983
+0.279966
+0.27995
+0.279934
+0.279919
+0.279904
+0.279889
+0.279875
+0.279861
+0.279848
+0.279835
+0.279823
+0.279811
+0.279799
+0.279788
+0.279777
+0.279767
+0.279757
+0.279748
+0.279739
+0.27973
+0.279722
+0.279714
+0.279707
+0.2797
+0.279694
+0.279688
+0.279682
+0.279677
+0.279672
+0.279668
+0.279664
+0.279661
+0.279658
+0.279655
+0.279653
+0.279651
+0.27965
+0.279649
+0.279649
+0.279649
+0.279649
+0.27965
+0.279651
+0.279653
+0.279655
+0.279658
+0.279661
+0.279664
+0.279668
+0.279672
+0.279677
+0.279682
+0.279688
+0.279694
+0.2797
+0.279707
+0.279714
+0.279722
+0.27973
+0.279739
+0.279748
+0.279757
+0.279767
+0.279777
+0.279788
+0.279799
+0.279811
+0.279823
+0.279835
+0.279848
+0.279861
+0.279875
+0.279889
+0.279904
+0.279919
+0.279934
+0.27995
+0.279966
+0.279983
+0.28
+0.280017
+0.280035
+0.280054
+0.280073
+0.280092
+0.280112
+0.280132
+0.280152
+0.280173
+0.280195
+0.280217
+0.280239
+0.280262
+0.280285
+0.280308
+0.280332
+0.280357
+0.280381
+0.280407
+0.280432
+0.280459
+0.280485
+0.280512
+0.28054
+0.280567
+0.280596
+0.280624
+0.280653
+0.280683
+0.280713
+0.280743
+0.280774
+0.280806
+0.280837
+0.280869
+0.280902
+0.280935
+0.280968
+0.281002
+0.281036
+0.281071
+0.281106
+0.281142
+0.281178
+0.281214
+0.281251
+0.281288
+0.281326
+0.281364
+0.281403
+0.281442
+0.281481
+0.281521
+0.281562
+0.281602
+0.281643
+0.281685
+0.281727
+0.28177
+0.281812
+0.281856
+0.2819
+0.281944
+0.281988
+0.282034
+0.282079
+0.282125
+0.282171
+0.282218
+0.282265
+0.282313
+0.282361
+0.28241
+0.282459
+0.282508
+0.282558
+0.282608
+0.282659
+0.28271
+0.282762
+0.282814
+0.282866
+0.282919
+0.282973
+0.283027
+0.283081
+0.283136
+0.283191
+0.283246
+0.283302
+0.283359
+0.283416
+0.283473
+0.283531
+0.283589
+0.283648
+0.283707
+0.283767
+0.283827
+0.283887
+0.283948
+0.284009
+0.284071
+0.284134
+0.284196
+0.28426
+0.284323
+0.284387
+0.284452
+0.284517
+0.284582
+0.284648
+0.284715
+0.284782
+0.284849
+0.284917
+0.284985
+0.285054
+0.285123
+0.285193
+0.285263
+0.285333
+0.285404
+0.285476
+0.285548
+0.28562
+0.285693
+0.285767
+0.285841
+0.285915
+0.28599
+0.286065
+0.286141
+0.286217
+0.286294
+0.286371
+0.286449
+0.286527
+0.286606
+0.286685
+0.286765
+0.286845
+0.286926
+0.287007
+0.287089
+0.287171
+0.287254
+0.287337
+0.287421
+0.287505
+0.28759
+0.287675
+0.287761
+0.287847
+0.287934
+0.288022
+0.28811
+0.288198
+0.288287
+0.288376
+0.288466
+0.288557
+0.288648
+0.28874
+0.288832
+0.288925
+0.289018
+0.289112
+0.289206
+0.289301
+0.289397
+0.289493
+0.289589
+0.289687
+0.289784
+0.289883
+0.289982
+0.290081
+0.290181
+0.290282
+0.290383
+0.290485
+0.290587
+0.290691
+0.290794
+0.290899
+0.291003
+0.291109
+0.291215
+0.291322
+0.291429
+0.291537
+0.291646
+0.291755
+0.291865
+0.291976
+0.292087
+0.292199
+0.292312
+0.292425
+0.292539
+0.292653
+0.292769
+0.292885
+0.293001
+0.293119
+0.293237
+0.293355
+0.293475
+0.293595
+0.293716
+0.293838
+0.29396
+0.294083
+0.294207
+0.294331
+0.294457
+0.294583
+0.29471
+0.294838
+0.294966
+0.295095
+0.295225
+0.295356
+0.295488
+0.29562
+0.295753
+0.295888
+0.296022
+0.296158
+0.296295
+0.296432
+0.296571
+0.29671
+0.29685
+0.296991
+0.297133
+0.297276
+0.297419
+0.297564
+0.29771
+0.297856
+0.298003
+0.298152
+0.298301
+0.298451
+0.298603
+0.298755
+0.298908
+0.299063
+0.299218
+0.299374
+0.299532
+0.29969
+0.29985
+0.30001
+0.300172
+0.300334
+0.300498
+0.300663
+0.300829
+0.300996
+0.301165
+0.301334
+0.301505
+0.301676
+0.301849
+0.302024
+0.302199
+0.302376
+0.302553
+0.302733
+0.302913
+0.303095
+0.303277
+0.303462
+0.303647
+0.303834
+0.304022
+0.304212
+0.304403
+0.304595
+0.304789
+0.304984
+0.30518
+0.305378
+0.305578
+0.305779
+0.305981
+0.306185
+0.30639
+0.306597
+0.306806
+0.307016
+0.307228
+0.307441
+0.307656
+0.307872
+0.30809
+0.30831
+0.308532
+0.308755
+0.30898
+0.309206
+0.309435
+0.309665
+0.309897
+0.310131
+0.310366
+0.310604
+0.310843
+0.311084
+0.311328
+0.311573
+0.31182
+0.312069
+0.31232
+0.312573
+0.312828
+0.313085
+0.313344
+0.313606
+0.313869
+0.314135
+0.314402
+0.314672
+0.314945
+0.315219
+0.315496
+0.315775
+0.316056
+0.31634
+0.316626
+0.316914
+0.317205
+0.317498
+0.317794
+0.318092
+0.318393
+0.318696
+0.319002
+0.31931
+0.319622
+0.319935
+0.320252
+0.320571
+0.320893
+0.321218
+0.321545
+0.321875
+0.322209
+0.322545
+0.322884
+0.323226
+0.32357
+0.323918
+0.324269
+0.324624
+0.324981
+0.325341
+0.325705
+0.326071
+0.326441
+0.326815
+0.327191
+0.327571
+0.327955
+0.328341
+0.328732
+0.329125
+0.329523
+0.329923
+0.330328
+0.330736
+0.331148
+0.331563
+0.331982
+0.332405
+0.332832
+0.333263
+0.333697
+0.334136
+0.334579
+0.335025
+0.335476
+0.335931
+0.33639
+0.336853
+0.33732
+0.337792
+0.338268
+0.338748
+0.339233
+0.339723
+0.340216
+0.340715
+0.341218
+0.341725
+0.342238
+0.342755
+0.343277
+0.343804
+0.344335
+0.344872
+0.345413
+0.34596
+0.346512
+0.347068
+0.34763
+0.348198
+0.34877
+0.349348
+0.349931
+0.35052
+0.351114
+0.351714
+0.35232
+0.352931
+0.353547
+0.35417
+0.354799
+0.355433
+0.356073
+0.356719
+0.357372
+0.35803
+0.358695
+0.359366
+0.360043
+0.360727
+0.361417
+0.362113
+0.362816
+0.363526
+0.364242
+0.364965
+0.365695
+0.366432
+0.367176
+0.367927
+0.368684
+0.369449
+0.370221
+0.371001
+0.371787
+0.372582
+0.373383
+0.374192
+0.375009
+0.375833
+0.376666
+0.377506
+0.378354
+0.379209
+0.380073
+0.380945
+0.381826
+0.382714
+0.383611
+0.384516
+0.38543
+0.386352
+0.387283
+0.388223
+0.389172
+0.390129
+0.391095
+0.392071
+0.393055
+0.394049
+0.395052
+0.396064
+0.397086
+0.398118
+0.399159
+0.40021
+0.40127
+0.40234
+0.403421
+0.404511
+0.405612
+0.406723
+0.407844
+0.408975
+0.410117
+0.41127
+0.412433
+0.413607
+0.414792
+0.415988
+0.417195
+0.418413
+0.419643
+0.420883
+0.422135
+0.423399
+0.424674
+0.425961
+0.42726
+0.428571
+0.429893
+0.431228
+0.432575
+0.433934
+0.435306
+0.43669
+0.438087
+0.439497
+0.440919
+0.442355
+0.443803
+0.445265
+0.44674
+0.448228
+0.449729
+0.451245
+0.452773
+0.454316
+0.455873
+0.457443
+0.459028
+0.460627
+0.46224
+0.463868
+0.46551
+0.467167
+0.468839
+0.470525
+0.472227
+0.473944
+0.475676
+0.477423
+0.479186
+0.480965
+0.482759
+0.484569
+0.486395
+0.488237
+0.490096
+0.49197
+0.493861
+0.495769
+0.497693
+0.499634
+0.501592
+0.503567
+0.50556
+0.507569
+0.509596
+0.511641
+0.513703
+0.515783
+0.517881
+0.519997
+0.522131
+0.524283
+0.526454
+0.528644
+0.530852
+0.533079
+0.535325
+0.53759
+0.539874
+0.542177
+0.5445
+0.546843
+0.549205
+0.551587
+0.55399
+0.556412
+0.558854
+0.561317
+0.563801
+0.566305
+0.56883
+0.571376
+0.573942
+0.57653
+0.57914
+0.581771
+0.584423
+0.587097
+0.589793
+0.592511
+0.595251
+0.598013
+0.600798
+0.603605
+0.606435
+0.609288
+0.612164
+0.615062
+0.617984
+0.62093
+0.623898
+0.626891
+0.629907
+0.632947
+0.63601
+0.639099
+0.642211
+0.645348
+0.648509
+0.651695
+0.654905
+0.658141
+0.661402
+0.664688
+0.667999
+0.671336
+0.674698
+0.678086
+0.6815
+0.68494
+0.688406
+0.691898
+0.695416
+0.698961
+0.702533
+0.706131
+0.709757
+0.713409
+0.717088
+0.720795
+0.724529
+0.72829
+0.732079
+0.735896
+0.739741
+0.743614
+0.747514
+0.751443
+0.755401
+0.759386
+0.763401
+0.767444
+0.771516
+0.775616
+0.779746
+0.783905
+0.788093
+0.792311
+0.796557
+0.800834
+0.80514
+0.809476
+0.813842
+0.818238
+0.822664
+0.82712
+0.831607
+0.836124
+0.840671
+0.845249
+0.849858
+0.854497
+0.859168
+0.863869
+0.868601
+0.873365
+0.87816
+0.882986
+0.887844
+0.892733
+0.897654
+0.902606
+0.90759
+0.912607
+0.917655
+0.922735
+0.927847
+0.932991
+0.938168
+0.943377
+0.948618
+0.953892
+0.959199
+0.964538
+0.969909
+0.975314
+0.980751
+0.986222
+0.991725
+0.997261
+1.00283
+1.00843
+1.01407
+1.01974
+1.02544
+1.03118
+1.03695
+1.04275
+1.04859
+1.05446
+1.06036
+1.0663
+1.07227
+1.07828
+1.08432
+1.09039
+1.0965
+1.10264
+1.10882
+1.11503
+1.12127
+1.12755
+1.13387
+1.14022
+1.1466
+1.15302
+1.15947
+1.16596
+1.17248
+1.17904
+1.18563
+1.19226
+1.19893
+1.20563
+1.21236
+1.21913
+1.22594
+1.23278
+1.23965
+1.24657
+1.25352
+1.2605
+1.26752
+1.27458
+1.28167
+1.28881
+1.29597
+1.30318
+1.31042
+1.3177
+1.32501
+1.33237
+1.33976
+1.34719
+1.35465
+1.36216
+1.3697
+1.37728
+1.3849
+1.39256
+1.40026
+1.408
+1.41577
+1.42359
+1.43144
+1.43934
+1.44728
+1.45526
+1.46327
+1.47133
+1.47943
+1.48758
+1.49576
+1.50399
+1.51226
+1.52057
+1.52893
+1.53733
+1.54577
+1.55426
+1.56279
+1.57137
+1.57999
+1.58866
+1.59738
+1.60614
+1.61495
+1.6238
+1.63271
+1.64166
+1.65067
+1.65972
+1.66882
+1.67797
+1.68718
+1.69643
+1.70574
+1.7151
+1.72451
+1.73397
+1.74349
+1.75307
+1.7627
+1.77238
+1.78212
+1.79192
+1.80178
+1.8117
+1.82167
+1.8317
+1.8418
+1.85195
+1.86217
+1.87245
+1.88279
+1.89319
+1.90366
+1.91419
+1.92479
+1.93545
+1.94619
+1.95698
+1.96785
+1.97879
+1.98979
+2.00086
+2.01201
+2.02323
+2.03451
+2.04587
+2.05731
+2.06882
+2.0804
+2.09206
+2.10379
+2.1156
+2.12748
+2.13945
+2.15149
+2.16361
+2.17581
+2.18809
+2.20044
+2.21288
+2.2254
+2.238
+2.25068
+2.26344
+2.27629
+2.28922
+2.30223
+2.31532
+2.3285
+2.34176
+2.3551
+2.36853
+2.38204
+2.39563
+2.40931
+2.42307
+2.43691
+2.45084
+2.46485
+2.47894
+2.49312
+2.50738
+2.52172
+2.53614
+2.55064
+2.56522
+2.57988
+2.59462
+2.60944
+2.62434
+2.63931
+2.65436
+2.66948
+2.68468
+2.69995
+2.71529
+2.7307
+2.74618
+2.76173
+2.77734
+2.79302
+2.80876
+2.82456
+2.84043
+2.85634
+2.87232
+2.88835
+2.90443
+2.92056
+2.93674
+2.95296
+2.96923
+2.98553
+3.00188
+3.01826
+3.03467
+3.05111
+3.06758
+3.08407
+3.10059
+3.11712
+3.13367
+3.15023
+3.1668
+3.18338
+3.19995
+3.21653
+3.2331
+3.24966
+3.26621
+3.28274
+3.29926
+3.31575
+3.33221
+3.34864
+3.36504
+3.38139
+3.3977
+3.41397
+3.43018
+3.44633
+3.46242
+3.47845
+3.4944
+3.51028
+3.52609
+3.54181
+3.55743
+3.57297
+3.58841
+3.60375
+3.61898
+3.63409
+3.64909
+3.66397
+3.67873
+3.69335
+3.70783
+3.72218
+3.73638
+3.75042
+3.76432
+3.77805
+3.79162
+3.80502
+3.81824
+3.83129
+3.84415
+3.85682
+3.8693
+3.88158
+3.89366
+3.90554
+3.9172
+3.92865
+3.93987
+3.95088
+3.96165
+3.97219
+3.9825
+3.99256
+4.00238
+4.01195
+4.02127
+4.03033
+4.03913
+4.04767
+4.05595
+4.06395
+4.07168
+4.07913
+4.0863
+4.09319
+4.09979
+4.10611
+4.11213
+4.11786
+4.1233
+4.12843
+4.13327
+4.13781
+4.14204
+4.14596
+4.14958
+4.15288
+4.15588
+4.15856
+4.16094
+4.16299
+4.16474
+4.16616
+4.16727
+4.16807
+4.16854
+4.1687
diff --git a/examples/DDG4/data/Si220pl/Si220pl_pot.txt b/examples/DDG4/data/Si220pl/Si220pl_pot.txt
new file mode 100644
index 000000000..474e01151
--- /dev/null
+++ b/examples/DDG4/data/Si220pl/Si220pl_pot.txt
@@ -0,0 +1,2050 @@
+2048 1 1
+1.91979e-10 2.21679e-10 3.13501e-10
+15.3567
+15.3563
+15.355
+15.3527
+15.3496
+15.3456
+15.3407
+15.3349
+15.3282
+15.3207
+15.3123
+15.3029
+15.2928
+15.2817
+15.2698
+15.257
+15.2433
+15.2288
+15.2134
+15.1972
+15.1801
+15.1622
+15.1435
+15.1239
+15.1035
+15.0823
+15.0602
+15.0374
+15.0138
+14.9894
+14.9641
+14.9382
+14.9114
+14.8839
+14.8556
+14.8266
+14.7968
+14.7663
+14.7351
+14.7032
+14.6705
+14.6372
+14.6032
+14.5685
+14.5331
+14.497
+14.4603
+14.423
+14.385
+14.3465
+14.3073
+14.2675
+14.2271
+14.1861
+14.1445
+14.1024
+14.0597
+14.0165
+13.9728
+13.9285
+13.8837
+13.8385
+13.7927
+13.7464
+13.6997
+13.6525
+13.6049
+13.5569
+13.5084
+13.4595
+13.4102
+13.3604
+13.3104
+13.2599
+13.2091
+13.1579
+13.1063
+13.0545
+13.0023
+12.9498
+12.897
+12.8439
+12.7905
+12.7369
+12.683
+12.6288
+12.5744
+12.5198
+12.4649
+12.4098
+12.3545
+12.299
+12.2434
+12.1875
+12.1315
+12.0753
+12.019
+11.9625
+11.9059
+11.8492
+11.7923
+11.7354
+11.6783
+11.6212
+11.5639
+11.5066
+11.4492
+11.3917
+11.3342
+11.2767
+11.2191
+11.1614
+11.1038
+11.0461
+10.9884
+10.9307
+10.8729
+10.8152
+10.7575
+10.6999
+10.6422
+10.5846
+10.527
+10.4694
+10.4119
+10.3544
+10.297
+10.2396
+10.1823
+10.1251
+10.0679
+10.0109
+9.95386
+9.89695
+9.84012
+9.78339
+9.72676
+9.67023
+9.61381
+9.55749
+9.50128
+9.44519
+9.38922
+9.33337
+9.27764
+9.22204
+9.16657
+9.11123
+9.05603
+9.00096
+8.94604
+8.89126
+8.83662
+8.78213
+8.72779
+8.67359
+8.61956
+8.56567
+8.51194
+8.45837
+8.40496
+8.35171
+8.29863
+8.2457
+8.19294
+8.14035
+8.08792
+8.03567
+7.98358
+7.93166
+7.87992
+7.82834
+7.77694
+7.72571
+7.67466
+7.62378
+7.57308
+7.52255
+7.4722
+7.42202
+7.37203
+7.32221
+7.27256
+7.2231
+7.17381
+7.1247
+7.07577
+7.02702
+6.97845
+6.93005
+6.88184
+6.8338
+6.78594
+6.73826
+6.69076
+6.64343
+6.59628
+6.54931
+6.50252
+6.4559
+6.40946
+6.3632
+6.31711
+6.2712
+6.22546
+6.17989
+6.13451
+6.08929
+6.04425
+5.99938
+5.95468
+5.91016
+5.8658
+5.82162
+5.77761
+5.73376
+5.69009
+5.64658
+5.60324
+5.56007
+5.51706
+5.47422
+5.43155
+5.38904
+5.3467
+5.30451
+5.26249
+5.22064
+5.17894
+5.1374
+5.09603
+5.05481
+5.01375
+4.97285
+4.93211
+4.89152
+4.85109
+4.81081
+4.77069
+4.73072
+4.6909
+4.65124
+4.61172
+4.57236
+4.53315
+4.49408
+4.45517
+4.4164
+4.37778
+4.3393
+4.30097
+4.26278
+4.22474
+4.18684
+4.14909
+4.11147
+4.074
+4.03667
+3.99947
+3.96242
+3.9255
+3.88872
+3.85208
+3.81557
+3.7792
+3.74296
+3.70685
+3.67088
+3.63504
+3.59934
+3.56376
+3.52831
+3.493
+3.45781
+3.42275
+3.38781
+3.35301
+3.31833
+3.28377
+3.24934
+3.21504
+3.18085
+3.14679
+3.11285
+3.07903
+3.04534
+3.01176
+2.9783
+2.94496
+2.91174
+2.87864
+2.84565
+2.81278
+2.78002
+2.74738
+2.71485
+2.68244
+2.65014
+2.61795
+2.58587
+2.5539
+2.52205
+2.4903
+2.45867
+2.42714
+2.39572
+2.3644
+2.3332
+2.3021
+2.2711
+2.24021
+2.20943
+2.17874
+2.14817
+2.11769
+2.08732
+2.05705
+2.02687
+1.9968
+1.96683
+1.93696
+1.90719
+1.87752
+1.84794
+1.81846
+1.78908
+1.75979
+1.7306
+1.70151
+1.67251
+1.6436
+1.61479
+1.58607
+1.55744
+1.52891
+1.50046
+1.47211
+1.44385
+1.41568
+1.3876
+1.3596
+1.3317
+1.30388
+1.27615
+1.24851
+1.22096
+1.19349
+1.16611
+1.13881
+1.1116
+1.08447
+1.05743
+1.03047
+1.0036
+0.976803
+0.950092
+0.923464
+0.896916
+0.87045
+0.844065
+0.817759
+0.791534
+0.765388
+0.739321
+0.713332
+0.687422
+0.66159
+0.635835
+0.610157
+0.584556
+0.559031
+0.533582
+0.508208
+0.48291
+0.457686
+0.432537
+0.407461
+0.38246
+0.357531
+0.332676
+0.307893
+0.283182
+0.258543
+0.233976
+0.20948
+0.185054
+0.160699
+0.136414
+0.112199
+0.0880539
+0.0639775
+0.0399699
+0.0160308
+-0.00784012
+-0.0316432
+-0.0553788
+-0.0790472
+-0.102649
+-0.126184
+-0.149652
+-0.173055
+-0.196392
+-0.219664
+-0.24287
+-0.266012
+-0.289089
+-0.312102
+-0.335051
+-0.357937
+-0.380759
+-0.403519
+-0.426215
+-0.448849
+-0.471421
+-0.49393
+-0.516378
+-0.538765
+-0.561091
+-0.583355
+-0.605559
+-0.627703
+-0.649786
+-0.67181
+-0.693774
+-0.715679
+-0.737524
+-0.759311
+-0.781039
+-0.802709
+-0.824321
+-0.845874
+-0.86737
+-0.888809
+-0.91019
+-0.931515
+-0.952783
+-0.973994
+-0.995149
+-1.01625
+-1.03729
+-1.05828
+-1.07921
+-1.10009
+-1.12091
+-1.14168
+-1.16239
+-1.18305
+-1.20365
+-1.2242
+-1.2447
+-1.26514
+-1.28553
+-1.30586
+-1.32615
+-1.34638
+-1.36656
+-1.38668
+-1.40675
+-1.42678
+-1.44675
+-1.46666
+-1.48653
+-1.50634
+-1.52611
+-1.54582
+-1.56548
+-1.58509
+-1.60466
+-1.62417
+-1.64363
+-1.66304
+-1.6824
+-1.70171
+-1.72097
+-1.74019
+-1.75935
+-1.77847
+-1.79753
+-1.81655
+-1.83552
+-1.85444
+-1.87331
+-1.89214
+-1.91091
+-1.92964
+-1.94832
+-1.96696
+-1.98555
+-2.00408
+-2.02258
+-2.04102
+-2.05942
+-2.07777
+-2.09608
+-2.11434
+-2.13255
+-2.15072
+-2.16884
+-2.18692
+-2.20495
+-2.22293
+-2.24087
+-2.25877
+-2.27661
+-2.29442
+-2.31218
+-2.32989
+-2.34756
+-2.36518
+-2.38277
+-2.4003
+-2.41779
+-2.43524
+-2.45264
+-2.47
+-2.48732
+-2.50459
+-2.52182
+-2.53901
+-2.55615
+-2.57325
+-2.5903
+-2.60732
+-2.62429
+-2.64121
+-2.6581
+-2.67494
+-2.69174
+-2.70849
+-2.72521
+-2.74188
+-2.75851
+-2.7751
+-2.79164
+-2.80815
+-2.82461
+-2.84103
+-2.85741
+-2.87375
+-2.89004
+-2.9063
+-2.92251
+-2.93869
+-2.95482
+-2.97091
+-2.98696
+-3.00297
+-3.01894
+-3.03486
+-3.05075
+-3.0666
+-3.0824
+-3.09817
+-3.11389
+-3.12958
+-3.14523
+-3.16083
+-3.1764
+-3.19192
+-3.20741
+-3.22285
+-3.23826
+-3.25363
+-3.26895
+-3.28424
+-3.29949
+-3.3147
+-3.32987
+-3.345
+-3.36009
+-3.37515
+-3.39016
+-3.40513
+-3.42007
+-3.43497
+-3.44983
+-3.46465
+-3.47943
+-3.49417
+-3.50887
+-3.52354
+-3.53817
+-3.55276
+-3.56731
+-3.58182
+-3.5963
+-3.61073
+-3.62513
+-3.63949
+-3.65382
+-3.6681
+-3.68235
+-3.69656
+-3.71073
+-3.72486
+-3.73896
+-3.75302
+-3.76704
+-3.78103
+-3.79497
+-3.80888
+-3.82276
+-3.83659
+-3.85039
+-3.86415
+-3.87788
+-3.89156
+-3.90521
+-3.91883
+-3.9324
+-3.94594
+-3.95944
+-3.97291
+-3.98634
+-3.99973
+-4.01309
+-4.02641
+-4.03969
+-4.05294
+-4.06615
+-4.07932
+-4.09246
+-4.10556
+-4.11862
+-4.13165
+-4.14465
+-4.1576
+-4.17052
+-4.18341
+-4.19625
+-4.20907
+-4.22184
+-4.23458
+-4.24729
+-4.25995
+-4.27259
+-4.28518
+-4.29775
+-4.31027
+-4.32276
+-4.33521
+-4.34763
+-4.36002
+-4.37236
+-4.38468
+-4.39695
+-4.40919
+-4.4214
+-4.43357
+-4.44571
+-4.4578
+-4.46987
+-4.4819
+-4.49389
+-4.50585
+-4.51778
+-4.52966
+-4.54152
+-4.55334
+-4.56512
+-4.57687
+-4.58858
+-4.60026
+-4.61191
+-4.62352
+-4.63509
+-4.64663
+-4.65813
+-4.6696
+-4.68104
+-4.69244
+-4.70381
+-4.71514
+-4.72643
+-4.7377
+-4.74892
+-4.76012
+-4.77127
+-4.7824
+-4.79349
+-4.80454
+-4.81556
+-4.82655
+-4.8375
+-4.84842
+-4.8593
+-4.87015
+-4.88097
+-4.89175
+-4.90249
+-4.91321
+-4.92388
+-4.93453
+-4.94514
+-4.95571
+-4.96625
+-4.97676
+-4.98724
+-4.99768
+-5.00808
+-5.01845
+-5.02879
+-5.03909
+-5.04936
+-5.0596
+-5.0698
+-5.07997
+-5.0901
+-5.10021
+-5.11027
+-5.12031
+-5.1303
+-5.14027
+-5.1502
+-5.1601
+-5.16997
+-5.1798
+-5.18959
+-5.19936
+-5.20909
+-5.21879
+-5.22845
+-5.23808
+-5.24768
+-5.25724
+-5.26677
+-5.27626
+-5.28573
+-5.29515
+-5.30455
+-5.31391
+-5.32324
+-5.33254
+-5.3418
+-5.35103
+-5.36022
+-5.36939
+-5.37852
+-5.38761
+-5.39668
+-5.4057
+-5.4147
+-5.42366
+-5.43259
+-5.44149
+-5.45036
+-5.45919
+-5.46798
+-5.47675
+-5.48548
+-5.49418
+-5.50284
+-5.51148
+-5.52008
+-5.52864
+-5.53718
+-5.54568
+-5.55414
+-5.56258
+-5.57098
+-5.57935
+-5.58769
+-5.59599
+-5.60426
+-5.6125
+-5.6207
+-5.62887
+-5.63701
+-5.64512
+-5.65319
+-5.66124
+-5.66924
+-5.67722
+-5.68516
+-5.69307
+-5.70095
+-5.70879
+-5.71661
+-5.72439
+-5.73213
+-5.73985
+-5.74753
+-5.75518
+-5.76279
+-5.77038
+-5.77793
+-5.78545
+-5.79293
+-5.80039
+-5.80781
+-5.8152
+-5.82255
+-5.82988
+-5.83717
+-5.84443
+-5.85166
+-5.85885
+-5.86601
+-5.87314
+-5.88024
+-5.8873
+-5.89434
+-5.90134
+-5.9083
+-5.91524
+-5.92214
+-5.92901
+-5.93585
+-5.94266
+-5.94943
+-5.95617
+-5.96288
+-5.96956
+-5.9762
+-5.98282
+-5.9894
+-5.99595
+-6.00246
+-6.00895
+-6.0154
+-6.02182
+-6.0282
+-6.03456
+-6.04088
+-6.04717
+-6.05343
+-6.05966
+-6.06586
+-6.07202
+-6.07815
+-6.08425
+-6.09031
+-6.09635
+-6.10235
+-6.10832
+-6.11426
+-6.12017
+-6.12604
+-6.13188
+-6.1377
+-6.14347
+-6.14922
+-6.15494
+-6.16062
+-6.16627
+-6.17189
+-6.17748
+-6.18303
+-6.18856
+-6.19405
+-6.19951
+-6.20493
+-6.21033
+-6.21569
+-6.22103
+-6.22633
+-6.2316
+-6.23683
+-6.24204
+-6.24721
+-6.25235
+-6.25746
+-6.26254
+-6.26759
+-6.2726
+-6.27758
+-6.28253
+-6.28745
+-6.29234
+-6.2972
+-6.30202
+-6.30681
+-6.31157
+-6.3163
+-6.321
+-6.32567
+-6.3303
+-6.3349
+-6.33947
+-6.34401
+-6.34852
+-6.35299
+-6.35744
+-6.36185
+-6.36623
+-6.37058
+-6.3749
+-6.37919
+-6.38344
+-6.38766
+-6.39185
+-6.39601
+-6.40014
+-6.40424
+-6.4083
+-6.41234
+-6.41634
+-6.42031
+-6.42425
+-6.42816
+-6.43203
+-6.43588
+-6.43969
+-6.44347
+-6.44722
+-6.45094
+-6.45463
+-6.45828
+-6.46191
+-6.4655
+-6.46906
+-6.47259
+-6.47609
+-6.47956
+-6.48299
+-6.4864
+-6.48977
+-6.49311
+-6.49642
+-6.4997
+-6.50295
+-6.50616
+-6.50935
+-6.5125
+-6.51562
+-6.51871
+-6.52177
+-6.5248
+-6.52779
+-6.53076
+-6.53369
+-6.53659
+-6.53947
+-6.5423
+-6.54511
+-6.54789
+-6.55064
+-6.55335
+-6.55603
+-6.55868
+-6.5613
+-6.56389
+-6.56645
+-6.56898
+-6.57147
+-6.57394
+-6.57637
+-6.57877
+-6.58114
+-6.58348
+-6.58578
+-6.58806
+-6.59031
+-6.59252
+-6.5947
+-6.59685
+-6.59897
+-6.60106
+-6.60312
+-6.60514
+-6.60714
+-6.6091
+-6.61103
+-6.61294
+-6.61481
+-6.61664
+-6.61845
+-6.62023
+-6.62197
+-6.62369
+-6.62537
+-6.62702
+-6.62864
+-6.63023
+-6.63179
+-6.63331
+-6.63481
+-6.63627
+-6.63771
+-6.63911
+-6.64048
+-6.64182
+-6.64313
+-6.6444
+-6.64565
+-6.64687
+-6.64805
+-6.6492
+-6.65032
+-6.65141
+-6.65247
+-6.6535
+-6.6545
+-6.65546
+-6.6564
+-6.6573
+-6.65817
+-6.65901
+-6.65982
+-6.6606
+-6.66135
+-6.66206
+-6.66275
+-6.6634
+-6.66403
+-6.66462
+-6.66518
+-6.66571
+-6.6662
+-6.66667
+-6.66711
+-6.66751
+-6.66789
+-6.66823
+-6.66854
+-6.66882
+-6.66907
+-6.66929
+-6.66947
+-6.66963
+-6.66975
+-6.66985
+-6.66991
+-6.66994
+-6.66994
+-6.66991
+-6.66985
+-6.66975
+-6.66963
+-6.66947
+-6.66929
+-6.66907
+-6.66882
+-6.66854
+-6.66823
+-6.66789
+-6.66751
+-6.66711
+-6.66667
+-6.6662
+-6.66571
+-6.66518
+-6.66462
+-6.66403
+-6.6634
+-6.66275
+-6.66206
+-6.66135
+-6.6606
+-6.65982
+-6.65901
+-6.65817
+-6.6573
+-6.6564
+-6.65546
+-6.6545
+-6.6535
+-6.65247
+-6.65141
+-6.65032
+-6.6492
+-6.64805
+-6.64687
+-6.64565
+-6.6444
+-6.64313
+-6.64182
+-6.64048
+-6.63911
+-6.63771
+-6.63627
+-6.63481
+-6.63331
+-6.63179
+-6.63023
+-6.62864
+-6.62702
+-6.62537
+-6.62369
+-6.62197
+-6.62023
+-6.61845
+-6.61664
+-6.61481
+-6.61294
+-6.61103
+-6.6091
+-6.60714
+-6.60514
+-6.60312
+-6.60106
+-6.59897
+-6.59685
+-6.5947
+-6.59252
+-6.59031
+-6.58806
+-6.58578
+-6.58348
+-6.58114
+-6.57877
+-6.57637
+-6.57394
+-6.57147
+-6.56898
+-6.56645
+-6.56389
+-6.5613
+-6.55868
+-6.55603
+-6.55335
+-6.55064
+-6.54789
+-6.54511
+-6.5423
+-6.53947
+-6.53659
+-6.53369
+-6.53076
+-6.52779
+-6.5248
+-6.52177
+-6.51871
+-6.51562
+-6.5125
+-6.50935
+-6.50616
+-6.50295
+-6.4997
+-6.49642
+-6.49311
+-6.48977
+-6.4864
+-6.48299
+-6.47956
+-6.47609
+-6.47259
+-6.46906
+-6.4655
+-6.46191
+-6.45828
+-6.45463
+-6.45094
+-6.44722
+-6.44347
+-6.43969
+-6.43588
+-6.43203
+-6.42816
+-6.42425
+-6.42031
+-6.41634
+-6.41234
+-6.4083
+-6.40424
+-6.40014
+-6.39601
+-6.39185
+-6.38766
+-6.38344
+-6.37919
+-6.3749
+-6.37058
+-6.36623
+-6.36185
+-6.35744
+-6.35299
+-6.34852
+-6.34401
+-6.33947
+-6.3349
+-6.3303
+-6.32567
+-6.321
+-6.3163
+-6.31157
+-6.30681
+-6.30202
+-6.2972
+-6.29234
+-6.28745
+-6.28253
+-6.27758
+-6.2726
+-6.26759
+-6.26254
+-6.25746
+-6.25235
+-6.24721
+-6.24204
+-6.23683
+-6.2316
+-6.22633
+-6.22103
+-6.21569
+-6.21033
+-6.20493
+-6.19951
+-6.19405
+-6.18856
+-6.18303
+-6.17748
+-6.17189
+-6.16627
+-6.16062
+-6.15494
+-6.14922
+-6.14347
+-6.1377
+-6.13188
+-6.12604
+-6.12017
+-6.11426
+-6.10832
+-6.10235
+-6.09635
+-6.09031
+-6.08425
+-6.07815
+-6.07202
+-6.06586
+-6.05966
+-6.05343
+-6.04717
+-6.04088
+-6.03456
+-6.0282
+-6.02182
+-6.0154
+-6.00895
+-6.00246
+-5.99595
+-5.9894
+-5.98282
+-5.9762
+-5.96956
+-5.96288
+-5.95617
+-5.94943
+-5.94266
+-5.93585
+-5.92901
+-5.92214
+-5.91524
+-5.9083
+-5.90134
+-5.89434
+-5.8873
+-5.88024
+-5.87314
+-5.86601
+-5.85885
+-5.85166
+-5.84443
+-5.83717
+-5.82988
+-5.82255
+-5.8152
+-5.80781
+-5.80039
+-5.79293
+-5.78545
+-5.77793
+-5.77038
+-5.76279
+-5.75518
+-5.74753
+-5.73985
+-5.73213
+-5.72439
+-5.71661
+-5.70879
+-5.70095
+-5.69307
+-5.68516
+-5.67722
+-5.66924
+-5.66124
+-5.65319
+-5.64512
+-5.63701
+-5.62887
+-5.6207
+-5.6125
+-5.60426
+-5.59599
+-5.58769
+-5.57935
+-5.57098
+-5.56258
+-5.55414
+-5.54568
+-5.53718
+-5.52864
+-5.52008
+-5.51148
+-5.50284
+-5.49418
+-5.48548
+-5.47675
+-5.46798
+-5.45919
+-5.45036
+-5.44149
+-5.43259
+-5.42366
+-5.4147
+-5.4057
+-5.39668
+-5.38761
+-5.37852
+-5.36939
+-5.36022
+-5.35103
+-5.3418
+-5.33254
+-5.32324
+-5.31391
+-5.30455
+-5.29515
+-5.28573
+-5.27626
+-5.26677
+-5.25724
+-5.24768
+-5.23808
+-5.22845
+-5.21879
+-5.20909
+-5.19936
+-5.18959
+-5.1798
+-5.16997
+-5.1601
+-5.1502
+-5.14027
+-5.1303
+-5.12031
+-5.11027
+-5.10021
+-5.0901
+-5.07997
+-5.0698
+-5.0596
+-5.04936
+-5.03909
+-5.02879
+-5.01845
+-5.00808
+-4.99768
+-4.98724
+-4.97676
+-4.96625
+-4.95571
+-4.94514
+-4.93453
+-4.92388
+-4.91321
+-4.90249
+-4.89175
+-4.88097
+-4.87015
+-4.8593
+-4.84842
+-4.8375
+-4.82655
+-4.81556
+-4.80454
+-4.79349
+-4.7824
+-4.77127
+-4.76012
+-4.74892
+-4.7377
+-4.72643
+-4.71514
+-4.70381
+-4.69244
+-4.68104
+-4.6696
+-4.65813
+-4.64663
+-4.63509
+-4.62352
+-4.61191
+-4.60026
+-4.58858
+-4.57687
+-4.56512
+-4.55334
+-4.54152
+-4.52966
+-4.51778
+-4.50585
+-4.49389
+-4.4819
+-4.46987
+-4.4578
+-4.44571
+-4.43357
+-4.4214
+-4.40919
+-4.39695
+-4.38468
+-4.37236
+-4.36002
+-4.34763
+-4.33521
+-4.32276
+-4.31027
+-4.29775
+-4.28518
+-4.27259
+-4.25995
+-4.24729
+-4.23458
+-4.22184
+-4.20907
+-4.19625
+-4.18341
+-4.17052
+-4.1576
+-4.14465
+-4.13165
+-4.11862
+-4.10556
+-4.09246
+-4.07932
+-4.06615
+-4.05294
+-4.03969
+-4.02641
+-4.01309
+-3.99973
+-3.98634
+-3.97291
+-3.95944
+-3.94594
+-3.9324
+-3.91883
+-3.90521
+-3.89156
+-3.87788
+-3.86415
+-3.85039
+-3.83659
+-3.82276
+-3.80888
+-3.79497
+-3.78103
+-3.76704
+-3.75302
+-3.73896
+-3.72486
+-3.71073
+-3.69656
+-3.68235
+-3.6681
+-3.65382
+-3.63949
+-3.62513
+-3.61073
+-3.5963
+-3.58182
+-3.56731
+-3.55276
+-3.53817
+-3.52354
+-3.50887
+-3.49417
+-3.47943
+-3.46465
+-3.44983
+-3.43497
+-3.42007
+-3.40513
+-3.39016
+-3.37515
+-3.36009
+-3.345
+-3.32987
+-3.3147
+-3.29949
+-3.28424
+-3.26895
+-3.25363
+-3.23826
+-3.22285
+-3.20741
+-3.19192
+-3.1764
+-3.16083
+-3.14523
+-3.12958
+-3.11389
+-3.09817
+-3.0824
+-3.0666
+-3.05075
+-3.03486
+-3.01894
+-3.00297
+-2.98696
+-2.97091
+-2.95482
+-2.93869
+-2.92251
+-2.9063
+-2.89004
+-2.87375
+-2.85741
+-2.84103
+-2.82461
+-2.80815
+-2.79164
+-2.7751
+-2.75851
+-2.74188
+-2.72521
+-2.70849
+-2.69174
+-2.67494
+-2.6581
+-2.64121
+-2.62429
+-2.60732
+-2.5903
+-2.57325
+-2.55615
+-2.53901
+-2.52182
+-2.50459
+-2.48732
+-2.47
+-2.45264
+-2.43524
+-2.41779
+-2.4003
+-2.38277
+-2.36518
+-2.34756
+-2.32989
+-2.31218
+-2.29442
+-2.27661
+-2.25877
+-2.24087
+-2.22293
+-2.20495
+-2.18692
+-2.16884
+-2.15072
+-2.13255
+-2.11434
+-2.09608
+-2.07777
+-2.05942
+-2.04102
+-2.02258
+-2.00408
+-1.98555
+-1.96696
+-1.94832
+-1.92964
+-1.91091
+-1.89214
+-1.87331
+-1.85444
+-1.83552
+-1.81655
+-1.79753
+-1.77847
+-1.75935
+-1.74019
+-1.72097
+-1.70171
+-1.6824
+-1.66304
+-1.64363
+-1.62417
+-1.60466
+-1.58509
+-1.56548
+-1.54582
+-1.52611
+-1.50634
+-1.48653
+-1.46666
+-1.44675
+-1.42678
+-1.40675
+-1.38668
+-1.36656
+-1.34638
+-1.32615
+-1.30586
+-1.28553
+-1.26514
+-1.2447
+-1.2242
+-1.20365
+-1.18305
+-1.16239
+-1.14168
+-1.12091
+-1.10009
+-1.07921
+-1.05828
+-1.03729
+-1.01625
+-0.995149
+-0.973994
+-0.952783
+-0.931515
+-0.91019
+-0.888809
+-0.86737
+-0.845874
+-0.824321
+-0.802709
+-0.781039
+-0.759311
+-0.737524
+-0.715679
+-0.693774
+-0.67181
+-0.649786
+-0.627703
+-0.605559
+-0.583355
+-0.561091
+-0.538765
+-0.516378
+-0.49393
+-0.471421
+-0.448849
+-0.426215
+-0.403519
+-0.380759
+-0.357937
+-0.335051
+-0.312102
+-0.289089
+-0.266012
+-0.24287
+-0.219664
+-0.196392
+-0.173055
+-0.149652
+-0.126184
+-0.102649
+-0.0790472
+-0.0553788
+-0.0316432
+-0.00784012
+0.0160308
+0.0399699
+0.0639775
+0.0880539
+0.112199
+0.136414
+0.160699
+0.185054
+0.20948
+0.233976
+0.258543
+0.283182
+0.307893
+0.332676
+0.357531
+0.38246
+0.407461
+0.432537
+0.457686
+0.48291
+0.508208
+0.533582
+0.559031
+0.584556
+0.610157
+0.635835
+0.66159
+0.687422
+0.713332
+0.739321
+0.765388
+0.791534
+0.817759
+0.844065
+0.87045
+0.896916
+0.923464
+0.950092
+0.976803
+1.0036
+1.03047
+1.05743
+1.08447
+1.1116
+1.13881
+1.16611
+1.19349
+1.22096
+1.24851
+1.27615
+1.30388
+1.3317
+1.3596
+1.3876
+1.41568
+1.44385
+1.47211
+1.50046
+1.52891
+1.55744
+1.58607
+1.61479
+1.6436
+1.67251
+1.70151
+1.7306
+1.75979
+1.78908
+1.81846
+1.84794
+1.87752
+1.90719
+1.93696
+1.96683
+1.9968
+2.02687
+2.05705
+2.08732
+2.11769
+2.14817
+2.17874
+2.20943
+2.24021
+2.2711
+2.3021
+2.3332
+2.3644
+2.39572
+2.42714
+2.45867
+2.4903
+2.52205
+2.5539
+2.58587
+2.61795
+2.65014
+2.68244
+2.71485
+2.74738
+2.78002
+2.81278
+2.84565
+2.87864
+2.91174
+2.94496
+2.9783
+3.01176
+3.04534
+3.07903
+3.11285
+3.14679
+3.18085
+3.21504
+3.24934
+3.28377
+3.31833
+3.35301
+3.38781
+3.42275
+3.45781
+3.493
+3.52831
+3.56376
+3.59934
+3.63504
+3.67088
+3.70685
+3.74296
+3.7792
+3.81557
+3.85208
+3.88872
+3.9255
+3.96242
+3.99947
+4.03667
+4.074
+4.11147
+4.14909
+4.18684
+4.22474
+4.26278
+4.30097
+4.3393
+4.37778
+4.4164
+4.45517
+4.49408
+4.53315
+4.57236
+4.61172
+4.65124
+4.6909
+4.73072
+4.77069
+4.81081
+4.85109
+4.89152
+4.93211
+4.97285
+5.01375
+5.05481
+5.09603
+5.1374
+5.17894
+5.22064
+5.26249
+5.30451
+5.3467
+5.38904
+5.43155
+5.47422
+5.51706
+5.56007
+5.60324
+5.64658
+5.69009
+5.73376
+5.77761
+5.82162
+5.8658
+5.91016
+5.95468
+5.99938
+6.04425
+6.08929
+6.13451
+6.17989
+6.22546
+6.2712
+6.31711
+6.3632
+6.40946
+6.4559
+6.50252
+6.54931
+6.59628
+6.64343
+6.69076
+6.73826
+6.78594
+6.8338
+6.88184
+6.93005
+6.97845
+7.02702
+7.07577
+7.1247
+7.17381
+7.2231
+7.27256
+7.32221
+7.37203
+7.42202
+7.4722
+7.52255
+7.57308
+7.62378
+7.67466
+7.72571
+7.77694
+7.82834
+7.87992
+7.93166
+7.98358
+8.03567
+8.08792
+8.14035
+8.19294
+8.2457
+8.29863
+8.35171
+8.40496
+8.45837
+8.51194
+8.56567
+8.61956
+8.67359
+8.72779
+8.78213
+8.83662
+8.89126
+8.94604
+9.00096
+9.05603
+9.11123
+9.16657
+9.22204
+9.27764
+9.33337
+9.38922
+9.44519
+9.50128
+9.55749
+9.61381
+9.67023
+9.72676
+9.78339
+9.84012
+9.89695
+9.95386
+10.0109
+10.0679
+10.1251
+10.1823
+10.2396
+10.297
+10.3544
+10.4119
+10.4694
+10.527
+10.5846
+10.6422
+10.6999
+10.7575
+10.8152
+10.8729
+10.9307
+10.9884
+11.0461
+11.1038
+11.1614
+11.2191
+11.2767
+11.3342
+11.3917
+11.4492
+11.5066
+11.5639
+11.6212
+11.6783
+11.7354
+11.7923
+11.8492
+11.9059
+11.9625
+12.019
+12.0753
+12.1315
+12.1875
+12.2434
+12.299
+12.3545
+12.4098
+12.4649
+12.5198
+12.5744
+12.6288
+12.683
+12.7369
+12.7905
+12.8439
+12.897
+12.9498
+13.0023
+13.0545
+13.1063
+13.1579
+13.2091
+13.2599
+13.3104
+13.3604
+13.4102
+13.4595
+13.5084
+13.5569
+13.6049
+13.6525
+13.6997
+13.7464
+13.7927
+13.8385
+13.8837
+13.9285
+13.9728
+14.0165
+14.0597
+14.1024
+14.1445
+14.1861
+14.2271
+14.2675
+14.3073
+14.3465
+14.385
+14.423
+14.4603
+14.497
+14.5331
+14.5685
+14.6032
+14.6372
+14.6705
+14.7032
+14.7351
+14.7663
+14.7968
+14.8266
+14.8556
+14.8839
+14.9114
+14.9382
+14.9641
+14.9894
+15.0138
+15.0374
+15.0602
+15.0823
+15.1035
+15.1239
+15.1435
+15.1622
+15.1801
+15.1972
+15.2134
+15.2288
+15.2433
+15.257
+15.2698
+15.2817
+15.2928
+15.3029
+15.3123
+15.3207
+15.3282
+15.3349
+15.3407
+15.3456
+15.3496
+15.3527
+15.355
+15.3563
+15.3567
diff --git a/examples/DDG4/scripts/Channeling.py b/examples/DDG4/scripts/Channeling.py
new file mode 100644
index 000000000..fff518cc9
--- /dev/null
+++ b/examples/DDG4/scripts/Channeling.py
@@ -0,0 +1,81 @@
+# ==========================================================================
+#  AIDA Detector description implementation
+# --------------------------------------------------------------------------
+# Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
+# All rights reserved.
+#
+# For the licensing terms see $DD4hepINSTALL/LICENSE.
+# For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
+#
+# ==========================================================================
+#
+from __future__ import absolute_import, unicode_literals
+import os
+import sys
+import DDG4
+from DDG4 import OutputLevel as Output
+from g4units import keV
+#
+#
+"""
+
+   dd4hep simulation example setup using the python configuration
+
+   @author  M.Frank
+   @version 1.0
+
+"""
+
+
+def run():
+  kernel = DDG4.Kernel()
+  install_dir = os.environ['DD4hepExamplesINSTALL']
+  kernel.loadGeometry(str("file:" + install_dir + "/examples/DDG4/compact/Channeling.xml"))
+
+  DDG4.importConstants(kernel.detectorDescription(), debug=False)
+  geant4 = DDG4.Geant4(kernel, tracker='Geant4TrackerCombineAction')
+  geant4.printDetectors()
+  # Configure UI
+  if len(sys.argv) > 1:
+    geant4.setupCshUI(macro=sys.argv[1])
+  else:
+    geant4.setupCshUI()
+
+  # Configure field
+  geant4.setupTrackingField(prt=True)
+  # Configure Event actions
+  prt = DDG4.EventAction(kernel, 'Geant4ParticlePrint/ParticlePrint')
+  prt.OutputLevel = Output.DEBUG
+  prt.OutputType = 3  # Print both: table and tree
+  kernel.eventAction().adopt(prt)
+
+  generator_output_level = Output.INFO
+
+  # Configure G4 geometry setup
+  seq, act = geant4.addDetectorConstruction("Geant4DetectorGeometryConstruction/ConstructGeo")
+  act.DebugMaterials = True
+  act.DebugElements = False
+  act.DebugVolumes = True
+  act.DebugShapes = True
+  act.DebugSurfaces = True
+
+  # Setup particle gun
+  gun = geant4.setupGun("Gun", particle='gamma', energy=5 * keV, multiplicity=1)
+  gun.OutputLevel = generator_output_level
+
+  geant4.setupTracker('ChannelingDevice')
+
+  # Now build the physics list:
+  phys = geant4.setupPhysics('QGSP_BERT')
+  ph = DDG4.PhysicsList(kernel, 'Channeling')
+  ph.addPhysicsConstructor('Geant4ChannelingPhysics')
+  ph.enableUI()
+  phys.adopt(ph)
+
+  phys.dump()
+
+  geant4.execute()
+
+
+if __name__ == "__main__":
+  run()
diff --git a/examples/DDG4/src/ChannelingCrystalMaterial.cpp b/examples/DDG4/src/ChannelingCrystalMaterial.cpp
new file mode 100644
index 000000000..fe8a134d9
--- /dev/null
+++ b/examples/DDG4/src/ChannelingCrystalMaterial.cpp
@@ -0,0 +1,81 @@
+//==========================================================================
+//  AIDA Detector description implementation 
+//--------------------------------------------------------------------------
+// Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
+// All rights reserved.
+//
+// For the licensing terms see $DD4hepINSTALL/LICENSE.
+// For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
+//
+// Author     : M.Frank
+//
+//==========================================================================
+
+// Framework include files
+#include <XML/XMLElements.h>
+#include <DDG4/Factories.h>
+#include <DD4hep/ComponentProperties.h>
+
+// Geant4 include files
+#include <G4SystemOfUnits.hh>
+#include <G4ExtendedMaterial.hh>
+#include <G4CrystalExtension.hh>
+#include <G4ChannelingMaterialData.hh>
+#include <G4ChannelingOptrMultiParticleChangeCrossSection.hh>
+
+// Forward declarations
+namespace {  class ChannelingCrystalMaterial;  }
+
+/// Namespace example name of the user
+namespace dd4hep   {
+
+  /// Class member specialization to create G4ExtendedMaterial
+  /** Class member specialization to create G4ExtendedMaterial
+   *
+   *  \author  M.Frank
+   *  \version 1.0
+   *  \ingroup DD4HEP_SIMULATION
+   */
+  template <>  G4Material* 
+  Geant4MaterialFactory<ChannelingCrystalMaterial>::create(dd4hep::Detector& /* description */,
+							   Material    mat,
+							   G4Material* base_material)
+  {
+    std::string crystal_data, orientation, name = base_material->GetName();
+    auto* material = new G4ExtendedMaterial(name+".extended", base_material);
+    material->RegisterExtension(std::unique_ptr<G4CrystalExtension>(new G4CrystalExtension(material)));
+    auto* ext = (G4CrystalExtension*)material->RetrieveExtension("crystal");
+    ext->SetUnitCell(new G4CrystalUnitCell(5.43 * CLHEP::angstrom,
+					   5.43 * CLHEP::angstrom,
+					   5.43 * CLHEP::angstrom,
+					   CLHEP::halfpi,
+					   CLHEP::halfpi,
+					   CLHEP::halfpi,
+					   227));
+
+    material->RegisterExtension(std::unique_ptr<G4ChannelingMaterialData>(new G4ChannelingMaterialData("channeling")));
+    auto* channelingData =(G4ChannelingMaterialData*)material->RetrieveExtension("channeling");
+    crystal_data = mat.constPropertyRef("Geant4-ignore:crystal_data");
+    if ( crystal_data.empty() )   {
+      except("ChannelingCrystalMaterial",
+	     "====> Extended material: %s misses channeling data file.",
+	     name.c_str());
+    }
+    channelingData->SetFilename(crystal_data);
+
+    /// Set crystal orientation
+    orientation = mat.propertyRef("Geant4-ignore:crystal_orientation");
+    if ( !orientation.empty() )   {
+      Direction direction;
+      Property  prop(direction);
+      prop.str(orientation);
+      channelingData->SetBR(direction.X());
+    }
+
+    printout(ALWAYS,"ChannelingCrystalMaterial",
+	     "====> Created extended material: %s Data: %s Orientation: %s",
+	     name.c_str(), crystal_data.c_str(), orientation.c_str());
+    return material;
+  }
+}      // End namespace dd4hep
+DECLARE_GEANT4EXTENDEDMATERIAL(ChannelingCrystalMaterial)
diff --git a/examples/DDG4/src/ChannelingCrystalVolume.cpp b/examples/DDG4/src/ChannelingCrystalVolume.cpp
new file mode 100644
index 000000000..41be0d580
--- /dev/null
+++ b/examples/DDG4/src/ChannelingCrystalVolume.cpp
@@ -0,0 +1,52 @@
+//==========================================================================
+//  AIDA Detector description implementation 
+//--------------------------------------------------------------------------
+// Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
+// All rights reserved.
+//
+// For the licensing terms see $DD4hepINSTALL/LICENSE.
+// For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
+//
+// Author     : M.Frank
+//
+//==========================================================================
+
+// Framework include files
+#include <DDG4/Factories.h>
+// Geant4 include files
+#include <G4ExtendedMaterial.hh>
+#include <G4LogicalCrystalVolume.hh>
+
+// Forward declarations
+namespace  {  class ChannelingCrystalVolume;  }
+
+/// Namespace example name of the user
+namespace dd4hep   {
+
+  /// Class member specialization to create a G4LogicalCrystalVolume
+  /** Class member specialization to create a G4LogicalCrystalVolume
+   *  
+   *  \author  M.Frank
+   *  \version 1.0
+   *  \ingroup DD4HEP_SIMULATION
+   */
+  template <> G4LogicalVolume* 
+  Geant4LogicalVolumeFactory<ChannelingCrystalVolume>::create(dd4hep::Detector& /* description */,
+							      Volume      volume,
+							      G4VSolid*   solid,
+							      G4Material* material)
+  {
+    G4ExtendedMaterial* mat = dynamic_cast<G4ExtendedMaterial*>(material);
+    if ( !mat )   {
+      except("ChannelingCrystalVolume", 
+	     "====> Material %s is no G4ExtendedMaterial!", material->GetName().c_str());
+    }
+    auto* ptr = new G4LogicalCrystalVolume(solid, mat, volume.name());
+    printout(ALWAYS,"ChannelingCrystalVolume",
+	     "====> Created specialize logical volume [G4LogicalCrystalVolume]: %s",
+	     volume.name());
+    return ptr;
+  }
+}      // End namespace dd4hep
+
+DECLARE_GEANT4LOGICALVOLUME(ChannelingCrystalVolume)
-- 
GitLab