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 |
|
|
}
|