UmassNtupleAnalysis

UMAnalysis

The newest version of the ntuple maker has had its name changed to UMAnalysis. The framework has a different method of filling variables from a tree, the rest of the framework is the same.

The new makeSelector class is called anyTree.C, and now only loads branches which you specify in a configuration file, such as oxfordConfig.cxx. For each branch, it checks for the branches existance and type, and based on this places the branch data in a map, with a string specifying what variable the data corresponds with. The advantage of this is that you only need to change the configuration file to change the ntuple you use, and instead of getting a crash when your ntuple changes, you will just get a list of which variable failed to load.

It is functional currently for oxford ntuples, which exist at the umass2 machine, and is being adapted for D3PDMaker ntuples. More info to come!

UMAnalysis

This twiki is under continuous development. Please email me if you have any questions/problems (ameade@physicsNOSPAMPLEASE.umass.edu).

UMAnalysis is an Ntuple analysis framework framework. It has two basic features:

  • Generalized Ntuple Reader: Takes a list of branches as input, loading those branches only. It performs the following checks: the branch exists, the type is expect, and that the variable is loaded into the analysis correctly. Errors are stored and printed. Should work with any ntuple, and should not crash if the ntuple used changes.
  • Internal Event Data Model (EDM): Its own simple particle classes, particle selection, hlstograming, etc. The goal is to have tools/structures to build your analysis that can be independent of changes in the data/ntuple format used.

The next two sections will be a howto, describing how to start using the package, and then a more thorough description of the analysis package and its features.

How To Download and Use the Framework

The package can be downloaded from anywhere you have access to SVN, probably hal or lxplus. This tutorial assumes you are at lxatlasumass2, where oxford analysis ntuples are already present to be used. See below for other ntuple types. The SVN Group Area ($SVNGRP) is set when you setup athena, so you must do this first. If you are not at lxplus, you will also need to kinit. For more information on setting up cmt and you work area, see the WorkBook.

source ~/cmthome/setup.sh -tag=15.5.1
svn co $SVNGRP/Institutes/UMass/UMAnalysis/trunk UMAnalysis
cd UMAnalysis/run
root -n -b -q oxLoadAnalysis.C

When you check out the package, the last argument is the name of the directory you create. The oxLoadAnalysis.C macro compiles and executes all classes necessary for the analysis, and by default runs 10 events on the included sample file, writing output histograms to testoutput.root. To checkout a specific version of the framework, and run on your own file, do:

svn co $SVNGRP/Institutes/UMass/UMAnalysis/tag/UMAnalysis-XX-YY-ZZ UMAnalysis
cd UMAnalysis/run
root -n -b -q 'oxLoadAnalysis.C("/path/to/my/file/EWPA.D3PD.root","outputName.root", eventsToProcess)'

The option -n tell root not to load the log on script, this can be necessary, because some log on scripts interfere with the framework.

Structure of the Analysis Package

The following are the files you will want to look at and alter, in order of importance to the end user:

  • myAnalysis.cxx: This is where the final analysis is done on internal particle classes (muons,jets,etc), which are accessed through m_data, for example m_data->preMuons gives you a pointer to the vector of muons before selection.

  • oxfordConfig.cxx: Here is where the configuration of the analysis is done. So you can change for example the minimum p_T of muons you will select, or what jet algorithms to consider. Note that all the branch name information is set here, so there is a different config file for each ntuple type. Each of them set a number of static variables in myConfig.h, that are then picked up by the other classes. Example:
     myConfig::MuonBranchName 

  • dataManager.cxx: This class stores the particle vectors in the struct particleData. It fills the internal particle classes from the branch information, and selects a subset of these particles that you actually want to use in your analysis (for example m_data->selMuons), based on a number of cuts on momentum, angular distribution and quality criteria (for example isolation).

  • particle.cxx: These classes define particles (muon, egamma, jet, etc). Each contains a momentum TLorenzVector and some amount of specific info such as vertex info and quality info. Each cxx file defines the ranges for all the cuts on that particle (though which cuts are actually used are defined in dataManager::selectParticles() ). Each particle has a method selectParticle() which returns a particleSelection object.

  • particleSelection.h: These classes define an integer in which each bit corresponds to a cut. There is an enumeration for each cut. So when selectParticle() is called, the particleSelection object contains which cuts the particle passes. This object also contains a number of methods to manipulate and evaluate these cuts, most important being passesCuts(bits), where bits is an integer specifying which cuts you want the particle to pass. If the particle passes all the cuts, it returns true. Examples are in dataManager::selectParticles(). Note that in the bitword, 0 specifies a bit that IS passed.

  • anyTree.C/.h: This is a custom TSelector class: it allocates the data vectors to which the branches are writen, and then defines an event loop in which all the other classes are called. It is very different from your typical makeSelector() class, in that it only reads the branches specified in the config file, and all other branches are turned off, which improves performance. In addition, it checks for the existence of the branches before reading them/allocating memory, and as such avoids the crashes that often plague this type of ntuple analysis.

Advanced Features and Utilities

The above functionality should remain constant, however, new features and utilities are expected to be added as time goes on

  • HistoMgr: This class allows you to avoid having to declare, initialize, and fill histograms in different files. Instead, you just call its "fill" function, and it does the declaration/initialization for you. Because of its implementation (as a singleton) it can be used across any number of files, and is already setup in myAnalysis, dataManager, and WZanalyMod, using the following:
       HistoMgr * m_histoMgr = HistoMgr::Instance();
       //...
       m_histoMgr->fillHisto(histoName, histoTitle, nBins, xLow, xHigh, fillValue, weight);
       

  • analysisModules: These are classes to contain blocks of self-contained analysis code, that can be then included and called from myAnalysis. They only take particleData (m_data) as input. This allows us to save our own code to the EWAnalysis package without changing myAnalysis, and also allows us to easily turn-on/turn-off blocks of analysis code. Right now only the following exist: myModule.cxx, a dummy analyisis, and WZacceptMod.cxx, a module I wrote to look at systematics for W/Z studies.

  • CutflowMgr: This class keeps track of particle and event cuts you make in your analysis, and then prints a summary of the efficiencies of each of these cuts. It is configured automatically for particle selection cuts. For event level cuts, it is called using the following:
       m_cutMgr->applyCut("objectBeingCut","nameOfCut",passesCut);
       //... more cuts ...
       m_cutMgr->cutsDoneFor("objectBeingCut");
       
    This last command lets the tool know the event/particle/other-object is done being cut. After the analysis is over, myAnalysis calls m_cutMgr->printCuts(). A sample of the output is attached here: cutflowOutput.txt. The sequential efficiency is out of events passing the cuts up till that point, the raw efficiency is out of all events, and the exclusive efficiency is out of events where all other cuts are passed.

Supported Ntuples

Right now only Oxford Analysis Ntuples are supported. However Tulin and I are in the process of making a configuration file for D3PDMaker Ntuples. In addition, I have all the information to support EWPA ntuples, let me know if you need this functionality.

Adding Variables

If you want a variable added, say muon chi2 information, you just need to do the following: (1) add it to yourConfig.cxx initMuonMap (where the postfixes are added), (2) add the variable to muon.h, and then go (3) dataManager::makeMuon() and add a line to fill the variable.

(1)  muonMap["fitChi2"] = "trackfitchi2"; // in your config file
(3)   aMuon->fitChi2 =  Get<double>("muon_fitChi2",i);    //in dataManager::makeMuon()
Make sure that the type of the variable you add is the same type as the branch. Get will attempt to cast the branch found to that type, which in some cases could cause problems.

Other Ntuple Types

It should be very easy to alter this for your favorate ntuple type. Just change the branch name prefixes (ex MuidMuon) and postfixes (ex. _px) listed in oxfordConfig.cxx to your own! That should be all. If you do, let me know, and put the config file in the package to share! Contact me if there are any problems.

Available Ntuples

As mentioned on lxatlasumass2 a selection of oxford ntuples exist at /data/atlas/data/oxNtuple. D3PD ntuples will be stored somewhere soon... To make your own ntuples, see one of the following pages:

-- AndrewMeade - 18-Jan-2010

Topic attachments
I Attachment History Action Size Date Who Comment
Texttxt cutflowOutput.txt r1 manage 0.6 K 2010-02-03 - 17:21 AndrewMeade sample output from CutflowMgr
Edit | Attach | Watch | Print version | History: r16 | r10 < r9 < r8 < r7 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r8 - 2010-02-26 - AndrewMeade
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback