Skip to content

Commit 345aacb

Browse files
clementhelsensclementhelsens
committed
udapte
1 parent 81f6e65 commit 345aacb

File tree

10 files changed

+109
-67
lines changed

10 files changed

+109
-67
lines changed

‎analyzers/dataframe/myUtils.cc‎

Lines changed: 8 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -279,8 +279,9 @@ myUtils::get_VertexObject(ROOT::VecOps::RVec<VertexingUtils::FCCAnalysesVertexMC
279279
if (recoparticles.size()<2)continue;
280280

281281
VertexingUtils::FCCAnalysesVertex TheVertex;
282-
if (v==0) TheVertex = VertexFitterSimple::VertexFitter(1,recoparticles, tracks );
283-
else TheVertex = VertexFitterSimple::VertexFitter(0,recoparticles, tracks );
282+
if (v==0) TheVertex = VertexFitterSimple::VertexFitter(1,recoparticles, tracks, true, 4.5, 20e-3, 300 );
283+
//if (v==0) TheVertex = VertexFitterSimple::VertexFitter(1,recoparticles, tracks);
284+
else TheVertex = VertexFitterSimple::VertexFitter(0,recoparticles, tracks);
284285
if (std::isnan(TheVertex.vertex.chi2))continue;
285286
TheVertex.mc_ind=v;
286287
result.push_back(TheVertex);
@@ -1441,7 +1442,7 @@ myUtils::build_composite_vertex::operator() (ROOT::VecOps::RVec<edm4hep::Reconst
14411442
recoparticles.push_back(recop.at(tmpvec.at(k)));
14421443
}
14431444

1444-
VertexingUtils::FCCAnalysesVertex TheVertex = VertexFitterSimple::VertexFitter(0,recoparticles, tracks );
1445+
VertexingUtils::FCCAnalysesVertex TheVertex = VertexFitterSimple::VertexFitter(0,recoparticles, tracks);
14451446
float chi2 = TheVertex.vertex.chi2;
14461447

14471448
if (chi2<0.01)continue;
@@ -1520,7 +1521,7 @@ myUtils::build_D0::operator() (ROOT::VecOps::RVec<edm4hep::ReconstructedParticle
15201521
ROOT::VecOps::RVec<edm4hep::ReconstructedParticleData> seltracks;
15211522
seltracks.push_back(recop.at(pions.at(i)));
15221523
seltracks.push_back(recop.at(kaons.at(j)));
1523-
VertexingUtils::FCCAnalysesVertex TheVertex = VertexFitterSimple::VertexFitter(0,seltracks, tracks );
1524+
VertexingUtils::FCCAnalysesVertex TheVertex = VertexFitterSimple::VertexFitter(0,seltracks, tracks);
15241525
float chi2 = TheVertex.vertex.chi2;
15251526

15261527
if (chi2<0.01)continue;
@@ -1630,7 +1631,7 @@ ROOT::VecOps::RVec<float> myUtils::awkwardtest(ROOT::VecOps::RVec<edm4hep::Recon
16301631
ROOT::VecOps::RVec<int> tk_ind;
16311632

16321633
ROOT::VecOps::RVec<edm4hep::ReconstructedParticleData> seltracks = VertexingUtils::selTracks(0.,3.,0.,3.)( recop, tracks);
1633-
VertexingUtils::FCCAnalysesVertex ThePVertex = VertexFitterSimple::VertexFitter(0,seltracks, tracks );
1634+
VertexingUtils::FCCAnalysesVertex ThePVertex = VertexFitterSimple::VertexFitter(0,seltracks, tracks);
16341635

16351636
int PV_ntrk = ThePVertex.ntracks;
16361637
float PV_chi2 = ThePVertex.vertex.chi2;
@@ -1750,7 +1751,7 @@ ROOT::VecOps::RVec<float> myUtils::awkwardtest(ROOT::VecOps::RVec<edm4hep::Recon
17501751
}
17511752

17521753
//VertexingUtils::FCCAnalysesVertex TheVertexActs = VertexFitterActs::VertexFitterFullBilloir(recoparticles, tracks );
1753-
VertexingUtils::FCCAnalysesVertex TheVertex = VertexFitterSimple::VertexFitter(0,recoparticles, tracks );
1754+
VertexingUtils::FCCAnalysesVertex TheVertex = VertexFitterSimple::VertexFitter(0,recoparticles, tracks);
17541755
float chi2 = TheVertex.vertex.chi2;
17551756

17561757
if (chi2<0.01)continue;
@@ -1991,7 +1992,7 @@ myUtils::build_tau23pi_vertexing::operator() (ROOT::VecOps::RVec<edm4hep::Recons
19911992
recoparticles.push_back(recop.at(tmpvec.at(k)));
19921993
}
19931994

1994-
VertexingUtils::FCCAnalysesVertex TheVertex = VertexFitterSimple::VertexFitter(0,recoparticles, tracks );
1995+
VertexingUtils::FCCAnalysesVertex TheVertex = VertexFitterSimple::VertexFitter(0,recoparticles, tracks);
19951996
float chi2 = TheVertex.vertex.chi2;
19961997

19971998
if (chi2<0.01)continue;

‎examples/FCCee/flavour/Bc2TauNu/analysis_stage1.py‎

Lines changed: 2 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -23,7 +23,6 @@
2323

2424
ROOT.gInterpreter.ProcessLine('''
2525
TMVA::Experimental::RBDT<> bdt("Bc2TauNu_BDT", "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/xgb_bdt_vtx.root");
26-
2726
computeModel = TMVA::Experimental::Compute<18, float>(bdt);
2827
''')
2928

@@ -42,7 +41,7 @@ def __init__(self, inputlist, outname, ncpu):
4241
print (" init done, about to run")
4342
#__________________________________________________________
4443
def run(self):
45-
#df2 = (self.df.Range(10)
44+
#df2 = (self.df.Range(100)
4645
df2 = (self.df
4746
#############################################
4847
## Aliases for # in python ##
@@ -376,7 +375,7 @@ def run(self):
376375
df2.Snapshot("events", self.outname, branchList)
377376

378377
# example call for standalone file
379-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage1.py p8_ee_Zbb_Bc2TauNu_stage1.root /eos/experiment/fcc/ee/generation/DelphesEvents/spring2021/IDEA/p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU/events_003834121.root
378+
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage1.py p8_ee_Zbb_Bc2TauNu_stage1.root /eos/experiment/fcc/ee/generation/DelphesEvents/spring2021/IDEA/p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU/events_003907469.root
380379

381380
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage1.py p8_ee_Zbb_Bu2TauNu_stage1.root /eos/experiment/fcc/ee/generation/DelphesEvents/fcc_tmp_v03/p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU/events_026079857.root
382381

‎examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py‎

Lines changed: 30 additions & 24 deletions
Original file line numberDiff line numberDiff line change
@@ -22,7 +22,12 @@
2222

2323
ROOT.gInterpreter.ProcessLine('''
2424
TMVA::Experimental::RBDT<> bdt("Bc2TauNu_BDT2", "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/xgb_bdt_stage2.root");
25-
computeModel = TMVA::Experimental::Compute<20, float>(bdt);
25+
computeModel = TMVA::Experimental::Compute<21, float>(bdt);
26+
''')
27+
28+
ROOT.gInterpreter.ProcessLine('''
29+
TMVA::Experimental::RBDT<> bdt("Bc2TauNu_BDT", "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/xgb_bdt_vtx.root");
30+
computeModelBis = TMVA::Experimental::Compute<18, float>(bdt);
2631
''')
2732

2833

@@ -58,7 +63,7 @@ def run(self):
5863

5964
.Define("LOCAL_CandVtxInd", "Tau23PiCandidates_vertex.at(CUT_CandInd)")
6065
##LOCAL INDEX screwed up in prod02!!!! need -1 because PV is removed
61-
.Define("CUT_CandVtxThrustEmin", "Vertex_thrusthemis_emin.at(LOCAL_CandVtxInd-1)")
66+
.Define("CUT_CandVtxThrustEmin", "Vertex_thrusthemis_emin.at(LOCAL_CandVtxInd)")
6267
.Filter("CUT_CandVtxThrustEmin>0")
6368

6469
.Define("EVT_CandN", "float(Tau23PiCandidates_vertex.size())")
@@ -103,12 +108,30 @@ def run(self):
103108
.Define("EVT_DVmass_ave", "myFinalSel::get_ave(Vertex_mass, Vertex_isPV, LOCAL_CandVtxInd)")
104109
.Define("EVT_PVmass", "Vertex_mass.at(0)")
105110

106-
.Define("MVAVec", ROOT.computeModel, ("EVT_CandMass","EVT_CandRho1Mass","EVT_CandRho2Mass","EVT_CandN","EVT_CandVtxFD","EVT_CandVtxChi2","EVT_CandPx","EVT_CandPy","EVT_CandPz","EVT_CandP","EVT_CandD0","EVT_CandZ0","EVT_CandAngleThrust","EVT_DVd0_min","EVT_DVd0_max","EVT_DVd0_ave","EVT_DVz0_min","EVT_DVz0_max","EVT_DVz0_ave","EVT_Nominal_B_E"))
107-
111+
.Define("MVAVec", ROOT.computeModel, ("EVT_CandMass", "EVT_CandRho1Mass", "EVT_CandRho2Mass",
112+
"EVT_CandN", "EVT_CandVtxFD", "EVT_CandVtxChi2",
113+
"EVT_CandPx", "EVT_CandPy", "EVT_CandPz",
114+
"EVT_CandP", "EVT_CandD0", "EVT_CandZ0",
115+
"EVT_CandAngleThrust", "EVT_DVd0_min", "EVT_DVd0_max",
116+
"EVT_DVd0_ave", "EVT_DVz0_min", "EVT_DVz0_max",
117+
"EVT_DVz0_ave", "EVT_PVmass", "EVT_Nominal_B_E"))
118+
119+
108120
.Define("EVT_MVA2", "MVAVec.at(0)")
109121
.Filter(MVAFilter)
110122

111-
123+
.Define("MVAVecBis", ROOT.computeModelBis, ("EVT_ThrustEmin_E", "EVT_ThrustEmax_E",
124+
"EVT_ThrustEmin_Echarged", "EVT_ThrustEmax_Echarged",
125+
"EVT_ThrustEmin_Eneutral", "EVT_ThrustEmax_Eneutral",
126+
"EVT_ThrustEmin_Ncharged", "EVT_ThrustEmax_Ncharged",
127+
"EVT_ThrustEmin_Nneutral", "EVT_ThrustEmax_Nneutral",
128+
"EVT_NtracksPV", "EVT_NVertex",
129+
"EVT_NTau23Pi", "EVT_ThrustEmin_NDV",
130+
"EVT_ThrustEmax_NDV", "EVT_dPV2DVmin",
131+
"EVT_dPV2DVmax", "EVT_dPV2DVave"))
132+
.Define("EVT_MVA1Bis", "MVAVecBis.at(0)")
133+
134+
112135
.Define("EVT_minRhoMass", "if (EVT_CandRho1Mass<EVT_CandRho2Mass) return EVT_CandRho1Mass; else return EVT_CandRho2Mass;")
113136
.Define("EVT_maxRhoMass", "if (EVT_CandRho1Mass>EVT_CandRho2Mass) return EVT_CandRho1Mass; else return EVT_CandRho2Mass;")
114137
.Define("EVT_ThrustDiff_E", "EVT_ThrustEmax_E-EVT_ThrustEmin_E")
@@ -156,7 +179,7 @@ def run(self):
156179
"EVT_CandPx","EVT_CandP","EVT_CandPz","EVT_CandPy",
157180
"EVT_CandD0","EVT_CandZ0","EVT_CandAngleThrust",
158181
"EVT_minRhoMass", "EVT_maxRhoMass",
159-
"EVT_MVA1", "EVT_MVA2",
182+
"EVT_MVA1", "EVT_MVA2","EVT_MVA1Bis",
160183

161184
"EVT_CandPion1P","EVT_CandPion1D0","EVT_CandPion1Z0",
162185
"EVT_CandPion2P","EVT_CandPion2D0","EVT_CandPion2Z0",
@@ -181,28 +204,11 @@ def run(self):
181204
#df2.Snapshot("events", self.outname, branchList, opts)
182205
df2.Snapshot("events", self.outname, branchList)
183206

184-
# example call for standalone file
185-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Batch_Analysis_stage1/p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU/*.root"
186-
187-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Batch_Analysis_stage1/p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU/*.root"
188-
189-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Batch_Analysis_stage1/p8_ee_Zbb_ecm91/*.root"
190-
191-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zcc_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Batch_Analysis_stage1/p8_ee_Zcc_ecm91/*.root"
192-
193-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zuds_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Batch_Analysis_stage1/p8_ee_Zuds_ecm91/*.root"
194207

195208
#########################################
196209

197-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Batch_Training_4stage1/p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU/*.root"
198-
199-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Batch_Training_4stage1/p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU/*.root"
200-
201-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Batch_Training_4stage1/p8_ee_Zbb_ecm91/*.root"
202-
203-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zcc_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Batch_Training_4stage1/p8_ee_Zcc_ecm91/*.root"
210+
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU.root /eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_03/Batch_Training_4stage1/p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU/flat_chunk_0.root
204211

205-
# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zuds_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Batch_Training_4stage1/p8_ee_Zuds_ecm91/*.root"
206212

207213
if __name__ == "__main__":
208214

‎examples/FCCee/flavour/Bc2TauNu/finalSel.py‎

Lines changed: 19 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -5,18 +5,15 @@
55
import ROOT
66

77
###Input directory where the files produced at the pre-selection level are
8-
baseDir = "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Analysis_stage2/"
8+
baseDir = "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_04/Analysis_stage2/"
99

1010
###Link to the dictonary that contains all the cross section informations etc...
1111
procDict = os.path.join(os.getenv('FCCDICTSDIR', deffccdicts), '') + "FCCee_procDict_spring2021_IDEA.json"
1212

1313

14-
#baseDir = "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Training_4stage2/"
15-
#procDict = os.path.join(os.getenv('FCCDICTSDIR', deffccdicts), '') + "FCCee_procDict_fcc_tmp_training.json"
16-
17-
1814
process_list=['p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU',
1915
'p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU',
16+
'p8_ee_Zbb_ecm91_EvtGen',
2017
'p8_ee_Zbb_ecm91',
2118
'p8_ee_Zcc_ecm91',
2219
'p8_ee_Zuds_ecm91',
@@ -54,7 +51,6 @@
5451
'p8_ee_Zbb_ecm91_EvtGen_Lb2LcstTauNu',
5552
]
5653

57-
#process_list=['p8_ee_Zbb_ecm91_EvtGen_Bs2DsstDsst']
5854

5955
define_list={
6056
# "EVT_minRhoMass":"if (EVT_CandRho1Mass<EVT_CandRho2Mass) return EVT_CandRho1Mass; else return EVT_CandRho2Mass;",
@@ -63,13 +59,15 @@
6359
}
6460

6561
###Dictionnay of the list of cuts. The key is the name of the selection that will be added to the output file
66-
cut_list = {"sel0":"EVT_MVA1>0.8 && EVT_MVA2>0.8",
62+
cut_list = {#"sel0":"EVT_MVA1>0.8 && EVT_MVA2>0.8",
6763
#"sel1":"EVT_MVA1>0.9 && EVT_MVA2>0.9",
6864
#"sel2":"EVT_MVA1>0.95 && EVT_MVA2>0.95",
6965
#"sel3":"EVT_MVA1>0.98 && EVT_MVA2>0.98",
7066
#"sel4":"EVT_MVA1>0.99 && EVT_MVA2>0.99",
67+
"sel5":"EVT_MVA1>0.9998979591836735 && EVT_MVA2>0.9956122448979592 && EVT_ThrustDiff_E > 10",
7168
}
7269

70+
7371
###Dictionary for the ouput variable/hitograms. The key is the name of the variable in the output files. "name" is the name of the variable in the input file, "title" is the x-axis label of the histogram, "bin" the number of bins of the histogram, "xmin" the minimum x-axis value and "xmax" the maximum x-axis value.
7472
variables = {
7573

@@ -97,11 +95,23 @@
9795
"EVT_CandAngleThrust_2" :{"name":"EVT_CandAngleThrust","title":"Angle between candidate and thrust","bin":50,"xmin":0.,"xmax":0.5},
9896
"EVT_MVA1" :{"name":"EVT_MVA1","title":"MVA1 score","bin":100,"xmin":0.6,"xmax":1.},
9997
"EVT_MVA2" :{"name":"EVT_MVA2","title":"MVA2 score","bin":100,"xmin":0.6,"xmax":1.},
100-
98+
"EVT_Nominal_B_E" :{"name":"EVT_Nominal_B_E","title":"Nominal B Energy [GeV]","bin":100,"xmin":0.,"xmax":50.},
99+
"EVT_PVmass" :{"name":"EVT_PVmass","title":"Primary Vertex Mass [GeV]","bin":100,"xmin":0.,"xmax":10.},
100+
"EVT_DVmass_min" :{"name":"EVT_DVmass_min","title":"Displaced Vertex Mass Min [GeV]","bin":100,"xmin":0.,"xmax":20.},
101+
"EVT_DVmass_max" :{"name":"EVT_DVmass_max","title":"Displaced Vertex Mass Max [GeV]","bin":100,"xmin":0.,"xmax":20.},
102+
"EVT_DVmass_ave" :{"name":"EVT_DVmass_ave","title":"Displaced Vertex Mass Ave [GeV]","bin":100,"xmin":0.,"xmax":20.},
103+
"EVT_DVd0_min" :{"name":"EVT_DVd0_min","title":"Displaced Vertex d0 Min [mm]","bin":100,"xmin":-10.,"xmax":10.},
104+
"EVT_DVd0_max" :{"name":"EVT_DVd0_max","title":"Displaced Vertex d0 Max [mm]","bin":100,"xmin":-10.,"xmax":10.},
105+
"EVT_DVd0_ave" :{"name":"EVT_DVd0_ave","title":"Displaced Vertex d0 Ave [mm]","bin":100,"xmin":-10.,"xmax":10.},
106+
"EVT_DVz0_min" :{"name":"EVT_DVz0_min","title":"Displaced Vertex z0 Min [mm]","bin":100,"xmin":-10.,"xmax":10.},
107+
"EVT_DVz0_max" :{"name":"EVT_DVz0_max","title":"Displaced Vertex z0 Max [mm]","bin":100,"xmin":-10.,"xmax":10.},
108+
"EVT_DVz0_ave" :{"name":"EVT_DVz0_ave","title":"Displaced Vertex z0 Ave [mm]","bin":100,"xmin":-10.,"xmax":10.},
109+
110+
101111
}
102112

103113
###Number of CPUs to use
104-
NUM_CPUS = 4
114+
NUM_CPUS = 8
105115

106116
###This part is standard to all analyses
107117
import config.runDataFrameFinal as rdf

‎examples/FCCee/flavour/Bc2TauNu/plots.py‎

Lines changed: 30 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -6,11 +6,11 @@
66
delphesVersion = "3.4.2"
77
energy = 91.0
88
collider = "FCC-ee"
9-
inputDir = "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Analysis_stage2/"
9+
inputDir = "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_04/Analysis_stage2/"
1010
formats = ['png','pdf']
1111
yaxis = ['lin','log']
1212
stacksig = ['nostack']
13-
outdir = '/eos/home-h/helsens/www/FCC/ee/flavour/Bc2TauNu/spring2021/prod_01/'
13+
outdir = '/eos/home-h/helsens/www/FCC/ee/flavour/Bc2TauNu/spring2021/prod_04/plots/'
1414
variables = [ "EVT_CandMass",
1515
"EVT_ThrustEmin_E","EVT_ThrustEmax_E",
1616
"EVT_ThrustEmin_Echarged", "EVT_ThrustEmax_Echarged",
@@ -25,28 +25,45 @@
2525
"EVT_CandVtxFD",
2626
"EVT_CandAngleThrust",
2727
"EVT_CandAngleThrust_2",
28-
"EVT_MVA1","EVT_MVA2"
28+
"EVT_MVA1","EVT_MVA2",
29+
30+
"EVT_Nominal_B_E",
31+
"EVT_PVmass",
32+
"EVT_DVmass_min",
33+
"EVT_DVmass_max",
34+
"EVT_DVmass_ave",
35+
"EVT_DVd0_min",
36+
"EVT_DVd0_max",
37+
"EVT_DVd0_ave",
38+
"EVT_DVz0_min",
39+
"EVT_DVz0_max",
40+
"EVT_DVz0_ave",
2941

42+
3043
]
3144

3245
scaleSig=1.
3346
###Dictonnary with the analysis name as a key, and the list of selections to be plotted for this analysis. The name of the selections should be the same than in the final selection
3447
selections = {}
3548
selections['Inclusive'] = ["sel1","sel2","sel3","sel4"]
3649
selections['ExclusiveMerged'] = ["sel1","sel2","sel3","sel4"]
37-
selections['Inclusive'] = ["sel0"]
38-
selections['ExclusiveMerged'] = ["sel0"]
50+
51+
selections['Inclusive'] = ["sel5"]
52+
selections['InclusiveEvtGen'] = ["sel5"]
53+
selections['ExclusiveMerged'] = ["sel5"]
3954

4055
extralabel = {}
4156
extralabel['sel0'] = "Selection: MVA1>0.8, MVA2>0.8"
4257
extralabel['sel1'] = "Selection: MVA1>0.9, MVA2>0.9"
4358
extralabel['sel2'] = "Selection: MVA1>0.95, MVA2>0.95"
4459
extralabel['sel3'] = "Selection: MVA1>0.98, MVA2>0.98"
4560
extralabel['sel4'] = "Selection: MVA1>0.99, MVA2>0.99"
61+
extralabel['sel5'] = "Selection: Donal's cuts"
4662

4763
colors = {}
4864
colors['Z_Bc'] = ROOT.kRed
4965
colors['Z_bb'] = ROOT.kBlue
66+
colors['Z_bb_EvtGen'] = ROOT.kBlue
5067
colors['Z_bb_B'] = ROOT.kBlue
5168
colors['Z_bb_Bu'] = ROOT.kBlue+1
5269
colors['Z_bb_Bd'] = ROOT.kBlue+2
@@ -65,6 +82,13 @@
6582
}
6683
}
6784

85+
plots['InclusiveEvtGen'] = {'signal':{'Z_Bc':['p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU']},
86+
'backgrounds':{'Z_bb_EvtGen':['p8_ee_Zbb_ecm91_EvtGen'],
87+
'Z_cc':['p8_ee_Zcc_ecm91'],
88+
'Z_uds':['p8_ee_Zuds_ecm91'],
89+
'Z_Bu':['p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU']
90+
}
91+
}
6892

6993

7094
plots['ExclusiveMerged'] = {'signal':{'Z_Bc':['p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU']},
@@ -109,6 +133,7 @@
109133
legend = {}
110134
legend['Z_Bc'] = 'B_{c}#rightarrow #tau#nu (#tau#rightarrow 3#pi)'
111135
legend['Z_bb'] = 'Z#rightarrow b#bar{b}'
136+
legend['Z_bb_EvtGen'] = 'Z#rightarrow b#bar{b} EvtGen'
112137
legend['Z_bb_B'] = 'Z#rightarrow b#bar{b} B_{u}+B_{d}+B_{s}+L_{b}'
113138
legend['Z_cc'] = 'Z#rightarrow c#bar{c}'
114139
legend['Z_uds'] = 'Z#rightarrow q#bar{q}'

0 commit comments

Comments
 (0)