Tracking, Computing & other Stuff

1. Tracking, Computing & other Stuff

2. Correlation of detector hits • The track segments of inner and outer MDCs are matched on Cluster level • Outer segments are matched with META hits • Inner segments are matched to RICH hits MAGNETIC FIELD BEAM RICH MDCI+II MDCIII+IV META

3. Tracking in MDCs • step 1 : projection of drift cells with respect to the target provides Cluster in inner segment • step 2 : projection of outer drift cells with respect to hit point of inner segment on the kickplane defines Cluster in outer segment • step 3 : Fitting of inner and outer Cluster by a straight line model function using the measured drift times (Least square minimization of the functional below)

4. Correlation of detector hits • HMetaMatchF task matches the detector hits of RICH,MDC and META • The matched objects provide input to the momentum reconstruction algorithms • Rings in the RICH and track segments in the inner MDCs are spatial correlated inside a matching window given by • Outer segments are matched with META hits using straight lines (neglecting the curvature of the track in magnetic field). From the difference of the propagated hit position of the segment on the META and the measured Hits position a quality factor (normalized by the errors) is calculated.

5. MAGNETIC FIELD BEAM MDCI+II META Momentum reconstruction Kickplane • KickPlane: • use inner segment + META hit to define the deflection angle and momentum • low resolution, very fast • Spline : • use 5th order spline to emulate a track in the field region • use inner and outer segments • the segment directions are not touched by the algorithm • good resolution, fast • Runge Kutta: • propagates particle through magnetic field • uses inner and outer segments • modifies the segments direction and position • best resolution, computing intensive MAGNETIC FIELD Spline BEAM MDCI+II MDCIII+IV META MAGNETIC FIELD Runge Kutta BEAM MDCI+II MDCIII+IV META MDCI+II MDCIII+IV META

6. HParticleTrackSorter • After the tracking procedure (hit reconstruction, Hit matching, momentum reconstruction still tracks are inside the sample of candidates which share the same detector hits • To get rid of such candidates a cleaning procedure is applied • All track candidates are inspected and flagged according certain properties • Two select functions (selectLeptons,selectHadrons) are performed to apply conditions on the candidates • Candidates which have passed the select functions are inserted in a list which is sorted by track quality • Starting from the best candidate a list of used detector hits is build and each following candidate which reuses an hit is discarded • Candidates which survives the selecting are marked with kIsUsed • Candidates selected by the selectLepton function are marked in addition as kIsLepton

7. HParticleTrackSorter

8. Helper functions static TObjArray* HParticleTool::getMdcCal1Cluster(Int_t segind) static TObjArray* HParticleTool::getMdcCal1Seg(Int_t segind) static HMdcClus* HParticleTool::getMdcClus(Int_t segind) static HMdcClusFit* HParticleTool::getMdcClusFit(Int_t segind) static HMdcClusInf* HParticleTool::getMdcClusInf(Int_t segind, Int_t nhit = 0) static HMdcHit* HParticleTool::getMdcHit(Int_t segind, Int_t nhit = 0) static HMdcSeg* HParticleTool::getMdcSeg(Int_t segind) static HMdcTrkCand* HParticleTool::getMdcTrkCand(Int_t metaind) static TObjArray* HParticleTool::getMdcWireFitSeg(Int_t segind) static HMetaMatch2* HParticleTool::getMetaMatch(Int_t metaind) static HRichHit* HParticleTool::getRichHit(Int_t richind) static HRpcCluster* HParticleTool::getRpcCluster(Int_t rpcind) static HShowerHit* HParticleTool::getShowerHit(Int_t showerind) static HTofCluster* HParticleTool::getTofCluster(Int_t tofind) static HTofHit* HParticleTool::getTofHit(Int_t tofind)

9. Things which are obvious how to use ... HMdcSeg HMdcHit HMdcClusInf HMdcClusFit HMdcWireFit TObjArray* HParticleTool::getMdcWireFitSeg(Int_t segind)

10. New Batch farm • In spring 2012 a new batch farm will go into operation (see Walter’s talk) • The new farm will have new features • The operating system will be 64bit • 32bit execution is not supported any more • No user file systems will be mounted, only /lustre  stability reason • The user have to copy their scripts , parameter files etc to /lustre before submitting the job • The software will be distributed experiment specific using CVMFS  avoid bottle neck of parallel access • LSF will be replaced by Sun grid engine (SGE): batch scripts have to be adapted

11. New Batch farm publish lxbuild02.gsi.de /cvmfs/hades.gsi.de • Installation of the HADES software: • The software is build on lxbuild02.gsi.de and stored local /cvmfs/hades.gsi.de • After installation the software has to be published to enable the user access it • The publish command runs basically a rsync to the cvmfs server • From the server the software will be distributed to all hosts and seen as /cvmfs/hades.gsi.de Cvmfs server lxb320.gsi.de lxb321.gsi.de distribute lxb322.gsi.de lxb323.gsi.de lxb324.gsi.de

12. Old Batch farm • Debian etch (current system) • Very old • Out of support  security issue • switched off at GSI January 2012 • Our situation: • Hydra2/HGeant2 will run on any new system • Old analyses of data taken 2005 and before requires old software versions • Solution: • Old batch farm will be switched to lenny64 • Old binaries (32bit compiled) will run without change • Programs can be recompiled under lenny32 (lennylust32.gsi.de) • Tests have been successful  we are ready to switch

13. Extended Event Display • Based on TEve of ROOT • OpenGL • Display Detector geometry taken from ROOT TGeomModeller • Helper classes implemented in libEventdisplay • User can easily modify displayed content by changing macros • Full event display for REAL and SIM Data • Display of HGEANT particles • Runge-Kutta tracking trough Detector using a magnetic Field map. Fit to reference points given by HGeant • TODO: • Pluto particles • ??? Requests welcome

14. As easy as this …edit nextEvent.C HParticleCandSim* cand; Int_t size = particleCandCat->getEntries(); for(Int_t i = 0; i < size; i ++){ cand = HCategoryManager::getObject(cand,particleCandCat,i); if(cand){ HEDParticleCand* edcand = 0; Int_t s = cand->getSector(); // Fill all objects. Group them into different lists // Groups get different colors or line styles //--------------------------------------------------------- // Candidates accepted as leptons by HParticleTrackSorter // Full tracks with RICH and accepted after sorting if(cand->isFlagBit(kIsLepton)){ edcand = new HEDParticleCand(cand); edcand->SetLineColor(kRed); particlecandLep->AddElement(s,edcand); } //--------------------------------------------------------- //--------------------------------------------------------- // Ghost tracks get dashed lines if(edcand && cand->isGhostTrack()) edcand->SetLineStyle(7); //--------------------------------------------------------- } } // end loop particlecand

18. Geometry display options can be modified via GUI • The event scene of the event display is defined inside a macro an be easily adopted to the user needs by non experts

19. Geometry display options can be modified via GUI • The event scene of the event display is defined inside a macro an be easily adopted to the user needs by non experts

20. Thanks to the flexible event scene the user can easily modify the displayed contents

21. Important Infos about LSF look to the man pages of: bsub bjobs bjobs -u all bjobs -u all -q hades bjobs -r bjobs -w bpeek bkill bkill id bkill 0 bhosts bqueues lsrun xlsbatch