Analysis of the channel H → WW → lνl jj

Technical aspects

Recipe for Ntuples production

setenv SCRAM_ARCH slc5_amd64_gcc462
cmsrel CMSSW_5_2_5
cd CMSSW_5_2_5/src
cmsenv
cvs co -d PhysicsTools/NtupleUtils UserCode/Bicocca/PhysicsTools/NtupleUtils
cvs co -d PhysicsTools/MiBiCommonPAT UserCode/DiMatteo/MiBiCommonPAT
cvs co -r V00-02-05 -d CMGTools/External UserCode/CMG/CMGTools/External
scramv1 b -j9

Set up the environment

The analysis package is independent from the CMS software. To set up the working area, do as following:

export CVS_RSH=ssh
export CVSROOT=:ext:USERNAME@cmscvs.cern.ch:/cvs_server/repositories/CMSSW
cvs co -d NtuplePackage UserCode/abenagli/Fall10/NtuplePackage
cvs co -d KinFitter UserCode/abenagli/Fall10/KinFitter
cvs co -d VBFAnalysisPackage UserCode/abenagli/Fall10/VBFAnalysisPackage

Here the link to the cvs repository of the NtuplePackage.
Here the link to the cvs repository of the VBFAnalysisPackage.

Instructions to compile the package

The !NtuplePackage contains general functions and libraries (the ntuple reader!), necessary to run the analysis code. To compile it, do as following:

cd NtuplePackage
source scripts/setup.(c)sh
make
In general, once these libraries have been compiled for the first, there is no need to access this directory anymore.

The !KinFitter contains general functions and libraries, necessary to use the KinFit algorithms. To compile it, do as following:

cd KinFitter
source scripts/setup.(c)sh
make
In general, once these libraries have been compiled for the first, there is no need to access this directory anymore.

The !VBFAnalysisPackage contains the actual analysis code. To compile the functions and the executables, do as following:

cd VBFAnalysisPackage
source scripts/setup.(c)sh
make
make exe

Analysis @ 7 TeV

The current version of the analysis is based on the official Fall10 MonteCarlo production.
For signal samples, the cross-section and branching-ratio values are taken from the LHC Higgs Cross Section Working Group

mH (GeV/c2) σNLO(pp → qqH) (pb) BR(H → WW → lν qq) σNLO(pp → qqH, H → WW → lν qq) (pb)
120 1.27E+00 4.18E-02 5.30E-02
130 1.15E+00 8.94E-02 1.03E-01
140 1.05E+00 1.48E-01 1.55E-01
150 9.62E-01 2.04E-01 1.96E-01
160 8.79E-01 2.66E-01 2.34E-01
170 8.17E-01 2.82E-01 2.30E-01
180 7.48E-01 2.74E-01 2.05E-01
190 6.93E-01 2.30E-01 1.59E-01
200 6.37E-01 2.16E-01 1.38E-01
250 4.30E-01 2.06E-01 8.87E-02
300 3.01E-01 2.02E-01 6.08E-02
350 0.00E-01 2.02E-01 0.00E-02
400 1.62E-01 1.98E-01 3.21E-02
450 1.24E-01 1.91E-01 2.35E-02
500 9.49E-02 1.84E-01 1.75E-02
550 7.36E-02 1.78E-01 1.31E-02
600 5.76E-02 1.74E-01 1.00E-02

Cross-section values for background samples are taken from ...

Sample Cross-section (pb) Number of events L (pb -1 )
Single boson + jets
/WJets-madgraph/Summer09-MC_31X_V3-v1/ 36 600 10 009 466 250.237
/ZJets-madgraph/Summer09-MC_31X_V3-v1/ 3 600 1 163 479 314.454
/AstarJets-madgraph/Summer09-MC_31X_V3-v1/ 376 114485 197.388
DiBoson + jets
/VVJets-madgraph/Summer09-MC_31X_V3-v1/ 8 97719 8281.27
/WW/Summer09-MC_31X_V3-v1/GEN-SIM-RECO 44.4 203591 2747.52
/ZZ/Summer09-MC_31X_V3-v1/ 7.3 200564 19101.3
/WZ/Summer09-MC_31X_V3-v1/ 17.6 246550 7609.57
Top + jets
/TTbarJets-madgraph/Summer09-MC_31X_V3-v2/ 288 1 028 322 2 742.19
/TTbar-mcatnlo/Summer09-MC_31X_V3-v1/ 375 1 028 322 2 742.19
/SingleTop_sChannel-madgraph/Summer09-MC_31X_V3-v2/ 1.7 11999 2399.8
/SingleTop_tChannel-madgraph/Summer09-MC_31X_V3-v1/ 41.6 122936 4239.17
/SingleTop_tWChannel-madgraph/Summer09-MC_31X_V3-v3/ 27.6 239949 2121.56
QCD
/QCD_Pt15/Summer09-MC_31X_V3-v1/ 1.458 e+09 603 910 0.00041
/QCD_Pt30/Summer09-MC_31X_V3-v1/ 1.090 e+08 603 910 0.00041
/QCD_Pt80/Summer09-MC_31X_V3-v1/ 1.936 e+06 603 910 0.00041
/QCD_Pt170/Summer09-MC_31X_V3-v1/ 6.251 e+04 603 910 0.00041
/QCD_Pt300/Summer09-MC_31X_V3-v1/ 3669 603 910 0.00041
/QCD_Pt470/Summer09-MC_31X_V3-v1/ 351.3 603 910 0.00041
/QCD_Pt800/Summer09-MC_31X_V3-v1/ 11.9 603 910 0.00041
/QCD_Pt1400/Summer09-MC_31X_V3-v1/ 0.2 603 910 0.00041
/QCD_Pt2200/Summer09-MC_31X_V3-v1/ 1.422 e-03 603 910 0.00041
/QCD_Pt3000/Summer09-MC_31X_V3-v1/ 8.595 e-06 603 910 0.00041

VBF Preselection

A first preselection of the events is performed on signal and background samples. This preselection is intended to be in common with all VBF channels. The event requirements are as follow:

  • At least one lepton (e/mu) with pT > 10 GeV/c, loosely isolated Sum pT trakcs / pT < 0.5 and, if electron, passing the RobustLoose eID
  • At least two jets with eT > 15 GeV (if not corrected) or eT > 30 GeV (if L2L3 corrected)

Events passing these first requirements are saved in a root ntuple. Here the instructions:

INSTRUCTIONS FOR SIGNAL SAMPLES:

Go to folder HiggsAnalysis/VBFHiggsToVV/scripts/signal.
Change parameters in file params_lxbatch.CFG according to your needs.
The production of the ntuples is sent to the CERN batch system typing perl launchJobs_lxbatch.pl params_lxbatch.CFG.

INSTRUCTIONS FOR BACKGROUND SAMPLES:

Since background samples are available on the grid, crab is used to launch the ntuple production. Go to folder HiggsAnalysis/VBFHiggsToVV/scripts/background.
Specify the samples you want to run over in file multicrab.cfg and edit the user-dependent parameters in file crab_background.cfg.
Launch the production of ntuples typing multicrab -create -submit -cfg multicrab.cfg

LOCATION:

In both cases, ntuples are named VBFPreselection_%d.root and are stored on hercules, in folder

/gwteray/users/benaglia/NTUPLES/Fall10/

VBFHiggsToWWTolnujjPreselection

A second, more stringent preselection of the events is performed on signal and background samples. This preselection is intended to be specific for the VBF H -> WW -> lnu jj channel. The event requirements are as follow:

  • at least 1 lepton (e/mu) with pT > 10 GeV/c, loosely isolated Sum pT trakcs / pT < 0.5 and, if electron, passing the RobustLoose eID
  • at least 4 jets with eT > 15 GeV (if not corrected) or eT > 30 GeV (if L2L3 corrected)
  • at least a couple of jets with Deta > 1.5 and mJJ > 250 GeV/c2

INSTRUCTIONS

This second preselection is run outside CMSSW. Just go to the base directory of the VBFAnalysisPackage and set the necessary environement variables by typing source scripts/setup.(c)sh.

Create the list of samples to run over by typing the command rfdir $CASTOR/NTUPLES/VBF/Summer09_10TeV | awk '{print $9}' > data/listOfSamples.txt. A file called listOfSamples.txt is created under data/ and contains the list of the names of all castor folders.

To launch the ntuple production, go to folder scripts/VBFHiggsToWWTolnujjPreselection, modify the file params_lxbatch.CFG as required and then perl launchJobs_lxbatch.pl params_lxbathc.CFG

In both cases, ntuples are named VBFHiggsToWWTolnujjPreselection_%d.root and are stored on castor, in folder

/castor/cern.ch/user/a/abenagli/NTUPLES/VBF/Summer09_10TeV/sample_name

Note that, at this stage, ntuple size should be small enough to allow for a disk copy of all samples. To automatically launch the copy of the ntuples in /tmp/abenagli just do perl launchCopy.pl params_lxbatch.CFG.

-- AndreaBenaglia - 21-Nov-2009

Edit | Attach | Watch | Print version | History: r17 < r16 < r15 < r14 < r13 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r17 - 2012-05-22 - LeonardoDiMatteo
 
    • 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-2024 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