Index: /trunk/Mars/fact/analysis/lightcurve.C =================================================================== --- /trunk/Mars/fact/analysis/lightcurve.C (revision 17103) +++ /trunk/Mars/fact/analysis/lightcurve.C (revision 17104) @@ -3,7 +3,8 @@ #include #include +#include #include #include - +#include #include #include @@ -31,5 +32,5 @@ -int lightcurve(Int_t sourcekey=1, Int_t nightmin=0, Int_t nightmax=0, Int_t minutes=20, TString table="AnalysisResultsRunLP", TString outfile="lc") +int lightcurve(Int_t sourcekey=1, Int_t nightmin=0, Int_t nightmax=0, Int_t minutes=20, TString table="AnalysisResultsRunLP", TString outfile="lc", Bool_t dch=kTRUE) { @@ -54,5 +55,22 @@ query+=Form(" AND fNight <=%d",nightmax); - cout << "Q: " << query << endl; + // datacheck + TString datacheck=" "; + //remove data with wrong settings + datacheck+=" AND fNight>20120420 AND NOT fNight IN (20120406,20120410,20120503) AND";//data with different bias voltage + datacheck+=" NOT fNight BETWEEN 20121206 AND 20130110"; // broken bias channel + //datacheck+=" AND NOT (fNight=20120608 AND fRunId=65) "; // something went wrong with tracking? + // 24.6. new coefficients + TString zdparam=" pow(0.753833 * cos(Radians(fZenithDistanceMean)), 7.647435) * exp(-5.753686*pow(Radians(fZenithDistanceMean),2.089609))"; + TString thparam=" pow((if(isnull(fThresholdMinSet),fThresholdMedian,fThresholdMinSet)-329.4203),2) * (-0.0000002044803) "; + TString param=" (fNumEvtsAfterBgCuts/5-fNumSigEvts)/fOnTimeAfterCuts - "+zdparam+" - "+thparam+" "; + datacheck+=" AND -0.085 < ("+param+") "; + datacheck+=" AND 0.25 > ("+param+") "; + + + if (dch) + query+=datacheck; + + //cout << "Q: " << query << endl; TSQLResult *res = serv.Query(query); @@ -69,5 +87,4 @@ query+=Form(" WHERE fSourceKey=%d", sourcekey); query+=Form(" AND fNight BETWEEN %s AND %s ", nightstart.Data(), nightstop.Data()); - //query+=" AND NOT ISNULL(fOnTime) ";//at the moment filled at ISDC query+=" GROUP BY fNight "; query+=" ORDER BY fNight "; @@ -80,6 +97,9 @@ return 1; - Int_t counter=0; - Int_t counter2=0; + //counter + Int_t counter=0;//counter for graphs + Int_t counter2=0;//counter for loop + + //graphs and histograms TGraphErrors excessrate; TGraphErrors backgroundrate; @@ -87,10 +107,18 @@ TGraphErrors significancerate; TGraphErrors excratevsbgrate; + TGraphErrors excvsbg; TGraphErrors excvszd; - TGraphErrors excvsbg; TGraphErrors bgvszd; + TGraphErrors excvscur; + TGraphErrors bgvscur; + TGraphErrors excvsth; + TGraphErrors bgvsth; + TGraphErrors zdvstime; + TGraphErrors curvstime; + TGraphErrors thvstime; TH1F exc("exc","rates", 43, -12.5, 202.5); TH1F bg("bg","rates", 43, -12.5, 202.5); + //variables for calculations and graphs Float_t significance; TString excevts; @@ -98,7 +126,8 @@ TString sigevts; TString ontime; - TString zd; TString zd2; TString zd3; + TString cur; + TString th; Float_t excevtssum=0; Float_t bgevtssum=0; @@ -109,25 +138,42 @@ Float_t signifrate=0; Float_t excerr=0; - Float_t zdsum=0; + Float_t bgerr=0; Float_t zdmean=0; + Float_t curmean=0; + Float_t thmean=0; + Float_t zddiff=0; + Float_t curdiff=0; + Float_t thdiff=0; Float_t zdmin=0; Float_t zdmax=0; + Float_t curmin=0; + Float_t curmax=0; + Float_t thmin=0; + Float_t thmax=0; MTime start; MTime stop; - Double_t mjdstart; - Double_t mjdstop; - Double_t mjd; + Double_t mjdstart=0; + Double_t mjdstop=0; + Double_t mjd=0; + Double_t mjddiff=0; + TSQLRow *row3=0; while ((row3=res3->Next())) { cout << "NIGHT " << (*row3)[0] << endl; - query =" SELECT IF(ISNULL(fOnTimeAfterCuts), IF(ISNULL(fOnTime), Time_to_sec(Timediff(fRunStop,fRunStart)), fOnTime), fOnTimeAfterCuts), "; - query+=" fNumBgEvts, fNumSigEvts, fNumExcEvts, "; - query+=" fRunStart, fRunStop, (fZenithDistanceMax-fZenithDistanceMin)/2, "; - query+=" fZenithDistanceMax, fZenithDistanceMin FROM "+table; + //query =" SELECT IF(ISNULL(fOnTimeAfterCuts), IF(ISNULL(fOnTime), Time_to_sec(Timediff(fRunStop,fRunStart)), fOnTime), fOnTimeAfterCuts), "; + query =" SELECT fOnTimeAfterCuts, fNumBgEvts, fNumSigEvts, fNumExcEvts, "; + query+=" fRunStart, fRunStop, "; + query+=" fZenithDistanceMax, fZenithDistanceMin, "; + query+=" fCurrentsMedMeanBeg, fThresholdMinSet FROM "+table; query+=" LEFT JOIN RunInfo USING(fNight, fRunId) "; query+=Form(" WHERE fSourceKey=%d", sourcekey); query+=Form(" AND fNight= %s ", (*row3)[0]); - //query+=" AND NOT ISNULL(fOnTime) ";//at the moment filled at isdc + query+=" AND fOnTimeAfterCuts < 1000 "; //exclude runs with wrong/too high ontime + query+=" AND NOT ISNULL(fNumExcEvts) ";// only where excess was extracted + + if (dch) + query+=datacheck; + query+=" ORDER BY fRunId "; @@ -147,122 +193,184 @@ { counter2++; - cout << "-----" << (*row2)[0] << " " << (*row2)[3] << endl; ontime=(*row2)[0]; bgevts=(*row2)[1]; sigevts=(*row2)[2]; excevts=(*row2)[3]; - zd=(*row2)[6]; - zd2=(*row2)[7]; - zd3=(*row2)[8]; - start.SetSqlDateTime((*row2)[4]); - stop.SetSqlDateTime((*row2)[5]); + zd2=(*row2)[6];//zdmin + zd3=(*row2)[7];//zdmax + cur=(*row2)[8];//currents + th=(*row2)[9];//threshold if (counter2==1) { - mjdstart=start.GetMjd(); - cout << "mjdstart: " << (*row2)[4] << endl; + start.SetSqlDateTime((*row2)[4]); zdmin=zd3.Atof(); zdmax=zd2.Atof(); + curmin=cur.Atof(); + curmax=cur.Atof(); + thmin=th.Atof(); + thmax=th.Atof(); } - cout << "# " << counter2 << " " << (*row2)[0] << " " << (*row2)[1] << " " << (*row2)[2] << " " << (*row2)[3] << endl; - - if (ontimesum+ontime.Atof()>minutes) + if (ontimesum+ontime.Atof()>minutes && minutes>0) { - cout << "reached limit" << endl; - significance=MMath::SignificanceLiMa(sigevtssum, bgevtssum, 0.2); - cout << " significance: " << significance << endl; - cout << " excess: " << excevtssum << " evts" << endl; - cout << " bg: " << bgevtssum << " evts" << endl; - cout << " ontime: " << ontimesum << " s " << endl; - cout << " excessrate: " << excevtssum/ontimesum*3600 << " evts/h" << endl; - cout << " bgrate: " << bgevtssum/ontimesum*3600 << " evts/h" << endl; - signifrate=significance/sqrt(ontimesum*3600); - bgrate=bgevtssum/ontimesum*3600; - excrate=excevtssum/ontimesum*3600; - excerr = MMath::ErrorExc(significance, bgevtssum, 0.2)/ontimesum*3600.; - //zdmean=zdsum/counter2; - zdmean=zdmin+(zdmax-zdmin)/2; - mjdstop=stop.GetMjd(); - mjd=mjdstart+(mjdstop-mjdstart)/2; - cout << " mjd: " << mjd << " (" << mjdstart << " - " << mjdstop << ") " << endl; - if (excrate>300) - cout << "==>" << (*row2)[0] << " " << (*row2)[1] << " " << (*row2)[2] - << " " << (*row2)[3] << " " << (*row2)[4] << " " << (*row2)[5] - << " " << (*row2)[6] << " " << (*row2)[7] << " " << (*row2)[8] - << endl; - + //calculate values + significance = MMath::SignificanceLiMaSigned(sigevtssum, bgevtssum*5, 0.2); + signifrate = significance/sqrt(ontimesum/3600); + bgrate = bgevtssum/ontimesum*3600; + excrate = excevtssum/ontimesum*3600; + excerr = MMath::ErrorExc(sigevtssum, bgevtssum*5, 0.2)/ontimesum*3600.; + bgerr = sqrt(bgevtssum)/ontimesum*3600; + zdmean = zdmin+(zdmax-zdmin)/2; + zddiff = (zdmax-zdmin)/2; + curmean = curmin+(curmax-curmin)/2; + curdiff = (curmax-curmin)/2; + thmean = thmin+(thmax-thmin)/2; + thdiff = (thmax-thmin)/2; + mjdstart = start.GetMjd(); + mjdstop = stop.GetMjd(); + mjd = mjdstart+(mjdstop-mjdstart)/2.; + mjddiff = ontimesum/3600./24/2; + + //fill histograms and graphs + exc.Fill(excrate); + bg.Fill(bgrate); + excessrate.SetPoint(counter, mjd, excrate); + excessrate.SetPointError(counter, mjddiff, excerr); + backgroundrate.SetPoint(counter, mjd, bgrate); + backgroundrate.SetPointError(counter, mjddiff, bgerr); + excvsbg.SetPoint(counter, bgrate, excrate); + excvsbg.SetPointError(counter, bgerr, excerr); + excvszd.SetPoint(counter, zdmean, excrate); + excvszd.SetPointError(counter, zddiff, excerr); + bgvszd.SetPoint(counter, zdmean, bgrate); + bgvszd.SetPointError(counter, zddiff, bgerr); + excvscur.SetPoint(counter, curmean, excrate); + excvscur.SetPointError(counter, curdiff, excerr); + bgvscur.SetPoint(counter, curmean, bgrate); + bgvscur.SetPointError(counter, curdiff, bgerr); + excvsth.SetPoint(counter, thmean, excrate); + excvsth.SetPointError(counter, thdiff, excerr); + bgvsth.SetPoint(counter, thmean, bgrate); + bgvsth.SetPointError(counter, thdiff, bgerr); + signif.SetPoint(counter, mjd, significance); + signif.SetPointError(counter, mjddiff, 0); + significancerate.SetPoint(counter, mjd, signifrate); + significancerate.SetPointError(counter, mjddiff, 0); + zdvstime.SetPoint(counter, mjd, zdmean); + zdvstime.SetPointError(counter, mjddiff, zddiff); + curvstime.SetPoint(counter, mjd, curmean); + curvstime.SetPointError(counter, mjddiff, curdiff); + thvstime.SetPoint(counter, mjd, thmean); + thvstime.SetPointError(counter, mjddiff, thdiff); + + //set couter for graphs + counter++; + //reset counter for calculation + counter2=0; + //reset variables + start.SetSqlDateTime((*row2)[4]); excevtssum=0; bgevtssum=0; sigevtssum=0; ontimesum=0; - zdsum=0; - exc.Fill(excrate); - bg.Fill(bgrate); - excessrate.SetPoint(counter, mjd, excrate); - excessrate.SetPointError(counter, ontimesum/3600./24/2, excerr); - backgroundrate.SetPoint(counter, mjd, bgrate); - backgroundrate.SetPointError(counter, ontimesum/3600./24/2, sqrt(bgevtssum)/ontimesum*3600); - excvszd.SetPoint(counter, zdmean, excrate); - excvszd.SetPointError(counter, (zdmax-zdmin)/2, excerr); - excvsbg.SetPoint(counter, bgrate, excrate); - excvsbg.SetPointError(counter, sqrt(bgevtssum)/ontimesum*3600, excerr); - bgvszd.SetPoint(counter, zdmean, bgrate); - bgvszd.SetPointError(counter, (zdmax-zdmin)/2, sqrt(bgevtssum)/ontimesum*3600); - signif.SetPoint(counter, mjd, significance); - signif.SetPointError(counter, ontimesum/3600./24/2, 0); - significancerate.SetPoint(counter, mjd, signifrate); - significancerate.SetPointError(counter, ontimesum/3600./24/2, 0); - counter++; - counter2=0; } + //set stop time + stop.SetSqlDateTime((*row2)[5]); + //add up variables ontimesum+=ontime.Atof(); excevtssum+=excevts.Atof(); bgevtssum+=bgevts.Atof(); sigevtssum+=sigevts.Atof(); - zdsum+=zd.Atof(); + + //update zdmin if (zdmin>zd3.Atof()) zdmin=zd3.Atof(); + //update zdmax if (zdmax on " << ontimesum << " exc " << excevtssum << " bg " << bgevtssum << " sig " << sigevtssum << endl; + + //update curmin + if (curmin>cur.Atof()) + curmin=cur.Atof(); + //update curmax + if (curmaxth.Atof()) + thmin=th.Atof(); + //update thmax + if (thmax(minutes-0.1*minutes)) + if (ontimesum>(minutes-0.1*minutes) || minutes<0) { - cout << "ontime > " << minutes-0.1*minutes << endl; - significance=MMath::SignificanceLiMa(sigevtssum, bgevtssum, 0.2); - cout << " significance: " << significance << endl; - cout << " excess: " << excevtssum << " evts" << endl; - cout << " bg: " << bgevtssum << " evts" << endl; - cout << " ontime: " << ontimesum << " s " << endl; - cout << " excessrate: " << excevtssum/ontimesum*3600 << " evts/h" << endl; - cout << " bgrate: " << bgevtssum/ontimesum*3600 << " evts/h" << endl; - signifrate=significance/sqrt(ontimesum*3600); - bgrate=bgevtssum/ontimesum*3600; - excrate=excevtssum/ontimesum*3600; - excerr = MMath::ErrorExc(significance, bgevtssum, 0.2)/ontimesum*3600.; - zdmean=zdmin+(zdmax-zdmin)/2; - mjdstop=stop.GetMjd(); - mjd=mjdstart+(mjdstop-mjdstart)/2; - cout << " mjd: " << mjd << " (" << mjdstart << " - " << mjdstop << ") " << endl; - - exc.Fill(excrate); - bg.Fill(bgrate); - excessrate.SetPoint(counter, mjd, excrate); - excessrate.SetPointError(counter, ontimesum/3600./24/2, excerr); - backgroundrate.SetPoint(counter, mjd, bgrate); - backgroundrate.SetPointError(counter, ontimesum/3600./24/2, sqrt(bgevtssum)/ontimesum*3600); - excvszd.SetPoint(counter, zdmean, excrate); - excvszd.SetPointError(counter, (zdmax-zdmin)/2, excerr); - excvsbg.SetPoint(counter, bgrate, excrate); - excvsbg.SetPointError(counter, sqrt(bgevtssum)/ontimesum*3600, excerr); - bgvszd.SetPoint(counter, zdmean, bgrate); - bgvszd.SetPointError(counter, (zdmax-zdmin)/2, sqrt(bgevtssum)/ontimesum*3600); - signif.SetPoint(counter, mjd, significance); - signif.SetPointError(counter, ontimesum/3600./24/2, 0); - significancerate.SetPoint(counter, mjd, signifrate); - significancerate.SetPointError(counter, ontimesum/3600./24/2, 0); + if (minutes<0 && ontimesum>20*60 && bgevtssum!=0) + { + //calculate values + significance = MMath::SignificanceLiMaSigned(sigevtssum, bgevtssum*5, 0.2); + signifrate = significance/sqrt(ontimesum/3600); + bgrate = bgevtssum/ontimesum*3600; + excrate = excevtssum/ontimesum*3600; + excerr = MMath::ErrorExc(sigevtssum, bgevtssum*5, 0.2)/ontimesum*3600.; + bgerr = sqrt(bgevtssum)/ontimesum*3600; + zdmean = zdmin+(zdmax-zdmin)/2; + zddiff = (zdmax-zdmin)/2; + curmean = curmin+(curmax-curmin)/2; + curdiff = (curmax-curmin)/2; + thmean = thmin+(thmax-thmin)/2; + thdiff = (thmax-thmin)/2; + mjdstart = start.GetMjd(); + mjdstop = stop.GetMjd(); + mjd = mjdstart+(mjdstop-mjdstart)/2.; + mjddiff = ontimesum/3600./24/2; + + //fill histograms and graphs + exc.Fill(excrate); + bg.Fill(bgrate); + excessrate.SetPoint(counter, mjd, excrate); + excessrate.SetPointError(counter, mjddiff, excerr); + backgroundrate.SetPoint(counter, mjd, bgrate); + backgroundrate.SetPointError(counter, mjddiff, bgerr); + excvsbg.SetPoint(counter, bgrate, excrate); + excvsbg.SetPointError(counter, bgerr, excerr); + excvszd.SetPoint(counter, zdmean, excrate); + excvszd.SetPointError(counter, zddiff, excerr); + bgvszd.SetPoint(counter, zdmean, bgrate); + bgvszd.SetPointError(counter, zddiff, bgerr); + excvscur.SetPoint(counter, curmean, excrate); + excvscur.SetPointError(counter, curdiff, excerr); + bgvscur.SetPoint(counter, curmean, bgrate); + bgvscur.SetPointError(counter, curdiff, bgerr); + excvsth.SetPoint(counter, thmean, excrate); + excvsth.SetPointError(counter, thdiff, excerr); + bgvsth.SetPoint(counter, thmean, bgrate); + bgvsth.SetPointError(counter, thdiff, bgerr); + signif.SetPoint(counter, mjd, significance); + signif.SetPointError(counter, mjddiff, 0); + significancerate.SetPoint(counter, mjd, signifrate); + significancerate.SetPointError(counter, mjddiff, 0); + zdvstime.SetPoint(counter, mjd, zdmean); + zdvstime.SetPointError(counter, mjddiff, zddiff); + curvstime.SetPoint(counter, mjd, curmean); + curvstime.SetPointError(counter, mjddiff, curdiff); + thvstime.SetPoint(counter, mjd, thmean); + thvstime.SetPointError(counter, mjddiff, thdiff); + + //set counter for graphs + counter++; + + //reset variables + excevtssum=0; + bgevtssum=0; + sigevtssum=0; + ontimesum=0; + } + + //reset for calculations + counter2=0; } delete res2; @@ -275,60 +383,130 @@ MStatusDisplay *d = new MStatusDisplay; +// if (dch) +// backgroundrate.SetTitle(Form("Rates vs MJD (dch: %s)", limit.Data())); +// else + + + Double_t x; + Double_t y; + Double_t min; + Double_t max; + + TText txt; + txt.SetTextColor(kGray+2); + txt.SetTextSize(0.07); + TString text="FACT Preliminary"; + + + + gROOT->SetSelectedPad(0); TCanvas &c1 = d->AddTab("Rates vs MJD", "Rates vs MJD"); gPad->SetGridy(); - excessrate.SetTitle("Rates vs MJD"); - excessrate.SetMarkerStyle(7); - excessrate.SetMarkerSize(2); - excessrate.GetYaxis()->SetTitle("evts/h"); - excessrate.GetXaxis()->SetTitle("MJD"); - excessrate.GetXaxis()->CenterTitle(); - excessrate.GetYaxis()->SetRangeUser(-50,400); - excessrate.DrawClone("AP"); - - backgroundrate.SetMarkerStyle(7); - backgroundrate.SetMarkerSize(2); + + excessrate.GetPoint(1,x,y); + max = excessrate.GetHistogram()->GetMaximum(); + min = excessrate.GetHistogram()->GetMinimum(); + + backgroundrate.SetTitle("Rates vs MJD"); + backgroundrate.SetMarkerStyle(8); + backgroundrate.SetMarkerSize(0.8); + backgroundrate.SetLineWidth(2); backgroundrate.SetMarkerColor(kBlue); backgroundrate.SetLineColor(kBlue); - backgroundrate.DrawClone("Psame"); - - TCanvas &c2 = d->AddTab("Excess Rate vs MJD", "Excess Rate vs MJD"); + backgroundrate.GetYaxis()->SetTitle("evts/h"); + backgroundrate.GetXaxis()->SetTitle("MJD"); + backgroundrate.GetXaxis()->CenterTitle(); + backgroundrate.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); + backgroundrate.DrawClone("AP"); + + excessrate.SetMarkerStyle(8); + excessrate.SetMarkerSize(0.8); + excessrate.SetLineWidth(2); + excessrate.DrawClone("Psame"); + + txt.DrawText(x, max, text); + + + + gROOT->SetSelectedPad(0); + TCanvas &c2 = d->AddTab("Excess vs MJD", "Excess vs MJD"); gPad->SetGridy(); excessrate.SetTitle("Excess Rate vs MJD"); excessrate.SetMarkerStyle(7); excessrate.SetMarkerSize(2); - excessrate.GetYaxis()->SetTitle("exc evts/h"); + excessrate.GetYaxis()->SetTitle("excess evts/h"); excessrate.GetXaxis()->SetTitle("MJD"); excessrate.GetXaxis()->CenterTitle(); + excessrate.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); excessrate.DrawClone("AP"); - + txt.DrawText(x, max, text); + + + + gROOT->SetSelectedPad(0); + TCanvas &c4 = d->AddTab("Exc vs Bg", "Exc vs Bg"); + + backgroundrate.GetPoint(1,x,y); + max = excvsbg.GetHistogram()->GetMaximum(); + min = excvsbg.GetHistogram()->GetMinimum(); + + excvsbg.SetTitle("Excess Rate vs Background Rate"); + excvsbg.SetMarkerStyle(7); + excvsbg.SetMarkerSize(2); + excvsbg.GetYaxis()->SetTitle("excess evts/h"); + excvsbg.GetXaxis()->SetTitle("background evts/h"); + excvsbg.GetXaxis()->CenterTitle(); + excvsbg.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); + excvsbg.DrawClone("AP"); + txt.DrawText(excvsbg.GetMean(1), max, text); + + + + gROOT->SetSelectedPad(0); TCanvas &c3 = d->AddTab("Exc vs Zd", "Exc vs Zd"); + + excvszd.GetPoint(1,x,y); + max = excvszd.GetHistogram()->GetMaximum(); + min = excvszd.GetHistogram()->GetMinimum(); + excvszd.SetTitle("Excess Rate vs Zd"); excvszd.SetMarkerStyle(7); excvszd.SetMarkerSize(2); - excvszd.GetYaxis()->SetTitle("exc evts/h"); - excvszd.GetXaxis()->SetTitle("zd"); + excvszd.GetYaxis()->SetTitle("excess evts/h"); + excvszd.GetXaxis()->SetTitle("deg"); excvszd.GetXaxis()->CenterTitle(); + excvszd.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); excvszd.DrawClone("AP"); - - TCanvas &c4 = d->AddTab("Exc vs Bg", "Exc vs Bg"); - excvsbg.SetTitle("Excess Rate vs Background Rate"); - excvsbg.SetMarkerStyle(7); - excvsbg.SetMarkerSize(2); - excvsbg.GetYaxis()->SetTitle("exc evts/h"); - excvsbg.GetXaxis()->SetTitle("bg evts/h"); - excvsbg.GetXaxis()->CenterTitle(); - excvsbg.DrawClone("AP"); - + txt.DrawText(x, max, text); + + + + gROOT->SetSelectedPad(0); TCanvas &c5 = d->AddTab("Bg vs Zd", "Bg vs Zd"); + + zdvstime.GetPoint(1,x,y); + max = bgvszd.GetHistogram()->GetMaximum(); + min = bgvszd.GetHistogram()->GetMinimum(); + bgvszd.SetTitle("Background Rate vs Zd"); bgvszd.SetMarkerStyle(7); bgvszd.SetMarkerSize(2); - bgvszd.GetYaxis()->SetTitle("bg evts/h"); - bgvszd.GetXaxis()->SetTitle("zd"); + bgvszd.GetYaxis()->SetTitle("background evts/h"); + bgvszd.GetXaxis()->SetTitle("deg"); bgvszd.GetXaxis()->CenterTitle(); + bgvszd.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); bgvszd.DrawClone("AP"); - + txt.DrawText(y, max, text); + + + + gROOT->SetSelectedPad(0); TCanvas &c6 = d->AddTab("Signifrate vs MJD", "Signifrate vs MJD"); gPad->SetGridy(); + + significancerate.GetPoint(1,x,y); + max = significancerate.GetHistogram()->GetMaximum(); + min = significancerate.GetHistogram()->GetMinimum(); + significancerate.SetTitle("Significance Rate vs MJD"); significancerate.SetMarkerStyle(7); @@ -337,7 +515,17 @@ significancerate.GetXaxis()->SetTitle("MJD"); significancerate.GetXaxis()->CenterTitle(); + significancerate.GetYaxis()->SetRangeUser(0,max*1.1); significancerate.DrawClone("AP"); - - TCanvas &c7 = d->AddTab("Significance vs MJD", "Significance vs MJD"); + txt.DrawText(x, max, text); + + + + gROOT->SetSelectedPad(0); + TCanvas &c7 = d->AddTab("Signif vs MJD", "Signif vs MJD"); + + signif.GetPoint(1,x,y); + max = signif.GetHistogram()->GetMaximum(); + min = signif.GetHistogram()->GetMinimum(); + signif.SetTitle("Significance vs MJD"); signif.SetMarkerStyle(7); @@ -346,10 +534,161 @@ signif.GetXaxis()->SetTitle("MJD"); signif.GetXaxis()->CenterTitle(); + signif.GetYaxis()->SetRangeUser(0,max*1.1); signif.DrawClone("AP"); - + txt.DrawText(x, max, text); + + + + gROOT->SetSelectedPad(0); TCanvas &c8 = d->AddTab("Rate Hist", "Rate Hist"); + + max = bg.GetMaximum()*1.1; + min = bg.GetMinimum(); + bg.SetLineColor(kBlue); + bg.GetYaxis()->SetRangeUser(0,max*1.1); + bg.GetXaxis()->SetTitle("evts/h"); bg.DrawCopy("hist"); exc.DrawCopy("histsame"); + txt.DrawText(min, max, text); + + + + gROOT->SetSelectedPad(0); + TCanvas &c9 = d->AddTab("Zd vs Time", "Zd vs Time"); + gPad->SetGridy(); + + zdvstime.GetPoint(1,x,y); + max = zdvstime.GetHistogram()->GetMaximum(); + min = zdvstime.GetHistogram()->GetMinimum(); + + zdvstime.SetTitle("Zd vs MJD"); + zdvstime.SetMarkerStyle(7); + zdvstime.SetMarkerSize(2); + zdvstime.GetYaxis()->SetTitle("deg"); + zdvstime.GetXaxis()->SetTitle("MJD"); + zdvstime.GetXaxis()->CenterTitle(); + zdvstime.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); + zdvstime.DrawClone("AP"); + txt.DrawText(x, max, text); + + + + gROOT->SetSelectedPad(0); + TCanvas &c10 = d->AddTab("Cur vs Time", "Cur vs Time"); + gPad->SetGridy(); + + curvstime.GetPoint(1,x,y); + max = curvstime.GetHistogram()->GetMaximum(); + min = curvstime.GetHistogram()->GetMinimum(); + + curvstime.SetTitle("Currents vs MJD"); + curvstime.SetMarkerStyle(7); + curvstime.SetMarkerSize(2); + curvstime.GetYaxis()->SetTitle("#muA"); + curvstime.GetXaxis()->SetTitle("MJD"); + curvstime.GetXaxis()->CenterTitle(); + curvstime.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); + curvstime.DrawClone("AP"); + txt.DrawText(x, max, text); + + + + gROOT->SetSelectedPad(0); + TCanvas &c11 = d->AddTab("Th vs Time", "Th vs Time"); + gPad->SetGridy(); + + thvstime.GetPoint(1,x,y); + max = thvstime.GetHistogram()->GetMaximum(); + min = thvstime.GetHistogram()->GetMinimum(); + + thvstime.SetTitle("Threshold vs MJD"); + thvstime.SetMarkerStyle(7); + thvstime.SetMarkerSize(2); + thvstime.GetYaxis()->SetTitle("DAC counts"); + thvstime.GetXaxis()->SetTitle("MJD"); + thvstime.GetXaxis()->CenterTitle(); + thvstime.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); + thvstime.DrawClone("AP"); + txt.DrawText(x, max, text); + + + + gROOT->SetSelectedPad(0); + TCanvas &c12 = d->AddTab("Exc vs Cur", "Exc vs Cur"); + + excvscur.GetPoint(1,x,y); + max = excvscur.GetHistogram()->GetMaximum(); + min = excvscur.GetHistogram()->GetMinimum(); + + excvscur.SetTitle("Excess Rate vs Currents"); + excvscur.SetMarkerStyle(7); + excvscur.SetMarkerSize(2); + excvscur.GetYaxis()->SetTitle("excess evts/h"); + excvscur.GetXaxis()->SetTitle("#muA"); + excvscur.GetXaxis()->CenterTitle(); + excvscur.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); + excvscur.DrawClone("AP"); + txt.DrawText(x, max, text); + + + + gROOT->SetSelectedPad(0); + TCanvas &c13 = d->AddTab("Bg vs Cur", "Bg vs Cur"); + + curvstime.GetPoint(1,x,y); + max = bgvscur.GetHistogram()->GetMaximum(); + min = bgvscur.GetHistogram()->GetMinimum(); + + bgvscur.SetTitle("Background Rate vs Current"); + bgvscur.SetMarkerStyle(7); + bgvscur.SetMarkerSize(2); + bgvscur.GetYaxis()->SetTitle("background evts/h"); + bgvscur.GetXaxis()->SetTitle("#muA"); + bgvscur.GetXaxis()->CenterTitle(); + bgvscur.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); + bgvscur.DrawClone("AP"); + txt.DrawText(y, max, text); + + + + gROOT->SetSelectedPad(0); + TCanvas &c14 = d->AddTab("Exc vs Th", "Exc vs Th"); + + excvsth.GetPoint(1,x,y); + max = excvsth.GetHistogram()->GetMaximum(); + min = excvsth.GetHistogram()->GetMinimum(); + + excvsth.SetTitle("Excess Rate vs Threshold"); + excvsth.SetMarkerStyle(7); + excvsth.SetMarkerSize(2); + excvsth.GetYaxis()->SetTitle("excess evts/h"); + excvsth.GetXaxis()->SetTitle("DAC counts"); + excvsth.GetXaxis()->CenterTitle(); + excvsth.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); + excvsth.DrawClone("AP"); + txt.DrawText(x, max, text); + + + + gROOT->SetSelectedPad(0); + TCanvas &c15 = d->AddTab("Bg vs Th", "Bg vs Th"); + + thvstime.GetPoint(1,x,y); + max = bgvsth.GetHistogram()->GetMaximum(); + min = bgvsth.GetHistogram()->GetMinimum(); + + bgvsth.SetTitle("Background Rate vs Threshold"); + bgvsth.SetMarkerStyle(7); + bgvsth.SetMarkerSize(2); + bgvsth.GetYaxis()->SetTitle("background evts/h"); + bgvsth.GetXaxis()->SetTitle("DAC counts"); + bgvsth.GetXaxis()->CenterTitle(); + bgvsth.GetYaxis()->SetRangeUser(min-min*0.1,max*1.1); + bgvsth.DrawClone("AP"); + txt.DrawText(y, max, text); + + + d->SaveAs(outfile);