Difference between revisions of "Reconstruction, June 9, 2009"
From GlueXWiki
m (Text replacement - "http://argus.phys.uregina.ca/gluex" to "https://halldweb.jlab.org") |
|||
(9 intermediate revisions by one other user not shown) | |||
Line 22: | Line 22: | ||
* [[Reconstruction References]] | * [[Reconstruction References]] | ||
− | # | + | # [https://halldweb.jlab.org/DocDB/0008/000856/001/pythia_data_v1.pdf Gluex-doc-856] |
− | # | + | # [https://halldweb.jlab.org/DocDB/0009/000985/002/GX_calorimetry.pdf GlueX Calorimeter Review Feb 2008] |
− | + | ||
− | + | ||
= Tentative Agenda = | = Tentative Agenda = | ||
# Review Reconstruction Plans in list above | # Review Reconstruction Plans in list above | ||
− | # | + | # a2_1320 |
+ | ##[[Image:Lg triangle yellow questionmark.gif|20px]]η's and π0's are not separate objects (as much as I understand) that I can apply cuts to when looking at the 4 photon invariant mass in the current code. The η doesn't know that I also made a π0. | ||
+ | ##proton reconstruction is now included via HDParSim | ||
# bggen and Pythia | # bggen and Pythia | ||
## HOW-TO: | ## HOW-TO: | ||
Line 36: | Line 36: | ||
##* filter_HighE: Filter out events with total neutral energy less than 6GeV. | ##* filter_HighE: Filter out events with total neutral energy less than 6GeV. | ||
##* hdgeant: run the background thrown data set through hdGeant. | ##* hdgeant: run the background thrown data set through hdGeant. | ||
− | ##* fCalTree: reconstruction possible pi0s and etas | + | ##* fCalTree: reconstruction of possible pi0s and etas.[[Image:Lg triangle yellow questionmark.gif|20px]] η's and π0's are not separate objects |
− | + | ||
# Discussion on Analysis | # Discussion on Analysis | ||
+ | |||
+ | == a2(1320) analysis == | ||
+ | |||
+ | |||
+ | [[Image:PhotEThrownE 2D.png|400px| reconstructed photon energy ratio as a function of thrown neutral energy|]][[Image:PhotEThrownE fit.png|400px| reconstructed photon energy ratio]] | ||
+ | |||
+ | Above:reconstructed energy ratio as a function generated energy | ||
+ | |||
+ | [[Image:Photonmultiplicity.png|400px]] | ||
+ | |||
+ | Above:Photon multiplicity | ||
+ | |||
+ | [[Image:4gammamass.png|400px]] | ||
+ | |||
+ | Above: invariant mass distribution selecting events with only 4 reconstructed photons | ||
+ | |||
+ | == pythia background == | ||
+ | |||
+ | [[Image:Bggen photEThrownE 2D.png|400px]] | ||
+ | |||
+ | Above: photon energy fraction as a function of Thrown energy. '''Full of Fail.''' | ||
+ | |||
+ | [[Image:Bggen photmultiplicity.png|400px]] | ||
+ | |||
+ | Above: Pythia photon multiplicity after neutral energy filtering and reconstruction. | ||
+ | |||
+ | |||
+ | [[Image:Bggen 4gammaEtot nPhot4.png||400px]] | ||
+ | |||
+ | Above: invariant mass distribution selecting events with only 4 reconstructed photons from Pythia background. |
Latest revision as of 12:38, 14 March 2017
Time: 11:00 EST/10:00 CST
- Ventrilo for sound. (Zisis will connect) See GlueX Communications for detailed instructions
- ESNET Number 8542553
Contents
Reconstruction To Do List
- Generation
- Use genr8 for a2 production
- Skim off (hd_filter) the wanted events by cutting on the total neutral energy (> 6 GeV) but record the total number of events before skimming. The resulting MC data set will be stored for general use.
- Use Pythia for background generation after removing a2 → η π0 from the MC data stream as well as events with total neutral energy <= 6GeV.
- Use a2 charge exchange production xsection vs pythia total xsection to estimate number of signal and noise events 124μb @ 9GeV ?
- Analysis
- After pruning the Pythia events, check the generated vs. reconstructed energy in a plot.
- Look at the a2 mass peak - adjust cuts on the eta and pi0 masses then look at the effect on the different angular distributions (i.e. Gottfried-Jackson angles, etc.)
- look at difference between reconstructed and generated z position for less than 65 cm
- look at single photon reconstruction efficiencies
- adjust the error matrix for the BCAL (look in PID/DPhoton_factory.cc)
- Gather all relevant DocDB documents and cross section publications at one spot on the Wiki.
Documents to Review
Tentative Agenda
- Review Reconstruction Plans in list above
- a2_1320
- bggen and Pythia
- HOW-TO:
- bggen: by default uses the photon beam spectra over 0.15 to 12 GeV. Cross-section??. This has to be modified in fort.15 to 9GeV for us to compare cross sections. This also now dramatically increases our background data set by a factor of 10 over the default beam energy. Why is there a limit on 2.2 million events?
- filter_HighE: Filter out events with total neutral energy less than 6GeV.
- hdgeant: run the background thrown data set through hdGeant.
- fCalTree: reconstruction of possible pi0s and etas. η's and π0's are not separate objects
- HOW-TO:
- Discussion on Analysis
a2(1320) analysis
Above:reconstructed energy ratio as a function generated energy
Above:Photon multiplicity
Above: invariant mass distribution selecting events with only 4 reconstructed photons
pythia background
Above: photon energy fraction as a function of Thrown energy. Full of Fail.
Above: Pythia photon multiplicity after neutral energy filtering and reconstruction.
Above: invariant mass distribution selecting events with only 4 reconstructed photons from Pythia background.