GoodDataSelTimingPeakStart.C

Go to the documentation of this file.

void DrawArrow(double x, double dx, double ymin, double ymax){

        TLine line;
        line.SetLineWidth(3);
        line.SetLineColor(2);
        line.DrawLine(x,ymin,x,ymax);
        TArrow arr;
        arr.SetLineWidth(3);
        arr.SetLineColor(2);
        arr.SetFillColor(2);
        arr.DrawArrow(x,ymax,x+dx,ymax,0.01,"|>");

}

void SetLeg(TLegend *leg){

        leg->SetFillStyle(0);
        leg->SetBorderSize(0);
        leg->SetTextSize(0.04);
        leg->SetTextFont(42);
        leg->SetY1(leg->GetY2()-leg->GetNRows()*0.05);

        return;

}

TLatex *MiscText(float x, float y, float size, TString text)
{

        TLatex *l = new TLatex(x,y,text);
        l->SetNDC();
        l->SetTextSize(size);
        l->SetTextColor(kBlack);
        l->Draw();

        return l;

}

void UpdateStamp(int ptype, int maxt){

        TDatime tnow(maxt);
        TString lasttime = TString::Format("Last updated on: %s", tnow.AsString());
        TLatex *ltx = MiscText(0.6,0.96,0.03,lasttime);
        //  ltx->SetTextAngle(-90);

}

void StripSpecialChar(string &instr){

        cout<<"\n"<<instr<<endl;
        cout<<"\n"<<endl;
        char str[50];
        strcpy(str,instr.c_str());
        // char str[50] = instr;
        //std::string str[50]=instr.c_str();
        char strStripped[50];
        int c = 0;

        for(int i=0; i<strlen(str); i++){

                if (isalnum(str[i])){

                        strStripped[c] = str[i];
                        c++;

                }

        }

        strStripped[c] = '\0';
        instr = strStripped;

}

void GoodDataSelTimingPeakStart(int ptype=0){
        TDatime t;
        int GMToffset = (t.Convert(true) - t.Convert());
        gSystem->Setenv("TZ","UTC");
        TString stdP1 = "1";
        t.Set();
        int maxt = t.Convert();
        int mint = TDatime(20140801,0).Convert();
        //  TDatime(maxt).Print();
        int dt = 3600*24;

        if(ptype==1) dt *= 7;

        if(ptype==0) dt *= 30;

        if(ptype==3) dt = maxt - mint;

        TString tsel = TString::Format("&& lastsec > %d && lastsec < %d",maxt - dt + 2*GMToffset, maxt + 2*GMToffset);
        stdP1 += tsel;
        int tdiv = 706;
        int tbin = 3600*24;
        TString tlab = "Date (dd/mm/yy - Central Time)";
        TString tform = "%d/%m/%y";
        TString suf = "Month";

        if(ptype==1){

                tbin = 3600*6;
                tdiv = 408;
                tlab = "Date (dd/mm - Central Time)";
                tform = "%d/%m";
                suf = "Week";

        }

        if(ptype==2){

                tbin = 3600;
                tdiv = 309;
                tlab = "Central Time (hour:min)";
                tform = "%H:%M";
                suf = "Day";

        }

        if(ptype==3){

                tbin = (maxt - mint)/30;
                tdiv = 706;
                tlab = "Date (dd/mm/yy - Central Time)";
                tform = "%d/%m/%y";
                suf = "All";

        }

        TString metricstr = "TMath::Max(TMath::Min(timingPeakStart * 1e-3,230.24),210.26)"; 
        TString sel = stdP1; 
        TString title = "Timing Peak Start (#mus)";
        int nbins = 40; 
        double xmin = 210.25; 
        double xmax = 230.25;
        TString filenames = "MetricsTree*.root";
        gStyle->SetTimeOffset(0);
        gStyle->SetOptStat(0);

        if(sel=="") sel = "1";

        if(title=="") title = metricstr;

        TString selpre = "((" + sel + ") && preliminary)";
        //  TString selpre = "((" + sel + ") && 0)";
        TChain *mytree = new TChain("mytree");
        mytree->Add(filenames);
        //  MakeLongCanvas();
        //*************************************************
        //  The code for the above function was copied below
        TCanvas *c1 = new TCanvas("c1","c1",1000,600);
        c1->SetBottomMargin(0.12);
        c1->SetLeftMargin(0.12);
        c1->SetRightMargin(0.07);
        c1->SetTopMargin(0.07);
        //*************************************************
        TLegend *leg = new TLegend(0.15,0.65,0.65,0.9);
        TLegend *legpre = new TLegend(0.67,0.65,0.95,0.9);
        mytree->Draw(metricstr,sel,"hist");
        TH1F* hAll = 0;
        hAll = (TH1F*)(mytree->GetHistogram())->Clone("");
        if(xmin == -9999 || xmax == -9999){

                xmin = hAll->GetXaxis()->GetXmin();
                xmax = hAll->GetXaxis()->GetXmax();

        }

        if(nbins<0){

                xmin = floor(xmin);
                xmax = ceil(xmax);
                nbins *= -1;

                if(int(xmax-xmin)%nbins!=0){

                        xmax += (nbins - int(xmax-xmin)%nbins)/2;
                        xmin -= (nbins - int(xmax-xmin)%nbins);

                }

                nbins = (xmax - xmin) / nbins;

        }

        delete hAll;
        hAll =  new TH1F("hAll","hAll",nbins,xmin,xmax);
        mytree->Draw(metricstr+">>hAll",sel,"hist");
        hAll->Clear();
        hAll->Draw("hist");
        hAll->SetLineColor(1);
        hAll->SetLineWidth(2);
        hAll->SetTitle(";"+title+";Subruns");
        hAll->GetXaxis()->SetTitleOffset(1);
        hAll->GetYaxis()->SetTitleOffset(1);
        hAll->GetXaxis()->CenterTitle();
        hAll->GetYaxis()->CenterTitle();
        hAll->GetXaxis()->SetNdivisions();
        hAll->GetYaxis()->SetNdivisions();
        hAll->SetMinimum(0.1);
        hAll->SetStats(kFALSE);
        double max = hAll->GetMaximum();
        hAll->SetMaximum(pow(10,1.5*(log10(max)+1)-1));
        leg->AddEntry(hAll,"Good Subruns","f");

        mytree->Draw(metricstr,"!(pass_all && !preliminary) && "+sel,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        mytree->GetHistogram()->SetFillColor(kOrange+7);
        leg->AddEntry(mytree->GetHistogram(),"Failed Slice","f");

        mytree->Draw(metricstr,"!(pass_db && pass_hits && passemptyspill && pass_timingpeak && pass_time && pass_duration && !preliminary) && "+sel,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        mytree->GetHistogram()->SetFillColor(kGreen+1);
        leg->AddEntry(mytree->GetHistogram(),"Failed Diblock","f");

        mytree->Draw(metricstr,"!(pass_hits && passemptyspill && pass_timingpeak && pass_time && pass_duration && !preliminary) && "+sel,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        mytree->GetHistogram()->SetFillColor(kBlue+2);
        leg->AddEntry(mytree->GetHistogram(),"Failed Hit Rate","f");

        mytree->Draw(metricstr,"!(passemptyspill && pass_timingpeak && pass_time && pass_duration && !preliminary) && "+sel,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        mytree->GetHistogram()->SetFillColor(kBlue-9);
        leg->AddEntry(mytree->GetHistogram(),"Failed Empty Spill","f");

        mytree->Draw(metricstr,"!(pass_timingpeak && pass_time && pass_duration && !preliminary) && "+sel,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        //mytree->GetHistogram()->SetFillColor(kMagenta-8);
        mytree->GetHistogram()->SetFillColor(kRed+1);
        leg->AddEntry(mytree->GetHistogram(),"Failed Timing Peak","f");

        mytree->Draw(metricstr,"!(pass_time && pass_duration && !preliminary) && "+sel,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        mytree->GetHistogram()->SetFillColor(kMagenta+2);
        leg->AddEntry(mytree->GetHistogram(),"Failed NuMI Live Time","f");

        mytree->Draw(metricstr,"!(pass_time && !preliminary) && "+sel,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        mytree->GetHistogram()->SetFillColor(kBlack);
        leg->AddEntry(mytree->GetHistogram(),"Failed Timestamp","f");

        mytree->Draw(metricstr,selpre,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        mytree->GetHistogram()->SetFillColor(18);
        legpre->AddEntry(mytree->GetHistogram()," ","f");

        mytree->Draw(metricstr,"!pass_all && "+selpre,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        mytree->GetHistogram()->SetFillColor(kGreen-9);
        legpre->AddEntry(mytree->GetHistogram()," ","f");

        mytree->Draw(metricstr,"!(pass_time && pass_duration) && "+selpre,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        mytree->GetHistogram()->SetFillColor(kMagenta-7);
        legpre->AddEntry(mytree->GetHistogram()," ","f");

        mytree->Draw(metricstr,"!(pass_time) && "+selpre,"hist same");
        mytree->GetHistogram()->SetLineColor(1);
        mytree->GetHistogram()->SetLineWidth(2);
        mytree->GetHistogram()->SetFillColor(kGray+2);
        legpre->AddEntry(mytree->GetHistogram()," Preliminary","f");
        SetLeg(leg);
        leg->SetNColumns(2);
        leg->SetY1(0.9-leg->GetNRows()*0.05);
        leg->Draw();
        legpre->SetNColumns(2);
        SetLeg(legpre);
        legpre->SetColumnSeparation(-0.05);
        legpre->Draw();
        c1->SetLogy();

        MiscText(0.15,0.94,0.055,"Near Detector");
        UpdateStamp(ptype,maxt);
        DrawArrow(216.75,-1,0.1,0.7*hAll->GetBinContent(hAll->GetMaximumBin()));
        DrawArrow(219.25,1,0.1,0.7*hAll->GetBinContent(hAll->GetMaximumBin()));
        c1->SaveAs("pngfiles/NearDet-t00-P1GoodDataSelTimingPeakStart"+suf+".png");
        c1->Close();
}