Skip to content

Commit

Permalink
Merge pull request cms-analysis#1 from gkrintir/topic_thq_v2
Browse files Browse the repository at this point in the history 
small update
  • Loading branch information
@hbakhshi
hbakhshi authored Feb 8, 2017
2 parents b04cd86 + a83a8f3 commit ff830cf
Show file tree
Hide file tree
Showing 9 changed files with 719 additions and 24 deletions.
34 changes: 32 additions & 2 deletions DataFormats/interface/THQLeptonicTag.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -28,7 +28,38 @@ namespace flashgg {
const std::vector<edm::Ptr<Jet> > bJets() const { return BJets_;}
const std::vector<edm::Ptr<Jet> > nonbJets() const { return nonBJets_;}
float thqleptonicMvaRes() const {return thqleptonicMvaRes_;}


float getElecAlpha(int eleIndex) const{
float eleta = electrons()[eleIndex]->eta();

//for isolation recalculation
float Aeff = 0;
//cmssw/RecoEgamma/ElectronIdentification/data/Summer16/effAreaElectrons_cone03_pfNeuHadronsAndPhotons_80X.txt
if( fabs(eleta) <= 1.0000 ){
Aeff = 0.1703;
} else if( fabs(eleta) <= 1.4790 ){
Aeff = 0.1715;
} else if( fabs(eleta) <= 2.0000 ){
Aeff = 0.1213;
} else if( fabs(eleta) <= 2.2000 ){
Aeff = 0.1230;
} else if( fabs(eleta) <= 2.3000 ){
Aeff = 0.1635;
} else if( fabs(eleta) <= 2.4000 ){
Aeff = 0.1937;
} else if( fabs(eleta) <= 5.0000 ){
Aeff = 0.2393;
}

return Aeff;

};
/*
float getMuoDz(int muIndex) const{
mouons()[muIndex]->muonBestTrack()->dz( diPhoton()->vtx()->position() ) ;
}
*/
// const reco::LeafCandidate getLepton() const{
// if( electrons().size() == 1 && muons().size() == 0 ){
// return ( *(electrons()[0]) );
Expand Down Expand Up @@ -119,12 +150,11 @@ namespace flashgg {
void setbJet( Ptr<Jet> bj ) { bJet.clear() ; bJet.push_back( bj ) ;}


void setValues( float fox , float aplan , float topMass , float met, float metPhi ){
void setValues( float fox , float aplan , float topMass , float met, float metPhi ){
FoxWolframMoment_ONE = fox ;
Aplanarity = aplan;

TopMass = topMass;

MET = met ;
MET_Phi = metPhi ;
}
Expand Down
1 change: 1 addition & 0 deletions MetaData/data/cross_sections.json
View file
Original file line number Diff line number Diff line change
Expand Up @@ -70,6 +70,7 @@
"TTGG_0Jets_TuneCUETP8M1_13TeV_amcatnlo_madspin_pythia8" : { "xs" : 0.017, "br" : 1.0, "kf" : 1.0, "itype":1 },
"TTGJets_TuneCUETP8M1_13TeV-amcatnloFXFX-madspin-pythia8" : { "xs" : 3.697, "br" : 1.0, "kf" : 1.0, "itype":1 },
"TGJets_TuneCUETP8M1_13TeV_amcatnlo_madspin_pythia8" : { "xs" : 2.967, "br" : 1.0, "kf" : 1.0, "itype":1 },
"ZGTo2LG_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8" : { "xs" : 117.864, "br" : 1.0, "kf" : 1.0, "itype":1 },
"VHToGG_M120_13TeV_amcatnloFXFX_madspin_pythia8" : { "xs" : 2.559, "br" : 0.00222, "itype":-120400 },
"VHToGG_M121_13TeV_amcatnloFXFX_madspin_pythia8" : { "xs" : 2.495, "br" : 0.00223, "itype":-121400 },
"VHToGG_M122_13TeV_amcatnloFXFX_madspin_pythia8" : { "xs" : 2.433, "br" : 0.00225, "itype":-122400 },
Expand Down
20 changes: 2 additions & 18 deletions Systematics/python/SystematicDumperDefaultVariables.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -18,24 +18,8 @@
"leadIDMVA :=leadingView.phoIdMvaWrtChosenVtx",
"subleadIDMVA :=subLeadingView.phoIdMvaWrtChosenVtx",
"maxEta :=max(abs(diPhoton().leadingPhoton.superCluster.eta),abs(diPhoton().leadingPhoton.superCluster.eta))",
"vtxZ :=diPhoton().vtx().z",
"bjet1_pt := ?bJets.size>0? bJets.at(0).pt : -1",
"bjet2_pt := ?bJets.size>1? bJets.at(1).pt : -1",
"bjet1_eta:= ?bJets.size>0? bJets.at(0).eta: -1",
"bjet2_eta:= ?bJets.size>1? bJets.at(1).eta:-1",
"bjet1_abseta := ?bJets.size>0? abs(bJets.at(0).eta): -1",
"bjet2_abseta := ?bJets.size>1? abs(bJets.at(1).eta): -1",
"jet1_pt := ?(jets.size>0)? jets.at(0).pt : -1",
"jet2_pt := ?(jets.size>1)? jets.at(1).pt : -1",
"jet1_eta := ?(jets.size>0)? jets.at(0).eta : -1",
"jet2_eta := ?(jets.size>1)? jets.at(1).eta: -1",
"jet1_abseta := ?(jets.size>0)? abs(jets.at(0).eta): -1",
"jet2_abseta := ?(jets.size>1)? abs(jets.at(1).eta): -1",
"electronjet1deta := ?(electrons.size>0)? abs(electrons.at(0).eta - jets.at(0).eta) : -1",
"electronjet2deta := ?(electrons.size>0)? abs(electrons.at(0).eta - jets.at(1).eta) : -1",
"muonjet1deta := ?(muons.size>0)? abs(muons.at(0).eta - jets.at(0).eta) : -1",
"muonjet2deta := ?(muons.size>0)? abs(muons.at(0).eta - jets.at(1).eta) : -1",
]
"vtxZ :=diPhoton().vtx().z"
]


defaultHistograms=["CMS_hgg_mass>>mass(160,100,180)",
Expand Down
2 changes: 1 addition & 1 deletion Systematics/test/bkg_jobs_forthq.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -5,4 +5,4 @@ outdir="/afs/cern.ch/user/h/hbakhshi/work/tHq/CMSSW_8_0_25/src/flashgg/Systemati
queue="1nd"
useAAA=1
version="1"
fggRunJobs.py --load bkg_jobs_forthq.json -d $outdir/bkg_jobs_$version -x cmsRun workspaceStd_testThq.py maxEvents=-1 -n 500 -q 8nh -H -D -P useAAA=1 --no-use-tarball puTarget=2.51e+05,1.15e+06,2.47e+06,3.72e+06,5.19e+06,6.79e+06,8.67e+06,2.31e+07,5.89e+07,1.38e+08,3.12e+08,5.71e+08,8.76e+08,1.21e+09,1.56e+09,1.87e+09,2.08e+09,2.19e+09,2.24e+09,2.28e+09,2.29e+09,2.24e+09,2.15e+09,2.03e+09,1.88e+09,1.71e+09,1.54e+09,1.36e+09,1.19e+09,1.01e+09,8.48e+08,6.94e+08,5.57e+08,4.38e+08,3.37e+08,2.53e+08,1.85e+08,1.31e+08,8.96e+07,5.87e+07,3.68e+07,2.2e+07,1.25e+07,6.75e+06,3.46e+06,1.68e+06,7.79e+05,3.44e+05,1.46e+05,6.13e+04,2.68e+04,1.33e+04,8.25e+03,6.3e+03,5.45e+03,4.97e+03,4.59e+03,4.25e+03,3.92e+03,3.58e+03,3.25e+03,2.93e+03,2.62e+03,2.33e+03,2.05e+03,1.79e+03,1.56e+03,1.34e+03,1.14e+03,969,815,680,563,463,378
fggRunJobs.py --load bkg_jobs_forthq_ZG.json -d $outdir/bkg_jobs_$version -x cmsRun workspaceStd_testThq.py maxEvents=-1 -n 500 -q 1nd -D -P useAAA=1 --no-use-tarball puTarget=2.51e+05,1.15e+06,2.47e+06,3.72e+06,5.19e+06,6.79e+06,8.67e+06,2.31e+07,5.89e+07,1.38e+08,3.12e+08,5.71e+08,8.76e+08,1.21e+09,1.56e+09,1.87e+09,2.08e+09,2.19e+09,2.24e+09,2.28e+09,2.29e+09,2.24e+09,2.15e+09,2.03e+09,1.88e+09,1.71e+09,1.54e+09,1.36e+09,1.19e+09,1.01e+09,8.48e+08,6.94e+08,5.57e+08,4.38e+08,3.37e+08,2.53e+08,1.85e+08,1.31e+08,8.96e+07,5.87e+07,3.68e+07,2.2e+07,1.25e+07,6.75e+06,3.46e+06,1.68e+06,7.79e+05,3.44e+05,1.46e+05,6.13e+04,2.68e+04,1.33e+04,8.25e+03,6.3e+03,5.45e+03,4.97e+03,4.59e+03,4.25e+03,3.92e+03,3.58e+03,3.25e+03,2.93e+03,2.62e+03,2.33e+03,2.05e+03,1.79e+03,1.56e+03,1.34e+03,1.14e+03,969,815,680,563,463,378
6 changes: 6 additions & 0 deletions Systematics/test/bkg_jobs_forthq_ZG.json
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,6 @@
{
"processes" : {
"DY" : [ "/ZGTo2LG_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8"]
},
"cmdLine" : "campaign=RunIISummer16-2_4_1-25ns_Moriond17"
}
13 changes: 10 additions & 3 deletions Systematics/test/workspaceStd_testThq.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -147,8 +147,10 @@
# Or use the official tool instead
useEGMTools(process)


import flashgg.Taggers.THQLeptonicTagVariables as var
'''
variablesToUse = defaultVariables
variablesToUse.append("topMass :=getTopMass()")
variablesToUse.append("fwdJetEta :=getFwdJet().eta")
Expand All @@ -171,6 +173,10 @@
variablesToUse.append("prompt_pho_1 := diPhoton.leadingPhoton.genMatchType()")
variablesToUse.append("prompt_pho_2 := diPhoton.subLeadingPhoton.genMatchType()")
variablesToUse.append("MuonsPt := getMuonsPt()" )
'''
variablesToUse = defaultVariables + var.vtx_variables + var.dipho_variables + var.photon_variables + var.lepton_variables + var.jet_variables
print "-------------------------------------------------"
print "--- Variables to be dumped, including systematic weights ---"
print variablesToUse
Expand All @@ -191,12 +197,13 @@

process.extraDumpers = cms.Sequence()
from flashgg.Taggers.tagsDumpers_cfi import *
#process.load("flashgg.Taggers.globalVariables_cff")

process.thqLeptonicTagDumper = createTagDumper("THQLeptonicTag")
process.thqLeptonicTagDumper.dumpTrees = customize.dumpTrees
process.thqLeptonicTagDumper.dumpWorkspace = customize.dumpWorkspace
process.thqLeptonicTagDumper.nameTemplate ="$PROCESS_$SQRTS_$LABEL_$SUBCAT"

#process.thqLeptonicTagDumper.globalVariables = process.globalVariables

import flashgg.Taggers.dumperConfigTools as cfgTools
cfgTools.addCategories(process.thqLeptonicTagDumper,
Expand Down Expand Up @@ -316,7 +323,7 @@
#print >> processDumpFile, process.dumpPython()

# set default options if needed
customize.setDefault("maxEvents",-1)
customize.setDefault("maxEvents",2000)
customize.setDefault("targetLumi",1.00e+3)
# call the customization
customize(process)
Loading

0 comments on commit ff830cf

Please sign in to comment.