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