MultivariateAnalysis/root/make_training_sample.C

{
<<<<<<< make_training_sample.C
  TChain ch_ttbar("ttljets");
  //TChain ch_ttbar("mva/ttljets");
  TChain ch_w("ttljets");
  //TChain ch_w("mva/ttljets");
=======
  //TChain ch_ttbar("ttljets");
  TChain ch_ttbar("ttljets");
  //TChain ch_w("ttljets");
  TChain ch_w("ttljets");
>>>>>>> 1.12
  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");
<<<<<<< make_training_sample.C
  //ch_ttbar . Add("TTJets_Fall08.root");
=======
// ch_ttbar . Add("TTJets_Fall08.root");
  ch_ttbar . Add("TTJets-madgraph_Fall08_IDEAL_jes_down_5_12Jun2009.root");
>>>>>>> 1.12
  //
  ch_w . Add("Wjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
  //ch_w . Add("Wjets_jes_up_21may2009.root");
  //ch_w . Add("Wjets_jes_down_21may2009.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_Threshold20GeV_22may2009.root");
  //ch_w . Add("WJets_Threshold5GeV_22may2009.root");
<<<<<<< make_training_sample.C
  //ch_w . Add("Wjets_Fullsim.root");
  //
  ch_z . Add("Zjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
=======
// ch_w . Add("Wjets_Fullsim.root");
 ch_w . Add("Wjets-madgraph_Winter09_IDEAL_V11_FastSim_jes_down_5_12Jun2009.root");
         //
//  ch_z . Add("Zjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
>>>>>>> 1.12
  //ch_z . Add("Zjets_jes_up_21may2009.root");
  //ch_z . Add("Zjets_jes_down_21may2009.root");
  ch_z . Add("Zjets-madgraph_Winter09_IDEAL_V11_FastSim_jes_down_5_12Jun2009.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_Summer08_IDEAL_jes_down_5_12Jun2009.root");
  //ch_qcd . Add("InclusiveMuPt15_jes_up_21may2009.root");
  //ch_qcd . Add("InclusiveMuPt15_jes_down_21may2009.root");
  /**************************
  ch_ttbar . Add("TTJets-madgraph_Fall08_IDEAL_V9_v2_14may2009.root");
  ch_w . Add("Wjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
  ch_z . Add("Zjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.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");
  ***************************/

  // selection
  #include "cuts.C"
  the_cut = cut4;
  TCut the_qcd_cut = cut4_qcd;

  // when copying, include all branches
  ch_ttbar . SetBranchStatus("*",1);
  ch_w . SetBranchStatus("*",1);
  ch_z . SetBranchStatus("*",1);
  ch_qcd . SetBranchStatus("*",1);

  // create new TTrees from the input
  TTree * t_sig = ch_ttbar . CopyTree(the_cut);
  TTree * t_wjets = ch_w . CopyTree(the_cut);
  TTree * t_zjets = ch_z . CopyTree(the_cut);
  TTree * t_t  = ch_t  . CopyTree(the_cut);
  TTree * t_tW = ch_tW . CopyTree(the_cut);
  TTree * t_s  = ch_s  . CopyTree(the_cut);
  TTree * t_qcd = ch_qcd . CopyTree(the_qcd_cut);
  //
  //TTree * _sig = t_sig->CopyTree("","",47831,392); // nominal, leave events for mock data
  //TTree * _sig = t_sig->CopyTree("","");             // systematics - take all events
  //
  //TTree * _wjets = t_wjets->CopyTree("","");         // systematics - take all events
  //TTree * _wjets = t_wjets->CopyTree("","",13177,0);
  //TTree * _zjets = t_zjets->CopyTree("","",2474,0);
  TTree * _t  = t_t ->CopyTree("","",523,0);
  TTree * _tW = t_tW->CopyTree("","",1092,0);
  TTree * _s  = t_s ->CopyTree("","",18,0);
  //
  //_____ WJets threshold 20 GeV
  //
  /****************
  TTree * _wjets = t_wjets->CopyTree("","",3592,0);
  TTree * _zjets = t_zjets->CopyTree("","",508,0);
  TTree * _t  = t_t ->CopyTree("","",145,0);
  TTree * _tW = t_tW->CopyTree("","",99,0);
  TTree * _s  = t_s ->CopyTree("","",5,0);
  *****************/
  //
  //_____ WJets scale up
  //
  /****************
  TTree * _wjets = t_wjets->CopyTree("","",1115,0);
  TTree * _zjets = t_zjets->CopyTree("","",331,0);
  TTree * _t  = t_t ->CopyTree("","",95,0);
  TTree * _tW = t_tW->CopyTree("","",64,0);
  TTree * _s  = t_s ->CopyTree("","",3,0);
  *****************/
  //
  //_____ WJets threshold 5 GeV
  //
  /****************
  TTree * _wjets = t_wjets->CopyTree("","",2014,0);
  TTree * _zjets = t_zjets->CopyTree("","",359,0);
  TTree * _t  = t_t ->CopyTree("","",103,0);
  TTree * _tW = t_tW->CopyTree("","",70,0);
  TTree * _s  = t_s ->CopyTree("","",4,0);
  *****************/
  //
  //_____ JES up
  //
  /****************
  TTree * _sig = t_sig->CopyTree("","");
//  TTree * _wjets = t_wjets->CopyTree("","",18324,0);
//   TTree * _wjets = t_wjets->CopyTree("","",15572,0); //jes up 5
   TTree * _wjets = t_wjets->CopyTree("","",14093,0); //jes up 2
//  TTree * _zjets = t_zjets->CopyTree("","",3321,0);
//  TTree * _zjets = t_zjets->CopyTree("","",2952,0);//jes up 5
  TTree * _zjets = t_zjets->CopyTree("","",2657,0);//jes up 2
  *****************/
  //
  //_____ JES down
  //
// /******************** nominal
  TTree * _sig = t_sig->CopyTree("","");
//  TTree * _wjets = t_wjets->CopyTree("","",9231,0);
  TTree * _wjets = t_wjets->CopyTree("","",11126,0);//jes down 5
//  TTree * _wjets = t_wjets->CopyTree("","",12360,0);//jes down 2
//  TTree * _zjets = t_zjets->CopyTree("","",1635,0);
  TTree * _zjets = t_zjets->CopyTree("","",2077,0);//jes down 5
//  TTree * _zjets = t_zjets->CopyTree("","",2314,0);//jes down 2
//  *********************/
  //
  //_wjets -> CopyAddresses(_zjets);
  //_wjets -> CopyEntries(_zjets);
  //_wjets -> CopyAddresses(_t);
  //_wjets -> CopyEntries(_t);
  //_wjets -> CopyAddresses(_tW);
  //_wjets -> CopyEntries(_tW);
  //_wjets -> CopyAddresses(_s);
  //_wjets -> CopyEntries(_s);
  //
  TTree * _qcd = t_qcd->CopyTree("","");


  // define new TTree names
  _sig.SetName("ttbar");
  _wjets.SetName("wjets");
  _qcd.SetName("qcd");
  //_sig->Print();
  //_wjets->Print();
  //_qcd->Print();

  // output file
  TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training_jes_down_5_15jun2009.root", "RECREATE");
  //TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-wfastsim-summer08-30may2009.root", "RECREATE");
  //TFile * out_file = new TFile("./tmva_training-summer08-25feb2009.root", "RECREATE");
  //TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-jes-up-summer08-22may2009.root", "RECREATE");
  //TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-jes-down-summer08-22may2009.root", "RECREATE");
  //TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-tauola-summer08-23may2009.root", "RECREATE");
  //TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-ISR_FSR_larger-summer08-23may2009.root", "RECREATE");
  //TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-ISR_FSR_smaller-summer08-23may2009.root", "RECREATE");
  //TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training_WJets_threshold20-summer08-23may2009.root", "RECREATE");
  //TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training_WJets_threshold5-summer08-23may2009.root", "RECREATE");
  //TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training_WJets_scaleUp-summer08-28may2009.root", "RECREATE");
  out_file->cd();
  _sig->Write();
  _wjets->Write();
  _qcd->Write();
  
  out_file->Write();
  delete out_file;
}
 
Revision:	1.13
Committed:	Thu Jun 18 16:46:55 2009 UTC (15 years, 10 months ago) by kukartse
Content type:	text/plain
Branch:	MAIN
CVS Tags:	V00-02-02, gak011410, gak010310, ejterm2010_25nov2009, V00-02-01, V00-02-00, gak112409, CMSSW_22X_branch_base, segala101609, V00-01-15
Branch point for:	CMSSW_22X_branch
Changes since 1.12:	+19 -2 lines
Log Message:	* empty log message *
#	User	Rev	Content
1	kukartse	1.1	{
2	kukartse	1.13	<<<<<<< make_training_sample.C
3			TChain ch_ttbar("ttljets");
4			//TChain ch_ttbar("mva/ttljets");
5			TChain ch_w("ttljets");
6			//TChain ch_w("mva/ttljets");
7			=======
8	kukartse	1.11	//TChain ch_ttbar("ttljets");
9	jindal	1.12	TChain ch_ttbar("ttljets");
10	kukartse	1.11	//TChain ch_w("ttljets");
11	jindal	1.12	TChain ch_w("ttljets");
12	kukartse	1.13	>>>>>>> 1.12
13	kukartse	1.7	TChain ch_z("ttljets");
14			TChain ch_tW("ttljets"); // single top
15			TChain ch_t("ttljets"); // single top
16			TChain ch_s("ttljets"); // single top
17			TChain ch_qcd("ttljets");
18	kukartse	1.9	//
19	kukartse	1.13	ch_ttbar . Add("TTJets-madgraph_Fall08_IDEAL_V9_v2_14may2009.root");
20	kukartse	1.9	//ch_ttbar . Add("TTJets_jes_up_21may2009.root");
21			//ch_ttbar . Add("TTJets_jes_down_21may2009.root");
22			//ch_ttbar . Add("Tauola_23mar2009.root");
23			//ch_ttbar . Add("TTJets_LargerISRFSR_23may2009.root");
24			//ch_ttbar . Add("TTJets_SmallerISRFSR_23may2009.root");
25	kukartse	1.13	<<<<<<< make_training_sample.C
26			//ch_ttbar . Add("TTJets_Fall08.root");
27			=======
28	jindal	1.12	// ch_ttbar . Add("TTJets_Fall08.root");
29			ch_ttbar . Add("TTJets-madgraph_Fall08_IDEAL_jes_down_5_12Jun2009.root");
30	kukartse	1.13	>>>>>>> 1.12
31	kukartse	1.9	//
32	kukartse	1.13	ch_w . Add("Wjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
33	kukartse	1.9	//ch_w . Add("Wjets_jes_up_21may2009.root");
34			//ch_w . Add("Wjets_jes_down_21may2009.root");
35	kukartse	1.10	//ch_w . Add("/uscms_data/d2/lpcljm/MVA/Summer08/May_28_2009/WJets_ScaleUp-madgraph_Winter09/Muons/all.root");
36	kukartse	1.9	//ch_w . Add("/uscms_data/d2/lpcljm/MVA/Summer08/May_28_2009/WJets_ScaleDown-madgraph_Winter09/Muons/all.root");
37			//ch_w . Add("WJets_Threshold20GeV_22may2009.root");
38			//ch_w . Add("WJets_Threshold5GeV_22may2009.root");
39	kukartse	1.13	<<<<<<< make_training_sample.C
40			//ch_w . Add("Wjets_Fullsim.root");
41			//
42			ch_z . Add("Zjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
43			=======
44	jindal	1.12	// ch_w . Add("Wjets_Fullsim.root");
45			ch_w . Add("Wjets-madgraph_Winter09_IDEAL_V11_FastSim_jes_down_5_12Jun2009.root");
46			//
47			// ch_z . Add("Zjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
48	kukartse	1.13	>>>>>>> 1.12
49	kukartse	1.9	//ch_z . Add("Zjets_jes_up_21may2009.root");
50			//ch_z . Add("Zjets_jes_down_21may2009.root");
51	jindal	1.12	ch_z . Add("Zjets-madgraph_Winter09_IDEAL_V11_FastSim_jes_down_5_12Jun2009.root");
52			//
53	kukartse	1.9	ch_t . Add("SingleTop_tChannel_Summer08_IDEAL_V9_v1.root");
54			ch_tW. Add("SingleTop_tWChannel_Summer08_IDEAL_V9_v1.root");
55			ch_s . Add("SingleTop_sChannel_Summer08_IDEAL_V9_v1.root");
56			//
57	jindal	1.12	// ch_qcd . Add("InclusiveMuPt15_Summer08_IDEAL_V9_v1_14may2009.root");
58			ch_qcd . Add("InclusiveMuPt15_Summer08_IDEAL_jes_down_5_12Jun2009.root");
59	kukartse	1.9	//ch_qcd . Add("InclusiveMuPt15_jes_up_21may2009.root");
60			//ch_qcd . Add("InclusiveMuPt15_jes_down_21may2009.root");
61			/**************************
62	kukartse	1.7	ch_ttbar . Add("TTJets-madgraph_Fall08_IDEAL_V9_v2_14may2009.root");
63			ch_w . Add("Wjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
64			ch_z . Add("Zjets-madgraph_Winter09_IDEAL_V11_FastSim_v1.root");
65			ch_t . Add("SingleTop_tChannel_Summer08_IDEAL_V9_v1.root");
66			ch_tW. Add("SingleTop_tWChannel_Summer08_IDEAL_V9_v1.root");
67			ch_s . Add("SingleTop_sChannel_Summer08_IDEAL_V9_v1.root");
68			ch_qcd . Add("InclusiveMuPt15_Summer08_IDEAL_V9_v1_14may2009.root");
69	kukartse	1.9	***************************/
70	kukartse	1.1
71			// selection
72	kukartse	1.2	#include "cuts.C"
73	kukartse	1.7	the_cut = cut4;
74			TCut the_qcd_cut = cut4_qcd;
75	kukartse	1.1
76			// when copying, include all branches
77	kukartse	1.7	ch_ttbar . SetBranchStatus("*",1);
78			ch_w . SetBranchStatus("*",1);
79			ch_z . SetBranchStatus("*",1);
80			ch_qcd . SetBranchStatus("*",1);
81	kukartse	1.1
82			// create new TTrees from the input
83	kukartse	1.7	TTree * t_sig = ch_ttbar . CopyTree(the_cut);
84			TTree * t_wjets = ch_w . CopyTree(the_cut);
85			TTree * t_zjets = ch_z . CopyTree(the_cut);
86			TTree * t_t = ch_t . CopyTree(the_cut);
87			TTree * t_tW = ch_tW . CopyTree(the_cut);
88			TTree * t_s = ch_s . CopyTree(the_cut);
89			TTree * t_qcd = ch_qcd . CopyTree(the_qcd_cut);
90			//
91	kukartse	1.11	//TTree * _sig = t_sig->CopyTree("","",47831,392); // nominal, leave events for mock data
92	jindal	1.12	//TTree * _sig = t_sig->CopyTree("",""); // systematics - take all events
93	kukartse	1.9	//
94	jindal	1.12	//TTree * _wjets = t_wjets->CopyTree("",""); // systematics - take all events
95	kukartse	1.11	//TTree * _wjets = t_wjets->CopyTree("","",13177,0);
96			//TTree * _zjets = t_zjets->CopyTree("","",2474,0);
97	jindal	1.12	TTree * _t = t_t ->CopyTree("","",523,0);
98			TTree * _tW = t_tW->CopyTree("","",1092,0);
99			TTree * _s = t_s ->CopyTree("","",18,0);
100	kukartse	1.9	//
101			//_____ WJets threshold 20 GeV
102			//
103			/****************
104			TTree * _wjets = t_wjets->CopyTree("","",3592,0);
105			TTree * _zjets = t_zjets->CopyTree("","",508,0);
106			TTree * _t = t_t ->CopyTree("","",145,0);
107			TTree * _tW = t_tW->CopyTree("","",99,0);
108			TTree * _s = t_s ->CopyTree("","",5,0);
109			*****************/
110			//
111			//_____ WJets scale up
112			//
113	kukartse	1.10	/****************
114	kukartse	1.9	TTree * _wjets = t_wjets->CopyTree("","",1115,0);
115			TTree * _zjets = t_zjets->CopyTree("","",331,0);
116			TTree * _t = t_t ->CopyTree("","",95,0);
117			TTree * _tW = t_tW->CopyTree("","",64,0);
118			TTree * _s = t_s ->CopyTree("","",3,0);
119	kukartse	1.10	*****************/
120	kukartse	1.9	//
121			//_____ WJets threshold 5 GeV
122			//
123			/****************
124			TTree * _wjets = t_wjets->CopyTree("","",2014,0);
125			TTree * _zjets = t_zjets->CopyTree("","",359,0);
126			TTree * _t = t_t ->CopyTree("","",103,0);
127			TTree * _tW = t_tW->CopyTree("","",70,0);
128			TTree * _s = t_s ->CopyTree("","",4,0);
129			*****************/
130			//
131			//_____ JES up
132			//
133			/****************
134			TTree * _sig = t_sig->CopyTree("","");
135	jindal	1.12	// TTree * _wjets = t_wjets->CopyTree("","",18324,0);
136			// TTree * _wjets = t_wjets->CopyTree("","",15572,0); //jes up 5
137			TTree * _wjets = t_wjets->CopyTree("","",14093,0); //jes up 2
138			// TTree * _zjets = t_zjets->CopyTree("","",3321,0);
139			// TTree * _zjets = t_zjets->CopyTree("","",2952,0);//jes up 5
140			TTree * _zjets = t_zjets->CopyTree("","",2657,0);//jes up 2
141	kukartse	1.9	*****************/
142			//
143			//_____ JES down
144			//
145	jindal	1.12	// /******************** nominal
146	kukartse	1.9	TTree * _sig = t_sig->CopyTree("","");
147	jindal	1.12	// TTree * _wjets = t_wjets->CopyTree("","",9231,0);
148			TTree * _wjets = t_wjets->CopyTree("","",11126,0);//jes down 5
149			// TTree * _wjets = t_wjets->CopyTree("","",12360,0);//jes down 2
150			// TTree * _zjets = t_zjets->CopyTree("","",1635,0);
151			TTree * _zjets = t_zjets->CopyTree("","",2077,0);//jes down 5
152			// TTree * _zjets = t_zjets->CopyTree("","",2314,0);//jes down 2
153			// *********************/
154	kukartse	1.9	//
155	kukartse	1.11	//_wjets -> CopyAddresses(_zjets);
156			//_wjets -> CopyEntries(_zjets);
157			//_wjets -> CopyAddresses(_t);
158			//_wjets -> CopyEntries(_t);
159			//_wjets -> CopyAddresses(_tW);
160			//_wjets -> CopyEntries(_tW);
161			//_wjets -> CopyAddresses(_s);
162			//_wjets -> CopyEntries(_s);
163			//
164	kukartse	1.7	TTree * _qcd = t_qcd->CopyTree("","");
165	kukartse	1.1
166
167			// define new TTree names
168			_sig.SetName("ttbar");
169			_wjets.SetName("wjets");
170	kukartse	1.3	_qcd.SetName("qcd");
171	kukartse	1.1	//_sig->Print();
172			//_wjets->Print();
173			//_qcd->Print();
174
175			// output file
176	jindal	1.12	TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training_jes_down_5_15jun2009.root", "RECREATE");
177	kukartse	1.11	//TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-wfastsim-summer08-30may2009.root", "RECREATE");
178	kukartse	1.3	//TFile * out_file = new TFile("./tmva_training-summer08-25feb2009.root", "RECREATE");
179	kukartse	1.9	//TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-jes-up-summer08-22may2009.root", "RECREATE");
180			//TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-jes-down-summer08-22may2009.root", "RECREATE");
181			//TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-tauola-summer08-23may2009.root", "RECREATE");
182			//TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-ISR_FSR_larger-summer08-23may2009.root", "RECREATE");
183			//TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training-ISR_FSR_smaller-summer08-23may2009.root", "RECREATE");
184			//TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training_WJets_threshold20-summer08-23may2009.root", "RECREATE");
185			//TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training_WJets_threshold5-summer08-23may2009.root", "RECREATE");
186	kukartse	1.10	//TFile * out_file = new TFile("/uscms_data/d1/lpcljm/MVA/Summer08/training/muon_jets_training_WJets_scaleUp-summer08-28may2009.root", "RECREATE");
187	kukartse	1.1	out_file->cd();
188			_sig->Write();
189			_wjets->Write();
190	kukartse	1.3	_qcd->Write();
191	kukartse	1.1
192			out_file->Write();
193			delete out_file;
194			}
195