Skip to content

Commit

Permalink
Merge pull request cms-sw#56 from amarini/tt
Browse files Browse the repository at this point in the history 
moving TT sample
  • Loading branch information
@amarini
amarini committed May 20, 2016
2 parents 437cbdf + a77e4c1 commit b4851c6
Show file tree
Hide file tree
Showing 5 changed files with 17 additions and 59 deletions.
3 changes: 2 additions & 1 deletion dat/config.dat
View file
Original file line number Diff line number Diff line change
Expand Up @@ -12,12 +12,13 @@
#
sub=tag|v1.4
sub=amarini|/store/user/amarini/Nero
sub=klute|/store/user/klute/Nero

#### BKG ###
Files=%(amarini)s/%(tag)s/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8
addfiles=%(amarini)s/%(tag)s/WJetsToLNu_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8
#addfiles=%(amarini)s/%(tag)s/TT_TuneCUETP8M1_13TeV-amcatnlo-pythia8
addfiles=%(amarini)s/%(tag)s/TT_TuneCUETP8M1_13TeV-powheg-pythia8
addfiles=%(klute)s/%(tag)s/TT_TuneCUETP8M1_13TeV-powheg-pythia8

##powehg diboson
addfiles=%(amarini)s/%(tag)s/WWTo2L2Nu_13TeV-powheg
Expand Down
46 changes: 0 additions & 46 deletions dat/configQCDPurity.dat
View file
Original file line number Diff line number Diff line change
Expand Up @@ -3,52 +3,6 @@
##### if you want to make a config file that overwrite some changes just create a new one with the line

include=dat/config.dat
##### Files=....,....,...
### Files=%(amarini)s/%(tag)s/Tau
### addfiles=%(amarini)s/%(tag)s/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8
### addfiles=%(amarini)s/%(tag)s/WJetsToLNu_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8
### #addfiles=/store/user/amarini/Nero/%(tag)s/TT_TuneCUETP8M1_13TeV-amcatnlo-pythia8
### addfiles=%(amarini)s/%(tag)s/TT_TuneCUETP8M1_13TeV-powheg-pythia8
### ### WW/WZ/ZZ PYTHIA ###
### ##addfiles=/store/user/amarini/Nero/%(tag)s/WW_TuneCUETP8M1_13TeV-pythia8
### ##addfiles=/store/user/amarini/Nero/%(tag)s/WZ_TuneCUETP8M1_13TeV-pythia8
### ##addfiles=/store/user/amarini/Nero/%(tag)s/ZZ_TuneCUETP8M1_13TeV-pythia8
### ### POWHEG ###
### addfiles=%(amarini)s/%(tag)s/WWTo2L2Nu_13TeV-powheg
### addfiles=%(amarini)s/%(tag)s/WWToLNuQQ_13TeV-powheg
### addfiles=%(amarini)s/%(tag)s/WZTo1L1Nu2Q_13TeV_amcatnloFXFX_madspin_pythia8
### addfiles=%(amarini)s/%(tag)s/WZTo1L3Nu_13TeV_amcatnloFXFX_madspin_pythia8
### addfiles=%(amarini)s/%(tag)s/WZTo2L2Q_13TeV_amcatnloFXFX_madspin_pythia8
### addfiles=%(amarini)s/%(tag)s/WZTo3LNu_TuneCUETP8M1_13TeV-powheg-pythia8
### addfiles=%(amarini)s/%(tag)s/ZZTo2L2Nu_13TeV_powheg_pythia8
### addfiles=%(amarini)s/%(tag)s/ZZTo4L_13TeV-amcatnloFXFX-pythia8
### ### SIGNAL ###
### #addfiles=/store/user/jsalfeld/ChargedHiggs/NeroTuples/Hplus_Wh_18_5_2015
### addfiles=%(amarini)s/%(tag)s/ChargedHiggs_HplusTB_HplusToTauNu_M-200_13TeV_amcatnlo_pythia8
### addfiles=%(amarini)s/%(tag)s/ChargedHiggs_HplusTB_HplusToTauNu_M-250_13TeV_amcatnlo_pythia8
### addfiles=%(amarini)s/%(tag)s/ChargedHiggs_HplusTB_HplusToTauNu_M-300_13TeV_amcatnlo_pythia8
### addfiles=%(amarini)s/%(tag)s/ChargedHiggs_HplusTB_HplusToTauNu_M-350_13TeV_amcatnlo_pythia8
### addfiles=%(amarini)s/%(tag)s/ChargedHiggs_HplusTB_HplusToTauNu_M-400_13TeV_amcatnlo_pythia8
### addfiles=%(amarini)s/%(tag)s/ChargedHiggs_HplusTB_HplusToTauNu_M-500_13TeV_amcatnlo_pythia8
### #Signal Tau Nu ; Private production
### ## signal amc@nlo
### ### QCD
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_15to30_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_30to50_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_50to80_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_80to120_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_120to170_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_170to300_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_300to470_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_470to600_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_600to800_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_800to1000_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_1000to1400_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_1400to1800_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_1800to2400_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_2400to3200_TuneCUETP8M1_13TeV_pythia8
### addfiles=%(amarini)s/%(tag)s/QCD_Pt_3200toInf_TuneCUETP8M1_13TeV_pythia8
#

####

Expand Down
3 changes: 2 additions & 1 deletion dat/mc_database.txt
View file
Original file line number Diff line number Diff line change
Expand Up @@ -20,7 +20,8 @@ QCD-50-80 /store/user/amarini/Nero/v1.4/QCD_Pt_50to80_TuneCUETP8M1_13TeV_pythia8
QCD-600-800 /store/user/amarini/Nero/v1.4/QCD_Pt_600to800_TuneCUETP8M1_13TeV_pythia8/QCD-600-800/160513_175708/0000 3534700.0 183.354119673
QCD-800-1000 /store/user/amarini/Nero/v1.4/QCD_Pt_800to1000_TuneCUETP8M1_13TeV_pythia8/QCD-800-1000/160513_175732/0000 2784248.23049 33.2672244814
QCD-80-120 /store/user/amarini/Nero/v1.4/QCD_Pt_80to120_TuneCUETP8M1_13TeV_pythia8/QCD-80-120/160513_175519/0000 6131452.68331 2700760.35761
TT-powheg-ext3 /store/user/amarini/Nero/v1.4/TT_TuneCUETP8M1_13TeV-powheg-pythia8/TT-powheg-ext3/160513_180234/0000 94495075.0 831
#TT-powheg-ext3 /store/user/amarini/Nero/v1.4/TT_TuneCUETP8M1_13TeV-powheg-pythia8/TT-powheg-ext3/160513_180234/0000 94495075.0 831
TT-powheg-ext3 /store/user/klute/Nero/v1.4/TT_TuneCUETP8M1_13TeV-powheg-pythia8/TT-powheg-ext3/160513_180234/0000 94495075.0 831
W3Jets-madgraph /store/user/amarini/Nero/v1.4/W3JetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/W3Jets-madgraph/160513_180322/0000 7205063.0 939.461425781
#WJets-amcatnlo /store/user/amarini/Nero/v1.4/WJetsToLNu_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/WJets-amcatnlo/160513_180301/0000 3.73197735021e+12 40472.4644846
## NNLO
Expand Down
1 change: 1 addition & 0 deletions interface/PurityFitAnalytic.hpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -34,6 +34,7 @@ using namespace std;
#include "RooExponential.h"
#include "RooGaussian.h"
#include "RooGenericPdf.h"
#include "RooFFTConvPdf.h"

#include "interface/BaseFitter.hpp"

Expand Down
23 changes: 12 additions & 11 deletions src/PurityFitAnalytic.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -334,28 +334,29 @@ float PurityFitAnalytic::fit_specific( const TH1* h_, const TH1* sig_, const TH1

// ----------------------- FIT INV ISO TO EXTRACT QCD MODEL ---------------------------
if (initvalues.find("sigmaR_QCD") == initvalues.end())initvalues["sigmaR_QCD"] = 20;
if (initvalues.find("cE_QCD") == initvalues.end())initvalues["cE_QCD"] = 0.001;
if (initvalues.find("sigmaG_QCD") == initvalues.end())initvalues["sigmaG_QCD"] = 20;
if (initvalues.find("muG_QCD") == initvalues.end())initvalues["muG_QCD"] = 5;
if (initvalues.find("f1QCD") == initvalues.end())initvalues["f1QCD"] = .3;
if (initvalues.find("cE_QCD") == initvalues.end())initvalues["cE_QCD"] = -0.001;
if (initvalues.find("sigmaG_QCD") == initvalues.end())initvalues["sigmaG_QCD"] = 5;
if (initvalues.find("muG_QCD") == initvalues.end())initvalues["muG_QCD"] = 0;
if (initvalues.find("f1QCD") == initvalues.end())initvalues["f1QCD"] = .95;
if (initvalues.find("f2QCD") == initvalues.end())initvalues["f2QCD"] = .3;
if (initvalues.find("fInvIso") == initvalues.end())initvalues["fInvIso"] = .9;
if (initvalues.find("f") == initvalues.end())initvalues["f"] = .9;
//if (initvalues.find("XXX") == initvalues.end())initvalues["XXX"] = 40;
// x) create QCD model
RooRealVar sR("sigmaR_QCD","sigmaR",initvalues["sigmaR_QCD"]);
RooGenericPdf rayleighQCD ("rayleigh_QCD","@0/(@1*@1)*TMath::Exp(-@0*@0/(2*@1*@1))",RooArgList(x,sR) );
RooRealVar cE("cE_QCD","cE",initvalues["cE_QCD"]);
RooRealVar cE("cE_QCD","cE",initvalues["cE_QCD"],-1.,0.);
RooExponential expoQCD("expo_QCD","expoQCD",x,cE);
RooRealVar sG("sigmaG_QCD","sigmaG",initvalues["sigmaG_QCD"]);
RooRealVar mG("muG_QCD","muG",initvalues["muG_QCD"]);
RooGaussian gausQCD("gauss_QCD","gauss_QCD",x,mG,sG);
RooRealVar sG("sigmaG_QCD","sigmaG",initvalues["sigmaG_QCD"],0.00001,30);
RooRealVar mG("muG_QCD","muG",initvalues["muG_QCD"],-5,5);
RooGaussian gaussQCD("gauss_QCD","gauss_QCD",x,mG,sG);

RooRealVar f1QCD("f1QCD","f1QCD",initvalues["f1QCD"]);
RooRealVar f1QCD("f1QCD","f1QCD",initvalues["f1QCD"],.9,1.0);
RooRealVar f2QCD("f2QCD","f2QCD",initvalues["f2QCD"]);
//RooAddPdf modelQCD("modelQCD","modelQCD",RooArgList(rayleighQCD,expoQCD,gausQCD) , RooArgList(f1QCD,f2QCD) );
//RooAddPdf modelQCD("modelQCD","modelQCD",RooArgList(rayleighQCD,expoQCD,gaussQCD) , RooArgList(f1QCD,f2QCD) );
//RooGenericPdf modelQCD ("rayleigh_QCD","@0/(@1*@1)*TMath::Exp(-@0*@0/(2*@1*@1))",RooArgList(x,sR) );
RooAddPdf modelQCD("modelQCD","modelQCD",RooArgList(rayleighQCD,expoQCD) , RooArgList(f1QCD) );
//RooAddPdf modelQCD("modelQCD","modelQCD",RooArgList(rayleighQCD,expoQCD) , RooArgList(f1QCD) );
RooFFTConvPdf modelQCD("rxg","rayleigh (X) gauss",x,rayleighQCD,gaussQCD);

RooRealVar fInvIso("fInvIso","fInvIso",initvalues["fInvIso"]);
RooAddPdf modelInvIso("modelInvIso","modelInvIso",RooArgList(modelQCD,modelEwk) , RooArgList(fInvIso) );
Expand Down

0 comments on commit b4851c6

Please sign in to comment.