Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Adding badGlobalMuonTagger to ICMuonProducer #167

Open
wants to merge 78 commits into
base: master
Choose a base branch
from
Open
Changes from 1 commit
Commits
Show all changes
78 commits
Select commit Hold shift + click to select a range
97a3a89
Changes for eta 3 spikes, coming from Raffaele. TO BE CHECKED by AM
Dec 1, 2016
3f20049
fix bug on: Changes for eta 3 spikes, coming from Raffaele. TO BE CHE…
Dec 1, 2016
079db6e
Saving trigger bits for all METMHT. TO BE CHECKED by AM
Dec 1, 2016
1617c01
Adding trigger bits for new METMHTs
Dec 1, 2016
1f998be
changes to leptons and vetos for synch with monojet
amagnan Dec 1, 2016
8fa0af1
Removing redundancy as suggested by AM
Dec 1, 2016
fa04846
fix conflicts
amagnan Dec 1, 2016
7776cb7
add trigger weights for metmht
amagnan Dec 2, 2016
7d860a3
changes to ltana
amagnan Dec 2, 2016
012258b
LTs for Maria
Dec 2, 2016
7cca73a
Merge branch 'HToInv_RDM' of github.com:riccardodimaria/ICHiggsTauTau…
Dec 2, 2016
d11cb5b
Addressing first 2 point in the TO-DO list from AM + LTs for Maria
Dec 2, 2016
9fff0eb
adding file
amagnan Dec 2, 2016
639eaf4
Merge branch 'HToInv_RDM' of https://github.com/riccardodimaria/ICHig…
amagnan Dec 2, 2016
4ba1448
Adding charged_had_energy_frac
Dec 2, 2016
dcb1d95
correcting idweights to include also loose leptons
amagnan Dec 2, 2016
b35adaa
Merge branch 'HToInv_RDM' of https://github.com/riccardodimaria/ICHig…
amagnan Dec 2, 2016
a729646
Adding 161031 filelist up-to-date
Dec 2, 2016
2cda499
adding dielectron CR
amagnan Dec 4, 2016
5a5d01f
TEST submission LTAna
Dec 4, 2016
1790ec8
adding single to filelist
amagnan Dec 4, 2016
b8f8f06
TEST submission LTAna
Dec 4, 2016
dec39f7
TEST submission LTAna
Dec 4, 2016
430449c
Adding Zee to make datacards
Dec 4, 2016
08a71de
Adding Zee to make datacards
Dec 4, 2016
be0b4af
Fixing possible bugs: TO BE CHECKED by AM
Dec 4, 2016
d68c372
separating trig weight properly
amagnan Dec 5, 2016
b9504f7
fix conflict
amagnan Dec 5, 2016
19dd2d6
Changes before HExo
Dec 5, 2016
f3ee577
Merge branch 'HToInv_RDM' of github.com:riccardodimaria/ICHiggsTauTau…
Dec 5, 2016
52bbd35
propagate change of trigger weight to e regions
amagnan Dec 5, 2016
de22f51
Merge branch 'HToInv_RDM' of https://github.com/riccardodimaria/ICHig…
amagnan Dec 5, 2016
46344e5
enu CR still not working
Dec 6, 2016
18431be
cout working for Powheg-VBFHtoinv-mH125 but not for WJetsToLNu-mg_spl…
Dec 6, 2016
5645883
Fixing issue with ee and enu regions, adding tauveto and bveto
Dec 7, 2016
8d23045
run on AM folder
amagnan Dec 8, 2016
f7f33bd
Merge branch 'HToInv_RDM' of https://github.com/riccardodimaria/ICHig…
amagnan Dec 8, 2016
fdd6fea
Fixed issues due to change variable name
Dec 8, 2016
f42b5ef
Single Electron LTs (B C D run) failed, others (E F run) are fine; [N…
Dec 9, 2016
8820f21
Fixing bugs and disabling type1cor_ in JetMETModifier.cc
Dec 12, 2016
62ae329
adding filelist for DATA
Dec 21, 2016
a9316f0
fixcards bash script
Jan 4, 2017
7ebaaac
fixcards bash script
Jan 4, 2017
f7ddc26
adding all jet frac variables to tree
amagnan Jan 6, 2017
4430f36
fixing conflicts
amagnan Jan 6, 2017
8ad8423
upadting to 2016 data electron id
amagnan Jan 9, 2017
7de1cc9
fixing bug in jet smearing for unmatched jets
amagnan Jan 11, 2017
864db5b
Trying to reweight W/Z from kfactors.root
Jan 11, 2017
903b2b8
Adding kfactors from Zeynep
Jan 11, 2017
0e4bf8a
Merge remote-tracking branch 'origin' into HToInv_RDM
Jan 11, 2017
94ab1e8
Trying to reweight W/Z from kfactors.root
Jan 11, 2017
94656e1
Trying to reweight W/Z from kfactors.root
Jan 11, 2017
7354de6
Up to Final run
Jan 13, 2017
8729e3b
New LTs after latest AM chages, 170113_ICHEP
Jan 13, 2017
8e310c9
Latest plots and limit, 170116
Jan 16, 2017
4cbc62d
fixing histos boundaries, adding variable metcut for enu region
amagnan Jan 17, 2017
ab1d589
merge conflicts fixed
amagnan Jan 17, 2017
ee54d0e
adding 3D hist for datacards
amagnan Jan 18, 2017
0444fef
adding code to optimise on 3D histo
amagnan Jan 19, 2017
4c57589
Pointless changes for compariso table
Jan 22, 2017
68612d3
Merge branch 'HToInv_RDM' of github.com:riccardodimaria/ICHiggsTauTau…
Jan 22, 2017
0924b52
Bugs after AM commit
Jan 23, 2017
ac8bafd
clean .hists files and correct log cout of integral adding overflows
amagnan Jan 24, 2017
126193a
NLO reweighting should be validated; right now it is giving silly res…
Jan 24, 2017
ec8c3de
Merge branch 'HToInv_RDM' of github.com:riccardodimaria/ICHiggsTauTau…
Jan 24, 2017
bd555a4
Comparison table ultimated
Jan 24, 2017
2c75769
Adding some private stuff in order to switch CMSSW version (CMSSW_8_0…
Jan 26, 2017
4f7188a
Merge remote-tracking branch 'origin' into HToInv_RDM
Jan 26, 2017
92ba7aa
Changing ICMuonProducer to select badMuons from Giovanni's patch
Jan 27, 2017
e0b631c
Removing redundant cuts from SR and QCD CR at LTAna step
Jan 27, 2017
5438581
Fixing bugs in producer
ahinzmann Jan 19, 2017
1173cb1
Try if everything is ok
Jan 30, 2017
03bb2a9
Good before cherry-pick
Jan 30, 2017
1ca223b
After cherry-pick: TO BE CHECKED BY AM
Jan 30, 2017
b6bb34e
Merge pull request #1 from riccardodimaria/HToInv_RDM
Jan 30, 2017
30ad750
Removing single files from filelists
Jan 30, 2017
3b317d3
Adding 161031.tar.gz
Jan 30, 2017
b164d94
Merge pull request #2 from riccardodimaria/HToInv_RDM
Jan 30, 2017
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
adding dielectron CR
amagnan committed Dec 4, 2016
commit 2cda49957ac6c3e554d141646725caa70eba664f
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
output_name=tmp.root
#input_folder=/vols/cms/rd1715/HiggsToInv/output_lighttree_161121_ICHEP/
input_folder=/vols/cms/magnan/Hinvisible/RunIILT/output_lighttree_161201/
input_folder=/vols/cms/magnan/Hinvisible/RunIILT/output_lighttree_161202/
#syst=
input_params=../filelists/161031/Params161031.dat
filelist=filelists/run2filelist80X.dat
@@ -12,6 +12,7 @@ listset=""
#uncomment for QCD mindphi plot in signal region
#basesel=jet1_eta*jet2_eta<0 && abs(jet1_eta)<4.7 && abs(jet2_eta)<4.7 && dijet_M>800&&jet1_pt>80&&dijet_deta>3.6&& jet2_pt>70&&metnomuons>0
basesel=jet1_eta*jet2_eta<0 && abs(jet1_eta)<4.7 && abs(jet2_eta)<4.7 && dijet_M>1100&&jet1_pt>80&&dijet_deta>3.6&& jet2_pt>40&&metnomuons>200 && nloosephotons==0
baseselele=jet1_eta*jet2_eta<0 && abs(jet1_eta)<4.7 && abs(jet2_eta)<4.7 && dijet_M>1100&&jet1_pt>80&&dijet_deta>3.6&& jet2_pt>40&&metnoelectrons>200
channel=munu
runblind=true
#runblind=false #to run unblind
36 changes: 36 additions & 0 deletions Analysis/HiggsNuNu/LightTreeAna/scripts/ee.hists
Original file line number Diff line number Diff line change
@@ -0,0 +1,36 @@
m_ee(12,60.,120.) ;m_{ee};Events
pt_ee(15,0.,400.) ;p_{T}^{ee};Events
ele1_pt(22,20.,360.) ;p_{T}(lead electron) (GeV);Events
ele1_eta(30,-2.5,2.5) ;#eta(lead electron) (GeV);Events
ele2_pt(17,20.,275.) ;p_{T}(sublead electron) (GeV);Events
ele2_eta(30,-2.5,2.5) ;#eta(sublead electron) (GeV);Events
metnoelectrons*pow(sqrt(sumet-ele1_pt-ele2_pt),-1)(15,0.,15.) ;S;Events
n_vertices(20,0,40) ;n_{vtx};Events
jet1_pt(28,80.,500.) ;p_{T}^{j1} (GeV);Events
jet2_pt(12,40.,250.) ;p_{T}^{j2} (GeV);Events
jet1_uncorpt(28,80.,500.) ;uncor p_{T}^{j1} (GeV);Events
jet2_uncorpt(12,40.,250.) ;uncor p_{T}^{j2} (GeV);Events
jet1_eta(25,-5,5) ;#eta^{j1} (GeV);Events
jet2_eta(25,-5,5) ;#eta^{j2} (GeV);Events
cjvjetpt(19,15.,205.) ;p_{T}^{CJV jet} (GeV);Events
met(20,0.,550.) ;MET (GeV);Events
l1met(20,0.,550.) ;L1MET (GeV);Events
ht(20,100,1000) ;H_{T}^{p_{T}^{j}>15 GeV} (GeV);Events
ht30(20,100,1000) ;H_{T}^{p_{T}^{j}>30 GeV} (GeV);Events
metnoelectrons(10,200.,450.) ;E_{T,no-#mu}^{miss} (GeV);Events")
dijet_M(10,1100.,3000.) ;M_{jj} (GeV);Events
dijet_dphi(15,0.,3.1416) ;#Delta#phi_{jj};Events
dijet_deta(10,3.6,5.5) ;#Delta#eta_{jj};Events
dijet_sumeta(20,-5,5) ;#eta_{j1}+#eta_{j2};Events
jetmetnoel_mindphi(8,2.3,3.1416) ;#Delta#phi(E_{T}^{miss},j_{1,2});Events
alljetsmetnoel_mindphi(8,2.3,3.1416) ;#Delta#phi(E_{T,no-#el}^{miss},j);Events
n_jets_cjv_30(5,0,5) ;CJV jets (30 GeV);Events
n_jets_cjv_20EB_30EE(5,0,5) ;CJV jets (20 GeV EB, 30 GeV EE);Events
jet1_csv(21,0.,1.) ;Jet 1 CSV;Events
jet2_csv(21,0.,1.) ;Jet 2 CSV;Events
jet3_csv(21,0.,1.) ;Jet 3 CSV;Events
n_jets_30(10,0.,10.) ;N jets pt>30 GeV;Events
n_jets_15(15,0.,15.) ;N jets pt>15 GeV;Events
central_tag_eta(25,-5.,5.) ;Central tag jet #eta;Events
forward_tag_eta(25,-5.,5.) ;Forward tag jet #eta;Events
alljetsmetnoel_mindphi:dijet_M(24,800.,2000.,20,1.,3.) ;M_{jj} (GeV);#Delta#phi(E_{T}^{miss},j);Events
4 changes: 4 additions & 0 deletions Analysis/HiggsNuNu/LightTreeAna/scripts/ee_sig.hists
Original file line number Diff line number Diff line change
@@ -0,0 +1,4 @@
metnoelectrons(10,200.,450.) ;E_{T,no-e}^{miss} (GeV);Events
dijet_M(12,1100.,4000.) ;M_{jj} (GeV);Events
dijet_deta(10,3.6,5.5) ;#Delta#eta_{jj};Events
alljetsmetnomu_mindphi(8,2.3,3.1416) ;#Delta#phi(E_{T,no-#mu}^{miss},j);Events
10 changes: 5 additions & 5 deletions Analysis/HiggsNuNu/LightTreeAna/scripts/enu.hists
Original file line number Diff line number Diff line change
@@ -1,16 +1,16 @@
lep_mt(18,0.,180.) ;m_{T}(lepton+MET) (GeV);Events
ele1_pt(22,20.,350.) ;p_{T}(electron) (GeV);Events
ele1_eta(20,-2.4,2.4) ;#eta(electron) (GeV);Events
metnomuons*pow(sqrt(sumet),-1)(30,0.,15.) ;S;Events
ele1_eta(20,-2.5,2.5) ;#eta(electron) (GeV);Events
metnoelectrons*pow(sqrt(sumet),-1)(30,0.,15.) ;S;Events
n_vertices(40,0,40) ;n_{vtx};Events
jet1_pt(32,80.,560.) ;p_{T}^{j1} (GeV);Events
jet2_pt(16,70.,310.) ;p_{T}^{j2} (GeV);Events
jet2_pt(16,40.,310.) ;p_{T}^{j2} (GeV);Events
jet1_uncorpt(32,80.,560.) ;uncor p_{T}^{j1} (GeV);Events
jet2_uncorpt(16,70.,310.) ;uncor p_{T}^{j2} (GeV);Events
jet2_uncorpt(16,40.,310.) ;uncor p_{T}^{j2} (GeV);Events
jet1_eta(47,-4.7,4.7) ;#eta^{j1} (GeV);Events
jet2_eta(47,-4.7,4.7) ;#eta^{j2} (GeV);Events
cjvjetpt(19,15.,205.) ;p_{T}^{CJV jet} (GeV);Events
metnomuons(25,200.,575.) ;E_{T,no-#mu}^{miss} (GeV);Events
metnoelectrons(25,200.,575.) ;E_{T,no-#mu}^{miss} (GeV);Events
met(30,0.,550.) ;MET (GeV);Events
l1met(30,0.,550.) ;L1MET (GeV);Events
ht(30,100,1000) ;H_{T}^{p_{T}^{j}>15 GeV} (GeV);Events
8 changes: 4 additions & 4 deletions Analysis/HiggsNuNu/LightTreeAna/scripts/mumu.hists
Original file line number Diff line number Diff line change
@@ -1,15 +1,15 @@
m_mumu(12,60.,120.) ;m_{#mu#mu};Events
pt_mumu(15,0.,400.) ;p_{T}^{#mu#mu};Events
mu1_pt(22,20.,360.) ;p_{T}(lead muon) (GeV);Events
mu1_eta(30,-2.1,2.1) ;#eta(lead muon) (GeV);Events
mu1_eta(30,-2.5,2.5) ;#eta(lead muon) (GeV);Events
mu2_pt(17,20.,275.) ;p_{T}(sublead muon) (GeV);Events
mu2_eta(30,-2.1,2.1) ;#eta(sublead muon) (GeV);Events
mu2_eta(30,-2.5,2.5) ;#eta(sublead muon) (GeV);Events
metnomuons*pow(sqrt(sumet-mu1_pt-mu2_pt),-1)(15,0.,15.) ;S;Events
n_vertices(20,0,40) ;n_{vtx};Events
jet1_pt(28,80.,500.) ;p_{T}^{j1} (GeV);Events
jet2_pt(12,70.,250.) ;p_{T}^{j2} (GeV);Events
jet2_pt(12,40.,250.) ;p_{T}^{j2} (GeV);Events
jet1_uncorpt(28,80.,500.) ;uncor p_{T}^{j1} (GeV);Events
jet2_uncorpt(12,70.,250.) ;uncor p_{T}^{j2} (GeV);Events
jet2_uncorpt(12,40.,250.) ;uncor p_{T}^{j2} (GeV);Events
jet1_eta(25,-5,5) ;#eta^{j1} (GeV);Events
jet2_eta(25,-5,5) ;#eta^{j2} (GeV);Events
cjvjetpt(19,15.,205.) ;p_{T}^{CJV jet} (GeV);Events
6 changes: 3 additions & 3 deletions Analysis/HiggsNuNu/LightTreeAna/scripts/munu.hists
Original file line number Diff line number Diff line change
@@ -1,12 +1,12 @@
lep_mt(18,0.,180.) ;m_{T}(lepton+MET) (GeV);Events
mu1_pt(22,20.,350.) ;p_{T}(muon) (GeV);Events
mu1_eta(20,-2.1,2.1) ;#eta(muon) (GeV);Events
mu1_eta(20,-2.5,2.5) ;#eta(muon) (GeV);Events
metnomuons*pow(sqrt(sumet-mu1_pt),-1)(30,0.,15.) ;S;Events
n_vertices(40,0,40) ;n_{vtx};Events
jet1_pt(32,80.,560.) ;p_{T}^{j1} (GeV);Events
jet2_pt(16,70.,310.) ;p_{T}^{j2} (GeV);Events
jet2_pt(16,40.,310.) ;p_{T}^{j2} (GeV);Events
jet1_uncorpt(32,80.,560.) ;uncor p_{T}^{j1} (GeV);Events
jet2_uncorpt(16,70.,310.) ;uncor p_{T}^{j2} (GeV);Events
jet2_uncorpt(16,40.,310.) ;uncor p_{T}^{j2} (GeV);Events
jet1_eta(47,-4.7,4.7) ;#eta^{j1} (GeV);Events
jet2_eta(47,-4.7,4.7) ;#eta^{j2} (GeV);Events
cjvjetpt(19,15.,205.) ;p_{T}^{CJV jet} (GeV);Events
4 changes: 2 additions & 2 deletions Analysis/HiggsNuNu/LightTreeAna/scripts/nunu.hists
Original file line number Diff line number Diff line change
@@ -1,9 +1,9 @@
metnomuons*pow(sqrt(sumet),-1)(30,0.,15.) ;S;Events
n_vertices(40,0,40) ;n_{vtx};Events
jet1_pt(47,80.,550.) ;p_{T}^{j1} (GeV);Events
jet2_pt(28,70.,350.) ;p_{T}^{j2} (GeV);Events
jet2_pt(28,40.,350.) ;p_{T}^{j2} (GeV);Events
jet1_uncorpt(47,80.,550.) ;uncor p_{T}^{j1} (GeV);Events
jet2_uncorpt(28,70.,350.) ;uncor p_{T}^{j2} (GeV);Events
jet2_uncorpt(28,40.,350.) ;uncor p_{T}^{j2} (GeV);Events
jet1_eta(94,-4.7,4.7) ;#eta^{j1} (GeV);Events
jet2_eta(94,-4.7,4.7) ;#eta^{j2} (GeV);Events
cjvjetpt(19,15.,205.) ;p_{T}^{CJV jet} (GeV);Events
Original file line number Diff line number Diff line change
@@ -71,7 +71,7 @@ for syst in "" #JESUP JESDOWN JERBETTER JERWORSE ELEEFFUP ELEEFFDOWN MUEFFUP MUE
do
mkdir -p $JOBDIR$syst
mkdir -p $OUTPUTDIR$syst
for channels in qcd enu munu taunu mumu nunu #qcd #topl topb top gamma
for channels in qcd enu munu taunu mumu ee nunu #qcd #topl topb top gamma
do
JOB=$channels
#executable expect strings separated by "!"
@@ -103,6 +103,12 @@ do
MINDPHICUT="alljetsmetnomu_mindphi\<0.5"
#MINDPHICUT="alljetsmetnomu_mindphi\>=0"
fi
if [ "$channels" = "ee" ]; then
MINDPHICUT="alljetsmetnoel_mindphi\>2.3"
fi
if [ "$channels" = "enu" ]; then
MINDPHICUT="alljetsmetnoel_mindphi\>2.3"
fi
if [ "$syst" = "" ]
then
$JOBWRAPPER "./bin/LTAnalysisRun2_2016 --cfg=$CONFIG --channel=$channels --histTitlePar='$HISTSTRING' --shapePar='$SHAPESTRING' -o $OUTPUTDIR$syst/$OUTPUTNAME --jetmetdphicut=$MINDPHICUT &> $JOBDIR$syst/$JOB.log" $JOBDIR$syst/$JOB.sh $GRIDSETUP
4 changes: 2 additions & 2 deletions Analysis/HiggsNuNu/LightTreeAna/scripts/taunu.hists
Original file line number Diff line number Diff line change
@@ -4,9 +4,9 @@ tau1_eta(15,-2.3,2.3) ;#eta(tau) (GeV);Events
metnomuons*pow(sqrt(sumet),-1)(10,0.,15.) ;S;Events
n_vertices(20,0,40) ;n_{vtx};Events
jet1_pt(16,80.,400.) ;p_{T}^{j1} (GeV);Events
jet2_pt(9,70.,250.) ;p_{T}^{j2} (GeV);Events
jet2_pt(9,40.,250.) ;p_{T}^{j2} (GeV);Events
jet1_uncorpt(16,80.,400.) ;uncor p_{T}^{j1} (GeV);Events
jet2_uncorpt(9,70.,250.) ;uncor p_{T}^{j2} (GeV);Events
jet2_uncorpt(9,40.,250.) ;uncor p_{T}^{j2} (GeV);Events
jet1_eta(20,-4.7,4.7) ;#eta^{j1} (GeV);Events
jet2_eta(20,-4.7,4.7) ;#eta^{j2} (GeV);Events
cjvjetpt(19,15.,205.) ;p_{T}^{CJV jet} (GeV);Events
89 changes: 84 additions & 5 deletions Analysis/HiggsNuNu/LightTreeAna/test/LTAnalysisRun2_2016.cpp
Original file line number Diff line number Diff line change
@@ -107,6 +107,7 @@ int main(int argc, char* argv[]){
std::string inputparams;
std::string filelist;
std::string basesel;
std::string baseselele;

std::string channel;
std::string syst;
@@ -162,6 +163,7 @@ int main(int argc, char* argv[]){
("do_list", po::value<bool>(&do_list)->default_value(false))
("listset", po::value<std::string>(&listset)->default_value("Top"))
("basesel", po::value<std::string>(&basesel)->default_value("jet1_eta*jet2_eta<0 && jet1_eta<4.7 && jet2_eta<4.7 && dijet_M>=1100&&jet1_pt>80&&dijet_deta>3.6&& jet2_pt>70&&metnomuons>200"))
("baseselele", po::value<std::string>(&baseselele)->default_value("jet1_eta*jet2_eta<0 && jet1_eta<4.7 && jet2_eta<4.7 && dijet_M>=1100&&jet1_pt>80&&dijet_deta>3.6&& jet2_pt>70&&metnoelectrons>200"))
("channel", po::value<std::string>(&channel)->default_value("nunu"))

("runblind", po::value<bool>(&runblind)->default_value(true))
@@ -228,14 +230,16 @@ int main(int argc, char* argv[]){
analysis->SetInputParams(inputparams);

std::cout<<"Base selection: "<<basesel<<std::endl;
std::cout<<"Base selection for electrons: "<<baseselele<<std::endl;

//add metsig cut
std::string metsigsel;
if(channel=="mumu") metsigsel="&& metnomuons/sqrt(sumet-mu1_pt-mu2_pt)>"+metsigcut;
else if (channel=="munu") metsigsel="&& metnomuons/sqrt(sumet-mu1_pt)>"+metsigcut;
else metsigsel="&& metnomuons/sqrt(sumet)>"+metsigcut;

basesel = basesel;//+metsigsel;
if (channel=="ee" || channel=="enu") basesel = baseselele;//+metsigsel;
else basesel = basesel;//+metsigsel;

analysis->set_baseselection(basesel);

@@ -375,7 +379,10 @@ int main(int argc, char* argv[]){
std::string dataextrasel;
std::string mcextrasel;

if (channel!="gamma"){
if (channel=="ee" || channel=="enu"){
dataextrasel="&&(pass_singleEltrigger==1)";
}
else if (channel!="gamma"){
//if(!do_mettrig) dataextrasel="&&(pass_sigtrigger==1)";
//for metmht trigger
if(!do_mettrig) dataextrasel="&&(pass_metmht90trigger==1 || pass_metmht100trigger==1 || pass_metmht110trigger==1 || pass_metmht120trigger==1)";
@@ -396,9 +403,10 @@ int main(int argc, char* argv[]){
//std::string nunuqcdcat="nvetomuons==0&&nvetoelectrons==0&&alljetsmetnomu_mindphi>1";
//std::string nunuqcdcat="nvetomuons==0&&nvetoelectrons==0&&"+jetmetdphicut;

std::string eecat="nselelectrons>=1&&nvetoelectrons==2&&nvetomuons==0&&m_ee>60&&m_ee<120&&oppsign_ee&&ele1_pt>40&&"+jetmetdphicut;
std::string mumucat="nselmuons>=1&&nvetomuons==2&&nvetoelectrons==0&&m_mumu>60&&m_mumu<120&&oppsign_mumu&&"+jetmetdphicut;
std::string munucat="nselmuons==1&&nvetomuons==1&&nvetoelectrons==0&&lep_mt>=0&&"+jetmetdphicut;
std::string enucat="nselelectrons==1&&nvetomuons==0&&nvetoelectrons==1&&"+jetmetdphicut;
std::string enucat="nselelectrons==1&&nvetomuons==0&&nvetoelectrons==1&&ele1_pt>40&&"+jetmetdphicut;
std::string taunucat="ntaus==1&&nvetomuons==0&&nvetoelectrons==0&&"+jetmetdphicut;
std::string gammacat="ntightphotons==1&&nvetomuons==0&&nvetoelectrons==0&&"+jetmetdphicut;
std::string toplcat="nvetomuons==1&&nvetoelectrons==1&&nselmuons==1&&nselelectrons==1";
@@ -411,6 +419,9 @@ int main(int argc, char* argv[]){
if(channel=="mumu"){//zmumu
sigcat=mumucat;
}
else if(channel=="ee"){//zee
sigcat=eecat;
}
else if(channel=="munu"){//wmu
sigcat=munucat;
}
@@ -452,9 +463,13 @@ int main(int argc, char* argv[]){
if (syst=="TRIG2DOWN") mcweightpufactor<<"*weight_trig_6/weight_trig_0";

if(channel=="taunu"||channel=="gamma"||channel=="nunu"||channel=="qcd") sigmcweight="total_weight_lepveto"+mcweightpufactor.str();//+mcweightpufactordebug;
//remove trigger weight for e channels which do not use signal trigger
else if (channel=="ee" || channel=="enu") sigmcweight="total_weight_leptight/weight_trig_0"+mcweightpufactor.str();//+mcweightpufactordebug;
else sigmcweight="total_weight_leptight"+mcweightpufactor.str();//+mcweightpufactordebug;
sig125mcweight="total_weight_lepveto"+mcweightpufactor.str();

if (channel=="ee" || channel == "enu") dataset="SingleElectron";

std::string bothcentral="TMath::Abs(jet1_eta)<3&&TMath::Abs(jet2_eta)<3";
std::string bothforward="TMath::Abs(jet1_eta)>=3&&TMath::Abs(jet2_eta)>=3";
std::string j2forwardj1central="TMath::Abs(jet1_eta)<3&&TMath::Abs(jet2_eta)>=3";
@@ -693,6 +708,16 @@ int main(int argc, char* argv[]){

if (use_nlo) zmumuraw.set_dataset("ZJets_ll_nlo");

DataShape zeeraw("zeeraw");
zeeraw.set_dataset("ZJets_ll")
.set_dirname("zee")
.set_shape(shape)
.set_dataweight(sigmcweight)
.set_basesel(analysis->baseselection())
.set_cat(sigcat+mcextrasel);

if (use_nlo) zeeraw.set_dataset("ZJets_ll_nlo");

DataShape qcdzmumuraw("qcdzmumuraw");
qcdzmumuraw.set_dataset("ZJets_ll")
.set_dirname("zmumuqcd")
@@ -711,6 +736,24 @@ int main(int argc, char* argv[]){
.set_basesel(analysis->baseselection())
.set_cat(sigcat+mcextrasel);

DataShape qcdzeeraw("qcdzeeraw");
qcdzeeraw.set_dataset("ZJets_ll")
.set_dirname("zeeqcd")
.set_shape(shape)
.set_dataweight(sigmcweight)
.set_basesel(analysis->baseselection())
.set_cat(sigcat+mcextrasel);

if (use_nlo) qcdzeeraw.set_dataset("ZJets_ll_nlo");

DataShape ewkzeeraw("ewkzeeraw");
ewkzeeraw.set_dataset("EWK_ZJets_ll")
.set_dirname("zeeewk")
.set_shape(shape)
.set_dataweight(sigmcweight)
.set_basesel(analysis->baseselection())
.set_cat(sigcat+mcextrasel);

DataShape wmunuraw("wmunuraw");
wmunuraw.set_dataset("WJets_munu")
.set_dirname("wmu")
@@ -1035,6 +1078,34 @@ int main(int argc, char* argv[]){
.set_sample("zmumuewk");
if(!do_mcbkg)ewkzmumuele.set_has_dderrors(1);

LTPlotElement zeeele;
zeeele.set_is_data(false)
.set_scale(1)
.set_color(kAzure + 2)
.set_in_stack(true)
.set_is_inratioden(true)
.set_legname("Z#rightarrow#mu#mu")
.set_sample("zee");
if(!do_mcbkg)zeeele.set_has_dderrors(1);
LTPlotElement qcdzeeele;
qcdzeeele.set_is_data(false)
.set_scale(getPostFitSF(channel,"zqcd"))
.set_color(kAzure + 2)
.set_in_stack(true)
.set_is_inratioden(true)
.set_legname("QCD Z#rightarrow#mu#mu")
.set_sample("zeeqcd");
if(!do_mcbkg)qcdzeeele.set_has_dderrors(1);
LTPlotElement ewkzeeele;
ewkzeeele.set_is_data(false)
.set_scale(getPostFitSF(channel,"zewk"))
.set_color(kAzure + 4)
.set_in_stack(true)
.set_is_inratioden(true)
.set_legname("EWK Z#rightarrow#mu#mu")
.set_sample("zeeewk");
if(!do_mcbkg)ewkzeeele.set_has_dderrors(1);

LTPlotElement znunuele;
znunuele.set_is_data(false)
.set_scale(1)
@@ -1129,7 +1200,7 @@ int main(int argc, char* argv[]){
elementvec.push_back(qcdwele);
elementvec.push_back(ewkwele);
}
else if(channel!="mumu"){
else if(channel!="mumu" && channel != "ee"){
if (channel!="enu" && channel!="taunu"){
elementvec.push_back(qcdwmunuele);
elementvec.push_back(ewkwmunuele);
@@ -1147,11 +1218,15 @@ int main(int argc, char* argv[]){
elementvec.push_back(qcdzmumuele);
elementvec.push_back(ewkzmumuele);
}
else if(channel=="ee") {
elementvec.push_back(qcdzeeele);
elementvec.push_back(ewkzeeele);
}
//}
elementvec.push_back(topele);
elementvec.push_back(qcdele);
elementvec.push_back(vvele);
if(channel!="mumu"&&channel!="enu"&&channel!="munu"){
if(channel!="mumu"&&channel!="ee"&&channel!="enu"&&channel!="munu"){
if (channel!="taunu"){
elementvec.push_back(qcdznunuele);
elementvec.push_back(ewkznunuele);
@@ -1193,6 +1268,8 @@ int main(int argc, char* argv[]){
dirvec.push_back("wtauewk");
dirvec.push_back("zmumuqcd");
dirvec.push_back("zmumuewk");
dirvec.push_back("zeeqcd");
dirvec.push_back("zeeewk");
dirvec.push_back("wqcd");
dirvec.push_back("wewk");
dirvec.push_back("zvvqcd");
@@ -1234,6 +1311,8 @@ int main(int argc, char* argv[]){
analysis->AddModule(&gjetsraw);
analysis->AddModule(&qcdzmumuraw);
analysis->AddModule(&ewkzmumuraw);
analysis->AddModule(&qcdzeeraw);
analysis->AddModule(&ewkzeeraw);
analysis->AddModule(&qcdznunuraw);
analysis->AddModule(&ewkznunuraw);
}
Binary file not shown.
4 changes: 3 additions & 1 deletion Analysis/HiggsNuNu/interface/JetMETModifier.h
Original file line number Diff line number Diff line change
@@ -31,7 +31,8 @@ namespace ic {
enum jetmetCor {
jecData=1,
jecMC=2,
smearMC=3
smearMC=3,
type1cor=4
};
enum jetmetSyst {
none=0,
@@ -107,6 +108,7 @@ namespace ic {
bool reapplyJecData_;
bool reapplyJecMC_;
bool smear_;
bool type1cor_;

unsigned nRho_;

Loading