MultivariateAnalysis/root/TreeReader_New.C

#include <fstream>
#include <iomanip>
#include <iostream>
#include <string.h>
#include <fstream>

#include <TFile.h>
#include <TH1.h>
#include <TTree.h>
#include <TCanvas.h>
#include <TF1.h>
#include <TStyle.h>
#include <TCut.h>
#include "cuts.C"

int TreeReader(char *infile = "/uscms/home/kukarzev/work/disk_ljmet/kukarzev/tmp/ljmet_tmva_InclusiveMuPt15-12feb2009.root",char *outfile = "test.root")
{
  TFile *fin = TFile::Open(infile);
  TFile *outf = new TFile(outfile,"RECREATE");


   TH1D * hMt_nominal  = new TH1D( "hMt_nominal", "Mt (nominal selection)", 100, 0., 90.); 
   TH1D * hMt_reversed = new TH1D( "hMt_reversed", "Mt (reversed selection)", 100, 0., 90.); 
  
 
   TH1D * hht_nominal  = new TH1D( "hht_nominal", "ht  (nominal) ", 500, 0., 2000.); 
   TH1D * hht_reversed = new TH1D( "hht_reversed", "ht  (reversed) ", 500, 0., 2000.); 
   

   TH1D * hM123inv_nominal  = new TH1D( "hM123inv_nominal", "M123inv (nominal selection)", 500, 0., 4500.); 
   TH1D * hM123inv_reversed = new TH1D( "hM123inv_reversed", "M123inv (reversed selection)", 500, 0., 4500.); 
  

   TH1D * hgetHt3_nominal  = new TH1D( "hgetHt3_nominal", "Ht3 (nominal selection)", 100, 0., 900.); 
   TH1D * hgetHt3_reversed = new TH1D( "hgetHt3_reversed", "Ht3 (reversed selection)", 100, 0., 900.); 
 

   TH1D * hktMinPrime_nominal  = new TH1D( "hktMinPrime_nominal", "ktMinPrime (nominal selection)", 100, 0., 25.); 
   TH1D * hktMinPrime_reversed = new TH1D( "hktMinPrime_reversed", "ktMinPrime (reversed selection)", 100, 0., 25.); 

   
   TH1D * hht2p_nominal  = new TH1D( "hht2p_nominal", "ht2p (nominal selection)", 100, 0., 4.); 
   TH1D * hht2p_reversed = new TH1D( "hht2p_reversed", "ht2p (reversed selection)", 100, 0., 4.); 


   TH1D * hcentrality_nominal  = new TH1D( "hcentrality_nominal", "centrality (nominal selection)", 100, 0., 1.); 
   TH1D * hcentrality_reversed = new TH1D( "hcentrality_reversed", "centrality (reversed selection)", 100, 0., 1.); 


   TH1D * hsphericity_nominal  = new TH1D( "hsphericity_nominal", "sphericity (nominal selection)", 100, 0., 1.); 
   TH1D * hsphericity_reversed = new TH1D( "hsphericity_reversed", "sphericity (reversed selection)", 100, 0., 1.); 


   TH1D * hhtPlusLepton_nominal  = new TH1D( "hhtPlusLepton_nominal", "htPlusLepton (nominal selection)", 500, 0., 3000.); 
   TH1D * hhtPlusLepton_reversed = new TH1D( "hhtPlusLepton_reversed", "htPlusLepton (reversed selection)", 500, 0., 3000.); 


   TH1D * hdPhiLMet_nominal  = new TH1D( "hdPhiLMet_nominal", "dPhiLMet (nominal selection)", 100, 0., 3.5); 
   TH1D * hdPhiLMet_reversed = new TH1D( "hdPhiLMet_reversed", "dPhiLMet (reversed selection)", 100, 0., 3.5); 


   TH1D * hDphiJMET_nominal  = new TH1D( "hDphiJMET_nominal", "DphiJMET (nominal selection)", 100, 0., 3.5); 
   TH1D * hDphiJMET_reversed = new TH1D( "hDphiJMET_reversed", "DphiJMET (reversed selection)", 100, 0., 3.5); 


   TH1D * haplanarity_nominal  = new TH1D( "haplanarity_nominal", "aplanarity (nominal selection)", 100, 0., 0.5); 
   TH1D * haplanarity_reversed = new TH1D( "haplanarity_reversed", "aplanarity (reversed selection)", 100, 0., 0.5); 


   TH1D * hW_MT_nominal  = new TH1D( "hW_MT_nominal", "W_MT (nominal selection)", 500, 0., 5000); 
   TH1D * hW_MT_reversed = new TH1D( "hW_MT_reversed", "W_MT (reversed selection)", 500, 0., 5000); 


   TH1D * hmetHtPlusLepton_nominal  = new TH1D( "hmetHtPlusLepton_nominal", "metHtPlusLepton (nominal selection)", 500, 0., 4000.); 
   TH1D * hmetHtPlusLepton_reversed = new TH1D( "hmetHtPlusLepton_reversed", "metHtPlusLepton (reversed selection)", 500, 0., 4000.);


   TTree * t1 = (TTree*)fin->Get("ttljets");

   t1->Draw("ht>> hht_nominal",semilep_cut && reversible_cut);
   t1->Draw("getMt>> hMt_nominal",semilep_cut && reversible_cut);
   t1->Draw("getM123inv>> hM123inv_nominal",semilep_cut && reversible_cut); 
   t1->Draw("getHt3>> hgetHt3_nominal",semilep_cut && reversible_cut);
   t1->Draw("ktMinPrime>> hktMinPrime_nominal",semilep_cut && reversible_cut);
   t1->Draw("ht2p>> hht2p_nominal",semilep_cut && reversible_cut);
   t1->Draw("aplanarity>> haplanarity_nominal",semilep_cut && reversible_cut);
   t1->Draw("centrality>> hcentrality_nominal",semilep_cut && reversible_cut);
   t1->Draw("sphericity>> hsphericity_nominal",semilep_cut && reversible_cut);
   t1->Draw("dPhiLMet>> hdPhiLMet_nominal",semilep_cut && reversible_cut);
   t1->Draw("DphiJMET>> hDphiJMET_nominal",semilep_cut && reversible_cut);
   t1->Draw("W_MT>> hW_MT_nominal",semilep_cut && reversible_cut);
   t1->Draw("htPlusLepton>> hhtPlusLepton_nominal",semilep_cut && reversible_cut);
   t1->Draw("metHtPlusLepton>> hmetHtPlusLepton_nominal",semilep_cut && reversible_cut);
 

   t1->Draw("ht>> hht_reversed",semilep_cut && !reversible_cut);
   t1->Draw("getMt>> hMt_reversed",semilep_cut && !reversible_cut);
   t1->Draw("getM123inv>> hM123inv_reversed",semilep_cut && !reversible_cut); 
   t1->Draw("getHt3>> hgetHt3_reversed",semilep_cut && !reversible_cut);
   t1->Draw("ktMinPrime>> hktMinPrime_reversed",semilep_cut && !reversible_cut);
   t1->Draw("ht2p>> hht2p_reversed",semilep_cut && !reversible_cut);
   t1->Draw("aplanarity>> haplanarity_reversed",semilep_cut && !reversible_cut);
   t1->Draw("centrality>> hcentrality_reversed",semilep_cut && !reversible_cut);
   t1->Draw("sphericity>> hsphericity_reversed",semilep_cut && !reversible_cut);
   t1->Draw("dPhiLMet>> hdPhiLMet_reversed",semilep_cut && !reversible_cut);
   t1->Draw("DphiJMET>> hDphiJMET_reversed",semilep_cut && !reversible_cut);
   t1->Draw("W_MT>> hW_MT_reversed",semilep_cut && !reversible_cut);
   t1->Draw("htPlusLepton>> hhtPlusLepton_reversed",semilep_cut && !reversible_cut);
   t1->Draw("metHtPlusLepton>> hmetHtPlusLepton_reversed",semilep_cut && !reversible_cut);

  
  outf->Write();
  outf->Close();


}
Revision:	1.1
Committed:	Fri Feb 13 00:53:03 2009 UTC (16 years, 2 months ago) by kukartse
Content type:	text/plain
Branch:	MAIN
CVS Tags:	V00-03-01, ZMorph_BASE_20100408, gak040610_morphing, V00-02-02, gak011410, gak010310, ejterm2010_25nov2009, V00-02-01, V00-02-00, gak112409, CMSSW_22X_branch_base, segala101609, V00-01-15, V00-01-14, V00-01-13, V00-01-12, V00-01-11, V00-01-10, gak031009, gak030509, gak022309, gak021209, HEAD
Branch point for:	ZMorph-V00-03-01, CMSSW_22X_branch
Log Message:	* empty log message *
#	User	Rev	Content
1	kukartse	1.1	#include <fstream>
2			#include <iomanip>
3			#include <iostream>
4			#include <string.h>
5			#include <fstream>
6
7			#include <TFile.h>
8			#include <TH1.h>
9			#include <TTree.h>
10			#include <TCanvas.h>
11			#include <TF1.h>
12			#include <TStyle.h>
13			#include <TCut.h>
14			#include "cuts.C"
15
16			int TreeReader(char infile = "/uscms/home/kukarzev/work/disk_ljmet/kukarzev/tmp/ljmet_tmva_InclusiveMuPt15-12feb2009.root",char outfile = "test.root")
17			{
18			TFile *fin = TFile::Open(infile);
19			TFile *outf = new TFile(outfile,"RECREATE");
20
21
22
23			TH1D * hMt_nominal = new TH1D( "hMt_nominal", "Mt (nominal selection)", 100, 0., 90.);
24			TH1D * hMt_reversed = new TH1D( "hMt_reversed", "Mt (reversed selection)", 100, 0., 90.);
25
26
27
28			TH1D * hht_nominal = new TH1D( "hht_nominal", "ht (nominal) ", 500, 0., 2000.);
29			TH1D * hht_reversed = new TH1D( "hht_reversed", "ht (reversed) ", 500, 0., 2000.);
30
31
32			TH1D * hM123inv_nominal = new TH1D( "hM123inv_nominal", "M123inv (nominal selection)", 500, 0., 4500.);
33			TH1D * hM123inv_reversed = new TH1D( "hM123inv_reversed", "M123inv (reversed selection)", 500, 0., 4500.);
34
35
36			TH1D * hgetHt3_nominal = new TH1D( "hgetHt3_nominal", "Ht3 (nominal selection)", 100, 0., 900.);
37			TH1D * hgetHt3_reversed = new TH1D( "hgetHt3_reversed", "Ht3 (reversed selection)", 100, 0., 900.);
38
39
40
41			TH1D * hktMinPrime_nominal = new TH1D( "hktMinPrime_nominal", "ktMinPrime (nominal selection)", 100, 0., 25.);
42			TH1D * hktMinPrime_reversed = new TH1D( "hktMinPrime_reversed", "ktMinPrime (reversed selection)", 100, 0., 25.);
43
44
45
46			TH1D * hht2p_nominal = new TH1D( "hht2p_nominal", "ht2p (nominal selection)", 100, 0., 4.);
47			TH1D * hht2p_reversed = new TH1D( "hht2p_reversed", "ht2p (reversed selection)", 100, 0., 4.);
48
49
50			TH1D * hcentrality_nominal = new TH1D( "hcentrality_nominal", "centrality (nominal selection)", 100, 0., 1.);
51			TH1D * hcentrality_reversed = new TH1D( "hcentrality_reversed", "centrality (reversed selection)", 100, 0., 1.);
52
53
54			TH1D * hsphericity_nominal = new TH1D( "hsphericity_nominal", "sphericity (nominal selection)", 100, 0., 1.);
55			TH1D * hsphericity_reversed = new TH1D( "hsphericity_reversed", "sphericity (reversed selection)", 100, 0., 1.);
56
57
58			TH1D * hhtPlusLepton_nominal = new TH1D( "hhtPlusLepton_nominal", "htPlusLepton (nominal selection)", 500, 0., 3000.);
59			TH1D * hhtPlusLepton_reversed = new TH1D( "hhtPlusLepton_reversed", "htPlusLepton (reversed selection)", 500, 0., 3000.);
60
61
62			TH1D * hdPhiLMet_nominal = new TH1D( "hdPhiLMet_nominal", "dPhiLMet (nominal selection)", 100, 0., 3.5);
63			TH1D * hdPhiLMet_reversed = new TH1D( "hdPhiLMet_reversed", "dPhiLMet (reversed selection)", 100, 0., 3.5);
64
65
66			TH1D * hDphiJMET_nominal = new TH1D( "hDphiJMET_nominal", "DphiJMET (nominal selection)", 100, 0., 3.5);
67			TH1D * hDphiJMET_reversed = new TH1D( "hDphiJMET_reversed", "DphiJMET (reversed selection)", 100, 0., 3.5);
68
69
70
71			TH1D * haplanarity_nominal = new TH1D( "haplanarity_nominal", "aplanarity (nominal selection)", 100, 0., 0.5);
72			TH1D * haplanarity_reversed = new TH1D( "haplanarity_reversed", "aplanarity (reversed selection)", 100, 0., 0.5);
73
74
75			TH1D * hW_MT_nominal = new TH1D( "hW_MT_nominal", "W_MT (nominal selection)", 500, 0., 5000);
76			TH1D * hW_MT_reversed = new TH1D( "hW_MT_reversed", "W_MT (reversed selection)", 500, 0., 5000);
77
78
79			TH1D * hmetHtPlusLepton_nominal = new TH1D( "hmetHtPlusLepton_nominal", "metHtPlusLepton (nominal selection)", 500, 0., 4000.);
80			TH1D * hmetHtPlusLepton_reversed = new TH1D( "hmetHtPlusLepton_reversed", "metHtPlusLepton (reversed selection)", 500, 0., 4000.);
81
82
83
84
85			TTree * t1 = (TTree*)fin->Get("ttljets");
86
87			t1->Draw("ht>> hht_nominal",semilep_cut && reversible_cut);
88			t1->Draw("getMt>> hMt_nominal",semilep_cut && reversible_cut);
89			t1->Draw("getM123inv>> hM123inv_nominal",semilep_cut && reversible_cut);
90			t1->Draw("getHt3>> hgetHt3_nominal",semilep_cut && reversible_cut);
91			t1->Draw("ktMinPrime>> hktMinPrime_nominal",semilep_cut && reversible_cut);
92			t1->Draw("ht2p>> hht2p_nominal",semilep_cut && reversible_cut);
93			t1->Draw("aplanarity>> haplanarity_nominal",semilep_cut && reversible_cut);
94			t1->Draw("centrality>> hcentrality_nominal",semilep_cut && reversible_cut);
95			t1->Draw("sphericity>> hsphericity_nominal",semilep_cut && reversible_cut);
96			t1->Draw("dPhiLMet>> hdPhiLMet_nominal",semilep_cut && reversible_cut);
97			t1->Draw("DphiJMET>> hDphiJMET_nominal",semilep_cut && reversible_cut);
98			t1->Draw("W_MT>> hW_MT_nominal",semilep_cut && reversible_cut);
99			t1->Draw("htPlusLepton>> hhtPlusLepton_nominal",semilep_cut && reversible_cut);
100			t1->Draw("metHtPlusLepton>> hmetHtPlusLepton_nominal",semilep_cut && reversible_cut);
101
102
103
104			t1->Draw("ht>> hht_reversed",semilep_cut && !reversible_cut);
105			t1->Draw("getMt>> hMt_reversed",semilep_cut && !reversible_cut);
106			t1->Draw("getM123inv>> hM123inv_reversed",semilep_cut && !reversible_cut);
107			t1->Draw("getHt3>> hgetHt3_reversed",semilep_cut && !reversible_cut);
108			t1->Draw("ktMinPrime>> hktMinPrime_reversed",semilep_cut && !reversible_cut);
109			t1->Draw("ht2p>> hht2p_reversed",semilep_cut && !reversible_cut);
110			t1->Draw("aplanarity>> haplanarity_reversed",semilep_cut && !reversible_cut);
111			t1->Draw("centrality>> hcentrality_reversed",semilep_cut && !reversible_cut);
112			t1->Draw("sphericity>> hsphericity_reversed",semilep_cut && !reversible_cut);
113			t1->Draw("dPhiLMet>> hdPhiLMet_reversed",semilep_cut && !reversible_cut);
114			t1->Draw("DphiJMET>> hDphiJMET_reversed",semilep_cut && !reversible_cut);
115			t1->Draw("W_MT>> hW_MT_reversed",semilep_cut && !reversible_cut);
116			t1->Draw("htPlusLepton>> hhtPlusLepton_reversed",semilep_cut && !reversible_cut);
117			t1->Draw("metHtPlusLepton>> hmetHtPlusLepton_reversed",semilep_cut && !reversible_cut);
118
119
120			outf->Write();
121			outf->Close();
122
123
124			}