Difference between revisions of "GlueX Analysis Workshop 2013"

From GlueXWiki
Jump to: navigation, search
(Agenda: Thursday July 11, 2013)
Line 85: Line 85:
 
*  '''5:05  Adjourn'''
 
*  '''5:05  Adjourn'''
  
== Schedule Details ==
+
== Schedule Work Area ==
 
+
* [[July_2013_Tutorial_Schedule_WorkArea | Schedule Work Area]]
=== 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: &gamma;p&rarr;&pi;<sup>+</sup>&pi;<sup>+</sup>&pi;<sup>-</sup>n
+
* Exercise (15 min):  Generate &gamma;p&rarr;&pi;<sup>+</sup>&pi;<sup>+</sup>&pi;<sup>-</sup>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 (15 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'''
+
* Talk (15 mins):
+
** What it is, how it works, constraints, how neutrals & magnetic field are handled
+
** Give link for setting it up manually
+
** Contents of DParticleCombo
+
* Exercise (30 mins):
+
** Auto-kinematic fit 3pi_n, histogram kinematic fit confidence level, missing mass after kinematic fit (conlev > 1%)
+
** Auto-kinematic fit decay channel, histogram kinematic fit confidence level, missing mass after kinematic fit (conlev > 1%)
+
 
+
=== Break, Wed 2:50pm - 3:10pm ===
+
 
+
'''Session 6, Wed 3:10pm - 4:00pm: Analysis Details'''
+
* Talk (25 mins):
+
** Thrown/Reconstructed track matching, PID FOM Calculation, new DTrackTimeBased for new PIDs (q+/q- subtlety), Event RF Bunch Selection, Track Preselection
+
** Different object versions (DParticleCombo, hypotheses)
+
** REST skims
+
* Exercise (25 mins):
+
** Make skims
+
** Manually histogram: pid conlevs (dE/dx, timing, total)
+
** Stretch goal: histogram pid confidence level vs. projected track start time for different PIDs for BCAL/TOF
+
 
+
'''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: &gamma;p&rarr;&eta;'p,  &eta;'&rarr;&eta;&pi;<sup>+</sup>&pi;<sup>-</sup>,  &eta;&rarr;&gamma;&gamma;
+
*** Do kinematic fitting, output to ROOT, histogram PID, kinfitresults, mass peaks, etc.
+
 
+
=== 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
+

Revision as of 18:04, 29 June 2013

Purpose

The purpose of this tutorial/workshop session will be to give an introduction to using the Hall-D analysis tools. It is open to anyone either inside or outside the GlueX collaboration who has an interest in doing physics inside Hall-D. The tutorial will cover the basics of:

  • Physics Analysis Framework
  • Using Boosted Decision Trees
  • Amplitude Analysis using AMPtools

The tutorial will consist of a few talks with room for questions/discussion.

Location and Time

The tutorial will take place:

DATE: July. 10-11, 2013

TIME: 9:00am to 4:30pm EDT (i.e. JLab time).

LOCATION: CEBAF Center L102 - L104

Registration

Registration Page Participants

Remote Participation

To connect from the outside by phone please do the following

1.) ESNET: 8542553

2.) Polycom RealPrescence App for iOS and Android. Call 8542553@gk1.es.net

3.) Telephone:

Dial +1-866-740-1260 : US and Canada

Dial +1-303-248-0285 : International

Then use participant code: 3421244# (the # is needed when using the phone) or www.readytalk.com then type access code 3421244 into "join a meeting" (you need java plugin)

4.) As of January 1, 2013 EVO became SeeVogh and Jefferson Lab has a contract to use it.

   To connect via SeeVogh, vist their website at research.seevogh.com. Your old EVO credentials are supposed to work. There will be a room booked for the meeting.
   To phone into an EVO meeting, from the U.S. call (626) 395-2112 and then enter the EVO meeting code. Instructions for the Phone Bridge to EVO.
   Skype Bridge to EVO 

5.) Please visit the Communications page for some options.

Agenda: Wednesday July 10, 2013

  • 8:30 Introduction and Monte Carlo (75) --- Chair: Paul Mattione
    • 8:30 Session 1a (15) --- Welcome and Getting Started --- Paul Mattione
    • 8:45 Session 1b (60) --- Overview and Signal Generation --- Matt Shepherd
    • 9:45 Session 1c (50) --- Monte Carlo Generators and Detector Simulation --- Mark Ito
  • 10:35 Coffee (25)
  • 11:00 JANA/DANA Overview (60) --- Chair: Justin Stevens
    • 11:00 Session 2a (30) --- Introduction to JANA --- Paul Mattione
    • 11:30 Session 2b (35) --- Reconstruction and Particle Data --- Paul Mattione
  • 12:05 Lunch (75)
  • 1:20 Introduction to Analyses in DANA (115) --- Chair: Matt Shepherd
    • 1:20 Session 3a (70) --- Analysis Basics --- Paul Mattione
    • 2:30 Session 3b (45) --- Kinematic Fitting --- Paul Mattione
  • 3:15 Coffee (25)
  • 3:40 DANA Analysis Details and Actions (135) --- Chair: Justin Stevens
    • 3:40 Session 4a (45) --- Analysis Details --- Paul Mattione
    • 4:25 Session 4b (50) --- Analysis Actions --- Paul Mattione
    • 5:15 Session 4c (40) --- Custom Actions --- Paul Mattione
  • 5:55 Adjourn

Agenda: Thursday July 11, 2013

  • 8:30 ROOT and Full Exercise (120) --- Chair: Matt Shepherd
    • 8:30 Session 5a (60) --- Analysis in ROOT --- Paul Mattione
    • 9:30 Session 5b (60) --- Exercise: Analyze γp→η'p, η'→ηπ+π-, η→γγ --- Paul Mattione
  • 10:30 Coffee (20)
  • 10:50 Boosted Decision Tree Introduction (90) --- Chair: Paul Mattione
    • 10:50 Session 6a (45) --- Multivariate Analysis Overview --- Mike Williams
    • 11:35 Session 6b (45) --- Introduction to TMVA --- Justin Stevens
  • 12:20 Lunch (70)
  • 1:30 Boosted Decision Tree Details (75) --- Chair: Matt Shepherd
    • 1:30 Session 7a (45) --- BDT Details --- Justin Stevens
    • 2:15 Session 7b (30) --- Other Multivariate Classifiers --- Justin Stevens
  • 2:45 Coffee (25)
  • 3:10 Amplitude Analysis (115) --- Chair: Justin Stevens
    • 3:10 Session 8a (65) --- Amplitude Analysis Introduction --- Matt Shepherd
    • 4:15 Session 8b (50) --- Amplitude Analysis Details --- Matt Shepherd
  • 5:05 Adjourn

Schedule Work Area