Draw_FieldApprox.C

Go to the documentation of this file.

/*
 *====================================================================
 *
 *  Field approximation quality check
 *
 *  Authors: M.Zyzak
 *
 *  e-mail : M.Zyzak@gsi.de
 *
 *====================================================================
 */
 
#include <unistd.h>   // for dir navigation
 
void Draw_FieldApprox()
{
 
    TString fileName = "FieldApprox.root";
    const int NStations = 8;
 
    TStopwatch timer;
    timer.Start();
    // ------------------------------------------------------------------------
 
    TStyle* plain = new TStyle("Plain", "Plain Style(no colors/fill areas)");
    // plain->SetCanvasBorderMode(0);
    // plain->SetPadBorderMode(0);
    plain->SetPadColor(0);
    plain->SetCanvasColor(0);
    plain->SetTitleColor(0);
    plain->SetStatColor(0);
    plain->SetOptFit(1111);
    // plain->SetStatX(0.7);
    plain->SetStatW(0.225);
    plain->SetOptStat(0);
    plain->cd();
 
    TFile* f = new TFile(fileName.Data(), "read");
 
    TH2F* stB[20];
    TH2F* stBx[20];
    TH2F* stBy[20];
    TH2F* stBz[20];
 
    TCanvas* c[20];
 
    gStyle->SetPalette(1);
 
    FILE* fff;
    fff = fopen("maxErr.txt", "w");
 
    system("mkdir FieldApprox -p");
    chdir("FieldApprox");
    for (int ist = 0; ist < NStations; ist++) {
        TString nameCanv;
 
        //     if(ist<2)
        //       nameCanv = Form("MVD %i",ist+1);
        //     else
        //       nameCanv = Form("STS %i",ist-1);
        nameCanv = Form("STS %i", ist + 1);
 
        c[ist] = new TCanvas(nameCanv.Data(), nameCanv.Data(), 0, 0, 1000, 900);
 
        c[ist]->Divide(2, 2);
        int ff;
 
        cout << Form("station %i, dB, z = %i", ist + 1, ff) << endl;
        stB[ist] = (TH2F*)f->Get(Form("station %i, dB", ist + 1));
        stBx[ist] = (TH2F*)f->Get(Form("station %i, dBx", ist + 1));
        stBy[ist] = (TH2F*)f->Get(Form("station %i, dBy", ist + 1));
        stBz[ist] = (TH2F*)f->Get(Form("station %i, dBz", ist + 1));
 
        TPaletteAxis* palette;
 
        c[ist]->cd(1);
        stB[ist]->GetZaxis()->SetTitleOffset(2.1);
        stB[ist]->SetTitle(Form("z = %i cm", ff));
        stB[ist]->Draw("colz");
 
        gPad->SetRightMargin(0.2);
        gPad->Update();
 
        c[ist]->cd(2);
        stBx[ist]->GetZaxis()->SetTitleOffset(2.1);
        stBx[ist]->SetTitle(Form("z = %i cm", ff));
        stBx[ist]->Draw("colz");
 
        gPad->SetRightMargin(0.2);
        gPad->Update();
 
        c[ist]->cd(3);
        stBy[ist]->GetZaxis()->SetTitleOffset(2.1);
        stBy[ist]->SetTitle(Form("z = %i cm", ff));
        stBy[ist]->Draw("colz");
 
        gPad->SetRightMargin(0.2);
        gPad->Update();
 
        c[ist]->cd(4);
        stBz[ist]->GetZaxis()->SetTitleOffset(2.1);
        stBz[ist]->SetTitle(Form("z = %i cm", ff));
        stBz[ist]->Draw("colz");
 
        gPad->SetRightMargin(0.2);
        gPad->Update();
 
        double XmaxB = fabs(stB[ist]->GetMaximum());
        double XminB = fabs(stB[ist]->GetMinimum());
        if (XmaxB < XminB)
            XmaxB = XminB;
 
        double XmaxBx = fabs(stBx[ist]->GetMaximum());
        double XminBx = fabs(stBx[ist]->GetMinimum());
        if (XmaxBx < XminBx)
            XmaxBx = XminBx;
 
        double XmaxBy = fabs(stBy[ist]->GetMaximum());
        double XminBy = fabs(stBy[ist]->GetMinimum());
        if (XmaxBy < XminBy)
            XmaxBy = XminBy;
 
        double XmaxBz = fabs(stBz[ist]->GetMaximum());
        double XminBz = fabs(stBz[ist]->GetMinimum());
        if (XmaxBz < XminBz)
            XmaxBz = XminBz;
 
        fprintf(fff, "%.2g %.2g  %.2g  %.2g\n", XmaxB, XmaxBx, XmaxBy, XmaxBz);
 
        nameCanv += ".pdf";
        c[ist]->SaveAs(nameCanv.Data());
    }
 
    timer.Stop();
    Double_t rtime = timer.RealTime();
    Double_t ctime = timer.CpuTime();
    printf("RealTime=%f seconds, CpuTime=%f seconds\n", rtime, ctime);
 
    chdir("..");
}