Difference between revisions of "FSRoot"

From GlueXWiki
Jump to: navigation, search
(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.


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