1 / 13

Simulation

GlueX collaboration mtg Bloomington, May 20-22, 2004. R.T. Jones, UConn. Simulation. Monte Carlo. detector simulation. analysis. Detector simulation tools. MCFast reliable for acceptance, resolution sufficient for many detector design decisions, PWA studies 10 3 events/cpu/s

delu
Download Presentation

Simulation

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. GlueX collaboration mtg Bloomington, May 20-22, 2004 R.T. Jones, UConn Simulation Monte Carlo detector simulation analysis

  2. Detector simulation tools • MCFast • reliable for acceptance, resolution • sufficient for many detector design decisions, PWA studies • 103 events/cpu/s • in production • Geant 3 • reliable for e-m interactions, backgrounds • needed for development of reconstruction software • 10 events/cpu/s • in production standard interfaces need are needed

  3. Detector simulation interfaces Monte Carlo geometry Simulation engine hits

  4. Detector simulation interfaces geometry • Decided collaboration mtg 3/2001: • standard geometry interface is needed • should be expressed in xml • Completed 5/2001: version 1.0 • entire detector exists in prototype • interface to MCfast completed and tested • interface to Geant in completed and tested • project web site http://zeus.phys.uconn.edu/halld/geometry

  5. Detector simulation interfaces • Developments since Dec. 2003: • creation of schema for geometry package • upgrade to new XERCES-C v2 from Apache.org • new visualization tool hdds-root by Ed Brash • refinements to geometry description in several areas by Ed Brash, Curtis Meyer • overhaul of hdds-mcfast tool for HDFast db geometry • Plans for next six months: • serious work on understanding backgrounds • user support – please use the online forums!

  6. Detector simulation interfaces • Decided at workfest 5/2001: • standard data model is needed • should be expressed in xml Monte Carlo hits • Completed 6/2001: version 1.0 • data model for Monte Carlo generation • interface to genr8 completed and tested • interface to Mcfast completed • interface to Geant completed • project web site http://zeus.phys.uconn.edu/halld/datamodel/doc

  7. Detector simulation interfaces • Work since Dec. 2003: • upgrade to XERCES-C v2 • introduction of schema-based spec • new tools to support schemas • hddm-schema • schema-hddm • new general-purpose tools • hddmcat • xml-xml Monte Carlo hits • Plans for the next six months • begin draft data model for reconstruction • improved hddm library C++ interface

  8. Monte Carlo Data Acquisition digitized signals in buffers reaction specification generation event building final state raw event record simulation conversion hits geometric reconstruction tracks /clusters kinematic reconstruction final state http://zeus.phys.uconn.edu/halld/datamodel/doc

  9. Data Model: Template example #1: define hits in forward TOF <forwardTOF> <slab y=“float” minOccurs=“0” maxOccurs=“unbounded”> <side end=“int” minOccurs=“0” maxOccurs=“unbounded”> <hit t=“float” dE=“float” minOccurs=“0” maxOccurs=“unbounded” /> </side> </slab> </forwardTOF> http://zeus.phys.uconn.edu/halld/datamodel/doc

  10. Data Model: Example #1: TOF hits structure in memory typedef struct { t_ForwardTOF_t* forwardTOF; } t_HDDM_t; typedef struct { t_Slabs_t* slabs; } t_ForwardTOF_t; typedef struct { float y; t_Sides_t* sides; } t_Slab_t; typedef struct { int end; t_Hits_t* hits; } t_Side_t; typedef struct { float t; float dE; } t_Hit_t; http://zeus.phys.uconn.edu/halld/datamodel/doc

  11. Data Model: Example #2: define Monte Carlo event <HDDM class="s" version="1.0"> <physicsEvent eventNo="int" runNo="int"> <reaction type="int" weight="float" minOccurs=“1” maxOccurs=“unbounded”> <beam type="Particle_t"> <momentum px="float" py="float" pz="float" E="float" /> <properties charge="int" mass="float" /> </beam> <target type="Particle_t"> <momentum px="float" py="float" pz="float" E="float" /> <properties charge="int" mass="float" /> </target> <vertex minOccurs=“1” maxOccurs=“unbounded”> <product type="Particle_t" decayVertex="int" minOccurs=“1” maxOccurs=“unbounded”> <momentum px="float" py="float" pz="float" E="float" /> <properties charge="int" mass="float" /> </product> <origin vx="float" vy="float" vz="float" t="float" /> </vertex> </reaction> </physicsEvent> </HDDM> http://zeus.phys.uconn.edu/halld/datamodel/doc

  12. Data Model: Existing HDDM tools: • hddm-c: translates hddm data model into c • hddm-xml: converts hddm stream into xml listing • xml-hddm: converts xml listing to hddm stream • stdhep-hddm: converts stdhep files into hddm stream • hddm-schema: reads hddm template, writes schema • schema-hddm: reads schema, writes hddm template • hddmcat: concantenate multiple hddm streams • xml-xml: pretty-print free-format xml for viewing http://zeus.phys.uconn.edu/halld/datamodel/doc

  13. Data Model: Example: output from genr8 MC generator processed by stdhep-hddm <?xml version="1.0" encoding="UTF-8" ?> <HDDM class="s" version="1.0” version=“1.0” xmlns=“http://www.gluex.org/hddm”> <physicsEvent runNo="-9000" eventNo="1"> <reaction type="0" weight="0.000000"> <vertex> <product type="pi-" decayVertex="0"> <momentum E="5.937384" px="-0.197764" py="0.586868" pz="5.903338" /> <properties mass="0.140000" charge="-1" /> </product> <product type="pi+" decayVertex="0"> <momentum E="1.947875" px="0.001550" py="-0.182470" pz="1.934249" /> <properties mass="0.140000" charge="1" /> </product> <product type="proton" decayVertex="0"> <momentum E="1.052739" px="0.196214" py="-0.404399" pz="0.162411" /> <properties mass="0.938000" charge="1" /> </product> <origin t="0.000000" vx="0.000000" vy="0.000000" vz="0.000000" /> </vertex> </reaction> </physicsEvent> <physicsEvent runNo="-9000" eventNo="2"> . . . http://zeus.phys.uconn.edu/halld/datamodel/doc

More Related