HZZ Exercise

HZZ->4l Search Exercise

2 Oct 2014

Doing Exercise 1, looking at plots for the Z1 peak, File PlotStack4l-davec.C
Input files: inputfile="filelist_4mu_2012_paper_CERN-davec.txt";
std::string histolabel = "hMZ1_3"; // Z1 mass.
Output to screen: PlotStack4l-davec-output-text.txt

Files actually used written out to "test-data-files-used.txt":

Data files used, *f1
./log/histos4mu/output_DoubleMu_Run2012A-22Jan2013-v1.root
./log/histos4mu/output_DoubleMu_Run2012B-22Jan2013-v1.root
./log/histos4mu/output_DoubleMu_Run2012B-22Jan2013-v1_1.root
./log/histos4mu/output_DoubleMu_Run2012C-22Jan2013-v1_1.root
./log/histos4mu/output_DoubleMu_Run2012D-22Jan2013-v1.root
./log/histos4mu/output_DoubleMu_Run2012D-22Jan2013-v1_2.root
Background files used, *f2
./log/histos4mu/output_GluGluToZZTo4L_8TeV-gg2zz-pythia6.root
./log/histos4mu/output_ZZTo4mu_8TeV-powheg-pythia6.root
./log/histos4mu/output_ZZTo4e_8TeV-powheg-pythia6.root
./log/histos4mu/output_ZZTo2e2tau_8TeV-powheg-pythia6.root
./log/histos4mu/output_WWJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola.root
./log/histos4mu/output_Tbar_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola.root
./log/histos4mu/output_T_s-channel_TuneZ2star_8TeV-powheg-tauola.root
Files for Signal, *f0
./log/histos4mu/output_WH_HToZZTo4L_M-200_8TeV-pythia6.root
./log/histos4mu/output_WH_HToZZTo4L_M-126_8TeV-pythia6.root
./log/histos4mu/output_WH_HToZZTo4L_M-125_8TeV-pythia6.root
./log/histos4mu/output_ZH_HToZZTo4L_M-200_8TeV-pythia6.root
./log/histos4mu/output_VBF_HToZZTo4L_M-126_8TeV-powheg-pythia6.root
./log/histos4mu/output_VBF_HToZZTo4L_M-125_8TeV-powheg-pythia6.root
./log/histos4mu/output_SMHiggsToZZTo4L_M-126_8TeV-powheg15-JHUgenV3-pythia6.root
./log/histos4mu/output_SMHiggsToZZTo4L_M-125_8TeV-powheg15-JHUgenV3-pythia6.root
./log/histos4mu/output_GluGluToHToZZTo4L_M-800_8TeV-powheg15-pythia6.root

h_hMZ1_3_4mu_8TeV_ratio_log.png:

Data:
Data filelist from vector: "for (unsigned int datasetIdData=0; datasetIdData<Vdatasetnamedata.size(); datasetIdData++)"
Data accumulated through "htotaldata->Add(hfourlepbestmass_4l_afterSel_new_new);"
Data loaded into graph "gr" and new error bars calculated with "for (int nbins=1;nbins<=htotaldata->GetNbinsX(); nbins++){"
and with TGraphAsymmErrors *gr = new TGraphAsymmErrors(Nbins,x,y,exl,exh,eyl,eyh);
gr->SetMarkerColor(1); gr->SetMarkerStyle(20); gr->SetMarkerSize(0.95); at ~line 577
gr->Draw("EPsame"); at ~line 1375

Higgs at 125:
Signal files (histos), with titles ...HToZZTo4L... loaded one by one into hfourlepbestmass_4l_afterSel_new
with "for ( int datasetIdSig=Vdatasetnamesig.size()-1; datasetIdSig >=0; datasetIdSig--){"
Load up working histo: "hfourlepbestmass_4l_afterSel_new_new" after nRebin
Accumulate bincont125 from hfourlepbestmass_4l_afterSel_new_new if datasetname has M-125 (similarly for 126, 200, 350, 500, 800)
Load bincontgg125 and bincontvbf125 separately
Load up histo: hfourlepbestmass_4l_afterSel_new_signal125->SetBinContent(nbins,double(bincont125[nbins-1]));, line 1306
Draw with "hfourlepbestmass_4l_afterSel_new_signal125->Draw("same");" ~line 1429

Background:
Uses htotal which is a THStack collection of TH1 (or derived) objects. Each TH1 shown with its own colour etc.
Datasets loaded from: for ( int datasetId=Vdatasetnamebkg.size()-1; datasetId >=0; datasetId--)
Make hfourlepbestmass_4l_afterSel_new_ histos with DY, DYlight, DYbb, DYcc, WW, WZ, TT, Wj, ggZZ, ZZ, qqZZ, qcdMu, qcdEM, qcdBC, singlet, in title

WW
Load hfourlepbestmass_4l_afterSel_new_WW if WW in database name
if(datasetnamebkg.find("WW") < 80){hfourlepbestmass_4l_afterSel_new_WW->Add(hfourlepbestmass_4l_afterSel_new_new);, line 801
Accumulate into htotal: htotal->Add(hfourlepbestmass_4l_afterSel_new_WW);, line 1072
TOP
hfourlepbestmass_4l_afterSel_new_singlet->Add(hfourlepbestmass_4l_afterSel_new_new);, line 1003
Accumulate htotal: htotal->Add(hfourlepbestmass_4l_afterSel_new_singlet);
htotal->Draw("same");, line 1408

9 Sep 2014

New folders for root output files:

• /afs/cern.ch/work/d/davec/HZZ/root-output/2e2mu
• /afs/cern.ch/work/d/davec/HZZ/root-output/4e
• /afs/cern.ch/work/d/davec/HZZ/root-output/4mu

Folder defined in compilereference_4mu_data-davec.C (and similarly for 4e, 2e2mu)

// place root file in work area:
      path="/afs/cern.ch/work/d/davec/HZZ/root-output/4mu";
      cout << "Output directory for root file is " << path << endl;

Define type of analysis in HZZ4LeptonsAnalysis_4e.C (and similarly 4mu and 2e2mu) in the histo title:
TH1D *nEvent_4l = new TH1D("nEvent_4l", "nEventComplete - HZZ4LeptonsAnalysis_4e", 20, 0., 20.);

15 Aug 2014

2011 datasets from:

• data_input_4e_2011_paper.txt
• data_input_4mu_2011_paper.txt
• data_input_2e2mu_2011_paper.txt

Analysis	Data	Rootuble size	Output size	Events from roottree	Events after HZZ analysis	Comment
4mu	DoubleMu_Run2011A-16Jan2012-v1	48GB	120K	3,946,595	10	histos4mu, also done by me
4mu	DoubleMu_Run2011B-16Jan2012-v1	63GB	120K	4,449,454	9	histos4mu, also done by me
4e	DoubleElectron_Run2011A-16Jan2012-v1	53GB	120K	4,399,347	7	histos4e, 3Dec2013
4e	DoubleElectron_Run2011B-16Jan2012-v1	66GB	120K	4,729,123	8	histos4e, also done by me
2e2mu	DoubleMu_Run2011A-16Jan2012-v1	48GB	119K	3,946,595	3	Done by me on 8Sep2014
2e2mu	DoubleMu_Run2011B-16Jan2012-v1	63GB	119K	4,449,454	3	Done by me on 8Sep2014
2e2mu	DoubleElectron_Run2011A-16Jan2012-v1	53GB	118K	4,399,347	11	histos2e2mu, 3Dec2013
2e2mu	DoubleElectron_Run2011B-16Jan2012-v1	66GB	125K	4,729,123	27	Done by me on 8Sep2014
2e2mu	MuEG_Run2011A-13Dec2012-v1	100GB	103K	8,043,543	0	histos2e2mu, 3Dec2013
2e2mu	MuEG_Run2011A-13Dec2012-v1_1	100GB	103K	7,998,750	0	histos2e2mu, 3Dec2013
2e2mu	MuEG_Run2011A-13Dec2012-v1_2	28GB	112K	2,292,001	1	Done by me, 2011 OK
2e2mu	MuEG_Run2011B-13Dec2012-v1	100GB	114K	6,997,788	2	histos2e2mu, 3Dec2013
2e2mu	MuEG_Run2011B-13Dec2012-v1_1	58GB	115K	4,108,480	1	Done by me, 2011 OK

• histos4mu in /afs/cern.ch/user/d/davec/HZZ/CMSSW_5_3_9/src/HiggsAnalysis/HiggsToZZ4Leptons/test/macros/FSR2012/INCLUSIVE_MACROS/log/lustre/cms/store/user/defilip/histos4mu_nobnn
• histos4e in /afs/cern.ch/user/d/davec/HZZ/CMSSW_5_3_9/src/HiggsAnalysis/HiggsToZZ4Leptons/test/macros/FSR2012/INCLUSIVE_MACROS/log/lustre/cms/store/user/defilip/histos4e_nobnn
• histos2e2mu in /afs/cern.ch/user/d/davec/HZZ/CMSSW_5_3_9/src/HiggsAnalysis/HiggsToZZ4Leptons/test/macros/FSR2012/INCLUSIVE_MACROS/log/lustre/cms/store/user/defilip/histos2e2mu_nobnn
• My 4e histos in /afs/cern.ch/work/d/davec/HZZ/root-output/4e
• My 4mu histos in /afs/cern.ch/work/d/davec/HZZ/root-output/4mu
• My 2e2mu histos in /afs/cern.ch/work/d/davec/HZZ/root-output/2e2mu

14 Aug 2014

Using HZZ-root.C to analyse output files to give events from GB roottree file and events after HZZ analysis.

2012 datasets

• data_input_4mu_2012_paper.txt

Analysis	Data	Rootuble size	Output size	Events from roottree	Events after HZZ analysis	Comment
4mu	DoubleMu_Run2012A-22Jan2013-v1	29GB	117K	1,928,938	5	histos4mu.tgz, 3Dec2013
4mu	DoubleMu_Run2012B-22Jan2013-v1	100GB	123K	6,561,216	15	histos4mu.tgz, 3Dec2013
4mu	DoubleMu_Run2012B-22Jan2013-v1_1	100GB	119K	6,587,958	7	histos4mu.tgz, 3Dec2013
4mu	DoubleMu_Run2012B-22Jan2013-v1_2	18GB	106K	1,180,848	4	done by me, 9Sep2014
4mu	DoubleMu_Run2012C-22Jan2013-v1	100GB	125K	6,534,272	17	done by me, 9Sep2014
4mu	DoubleMu_Run2012C-22Jan2013-v1_1	100GB	126 (old 128K??)	6,565,407	24	histos4mu.tgz, 3Dec2013
4mu	DoubleMu_Run2012C-22Jan2013-v1_2	80GB	120K	5,242,375	8	done by me, 2Sep2014
4mu	DoubleMu_Run2012D-22Jan2013-v1	100GB	122 (old 124K??)	6,457,460	12	histos4mu.tgz, 3Dec2013
4mu	DoubleMu_Run2012D-22Jan2013-v1_1	100GB	123K	6,428,330	15	done by me, 2Sep2014
4mu	DoubleMu_Run2012D-22Jan2013-v1_2	87GB	120 (old 122K??)	5,614,444	8	histos4mu.tgz, 3Dec2013
4mu	DoubleMu_Run2012D-22Jan2013-v1_3	1GB	101K	59,617	0	done by me, 9Sep2014

• histos4mu.tgz expanded into
  
• /afs/cern.ch/user/d/davec/HZZ/CMSSW_5_3_9/src/HiggsAnalysis/HiggsToZZ4Leptons/test/macros/FSR2012/INCLUSIVE_MACROS/log/lustre/cms/store/user/defilip/histos4mu_nobnn
• My files in /afs/cern.ch/work/d/davec/HZZ/root-output/4mu

• data_input_2e2mu_2012_paper.txt

Analysis	Data	Rootuble size	Output size	Events from roottree	Events after HZZ analysis	Comment
2e2mu	DoubleElectron_Run2012A-22Jan2013-v1	23GB	116K	1,621,779	9	Done by me, 9Sep2014
2e2mu	DoubleElectron_Run2012B-22Jan2013-v1	100GB	126K	6,542,925	31	histos2e2mu, 3Dec2013
2e2mu	DoubleElectron_Run2012B-22Jan2013-v1_1	24GB	113K	1,647,962	2	Done by me, 9Sep2014
2e2mu	DoubleElectron_Run2012C-22Jan2013-v1	100GB	124K	6,540,600	23	Done by me, 9Sep2014
2e2mu	DoubleElectron_Run2012C-22Jan2013-v1_1	94GB	126K	6,175,001	27	Done by me, 9Sep2014
2e2mu	DoubleElectron_Run2012D-22Jan2013-v1	100GB	128K	6,378,311	35	histos2e2mu, 3Dec2013
2e2mu	DoubleElectron_Run2012D-22Jan2013-v1_1	50GB	123K	3,379,685	22	Done by me, 9Sep2014
2e2mu	DoubleMu_Run2012A-22Jan2013-v1	29GB	116K	1,928,938	1	histos2e2mu, 3Dec2013
2e2mu	DoubleMu_Run2012B-22Jan2013-v1	100GB	119K	6,561,216	2	Done by me, 9Sep2014
2e2mu	DoubleMu_Run2012B-22Jan2013-v1_1	100GB	121K	6,587,958	6	Done by me, 9Sep2014
2e2mu	DoubleMu_Run2012B-22Jan2013-v1_2	18GB	104K	1,180,848	0	histos2e2mu, 3Dec2013
2e2mu	DoubleMu_Run2012C-22Jan2013-v1	100GB	121K	6,534,272	6	histos2e2mu, 3Dec2013
2e2mu	DoubleMu_Run2012C-22Jan2013-v1_1	100GB	122K	6,565,407	9	histos2e2mu, 3Dec2013
2e2mu	DoubleMu_Run2012C-22Jan2013-v1_2	80GB	120K	5,242,375	5	Done by me, 9Sep2014
2e2mu	DoubleMu_Run2012D-22Jan2013-v1	100GB	122K	6,457,460	7	histos2e2mu, 3Dec2013
2e2mu	DoubleMu_Run2012D-22Jan2013-v1_1	100GB	120K	6,428,330	4	histos2e2mu, 3Dec2013
2e2mu	DoubleMu_Run2012D-22Jan2013-v1_2	87GB	121K	5,614,444	5	histos2e2mu, 3Dec2013
2e2mu	DoubleMu_Run2012D-22Jan2013-v1_3	1GB	99K	59,617	0	histos2e2mu, 3Dec2013
2e2mu	MuEG_Run2012A-22Jan2013-v1	11GB	100K	728,793	0	Done by me, 2Sep2014
2e2mu	MuEG_Run2012B-22Jan2013-v1	100GB	115K	6,523,647	3	Done by me, 2Sep2014
2e2mu	MuEG_Run2012B-22Jan2013-v1_1	29GB	101K	1,864,540	0	Done by me, , 2Sep2014, not in 2e2mu paper
2e2mu	MuEG_Run2012C-22Jan2013-v1	100GB	104K	6,514,014	0	histos2e2mu, 3Dec2013
2e2mu	MuEG_Run2012C-22Jan2013-v1_1	95GB	103K	6,171,564	0	Done by me, 2Sep2014
2e2mu	MuEG_Run2012D-22Jan2013-v1	100GB	105K	6,405,573	0	histos2e2mu, 3Dec2013
2e2mu	MuEG_Run2012D-22Jan2013-v1_1	100GB	92K	6,386,009	2	Done by me, 2Sep2014
2e2mu	MuEG_Run2012D-22Jan2013-v1_2	8GB	101K	527,859	0	histos2e2mu, 3Dec2013

• histos2e2mu in /afs/cern.ch/user/d/davec/HZZ/CMSSW_5_3_9/src/HiggsAnalysis/HiggsToZZ4Leptons/test/macros/FSR2012/INCLUSIVE_MACROS/log/lustre/cms/store/user/defilip/histos2e2mu_nobnn
• My files in /afs/cern.ch/work/d/davec/HZZ/root-output/2e2mu

9 Aug 2014

• data_input_4e_2012_paper.txt
• output root files in ~/HZZ/CMSSW_5_3_9/src/HiggsAnalysis/HiggsToZZ4Leptons/test/macros/FSR2012/INCLUSIVE_MACROS/log/lustre/cms/store/user/defilip/4e_nobnn/

Dataset	Castor	First level processed data source	Filename	Size	Events from roottree	Events after HZZ analysis
DoubleElectron_Run2012A-22Jan2013-v1	23GB	rootuple copied from Castor to "work", batch: run-4e-davec.sh	output_DoubleElectron_Run2012A-22Jan2013-v1.root	117K	1,621,779	3
DoubleElectron_Run2012B-22Jan2013-v1	100GB	From histos4e.tgz	output_DoubleElectron_Run2012B-22Jan2013-v1.root	124K	6,542,925	13
DoubleElectron_Run2012B-22Jan2013-v1_1	24GB	rootuple copied from Castor to "work", batch: run-4e-davec.sh	output_DoubleElectron_Run2012B-22Jan2013-v1_1.root	117K	1,647,962	4
DoubleElectron_Run2012C-22Jan2013-v1	100GB	rootuple analysed directly from Castor with batch: run-4e-davec.sh, using dirInput="root://castorcms//castor/cern.ch/user/n/ndefilip/Paper/data";	output_DoubleElectron_Run2012C-22Jan2013-v1.root	126K	6,540,600	19
DoubleElectron_Run2012C-22Jan2013-v1_1	94GB	From histos4e.tgz	output_DoubleElectron_Run2012C-22Jan2013-v1_1.root	124K	6,175,001	14
DoubleElectron_Run2012D-22Jan2013-v1	100GB	rootuple analysed directly from Castor with batch: run-4e-davec.sh, using dirInput="root://castorcms//castor/cern.ch/user/n/ndefilip/Paper/data";	output_DoubleElectron_Run2012D-22Jan2013-v1.root	124K	6,378,311	13
DoubleElectron_Run2012D-22Jan2013-v1_1	50GB	From histos4e.tgz	output_DoubleElectron_Run2012D-22Jan2013-v1_1.root	121K	3,379,685	8
					Total 4e	74

Processed data from histos4e.tgz, and from my analysis, are in:
~/HZZ/CMSSW_5_3_9/src/HiggsAnalysis/HiggsToZZ4Leptons/test/macros/FSR2012/INCLUSIVE_MACROS/log/lustre/cms/store/user/defilip/histos4e_nobnn

Castor files are:
-rw-r--r-- 1 cms003 zh 24476683419 23GB Jun 25 2013 /castor/cern.ch/user/n/ndefilip/Paper/data/roottree_leptons_DoubleElectron_Run2012A-22Jan2013-v1.root
-rw-r--r-- 1 cms003 zh 100051873173 100GB Jun 25 2013 /castor/cern.ch/user/n/ndefilip/Paper/data/roottree_leptons_DoubleElectron_Run2012B-22Jan2013-v1.root
-rw-r--r-- 1 cms003 zh 25194986104 24GB Jun 25 2013 /castor/cern.ch/user/n/ndefilip/Paper/data/roottree_leptons_DoubleElectron_Run2012B-22Jan2013-v1_1.root
-rw-r--r-- 1 cms003 zh 100047346417 100GB Jun 25 2013 /castor/cern.ch/user/n/ndefilip/Paper/data/roottree_leptons_DoubleElectron_Run2012C-22Jan2013-v1.root
-rw-r--r-- 1 cms003 zh 94947076917 94GB Jun 25 2013 /castor/cern.ch/user/n/ndefilip/Paper/data/roottree_leptons_DoubleElectron_Run2012C-22Jan2013-v1_1.root
-rw-r--r-- 1 cms003 zh 100053979480 100GB Jun 27 2013 /castor/cern.ch/user/n/ndefilip/Paper/data/roottree_leptons_DoubleElectron_Run2012D-22Jan2013-v1.root
-rw-r--r-- 1 cms003 zh 52997023923 50GB Jun 25 2013 /castor/cern.ch/user/n/ndefilip/Paper/data/roottree_leptons_DoubleElectron_Run2012D-22Jan2013-v1_1.root

The full list of roottree datasets on castor

8 Aug 2014

Running on the HZZ datasets:

• directory is ~/HZZ/CMSSW_5_3_9/src/HiggsAnalysis/HiggsToZZ4Leptons/test/macros/FSR2012/INCLUSIVE_MACROS
• Do ./HZZload.sh to get CMS environment
• edit run-4e-davec.sh or submit_HZZ4LeptonsAnalysis_CERN-davec.sh for inputs/outputs etc
• edit compilereference_4e_data-davec.C or compilereference_4e_data.C which have the main() loop
• edit HZZ4LeptonsAnalysis_4e.h and HZZ4LeptonsAnalysis_4e.C for analysis
• Compile with:
• bash compilereference-4e-davec.sh, or bash compilereference.sh, to compile things
• Makes executables:
• RunReference4e_data-davec RunReference4e
• ./run-4e-davec.sh to run job locally or ./submit_HZZ4LeptonsAnalysis_CERN-davec.sh
• job submitted from ./loopcheck_data_4e-davec.sh
• loopcheck picks up datasets to loop over from:
• data_input_4e_2012_paper.txt
• loopcheck writes a one line file for each of the datasets to /DataCards4e$2/data_input_${n}.txt where $2 = 2012
• loopcheck modifies submit_HZZ4LeptonsAnalysis_CERN-davec.sh with
• cat submit_HZZ4LeptonsAnalysis_CERN-davec.sh to make the final sh file jobs4e/submit_ReferenceAnalysis_${samplename}_4e.sh
• job submitted inside loopcheck is: bsub -q 8nh jobs4e/submit_ReferenceAnalysis_${samplename}_4e.sh

data_input_4e_2012_paper.txt
DoubleElectron_Run2012A-22Jan2013-v1
DoubleElectron_Run2012B-22Jan2013-v1
DoubleElectron_Run2012B-22Jan2013-v1_1
DoubleElectron_Run2012C-22Jan2013-v1
DoubleElectron_Run2012C-22Jan2013-v1_1
DoubleElectron_Run2012D-22Jan2013-v1
DoubleElectron_Run2012D-22Jan2013-v1_1

5 Aug 2014

Root files output from jobs at
~/HZZ/CMSSW_5_3_9/src/HiggsAnalysis/HiggsToZZ4Leptons/test/macros/FSR2012/INCLUSIVE_MACROS/log/lustre/cms/store/user/defilip/histos2e2mu_nobnn $
Files like:

-rw-r--r-- 1 davec zh  1.5M Dec  3  2013 output_VBF_ToHToZZTo4L_M-210_7TeV-powheg-pythia6.root
-rw-r--r-- 1 davec zh     0 Dec  3  2013 output_VBF_ToHToZZTo4L_M-220_7TeV-powheg-pythia6.root
-rwxr-xr-x 1 davec zh  8.3M Dec  3  2013 output_VBFToHToZZTo4L_M-225_8TeV-powheg15-pythia6.root
-rw-r--r-- 1 davec zh     0 Dec  3  2013 output_VBF_ToHToZZTo4L_M-230_7TeV-powheg-pythia6.root
-rw-r--r-- 1 davec zh  1.5M Dec  3  2013 output_VBF_ToHToZZTo4L_M-250_7TeV-powheg-pythia6.root
-rwxr-xr-x 1 davec zh     0 Dec  3  2013 output_VBFToHToZZTo4L_M-250_8TeV-powheg15-pythia6.root
-rw-r--r-- 1 davec zh  1.5M Dec  3  2013 output_VBF_ToHToZZTo4L_M-275_7TeV-powheg-pythia6.root
-rwxr-xr-x 1 davec zh  8.8M Dec  3  2013 output_VBFToHToZZTo4L_M-275_8TeV-powheg15-pythia6.root
-rw-r--r-- 1 davec zh  1.5M Dec  3  2013 output_VBF_ToHToZZTo4L_M-300_7TeV-powheg-pythia6.root
-rwxr-xr-x 1 davec zh     0 Dec  3  2013 output_VBFToHToZZTo4L_M-300_8TeV-powheg15-pythia6.root
-rw-r--r-- 1 davec zh     0 Dec  3  2013 output_VBF_ToHToZZTo4L_M-325_7TeV-powheg-pythia6.root
-rw-r--r-- 1 davec zh     0 Dec  3  2013 output_VBF_ToHToZZTo4L_M-350_7TeV-powheg-pythia6.root
-rwxr-xr-x 1 davec zh     0 Dec  3  2013 output_VBFToHToZZTo4L_M-350_8TeV-powheg15-pythia6.root
-rw-r--r-- 1 davec zh  1.5M Dec  3  2013 output_VBF_ToHToZZTo4L_M-400_7TeV-powheg-pythia6.root
-rwxr-xr-x 1 davec zh  9.0M Dec  3  2013 output_VBFToHToZZTo4L_M-400_8TeV-powheg15-pythia6.root
-rw-r--r-- 1 davec zh  1.6M Dec  3  2013 output_VBF_ToHToZZTo4L_M-425_7TeV-powheg-pythia6.root
-rw-r--r-- 1 davec zh  1.6M Dec  3  2013 output_VBF_ToHToZZTo4L_M-450_7TeV-powheg-pythia6.root
-rwxr-xr-x 1 davec zh  9.3M Dec  3  2013 output_VBFToHToZZTo4L_M-450_8TeV-powheg15-pythia6.root

30 Jul 2014

Running Step 2 of the exercise. .L Sgn_Bkg_superimpose2.C Sgn_Bkg_superimpose()
Challenges to set up pad size within c canvas. Can now print latex on canvas to give data details.

Particle flow isolation is on line 2942 of HZZ4LeptonsAnalysis_4mu.C. The summed energy in a cone of deltaR = 0.3 divided by the muon pT. The summed energy is from neutral and charged hadrons, photons but pileup (charged??) is subtracted:

Old:
//RECOMU_PFX_dB_new[0]=(RECOMU_PFchHad[iL[i]]+
      //          max( (RECOMU_PFneuHad[iL[i]]+
      //           (RECOMU_PFphoton[iL[i]]-RECOPFPHOT_PT[iLp[pi]])-
      //           max(RHO_mu,0.0)*(EffectiveArea)), 0.0) ) / 
                //                  RECOMU_PT[iL[i]];

New:
 RECOMU_PFX_dB_new[1]=(RECOMU_PFchHad[iL[j]] +
                          max (   0.,RECOMU_PFneuHad[iL[j]]+ 
                              (RECOMU_PFphoton[iL[j]]-RECOPFPHOT_PT[iLp[pj]]) - 
                              0.5*RECOMU_PFPUchAllPart[iL[j]])  )  /
                                 RECOMU_PT[iL[j]];

Observations:
Mass_Z1 low side cuts possible to remove ZZ->4mu and TT->2l+2nu+2b
hIso > 1.5 effective to remove TT->2l+2nu+2b
hSip > ~5 effective to remove TT->2l+2nu+2b
* 3datasets-3plots-step3.png:

28 Jul 2014

Going through Step 2. Compilation OK. However programme crashes. Step2-crash-dump.txt
Find code requests:
TFile *_filePUMC = TFile::Open(PU_MC); where PU_MC = PUvertices/PUvertices_%s.root = PUvertices/PUvertices_output.root.
PUvertices/PUvertices_output.root was missing. Copied one of the PUvertices_*.root files to PUvertices/PUvertices_output.root.
All OK now.

Note:
echo $ROOTSYS
/afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.19/x86_64-slc5-gcc47-opt/root/
after "cmsenv" becomes:
/cvmfs/cms.cern.ch/slc5_amd64_gcc462/lcg/root/5.32.00-cms21/

Look at the 3 GB file plots in work/HZZ/:
TFile *_file0 = TFile::Open("/afs/cern.ch/work/d/davec/HZZ/roottree_leptons_GluGluToHToZZTo4L_M-125_8TeV-powheg15-pythia6.root")
Tree: HZZ4LeptonsAnalysis
From command line, to avoid -999 initialisations do
HZZ4LeptonsAnalysis->Draw("RECOELE_ETA","RECO_ETA>-999.")
177530 events
Expect 177530 * 4 = 710120 leptons. In plot MC_LEPT_ETA have 677709 leptons, so lost 32411 leptons = 4.5% of generated events.

In output.root:
Have only 12411 events from 177530 events after cuts = 7%, from plots in output.root directory HZZ4LeptonsAnalysisReduced;1

Some plots in output.root:
* output-root-f-njets-pass.png:

Some plots in roottree_leptons_GluGluToHToZZTo4L_M-125_8TeV-powheg15-pythia6.root :

• GGHZZ-3GB-Run.png:
  

• GGHZZ-3GB-nvertices.png:
  

• GGHZZ-3GB-Events.png:
  

• GGHZZ-3GB-LEPT_ETA.png:
  

• GGHZZ-3GB-MC_eta.png:
  

• GGHZZ-3GB-MC_mass.png:
  

• GGHZZ-3GB-MC_theta.png:
  

• GGHZZ-3GB-RECOELE_ETA.png:
  

• GGHZZ-3GB-RECO_ZEE_MASS.png:
  

22 Jul 2014

First part of exercise.

Took from 10.00 to 01.00 next day to get root file via rfcp from CASTOR.
rfcp /castor/cern.ch/user/n/ndefilip/Paper/data/roottree_leptons_DoubleMu_Run2012D-22Jan2013-v1_3.root .

root -l roottree_leptons_DoubleMu_Run2012D-22Jan2013-v1_3.root

root> TBrowser browser
root> HZZ4LeptonsAnalysis->Draw("RECO_NMU","RECO_NMU>-999.") etc

How many events does it contain? 59617
How many muons per event? (check RECO_NMU): 2.06
What is their muon pT spectrum? (RECOMU_PT): Mean pT = 21.6 GeV
What is their isolation distribution? (RECOMU_PFX_dB): Mean = 1.75
Plot the dilepton mass opposite charge (RECO_ZMM_MASS[0]): 2nd peak at 90 GeV

• RECO_NMU.png:
  

• RECOMU_PT.png:
  

• RECOMU_PFX_dB.png:
  

• Reco_ZMM_Mass0.png:
  

-- DavidCockerill - 23 Jul 2014

• Step2-crash-dump.txt: Step2-crash-dump.txt

• GGHZZ-3GB-Events.png:
  

• GGHZZ-3GB-LEPT_ETA.png:
  

• GGHZZ-3GB-MC_eta.png:
  

• GGHZZ-3GB-MC_mass.png:
  

• GGHZZ-3GB-MC_theta.png:
  

• GGHZZ-3GB-RECOELE_ETA.png:
  

• GGHZZ-3GB-RECO_ZEE_MASS.png:
  

• GGHZZ-3GB-Run.png:
  

• GGHZZ-3GB-nvertices.png:
  

• output-root-f-njets-pass.png:
  

• 3datasets-3plots-step3.png:
  

• data_input_4e_2012_paper.txt: data_input_4e_2012_paper.txt

• data_input_4mu_2012_paper.txt: data_input_4mu_2012_paper.txt

• data_input_2e2mu_2012_paper.txt: data_input_2e2mu_2012_paper.txt

• roottree-data.txt: roottree-data.txt

• data_input_4e_2011_paper.txt: data_input_4e_2011_paper.txt

• data_input_4mu_2011_paper.txt: data_input_4mu_2011_paper.txt

• data_input_2e2mu_2011_paper.txt: data_input_2e2mu_2011_paper.txt

• PlotStack4l-davec-output-text.txt: PlotStack4l-davec-output-text.txt

• filelist_4mu_2012_paper_CERN-davec.txt: filelist_4mu_2012_paper_CERN-davec.txt

• h_hMZ1_3_4mu_8TeV_ratio_log.png: