Common Tools. G.Gavalian (JLAB)

Transcription

1 Common Tools G.Gavalian (JLAB) 1

2 Software components Data Input/Output EVIO data structure definitions. JAVA libraries for dictionary based I/O. Geometry Geometry primitives library. Abstract classes for detector description. Detector implementations. Reconstruction Framework Implementation of abstract classes for reconstruction. Unified CLARA interface. Scripting environment. Histograming Java Histograming package (plot/fit/export) Histogram I/O from JAVA Histogram ROOT import 2

3 EVIO Input/Output Standardized XML dictionaries for EVIO bank structures: Each detector component has an XML file describing data structures. Only requested banks are deserialized. Reading BOS (legacy) files is transparent to the program. 3

4 Geometry Point3D(double x, double y, double z) point in 3d space. convert to vector. (x,y,z) distance from given point. Line3D(Point3D origin, Point3D end) (x,y,z)1 line connecting two points. can be treated as ray or infinite line. line direction, line length. (x,y,z)0 Face3D(Point3D p0, Point3D p1, Point3D p2) P0 plane with three points. normal to the plane. P1 intersect with line. P2 4

5 Geometry shapes n Plane3D(Point3D origin, Vector3D dir) origin point and normal vector intersection with lines and segments intersection with planes Shape3D() collection of face 3D objects intersection with lines intersection with planes Arc3D(Point3D center, Vector3D A, Vector3D B) h origin, opening angle and normal. intersection with lines and segments. B n A 5

6 Geometry primitives hierarchy Point3D Vector3D Line3D Arc3D Face3D Plane3D Sector3D Cylindrical3D Path3D Shape3D Fundamental Derived Surface Container 6

7 Detector Hierarchy Properties Methods Detector DetectorId getdetectorid() String gettype() void show() int Sector List<Sector> getnumsectors() getsector(int sect) getallsectors() Sector int String void getsectorid() gettype() show() int getnumsuperlayers() Superlayer getsuperlayer(int supl) List<Superlayer> getallsuperlayers() Superlayer int getsuperlayerid() String gettype() void show() int Layer List<Layer> getnumlayers() getlayer(int layer) getalllayers() Layer int String void getlayerid() gettype() show() int getnumcomponents() Component getcomponent(int c) List<Component> getallsectors() Component int String void getcomponentid() gettype() show() Point3D double double getmidpoint() getlength() getwidth() 7

8 Detector Hierarchy FTOFDetector ECDetector DCDetector FTOFSector ECSector DCSector FTOFSuperlayer ECSuperlayer DCSuperlayer FTOFLayer ECLayer DCLayer Scintillator Paddle Prismatic Component DriftChamber Wire Shape3D 8

9 Detector: Usage ConstantProvider dcconstants = DataBaseLoader.getConstantsDC(); DCDetector dc = new DCFactory().createDetectorCLAS(dcConstants); DCSector sector = dc.getsector(0); DCSuperlayer superlayer = sector.getsuperlayer(5); DCLayer layer = superlayer.getlayer(5); DriftChamberWire wire = layer.getcomponent(12); Point3D wirestart = wire.getline().origin(); double length = dc.getsector(0).getsuperlayer(0).getlayer(0).getcomponent(0).getlength(); int numsectors = dc.getnumsectors(); List<Point3D> midpoints = new ArrayList(); for (DCSector sector : dc.getallsectors()) for (DCSuperlayer superlayer : sector.getallsuperlayers()) for (DCLayer layer : sup.getalllayers()) for (DriftChamberWire wire : layer.getallcomponents()) midpoints.add(wire.getmidpoint()); 9

10 Detector: Usage ConstantProvider constants = DataBaseLoader.getConstantsEC(); ECDetector detector; detector = new ECFactory().createDetectorCLAS(constants); detector = new ECFactory().createDetectorSector(constants); detector = new ECFactory().createDetectorTilted(constants); detector = new ECFactory().createDetectorLocal(constants); // Constructed in local coordinates (varies by detector) ECSector firstsector = factory.createsector(constants, 0); ECSuperlayer outerec = factory.createsuperlayer(constants, 0, 2); ECLayer wlayer = factory.createlayer(constants, 0, 2, 2); // Factories have no member variables or initialization procedure. // It is acceptable to instantiate a new factory every time one is needed. // However some detectors DO have non-negligible initialization procedures. // For example, the BST has 33,792(!) wires, the end points of which are // calculated using 67,584 line-plane intersections. This takes ~0.5 seconds. 10

11 Detector: coordinate Systems Local Tilted 11

12 Detector: coordinate Systems Sector 12

13 Detectors: coordinate System (CLAS) 13

14 Detector: Visualizer datadc = DataBaseLoader.getConstantsDC() factory = DCFactory() dcdetector = factory.createdetectorclas(datadc) dataec = DataBaseLoader.getConstantsEC() factory = ECFactory() ecdetector = factory.createdetectorclas(dataec) datasc = DataBaseLoader.getConstantsFTOF() factory = FTOFFactory() scdetector = factory.createdetectorclas(datasc) path = swimmer.getparticlepath(particle) hitsdc = dcdetector.gethits(path) hitsec = ecdetector.gethits(path) hitssc = scdetector.gethits(path) for hit in hitsdc: print '\t',hit.tostring() for hit in hitssc: print '\t',hit.tostring() for hit in hitsec: print '\t',hit.tostring() HIT : DC ( , ) XYZ = HIT : DC ( , ) XYZ = HIT : DC ( , ) XYZ = HIT : DC ( , ) XYZ = HIT : DC ( , ) XYZ = HIT : DC ( , ) XYZ = HIT : FTOF ( , ) XYZ = HIT : FTOF ( , ) XYZ = HIT : EC ( , ) XYZ = HIT : EC ( , ) XYZ = HIT : EC ( , ) XYZ = HIT : EC ( , ) XYZ = HIT : EC ( , ) XYZ = HIT : EC ( , ) XYZ =

15 Electron Acceptance Graph2D Graph2D Graph2D DVCS Generator

16 Fast Monte-Carlo Detector Response tracker. Fast monte-carlo module based on geometry library. Uses common magnetic swimmer. Very fast fiducial region simulation. Reconstruction: Combined generator and reconstruction module. Particle momentum smearing. Missing wire implementation. Dead region implementation. Analysis: Produced output is the same as for real reconstruction. 16

17 Reconstruction Suite Development of CLAS12 reconstruction software is implemented in JAVA. A common development framework was developed and packaged, and is available on github. The package comes with JYTHON and Groovy support, with scripting support. Package contains all necessary software for: EVIO input/output Embedded CLARA functionality CLAS12 detector geometry (ccdb sqlite is available) Plotting and fitting libraries (minuit) Plans to implement common validation framework. 17

18 Reconstruction BEFORE AFTER Reconstruction DC EVIO ALGORI THM Abstract Class GEOM GEOM CALIB CCDB CLARA CLARA EVIO FTOF EVIO ALGORI THM DC FTOF EC GEOM CLARA CCDB Significantly reduced the code base All major functionality is maintained in one class Easy to fix bugs and add features 18

19 Reconstruction Framework Generators: Physics process event generators (Pythia 6,HERWIG++) JAVA libraries for event selection. Simulation: GEMC 2.1 géant simulation. LUND input files. Reconstruction: JAVA based standalone reconstruction. Precompiled, included in COATJAVA package Service oriented reconstruction (multithreaded,in progress) Data Analysis: High Level DST readers, analysis visualization tools Interface to ROOT. 19

20 Event Generators Pythia 6 FORTRAN based event generator. Generated event files stored on JLAB network. Selecting events from generation: Usage : python PhysicsFilter.py -f filename.data -p [particle list ":" separated] [flags] Options: -f file : file name to filter -d dir : directory name where lund files are located (*.data files will be processed) -p list : particle list ":" separated (i.e. 11:2212:211:-211) -e : switch on exlusive flag (default is inclusive) Examples: To select events of inclusive pi0 production on proton use: >python PhysicsFilter.py -f file.data -p 11:2212:22:22 To select events of exclusive phi meson production on proton use: >python PhysicsFilter.py -f file.data -p 11:2212:321:-321 -e Event Generator output will be available through CLARA. 20

21 Reconstruction Reconstruction libraries for several detectors are implemented. Common suite can run in single threaded stand alone mode. Ability to dynamically change the reconstruction flow. Usage : runreconstruction inputfile.evio -[flags] Options: -o file : output file name -n num : number of events to reconstruct -e evt : reconstruct one event with EVT NUM = evt (bank 10 evio) -s LIST : list of reconstruction Modules to run (ALL=DC:FTOF:EC:HTCC) Flags: -r : run reconstruction with reversed field Examples: To run DC and FTOF reconstruction for 300 events on generated file use: >runreconstruction input.evio -s DC:FTOF -n

22 Reconstruction: Scripting Groovy scripting support Familiar JAVA syntax Interpreted at run time Can be compiled to distribute Full integration with CLAS tools Can run in CLARA Plotting/Fitting library included Easy prototyping 22

23 Histograming Histograming package is implemented in JAVA for reconstruction suite: H1D, H2D and H3D with output to EVIO Plotting, Fitting Libraries in Java (using jminuit) ROOT tools to read EVIO histograms 23

24 ROOT Tools Motivation for using ROOT as an analysis tool: Most of high analysis and prototyping are done in ROOT. Most users are familiar with ROOT. There is extensive documentation online. Common approach incorporating ROOT into analysis chain: Convert existing file format to ROOT tuples. Use custom format in ROOT. What we end up with: Several different versions of DSTs (redundant storage). Versions of analysis code that can not be used by other groups. 24

25 ROOT Tree TFile LoadTree TTree Branch Tree Player Leaf Leaf Selector Leaf Draw( x, x<1.0&&y>2.0 ) 25

26 ROOT EVIO Tree TFile Evio TTree Evio LoadTree Branch Tree Player Leaf Leaf Selector Leaf Draw( x, x<1.0&&y>2.0 ) 26

27 ROOT/EVIO Tree Browser root[10] tree->draw( FTOF1B::dgtz.paddles, FTOF1B::dgtz.sector==1 ); 27

28 Reading ROOT/EVIO Tree in CINT 28

29 Conclusion Framework: Existing software components are using the common-tools library. Unified I/O and Geometry Detector Geometry: Geometry primitives library Detector factories Detector initialization from mysql and sqlite Reconstruction: CLARA encapsulation Abstract classes to configure services and load geometry Reconstruction diagnostics Scripting language support Analysis: Java data analysis package (histograming/fitting) Unified physics analysis framework (includes corrections and cuts) ROOT evio Reader ROOT evio histogram reader 29

30 BACKUP-SLIDES 30

31 Conclusion Data Input/Output: Standard library for I/O in Java and ROOT. XML dictionaries for bank structures. Detector Geometry: Geometry primitives library Detector factories Detector initialization from mysql and sqlite Reconstruction: CLARA encapsulation Abstract classes to configure services and load geometry Reconstruction diagnostics Scripting language support Analysis: Java data analysis package (histograming/fitting) ROOT evio Reader ROOT evio histogram reader 31