Difference between revisions of "FSRoot"
From GlueXWiki
(→Here is a working example:) |
|||
Line 24: | Line 24: | ||
> flatten -in /u/home/remitche/FSRootExample/analysis.root -out flat.root | > flatten -in /u/home/remitche/FSRootExample/analysis.root -out flat.root | ||
+ | |||
+ | * where analysis.root was created by running /u/home/remitche/FSRootExample/MC.run | ||
+ | * which generates MC based on the config file /u/home/remitche/FSRootExample/MC.config | ||
+ | * which simulates the process gamma p --> eta' p --> eta pi+ pi- p --> (pi+ pi- pi0) pi+ pi- p using genr8 with the input file /u/home/remitche/FSRootExample/MC.genr8 | ||
2. Make a plot of the eta' mass: | 2. Make a plot of the eta' mass: | ||
root> FSHistogram::getTH1F("flat.root","ntFSGlueX_100_221","MASS(2,3,4,5,6)","(80,0.8,1.2)","")->Draw() | root> FSHistogram::getTH1F("flat.root","ntFSGlueX_100_221","MASS(2,3,4,5,6)","(80,0.8,1.2)","")->Draw() | ||
+ | |||
+ | * where particles 2-6 are the five pions from the eta' decay |
Revision as of 11:54, 10 December 2021
FSRoot is an alternative to DSelectors that uses flat trees.
There are two steps:
1. Use FlattenForFSRoot (hd_utilities/FlattenForFSRoot/flatten) to generate flat trees from the standard analysis trees.
> flatten -in [input file with analysis tree] -out [output file with flat tree]
- Run flatten with no flags to see further run options.
- This should work for data, reconstructed MC, and thrown MC.
- See the FlattenForFSRoot/Documentation directory for more information.
2. Use FSRoot (available at https://github.com/remitche66/FSRoot) to make histograms from the flat trees.
- See the FSRoot/Documentation directory for documentation.
Here is a working example:
1. Flatten an analysis tree:
> flatten -in /u/home/remitche/FSRootExample/analysis.root -out flat.root
- where analysis.root was created by running /u/home/remitche/FSRootExample/MC.run
- which generates MC based on the config file /u/home/remitche/FSRootExample/MC.config
- which simulates the process gamma p --> eta' p --> eta pi+ pi- p --> (pi+ pi- pi0) pi+ pi- p using genr8 with the input file /u/home/remitche/FSRootExample/MC.genr8
2. Make a plot of the eta' mass:
root> FSHistogram::getTH1F("flat.root","ntFSGlueX_100_221","MASS(2,3,4,5,6)","(80,0.8,1.2)","")->Draw()
- where particles 2-6 are the five pions from the eta' decay