Tags:
tag this topic
create new tag
view all tags
<!-- keep this as a security measure: #uncomment if the subject should only be modifiable by the listed groups # * Set ALLOWTOPICCHANGE = Main.TWikiAdminGroup,Main.CMSAdminGroup # * Set ALLOWTOPICRENAME = Main.TWikiAdminGroup,Main.CMSAdminGroup #uncomment this if you want the page only be viewable by the listed groups # * Set ALLOWTOPICVIEW = Main.TWikiAdminGroup,Main.CMSAdminGroup,Main.CMSAdminReaderGroup --> ---+ Analysis of the 'Ntuples' This page describes the setup and some tools for analyzing the 'ntuples'. ---++ Quick and dirty plots directly from the ntuple <pre>setenv LD_LIBRARY_PATH /swshare/glite/d-cache/dcap/lib/:${LD_LIBRARY_PATH} root gSystem->Load("..../lib/slc4_ia32_gcc345/libAna00.so"); TFile f("Bu2JpsiKp.root") T1->Draw("fCandidates.fMass", "abs(fCandidates.fType)==521443") T1->Draw("fCandidates.fMass", "abs(fCandidates.fType)==521&&fCandidates.fVertex.fChi2<2") </pre> Please note that you should find the proper relative path (replace .... with something appropriate above) ---++ Setup Login to t3ui02 <pre>kinit user@CERN.CH cmsrel CMSSW_3_5_6 cd CMSSW_3_5_6/src cmsenv setenv LD_LIBRARY_PATH /swshare/glite/d-cache/dcap/lib/:${LD_LIBRARY_PATH} cvs co HeavyFlavorAnalysis/Bs2MuMu cvs co AnalysisDataFormats/HeavyFlavorObjects cd AnalysisDataFormats/HeavyFlavorObjects/ make cd ../../HeavyFlavorAnalysis/Bs2MuMu/macros/ make bin/runTreeReaders -r bmmReader -c chains/100407.mc -C cuts/bmmReader.default.cuts (for MinBias MC) bin/runTreeReaders -r bmmReader -c chains/100407.data -C cuts/bmmReader.default.cuts (for prompt reco data) </pre> You can redirect the output rootfile with -D ~/somewhere/. If you are using a bash, you'll need to translate the 'setenv LD_LIBRART_PATH ...' command appropriately. runTreeReaders.cc is a wrapper executable that instantiates tree classes to be run and provides some command line parsing and other setup utilities. The class bmmReader is just an example for an analysis class. It inherits from treeReader01 (some of the stuff which needs virtually no editing is in treeReader01.icc). treeReader01::goodRun() returns true if the run is in a list of good runs as defined in the runregistry. Lumisections are not yet included. When you produce a class (presumably derived from treeReader01 or even from bmmReader), and if it is needed for the final analysis (validation, normalization, whatever), please add it to runTreeReaders.cc. ---++Running in batch When your chains are very large, it is much better to split them into chunks and run in parallel over all chunks in batch <pre> cd chains/batch1 run -c ../../tree.csh -t ../../../../../grid-100427.tar.gz -m batch -x 'bin/runTreeReaders -n 900 -r muCharmReader -C cuts/muCharm.default.cuts' -r 'STORAGE1 srm://t3se01.psi.ch:8443/srm/managerv2\?SFN=/pnfs/psi.ch/cms/trivcat/store/user/ursl/root/trees/batch1' chain* </pre> Note: The output filename is determined by the chain name. If you want to submit batches with different cut files, you have two options: * store the output into different SE locations * name the chain files according to the different cut scenarios Other useful stuff for splitting chains and then combining the results: * you can split the chains with <pre> split -d -l 10 100414.v05.data b100414.v05.data- </pre> * you can combine the split output files with <pre> $ROOTSYS/bin/hadd combined.root \ dcap://t3se01.psi.ch:22125/pnfs/psi.ch/cms/trivcat/store/user/starodumov/d0/b100414.v05.data-00.root \ dcap://t3se01.psi.ch:22125/pnfs/psi.ch/cms/trivcat/store/user/starodumov/d0/b100414.v05.data-01.root \ </pre> ---++Skimming the large ntuples If you want to 'skim' the very large ntuples (e.g. because you are looking at candidates that only appear in a small fraction of all events), have a look at HeavyFlavorAnalysis/Bs2MuMu/macros/copyReader.cc. It will copy the ENTIRE event, given the selection in eventProcessing(). Please modify that according to your needs. -- Main.UrsLangenegger - 2010-06-30
E
dit
|
A
ttach
|
Watch
|
P
rint version
|
H
istory
: r6
<
r5
<
r4
<
r3
<
r2
|
B
acklinks
|
V
iew topic
|
Ra
w
edit
|
M
ore topic actions
Topic revision: r6 - 2010-06-30
-
UrsLangenegger
CmsTier3
Log In
CmsTier3 Web
Create New Topic
Index
Search
Changes
Notifications
Statistics
Preferences
User Pages
Main Page
Policies
Monitoring Storage Space
Monitoring Slurm Usage
Physics Groups
Steering Board Meetings
Admin Pages
AdminArea
Cluster Specs
Home
Site map
CmsTier3 web
LCGTier2 web
PhaseC web
Main web
Sandbox web
TWiki web
CmsTier3 Web
Create New Topic
Index
Search
Changes
Notifications
RSS Feed
Statistics
Preferences
P
P
P
View
Raw View
Print version
Find backlinks
History
More topic actions
Edit
Raw edit
Attach file or image
Edit topic preference settings
Set new parent
More topic actions
Account
Log In
E
dit
A
ttach
Copyright © 2008-2024 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki?
Send feedback