manTree SFrame PileUpTool instructions


Start with a clean shell & setup sframe.

Then you need to checkout the package:

svn co svn+ssh:// PileupReweighting
You will need if your cern username is different to your manchester one. If you want to take a tagged version, e.g. 00-00-09, the command is like:
svn co svn+ssh:// PileupReweighting

Now you need to compile (after patching the makefile):

cd PileupReweighting/cmt
gmake -f Makefile.Standalone

This generates the shared library: PileupReweighting/StandAlone/

Now make a soft link of library to your SFrame library folder

ln -s $SFRAME_DIR/PileupReweighting/StandAlone/ $SFRAME_DIR/lib

What The tool does

More details can be found here. But the tool compares the average number of interactions per bunch crossing (mu) for the MC sample (mc10a (r2215_r2260)) to what we see in data.

* AvintprxperiodD:

  • Mu_mc10a.png:

Adding PileUpWeight to SFrame analysis code

This code is dealt with by topUtils. You need to check:

The following extract is from TobObjectSel.cxx

m_tPileUp = new Root::TPileupReweighting( (((std::string)"T")+"Pileup").c_str() );
int isGood = m_tPileUp->initialize( sframedir + "/topUtils/data/PileupMu/ilumicalc_histograms_EF_mu18_178044-180242.root",
                                      sframedir+ "/topUtils/data/PileupMu/mu_mc10a.root",

Since Mc10a mu is stored in eventinfo class as LumiBlock(). To use this In your analysis code when running on MC do something like:

float mu = (float)m_evinfo->LumiBlock();
const double pileupEventWeight = m_object_sel->getPileUpWeight(mu);
weight*= pileupEventWeight;

Add the library to the config XML file topSel/config/TopSingleEle_config.xml

<JobConfiguration JobName="TopSingleEleJob" OutputLevel="INFO">
  <Library Name="PileupReweightingLib"/>

-- JohnAlmond - 05-May-2011

Topic attachments
I Attachment History Action Size Date Who Comment
PNGpng Mu_mc10a.png r1 manage 9.1 K 2011-05-06 - 01:36 JohnAlmond  
PNGpng avintprxperiodD.png r1 manage 9.0 K 2011-05-06 - 01:13 JohnAlmond AvintprxperiodD
Edit | Attach | Watch | Print version | History: r6 < r5 < r4 < r3 < r2 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r6 - 2011-05-12 - MarkOwen
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Atlas All webs login

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