DvGraph Class Reference

Inheritance diagram for DvGraph:

Collaboration diagram for DvGraph:

Public Member Functions
	DvGraph (char *, string, const string &, int, int, bool)

Public Attributes
TGraph *	rtgraph

Detailed Description

Definition at line 1797 of file TRTCalib_makeplots.cxx.

Constructor & Destructor Documentation

DvGraph()

DvGraph::DvGraph	(	char *	infile,
		string	path,
		const string &	folder,
		int	det,
		int	lay,
		bool	isinverted
	)

Definition at line 1805 of file TRTCalib_makeplots.cxx.

{
 
    this->SetName(Form("Dv_%i_%i", det, lay));
 
    vector<string> levels;
 
    while (true)
    {
        if (const auto f = path.find(','); f==std::string::npos ){
            break;
        } else {
          levels.push_back(path.substr(0, f));
          path.erase(f+1,std::string::npos);
        }
    }
    levels.push_back(path.substr(0, path.find(",")));
    cout << " In DvGraph   folder: " << folder << " det " << det << " lay " << lay << endl;
    for (unsigned int i = 0; i < levels.size(); i++)
        cout << " " << levels.at(i);
    cout << endl;
 
    TFile *file = new TFile(infile);
    cout << infile << endl;
    file->cd();
 
    TDirectory *trt = (TDirectory *)file->FindKey("TRT_all")->ReadObj();
    if (folder != "")
    {
        TDirectory *thisdet = (TDirectory *)trt->FindKey(folder.c_str())->ReadObj();
        cout << "PLOT FOR " << folder << endl;
        if (thisdet->FindKey("rtgraph"))
        {
            rtgraph = (TGraphErrors *)thisdet->FindKey("rtgraph")->ReadObj();
            cout << "found rtgraph " << folder << endl;
        }
        else
        {
            rtgraph = nullptr;
            cout << "did not find rtgraph " << folder << endl;
        }
    }
    else if (trt->FindKey("rtgraph"))
    {
        rtgraph = (TGraphErrors *)trt->FindKey("rtgraph")->ReadObj();
        cout << "no folder, but found rtgraph " << folder << endl;
    }
    else
    {
        rtgraph = nullptr;
        cout << "no folder, no rtgraph " << folder << endl;
    }
 
    if (rtgraph)
    {
 
        TF1 *rtfunc = (TF1 *)rtgraph->GetListOfFunctions()->First();
 
        double p0 = rtfunc->GetParameter(0);
        double p1 = rtfunc->GetParameter(1);
        double p2 = rtfunc->GetParameter(2);
        double p3 = rtfunc->GetParameter(3);
 
        cout << "  R-t PARAMETERS: " << p0 << " " << p1 << " " << p2 << " " << p3 << endl;
 
        const bool isdines = false; //hardcoded flag
 
        if (isdines)
        {
            //coverity[dead_error_begin]
            string vrrelation;
            if (isinverted)
                vrrelation = "x*[2]/sqrt([0]*[0]+x*x)";
            else
                vrrelation = "sqrt([0]*[0]+x*x)/(x*[2])";
 
            TF1 *vr1 = new TF1("vr-relation", vrrelation.c_str(), 0.0001, 2.0);
            vr1->SetParameters(p0, p1, p2);
 
            for (int i = 0; i <= 100; i++)
            {
                float r = 0.1 + (2.0 - 0.1) * (float)i / 100;
                float v1 = vr1->Eval(r, 0, 0, 0);
                this->SetPoint(i, r, v1);
                // cout << i << " " << r << " " << v1 << endl;
            }
 
            vr1->Delete();
        }
        else
        {
            cout << "  USING POLYNOMIAL R-t RELATION" << endl;
 
            string rtrelation = "[0]+x*([1]+x*([2]+x*[3]))";
            string vtrelation;
            if (isinverted)
                vtrelation = "1/([0]+2*x*[1]+3*x*x*[2])";
            else
                vtrelation = "[0]+2*x*[1]+3*x*x*[2]";
 
            TF1 *rt1 = new TF1("rt-relation", rtrelation.c_str(), -10000, 10000);
            rt1->SetParameters(p0, p1, p2, p3);
            TF1 *vt1 = new TF1("vt-relation", vtrelation.c_str(), -10000, 10000);
            vt1->SetParameters(p1, p2, p3);
 
            for (int i = 0; i <= 100; i++)
            {
                float r = 0.1 + (2.0 - 0.1) * (float)i / 100;
                float t1 = rt1->GetX(r, -10000, 10000);
                float v1 = vt1->Eval(t1, 0, 0, 0);
                this->SetPoint(i, r, v1);
                // cout << i << " " << r << " " << t1 << " " << v1 << endl;
            }
 
            rt1->Delete();
            vt1->Delete();
        }
    }
}

Member Data Documentation

rtgraph

TGraph* DvGraph::rtgraph

Definition at line 1802 of file TRTCalib_makeplots.cxx.

---

The documentation for this class was generated from the following file:

• TRTCalib_makeplots.cxx