H to ZZ to 2L2NU

To Do List

Find which is the best Z+jets sample to use
Figure out the best discriminator values for the Track Counting algorithm (b-tagging)
Look into Nicola's 3D impact parameter stuff

Finished Items

Write Filter for only 2L2NU events
Run Some Data then test on filter
Running -- Create MonteCarlo Samples
Write HiggsToZZ2L2NU Package
--Create Pat Tuples to work with
--Reconstruct Z from Leptons
--Reconstruct Z from MET
https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideCandCombiner
--Write all analysis channels together
--Do MonteCarlo truth matching as we reconstruct
https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideCandCombiner

MonteCarlo

This file montecarlorpp.pdf is a PDF of pdgrid numbers and particle for use with the

pdgId() ==

command in c++.

List of Raw Data

115 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-d162eccdcfdb1619abbb00d2b4491240/USER
150 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-5b76a4d5ac2a2cd6e57a0f523c45e406/USER
200 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-42a0be463fd39c34af361ad3d112f85e/USER
250 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-0254735b2e2c307743d75b4d0c2f62b2/USER
300 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-c250a291024dcec21ed90be75d8781b0/USER
350 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-c1f12ad2dd1495d15bccd77798a9fc17/USER
400 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-5aaf0377d562f60d2832155e4dcd2517/USER
450 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-044047b2888a40223b7b58386a019f9a/USER
500 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-e0534a624e43895c0bd2c1775fed89dd/USER
550 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-395e480cb253c1e5855ff41bf03e274d/USER
600 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-ec4304b32e31125d5a2de76af031c2af/USER
650 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-af5479174a7832a91d27c8c62c42bb66/USER
700 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-4bad27a9b6a1d11820ddc53d559eb065/USER
750 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-2L2NU_GenSimRaw_CMSSW_3_7_0-aceba2100b53a5be308249da88a3f68b/USER

List of Reconstructed Data

115 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H115_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
150 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H150_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
200 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H200_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
250 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H250_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
300 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H300_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
350 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H350_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
400 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H400_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
450 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H450_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
500 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H500_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
550 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H550_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
600 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H600_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
650 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H650_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
700 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H700_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER
750 /2L2NU_GenSimRaw_CMSSW_3_7_0/mkress-H750_ZZ_2L2NU_7TeV_RECO_CMSSW_3_7_0-c267de12380b5c212f4e4c7c0229d1b5/USER

Interesting

| awk '{if($1>0 && $2~/Done/)printf "%s,",$1}'
awk '{printf("%s,",$1)}'

So the MonteCarlo is running and things seem to be going well. I have checked one file that finished with 6000 events this is the breakdown 2669 (44.4833%) passed the filter. Of those 875 Tau 899 Mu 895 e 916 Nu_Tau 895 Nu_mu 858 nu_e

For lxplus

source /afs/cern.ch/cms/LCG/LCG-2/UI/cms_ui_env.csh
source /afs/cern.ch/cms/ccs/wm/scripts/Crab/CRAB_2_7_2/crab.csh
cd /afs/cern.ch/user/m/mkress/scratch0/2L2NU/CMSSW_3_7_0/src/
cmsenv
cd GeneratorInterface/Pythia6Interface/test/

For cmslpc

source /uscmst1/prod/sw/cms/cshrc prod
source /uscmst1/prod/grid/gLite_SL5.csh
cd /uscms_data/d2/mkress/CMSSW/CMSSW_3_7_0/src
eval `scramv1 runtime -csh`
source /uscmst1/prod/grid/CRAB/crab.csh
cd HiggsAnalysis/HiggsToZZ4Leptons/test/
voms-proxy-init -voms cms

OLD
source /uscmst1/prod/sw/cms/cshrc prod
cd /uscms_data/d2/mkress/CMSSW/CMSSW_3_7_0/src
eval `scramv1 runtime -csh`
source /uscmst1/prod/grid/gLite_SL5_CRAB_27x.csh      
source /uscmst1/prod/grid/CRAB_2_7_3/crab.csh
cd GeneratorInterface/Pythia6Interface/test/
voms-proxy-init -voms cms

Log

Hi Matthew, yes, it means that you have now raw data in your hands and you need to run the reconstruction on top of it. You could either use the official cmsDriver command, like this:

cmsDriver.py step2_DT1_1 -s RAW2DIGI,RECO -n 1000 --eventcontent RECO --conditions FrontierConditions_GlobalTag,CRAFT09_R_V9::All --geometry Ideal --filein /your/path/your_input_file.root

where 1000 is the number of events you want to reconstruct in this file and CRAFT09_R_V9::All needs to be replaced with the correct global tag, which mainly depends on which release you are running in and in which release you simulated the input data. You can find all the different global tags here: https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideFrontierConditions make sure to take the right one for MC.

If this does not seem to work, please ask Nicola for what he did with the other files for reco.

Petra

good site to search for cmssw stuff http://cmslxr.fnal.gov