MultivariateAnalysis/root/qcd_purity.C

{ gROOT->SetStyle("Plain");
  gStyle->SetOptStat(0000); 
  gStyle->SetPalette(1);

  TChain ch_ttbar("ttljets");
  //TChain ch_ttbar("mva/ttljets");
  TChain ch_w("ttljets");
  TChain ch_z("ttljets");
  TChain ch_tW("ttljets"); // single top
  TChain ch_t("ttljets");  // single top
  TChain ch_s("ttljets");  // single top
  TChain ch_qcd("ttljets");
  //
  ch_ttbar . Add("TTJets-madgraph_Fall08_IDEAL_V9_v2_14may2009.root");
  //ch_ttbar . Add("TTJets_jes_up_21may2009.root");
  //ch_ttbar . Add("TTJets_jes_down_21may2009.root");
  //ch_ttbar . Add("Tauola_23mar2009.root");
  //ch_ttbar . Add("TTJets_LargerISRFSR_23may2009.root");
  //ch_ttbar . Add("TTJets_SmallerISRFSR_23may2009.root");
  //
   //ch_w . Add("Wjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
  ch_w . Add("WJets-madgraph_Fall08_IDEAL_V9_v1.root");
  //ch_w . Add("/uscms_data/d2/lpcljm/MVA/Summer08/May_28_2009/WJets_ScaleUp-madgraph_Winter09/Muons/all.root");
  //ch_w . Add("/uscms_data/d2/lpcljm/MVA/Summer08/May_28_2009/WJets_ScaleDown-madgraph_Winter09/Muons/all.root");
  //ch_w . Add("WJets_ScaleDown_27may2009.root");
  //ch_w . Add("WJets_Threshold20GeV_22may2009.root");
  //ch_w . Add("WJets_Threshold5GeV_22may2009.root");
  //
   //ch_z . Add("Zjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
  ch_z . Add("ZJets-madgraph_Fall08_IDEAL_V9_reco-v2.root");
  //ch_z . Add("Zjets_jes_up_21may2009.root");
  //ch_z . Add("Zjets_jes_down_21may2009.root");
  //
  ch_t . Add("SingleTop_tChannel_Summer08_IDEAL_V9_v1.root");
  ch_tW. Add("SingleTop_tWChannel_Summer08_IDEAL_V9_v1.root");
  ch_s . Add("SingleTop_sChannel_Summer08_IDEAL_V9_v1.root");
  //
  ch_qcd . Add("InclusiveMuPt15_Summer08_IDEAL_V9_v1_14may2009.root");
  //ch_qcd . Add("InclusiveMuPt15_jes_up_21may2009.root");
  //ch_qcd . Add("InclusiveMuPt15_jes_down_21may2009.root");
  //
  //_____ systematics samples
  //
  //ch_w . Add("WJets_ScaleUp_22may2009.root");
  //ch_w . Add("WJets_ScaleDown_22may2009.root");
  //ch_w . Add("WJets_Threshold20GeV_22may2009.root");
  //ch_w . Add("WJets_Threshold5GeV_22may2009.root");

#include "cuts.C"

  //the_cut = "";
  //the_cut = cut3;
  //the_cut = cut3_rev;
  //the_cut = cut3_qcd;
   the_cut = cut4;
   //the_cut = "getHt<200";
  //
  //the_cut = cut4_abcd;
   //the_cut = cut4_abcd && rel_iso_cut4 && lepton_d0sig_cut4;
   //the_cut = cut4_abcd && rel_iso_cut4 && lepton_d0sig_cut4_abcd;
   //the_cut = cut4_abcd && rel_iso_cut4_abcd && lepton_d0sig_cut4;
   //the_cut = cut4_abcd && rel_iso_cut4_abcd && lepton_d0sig_cut4_abcd;
  //
  //the_cut = cut4_qcd;
  //the_cut = muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && rel_iso_cut4 && lepton_IsoDeposit_em_cut4 && lepton_IsoDeposit_had_cut4 && lepton_d0sig_cut4;

  double lumi = 20.0;

  //double xsec_ttbar = 317.0;
  double xsec_ttbar = 414.0;
  double xsec_w   = 40000.0;
  double xsec_z   =  3700.0;
  double xsec_t   =   130.0;
  double xsec_tW  =    29.0;
  double xsec_s   =     5.0;
  double xsec_qcd =     0.5091e+9;

  double eff_gen_ttbar = 1.0;
  double eff_gen_w = 1.0;
  double eff_gen_z = 1.0;
  double eff_gen_tW = 1.0; //0.3257;
  double eff_gen_t  = 0.3257;
  double eff_gen_s  = 0.3257;
  double eff_gen_qcd = 0.000239;

  Long64_t n_processed_ttbar = 1018924; // TTJets madgraph
  //Long64_t n_processed_ttbar = 1334225; // TTJets ISRFSR high
  //Long64_t n_processed_ttbar = 1134190; // TTJets ISRFSR low
  //Long64_t n_processed_ttbar = 106540;  //Tauola
   //
   //Long64_t n_processed_w     = 101467739; // fast sim
  Long64_t n_processed_w       = 10110648; // full sim
  //Long64_t n_processed_w     = 14287301; // Wjets Q^2 scale up
  //Long64_t n_processed_w     = 11444008; // Wjets Q^2 scale down
  //Long64_t n_processed_w     = 21939653; //20 Gev Threshold
  //Long64_t n_processed_w     = 15480461; //5 Gev Threshold
  //
  //Long64_t n_processed_z     = 9679459; // fast sim PAT: 2000116
   Long64_t n_processed_z     = 1163479; // full sim
  Long64_t n_processed_tW    = 169048;
  Long64_t n_processed_t     = 281756;
  Long64_t n_processed_s     = 11999;
  //
  Long64_t n_processed_qcd   = 6142581; // PAT: 1282829

  cout << "Selecting events..." << endl;
  Long64_t n_selected_ttbar  = ch_ttbar  -> Draw("event", the_cut);
  Long64_t n_selected_w      = ch_w      -> Draw("event", the_cut);
  Long64_t n_selected_z      = ch_z      -> Draw("event", the_cut);
  Long64_t n_selected_t      = ch_t      -> Draw("event", the_cut);
  Long64_t n_selected_tW     = ch_tW     -> Draw("event", the_cut);
  Long64_t n_selected_s      = ch_s      -> Draw("event", the_cut);
  Long64_t n_selected_qcd    = ch_qcd    -> Draw("event", the_cut);

  //Selection efficiency, w events: 0.000166854
  //Selection efficiency, z events: 0.00033864
  double eff_sel_ttbar = (double)n_selected_ttbar/(double)n_processed_ttbar;
   double eff_sel_w = (double)n_selected_w/(double)n_processed_w;
   double eff_sel_z = (double)n_selected_z/(double)n_processed_z;
   //double eff_sel_w = 0.000166854;
  //double eff_sel_z = 0.00033864;
  double eff_sel_t = (double)n_selected_t/(double)n_processed_t;
  double eff_sel_tW= (double)n_selected_tW/(double)n_processed_tW;
  double eff_sel_s = (double)n_selected_s/(double)n_processed_s;
  double eff_sel_qcd = (double)n_selected_qcd/(double)n_processed_qcd;

  double n_expected_ttbar = xsec_ttbar * eff_gen_ttbar * eff_sel_ttbar * lumi;
  double n_expected_w = xsec_w * eff_gen_w * eff_sel_w * lumi;
  double n_expected_z = xsec_z * eff_gen_z * eff_sel_z * lumi;
  double n_expected_t = xsec_t * eff_gen_t * eff_sel_t * lumi;
  double n_expected_tW= xsec_tW* eff_gen_tW* eff_sel_tW* lumi;
  double n_expected_s = xsec_s * eff_gen_s * eff_sel_s * lumi;
  double n_expected_qcd = xsec_qcd * eff_gen_qcd * eff_sel_qcd * lumi;

  double n_expected_all = n_expected_ttbar + n_expected_w + n_expected_z + n_expected_t + n_expected_tW + n_expected_s + n_expected_qcd;
  double n_expected_phys = n_expected_ttbar + n_expected_w + n_expected_z + n_expected_t + n_expected_tW + n_expected_s;
  double n_expected_wzt = n_expected_w + n_expected_z + n_expected_t + n_expected_tW + n_expected_s;
  //double n_expected_all = n_expected_ttbar + n_expected_z + n_expected_qcd;
  //double n_expected_phys = n_expected_ttbar + n_expected_z;

  //
  // _____ mixing phys bg _______________________________________________
  //
  double eff_lumi_ttbar = n_processed_ttbar/xsec_ttbar/eff_gen_ttbar;
  double eff_lumi_w = n_processed_w/xsec_w/eff_gen_w;
  double eff_lumi_z = n_processed_z/xsec_z/eff_gen_z;
  double eff_lumi_t = n_processed_t/xsec_t/eff_gen_t;
  double eff_lumi_tW= n_processed_tW/xsec_tW/eff_gen_tW;
  double eff_lumi_s = n_processed_s/xsec_s/eff_gen_s;
  //
  double w_rate_corr = 133.483/103.891;
  double z_rate_corr = 25.0594/18.5393;
  //
  double n_template_w = (double)n_selected_w*eff_lumi_w*w_rate_corr/eff_lumi_w/w_rate_corr;
  double n_template_z = (double)n_selected_z*eff_lumi_w*z_rate_corr/eff_lumi_z/w_rate_corr;
  double n_template_t = (double)n_selected_t*eff_lumi_w/eff_lumi_t/w_rate_corr;
  double n_template_tW= (double)n_selected_tW*eff_lumi_w/eff_lumi_tW/w_rate_corr;
  double n_template_s = (double)n_selected_s*eff_lumi_w/eff_lumi_s/w_rate_corr;


  cout << endl << "Selection efficiency, ttbar events: " << eff_sel_ttbar << endl;
  cout << "Selected ttbar events: " << n_selected_ttbar << endl;
  cout << "Expected ttbar events: " << n_expected_ttbar << endl;
  cout << "Effective ttbar lumi: " << eff_lumi_ttbar << endl;
  //
  cout << endl << "Selection efficiency, w events: " << eff_sel_w << endl;
  cout << "Selected w events: " << n_selected_w << endl;
  cout << "Expected w events: " << n_expected_w << endl;
  cout << "Effective w lumi: " << eff_lumi_w << endl;
  cout << "Events for the template: " << n_template_w << endl;
  //
  cout << endl << "Selection efficiency, z events: " << eff_sel_z << endl;
  cout << "Selected z events: " << n_selected_z << endl;
  cout << "Expected z events: " << n_expected_z << endl;
  cout << "Effective z lumi: " << eff_lumi_z << endl;
  cout << "Events for the template: " << n_template_z << endl;
  //
  cout << endl << "Selection efficiency, single top t-channel events: " << eff_sel_t << endl;
  cout << "Selected t events: " << n_selected_t << endl;
  cout << "Expected t events: " << n_expected_t << endl;
  cout << "Effective t lumi: " << eff_lumi_t << endl;
  cout << "Events for the template: " << n_template_t << endl;
  //
  cout << endl << "Selection efficiency, single top tW-channel events: " << eff_sel_tW << endl;
  cout << "Selected tW events: " << n_selected_tW << endl;
  cout << "Expected tW events: " << n_expected_tW << endl;
  cout << "Effective tW lumi: " << eff_lumi_tW << endl;
  cout << "Events for the template: " << n_template_tW << endl;
  //
  cout << endl << "Selection efficiency, single top s-channel events: " << eff_sel_s << endl;
  cout << "Selected s events: " << n_selected_s << endl;
  cout << "Expected s events: " << n_expected_s << endl;
  cout << "Effective s lumi: " << eff_lumi_s << endl;
  cout << "Events for the template: " << n_template_s << endl;
  //
  cout << endl << "Selection efficiency, qcd events: " << eff_sel_qcd << endl;
  cout << "Selected qcd events: " << n_selected_qcd << endl;
  cout << "Expected qcd events: " << n_expected_qcd << endl;

  cout << endl << "Total events expected: " << n_expected_all << endl;
  cout << "Ttbar,W,Z+jets events expected: " << n_expected_phys << endl;
  cout << "ttbar events expected: " << n_expected_ttbar << endl;
  cout << "phys BG events expected: " << n_expected_wzt << endl;
  cout << "QCD events expected: " << n_expected_qcd << endl;


}

 
Revision:	1.13
Committed:	Sat Jul 4 13:29:09 2009 UTC (15 years, 10 months ago) by kukartse
Content type:	text/plain
Branch:	MAIN
CVS Tags:	gak010310, ejterm2010_25nov2009, V00-02-01, V00-02-00, gak112409, CMSSW_22X_branch_base, segala101609
Branch point for:	CMSSW_22X_branch
Changes since 1.12:	+1 -0 lines
Log Message:	* empty log message *
#	User	Rev	Content
1	kukartse	1.1	{ gROOT->SetStyle("Plain");
2			gStyle->SetOptStat(0000);
3			gStyle->SetPalette(1);
4
5	kukartse	1.7	TChain ch_ttbar("ttljets");
6	kukartse	1.9	//TChain ch_ttbar("mva/ttljets");
7	kukartse	1.6	TChain ch_w("ttljets");
8			TChain ch_z("ttljets");
9	kukartse	1.7	TChain ch_tW("ttljets"); // single top
10			TChain ch_t("ttljets"); // single top
11			TChain ch_s("ttljets"); // single top
12			TChain ch_qcd("ttljets");
13	kukartse	1.9	//
14	kukartse	1.7	ch_ttbar . Add("TTJets-madgraph_Fall08_IDEAL_V9_v2_14may2009.root");
15	kukartse	1.9	//ch_ttbar . Add("TTJets_jes_up_21may2009.root");
16			//ch_ttbar . Add("TTJets_jes_down_21may2009.root");
17			//ch_ttbar . Add("Tauola_23mar2009.root");
18			//ch_ttbar . Add("TTJets_LargerISRFSR_23may2009.root");
19			//ch_ttbar . Add("TTJets_SmallerISRFSR_23may2009.root");
20			//
21	kukartse	1.12	//ch_w . Add("Wjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
22			ch_w . Add("WJets-madgraph_Fall08_IDEAL_V9_v1.root");
23	kukartse	1.9	//ch_w . Add("/uscms_data/d2/lpcljm/MVA/Summer08/May_28_2009/WJets_ScaleUp-madgraph_Winter09/Muons/all.root");
24	kukartse	1.10	//ch_w . Add("/uscms_data/d2/lpcljm/MVA/Summer08/May_28_2009/WJets_ScaleDown-madgraph_Winter09/Muons/all.root");
25	kukartse	1.9	//ch_w . Add("WJets_ScaleDown_27may2009.root");
26			//ch_w . Add("WJets_Threshold20GeV_22may2009.root");
27			//ch_w . Add("WJets_Threshold5GeV_22may2009.root");
28			//
29	kukartse	1.12	//ch_z . Add("Zjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
30			ch_z . Add("ZJets-madgraph_Fall08_IDEAL_V9_reco-v2.root");
31	kukartse	1.9	//ch_z . Add("Zjets_jes_up_21may2009.root");
32			//ch_z . Add("Zjets_jes_down_21may2009.root");
33			//
34	kukartse	1.7	ch_t . Add("SingleTop_tChannel_Summer08_IDEAL_V9_v1.root");
35			ch_tW. Add("SingleTop_tWChannel_Summer08_IDEAL_V9_v1.root");
36			ch_s . Add("SingleTop_sChannel_Summer08_IDEAL_V9_v1.root");
37	kukartse	1.9	//
38	kukartse	1.7	ch_qcd . Add("InclusiveMuPt15_Summer08_IDEAL_V9_v1_14may2009.root");
39	kukartse	1.9	//ch_qcd . Add("InclusiveMuPt15_jes_up_21may2009.root");
40			//ch_qcd . Add("InclusiveMuPt15_jes_down_21may2009.root");
41			//
42			//_____ systematics samples
43			//
44			//ch_w . Add("WJets_ScaleUp_22may2009.root");
45			//ch_w . Add("WJets_ScaleDown_22may2009.root");
46			//ch_w . Add("WJets_Threshold20GeV_22may2009.root");
47			//ch_w . Add("WJets_Threshold5GeV_22may2009.root");
48	kukartse	1.1
49			#include "cuts.C"
50
51	kukartse	1.6	//the_cut = "";
52			//the_cut = cut3;
53			//the_cut = cut3_rev;
54			//the_cut = cut3_qcd;
55	kukartse	1.12	the_cut = cut4;
56			//the_cut = "getHt<200";
57	kukartse	1.8	//
58	kukartse	1.9	//the_cut = cut4_abcd;
59	kukartse	1.12	//the_cut = cut4_abcd && rel_iso_cut4 && lepton_d0sig_cut4;
60			//the_cut = cut4_abcd && rel_iso_cut4 && lepton_d0sig_cut4_abcd;
61			//the_cut = cut4_abcd && rel_iso_cut4_abcd && lepton_d0sig_cut4;
62			//the_cut = cut4_abcd && rel_iso_cut4_abcd && lepton_d0sig_cut4_abcd;
63	kukartse	1.8	//
64	kukartse	1.7	//the_cut = cut4_qcd;
65			//the_cut = muon_trigger_cut4 && muon_qual_cut4 && muon_chisq_cut4 && muon_nhits_cut4 && rel_iso_cut4 && lepton_IsoDeposit_em_cut4 && lepton_IsoDeposit_had_cut4 && lepton_d0sig_cut4;
66	kukartse	1.1
67	kukartse	1.2	double lumi = 20.0;
68	kukartse	1.1
69	kukartse	1.7	//double xsec_ttbar = 317.0;
70			double xsec_ttbar = 414.0;
71			double xsec_w = 40000.0;
72			double xsec_z = 3700.0;
73			double xsec_t = 130.0;
74			double xsec_tW = 29.0;
75			double xsec_s = 5.0;
76			double xsec_qcd = 0.5091e+9;
77	kukartse	1.1
78			double eff_gen_ttbar = 1.0;
79			double eff_gen_w = 1.0;
80			double eff_gen_z = 1.0;
81	kukartse	1.11	double eff_gen_tW = 1.0; //0.3257;
82	kukartse	1.7	double eff_gen_t = 0.3257;
83			double eff_gen_s = 0.3257;
84	kukartse	1.6	double eff_gen_qcd = 0.000239;
85	kukartse	1.1
86	kukartse	1.9	Long64_t n_processed_ttbar = 1018924; // TTJets madgraph
87			//Long64_t n_processed_ttbar = 1334225; // TTJets ISRFSR high
88			//Long64_t n_processed_ttbar = 1134190; // TTJets ISRFSR low
89			//Long64_t n_processed_ttbar = 106540; //Tauola
90	kukartse	1.12	//
91			//Long64_t n_processed_w = 101467739; // fast sim
92			Long64_t n_processed_w = 10110648; // full sim
93	kukartse	1.11	//Long64_t n_processed_w = 14287301; // Wjets Q^2 scale up
94	kukartse	1.9	//Long64_t n_processed_w = 11444008; // Wjets Q^2 scale down
95			//Long64_t n_processed_w = 21939653; //20 Gev Threshold
96			//Long64_t n_processed_w = 15480461; //5 Gev Threshold
97			//
98	kukartse	1.12	//Long64_t n_processed_z = 9679459; // fast sim PAT: 2000116
99			Long64_t n_processed_z = 1163479; // full sim
100	kukartse	1.7	Long64_t n_processed_tW = 169048;
101			Long64_t n_processed_t = 281756;
102			Long64_t n_processed_s = 11999;
103	kukartse	1.9	//
104	kukartse	1.6	Long64_t n_processed_qcd = 6142581; // PAT: 1282829
105	kukartse	1.1
106			cout << "Selecting events..." << endl;
107			Long64_t n_selected_ttbar = ch_ttbar -> Draw("event", the_cut);
108			Long64_t n_selected_w = ch_w -> Draw("event", the_cut);
109			Long64_t n_selected_z = ch_z -> Draw("event", the_cut);
110	kukartse	1.7	Long64_t n_selected_t = ch_t -> Draw("event", the_cut);
111			Long64_t n_selected_tW = ch_tW -> Draw("event", the_cut);
112			Long64_t n_selected_s = ch_s -> Draw("event", the_cut);
113	kukartse	1.1	Long64_t n_selected_qcd = ch_qcd -> Draw("event", the_cut);
114
115	kukartse	1.11	//Selection efficiency, w events: 0.000166854
116			//Selection efficiency, z events: 0.00033864
117	kukartse	1.1	double eff_sel_ttbar = (double)n_selected_ttbar/(double)n_processed_ttbar;
118	kukartse	1.12	double eff_sel_w = (double)n_selected_w/(double)n_processed_w;
119			double eff_sel_z = (double)n_selected_z/(double)n_processed_z;
120			//double eff_sel_w = 0.000166854;
121			//double eff_sel_z = 0.00033864;
122	kukartse	1.7	double eff_sel_t = (double)n_selected_t/(double)n_processed_t;
123			double eff_sel_tW= (double)n_selected_tW/(double)n_processed_tW;
124			double eff_sel_s = (double)n_selected_s/(double)n_processed_s;
125	kukartse	1.1	double eff_sel_qcd = (double)n_selected_qcd/(double)n_processed_qcd;
126
127			double n_expected_ttbar = xsec_ttbar * eff_gen_ttbar * eff_sel_ttbar * lumi;
128			double n_expected_w = xsec_w * eff_gen_w * eff_sel_w * lumi;
129			double n_expected_z = xsec_z * eff_gen_z * eff_sel_z * lumi;
130	kukartse	1.7	double n_expected_t = xsec_t * eff_gen_t * eff_sel_t * lumi;
131			double n_expected_tW= xsec_tW* eff_gen_tW* eff_sel_tW* lumi;
132			double n_expected_s = xsec_s * eff_gen_s * eff_sel_s * lumi;
133	kukartse	1.1	double n_expected_qcd = xsec_qcd * eff_gen_qcd * eff_sel_qcd * lumi;
134
135	kukartse	1.8	double n_expected_all = n_expected_ttbar + n_expected_w + n_expected_z + n_expected_t + n_expected_tW + n_expected_s + n_expected_qcd;
136	kukartse	1.7	double n_expected_phys = n_expected_ttbar + n_expected_w + n_expected_z + n_expected_t + n_expected_tW + n_expected_s;
137	kukartse	1.12	double n_expected_wzt = n_expected_w + n_expected_z + n_expected_t + n_expected_tW + n_expected_s;
138	kukartse	1.5	//double n_expected_all = n_expected_ttbar + n_expected_z + n_expected_qcd;
139			//double n_expected_phys = n_expected_ttbar + n_expected_z;
140	kukartse	1.1
141	kukartse	1.7	//
142			// _____ mixing phys bg _______________________________________________
143			//
144			double eff_lumi_ttbar = n_processed_ttbar/xsec_ttbar/eff_gen_ttbar;
145			double eff_lumi_w = n_processed_w/xsec_w/eff_gen_w;
146			double eff_lumi_z = n_processed_z/xsec_z/eff_gen_z;
147			double eff_lumi_t = n_processed_t/xsec_t/eff_gen_t;
148			double eff_lumi_tW= n_processed_tW/xsec_tW/eff_gen_tW;
149			double eff_lumi_s = n_processed_s/xsec_s/eff_gen_s;
150			//
151	kukartse	1.12	double w_rate_corr = 133.483/103.891;
152			double z_rate_corr = 25.0594/18.5393;
153	kukartse	1.11	//
154			double n_template_w = (double)n_selected_weff_lumi_ww_rate_corr/eff_lumi_w/w_rate_corr;
155			double n_template_z = (double)n_selected_zeff_lumi_wz_rate_corr/eff_lumi_z/w_rate_corr;
156			double n_template_t = (double)n_selected_t*eff_lumi_w/eff_lumi_t/w_rate_corr;
157			double n_template_tW= (double)n_selected_tW*eff_lumi_w/eff_lumi_tW/w_rate_corr;
158			double n_template_s = (double)n_selected_s*eff_lumi_w/eff_lumi_s/w_rate_corr;
159	kukartse	1.7
160
161	kukartse	1.1	cout << endl << "Selection efficiency, ttbar events: " << eff_sel_ttbar << endl;
162	kukartse	1.4	cout << "Selected ttbar events: " << n_selected_ttbar << endl;
163	kukartse	1.1	cout << "Expected ttbar events: " << n_expected_ttbar << endl;
164	kukartse	1.13	cout << "Effective ttbar lumi: " << eff_lumi_ttbar << endl;
165	kukartse	1.7	//
166	kukartse	1.1	cout << endl << "Selection efficiency, w events: " << eff_sel_w << endl;
167	kukartse	1.4	cout << "Selected w events: " << n_selected_w << endl;
168	kukartse	1.1	cout << "Expected w events: " << n_expected_w << endl;
169	kukartse	1.7	cout << "Effective w lumi: " << eff_lumi_w << endl;
170			cout << "Events for the template: " << n_template_w << endl;
171			//
172	kukartse	1.1	cout << endl << "Selection efficiency, z events: " << eff_sel_z << endl;
173	kukartse	1.4	cout << "Selected z events: " << n_selected_z << endl;
174	kukartse	1.1	cout << "Expected z events: " << n_expected_z << endl;
175	kukartse	1.7	cout << "Effective z lumi: " << eff_lumi_z << endl;
176			cout << "Events for the template: " << n_template_z << endl;
177			//
178			cout << endl << "Selection efficiency, single top t-channel events: " << eff_sel_t << endl;
179			cout << "Selected t events: " << n_selected_t << endl;
180			cout << "Expected t events: " << n_expected_t << endl;
181			cout << "Effective t lumi: " << eff_lumi_t << endl;
182			cout << "Events for the template: " << n_template_t << endl;
183			//
184			cout << endl << "Selection efficiency, single top tW-channel events: " << eff_sel_tW << endl;
185			cout << "Selected tW events: " << n_selected_tW << endl;
186			cout << "Expected tW events: " << n_expected_tW << endl;
187			cout << "Effective tW lumi: " << eff_lumi_tW << endl;
188			cout << "Events for the template: " << n_template_tW << endl;
189			//
190			cout << endl << "Selection efficiency, single top s-channel events: " << eff_sel_s << endl;
191			cout << "Selected s events: " << n_selected_s << endl;
192			cout << "Expected s events: " << n_expected_s << endl;
193			cout << "Effective s lumi: " << eff_lumi_s << endl;
194			cout << "Events for the template: " << n_template_s << endl;
195			//
196	kukartse	1.1	cout << endl << "Selection efficiency, qcd events: " << eff_sel_qcd << endl;
197	kukartse	1.4	cout << "Selected qcd events: " << n_selected_qcd << endl;
198	kukartse	1.1	cout << "Expected qcd events: " << n_expected_qcd << endl;
199
200			cout << endl << "Total events expected: " << n_expected_all << endl;
201			cout << "Ttbar,W,Z+jets events expected: " << n_expected_phys << endl;
202	kukartse	1.12	cout << "ttbar events expected: " << n_expected_ttbar << endl;
203			cout << "phys BG events expected: " << n_expected_wzt << endl;
204			cout << "QCD events expected: " << n_expected_qcd << endl;
205	kukartse	1.1
206
207
208			}
209
210