Difference between revisions of "Reconstruction, June 9, 2009"
From GlueXWiki
Line 22: | Line 22: | ||
* [[Reconstruction References]] | * [[Reconstruction References]] | ||
− | # | + | # [http://argus.phys.uregina.ca/gluex/DocDB/0008/000856/001/pythia_data_v1.pdf Gluex-doc-856] |
− | # | + | # [http://argus.phys.uregina.ca/gluex/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 |
# bggen and Pythia | # bggen and Pythia | ||
## HOW-TO: | ## HOW-TO: | ||
Line 39: | Line 37: | ||
# 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]] |
Revision as of 00:43, 9 June 2009
Time: 11:00 EST/10:00 CST
- Ventrilo for sound. (Zisis will connect) See GlueX Communications for detailed instructions
- ESNET Number 8542553
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 possible pi0s and etas and a2.
- HOW-TO:
- Discussion on Analysis
Above:reconstructed energy ratio as a function generated energy