ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/algomez/FourTop/TreeAnalyzer/test/runAnalysis.C
Revision: 1.3
Committed: Thu Feb 9 18:35:37 2012 UTC (13 years, 3 months ago) by algomez
Content type: text/plain
Branch: MAIN
Changes since 1.2: +127 -109 lines
Log Message:
new version with JEC and BDT incorporated

File Contents

# User Rev Content
1 algomez 1.1
2 algomez 1.3 void runAnalysis(TString sample="all", bool JEC=true, bool down=true, bool BDT=true, bool getLogs=false)
3 algomez 1.1 {
4     TString desdir = "/uscms/home/algomez/work/CMSSW_4_2_4/src/Yumiceva/TreeAnalyzer/test/";
5     TProof *p = TProof::Open("lite://", desdir ,desdir);
6    
7     //p->AddDynamicPath("");
8     p->Exec("gSystem->Load(\"/uscms/home/algomez/work/CMSSW_4_2_4/lib/slc5_amd64_gcc434/libYumicevaTop7TeV.so\")");
9     p->Exec("gSystem->Load(\"/uscms/home/algomez/work/CMSSW_4_2_4/lib/slc5_amd64_gcc434/libCondFormatsJetMETObjects.so\")");
10     p->AddIncludePath("/uscms/home/algomez/work/CMSSW_4_2_4/src/");
11    
12     p->Archive(" ",desdir);
13    
14     //p->AddInput(new TNamed("PROOF_OUTPUTFILE_LOCATION", "LOCAL"));
15 algomez 1.2
16 algomez 1.1 if (sample=="MC"||sample=="4Top_1000"||sample=="all")
17     {
18 algomez 1.2 TDSet *mc_4Top_1000 = new TDSet("top","*","/PATNtupleMaker");
19 algomez 1.3 mc_4Top_1000->Add("/uscms_data/d3/algomez/files/fourtop/Ntuples/4Top_1000_new_Tuple-PATskim.root");
20     if (!JEC) {
21     if (!BDT) mc_4Top_1000->Process("Analyzer.C+","sample=4Top1000");
22     else mc_4Top_1000->Process("Analyzer.C+","BDT sample=4Top1000");
23     } else {
24     if (!down) {
25     if (!BDT) mc_4Top_1000->Process("Analyzer.C+","JECUP sample=4Top1000");
26     else mc_4Top_1000->Process("Analyzer.C+","JECUP BDT sample=4Top1000");
27     } else {
28     if (!BDT) mc_4Top_1000->Process("Analyzer.C+","JECDOWN sample=4Top1000");
29     else mc_4Top_1000->Process("Analyzer.C+","JECDOWN BDT sample=4Top1000");
30     }
31     }
32 algomez 1.1 }
33    
34     if (sample=="MC"||sample=="4Top_500"||sample=="all")
35     {
36 algomez 1.2 TDSet *mc_4Top_500 = new TDSet("top","*","/PATNtupleMaker");
37 algomez 1.3 mc_4Top_500->Add("/uscms_data/d3/algomez/files/fourtop/Ntuples/4Top_500_new_Tuple-PATskim.root");
38     if (!JEC) {
39     if (!BDT) mc_4Top_500->Process("Analyzer.C+","sample=4Top500");
40     else mc_4Top_500->Process("Analyzer.C+","BDT sample=4Top500");
41     } else {
42     if (!down) {
43     if (!BDT) mc_4Top_500->Process("Analyzer.C+","JECUP sample=4Top500");
44     else mc_4Top_500->Process("Analyzer.C+","JECUP BDT sample=4Top500");
45     } else {
46     if (!BDT) mc_4Top_500->Process("Analyzer.C+","JECDOWN sample=4Top500");
47     else mc_4Top_500->Process("Analyzer.C+","JECDOWN BDT sample=4Top500");
48     }
49     }
50 algomez 1.1 }
51    
52 algomez 1.3 if (sample=="MC"||sample=="ttbar"||sample=="all"||sample=="4Top_1000"||sample=="4Top_500")
53 algomez 1.1 {
54     TDSet *mc_ttbar = new TDSet("top","*","/PATNtupleMaker");
55 algomez 1.3 mc_ttbar->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/TTbar_Mu.root");
56     if (!JEC) {
57     if (!BDT) mc_ttbar->Process("Analyzer.C+","sample=ttbar");
58     else mc_ttbar->Process("Analyzer.C+","BDT sample=ttbar");
59     } else {
60     if (!down) {
61     if (!BDT) mc_ttbar->Process("Analyzer.C+","JECUP sample=ttbar");
62     else mc_ttbar->Process("Analyzer.C+","JECUP BDT sample=ttbar");
63     } else {
64     if (!BDT) mc_ttbar->Process("Analyzer.C+","JECDOWN sample=ttbar");
65     else mc_ttbar->Process("Analyzer.C+","JECDOWN BDT sample=ttbar");
66     }
67     }
68 algomez 1.1 }
69    
70 algomez 1.3 if (!JEC){
71 algomez 1.1 if (sample=="MC"||sample=="WJets"||sample=="all")
72     {
73     TDSet *mc_Wjets = new TDSet("top","*","/PATNtupleMaker");
74 algomez 1.3 mc_Wjets->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/WJets_Mu.root");
75     if (!BDT) mc_Wjets->Process("Analyzer.C+","sample=WJets");
76     else mc_ttbar->Process("Analyzer.C+","BDT sample=WJets");
77 algomez 1.1 }
78 algomez 1.3
79 algomez 1.1 if (sample=="MC"||sample=="QCD"||sample=="all")
80     {
81     TDSet *mc_QCD = new TDSet("top","*","/PATNtupleMaker");
82 algomez 1.3 mc_QCD->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/QCD_Mu.root");
83     if (!BDT) mc_QCD->Process("Analyzer.C+","sample=QCD");
84     else mc_QCD->Process("Analyzer.C+","BDT sample=QCD");
85 algomez 1.1 }
86 algomez 1.3
87 algomez 1.1 if (sample=="MC"||sample=="STsch"||sample=="all")
88     {
89     TDSet *mc_STsch = new TDSet("top","*","/PATNtupleMaker");
90 algomez 1.3 mc_STsch->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STsch_Mu.root");
91     if (!BDT) mc_STsch->Process("Analyzer.C+","sample=STsch");
92     else mc_STsch->Process("Analyzer.C+","BDT sample=STsch");
93 algomez 1.1 }
94 algomez 1.3
95 algomez 1.1 if (sample=="MC"||sample=="STtch"||sample=="all")
96     {
97     TDSet *mc_STtch = new TDSet("top","*","/PATNtupleMaker");
98 algomez 1.3 mc_STtch->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STtch_Mu.root");
99     if (!BDT) mc_STtch->Process("Analyzer.C+","sample=STtch");
100     else mc_STtch->Process("Analyzer.C+","BDT sample=STtch");
101 algomez 1.1 }
102 algomez 1.3
103 algomez 1.1 if (sample=="MC"||sample=="STtWch"||sample=="all")
104     {
105     TDSet *mc_STtWch = new TDSet("top","*","/PATNtupleMaker");
106 algomez 1.3 mc_STtWch->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STtWch_Mu.root");
107     if (!BDT) mc_STtWch->Process("Analyzer.C+","sample=STtWch");
108     else mc_STtWch->Process("Analyzer.C+","BDT sample=STtWch");
109 algomez 1.1 }
110 algomez 1.3
111 algomez 1.1 if (sample=="MC"||sample=="STsch_bar"||sample=="all")
112     {
113     TDSet *mc_STsch_bar = new TDSet("top","*","/PATNtupleMaker");
114 algomez 1.3 mc_STsch_bar->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STsch_bar_Mu.root");
115     if (!BDT) mc_STsch_bar->Process("Analyzer.C+","sample=STsch_bar");
116     else mc_STsch_bar->Process("Analyzer.C+","BDT sample=STsch_bar");
117 algomez 1.1 }
118 algomez 1.3
119 algomez 1.1 if (sample=="MC"||sample=="STtch_bar"||sample=="all")
120     {
121     TDSet *mc_STtch_bar = new TDSet("top","*","/PATNtupleMaker");
122 algomez 1.3 mc_STtch_bar->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STtch_bar_Mu.root");
123     if (!BDT) mc_STtch_bar->Process("Analyzer.C+","sample=STtch_bar");
124     else mc_STtch_bar->Process("Analyzer.C+","BDT sample=STtch_bar");
125 algomez 1.1 }
126 algomez 1.3
127 algomez 1.1 if (sample=="MC"||sample=="STtWch_bar"||sample=="all")
128     {
129     TDSet *mc_STtWch_bar = new TDSet("top","*","/PATNtupleMaker");
130 algomez 1.3 mc_STtWch_bar->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/STtWch_bar_Mu.root");
131     if (!BDT) mc_STtWch_bar->Process("Analyzer.C+","sample=STtWch_bar");
132     else mc_STtWch_bar->Process("Analyzer.C+","BDT sample=STtWch_bar");
133 algomez 1.1 }
134 algomez 1.3
135 algomez 1.1 if (sample=="MC"||sample=="WW"||sample=="all")
136     {
137     TDSet *mc_WW = new TDSet("top","*","/PATNtupleMaker");
138 algomez 1.3 mc_WW->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/WW_Mu.root");
139     if (!BDT) mc_WW->Process("Analyzer.C+","sample=WW");
140     else mc_WW->Process("Analyzer.C+","BDT sample=WW");
141 algomez 1.1 }
142 algomez 1.3
143 algomez 1.1 if (sample=="MC"||sample=="WZ"||sample=="all")
144     {
145     TDSet *mc_WZ = new TDSet("top","*","/PATNtupleMaker");
146 algomez 1.3 mc_WZ->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/WZ_Mu.root");
147     if (!BDT) mc_WZ->Process("Analyzer.C+","sample=WZ");
148     else mc_WZ->Process("Analyzer.C+","BDT sample=WZ");
149 algomez 1.1 }
150 algomez 1.3
151 algomez 1.1 if (sample=="MC"||sample=="ZJets"||sample=="all")
152     {
153     TDSet *mc_ZJets = new TDSet("top","*","/PATNtupleMaker");
154 algomez 1.3 mc_ZJets->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/MC/v9_14/ZJets_Mu.root");
155     if (!BDT) mc_ZJets->Process("Analyzer.C+","sample=ZJets");
156     else mc_ZJets->Process("Analyzer.C+","BDT sample=ZJets");
157 algomez 1.2 }
158 algomez 1.3 }
159 algomez 1.1 if (sample=="data"||sample=="all")
160     {
161     TDSet *data = new TDSet("top","*","/PATNtupleMaker");
162 algomez 1.3 data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_May10ReReco_Mu.root");
163     data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v4_Run_165088_166967_Mu.root");
164     data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v4_Run_167039_167913_Mu.root");
165     //data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Aug05ReReco_Mu.root");
166     //data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v6_Run_172620_173198_Mu.root");
167     //data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v6_Run_173236_173692_Mu.root");
168     //data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Run2011B-PromptReco-v1_175860_178380_Mu.root");
169     //data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Run2011B-PromptReco-v1_178420_180252_Mu.root");
170     data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Aug05ReReco_Mu.root");
171     data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v6_Run_172620_173198_Mu.root");
172     data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_PromptReco-v6_Run_173236_173692_Mu.root");
173     data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Run2011B-PromptReco-v1_175860_177452_Mu.root");
174     data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Run2011B-PromptReco-v1_177718_178380_Mu.root");
175     data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/SingleMu_Run2011B-PromptReco-v1_178420_180252_Mu.root ");
176    
177     if (!BDT) data->Process("Analyzer.C+","sample=data");
178     else data->Process("Analyzer.C+","BDT sample=data");
179 algomez 1.1 // get log files
180     if (getLogs)
181     {
182     logList = p->GetManager()->GetSessionLogs()->GetListOfLogs();
183     for (int i=1; i< logList->GetSize(); ++i)
184     {
185     logElem = ( TProofLogElem* ) logList->At( i );
186     macro = logElem->GetMacro();
187     macro->SaveSource("data_muons_"+TString(Form("%i",i))+".stdout");
188     }
189     }
190     }
191 algomez 1.3
192     if (sample=="dataMissing"||sample=="all")
193     {
194     TDSet *data = new TDSet("top","*","/PATNtupleMaker");
195     data->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011Muon/42x_v9_12/missing.root ");
196    
197     if (!BDT) data->Process("Analyzer.C+","sample=dataMissing");
198     else data->Process("Analyzer.C+","BDT sample=dataMissing");
199     // get log files
200    
201     if (getLogs)
202     {
203     logList = p->GetManager()->GetSessionLogs()->GetListOfLogs();
204     for (int i=1; i< logList->GetSize(); ++i)
205     {
206     logElem = ( TProofLogElem* ) logList->At( i );
207     macro = logElem->GetMacro();
208     macro->SaveSource("data_missing_"+TString(Form("%i",i))+".stdout");
209     }
210     }
211     }
212    
213 algomez 1.1 if (sample=="dataQCD2"||sample=="all")
214     {
215     TDSet *dataQCD2 = new TDSet("top","*","/PATNtupleMaker");
216     dataQCD2->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011/42x/Run2011_May10ReReco_Mu.root");
217     dataQCD2->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011/42x/Run2011_PromptReco1_Mu.root");
218     dataQCD2->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011/42x/Run2011_PromptReco2_Mu.root");
219     dataQCD2->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011/42x/SingleMu_Aug05ReReco_Mu.root");
220     dataQCD2->Add("/uscms_data/d3/ttmuj/Documents/NtupleMaker/Data/2011/42x/SingleMu_PromptReco-v6_Mu.root");
221 algomez 1.3 if (!BDT) dataQCD2->Process("Analyzer.C+","QCD2 sample=data");
222     else dataQCD2->Process("Analyzer.C+","BDT QCD2 sample=data");
223 algomez 1.1 }
224    
225     }