Some useful notes about luminosity and pileup determinations.


<p><a href="http://cms.cern.ch/iCMS/analysisadmin/cadilines?line=LUM-12-002&tp=an&id=920&ancode=LUM-12-002">LUM-12-002 (unpublished)</a>

<p><h2>Pileup task</h2>


<p><i>With condor:</i>

<p><a href="http://batchdocs.web.cern.ch/batchdocs/local/quick.html">Quick start</a>

<p>condor_submit createPileupJSON.sub

<p><i>With old batch system:</i>

<p>~giamman/public/pu/myPileupScript.sh [label] [input JSON] [normtag] [xsection] [nBins] [local directory]
<br>~giamman/public/pu/createPileupJSON.sh

<p>bsub -q 1nw "bash createPileupJSON.sh"

<p>To do list:
<ul>
<li>lumipro area: establish backup policy (keep last two?)</li>
<li>PR on github</li>
<li>note: I need to edit PublicPlots/.git/config and "unset SSH_ASKPASS"</li>
<li>cron job</li>
<ol>
<li>in my area</li>
<li>then in lumipro</li>
<li>using "acrontab -e"</li>
<li><a href="https://kb.iu.edu/d/afiz">useful reference for syntax</a></li>
</ol>
<li>create public plot as in <a href="https://twiki.cern.ch/twiki/bin/view/CMSPublic/LumiPublicResults#Pileup_distribution">https://twiki.cern.ch/twiki/bin/view/CMSPublic/LumiPublicResults#Pileup_distribution</a></li>
<ul>
<li>cd /afs/cern.ch/user/l/lumipro/PublicPlots</li>
<li>ls -trl *pile*</li>
<li>Script: create_public_pileup_plots.py</li>
<li>Its config files: public_pileup_plots_pp_201*.cfg</li>
<li>python create_public_pileup_plots.py public_pileup_plots_pp_2017.cfg --ignore-cache</li>
</ul>
<li>Make script "incremental" to save time at each execution</li>
</ul>

<p><hr>

<p><h3>Create pileup JSON</h3>

<p><a href="http://cms-service-lumi.web.cern.ch/cms-service-lumi/brilwsdoc.html#brilcalclumi">brilcalclumi</a>

<p>Install:
<pre>
bash
export PATH=$HOME/.local/bin:/afs/cern.ch/cms/lumi/brilconda-1.0.3/bin:$PATH
pip install --install-option="--prefix=$HOME/.local" brilws
</pre>
<p>Create csv:
<pre>
brilcalc lumi --xing -i /afs/cern.ch/cms/CAF/CMSCOMM/COMM_DQM/certification/Collisions16/13TeV/DCSOnly/json_DCSONLY.txt -o /afs/cern.ch/work/g/giamman/pu/pileup.csv --normtag /afs/cern.ch/user/c/cmsbril/public/2016normtags/normtag_BRIL.json
</pre>

<p><a href="https://twiki.cern.ch/twiki/bin/viewauth/CMS/PileupJSONFileforData#Calculating_the_Relevant_Pileup">PU calculation</a>

<p>Create JSON (from inside a CMSSW area):
<pre>
estimatePileup_makeJSON_2015.py --csvInput /afs/cern.ch/work/g/giamman/pu/pileup.csv /afs/cern.ch/work/g/giamman/pu/pileup_JSON.txt
</pre>

<p><a href="https://github.com/cms-sw/cmssw/tree/CMSSW_8_0_X/RecoLuminosity/LumiDB/scripts">PU scripts in CMSSW</a>

<!---
<p><a href="https://cern.service-now.com/service-portal/article.do?n=KB0001998">EOS instructions</a>
<pre>
export EOS_MGM_URL=root://eoscms.cern.ch
eos ls -l /eos/cms/user/g/giamman
</pre>
--->

<p><hr>

<p><a href="https://twiki.cern.ch/twiki/bin/viewauth/CMS/PileupMCReweightingUtilities">https://twiki.cern.ch/twiki/bin/viewauth/CMS/PileupMCReweightingUtilities</a>
<br>
<pre>
Handle<std::vector< PileupSummaryInfo > >  PupInfo;
event.getByLabel(edm::InputTag("addPileupInfo"), PupInfo);

std::vector<PileupSummaryInfo>::const_iterator PVI;

float Tnpv = -1;
for(PVI = PupInfo->begin(); PVI != PupInfo->end(); ++PVI) {

   int BX = PVI->getBunchCrossing();

   if(BX == 0) { 
     Tnpv = PVI->getTrueNumInteractions();
     continue;
   }

}

double MyWeight = LumiWeights_.weight( Tnpv );
</pre>

<p>Use stand-alone: check PhysicsTools/Utilities/interface/LumiReweightingStandAlone.h
<br>See <a href="https://twiki.cern.ch/twiki/bin/viewauth/CMS/PileupMCReweightingUtilities">https://twiki.cern.ch/twiki/bin/viewauth/CMS/PileupMCReweightingUtilities</a>

<p>Use in EDAnalyzer (from <a href="https://raw.githubusercontent.com/nfaltermann/ST_RunII_EA/fca1a11ee568d814afdbf8a0e58df5add5196523/src/DMAnalysisTreeMaker.cc">here</a>):
<pre>
#include "PhysicsTools/Utilities/interface/LumiReWeighting.h"

private:
  virtual void beginJob() {
    std::string distr = "pileUp" + dataPUFile_ + ".root";
    LumiWeights_ = edm::LumiReWeighting(distr,"DataPileupHistogram_69mbMinBias.root",std::string("pileup"),std::string("pileup"));	
    LumiWeightsUp_ = edm::LumiReWeighting(distr,"DataPileupHistogram_69mbMinBias_up.root",std::string("pileup"),std::string("pileup"));
    LumiWeightsDown_ = edm::LumiReWeighting(distr,"DataPileupHistogram_69mbMinBias_down.root",std::string("pileup"),std::string("pileup"));
  }

  edm::LumiReWeighting LumiWeights_, LumiWeightsUp_, LumiWeightsDown_;
</pre>

<pre>
  edm::InputTag lumiBlock_ = iConfig.getParameter<edm::InputTag>("lumiBlock");
  t_lumiBlock_ = consumes< unsigned int >( lumiBlock_ );
  edm::InputTag runNumber_ = iConfig.getParameter<edm::InputTag>("runNumber");
  t_runNumber_ = consumes< unsigned int >( runNumber_ );
  edm::InputTag eventNumber_ = iConfig.getParameter<edm::InputTag>("eventNumber");
  t_eventNumber_ = consumes< ULong64_t >( eventNumber_ );
</pre>

<pre>

void DMAnalysisTreeMaker::getPUSF(){
  puWeight=(float) LumiWeights_.weight(nTruePU);
  puWeightUp = (float) LumiWeightsUp_.weight(nTruePU);
  puWeightDown = (float) LumiWeightsDown_.weight(nTruePU);
  float_values["Event_puWeight"]=puWeight;
  float_values["Event_puWeightUp"]=puWeightUp; 
  float_values["Event_puWeightDown"]=puWeightDown;
  float_values["Event_nTruePU"]=(float)nTruePU;
}
</pre>

<p>How to use "consumes": <a href="https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideEDMGetDataFromEvent#consumes">link</a>

<p>List of all available collections: <a href="https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideRecoDataTable">https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideRecoDataTable</a>
<br>To be used:
<ul>
<li>offlinePrimaryVertices</li>
<li>offlinePrimaryVerticesWithBS (as cross check)</li>
<li>pixelVertices</li>
<li>pixelTracks</li>
<!---
<li></li>
<li></li>
<li></li>
<li></li>
<li></li>
<li></li>
--->
</ul>


<p><hr>

<p><h2>Length scale calibration</h2>

<p>Section 6.3 of LUM-12-002 (see above)
<p><a href="https://indico.cern.ch/event/502071/contributions/1179003/attachments/1231672/1805862/LSC_reviewSalfeld.pdf">talk #1</a>
<p><a href="https://indico.cern.ch/event/449981/contributions/1945821/attachments/1162182/1673806/LumiDPG_Sep292015.pdf">talk #2</a>

<p><hr>

<p><h2>Brilcalc</h2>

<p><a href="http://fwyzard.web.cern.ch/fwyzard/lumi/brilcalc">Andrea Bocci's wrapper script to ensure correct environment</a>