macro/dirc_vis.mac

+# npsim --runType vis --compactFile athena.xml --random.seed 1 --enableGun --gun.energy 2*GeV --gun.thetaMin 90*deg --gun.thetaMax 90*deg --gun.distribution uniform --macro macro/dirc_vis.mac --outputFile dirc_vis.root
+
+/vis/list
+
+# Use this open statement to create an OpenGL view:
+/vis/open OGL 600x600-0+0
+
+
+#
+# Use this open statement to create a .prim file suitable for
+# viewing in DAWN:
+#/vis/open DAWNFILE
+#
+# Use this open statement to create a .heprep file suitable for
+# viewing in HepRApp:
+#/vis/open HepRepFile
+#
+# Use this open statement to create a .wrl file suitable for
+# viewing in a VRML viewer:
+#/vis/open VRML2FILE
+#
+# Disable auto refresh and quieten vis messages whilst scene and
+# trajectories are established:
+/vis/viewer/set/autoRefresh false
+/vis/viewer/set/rotationStyle freeRotation
+/vis/verbose errors
+#
+# Draw geometry:
+/vis/drawVolume
+
+# Specify zoom value:
+/vis/viewer/zoom 1.4
+
+# Specify style (surface or wireframe):
+#/vis/viewer/set/style wireframe
+
+#top
+# /vis/viewer/set/viewpointThetaPhi -90. 0.
+# /vis/viewer/set/viewpointVector 0 1 0.00000001
+# /vis/viewer/set/upVector 1 0 0
+
+#mcp
+#/vis/viewer/set/viewpointVector -1 0 0.00000001
+#/vis/viewer/set/upVector 0 1 0
+
+#side
+# /vis/viewer/set/viewpointVector 0 1 0.00000001
+# /vis/viewer/set/upVector 0 -1 0
+
+#lenses
+#/vis/viewer/set/viewpointThetaPhi 30 30
+#
+#/vis/viewer/set/projection orthogonal
+#/vis/viewer/zoomTo 4.5
+#/vis/viewer/set/targetPoint -110 20 -10 cm
+#/vis/viewer/set/background white
+
+
+# Draw coordinate axes:
+/vis/scene/add/axes 0 0 0 1 m
+#
+# Draw smooth trajectories at end of event, showing trajectory points
+# as markers 2 pixels wide:
+/vis/scene/add/trajectories smooth
+#/vis/modeling/trajectories/create/drawByCharge
+#/vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true
+#/vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 2
+
+/vis/modeling/trajectories/create/drawByParticleID
+/vis/modeling/trajectories/drawByParticleID-0/set opticalphoton yellow
+/vis/modeling/trajectories/drawByParticleID-0/setRGBA opticalphoton 1 0.8 0 0.3
+/vis/modeling/trajectories/drawByParticleID-0/setRGBA e+ 1 0 1 1
+/vis/modeling/trajectories/drawByParticleID-0/default/setDrawStepPts true
+/vis/modeling/trajectories/drawByParticleID-0/default/setStepPtsSize 2
+
+
+
+# (if too many tracks cause core dump => /tracking/storeTrajectory 0)
+#
+# Draw hits at end of event:
+#/vis/scene/add/hits
+#
+# To draw only gammas:
+#/vis/filtering/trajectories/create/particleFilter
+#/vis/filtering/trajectories/particleFilter-0/add gamma
+#
+# To invert the above, drawing all particles except gammas,
+# keep the above two lines but also add:
+#/vis/filtering/trajectories/particleFilter-0/invert true
+#
+# Many other options are available with /vis/modeling and /vis/filtering.
+# For example, to select colour by particle ID:
+#/vis/modeling/trajectories/create/drawByParticleID
+#/vis/modeling/trajectories/drawByParticleID-0/set e- blue
+#
+# To superimpose all of the events from a given run:
+/vis/scene/endOfEventAction accumulate
+
+# Re-establish auto refreshing and verbosity:
+/vis/viewer/set/autoRefresh true
+/vis/verbose warnings
+#
+# For file-based drivers, use this to create an empty detector view:
+#/vis/viewer/flush
+
+/vis/viewer/set/projection p 40 deg
+/vis/viewer/set/style s
+/vis/viewer/set/hiddenEdge true
+#/vis/viewer/set/background 1 1 1 0
+#/vis/viewer/set/background 0 0 0 0
+/vis/viewer/set/hiddenMarker 
+/vis/geometry/set/forceLineSegmentsPerCircle all 0 50
+# /vis/viewer/set/globalLineWidthScale 1.5
+# /vis/geometry/set/lineWidth all 0 1
+
+/vis/ogl/set/displayListLimit 500000
+
+
+#fp
+
+#/vis/viewer/save
+
+# #mcp
+# /vis/viewer/set/viewpointVector -1 0 0.00000001
+# /vis/viewer/set/upVector 0 1
+# /vis/viewer/set/viewpointThetaPhi 9  41  deg
+# /vis/viewer/zoomTo 6
+# /vis/viewer/set/targetPoint 88.5358 -17.0546 18.8692 cm 
+
+
+# lens
+# /vis/viewer/zoomTo 20
+# /vis/viewer/set/viewpointThetaPhi 129.792  1.39204  deg
+# /vis/viewer/set/targetPoint 100 20 210 cm
+
+# # #lens side
+# /vis/viewer/set/viewpointVector 1.31839e-16 -1 -6.93889e-18
+# /vis/viewer/set/upVector 0 1 0
+# /vis/viewer/set/projection orthogonal
+# /vis/viewer/zoomTo 11
+# /vis/viewer/scaleTo 1 1 1
+# /vis/viewer/set/targetPoint 1.01299 -0.569503 2.15227 m  
+
+
+# side wedge
+# /vis/viewer/set/viewpointVector -5.20417e-17 1 -3.46945e-18
+# /vis/viewer/set/upVector 0 1 0
+# /vis/viewer/set/projection orthogonal
+# /vis/viewer/zoomTo 13.4
+# /vis/viewer/scaleTo 1 1 1
+# /vis/viewer/set/targetPoint 1.0404 -0.128502 2.17344 m
+
+# whole length
+# /vis/viewer/zoomTo 7.0883
+# /vis/viewer/set/targetPoint 1.11348 -0.128502 0.164898 m
+#/vis/ogl/export ev_1.png
+
+# draw field
+#/vis/scene/add/magneticField 10
+
+#/vis/scene/create
+#/vis/scene/add/volume
+#
+#/vis/sceneHandler/create VRML2FILE
+#
+#/vis/viewer/create
+#
+#/tracking/storeTrajectory 1
+#/vis/scene/add/trajectories
+#/vis/scene/endOfEventAction accumulate
+
+/vis/viewer/refresh
+
+
+/gps/verbose 2
+/gps/particle pi-
+/gps/number 1
+
+#/gps/ene/type Gauss
+/gps/ene/mono  6.0 GeV
+#/gps/ene/sigma 3.0 GeV
+
+#/gps/pos/type Volume
+#/gps/pos/shape Cylinder
+#/gps/pos/centre 0.0 0.0 0.0 cm
+#/gps/pos/radius 0.01 cm
+#/gps/pos/halfz  10 cm
+/gps/position 0 0 0 cm
+#/gps/direction 1 0 0
+
+/gps/direction -1 0 0
+#/gps/ang/type iso
+#/gps/ang/mintheta  10 degree
+#/gps/ang/maxtheta  170 degree
+
+/vis/viewer/set/viewpointVector 0 -1 0
+/vis/viewer/set/upVector -1 0 0
+/vis/viewer/set/projection orthogonal
+/vis/viewer/zoomTo 6.8984
+/vis/viewer/scaleTo 1 1 1
+/vis/viewer/set/targetPoint -0.987804 0.579407 2.01758 m
+
+/run/beamOn 1
\ No newline at end of file

macro/gps.mac

+/control/verbose 2
+/run/initialize
+
+/gps/verbose 2
+/gps/particle pi-
+/gps/number 1
+
+/gps/ene/type Gauss
+/gps/ene/mono  10.0 GeV
+/gps/ene/sigma 3.0 GeV
+
+/gps/pos/type Volume
+/gps/pos/shape Cylinder
+/gps/pos/centre 0.0 0.0 0.0 cm
+/gps/pos/radius 0.01 cm
+/gps/pos/halfz  10 cm
+/gps/position 0 0 0 cm
+
+#/gps/direction 0 0.1 1.0 
+/gps/ang/type iso
+/gps/ang/mintheta  10 degree
+/gps/ang/maxtheta  170 degree
+
+/run/beamOn 100

macro/mrich_vis.mac

+/run/initialize
+
+#/vis/open OGL 800x800-0+0
+/vis/open OGLIQt 800x800-0+0
+#
+/vis/drawVolume
+/vis/viewer/set/viewpointThetaPhi 20 30
+/vis/viewer/zoom 1.1
+#/vis/viewer/set/style wireframe
+#/vis/scene/add/axes 0 0 0 1 m
+/vis/scene/add/trajectories rich smooth
+/vis/modeling/trajectories/create/drawByCharge
+#/vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true
+#/vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 2
+/vis/scene/add/hits
+/vis/scene/endOfEventAction accumulate 2 
+
+#/vis/filtering/trajectories/create/particleFilter
+#/vis/filtering/trajectories/particleFilter-0/add gamma
+#/vis/filtering/trajectories/particleFilter-0/invert true
+#/vis/filtering/trajectories/particleFilter-0/verbose true
+#/vis/filtering/trajectories/particleFilter-0/active true
+
+
+/vis/ogl/set/displayListLimit 50000
+
+/vis/viewer/flush
+
+
+/gps/verbose 2
+/gps/particle pi-
+/gps/number 1
+
+/gps/ene/type Gauss
+/gps/ene/mono  8.0 GeV
+/gps/ene/sigma 1.0 GeV
+
+/gps/pos/type Volume
+/gps/pos/shape Cylinder
+/gps/pos/centre 0.0 0.0 0.0 cm
+/gps/pos/radius 0.001 cm
+/gps/pos/halfz  1 cm
+/gps/position 0 0 0 cm
+
+#/gps/direction 0 0.1 1.0 
+/gps/ang/type iso
+/gps/ang/mintheta  10 degree
+/gps/ang/maxtheta  20 degree
+
+/run/beamOn 1

macro/vis.mac

+/vis/open OGL 800x800-0+0
+
+/vis/drawVolume
+/vis/viewer/set/viewpointThetaPhi 30 30
+#/vis/viewer/zoom 30.
+#/vis/viewer/set/style wireframe
+#/vis/scene/add/axes 0 0 0 1 m
+/vis/scene/add/trajectories rich smooth
+/vis/modeling/trajectories/create/drawByCharge
+#/vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true
+#/vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 2
+/vis/scene/add/hits
+/vis/scene/endOfEventAction accumulate 2 
+
+#/vis/filtering/trajectories/create/particleFilter
+#/vis/filtering/trajectories/particleFilter-0/add gamma
+#/vis/filtering/trajectories/particleFilter-0/invert true
+#/vis/filtering/trajectories/particleFilter-0/verbose true
+#/vis/filtering/trajectories/particleFilter-0/active true
+
+
+/vis/ogl/set/displayListLimit 500000
+
+/vis/viewer/flush
+
+/run/beamOn 1
+
+#/control/execute macro/gun.mac
+

macro/vis_OGL_along_z_beamOn_1.mac

+# Usage example:
+#   npsim --runType vis --compactFile athena.xml --random.seed 1 --enableGun --gun.energy 2*GeV --gun.thetaMin 90*deg --gun.thetaMax 90*deg --gun.distribution uniform --macro macro/vis_OGL_along_z_beamOn_1.mac --outputFile test.root
+/vis/open OGL 1920x1080
+/vis/scene/create
+/vis/scene/add/volume
+/vis/sceneHandler/attach
+/vis/viewer/zoom 10
+/vis/viewer/flush
+/vis/scene/add/trajectories
+/vis/scene/add/hits
+/run/beamOn 1

reports/templates/detector_views.md

+#AthenaDetector
+
+
+{% for image_path in images %}
+<img src="{{ base_url }}/{{ image_path }}?job=report" width="400px" />
+<br \>
+{% endfor %}

scripts/checkOverlaps.py

+from __future__ import absolute_import, unicode_literals
+import os
+import time
+import logging
+
+import argparse
+parser = argparse.ArgumentParser(
+     prog='checkOverlaps.py',
+     description='''Check for overlaps using Geant4''',
+     epilog='''
+     This program checks the compact detector file for overlaps using Geant4.
+         ''')
+parser.add_argument("-c", "--compact", help="compact detector file",default="athena.xml")
+parser.add_argument("-r", "--resolution", help="number of points on surface",default="10000")
+parser.add_argument("-t", "--tolerance", help="minimum distance (in mm) to report overlaps",default="0.1")
+parser.add_argument("-v", "--verbose", help="print output", action='store_true')
+
+args = parser.parse_args()
+
+import DDG4
+from g4units import keV, GeV, mm, ns, MeV
+
+def run():
+  kernel = DDG4.Kernel()
+  description = kernel.detectorDescription()
+  kernel.loadGeometry(str("file:" + args.compact))
+
+  DDG4.importConstants(description)
+
+  geant4 = DDG4.Geant4(kernel)
+  ui = geant4.setupCshUI(ui=None)
+  ui.Commands = [
+      '/geometry/test/resolution {}'.format(args.resolution),
+      '/geometry/test/tolerance {}'.format(args.tolerance),
+      '/geometry/test/verbosity {}'.format(1 if args.verbose else 0),
+      '/geometry/test/run'
+      ]
+  kernel.configure()
+  kernel.initialize()
+  kernel.run()
+  kernel.terminate()
+
+
+if __name__ == "__main__":
+  run()

scripts/convert_to_gdml.py

+from __future__ import absolute_import, unicode_literals
+import os
+import time
+import logging
+
+import argparse
+parser = argparse.ArgumentParser(
+     prog='convert_to_gdml.py',
+     description='''Convert DD4Hep description to GDML''',
+     epilog='''
+     This program converts the compact detector file to a single GDML file.
+         ''')
+parser.add_argument("-c", "--compact", help="compact detector file",default="athena.xml")
+parser.add_argument("-o", "--output", help="gdml detector file",default="athena.gdml")
+
+args = parser.parse_args()
+
+import DDG4
+from g4units import keV, GeV, mm, ns, MeV
+
+def run():
+  kernel = DDG4.Kernel()
+  description = kernel.detectorDescription()
+  kernel.loadGeometry(str("file:" + args.compact))
+
+  DDG4.importConstants(description)
+
+  geant4 = DDG4.Geant4(kernel)
+  ui = geant4.setupCshUI(ui=None)
+  #
+  # Setup the GDML writer action
+  writer = DDG4.Action(kernel, 'Geant4GDMLWriteAction/Writer')
+  writer.enableUI()
+  kernel.registerGlobalAction(writer)
+  ui.Commands = [
+      '/ddg4/Writer/Output {}'.format(args.output),
+      '/ddg4/Writer/OverWrite 1',
+      '/ddg4/Writer/write'
+      ]
+  kernel.configure()
+  kernel.initialize()
+  kernel.run()
+  kernel.terminate()
+
+
+if __name__ == "__main__":
+  run()

scripts/example_dd4hep.cxx

+#include "DD4hep/Detector.h"
+#include "DDG4/Geant4Data.h"
+#include "DDRec/CellIDPositionConverter.h"
+#include "DDRec/SurfaceManager.h"
+#include "DDRec/Surface.h"
+
+
+#include "fmt/core.h"
+
+
+/** Example using dd4hep
+ */
+void example_dd4hep(const char* fname = "test_tracker_disc.root"){
+
+  using namespace ROOT::Math;
+
+  // -------------------------
+  // Get the DD4hep instance
+  // Load the compact XML file
+  // Initialize the position converter tool
+  dd4hep::Detector& detector = dd4hep::Detector::getInstance();
+  detector.fromCompact("athena.xml");
+  dd4hep::rec::CellIDPositionConverter cellid_converter(detector);
+
+  fmt::print("Detector Types:\n");
+  for(const auto&  dtype : detector.detectorTypes() ) {
+    fmt::print("  {}\n", dtype); 
+  }
+
+  fmt::print("\n");
+  fmt::print("All detector subsystem names:\n");
+  for(const auto&  d : detector.detectors() ) {
+    fmt::print("  {}\n", d.first);
+  }
+
+}

scripts/input_data/central_events.hepmc

+HepMC::Version 3.02.02
+HepMC::Asciiv3-START_EVENT_LISTING
+E 0 4 9
+U GEV CM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 -9.9999999869440064e+00 1.0000000000000000e+01 5.1099888971089147e-04 4
+P 2 1 22 -5.8225640056623670e-03 3.7747924008719728e-02 -2.9249301763891555e+00 2.9248270767857321e+00 -4.5408321610621201e-02 3
+P 3 1 11 5.8225640056623670e-03 -3.7747924008719728e-02 -7.0750698105548508e+00 7.0751729232142679e+00 5.1101967722736764e-04 1
+P 4 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 9.9995598131387851e+01 1.0000000000000000e+02 9.3827199999876709e-01 4
+P 5 4 2212 0.0000000000000000e+00 0.0000000000000000e+00 9.9995598131482922e+01 1.0000000000009507e+02 9.3827199999973643e-01 3
+V -3 0 [2,5]
+P 6 -3 553 -1.4745608418116361e-02 -8.6828685571490272e-01 4.7980031183979008e+00 1.0642943648781511e+01 9.4603000000109834e+00 2
+P 7 -3 2212 8.9230444123625328e-03 9.0603477972421620e-01 9.2272664838178372e+01 9.2281883429672035e+01 9.3827200000070576e-01 1
+P 8 6 -13 4.4907597498767587e-01 2.0906320322010372e+00 -2.1757837479786910e+00 3.0524760533554351e+00 1.0565800000001382e-01 1
+P 9 6 13 -4.6382158340577534e-01 -2.9589188879149320e+00 6.9737868663710216e+00 7.5904675954137195e+00 1.0565800000002223e-01 1
+E 0 4 10
+U GEV CM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 -9.9999999869440064e+00 1.0000000000000000e+01 5.1099888971089147e-04 4
+P 2 1 22 -2.6839726606050791e-04 1.1264659797828775e-03 -6.3401069726794805e+00 6.3401067668601172e+00 -1.9876609922101014e-03 3
+P 3 1 11 2.6839726606050791e-04 -1.1264659797828775e-03 -3.6598930142645258e+00 3.6598932331398828e+00 5.1103633096365880e-04 1
+P 4 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 9.9995598131387851e+01 1.0000000000000000e+02 9.3827199999876709e-01 4
+P 5 4 2212 0.0000000000000000e+00 0.0000000000000000e+00 9.9995598131482922e+01 1.0000000000009507e+02 9.3827199999973643e-01 3
+V -3 0 [2,5]
+P 6 -3 553 1.9156252814767760e-01 2.1790735768546463e-01 -2.8073109612741973e+00 9.8723072682332713e+00 9.4603000029807429e+00 2
+P 7 -3 2212 -1.9183092541375127e-01 -2.1678089170561543e-01 9.6462802125352582e+01 9.6467799504704473e+01 9.3827199999779787e-01 1
+P 8 6 -13 -2.4852998465289872e-01 2.3091049977702758e+00 -5.7337843777899966e+00 6.1871786743193251e+00 1.0565799999995498e-01 1
+P 9 6 13 4.4342149179344215e-01 -2.1010576035711188e+00 2.9552377554469285e+00 3.6545395711486477e+00 1.0565800000003905e-01 1
+P 10 6 22 -3.3289790532426214e-03 9.8599634176282647e-03 -2.8764338046322176e-02 3.0589019653741657e-02 3.2927225399135965e-10 1
+E 0 4 9
+U GEV CM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 -9.9999999869440064e+00 1.0000000000000000e+01 5.1099888971089147e-04 4
+P 2 1 22 2.1549001953953920e-01 5.3446029602846556e-01 -6.2965867413439192e-01 6.1195545733789913e-01 -5.9503304168678661e-01 3
+P 3 1 11 -2.1549001953953920e-01 -5.3446029602846556e-01 -9.3703413128096145e+00 9.3880445426621009e+00 5.0981276742689016e-04 1
+P 4 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 9.9995598131387851e+01 1.0000000000000000e+02 9.3827199999876709e-01 4
+P 5 4 2212 0.0000000000000000e+00 0.0000000000000000e+00 9.9995598131482922e+01 1.0000000000009507e+02 9.3827199999973643e-01 3
+V -3 0 [2,5]
+P 6 -3 553 2.2825834022061950e-01 3.3334310888886032e-01 3.5569092662904950e+01 3.6807891116586688e+01 9.4603000000031017e+00 2
+P 7 -3 2212 -1.2768320678757170e-02 2.0111718714536578e-01 6.3796846795456545e+01 6.3804064341873378e+01 9.3827200000094813e-01 1
+P 8 6 -13 -3.2477759288923482e+00 3.3065214111492942e+00 1.3631832024299879e+01 1.4398581189769127e+01 1.0565800000146809e-01 1
+P 9 6 13 3.4760342691128936e+00 -2.9731783022605436e+00 2.1937260638593418e+01 2.2409309926805506e+01 1.0565799999797114e-01 1
+E 0 4 9
+U GEV CM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 -9.9999999869440064e+00 1.0000000000000000e+01 5.1099888971089147e-04 4
+P 2 1 22 8.1053930688133671e-03 -1.2356942914208306e-02 -2.8212751471844815e+00 2.8212599310500082e+00 -1.7442735092468440e-02 3
+P 3 1 11 -8.1053930688133671e-03 1.2356942914208306e-02 -7.1787248397595249e+00 7.1787400689499918e+00 5.1101148051165002e-04 1
+P 4 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 9.9995598131387851e+01 1.0000000000000000e+02 9.3827199999876709e-01 4
+P 5 4 2212 0.0000000000000000e+00 0.0000000000000000e+00 9.9995598131482922e+01 1.0000000000009507e+02 9.3827199999973643e-01 3
+V -3 0 [2,5]
+P 6 -3 553 -2.4074047134142244e-02 2.6879165673024504e-02 5.1101881387075991e+00 1.0752330024333787e+01 9.4602999999952058e+00 2
+P 7 -3 2212 3.2179440203073603e-02 -3.9236108587412617e-02 9.2064134846942238e+01 9.2068929908251562e+01 9.3827200000167510e-01 1
+P 8 6 -13 3.4101173607658830e+00 -1.3038459453885907e+00 5.9422927439678519e+00 6.9750212423676778e+00 1.0565800000106459e-01 1
+P 9 6 13 -3.4341914079000371e+00 1.3307251110616287e+00 -8.3210460525766239e-01 3.7773087819715596e+00 1.0565800000016515e-01 1
+E 0 1 3
+U GEV MM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 1.0000000000000000e+01 1.0000000000000000e+01 0.0000000000000000e+00 4
+P 2 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 9.3799999999999994e-01 9.3799999999999994e-01 4
+V -1 0 [1,2]
+P 3 -1 11 8.3077231638512800e-02 3.4688303864359871e-01 -9.3422177831854492e-01 1.0000001305604915e+00 5.1100000015238706e-04 1
+E 0 1 3
+U GEV MM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 1.0000000000000000e+01 1.0000000000000000e+01 0.0000000000000000e+00 4
+P 2 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 9.3799999999999994e-01 9.3799999999999994e-01 4
+V -1 0 [1,2]
+P 3 -1 11 7.2104820115711854e-01 1.5628742640734802e-01 -1.8589415622754393e+00 2.0000000652802488e+00 5.1100000036965185e-04 1
+E 0 1 3
+U GEV MM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 1.0000000000000000e+01 1.0000000000000000e+01 0.0000000000000000e+00 4
+P 2 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 9.3799999999999994e-01 9.3799999999999994e-01 4
+V -1 0 [1,2]
+P 3 -1 11 -1.4361173664326103e+00 -9.3836236801405426e-01 -2.4611060473140363e+00 3.0000000435201661e+00 5.1100000036965185e-04 1
+E 0 1 3
+U GEV MM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 1.0000000000000000e+01 1.0000000000000000e+01 0.0000000000000000e+00 4
+P 2 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 9.3799999999999994e-01 9.3799999999999994e-01 4
+V -1 0 [1,2]
+P 3 -1 11 -2.6126728458756499e-01 2.0189611701651842e+00 -3.4431867796228350e+00 4.0000000326401244e+00 5.1099999515529708e-04 1
+E 0 1 3
+U GEV MM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 1.0000000000000000e+01 1.0000000000000000e+01 0.0000000000000000e+00 4
+P 2 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 9.3799999999999994e-01 9.3799999999999994e-01 4
+V -1 0 [1,2]
+P 3 -1 11 8.5091281358231141e-01 -2.4919887833464918e+00 -4.2504046027827407e+00 5.0000000261121000e+00 5.1099999863153363e-04 1
+E 0 1 3
+U GEV MM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 1.0000000000000000e+01 1.0000000000000000e+01 0.0000000000000000e+00 4
+P 2 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 9.3799999999999994e-01 9.3799999999999994e-01 4
+V -1 0 [1,2]
+P 3 -1 11 -2.5637625613563038e+00 2.9774447329110156e+00 -4.5345280230083720e+00 6.0000000217600835e+00 5.1100000210777007e-04 1
+E 0 1 3
+U GEV MM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 1.0000000000000000e+01 1.0000000000000000e+01 0.0000000000000000e+00 4
+P 2 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 9.3799999999999994e-01 9.3799999999999994e-01 4
+V -1 0 [1,2]
+P 3 -1 11 8.5248676124549916e-01 1.3423464093622672e+00 -6.8169914507188132e+00 7.0000000186514999e+00 5.1100000210777007e-04 1
+E 0 1 3
+U GEV MM
+P 1 0 11 0.0000000000000000e+00 0.0000000000000000e+00 1.0000000000000000e+01 1.0000000000000000e+01 0.0000000000000000e+00 4
+P 2 0 2212 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 9.3799999999999994e-01 9.3799999999999994e-01 4
+V -1 0 [1,2]
+P 3 -1 11 2.7295844590149682e-01 -1.0660162875541499e+00 -7.9239575315293242e+00 8.0000000163200617e+00 5.1099999515529708e-04 1
+E 0 1 3
+U GEV MM
+HepMC::Asciiv3-END_EVENT_LISTING
+

scripts/plot_ce_ecal_placement.py

+'''
+    A script to calcualte placement of ecal endcap modules
+    lxml is not included in container, get it by simply typing 'pip install lxml'
+    Author: Chao Peng (ANL)
+    Date: 06/17/2021
+'''
+
+import os
+import numpy as np
+import argparse
+import DDG4
+from lxml import etree as ET
+from matplotlib import pyplot as plt
+from matplotlib.collections import PatchCollection
+from matplotlib.patches import Rectangle, Circle
+
+
+CRYSTAL_SIZE = (20., 20., 200.) # mm
+CRYSTAL_GAP = 0.5 # mm
+CRYSTAL_ALIGNMENT = [
+    (5, 21), (5, 21), (5, 21), (4, 22),
+    (3, 23), (0, 26), (0, 24), (0, 24),
+    (0, 24), (0, 24), (0, 24), (0, 24),
+    (0, 22), (0, 22), (0, 20), (0, 20),
+    (0, 18), (0, 18), (0, 16), (0, 16),
+    (0, 14), (0, 14), (0, 12), (0, 12),
+    (0,  6), (0,  6),
+]
+
+GLASS_SIZE = (40., 40., 400.) # mm
+GLASS_GAP = 1.0 # mm
+GLASS_ALIGNMENT = [
+    (13, 10), (13, 10), (13, 10), (12, 10),
+    (12, 10), (12, 10), (11, 11), (10, 11),
+    (9, 12),  (8, 13),  (7, 13),  (6, 14),
+    (3, 16),  (0, 18),  (0, 18),  (0, 16),
+    (0, 16),  (0, 14),  (0, 13),  (0, 11),
+    (0, 10),  (0, 7),   (0, 3),
+]
+
+# calculate positions of modules with a quad-alignment and module size
+def individual_placement(alignment, module_x, module_y, gap=0.):
+    placements = []
+    for row, (start, num) in enumerate(alignment):
+        for col in np.arange(start, start + num):
+            placements.append(((col + 0.5)*(module_y + gap), (row + 0.5)*(module_x + gap)))
+    placements = np.asarray(placements)
+    return np.vstack((placements,
+            np.vstack((placements.T[0]*-1., placements.T[1])).T,
+            np.vstack((placements.T[0], placements.T[1]*-1.)).T,
+            np.vstack((placements.T[0]*-1., placements.T[1]*-1.)).T))
+
+
+def draw_placement(axis, colors=('teal'), module_alignment=((CRYSTAL_SIZE, CRYSTAL_GAP, CRYSTAL_ALIGNMENT))):
+    xmin, ymin, xmax, ymax = 0., 0., 0., 0.
+    patches = []
+    numbers = []
+    for color, (mod_size, mod_gap, alignment) in zip(colors, module_alignment):
+        placements = individual_placement(alignment, *mod_size[:2], mod_gap)
+        boxes = [Rectangle((x - (mod_size[0] + mod_gap)/2., y - (mod_size[1] + mod_gap)/2.), mod_size[0], mod_size[1])
+                 for x, y in placements]
+        patches.append(Rectangle((0., 0.), *mod_size[:2], facecolor=color, alpha=0.5, edgecolor='k'))
+        numbers.append(len(placements))
+        pc = PatchCollection(boxes, facecolor=color, alpha=0.5, edgecolor='k')
+
+        xmin = min(xmin, placements.T[0].min() - 8.*(mod_size[0] + mod_gap))
+        ymin = min(ymin, placements.T[1].min() - 8.*(mod_size[1] + mod_gap))
+        xmax = max(xmax, placements.T[0].max() + 8.*(mod_size[0] + mod_gap))
+        ymax = max(ymax, placements.T[1].max() + 8.*(mod_size[1] + mod_gap))
+
+        # Add collection to axes
+        axis.add_collection(pc)
+    axis.set_xlim(xmin, xmax)
+    axis.set_ylim(ymin, ymax)
+    return axis, patches, numbers
+
+
+def compact_constants(path, names):
+    if not os.path.exists(path):
+        print('Cannot find compact file \"{}\".'.format(path))
+        return []
+    kernel = DDG4.Kernel()
+    description = kernel.detectorDescription()
+    kernel.loadGeometry("file:{}".format(path))
+    try:
+        vals = [description.constantAsDouble(n) for n in names]
+    except:
+        print('Fail to extract values from {}, return empty.'.format(names))
+        vals = []
+    kernel.terminate()
+    return vals
+
+
+if __name__ == '__main__':
+    parser = argparse.ArgumentParser()
+    parser.add_argument('-s', '--save', default='ce_ecal_placement_test.xml',
+            help='path to save compact file.')
+    parser.add_argument('-c', '--compact', default='',
+            help='compact file to get contant to plot')
+    parser.add_argument('--radii-constants', dest='radii', default='EcalBarrel_rmin',
+            help='constant names in compact file to plot, seprate by \",\"')
+    parser.add_argument('--individual', dest='indiv', action='store_true',
+            help='individual block placements instead of line placements')
+    args = parser.parse_args()
+
+    data = ET.Element('lccdd')
+    defines = ET.SubElement(data, 'define')
+
+    # constants: name, value
+    CONSTANTS = [
+        ('CrystalModule_sx', '{:.2f}*mm'.format(CRYSTAL_SIZE[0])),
+        ('CrystalModule_sy', '{:.2f}*mm'.format(CRYSTAL_SIZE[1])),
+        ('CrystalModule_sz', '{:.2f}*mm'.format(CRYSTAL_SIZE[2])),
+        ('CrystalModule_wrap', '{:.2f}*mm'.format(CRYSTAL_GAP)),
+
+        ('GlassModule_sx', '{:.2f}*mm'.format(GLASS_SIZE[0])),
+        ('GlassModule_sy', '{:.2f}*mm'.format(GLASS_SIZE[1])),
+        ('GlassModule_sz', '{:.2f}*mm'.format(GLASS_SIZE[2])),
+        ('GlassModule_wrap', '{:.2f}*mm'.format(GLASS_GAP)),
+
+        ('CrystalModule_z0', '10.*cm'),
+        ('GlassModule_z0', '0.0*cm'),
+        ('EcalEndcapN_z0', '-EcalEndcapN_zmin-max(CrystalModule_sz,GlassModule_sz)/2.'),
+
+        ('CrystalModule_dx', 'CrystalModule_sx + CrystalModule_wrap'),
+        ('CrystalModule_dy', 'CrystalModule_sy + CrystalModule_wrap'),
+        ('GlassModule_dx', 'GlassModule_sx + GlassModule_wrap'),
+        ('GlassModule_dy', 'GlassModule_sy + GlassModule_wrap'),
+    ]
+
+# line-by-line alignment start pos, total number of blocks
+    for name, value in CONSTANTS:
+        constant = ET.SubElement(defines, 'constant')
+        constant.set('name', name)
+        constant.set('value', value)
+
+    # this value will be used multiple times, so define it here
+    readout_name = 'EcalEndcapNHits'
+
+    # detector and its dimension/position/rotation
+    dets = ET.SubElement(data, 'detectors')
+    cmt = ET.SubElement(dets, 'comment')
+    cmt.text = ' Backwards Endcap EM Calorimeter, placements generated by script '
+
+    det = ET.SubElement(dets, 'detector')
+    det.set('id', 'ECalEndcapN_ID')
+    det.set('name', 'EcalEndcapN')
+    det.set('type', 'HomogeneousCalorimeter')
+    det.set('readout', readout_name)
+
+    pos = ET.SubElement(det, 'position')
+    pos.set('x', '0')
+    pos.set('y', '0')
+    pos.set('z', 'EcalEndcapN_z0')
+
+    rot = ET.SubElement(det, 'rotation')
+    rot.set('x', '0')
+    rot.set('y', '0')
+    rot.set('z', '0')
+
+    # placements of modules
+    plm = ET.SubElement(det, 'placements')
+    pltype = 'individuals' if args.indiv else 'lines'
+
+    # crystal
+    crystal = ET.SubElement(plm, pltype)
+    crystal.set('sector', '1')
+    crystal_mod = ET.SubElement(crystal, 'module')
+    crystal_mod.set('sizex', 'CrystalModule_sx')
+    crystal_mod.set('sizey', 'CrystalModule_sy')
+    crystal_mod.set('sizez', 'CrystalModule_sz')
+    crystal_mod.set('material', 'PbWO4')
+    crystal_mod.set('vis', 'AnlTeal')
+    crystal_wrap = ET.SubElement(crystal, 'wrapper')
+    crystal_wrap.set('thickness', 'CrystalModule_wrap')
+    crystal_wrap.set('material', 'Epoxy')
+    crystal_wrap.set('vis', 'WhiteVis')
+    # crystal placements (for individuals)
+    if args.indiv:
+        for m, (x, y) in enumerate(individual_placement(CRYSTAL_ALIGNMENT, *CRYSTAL_SIZE[:2], CRYSTAL_GAP)):
+            module = ET.SubElement(crystal, 'placement')
+            module.set('x', '{:.3f}*mm'.format(x))
+            module.set('y', '{:.3f}*mm'.format(y))
+            module.set('z', 'CrystalModule_z0')
+            module.set('id', '{:d}'.format(m))
+    # crystal placements (for lines)
+    else:
+        crystal.set('mirrorx', 'true')
+        crystal.set('mirrory', 'true')
+        for row, (begin, nmods) in enumerate(CRYSTAL_ALIGNMENT):
+            line = ET.SubElement(crystal, 'line')
+            line.set('axis', 'x')
+            line.set('x', 'CrystalModule_dx/2.')
+            line.set('y', 'CrystalModule_dy*{:d}/2.'.format(row*2 + 1))
+            line.set('z', 'CrystalModule_z0')
+            line.set('begin', '{:d}'.format(begin))
+            line.set('nmods', '{:d}'.format(nmods))
+
+
+    # glass
+    glass = ET.SubElement(plm, pltype)
+    glass.set('sector', '2')
+    glass_mod = ET.SubElement(glass, 'module')
+    glass_mod.set('sizex', 'GlassModule_sx')
+    glass_mod.set('sizey', 'GlassModule_sy')
+    glass_mod.set('sizez', 'GlassModule_sz')
+    # TODO: change glass material
+    glass_mod.set('material', 'PbGlass')
+    glass_mod.set('vis', 'AnlBlue')
+    glass_wrap = ET.SubElement(glass, 'wrapper')
+    glass_wrap.set('thickness', 'GlassModule_wrap')
+    glass_wrap.set('material', 'Epoxy')
+    glass_wrap.set('vis', 'WhiteVis')
+    # crystal placements (for individuals)
+    if args.indiv:
+        for m, (x, y) in enumerate(individual_placement(GLASS_ALIGNMENT, *GLASS_SIZE[:2], GLASS_GAP)):
+            module = ET.SubElement(glass, 'placement')
+            module.set('x', '{:.3f}*mm'.format(x))
+            module.set('y', '{:.3f}*mm'.format(y))
+            module.set('z', 'GlassModule_z0')
+            module.set('id', '{:d}'.format(m))
+    # crystal placements (for lines)
+    else:
+        glass.set('mirrorx', 'true')
+        glass.set('mirrory', 'true')
+        for row, (begin, nmods) in enumerate(GLASS_ALIGNMENT):
+            line = ET.SubElement(glass, 'line')
+            line.set('axis', 'x')
+            line.set('x', 'GlassModule_dx/2.')
+            line.set('y', 'GlassModule_dy*{:d}/2.'.format(row*2 + 1))
+            line.set('z', 'GlassModule_z0')
+            line.set('begin', '{:d}'.format(begin))
+            line.set('nmods', '{:d}'.format(nmods))
+
+
+    # readout
+    readouts = ET.SubElement(data, 'readouts')
+    cmt = ET.SubElement(readouts, 'comment')
+    cmt.text = 'Effectively no segmentation, the segmentation is used to provide cell dimension info'
+    readout = ET.SubElement(readouts, 'readout')
+    readout.set('name', readout_name)
+    seg = ET.SubElement(readout, 'segmentation')
+    # need segmentation to provide cell dimension info
+    # seg.set('type', 'NoSegmentation')
+    seg.set('type', 'MultiSegmentation')
+    seg.set('key', 'sector')
+    crystal_seg = ET.SubElement(seg, 'segmentation')
+    crystal_seg.set('name', 'CrystalSeg')
+    crystal_seg.set('key_value', '1')
+    crystal_seg.set('type', 'CartesianGridXY')
+    crystal_seg.set('grid_size_x', 'CrystalModule_dx')
+    crystal_seg.set('grid_size_y', 'CrystalModule_dy')
+    glass_seg = ET.SubElement(seg, 'segmentation')
+    glass_seg.set('name', 'GlassSeg')
+    glass_seg.set('key_value', '2')
+    glass_seg.set('type', 'CartesianGridXY')
+    glass_seg.set('grid_size_x', 'GlassModule_dx')
+    glass_seg.set('grid_size_y', 'GlassModule_dy')
+    rid = ET.SubElement(readout, 'id')
+    rid.text = 'system:8,sector:4,module:20,x:32:-16,y:-16'
+
+
+    text = ET.tostring(data, pretty_print=True)
+    with open(args.save, 'wb') as f:
+        f.write(text)
+
+
+    fig, ax = plt.subplots(figsize=(12, 12), dpi=160)
+    ax, patches, nblocks = draw_placement(ax, ['teal', 'royalblue'],
+            [(CRYSTAL_SIZE, CRYSTAL_GAP, CRYSTAL_ALIGNMENT), (GLASS_SIZE, GLASS_GAP, GLASS_ALIGNMENT)])
+    ax.set_xlabel('x (mm)', fontsize=24)
+    ax.set_ylabel('y (mm)', fontsize=24)
+    ax.tick_params(direction='in', labelsize=22, which='both')
+    ax.set_axisbelow(True)
+    ax.grid(linestyle=':', which='both')
+    ax.legend(patches, ['{} {}'.format(num, name) for num, name in zip(nblocks, ['PbWO$_4$', 'SciGlass'])], fontsize=24)
+
+    if args.compact and args.radii:
+        names = [c.strip() for c in args.radii.split(',')]
+        radii = compact_constants(args.compact, names)
+        for name, radius in zip(names, radii):
+            ax.add_patch(Circle((0, 0), radius*10., facecolor='none', edgecolor='k', linewidth=2))
+            ax.annotate(name, xy=(radius*10/1.4, radius*10/1.4), fontsize=22)
+    fig.savefig('ce_ecal_placement.png')
+

scripts/run_topside.py → scripts/run_detector_simulation.py

@@ -6,13 +6,13 @@ import logging
 
 import argparse
 parser = argparse.ArgumentParser(
-     prog='run_topside',
+     prog='run_detector_simulation',
      description='''This runs the simulation  the but that is okay''',
      epilog='''
      This program should be run in the "compact" directory.
          ''')
 parser.add_argument("-v","--verbose", help="increase output verbosity", type=int, default=0)
-parser.add_argument("--compact", help="compact detector file",default="topside.xml")
+parser.add_argument("--compact", help="compact detector file",default="athena.xml")
 parser.add_argument("--vis", help="vis true/false", action="store_true",default=False)
 parser.add_argument("--ui", help="ui setting tcsh or qt; default=qt", type=str,default="qt",dest="ui")
 parser.add_argument("-b","--batch", help="batch turns off vis/ui", action="store_true",default=False, dest="batch")
@@ -86,7 +86,7 @@ def run():
 
     outputfile = args.output
     if outputfile is None:
-        outputfile = 'data/topside_' + time.strftime('%Y-%m-%d_%H-%M')
+        outputfile = 'data/athena_' + time.strftime('%Y-%m-%d_%H-%M')
     #rootoutput = geant4.setupROOTOutput('RootOutput', outputfile)
     #rootoutput.HandleMCTruth = True
 
@@ -214,21 +214,21 @@ def run():
     #seq, act = geant4.setupDetector('HeavyGasCherenkov','PhotoMultiplierSDAction')
     #act.adopt(f3)
 
-    seq, act = geant4.setupTracker('SiTrackerBarrel_Layer1')
-    seq, act = geant4.setupTracker('SiTrackerEndcapP_Layer1')
-    seq, act = geant4.setupTracker('SiTrackerEndcapN_Layer1')
-    seq, act = geant4.setupTracker('SiTrackerBarrel_Layer2')
-    seq, act = geant4.setupTracker('SiTrackerEndcapP_Layer2')
-    seq, act = geant4.setupTracker('SiTrackerEndcapN_Layer2')
-    seq, act = geant4.setupTracker('SiTrackerBarrel_Layer3')
-    seq, act = geant4.setupTracker('SiTrackerEndcapP_Layer3')
-    seq, act = geant4.setupTracker('SiTrackerEndcapN_Layer3')
-    seq, act = geant4.setupTracker('SiTrackerBarrel_Layer4')
-    seq, act = geant4.setupTracker('SiTrackerEndcapP_Layer4')
-    seq, act = geant4.setupTracker('SiTrackerEndcapN_Layer4')
-    seq, act = geant4.setupTracker('SiTrackerBarrel_Layer5')
-    seq, act = geant4.setupTracker('SiTrackerEndcapP_Layer5')
-    seq, act = geant4.setupTracker('SiTrackerEndcapN_Layer5')
+    #seq, act = geant4.setupTracker('SiTrackerBarrel_Layer1')
+    #seq, act = geant4.setupTracker('SiTrackerEndcapP_Layer1')
+    #seq, act = geant4.setupTracker('SiTrackerEndcapN_Layer1')
+    #seq, act = geant4.setupTracker('SiTrackerBarrel_Layer2')
+    #seq, act = geant4.setupTracker('SiTrackerEndcapP_Layer2')
+    #seq, act = geant4.setupTracker('SiTrackerEndcapN_Layer2')
+    #seq, act = geant4.setupTracker('SiTrackerBarrel_Layer3')
+    #seq, act = geant4.setupTracker('SiTrackerEndcapP_Layer3')
+    #seq, act = geant4.setupTracker('SiTrackerEndcapN_Layer3')
+    #seq, act = geant4.setupTracker('SiTrackerBarrel_Layer4')
+    #seq, act = geant4.setupTracker('SiTrackerEndcapP_Layer4')
+    #seq, act = geant4.setupTracker('SiTrackerEndcapN_Layer4')
+    #seq, act = geant4.setupTracker('SiTrackerBarrel_Layer5')
+    #seq, act = geant4.setupTracker('SiTrackerEndcapP_Layer5')
+    #seq, act = geant4.setupTracker('SiTrackerEndcapN_Layer5')
     #seq, act = geant4.setupTracker('SiVertexBarrel')
     #seq, act = geant4.setupTracker('SiTrackerForward')
     #seq, act = geant4.setupCalorimeter('EcalBarrel')

scripts/run_simulation.sh

-#!/bin/bash
-
-npsim --runType batch --numberOfEvents 100000 \
-      --compactFile ./topside.xml \
-      --inputFiles  ../datasets/data/emcal_electron_0GeVto30GeV_100kEvt.hepmc \
-      --outputFile  ./sim_electron_0GeVto30GeV_100k_input.root

scripts/test_ACTS.cxx

+R__LOAD_LIBRARY(libDDCore.so)
+R__LOAD_LIBRARY(libActsPluginDD4hep.so)
+R__LOAD_LIBRARY(libDDG4.so)
+R__LOAD_LIBRARY(libDDG4IO.so)
+#include "DD4hep/Detector.h"
+#include "DD4hep/DetElement.h"
+#include "DD4hep/Objects.h"
+#include "DD4hep/Detector.h"
+#include "DDG4/Geant4Data.h"
+#include "DDRec/CellIDPositionConverter.h"
+#include "DDRec/SurfaceManager.h"
+#include "DDRec/Surface.h"
+
+#include "TCanvas.h"
+#include "TChain.h"
+
+#include "Acts/Geometry/TrackingGeometry.hpp"
+#include "Acts/Geometry/TrackingVolume.hpp"
+#include "Acts/Plugins/DD4hep/ConvertDD4hepDetector.hpp"
+
+/** Example loading ACTs.
+ *
+ *
+ */
+void test_ACTS(const char* compact = "athena.xml"){
+
+  using namespace ROOT::Math;
+  // -------------------------
+  // Get the DD4hep instance
+  // Load the compact XML file
+  // Initialize the position converter tool
+  dd4hep::Detector& detector = dd4hep::Detector::getInstance();
+  detector.fromCompact(compact);
+  dd4hep::rec::CellIDPositionConverter cellid_converter(detector);
+
+  //std::unique_ptr<const Acts::TrackingGeometry> 
+  auto acts_tracking_geometry  = Acts::convertDD4hepDetector (detector.world(),Acts::Logging::Level::VERBOSE);
+  //acts_tracking_geometry  = Acts::convertDD4hepDetector (detector.world(),Acts::Logging::Level::INFO);
+
+  if(acts_tracking_geometry) {
+    std::cout << "success?\n";
+  }
+  //  if(acts_tracking_geometry->highestTrackingVolume()) {
+  //    std::cout << " volume name \n ";
+  //    std::cout << acts_tracking_geometry->highestTrackingVolume()->volumeName() << std::endl;
+  //  } else {
+  //    std::cout << "derp\n";
+  //  }
+  //}
+
+}

scripts/view1/.DAWN_1.history

+1.34392e+07
+90
+180
+0
+0
+0
+0
+1
+1
+0.001
+0
+1
+1
+1
+0.5
+0.5
+0.5
+19
+71
+0.001
+0.001
+0.001
+3
+71
+0.001
+0
+0
+1
+evince
+0
+0

scripts/view1/generate_eps

@@ -9,8 +9,8 @@ function print_the_help {
   exit 
 }
 
-FILE_TAG="view1"
-INPUT_FILE="g4_0000.prim"
+FILE_TAG="view01"
+INPUT_FILE="../../g4_0000.prim"
 
 
 POSITIONAL=()
@@ -44,14 +44,30 @@ done
 set -- "${POSITIONAL[@]}" # restore positional parameters
 
 
-dawncut 1 0 0 1 ${INPUT_FILE} ${FILE_TAG}_temp0.prim 
-dawncut -1 0 0 1 ${FILE_TAG}_temp0.prim  ${FILE_TAG}.prim
-dawn -d ${FILE_TAG}.prim 
+# Side view (lines)
+dawncut -1 0 0 1 ${INPUT_FILE} ${FILE_TAG}.prim
+# dawncut 1 0 0 1 ${INPUT_FILE} ${FILE_TAG}_temp0.prim 
+# dawncut -1 0 0 1 ${FILE_TAG}_temp0.prim  ${FILE_TAG}.prim
+../../bin/dawn_tweak --mag 10 --draw 3 --theta 165 --phi 75 --light-theta 180 --light-phi 90
+
+dawn -d ${FILE_TAG}.prim
 ps2pdf ${FILE_TAG}.eps ${FILE_TAG}_full.pdf
 gs -o ${FILE_TAG}.pdf -sDEVICE=pdfwrite \
   -c "[/CropBox [51 250 550 590] /PAGES pdfmark" \
   -f ${FILE_TAG}_full.pdf
 
-pdftoppm ${FILE_TAG}.pdf ${FILE_TAG} -png -singlefile -cropbox
+pdftoppm ${FILE_TAG}.pdf ${FILE_TAG} -png -singlefile -cropbox -thinlinemode solid -aaVector yes -r 600
+
+# Top view
+#dawncut 0 -1 0 1 ${INPUT_FILE} ${FILE_TAG}.prim 
+dawncut 0 1 0 1 ${INPUT_FILE} ${FILE_TAG}_temp0.prim 
+dawncut 0 -1 0 1 ${FILE_TAG}_temp0.prim  ${FILE_TAG}.prim
+../../bin/dawn_tweak --mag 10 --draw 1 --theta 90 --phi 90
+dawn -d ${FILE_TAG}.prim 
+ps2pdf ${FILE_TAG}.eps ${FILE_TAG}_top_full.pdf
+gs -o ${FILE_TAG}_top.pdf -sDEVICE=pdfwrite \
+  -c "[/CropBox [51 250 550 590] /PAGES pdfmark" \
+  -f ${FILE_TAG}_top_full.pdf
 
+pdftoppm ${FILE_TAG}_top.pdf ${FILE_TAG}_top -png -singlefile -cropbox -thinlinemode solid -aaVector yes
 

scripts/view11/.DAWN_1.history

+1.34392e+07
+90
+180
+0
+0
+0
+25000
+2.0
+1
+0.001
+0
+1
+1
+1
+0.5
+0.5
+0.5
+19
+71
+0.001
+0.001
+0.001
+3
+71
+0.001
+0
+0
+1
+evince
+0
+0

scripts/view12/.DAWN_1.history

+1.34392e+07
+0
+180
+0
+0
+0
+0
+8
+1
+0.001
+0
+1
+1
+1
+0.5
+0.5
+0.5
+19
+71
+0.01
+0.01
+0.01
+3
+70
+0.01
+1
+1
+1
+evince
+0
+0

scripts/view12/generate_eps

@@ -18,7 +18,7 @@ tan ()
 
 
 function print_the_help {
-  echo "USAGE: $0 -i <PRIM_FILE>  "
+  echo "USAGE: $0 -i <PRIM_FILE> <slices ...> "
   echo "  OPTIONS: "
   echo "            -t,--tag           filename tag (default: view1)"
   exit 
@@ -48,12 +48,16 @@ do
       shift # past argument
       shift # past value
       ;;
-    *)    # unknown option
-      #POSITIONAL+=("$1") # save it in an array for later
+    -[a-zA-Z]*) # unknown option
+      POSITIONAL+=("$1") # save it in an array for later
       echo "unknown option $1"
       print_the_help
       shift # past argument
       ;;
+    *)     # positional options
+      POSITIONAL+=("$1") # save it in an array for later
+      shift # past argument
+      ;;
   esac
 done
 set -- "${POSITIONAL[@]}" # restore positional parameters
@@ -62,17 +66,6 @@ set -- "${POSITIONAL[@]}" # restore positional parameters
 
 echo "view12 produces a series of slightly rotated XY slices a different z locations. Along beamline"
 
-# slice at z = 2m
-# note the offset has to change with sign of the direction to cut in the opposite direction.
-dawncut 0  0  1  2005 ${INPUT_FILE} ${FILE_TAG}b_temp0.prim 
-dawncut 0  0 -1 -2000 ${FILE_TAG}b_temp0.prim  ${FILE_TAG}b.prim
-dawn -d ${FILE_TAG}b.prim 
-ps2pdf ${FILE_TAG}b.eps ${FILE_TAG}b_full.pdf
-gs -o ${FILE_TAG}b.pdf -sDEVICE=pdfwrite \
-  -c "[/CropBox [50 175 550 675] /PAGES pdfmark" \
-  -f ${FILE_TAG}b_full.pdf
-pdftoppm ${FILE_TAG}b.pdf ${FILE_TAG}b -png -singlefile -cropbox
-
 original_file_tag="${FILE_TAG}"
 
 make_slice(){
@@ -90,9 +83,10 @@ make_slice(){
   rm "${FILE_TAG}_temp0.prim"
   rm "${FILE_TAG}.prim"
 }
-for zzz in $(seq 50 50 2000) ;
+
+for zzz in $@ ;
 do
-  make_slice ${zzz} & 
+  make_slice ${zzz}
 done
 
 wait