|
|
| Line 1: |
Line 1: |
| − | 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.
| |
| | | | |
| − | == 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)
| |
| − |
| |
| − | * [http://zeus.phys.uconn.edu/wiki/index.php/Exotic_b1%CF%80_Channel_Simulation_and_Analysis Analysis of omega pi+ pi-] ([http://zeus.phys.uconn.edu/wiki/index.php/User:Senderovich Igor Senderovich])
| |
| − |
| |
| − | * 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)
| |
| − | ** [http://argus.phys.uregina.ca/cgi-bin/private/DocDB/ShowDocument?docid=1604 University of Regina PhD Thesis (see chapter 7)]
| |
| − |
| |
| − | == Physics at the Information Frontier Center ==
| |
| − | * [http://zeus.phys.uconn.edu/catwiki/ Collaborative Analysis Toolkit Wiki]
| |
| − |
| |
| − | == 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
| |
| − | ** [http://argus.phys.uregina.ca/cgi-bin/public/DocDB/ShowDocument?docid=1226 GlueX PAC30 Documentation]
| |
| − | ** [http://argus.phys.uregina.ca/cgi-bin/private/DocDB/ShowDocument?docid=621 Potential GlueX Physics During the First Two Years of 12 GeV Running]
| |
| − | ** [http://argus.phys.uregina.ca/cgi-bin/private/DocDB/ShowDocument?docid=624 Talk:The First Two Years of GlueX Physics]
| |
| − | * The experimental data input for the cross sections at 9 GeV come from SLAC, namely
| |
| − | ** [http://link.aps.org/doi/10.1103/PhysRevD.8.1277 H. H. Bingham et al. (SLAC) Phys. Rev. D8, 1277]
| |
| − | ** [http://link.aps.org/doi/10.1103/PhysRevD.7.3150 J. Ballam et al. (SLAC) Phys. Rev. D7, 3150]
| |
| − | ** and references therein
| |
| − | * SLAC photoproduction results
| |
| − | {| border="1"
| |
| − | |<math>\sigma_{tot}</math> || 124.0<math>\pm</math>2.5 <math>\mu</math>b
| |
| − | |-
| |
| − | |<math>p \pi^{+} \pi^{-}</math> || 14.7 <math>\pm</math>0.6 <math>\mu</math>b
| |
| − | |-
| |
| − | |<math>p K^{+} K^{-}</math> || 0.58 <math>\pm</math>0.05 <math>\mu</math>b
| |
| − | |-
| |
| − | |<math>n 2 \pi^{+} \pi^{-}</math> || 3.2 <math>\pm</math>0.7 <math>\mu</math>b
| |
| − | |}
| |
| − | * Basic formula to calculate expected events: <math> N_{raw} = \sigma \times N_{t} \times N_{\gamma} \times T</math> where
| |
| − | ** <math>N_{raw}</math>: expected raw events
| |
| − | ** <math>\sigma</math>: cross section for wanted reaction
| |
| − | ** <math>N_{t}</math>: density of target
| |
| − | ** <math>N_{\gamma}</math>: number of photons/s
| |
| − | ** <math>T</math>: integrated time
| |
| − | * The total number of events that we actually detect is then <math> n = N_{raw} \times \epsilon \times BR</math> where
| |
| − | ** <math>\epsilon</math>: expected construction efficiency for this channel
| |
| − | ** <math>BR</math>: the branching fraction of the particular decay to this channel
| |
| − | so that
| |
| − | <math> n = \sigma \times N_{t} \times N_{\gamma} \times T \times \epsilon \times BR </math>, and the number of events is the product of '''6''' factors.
| |
| − | * The target density, <math>N_{t}</math> is the easiest to determine. With a 30cm <math>LH_{2}</math> target, the density is <math>1.26 b^{-1}</math>
| |
| − | * The photon rate is expected to start at <math>10^{7} \gamma/s</math>, with a gradual increase to the maximum of <math>10^{8} \gamma/s</math>
| |
| − | * For the integrated time <math>T</math>, 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 <math>5 \times 10^{6}</math> s (this efficiency also takes into account the tagging ratio)
| |
| − | * The most difficult estimates are the factors of <math>\sigma</math> (the cross section), <math>\epsilon</math> (the efficiency of reconstruction), and <math>BR</math> (the branching ratio), all of which are dependent on the specific channel of interest. Here we take the channel <math>\pi \pi \pi (n)</math>, 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 <math>n 2 \pi^{+} \pi^{-}</math> is <math>\sigma=3.2 \mu b</math>, of which an unknown fraction will be our signal of interest, the <math>\pi_{1}(1600)</math>. Current (2011) efficiency studies by Jake Bennett at IU show that the reconstruction efficiency in the mass region of interest is around <math>20%</math>, with hopefully an increase as the kinematic fitter is further developed. Taking <math>BR=1</math>, we come to our final estimate of
| |
| − | <math> n = 3.2 [\mu b] \times 1.26 [b^{-1}] \times 10^{7} [s^{-1}] \times 5 \times 10^{6} [s] \times 0.20 \times 1 = 400 M</math> events.
| |
| − | * This number can be contrasted with the estimates given in the above references, for example, <math>p \pi^{+} \pi^{-}</math> assuming <math>\sigma = 20 \mu b</math> and <math>\epsilon=0.75</math> and the same flux and integrated time as above gives <math>940M</math> events in the 1st year.
| |
| − | * Also to be contrasted is the number of events in other analyses. These include:
| |
| − | ** [http://link.aps.org/doi/10.1103/PhysRevLett.102.102002 Nozar et al. (CLAS)] photoproduction of <math>\pi^{+} \pi^{+} \pi^{-} (n)</math>: 83k events for PWA
| |
| − | ** [http://link.aps.org/doi/10.1103/PhysRevD.73.072001 Dzierba et al. (E852-IU)] pionproduction of <math>\pi^{-} \pi^{-} \pi^{+} p</math>: 2.6M events
| |
| − | ** [http://link.aps.org/doi/10.1103/PhysRevLett.104.241803 Alekseev (COMPASS)] pionproduction of <math>\pi^{-} \pi^{-} \pi^{+} p</math>: 420k events
| |
| − | * In conclusion, even with the relatively pessimistic numbers of <math>T =</math> 2 months effective running, <math>\epsilon=0.20</math>, <math>N_{\gamma}=10^{7} \gamma/s</math>, 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 ==
| |
| − | * [http://argus.phys.uregina.ca/cgi-bin/private/DocDB/ShowDocument?docid=621 Potential GlueX Physics During the First Two Years of 12 GeV Running]
| |
| − | * [http://argus.phys.uregina.ca/cgi-bin/private/DocDB/ShowDocument?docid=624 Talk:The First Two Years of GlueX Physics]
| |
| − | * [http://argus.phys.uregina.ca/cgi-bin/private/DocDB/ShowDocument?docid=1232 Simulations of Exotic Mesons in GlueX]
| |
| − | * [http://argus.phys.uregina.ca/gluex/DocDB/0009/000989/001/simulation.pdf Calorimeter physics simulations and reconstruction for the Calorimetry Review (GlueX Portal)]
| |
| − | * [[PrimEx at GlueX]]
| |
| − | * [[Photoproduction & Lattice QCD]]
| |
| − | * [http://argus.phys.uregina.ca/cgi-bin/private/DocDB/ShowDocument?docid=380 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_Software | GlueX Analysis Library]]: JANA-based, details on GlueX reconstruction classes, event selection, etc.
| |
| − | # Partial Wave Analysis Tutorials
| |
| − | ## [http://www-meg.phys.cmu.edu/williams/wiki-ruby-pwa/ pwa-ruby wiki] : An installation and usage <i>howto</i> 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 [http://www-meg.phys.cmu.edu/williams/qft++/ qft++ Website].
| |
| − | ##A paper has also been submitted to Comp. Phys. Comm. (2008). [http://xxx.lanl.gov/abs/0805.2956 arXiv:0805.2956]
| |
| − | #Software Packages
| |
| − | ## [http://www-meg.phys.cmu.edu/williams/ruby-minuit/ 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.
| |
| − | ## [http://www-meg.phys.cmu.edu/williams/qft++/index.php 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 [http://argus.phys.uregina.ca/cgi-bin/private/DocDB/ShowDocument?docid=666 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), [http://xxx.lanl.gov/abs/0809.2548 arXiv:0809.2948].
| |
| − |
| |
| − | == Useful Articles ==
| |
| − |
| |
| − | * [http://www.curtismeyer.com/articles Article Archive from Curtis]
| |
| − | * [[Reconstruction References]]
| |
| − |
| |
| − | == Meetings ==
| |
| − |
| |
| − | [[Physics Working Group Meetings]]
| |
| − |
| |
| − | [[BCAL Reconstruction Work]]
| |
| − |
| |
| − | == Useful Links ==
| |
| − |
| |
| − | * [http://www.int.washington.edu/PROGRAMS/09-3.html Institute for Nuclear Theory (INT) workshop in November 2009]
| |
| − | * [http://conferences.jlab.org/php2008/index.html Workshop on photon-hadron physics with the GlueX detector]
| |
| − | * [https://wiki.jlab.org/pwawiki JLab PWA wiki]
| |
