00001
00002
00003
00004
00005
00006
00007
00008
00009 #include "AEMPlots.h"
00010 #include <iostream>
00011
00012 static int fNbins=24*(31+28+31+30+31+30+31+31+30+31+30);
00013
00014
00015
00016 AEMPlots::AEMPlots()
00017 {
00018 double fNbinsd = static_cast<double>(fNbins);
00019 fAEMWeekSpills = new TH1D("fAEMWeekSpills",
00020 "Spills per hour for the week",
00021 168, 0., 168.);
00022 fAEMWeekSpillTime = new TH1D("fAEMWeekSpillTime",
00023 "Spill-seconds per hour for the week",
00024 168, 0., 168.);
00025 fAEMWeekEvents = new TH1D("fAEMWeekEvents",
00026 "Events per hour for the week",
00027 168, 0., 168.);
00028 fAEMWeekRate = new TH1D("fAEMWeekRate",
00029 "Event rate for the week",
00030 168, 0., 168.);
00031 fAEMTotalSpills = new TH1D("fAEMTotalSpills",
00032 "Spills per hour for the run",
00033 fNbins, 0., fNbinsd);
00034 fAEMTotalSpillTime = new TH1D("fAEMTotalSpillTime",
00035 "Spill-seconds per hour for the run",
00036 fNbins, 0., fNbinsd);
00037 fAEMTotalEvents = new TH1D("fAEMTotalEvents",
00038 "Events per hour for the run",
00039 fNbins, 0., fNbinsd);
00040 fAEMTotalRate = new TH1D("fAEMTotalRate",
00041 "Event rate for the run",
00042 fNbins, 0., fNbinsd);
00043 fAEMWeekData = new TH2D("fAEMWeekData",
00044 "Data this week by target and momentum",
00045 2, 0., 2., 2, 0., 2.);
00046 fAEMTotalData = new TH2D("fAEMTotalData",
00047 "Integrated data by target and momentum",
00048 2, 0., 2., 2, 0., 2.);
00049 fAEMRateDist = new TH1D("fAEMRateDist",
00050 "DAQ event rate for runs this week",
00051 100, 0., 99.);
00052
00053 }
00054
00055
00056
00057 AEMPlots::~AEMPlots()
00058 {
00059 delete fAEMWeekSpills;
00060 delete fAEMWeekSpillTime;
00061 delete fAEMWeekEvents;
00062 delete fAEMWeekRate;
00063 delete fAEMTotalSpills;
00064 delete fAEMTotalSpillTime;
00065 delete fAEMTotalEvents;
00066 delete fAEMTotalRate;
00067 delete fAEMWeekData;
00068 delete fAEMTotalData;
00069 delete fAEMRateDist;
00070
00071 }
00072
00073
00074
00075 void AEMPlots::FillWeek(long fEpochStart, DBRunEntries& fRuns)
00076 {
00077 double fHourlySpills[168];
00078 double fHourlySpillTime[168];
00079 double fHourlyEvents[168];
00080 double fHourlyRate[168];
00081
00082 for (int i=0; i!=168; ++i) {
00083 fHourlySpills[i] = 0.;
00084 fHourlySpillTime[i] = 0.;
00085 fHourlyEvents[i] = 0.;
00086 fHourlyRate[i] = 0.;
00087 }
00088
00089 for (std::vector<DBRunEntry>::iterator it=fRuns.begin(); it !=fRuns.end(); ++it) {
00090 if (it->goodrun == 1) {
00091 int fRunStart = it->epochstart - fEpochStart;
00092 int fRunEnd = it->epochend - fEpochStart;
00093 double fRunLength = fRunEnd-fRunStart;
00094 if (fRunLength < 3600.) fRunLength = 3600.;
00095
00096 int fFirstHour;
00097 fFirstHour = (fRunStart<0 ? 0 : fRunStart/3600);
00098 for (int fHour=fFirstHour; (fHour!=(fRunEnd/3600)+1) && (fHour!=168); ++fHour) {
00099 double fRunSecThisHour, fHourStart, fHourEnd;
00100 fHourStart = fHour*3600.;
00101 fHourEnd = (fHour+1.)*3600. -1.;
00102
00103 fRunSecThisHour = 3600.;
00104 if (fRunStart>fHourStart) fRunSecThisHour -= (fRunStart-fHourStart);
00105 if (fRunEnd<fHourEnd) fRunSecThisHour -= (fHourEnd-fRunEnd);
00106
00107 fHourlyRate[fHour] += it->DAQRate() * fRunSecThisHour/3600.;
00108 fHourlySpillTime[fHour] += it->SpillTime() *fRunSecThisHour/fRunLength;
00109 fHourlySpills[fHour] += it->Spills() * fRunSecThisHour/fRunLength;
00110 fHourlyEvents[fHour] += it->SpillEvents() * fRunSecThisHour/fRunLength;
00111 }
00112
00113 fAEMRateDist->Fill(it->DAQRate());
00114 }
00115 }
00116
00117 double fWeekSpills, fWeekSpillTime, fWeekEvents;
00118 fWeekSpills = 0.;
00119 fWeekSpillTime = 0.;
00120 fWeekEvents = 0.;
00121
00122 for (int i=0; i!=168; ++i) {
00123 fAEMWeekSpills->Fill(i,fHourlySpills[i]);
00124 fAEMWeekSpillTime->Fill(i,fHourlySpillTime[i]);
00125 fAEMWeekEvents->Fill(i,fHourlyEvents[i]);
00126 fAEMWeekRate->Fill(i,fHourlyRate[i]);
00127 fWeekSpills += fHourlySpills[i];
00128 fWeekSpillTime += fHourlySpillTime[i];
00129 fWeekEvents += fHourlyEvents[i];
00130 }
00131
00132 std::vector<DBRunEntry>::iterator w=fRuns.begin();
00133 if ( w!=fRuns.end() && (w->epochstart - fEpochStart) < 0) ++w;
00134 std::string fTitle((w==fRuns.end()?"quiet week":w->datestart),0,10);
00135
00136 fAEMRateDist->SetTitle(fTitle.c_str());
00137 fAEMWeekSpills->SetTitle(fTitle.c_str());
00138 fAEMWeekSpillTime->SetTitle(fTitle.c_str());
00139 fAEMWeekEvents->SetTitle(fTitle.c_str());
00140 fAEMWeekRate->SetTitle(fTitle.c_str());
00141
00142 fAEMRateDist->SetFillColor(30);
00143 fAEMWeekSpills->SetFillColor(38);
00144 fAEMWeekSpillTime->SetFillColor(38);
00145 fAEMWeekEvents->SetFillColor(38);
00146 fAEMWeekRate->SetFillColor(30);
00147
00148 fAEMRateDist->SetXTitle("DAQ Rate");
00149 fAEMRateDist->SetYTitle("Entries");
00150 fAEMWeekSpills->SetXTitle("Hour");
00151 fAEMWeekSpills->SetYTitle("Spills per hour");
00152 fAEMWeekSpillTime->SetXTitle("Hour");
00153 fAEMWeekSpillTime->SetYTitle("Spill-seconds per hour");
00154 fAEMWeekEvents->SetXTitle("Hour");
00155 fAEMWeekEvents->SetYTitle("Events per hour");
00156 fAEMWeekRate->SetXTitle("Hour");
00157 fAEMWeekRate->SetYTitle("DAQ Rate");
00158
00159 fAEMWeekSpills->SetMinimum(0.);
00160 fAEMWeekSpillTime->SetMinimum(0.);
00161 fAEMWeekEvents->SetMinimum(0.);
00162 fAEMWeekRate->SetMinimum(0.);
00163
00164 fAEMWeekSpills->GetXaxis()->SetNdivisions(28,kFALSE);
00165 fAEMWeekSpillTime->GetXaxis()->SetNdivisions(28,kFALSE);
00166 fAEMWeekEvents->GetXaxis()->SetNdivisions(28,kFALSE);
00167 fAEMWeekRate->GetXaxis()->SetNdivisions(28,kFALSE);
00168
00169 fAEMWeekSpillTime->SetMaximum(360.);
00170 fAEMWeekRate->SetMaximum(60.);
00171
00172 std::cout << "Last week: \n" << fWeekEvents << " events in "
00173 << fWeekSpills << " spills with beam" << std::endl;
00174
00175 return;
00176 }
00177
00178
00179
00180 void AEMPlots::FillTotal(long fEpochStart, DBRunEntries& fRuns)
00181 {
00182 double fHourlySpills[fNbins];
00183 double fHourlySpillTime[fNbins];
00184 double fHourlyEvents[fNbins];
00185 double fHourlyRate[fNbins];
00186
00187 for (int i=0; i!=fNbins; ++i) {
00188 fHourlySpills[i] = 0.;
00189 fHourlySpillTime[i] = 0.;
00190 fHourlyEvents[i] = 0.;
00191 fHourlyRate[i] = 0.;
00192 }
00193
00194 for (std::vector<DBRunEntry>::iterator it=fRuns.begin(); it !=fRuns.end(); ++it) {
00195 if (it->goodrun == 1) {
00196 int fRunStart = it->epochstart - fEpochStart;
00197 int fRunEnd = it->epochend - fEpochStart;
00198 double fRunLength = fRunEnd-fRunStart;
00199 if (fRunLength < 3600.) fRunLength = 3600.;
00200
00201 int fFirstHour;
00202 fFirstHour = (fRunStart<0 ? 0 : fRunStart/3600);
00203 for (int fHour=fFirstHour; (fHour!=(fRunEnd/3600)+1) && (fHour!=fNbins); ++fHour) {
00204 double fRunSecThisHour, fHourStart, fHourEnd;
00205 fHourStart = fHour*3600.;
00206 fHourEnd = (fHour+1.)*3600. -1.;
00207
00208 fRunSecThisHour = 3600.;
00209 if (fRunStart>fHourStart) fRunSecThisHour -= (fRunStart-fHourStart);
00210 if (fRunEnd<fHourEnd) fRunSecThisHour -= (fHourEnd-fRunEnd);
00211
00212 fHourlyRate[fHour] += it->DAQRate() * fRunSecThisHour/3600.;
00213 fHourlySpillTime[fHour] += it->SpillTime() *fRunSecThisHour/fRunLength;
00214 fHourlySpills[fHour] += it->Spills() * fRunSecThisHour/fRunLength;
00215 fHourlyEvents[fHour] += it->SpillEvents() * fRunSecThisHour/fRunLength;
00216 }
00217
00218 fAEMRateDist->Fill(it->DAQRate());
00219 }
00220 }
00221
00222 double fTotalSpills, fTotalSpillTime, fTotalEvents;
00223 fTotalSpills = 0.;
00224 fTotalSpillTime = 0.;
00225 fTotalEvents = 0.;
00226
00227 for (int i=0; i!=fNbins; ++i) {
00228 fAEMTotalSpills->Fill(i,fHourlySpills[i]);
00229 fAEMTotalSpillTime->Fill(i,fHourlySpillTime[i]);
00230 fAEMTotalEvents->Fill(i,fHourlyEvents[i]);
00231 fAEMTotalRate->Fill(i,fHourlyRate[i]);
00232 fTotalSpills += fHourlySpills[i];
00233 fTotalSpillTime += fHourlySpillTime[i];
00234 fTotalEvents += fHourlyEvents[i];
00235 }
00236
00237 std::string fTitle("Total Run");
00238
00239 fAEMTotalSpills->SetTitle(fTitle.c_str());
00240 fAEMTotalSpillTime->SetTitle(fTitle.c_str());
00241 fAEMTotalEvents->SetTitle(fTitle.c_str());
00242 fAEMTotalRate->SetTitle(fTitle.c_str());
00243
00244 fAEMTotalSpills->SetFillColor(38);
00245 fAEMTotalSpillTime->SetFillColor(38);
00246 fAEMTotalEvents->SetFillColor(38);
00247 fAEMTotalRate->SetFillColor(30);
00248
00249 fAEMTotalSpills->SetXTitle("Hour");
00250 fAEMTotalSpills->SetYTitle("Spills per hour");
00251 fAEMTotalSpillTime->SetXTitle("Hour");
00252 fAEMTotalSpillTime->SetYTitle("Spill-seconds per hour");
00253 fAEMTotalEvents->SetXTitle("Hour");
00254 fAEMTotalEvents->SetYTitle("Events per hour");
00255 fAEMTotalRate->SetXTitle("Hour");
00256 fAEMTotalRate->SetYTitle("DAQ Rate");
00257
00258 fAEMTotalSpills->SetMinimum(0.);
00259 fAEMTotalSpillTime->SetMinimum(0.);
00260 fAEMTotalEvents->SetMinimum(0.);
00261 fAEMTotalRate->SetMinimum(0.);
00262
00263 fAEMTotalSpillTime->SetMaximum(360.);
00264 fAEMTotalRate->SetMaximum(60.);
00265
00266 std::cout << "Total: \n" << fTotalEvents << " events in "
00267 << fTotalSpills << " spills with beam" << std::endl;
00268
00269 return;
00270 }
