LPC Organization
Physics Topologies:
Lepton + Jets + MET = LJMET
Getting started with the Starter Kit
(instructions from Jia Fu Low)
================================
README for CMSSW_1_6_12
================================
1. Be sure that you have done the preliminary setup for CMSSW when you log in;
on the cmslpc cluster, this is `source /uscmst1/prod/sw/cms/setup/cshrc prod'
2. Create a new CMSSW working area for yourself via `cmsrel CMSSW_1_6_12', and
cd to the CMSSW_1_6_12/src area.
3. Type `cmsenv' to establish your environment variables
4. Follow the instructions at
https://twiki.cern.ch/twiki/bin/view/CMSWorkBookComputingConcepts#CvsCheckoutModules
to properly set up CVS access. Anonymous access will work if you only want to
check out code and not check in updates later.
5. Check out the codes with CVS from UserCode/jiafulow via
`cvs co -d PhysicsTools UserCode/jiafulow/PhysicsTools`. The code should now
be in the area CMSSW_1_6_12/src/PhysicsTools.
6. `cd PhysicsTools/StarterKit' and `scramv1 b' to build the executable code.
7. The main config file is the `topSemiLepMuon.cfg' (in PhysicsTools/StarterKit).
`cmsRun topSemiLepMuon.cfg >& log' will produce a `histos.root' which contains
the StarterKit histograms as well as additional histograms in `MyHistos'.
You may need to edit the config file to specify the input file by uncommenting
the include line for the file that specifies the input chowder sample, and you
may want to also adjust the maximum number of events processed.
(Note that the default input file is
/CSA07Muon/CMSSW_1_6_7-CSA07-Chowder-P1-PDMuon-Skims4-topSemiLepMuon/USER)
8. The second config file `topSemiLepMuonQCD.cfg' is used to run on QCD event
datasets that are not part of the CSA07 soups.
(Note that the default QCD input file is
/ppMuPt20-15/CMSSW_1_6_7-CSA07-1205151815/AODSIM)
9. In both config files, several paths are already defined. For each path,
selection cuts are applied according to the Analysis Note TOP08-05.
A complete set of plots are made in each path.
a. Pre : the default PAT outputs, without any cut
(Preselection in TOP08-05)
b. LooseMuon : Pre + Loose Muon
c. JetBin4 : Pre + Inclusive 4 Jets
(4 Jets pT > 65/40/40/40 GeV in TOP08-05)
d. Bin4Loose : Pre + Loose Muon + Inclusive 4 Jets
(1 Muon pT > 30 GeV in TOP08-05)
e. Bin1 : Pre + Tight Muon + Exclusive 1 Jet
f. Bin2 : Pre + Tight Muon + Exclusive 2 Jets
g. Bin3 : Pre + Tight Muon + Exclusive 3 Jets
h. Bin4 : Pre + Tight Muon + Inclusive 4 Jets
(dRmin > 0.3 & Ecalo < 1 GeV cut in TOP08-05)
-> Loose Muon: leading muon pT > 30 GeV, |eta| < 2.1, trackIso < 3 GeV,
caloIso < 5 GeV
other muon pT < 30 GeV
-> Tight Muon: Loose Muon
+ leading muon caloIso < 1 GeV
+ deltaR between leading muon and any jet > 0.3
-> Jet : leading jet pT > 65 GeV, other jet pT > 40 GeV. All jets |eta| < 2.4
10. The following modules are called directly:
a. CSA07EventWeightProducer
b. The standard patLayer0 and patLayer1 sequences (in PhysicsTools/PatAlgos/data)
c. CSA07ProcessIdFilter
d. TopSemiLepKit (extra histograms are made here)
The parameters given to these modules can be found in PhysicsTools/StarterKit/data
Using these parameters, the results obtained are similar to those in TOP08-05.
11. TopSemiLepKit inherits from the StarterKit so it gets all the PAT outputs and
the standard StarterKit histograms. In addition, it calls `MyHistos' class to
add more histograms (in the same way as StarterKit calling PhysicsHistograms
class to make histograms).
12. In `MyHistos' class, different kinds of histograms (PhysVarHistos) are made:
a. di-jet deltaR, invariant mass
b. jet maximum |eta|, multiplicity, minimum pT
c. deltaPhi between MET & leading jet, MET & leading muon
d. deltaR between leading muon & leading jet
e. W transverse mass from muon and MET
f. Z invariant mass from two muons
g. event shape variables, e.g. sphericity, aplanarity, centrality, HT
h. some extra plots which might or might not be useful
13. In calculating the event shape variables, three different cases are used:
a. exactly four jets
b. all jets with pT > 40 GeV & |eta| < 2.1
c. all jets with pT > 20 GeV
For each case, the variables are calculated using
i. jets only
ii. muon + jets (default)
iii. muon + MET + jets
14. The macros `plotMacro.C' and `plotMacro2.C' (more organized) is used to plot
the histograms and save them in .png format.
15. All suggestions and comments are welcome!
================================
Results
================================
| ttbar | W | Z | QCD
----------------------------+---------------+---------------+----------------+------------
Preselection | 1289.6 (1276) | 7674.8 (7474) | 1468.1 (1430) |
| | | |
4 Jet pT > 65/40/40/40 GeV | 401.2 (371) | 87.2 (83) | 17.4 (16) |
| | | |
1 Muon pT > 30 GeV | 200.6 (195) | 58.7 (57) | 7.0 (8) | 125.3 (125)
| | | |
Rmin > 0.3 & Ecalo < 1 GeV | 156.5 (153) | 45.6 (45) | 5.5 (7) | 11.5 (11)
* The numbers were obtained when I used the codes to run over all events in the
two datasets. The numbers in bracket are the numbers quoted in TOP08-05.
================================
README for CMSSW_2_1_0
================================
1. To get the code and compile it:
# Set up a new project area
scramv1 p CMSSW CMSSW_2_1_0_pre8
cd CMSSW_2_1_0_pre8/src
# Set up CVS
cmscvsroot
cvs login
(password 98passwd)
# Get the pat code
cvs co -r V03-08-00 DataFormats/PatCandidates
cvs co -r V04-03-02 PhysicsTools/PatAlgos
cvs co -r V03-01-00 PhysicsTools/PatUtils
# Get the starter kit code
cvs co -r V02-02-00 PhysicsTools/StarterKit
# Get the top semileptonic kit code
cvs co -d Analysis/TopSemiLepKit -A UserCode/jiafulow/Analysis/TopSemiLepKit
# compile
scramv1 b -j 4
# setup runtime paths
cmsenv
# run
cd Analysis/TopSemiLepKit/test
cmsRun topSemiLepMuon_cfg.py
2. The main config file is the `topSemiLepMuon.cfg' (in PhysicsTools/StarterKit).
`cmsRun topSemiLepMuon.cfg >& log' will produce a `histos.root' which contains
the StarterKit histograms as well as additional histograms from TopSemiLepKit.
You can change the PoolSource to be whatever you'd like, by default it is
/RelValTTbar/CMSSW_2_1_0_pre6-RelVal-1214239099-STARTUP_V1-2nd/GEN-SIM-DIGI-RAW-HLTDEBUG-RECO
3. Several paths are already defined. For each path,
selection cuts are applied according to the Analysis Note TOP08-05.
A complete set of plots are made in each path.
a. Pre : the default PAT outputs, without any cut
(Preselection in TOP08-05)
b. LooseMuon : Pre + Loose Muon
c. JetBin4 : Pre + Inclusive 4 Jets
(4 Jets pT > 65/40/40/40 GeV in TOP08-05)
d. Bin4Loose : Pre + Loose Muon + Inclusive 4 Jets
(1 Muon pT > 30 GeV in TOP08-05)
e. Bin1 : Pre + Tight Muon + Exclusive 1 Jet
f. Bin2 : Pre + Tight Muon + Exclusive 2 Jets
g. Bin3 : Pre + Tight Muon + Exclusive 3 Jets
h. Bin4 : Pre + Tight Muon + Inclusive 4 Jets
(dRmin > 0.3 & Ecalo < 1 GeV cut in TOP08-05)
-> Loose Muon: leading muon pT > 30 GeV, |eta| < 2.1, trackIso < 3 GeV,
caloIso < 5 GeV
other muon pT < 30 GeV
-> Tight Muon: Loose Muon
+ leading muon caloIso < 1 GeV
+ deltaR between leading muon and any jet > 0.3
-> Jet : leading jet pT > 65 GeV, other jet pT > 40 GeV. All jets |eta| < 2.4
4. TopSemiLepKit uses PatAnalyzerKitHelper, so it gets all the PAT outputs and
the standard PatAnalyzerKit histograms. In addition, it calls `TopSemiLepKitHelperHistos' class to
add more histograms (in the same way as PatAnalyzerKit calling PatAnalyzerKitHelper
class to make histograms).
5. In `TopSemiLepKitHelperHistos' class, different kinds of histograms (PhysVarHistos) are made:
a. di-jet deltaR, invariant mass
b. jet maximum |eta|, multiplicity, minimum pT
c. deltaPhi between MET & leading jet, MET & leading muon
d. deltaR between leading muon & leading jet
e. W transverse mass from muon and MET
f. Z invariant mass from two muons
g. event shape variables, e.g. sphericity, aplanarity, centrality, HT
h. some extra plots which might or might not be useful
6. In calculating the event shape variables, three different cases are used:
a. exactly four jets
b. all jets with pT > 40 GeV & |eta| < 2.1
c. all jets with pT > 20 GeV
For each case, the variables are calculated using
i. jets only
ii. muon + jets (default)
iii. muon + MET + jets
7. The macros `plotMacro.C' and `plotMacro2.C' (more organized) is used to plot
the histograms and save them in .png format.
