| 1 | /*void Fit(TGraph &g, double rate_start=3e5, double rate_kink=80, double rate_end=1)
|
|---|
| 2 | {
|
|---|
| 3 | double begin = -1;
|
|---|
| 4 | double end = -1;
|
|---|
| 5 |
|
|---|
| 6 | double kinkx = -1;
|
|---|
| 7 | double kinky = 0;
|
|---|
| 8 |
|
|---|
| 9 | TGraph gl;
|
|---|
| 10 | for (int j=0; j<g.GetN(); j++)
|
|---|
| 11 | {
|
|---|
| 12 | gl.SetPoint(gl.GetN(), g.GetX()[j], log10(g.GetY()[j]));
|
|---|
| 13 | if (g.GetY()[j]<rate_start && begin<0)
|
|---|
| 14 | begin = g.GetX()[j];
|
|---|
| 15 | if (g.GetY()[j]<rate_end && end<0)
|
|---|
| 16 | {
|
|---|
| 17 | end = g.GetX()[j];
|
|---|
| 18 | //if (end<550)
|
|---|
| 19 | // end = 750;
|
|---|
| 20 | }
|
|---|
| 21 | if (g.GetY()[j]<rate_kink && kinkx<0)
|
|---|
| 22 | {
|
|---|
| 23 | kinkx = g.GetX()[j];
|
|---|
| 24 | kinky = g.GetY()[j];
|
|---|
| 25 | }
|
|---|
| 26 | }
|
|---|
| 27 |
|
|---|
| 28 | if (end==-1)
|
|---|
| 29 | end = 600;//g.GetX()[g.GetN()-1];
|
|---|
| 30 |
|
|---|
| 31 | TF1 func2("f2", "log10("
|
|---|
| 32 | " [1]*exp([2]*(x-[0]))"
|
|---|
| 33 | "+ [1]*exp([3]*(x-[0]))"
|
|---|
| 34 | //"+ [4]*exp([5]*(x-700))"
|
|---|
| 35 | ")", begin, end);
|
|---|
| 36 |
|
|---|
| 37 | Double_t v00 = g.Eval(kinkx-150, 0, "S");
|
|---|
| 38 | Double_t v01 = g.Eval(kinkx- 50, 0, "S");
|
|---|
| 39 |
|
|---|
| 40 | Double_t v10 = g.Eval(kinkx+100, 0, "S");
|
|---|
| 41 | Double_t v11 = g.Eval(kinkx+200, 0, "S");
|
|---|
| 42 |
|
|---|
| 43 | //Double_t v20 = g.Eval(700-50, 0, "S");
|
|---|
| 44 | //Double_t v21 = g.Eval(700, 0, "S");
|
|---|
| 45 | //Double_t v22 = g.Eval(700+50, 0, "S");
|
|---|
| 46 |
|
|---|
| 47 | func2.SetParameter(0, kinkx);
|
|---|
| 48 | func2.SetParameter(1, kinky);
|
|---|
| 49 | func2.SetParameter(2, log(v01/v00)/100);
|
|---|
| 50 | func2.SetParameter(3, log(v11/v10)/100);
|
|---|
| 51 |
|
|---|
| 52 | func2.FixParameter(4, 0);
|
|---|
| 53 | func2.FixParameter(5, 0);
|
|---|
| 54 | func2.SetLineWidth(1);
|
|---|
| 55 | func2.SetLineColor(kBlue);
|
|---|
| 56 |
|
|---|
| 57 | gl.Fit(&func2, "N0", "", begin, end);
|
|---|
| 58 |
|
|---|
| 59 | TF1 func("f1",
|
|---|
| 60 | "[0]*exp([1]*(x-[6])) + "
|
|---|
| 61 | "[2]*exp([3]*(x-[6])) + "
|
|---|
| 62 | "[4]*exp([5]*(x-700))", begin, end);
|
|---|
| 63 | func.FixParameter(0, func2.GetParameter(1));
|
|---|
| 64 | func.FixParameter(1, func2.GetParameter(2));
|
|---|
| 65 | func.FixParameter(2, func2.GetParameter(1));
|
|---|
| 66 | func.FixParameter(3, func2.GetParameter(3));
|
|---|
| 67 | func.FixParameter(4, func2.GetParameter(4));
|
|---|
| 68 | func.FixParameter(5, func2.GetParameter(5));
|
|---|
| 69 | func.FixParameter(6, func2.GetParameter(0));
|
|---|
| 70 | func.SetLineWidth(2);
|
|---|
| 71 | func.SetLineColor(kBlue);
|
|---|
| 72 | func.DrawClone("same");
|
|---|
| 73 | func.SetLineWidth(1);
|
|---|
| 74 |
|
|---|
| 75 | cout << begin << " -- " << end << " : " << kinkx << endl;
|
|---|
| 76 |
|
|---|
| 77 | func.SetLineStyle(kDashed);
|
|---|
| 78 | func.SetRange(0, 1000);
|
|---|
| 79 |
|
|---|
| 80 | func.FixParameter(0, 0);
|
|---|
| 81 | func.FixParameter(2, 0);
|
|---|
| 82 | func.FixParameter(4, func2.GetParameter(4));
|
|---|
| 83 | func.DrawClone("same");
|
|---|
| 84 |
|
|---|
| 85 | func.FixParameter(0, func2.GetParameter(1));
|
|---|
| 86 | func.FixParameter(2, 0);
|
|---|
| 87 | func.FixParameter(4, 0);
|
|---|
| 88 | func.DrawClone("same");
|
|---|
| 89 |
|
|---|
| 90 | func.FixParameter(0, 0);
|
|---|
| 91 | func.FixParameter(2, func2.GetParameter(1));
|
|---|
| 92 | func.FixParameter(4, 0);
|
|---|
| 93 | func.DrawClone("same");
|
|---|
| 94 |
|
|---|
| 95 |
|
|---|
| 96 | TLatex txt;
|
|---|
| 97 | txt.SetTextColor(kBlue);
|
|---|
| 98 | txt.SetTextAlign(33);
|
|---|
| 99 | txt.DrawLatex(func2.GetParameter(0), func2.GetParameter(1),
|
|---|
| 100 | Form("%.1fHz @ %.0f", func2.GetParameter(1), func2.GetParameter(0)));
|
|---|
| 101 |
|
|---|
| 102 | // [1]*exp([2]*(x-[0])) / exp(1) = [1]*exp([3]*(x-[0]))
|
|---|
| 103 | // exp([2]*(x-[0])) / exp(1) = exp([3]*(x-[0]))
|
|---|
| 104 | // exp([2]*(x-[0])) / exp([3]*(x-[0])) = exp(1)
|
|---|
| 105 | // exp([2]*(x-[0]) - [3]*(x-[0])) = exp(1)
|
|---|
| 106 | // [2]*(x-[0]) - [3]*(x-[0]) = 1
|
|---|
| 107 | //
|
|---|
| 108 | // x = 1/([2]-[3]) + [0]
|
|---|
| 109 |
|
|---|
| 110 |
|
|---|
| 111 | Double_t th = 1. / (func2.GetParameter(3)-func2.GetParameter(2)) + func2.GetParameter(0);
|
|---|
| 112 | Double_t rate = pow(10, func2.Eval(th));
|
|---|
| 113 | Double_t rate_end = pow(10, func2.Eval(end)-0.5);
|
|---|
| 114 |
|
|---|
| 115 | txt.SetTextAlign(11);
|
|---|
| 116 | txt.DrawLatex(th+10, rate, Form("%.1fHz @ %.0f", rate, th));
|
|---|
| 117 |
|
|---|
| 118 | txt.SetTextAlign(12);
|
|---|
| 119 | txt.DrawLatex(180, 3e6, Form("%.1f", 1./func2.GetParameter(2)));
|
|---|
| 120 | txt.DrawLatex(end, rate_end, Form("%.0f", 1./func2.GetParameter(3)));
|
|---|
| 121 |
|
|---|
| 122 | TLine line;
|
|---|
| 123 | line.SetLineStyle(kDashed);
|
|---|
| 124 | line.SetLineColor(kBlue);
|
|---|
| 125 | line.DrawLine(th, 1e-2, th, 1e8);
|
|---|
| 126 | }*/
|
|---|
| 127 |
|
|---|
| 128 | Bool_t plot(MSQLMagic &serv, Long64_t search_id, TGraph &g, TGraph *gb)
|
|---|
| 129 | {
|
|---|
| 130 | TString query;
|
|---|
| 131 | query += "SELECT fTimeBegin, fTimeEnd, fVoltageIsOn, fOvervoltage, fCurrentMedMean, fNight ";
|
|---|
| 132 | query += "FROM Ratescan WHERE fRatescanID=";
|
|---|
| 133 | query += search_id;
|
|---|
| 134 |
|
|---|
| 135 | TSQLResult *res = serv.Query(query);
|
|---|
| 136 | if (!res)
|
|---|
| 137 | return kFALSE;
|
|---|
| 138 |
|
|---|
| 139 | if (res->GetRowCount()!=1)
|
|---|
| 140 | {
|
|---|
| 141 | cout << "ERROR - Row count is not 1." << endl;
|
|---|
| 142 | delete res;
|
|---|
| 143 | return kTRUE;
|
|---|
| 144 | }
|
|---|
| 145 |
|
|---|
| 146 | TSQLRow *row = res->Next();
|
|---|
| 147 |
|
|---|
| 148 | const char *time_beg = (*row)[0];
|
|---|
| 149 | const char *time_end = (*row)[1];
|
|---|
| 150 | const char *night = (*row)[5];
|
|---|
| 151 |
|
|---|
| 152 | int voltage_on = (*row)[2] ? atoi((*row)[2]) : -1;
|
|---|
| 153 | float overvoltage = (*row)[3] ? atof((*row)[3]) : -100;
|
|---|
| 154 | float current = (*row)[4] ? atof((*row)[4]) : -1;
|
|---|
| 155 |
|
|---|
| 156 | delete row;
|
|---|
| 157 |
|
|---|
| 158 | gROOT->SetSelectedPad(0);
|
|---|
| 159 | TCanvas *c = new TCanvas(time_beg+11, time_beg);
|
|---|
| 160 |
|
|---|
| 161 | c->SetGrid();
|
|---|
| 162 | c->SetLogy();
|
|---|
| 163 | c->SetTopMargin(0.01);
|
|---|
| 164 | c->SetRightMargin(0.005);
|
|---|
| 165 |
|
|---|
| 166 | TPaveText leg(600, 1e5, 1050, 8e8, "br");
|
|---|
| 167 | leg.SetBorderSize(1);
|
|---|
| 168 | leg.SetFillColor(kWhite);
|
|---|
| 169 | leg.AddText("Ratescan");
|
|---|
| 170 | leg.AddText("");
|
|---|
| 171 | leg.AddText(Form("Begin %s", time_beg));
|
|---|
| 172 | leg.AddText(Form("End %s", time_end));
|
|---|
| 173 | leg.AddText("");
|
|---|
| 174 | if (voltage_on==0)
|
|---|
| 175 | leg.AddText("Voltage off");
|
|---|
| 176 | else
|
|---|
| 177 | {
|
|---|
| 178 | if (current>=0)
|
|---|
| 179 | leg.AddText(Form("Current <I_{med}> = %.1f #muA", current));
|
|---|
| 180 | if (overvoltage>-70)
|
|---|
| 181 | leg.AddText(Form("Voltage #DeltaU = %+.2f V", overvoltage));
|
|---|
| 182 | }
|
|---|
| 183 |
|
|---|
| 184 | TH1D h0("frame", "", 1050, 0, 1050);
|
|---|
| 185 | h0.SetDirectory(0);
|
|---|
| 186 | h0.SetStats(kFALSE);
|
|---|
| 187 | h0.SetXTitle("Threshold [dac counts]");
|
|---|
| 188 | h0.SetYTitle("Trigger rate [Hz]");
|
|---|
| 189 | h0.GetXaxis()->SetLabelSize(0.05);
|
|---|
| 190 | h0.GetXaxis()->SetTitleSize(0.05);
|
|---|
| 191 | h0.GetYaxis()->SetLabelSize(0.05);
|
|---|
| 192 | h0.GetYaxis()->SetTitleSize(0.05);
|
|---|
| 193 | h0.SetMaximum(8e8);
|
|---|
| 194 | h0.SetMinimum(0.61);
|
|---|
| 195 | h0.DrawCopy();
|
|---|
| 196 |
|
|---|
| 197 | leg.DrawClone();
|
|---|
| 198 |
|
|---|
| 199 | for (int i=0; i<40; i++)
|
|---|
| 200 | {
|
|---|
| 201 | gb[i].SetNameTitle(Form("Board%02d", i), "Board trigger rate (time 40)");
|
|---|
| 202 | gb[i].DrawClone("P");
|
|---|
| 203 | }
|
|---|
| 204 |
|
|---|
| 205 | g.SetMarkerStyle(kFullDotMedium);
|
|---|
| 206 | g.SetMarkerColor(kBlue);
|
|---|
| 207 | g.SetLineColor(kBlue);
|
|---|
| 208 | g.SetNameTitle("Camera", "Camera trigger rate");
|
|---|
| 209 | g.DrawClone("PL");
|
|---|
| 210 |
|
|---|
| 211 | c->Write();
|
|---|
| 212 |
|
|---|
| 213 | TString name;
|
|---|
| 214 | name += night;
|
|---|
| 215 | name += "-ratescan ";
|
|---|
| 216 | name += time_beg;
|
|---|
| 217 |
|
|---|
| 218 | c->SaveAs(name+".pdf");
|
|---|
| 219 | c->SaveAs(name+".eps");
|
|---|
| 220 | c->SaveAs(name+".png");
|
|---|
| 221 |
|
|---|
| 222 | delete c;
|
|---|
| 223 | }
|
|---|
| 224 |
|
|---|
| 225 | Int_t plotid(MSQLMagic &serv, fits &file, const char *_search_id, Int_t row)
|
|---|
| 226 | {
|
|---|
| 227 | Long64_t search_id = atol(_search_id);
|
|---|
| 228 |
|
|---|
| 229 | cout << " " << search_id << endl;
|
|---|
| 230 |
|
|---|
| 231 | UInt_t offset = file.GetUInt("MJDREF");
|
|---|
| 232 |
|
|---|
| 233 | bool old = file.HasColumn("Data0");
|
|---|
| 234 |
|
|---|
| 235 | Double_t *ptime = file.SetPtrAddress("Time");
|
|---|
| 236 | UInt_t *pid = file.SetPtrAddress(old ? "Data0" : "Id");
|
|---|
| 237 | Float_t *rates = file.SetPtrAddress(old ? "Data5" : "BoardRate");
|
|---|
| 238 | UInt_t *pth = file.SetPtrAddress(old ? "Data1" : "Threshold");
|
|---|
| 239 | Float_t *ptrig = file.SetPtrAddress(old ? "Data4" : "TriggerRate");
|
|---|
| 240 | Float_t *pontime = file.SetPtrAddress(old ? "Data3" : "RelOnTime");
|
|---|
| 241 |
|
|---|
| 242 | if (!ptime || !pid || !rates || !pth || !ptrig || !pontime)
|
|---|
| 243 | return -1;
|
|---|
| 244 |
|
|---|
| 245 | TGraph g;
|
|---|
| 246 | TGraph gb[40];
|
|---|
| 247 |
|
|---|
| 248 | while (file.GetRow(row++))
|
|---|
| 249 | {
|
|---|
| 250 | if (pid[0]<search_id)
|
|---|
| 251 | continue;
|
|---|
| 252 |
|
|---|
| 253 | if (pid[0]!=search_id)
|
|---|
| 254 | break;
|
|---|
| 255 |
|
|---|
| 256 | g.SetPoint(g.GetN(), pth[0], ptrig[0]/pontime[0]);
|
|---|
| 257 | for (int i=0; i<40; i++)
|
|---|
| 258 | gb[i].SetPoint(gb[i].GetN(), pth[0], rates[i]*40);
|
|---|
| 259 | }
|
|---|
| 260 |
|
|---|
| 261 | if (g.GetN()>0)
|
|---|
| 262 | plot(serv, search_id, g, gb);
|
|---|
| 263 |
|
|---|
| 264 | return row;
|
|---|
| 265 |
|
|---|
| 266 |
|
|---|
| 267 | }
|
|---|
| 268 |
|
|---|
| 269 | Bool_t plotratescan(MSQLMagic &serv, const char *_night)
|
|---|
| 270 | {
|
|---|
| 271 | Long64_t night = atol(_night);
|
|---|
| 272 |
|
|---|
| 273 | TString fname = Form("/fact/aux/%04d/%02d/%02d/%06d.RATE_SCAN_DATA.fits",
|
|---|
| 274 | night/10000, (night/100)%100, night%100, night);
|
|---|
| 275 |
|
|---|
| 276 | cout << " " << fname << endl;
|
|---|
| 277 |
|
|---|
| 278 | fits file(fname.Data());
|
|---|
| 279 | if (!file)
|
|---|
| 280 | {
|
|---|
| 281 | cout << "ERROR - Cannot access " << fname << ": " << gSystem->GetError() << endl;
|
|---|
| 282 | return kFALSE;
|
|---|
| 283 | }
|
|---|
| 284 |
|
|---|
| 285 | TString query;
|
|---|
| 286 | query += "SELECT fRatescanID FROM Ratescan WHERE fNight=";
|
|---|
| 287 | query += night;
|
|---|
| 288 | query += " ORDER BY fRatescanID";
|
|---|
| 289 | /*
|
|---|
| 290 | query += "SELECT fRatescanID, fTimeBeg, fTimeEnd, fVoltageIsOn, ";
|
|---|
| 291 | query += " fOvervoltage, fCurrentMedMean, fRateBegin, fRateEnd, ";
|
|---|
| 292 | query += " fThresholdBegin, fThresholdEnd, fNumPoints ";
|
|---|
| 293 | query += "WHERE fNight=";
|
|---|
| 294 | query += night;
|
|---|
| 295 | */
|
|---|
| 296 | TSQLResult *res = serv.Query(query);
|
|---|
| 297 | if (!res)
|
|---|
| 298 | return kFALSE;
|
|---|
| 299 |
|
|---|
| 300 | if (res->GetRowCount()==0)
|
|---|
| 301 | {
|
|---|
| 302 | delete res;
|
|---|
| 303 | return kTRUE;
|
|---|
| 304 | }
|
|---|
| 305 |
|
|---|
| 306 | TString oname = Form("%06d-ratescan.root", night);
|
|---|
| 307 |
|
|---|
| 308 | cout << " " << oname << '\n' << endl;
|
|---|
| 309 |
|
|---|
| 310 | TFile rootfile(oname.Data(), "recreate");
|
|---|
| 311 |
|
|---|
| 312 | TSQLRow *row = 0;
|
|---|
| 313 |
|
|---|
| 314 | Int_t cnt = 0;
|
|---|
| 315 | while ((row=res->Next()) && cnt>=0)
|
|---|
| 316 | {
|
|---|
| 317 | const char *id = (*row)[0];
|
|---|
| 318 | cnt = plotid(serv, file, id, cnt);
|
|---|
| 319 | delete row;
|
|---|
| 320 | }
|
|---|
| 321 |
|
|---|
| 322 | delete res;
|
|---|
| 323 |
|
|---|
| 324 | cout << endl;
|
|---|
| 325 |
|
|---|
| 326 | return cnt>=0;
|
|---|
| 327 | }
|
|---|
| 328 |
|
|---|
| 329 | Bool_t plotratescan(const char *night)
|
|---|
| 330 | {
|
|---|
| 331 | MSQLMagic serv("sql.rc");
|
|---|
| 332 | if (!serv.IsConnected())
|
|---|
| 333 | {
|
|---|
| 334 | cout << "ERROR - Connection to database failed." << endl;
|
|---|
| 335 | return 0;
|
|---|
| 336 | }
|
|---|
| 337 |
|
|---|
| 338 | cout << "plotratescan" << endl;
|
|---|
| 339 | cout << "------------" << endl;
|
|---|
| 340 | cout << endl;
|
|---|
| 341 | cout << "Connected to " << serv.GetName() << endl;
|
|---|
| 342 | cout << "Night: " << night << endl;
|
|---|
| 343 | cout << endl;
|
|---|
| 344 |
|
|---|
| 345 | return plotratescan(serv, night);
|
|---|
| 346 | }
|
|---|
| 347 |
|
|---|
| 348 | Bool_t plotratescan()
|
|---|
| 349 | {
|
|---|
| 350 | MSQLMagic serv("sql.rc");
|
|---|
| 351 | if (!serv.IsConnected())
|
|---|
| 352 | {
|
|---|
| 353 | cout << "ERROR - Connection to database failed." << endl;
|
|---|
| 354 | return 0;
|
|---|
| 355 | }
|
|---|
| 356 |
|
|---|
| 357 | cout << "plotratescan" << endl;
|
|---|
| 358 | cout << "------------" << endl;
|
|---|
| 359 | cout << endl;
|
|---|
| 360 | cout << "Connected to " << serv.GetName() << endl;
|
|---|
| 361 | cout << endl;
|
|---|
| 362 |
|
|---|
| 363 | TSQLResult *res = serv.Query("SELECT fNight FROM Ratescan GROUP BY fNight ORDER BY fNight");
|
|---|
| 364 | if (!res)
|
|---|
| 365 | return kFALSE;
|
|---|
| 366 |
|
|---|
| 367 | TSQLRow *row = 0;
|
|---|
| 368 | while ((row=res->Next()))
|
|---|
| 369 | {
|
|---|
| 370 | const char *night = (*row)[0];
|
|---|
| 371 | plotratescan(serv, night);
|
|---|
| 372 | delete row;
|
|---|
| 373 | }
|
|---|
| 374 |
|
|---|
| 375 | delete res;
|
|---|
| 376 |
|
|---|
| 377 | return kTRUE;
|
|---|
| 378 | }
|
|---|
