Difference between revisions of "July 2013 Tutorial Schedule WorkArea"
From GlueXWiki
(→ANALYSIS Library (Paul)) |
(→ANALYSIS Library (Paul)) |
||
Line 50: | Line 50: | ||
** Give link for setting it up manually | ** Give link for setting it up manually | ||
** Track Preselection, PID FOM Calculation, new DTrackTimeBased for new PIDs (q+/q- subtlety), Event RF Bunch Selection, Thrown/Reconstructed track matching | ** Track Preselection, PID FOM Calculation, new DTrackTimeBased for new PIDs (q+/q- subtlety), Event RF Bunch Selection, Thrown/Reconstructed track matching | ||
− | * Exercise ( | + | * Exercise (45 mins): |
** Auto-kinematic fit 3pi_n, histogram kinematic fit confidence level, missing mass after kinematic fit (conlev > 1%) | ** Auto-kinematic fit 3pi_n, histogram kinematic fit confidence level, missing mass after kinematic fit (conlev > 1%) | ||
** Manually histogram: pid conlevs (dE/dx, timing, total) | ** Manually histogram: pid conlevs (dE/dx, timing, total) |
Revision as of 19:30, 7 July 2013
Contents
- 1 Session 1, Wed 8:30am - 9:30am: Overview (Matt)
- 2 Session 2, Wed 9:30am - 10:20am: Monte Carlo Generators and Detector Simulation (Mark)
- 3 Break, Wed 10:20am - 10:45am
- 4 Session 3, Wed 10:45am - 11:45am: JANA/DANA (Paul)
- 5 ANALYSIS Library (Paul)
- 6 Break, Wed 2:50pm - 3:10pm
- 7 Break, Thur 10:30am - 10:50am
- 8 Boosted Decision Trees (Mike & Justin)
- 9 Break, Thur 2:45pm - 3:10pm
- 10 Amplitude Analysis (Matt)
Session 1, Wed 8:30am - 9:30am: Overview (Matt)
- Talk (30 min): Physics Overview and Analysis Strategy
- Define the physics goals
- Discuss detector
- Discuss general analysis strategy
- Talk (15 min): A specific physics example: γp→π+π+π-n
- Exercise (15 min): Generate γp→π+π+π-n with intermediate resonances
Session 2, Wed 9:30am - 10:20am: Monte Carlo Generators and Detector Simulation (Mark)
- Talk (30 min): Other standard monte carlo tools (bggen, genr8, hdgeant, mcsmear)
- Exercise (20 min): Generate pythia data, run HDGeant and mcsmear
Break, Wed 10:20am - 10:45am
Session 3, Wed 10:45am - 11:45am: JANA/DANA (Paul)
- Talk (20 mins): Intro to JANA/DANA (factories & plugins)
- JANA is a framework: multithreaded, factory-based (describe concept of factories, tags (e.g. BCAL)), event-processors control the work flow
- The primary executable (hd_root) does virtually nothing by itself, the processor in your plugin tells it what to do
- JANA automatically loops over the events in your input files, executing the primary factory & plugin methods as needed: init(), brun(), evnt(), erun(), fini()
- hd_dump: Prints data to screen
- DMCThrown and DKinematicData
- Exercise (10 mins):
- Print the kinematics of the thrown particles to screen via hd_dump and inside of plugin (all three kinds)
- Talk (10 mins): Reconstruction and Data Classes
- Track/Shower reconstruction, Hypotheses, REST format, danarest plugin
- Basic Analysis Classes: DKinematicData, DChargedTrackHypothesis, DChargedTrack, DNeutralShower, DMCThrown, DBeamPhoton
- Exercise (20 mins):
- Use danarest plugin to do track reconstruction.
- Print basic analysis objects (DChargedTrackHypothesis, etc.) to screen via hd_dump and inside of plugin
Lunch, Wed 11:45am - 1:00pm
ANALYSIS Library (Paul)
Session 4, Wed 1:00pm - 1:50pm: Analysis Basics
- Talk (25 mins):
- User specifies reaction & decay chain, library generates all particle combinations
- Particle Combinations: Why, and why so many? (extra tracks, any PID of same q, decay products of different particles, beware double-counting), DParticleCombo contents
- DParticleCombo contents
- Basic example plugin: 3pi_n (go over line by line)
- Exercise (45 mins):
- Generate 3pi_n plugin, manually histogram: #combos, #events with >= 1 combo, missing mass (calc manually), p-vs-theta & vertex-z for pi+'s and pi-'s
- Generate decay plugin, manually histogram: #combos, #events with >= 1 combo, missing mass (calc manually), p-vs-theta & vertex-z for pi+'s and pi-'s
Session 5, Wed 1:50pm - 2:50pm: Kinematic Fitting & Analysis Details
- Talk (30 mins):
- What it is, how it works, constraints, how neutrals & magnetic field are handled
- Give link for setting it up manually
- Track Preselection, PID FOM Calculation, new DTrackTimeBased for new PIDs (q+/q- subtlety), Event RF Bunch Selection, Thrown/Reconstructed track matching
- Exercise (45 mins):
- Auto-kinematic fit 3pi_n, histogram kinematic fit confidence level, missing mass after kinematic fit (conlev > 1%)
- Manually histogram: pid conlevs (dE/dx, timing, total)
Break, Wed 2:50pm - 3:10pm
Session 7, Wed 4:00pm - 4:50pm: Analysis Actions - Pre-existing
- Talk (25 mins):
- Reaction-Independent (can be called manually): thrown, reconstructed track distros, gen-recon comparison (go over line-by-line on one)
- Reaction-Dependent Histograms (executed for each combo): PID, kinfit results, missing mass, invariant mass (go over line-by-line on one)
- Reaction-Dependent Cuts (executed for each combo): PID, kinfit conlev, missing mass, invariant mass (go over line-by-line on one)
- #events/combos pass cuts histograms
- Exercise (25 mins):
- Add reaction-independent analysis actions: fill histograms and look at kinematics in ROOT
- Add reaction-dependent analysis actions: hist/cut pid, hist/cut kinematic fit results, hist missing mass before/after kinfit conlev cut
Session 8, Wed 4:50pm - 5:30pm: Analysis Actions - Custom
- Talk (15 mins):
- Creating your own custom actions, remind to beware double-counting, warn about multi-threading issues
- Exercise (25 mins):
- Generate custom action to make dalitz plot of 3 pions.
Session 9, Thur 8:30am - 9:30am: REST & ROOT
- Talk (25 mins):
- REST Skims, TTree format, Making/Using a TSelector, beware double-counting when filling ROOT trees.
- Exercise (35 mins):
- Create a skim of 3pi_n events, analyze the skim data, save the results to a root ttree, analyze the data in ROOT (make some plots)
Session 10, Thur 9:30am - 10:30am: Analyzing Decays
- Talk (5 mins):
- Setting up DReaction
- Exercise (55 mins):
- Make a plugin to analyze: γp→η'p, η'→ηπ+π-, η→γγ
- Do kinematic fitting, output to ROOT, histogram PID, kinfitresults, mass peaks, etc.
- Make a plugin to analyze: γp→η'p, η'→ηπ+π-, η→γγ
Break, Thur 10:30am - 10:50am
Boosted Decision Trees (Mike & Justin)
Session 11, Thur 10:50am - 12:20pm: BDT Introduction
- Talk (45 mins): Multivariate Analysis Overview
- Talk (20 mins):
- Introduce TMVA, setup factory, select input variables, book methods, and training samples
- Exercise (25 mins):
- Train BDT using 3pi_n samples, make plots of input variables and classifier response with TMVAGui
Lunch, Thur 12:20pm - 1:30pm
Session 12, Thur 1:30pm - 2:15pm: BDT Details
- Talk (20 mins):
- Describe weight files, setup reader, evaluation on larger samples, and storing classifier response in TTree
- Exercise (25 mins):
- Apply BDT to large background sample, compute efficiency and purity
Session 13, Thur 2:15pm - 2:45pm: BDT Details
- Talk (10 mins):
- Describe setup for other Multivariate Classifiers
- Exercise (20 mins):
- Try many different classifiers for 3pi_n channel and evaluate with TMVAGui
Break, Thur 2:45pm - 3:10pm
Amplitude Analysis (Matt)
Session 14, Thur 3:10pm - 4:15pm: Amplitude Analysis Introduction
- Talk (45 mins): Amplitude Analysis Fitting Strategy
- Exercise (20 mins): performing a single fit and analyzing projections
Session 15, Thur 4:15pm - 5:05pm: Amplitude Analysis Details
- Talk (30 mins): Multiple fits, opportunities for extensions
- Exercise (20 mins): performing multiple fits in different bins of a mass spectrum