Old Physics Analysis Old Physics Analysis
From GlueXWiki
Revision as of 21:53, 29 July 2013 by Pmatt (Talk | contribs) (Created page with "This page is intended to provide access to the various physics analysis topics that are being pursued with the GlueX detector. At this stage these typically consist of Monte Car...")
This page is intended to provide access to the various physics analysis topics that are being pursued with the GlueX detector. At this stage these typically consist of Monte Carlo studies which help demonstrate the feasibility of a particular analysis.
Contents
Ongoing Analyses
- Grid issues, reconstruction, and PWA of pi+pi+pi-n (Jake Bennett)
- Particle ID and tracking issues using K+K-p and pi+pi-p (Kei Moriya)
- Kinematic fitting in b1pi (Will Levine)
- Simulations of cascade baryons (Nathan Sparks)
- Analysis of gamma p --> eta pi0 p (Irina Semenova and Andrei Semenov)
Completed Analyses
- Analysis of gamma p --> eta pi0 p (Blake Leverington)
Physics at the Information Frontier Center
Expected Statistics for GlueX
- The expected statistics for GlueX in 1 year of running will depend on several crucial factors, summarized below.
- Initial studies for the expected data rates for GlueX were done in
- The experimental data input for the cross sections at 9 GeV come from SLAC, namely
- H. H. Bingham et al. (SLAC) Phys. Rev. D8, 1277
- J. Ballam et al. (SLAC) Phys. Rev. D7, 3150
- and references therein
- SLAC photoproduction results
124.02.5 b | |
14.7 0.6 b | |
0.58 0.05 b | |
3.2 0.7 b |
- Basic formula to calculate expected events: where
- : expected raw events
- : cross section for wanted reaction
- : density of target
- : number of photons/s
- : integrated time
- The total number of events that we actually detect is then where
- : expected construction efficiency for this channel
- : the branching fraction of the particular decay to this channel
so that , and the number of events is the product of 6 factors.
- The target density, is the easiest to determine. With a 30cm target, the density is
- The photon rate is expected to start at , with a gradual increase to the maximum of
- For the integrated time , we assume, as in the previous studies listed above, that we run for 1 year, of which 6 months is datataking, with 1/3 efficiency (effectively 2 months out of the first year), which gives us s (this efficiency also takes into account the tagging ratio)
- The most difficult estimates are the factors of (the cross section), (the efficiency of reconstruction), and (the branching ratio), all of which are dependent on the specific channel of interest. Here we take the channel , since this is one of the prime channels of interest, and Monte Carlo studies have been done on this channel. According to the SLAC data above, the cross section for is , of which an unknown fraction will be our signal of interest, the . Current (2011) efficiency studies by Jake Bennett at IU show that the reconstruction efficiency in the mass region of interest is around , with hopefully an increase as the kinematic fitter is further developed. Taking , we come to our final estimate of
events.
- This number can be contrasted with the estimates given in the above references, for example, assuming and and the same flux and integrated time as above gives events in the 1st year.
- Also to be contrasted is the number of events in other analyses. These include:
- Nozar et al. (CLAS) photoproduction of : 83k events for PWA
- Dzierba et al. (E852-IU) pionproduction of : 2.6M events
- Alekseev (COMPASS) pionproduction of : 420k events
- In conclusion, even with the relatively pessimistic numbers of 2 months effective running, , , GlueX should accumulate on the order of several hundred million events in the channels of interest, which is at least 2-3 orders of magnitude higher than other experiments.
Documentation/Studies
- Potential GlueX Physics During the First Two Years of 12 GeV Running
- Talk:The First Two Years of GlueX Physics
- Simulations of Exotic Mesons in GlueX
- Calorimeter physics simulations and reconstruction for the Calorimetry Review (GlueX Portal)
- PrimEx at GlueX
- Photoproduction & Lattice QCD
- A Summary of Hybrid Masses and Decays
- Radiation length scans
- b1pi with GlueX Monte Carlo
Software Analysis Tools
- JANA-based GlueX Analysis Framework:
- GlueX Analysis Library: JANA-based, details on GlueX reconstruction classes, event selection, etc.
- Partial Wave Analysis Tutorials
- pwa-ruby wiki : An installation and usage howto for the ruby-pwa code. The site also includes a fully worked out simple example.
- The qft++ code
- The qft++ code was developed by Mike Williams at Carnegie Mellon University to facilitate analytic calculation of tree-level diagrams for partial wave analysis.
- Details can be found on the qft++ Website.
- A paper has also been submitted to Comp. Phys. Comm. (2008). arXiv:0805.2956
- Software Packages
- The Ruby-Minuit Package Ruby-Minuit is a Ruby binding to CERNLIB's MINUIT package (which is written in FORTRAN). It provides access to MINUIT's powerfull minimization routines from the infinitely flexible Ruby language.
- The qft++ package was written (in C++) to numerically calculate quantum field theory expressions. The API was designed such that the code strongly resembles the hand-written expressions.
- Maximum Likelihood Fitting.
- A short note on maximum Likelihood Fitting.
- A generalized background subtraction method was developed to facilitate the baryon PWA done at CMU. This has been published as JINST 4 P10003, (2009), arXiv:0809.2948.
Useful Articles
Meetings
Physics Working Group Meetings