//************************************************
// Author: Federica Lionetto
// Created on: 25/02/2015
//************************************************
/*
SCurve contains some useful functions related to s-curves:
- fLeft describes the falling s-curve;
- fRight describes the rising s-curve;
- fSigma describes the sigma as a function of the z position;
- assignParInfoSCurves sets the starting values and the limits of the parameters of the two erf functions describing the s-curves according to the filename of the measurement;
- assignParInfoSigma sets the starting values and the limits of the parameters of the two sigma functions according to the filename of the measurement;
- assignParInfoEta sets the starting values and the limits of the parameters of the erf function describing the charge sharing coefficient according to the filename of the measurement;
- fitSCurve fits the s-curve and returns the result of the fit;
- fitMu fits the mu as a function of the z position and returns the result of the fit.
- fitSigma fits the sigma as a function of the z position and returns the result of the fit.
*/
// Header guard.
#ifndef __SCURVE_C_INCLUDED__
#define __SCURVE_C_INCLUDED__
#include "Lib.C"
#include "lhcbStyle.C"
#include "Style.C"
#include "Par.C"
// Declarations.
Double_t fLeft(Double_t *x, Double_t *par);
Double_t fRight(Double_t *x, Double_t *par);
Double_t fSigma(Double_t *x, Double_t *par);
void assignParInfoSCurves(string filename, Float_t zvalue, Float_t *minRight, Float_t *maxRight, Float_t *minLeft, Float_t *maxLeft, Float_t parSCurves[], Int_t lengthParSCurves);
void assignParInfoSigma(string filename, Float_t parSigma[], Int_t lengthParSigma);
void assignParInfoEta(string filename, Float_t parEta[], Int_t lengthParEta);
void fitSCurve(string edge, TGraph *graph, TF1 *f, Float_t parSCurves[], Int_t lengthParSCurves, Float_t *a, Float_t *mu, Float_t *sigma, Float_t *b, Float_t *ua, Float_t *umu, Float_t *usigma, Float_t *ub);
void fitMu(string edge, TGraph *graph, TF1 *f, Float_t *p0, Float_t *p1, Float_t *up0, Float_t *up1);
void fitSigma(string edge, TGraph *graph, TF1 *f, Float_t parSigma[], Int_t lengthParSigma, Float_t *q0, Float_t *q1, Float_t *q2, Float_t *uq0, Float_t *uq1, Float_t *uq2);
// Definitions.
// fLeft has four parameters:
// - par[0] is the y-spread around the middle of the curve;
// - par[1] is the x-coordinate of the middle of the curve;
// - par[2] is the steepness of the curve;
// - par[3] is the y-coordinate of the middle of the curve.
Double_t fLeft(Double_t *x, Double_t *par)
{
return par[0]*(1.-erf((x[0]-par[1])/(sqrt(2.)*par[2])))+par[3];
}
// fRight has four parameters:
// - par[0] is the y-spread around the middle of the curve;
// - par[1] is the x-coordinate of the middle of the curve;
// - par[2] is the steepness of the curve;
// - par[3] is the y-coordinate of the middle of the curve.
Double_t fRight(Double_t *x, Double_t *par)
{
return par[0]*erf((x[0]-par[1])/(sqrt(2.)*par[2]))+par[3];
}
// fSigma is used to fit the sigma as a function of the z position and has three parameters:
// - par[0] is the angular coefficient of the two straight lines;
// - par[1] is the x-coordinate of the lowest point;
// - par[2] is the y-coordinate of the lowest point.
Double_t fSigma(Double_t *x, Double_t *par)
{
return par[0]*abs(x[0]-par[1])+par[2];
}
// assignParInfoSCurves has the following parameters:
// - string filename, the filename of the measurement;
// - Float_t zvalue, the position along z;
// - Float_t *minRight, the lower limit of the fitting range for the fRight function;
// - Float_t *maxRight, the upper limit of the fitting range for the fRight function;
// - Float_t *minLeft, the lower limit of the fitting range for the fLeft function;
// - Float_t *maxLeft, the upper limit of the fitting range for the fLeft function;
// - Float_t parSCurves[], the array containing the parameters of the s-curves;
// - Float_t lengthSCurves, the length of the array containing the parameters of the s-curves.
void assignParInfoSCurves(string filename, Float_t zvalue, Float_t *minRight, Float_t *maxRight, Float_t *minLeft, Float_t *maxLeft, Float_t parSCurves[], Int_t lengthParSCurves)
{
if (filename == "ProcessRawData-201501093-216") {
// Right.
if (zvalue < 120.)
{
*minRight = 70.;
*maxRight = 180.;
}
else
{
*minRight = 60.;
*maxRight = 140.;
}
// Left.
if (zvalue < 40.)
{
*minLeft = 180.;
*maxLeft = 280.;
}
else if (zvalue < 120.)
{
*minLeft = 180.;
*maxLeft = 270.;
}
else
{
*minLeft = 180.;
*maxLeft = 260.;
}
// Right.
if (zvalue < 120.)
parSCurves[0] = 125.;
else
parSCurves[0] = 140.;
parSCurves[1] = 95.;
parSCurves[2] = 10.;
parSCurves[3] = 200.;
parSCurves[4] = 90.;
parSCurves[5] = 170.;
parSCurves[6] = 70.;
parSCurves[7] = 120.;
parSCurves[8] = 3.;
parSCurves[9] = 20.;
parSCurves[10] = 150.;
parSCurves[11] = 220.;
// Left.
if (zvalue < 120.)
parSCurves[12] = 125.;
else
parSCurves[12] = 140.;
parSCurves[13] = 230.;
parSCurves[14] = 10.;
parSCurves[15] = 200.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = 215.;
parSCurves[19] = 250.;
parSCurves[20] = 3.;
parSCurves[21] = 20.;
// fitLeft->SetParLimits(3,150.,220.);
}
// It does not look very nice.
else if (filename == "ProcessRawData-20150111-216") {
// Right.
if (zvalue < 105.)
{
*minRight = 60.;
*maxRight = 120.;
}
else
{
*minRight = 50.;
*maxRight = 120.;
}
// Left.
*minLeft = 180.;
*maxLeft = 250.;
// Right.
parSCurves[0] = 140.;
parSCurves[1] = 85.;
parSCurves[2] = 6.;
parSCurves[3] = 175.;
parSCurves[4] = 110.;
parSCurves[5] = 170.;
parSCurves[6] = 70.;
parSCurves[7] = 100.;
parSCurves[8] = 3.;
parSCurves[9] = 10.;
parSCurves[10] = 150.;
parSCurves[11] = 200.;
// Left.
parSCurves[12] = 140.;
parSCurves[13] = 220.;
parSCurves[14] = 6.;
parSCurves[15] = 175.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = 210.;
parSCurves[19] = 230.;
parSCurves[20] = 3.;
parSCurves[21] = 10.;
// fitLeft->SetParLimits(3,150.,220.);
}
else if (filename == "ProcessRawData-201501124-216") {
// Right.
if (zvalue < 80.)
{
*minRight = -8440.;
*maxRight = -8350.;
}
else if (zvalue < 180.)
{
*minRight = -8450.;
*maxRight = -8370.;
}
else
{
*minRight = -8460.;
*maxRight = -8370.;
}
// Left.
if (zvalue == 100.)
{
*minLeft = -8320.;
*maxLeft = -8240.;
}
else if (zvalue < 120.)
{
*minLeft = -8350.;
*maxLeft = -8240.;
}
else if (zvalue < 180.)
{
*minLeft = -8350.;
*maxLeft = -8250.;
}
else
{
*minLeft = -8360.;
*maxLeft = -8260.;
}
// Right.
if (zvalue < 100.)
parSCurves[0] = 100.;
else
parSCurves[0] = 130.;
if (zvalue < 100.)
parSCurves[1] = -8400.;
else
parSCurves[1] = -8420.;
parSCurves[2] = 10.;
if (zvalue < 100.)
parSCurves[3] = 200.;
else
parSCurves[3] = 170.;
parSCurves[4] = 80.;
parSCurves[5] = 140.;
parSCurves[6] = -8450.;
parSCurves[7] = -8370.;
parSCurves[8] = 3.;
parSCurves[9] = 20.;
parSCurves[10] = 150.;
parSCurves[11] = 250.;
// Left.
if (zvalue < 100.)
parSCurves[12] = 100.;
else
parSCurves[12] = 135.;
parSCurves[13] = -8280.;
parSCurves[14] = 10.;
if (zvalue < 100.)
parSCurves[15] = 200.;
else
parSCurves[15] = 180.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = -8300.;
parSCurves[19] = -8250.;
parSCurves[20] = 3.;
parSCurves[21] = 20.;
// fitLeft->SetParLimits(3,150.,220.);
}
else if (filename == "ProcessRawData-20150221-216") {
// Right.
if (zvalue < 180.)
{
*minRight = 60.;
*maxRight = 160.;
}
else
{
*minRight = 50.;
*maxRight = 160.;
}
// Left.
if (zvalue < 180.)
{
*minLeft = 160.;
*maxLeft = 260.;
}
else
{
*minLeft = 160.;
*maxLeft = 250.;
}
// Right.
parSCurves[0] = 15.;
parSCurves[1] = 100.;
parSCurves[2] = 6.;
parSCurves[3] = 12.;
parSCurves[4] = -999.;
parSCurves[5] = -999.;
parSCurves[6] = 70.;
parSCurves[7] = 110.;
parSCurves[8] = 3.;
parSCurves[9] = 15.;
parSCurves[10] = -999.;
parSCurves[11] = -999.;
// Left.
parSCurves[12] = 17.;
parSCurves[13] = 220.;
parSCurves[14] = 6.;
parSCurves[15] = 10.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = 200.;
parSCurves[19] = 230.;
parSCurves[20] = 3.;
parSCurves[21] = 15.;
// fitLeft->SetParLimits(3,150.,220.);
}
else if (filename == "ProcessRawData-201502235-216") {
// Right.
if (zvalue < 80.)
{
*minRight = 120.;
*maxRight = 240.;
}
else
{
*minRight = 100.;
*maxRight = 230.;
}
// Left.
if (zvalue < 80.)
{
*minLeft = 240.;
*maxLeft = 330.;
}
else
{
*minLeft = 230.;
*maxLeft = 320.;
}
// Right.
parSCurves[0] = 16.;
parSCurves[1] = 160.;
parSCurves[2] = 6.;
if (zvalue < 60.)
parSCurves[3] = 12.;
else
parSCurves[3] = 10.;
parSCurves[4] = -999.;
parSCurves[5] = -999.;
parSCurves[6] = 120.;
parSCurves[7] = 200.;
parSCurves[8] = 3.;
parSCurves[9] = 15.;
parSCurves[10] = -999.;
parSCurves[11] = -999.;
// Left.
parSCurves[12] = 16.;
parSCurves[13] = 280.;
parSCurves[14] = 6.;
parSCurves[15] = 12.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = 260.;
parSCurves[19] = 300.;
parSCurves[20] = 3.;
parSCurves[21] = 15.;
// fitLeft->SetParLimits(3,150.,220.);
}
else if ((filename == "ProcessRawData-201502253-216") || (filename == "ProcessRawData-20150226-216") || (filename == "ProcessRawData-201502262-216") || (filename == "ProcessRawData-20150226-100V-216") || (filename == "ProcessRawData-20150226-120V-216") || (filename == "ProcessRawData-20150226-140V-216") || (filename == "ProcessRawData-20150226-160V-216") || (filename == "ProcessRawData-20150226-180V-216") || (filename == "ProcessRawData-20150226-200V-216") || (filename == "ProcessRawData-20150227-21.0dB-216") || (filename == "ProcessRawData-20150227-21.5dB-216") || (filename == "ProcessRawData-20150227-22.0dB-216") || (filename == "ProcessRawData-20150227-22.5dB-216") || (filename == "ProcessRawData-20150227-23.0dB-216") || (filename == "ProcessRawData-20150227-23.5dB-216") || (filename == "ProcessRawData-20150227-24.0dB-216") || (filename == "ProcessRawData-20150227-24.5dB-216") || (filename == "ProcessRawData-20150227-25.0dB-216") || (filename == "ProcessRawData-20150309-67ns-216") || (filename == "ProcessRawData-20150309-68ns-216") || (filename == "ProcessRawData-20150309-69ns-216") ||(filename == "ProcessRawData-20150309-70ns-216") || (filename == "ProcessRawData-20150309-71ns-216") || (filename == "ProcessRawData-20150309-72ns-216") || (filename == "ProcessRawData-20150309-73ns-216") || (filename == "ProcessRawData-20150309-74ns-216") || (filename == "ProcessRawData-20150309-75ns-216") || (filename == "ProcessRawData-20150315-67ns-216") || (filename == "ProcessRawData-20150315-68ns-216") || (filename == "ProcessRawData-20150315-69ns-216") ||(filename == "ProcessRawData-20150315-70ns-216") || (filename == "ProcessRawData-20150315-71ns-216") || (filename == "ProcessRawData-20150315-72ns-216") || (filename == "ProcessRawData-20150315-73ns-216") || (filename == "ProcessRawData-20150315-74ns-216") || (filename == "ProcessRawData-20150315-75ns-216") || (filename == "ProcessRawData-20150317-216")) {
// Right.
*minRight = 30.;
*maxRight = 140.;
// Left.
if (zvalue < 150.) {
*minLeft = 140.;
*maxLeft = 235;
}
else if (zvalue < 230.) {
*minLeft = 140.;
*maxLeft = 240.;
}
else {
*minLeft = 140.;
*maxLeft = 230.;
}
// Right.
if (zvalue < 130.)
parSCurves[0] = 11.;
else
parSCurves[0] = 16.;
if (zvalue < 110.)
parSCurves[1] = 80.;
else if (zvalue < 150.)
parSCurves[1] = 70.;
else
parSCurves[1] = 60.;
if (zvalue<130.)
parSCurves[2] = 10.;
else
parSCurves[2] = 6.;
parSCurves[3] = 10.;
parSCurves[4] = -999.;
parSCurves[5] = -999.;
parSCurves[6] = 20.;
parSCurves[7] = 100.;
parSCurves[8] = 3.;
parSCurves[9] = 20.;
parSCurves[10] = -999.;
parSCurves[11] = -999.;
// Left.
if (zvalue < 130.)
parSCurves[12] = 11.;
else
parSCurves[12] = 16.;
parSCurves[13] = 190.;
if (zvalue < 130.)
parSCurves[14] = 10.;
else
parSCurves[14] = 6.;
parSCurves[15] = 10.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = 170.;
parSCurves[19] = 210.;
parSCurves[20] = 3.;
parSCurves[21] = 20.;
// fitLeft->SetParLimits(3,150.,220.);
}
else if ((filename == "ProcessRawData-201503232-216") || (filename == "ProcessRawData-20150325-216") || (filename == "ProcessRawData-201503252-67ns-216") || (filename == "ProcessRawData-201503252-68ns-216") || (filename == "ProcessRawData-201503252-69ns-216") || (filename == "ProcessRawData-201503252-70ns-216") || (filename == "ProcessRawData-201503252-71ns-216") || (filename == "ProcessRawData-201503252-72ns-216") || (filename == "ProcessRawData-201503252-73ns-216") || (filename == "ProcessRawData-201503252-74ns-216") || (filename == "ProcessRawData-201503252-75ns-216") || (filename == "ProcessRawData-201503253-150V-216") || (filename == "ProcessRawData-201503253-170V-216") || (filename == "ProcessRawData-201503253-190V-216") || (filename == "ProcessRawData-201503253-210V-216") || (filename == "ProcessRawData-201503253-230V-216") || (filename == "ProcessRawData-201503253-250V-216") || (filename == "ProcessRawData-20150326-24.5dB-216") || (filename == "ProcessRawData-20150326-25.0dB-216") || (filename == "ProcessRawData-20150326-25.5dB-216") || (filename == "ProcessRawData-20150326-26.0dB-216") || (filename == "ProcessRawData-20150326-26.5dB-216") || (filename == "ProcessRawData-20150326-27.0dB-216") || (filename == "ProcessRawData-20150326-27.5dB-216") || (filename == "ProcessRawData-20150326-28.0dB-216") || (filename == "ProcessRawData-20150326-28.5dB-216")) {
// Right.
if (zvalue < 150.) {
*minRight = 120.;
*maxRight = 200.;
}
else {
*minRight = 110.;
*maxRight = 200.;
}
// Left.
if (zvalue < 170.) {
*minLeft = 240.;
*maxLeft = 320.;
}
else if (zvalue < 250.) {
*minLeft = 220.;
*maxLeft = 310.;
}
else {
*minLeft = 220.;
*maxLeft = 300.;
}
// Right.
parSCurves[0] = 2.5;
parSCurves[1] = 150.;
parSCurves[2] = 10.;
parSCurves[3] = 3.;
parSCurves[4] = -999.;
parSCurves[5] = -999.;
parSCurves[6] = 120.;
parSCurves[7] = 180.;
parSCurves[8] = 3.;
parSCurves[9] = 20.;
parSCurves[10] = -999.;
parSCurves[11] = -999.;
// Left.
parSCurves[12] = 2.5;
parSCurves[13] = 275.;
parSCurves[14] = 10.;
parSCurves[15] = 3.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = 250.;
parSCurves[19] = 300.;
parSCurves[20] = 3.;
parSCurves[21] = 20.;
// fitLeft->SetParLimits(3,150.,220.);
}
else if (filename == "ProcessRawData-201503273-216"){
// Right.
if (zvalue < 210.) {
*minRight = 90.;
*maxRight = 170.;
}
else {
*minRight = 80.;
*maxRight = 170.;
}
// Left.
if (zvalue < 90.) {
*minLeft = 200.;
*maxLeft = 300.;
}
else if (zvalue < 210.) {
*minLeft = 200.;
*maxLeft = 290.;
}
else {
*minLeft = 200.;
*maxLeft = 280.;
}
// Right.
parSCurves[0] = 2.5;
parSCurves[1] = 120.;
parSCurves[2] = 10.;
parSCurves[3] = 3.;
parSCurves[4] = -999.;
parSCurves[5] = -999.;
parSCurves[6] = 100.;
parSCurves[7] = 150.;
parSCurves[8] = 3.;
parSCurves[9] = 20.;
parSCurves[10] = -999.;
parSCurves[11] = -999.;
// Left.
parSCurves[12] = 2.5;
parSCurves[13] = 275.;
parSCurves[14] = 10.;
parSCurves[15] = 3.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = 250.;
parSCurves[19] = 300.;
parSCurves[20] = 3.;
parSCurves[21] = 20.;
// fitLeft->SetParLimits(3,150.,220.);
}
else if (filename == "ProcessRawData-201507248-216"){
// Right.
if (zvalue < -20.) {
*minRight = 35.;
*maxRight = 75.;
}
else if (zvalue < 0.) {
*minRight = 40.;
*maxRight = 70.;
}
else if (zvalue < 40.) {
*minRight = 30.;
*maxRight = 70.;
}
else {
*minRight = 30.;
*maxRight = 65.;
}
// Left.
if (zvalue < -40.) {
*minLeft = 0.;
*maxLeft = 40.;
}
else if (zvalue < 0.) {
*minLeft = -5.;
*maxLeft = 35.;
}
else if (zvalue < 20.) {
*minLeft = -5.;
*maxLeft = 30.;
}
else if (zvalue < 80.) {
*minLeft = -10.;
*maxLeft = 31.;
}
else {
*minLeft = -15.;
*maxLeft = 35.;
}
// Right.
parSCurves[0] = 15;
parSCurves[1] = 60.;
parSCurves[2] = 10.;
parSCurves[3] = -45.;
parSCurves[4] = -999.;
parSCurves[5] = -999.;
parSCurves[6] = 50.;
parSCurves[7] = 60.;
parSCurves[8] = 3.;
parSCurves[9] = 20.;
parSCurves[10] = -999.;
parSCurves[11] = -999.;
// Left.
parSCurves[12] = 20.; // 15.
parSCurves[13] = 5.;
parSCurves[14] = 10.;
if (zvalue < 80.)
parSCurves[15] = -45.;
else
parSCurves[15] = -53.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = -10.;
parSCurves[19] = 20.;
parSCurves[20] = 3.;
parSCurves[21] = 20.;
// fitLeft->SetParLimits(3,150.,220.);
}
else if (filename == "ProcessRawData-20150725-216"){
// Right.
if (zvalue < -70.) {
*minRight = 25.;
*maxRight = 65.;
}
else if (zvalue < -50.) {
*minRight = 20.;
*maxRight = 60.;
}
else if (zvalue < -40.) {
*minRight = 15.;
*maxRight = 55.;
}
else if (zvalue < -30.) {
*minRight = 20.;
*maxRight = 55.;
}
else if (zvalue < -20.) {
*minRight = 25.;
*maxRight = 55.;
}
else if (zvalue < 50.) {
*minRight = 30.;
*maxRight = 55.;
}
else if (zvalue < 60.) {
*minRight = 25.;
*maxRight = 60.;
}
else if (zvalue < 80.) {
*minRight = 20.;
*maxRight = 60.;
}
else {
*minRight = 25.;
*maxRight = 55.;
}
// Left.
if (zvalue < -70.) {
*minLeft = -10.;
*maxLeft = 25.;
}
else if (zvalue < -50.) {
*minLeft = -15.;
*maxLeft = 20.;
}
else if (zvalue < -40.) {
*minLeft = -15.;
*maxLeft = 15.;
}
else if (zvalue < 40.) {
*minLeft = -15.;
*maxLeft = 10.;
}
else if (zvalue < 60.) {
*minLeft = -15.;
*maxLeft = 11.;
}
else if (zvalue < 80.) {
*minLeft = -15.;
*maxLeft = 30.;
}
else {
*minLeft = -15.;
*maxLeft = 25.;
}
// Right.
parSCurves[0] = 6.;
parSCurves[1] = 50.;
parSCurves[2] = 10.;
parSCurves[3] = -29.;
parSCurves[4] = -999.;
parSCurves[5] = -999.;
parSCurves[6] = 30.;
parSCurves[7] = 60.;
parSCurves[8] = 3.;
parSCurves[9] = 20.;
parSCurves[10] = -999.;
parSCurves[11] = -999.;
// Left.
parSCurves[12] = 6.;
parSCurves[13] = 5.;
parSCurves[14] = 10.;
parSCurves[15] = -29.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = -10.;
parSCurves[19] = 20.;
parSCurves[20] = 3.;
parSCurves[21] = 20.;
// fitLeft->SetParLimits(3,150.,220.);
}
else if ((filename == "ProcessRawData-201507252-216") || (filename == "ProcessRawData-201507253-56ns-216") || (filename == "ProcessRawData-201507253-57ns-216") || (filename == "ProcessRawData-201507253-58ns-216") || (filename == "ProcessRawData-201507253-59ns-216") || (filename == "ProcessRawData-201507253-60ns-216") || (filename == "ProcessRawData-201507253-61ns-216") || (filename == "ProcessRawData-201507253-62ns-216") || (filename == "ProcessRawData-201507253-63ns-216") || (filename == "ProcessRawData-201507253-64ns-216") || (filename == "ProcessRawData-20150726--250V-216") || (filename == "ProcessRawData-20150726--230V-216") || (filename == "ProcessRawData-20150726--210V-216") || (filename == "ProcessRawData-20150726--190V-216") || (filename == "ProcessRawData-20150726--170V-216") || (filename == "ProcessRawData-20150726--150V-216") || (filename == "ProcessRawData-201507262-24.5dB-216") || (filename == "ProcessRawData-201507262-25.0dB-216") || (filename == "ProcessRawData-201507262-25.5dB-216") || (filename == "ProcessRawData-201507262-26.0dB-216") || (filename == "ProcessRawData-201507262-26.5dB-216") || (filename == "ProcessRawData-201507262-27.0dB-216") || (filename == "ProcessRawData-201507262-27.5dB-216") || (filename == "ProcessRawData-201507262-28.0dB-216") || (filename == "ProcessRawData-201507262-28.5dB-216") || (filename == "ProcessRawData-201507263-56ns-216") || (filename == "ProcessRawData-201507263-57ns-216") || (filename == "ProcessRawData-201507263-58ns-216") || (filename == "ProcessRawData-201507263-59ns-216") || (filename == "ProcessRawData-201507263-60ns-216") || (filename == "ProcessRawData-201507263-61ns-216") || (filename == "ProcessRawData-201507263-62ns-216") || (filename == "ProcessRawData-201507263-63ns-216") || (filename == "ProcessRawData-201507263-64ns-216")){
// Right.
*minRight = 30.;
*maxRight = 55.;
// Left.
*minLeft = -20.;
*maxLeft = 5.;
// Right.
parSCurves[0] = 15.;
parSCurves[1] = 50.;
parSCurves[2] = 10.;
parSCurves[3] = -23.;
parSCurves[4] = -999.;
parSCurves[5] = -999.;
parSCurves[6] = 30.;
parSCurves[7] = 60.;
parSCurves[8] = 2.;
parSCurves[9] = 20.;
parSCurves[10] = -999.;
parSCurves[11] = -999.;
if (filename == "ProcessRawData-201507262-27.5dB-216") {
*maxRight = 60.;
*maxLeft = 10.;
parSCurves[0] = 11.;
parSCurves[1] = 45.;
parSCurves[2] = 4.;
parSCurves[3] = -17.;
}
else if (filename == "ProcessRawData-201507262-28.0dB-216") {
*maxRight = 60.;
*maxLeft = 10.;
parSCurves[0] = 11.;
parSCurves[1] = 49.;
parSCurves[2] = 13.;
parSCurves[3] = -15.;
parSCurves[4] = 10.;
parSCurves[5] = 12.;
parSCurves[6] = 46.;
parSCurves[7] = 50.;
parSCurves[8] = 4.;
parSCurves[9] = 20.;
parSCurves[10] = -20.;
parSCurves[11] = -10.;
}
if (filename == "ProcessRawData-201507263-58ns-216") {
parSCurves[1] = 45.;
parSCurves[2] = 4.;
parSCurves[4] = 10.;
parSCurves[5] = 18.;
parSCurves[10] = -25.;
parSCurves[11] = -15.;
}
// Left.
parSCurves[12] = 15.;
parSCurves[13] = -10.;
parSCurves[14] = 10.;
parSCurves[15] = -23.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = -20.;
parSCurves[19] = 10.;
parSCurves[20] = 2.;
parSCurves[21] = 20.;
// fitLeft->SetParLimits(3,150.,220.);
/*
if ((filename == "ProcessRawData-201507262-27.5dB-216") || (filename == "ProcessRawData-201507262-28.0dB-216")) {
parSCurves[12] = 11.;
parSCurves[13] = -10.;
parSCurves[14] = 4.;
parSCurves[15] = -17.;
// parSCurves[18] = -7.;
// parSCurves[19] = 1.;
}
*/
}
else if (filename == "ProcessRawData-201507273-216"){
// Right.
// z = -80, -70, -60, -50, -40, -30, -20, -10, 0
if (zvalue < 10.) {
*minRight = 95.;
*maxRight = 135.;
}
// z = 10, 20, 30, 40, 50, 60, 70, 80
else {
*minRight = 95.;
*maxRight = 130.;
}
// Left.
// z = -80, -70, -60, -50, -40
if (zvalue < -30.) {
*minLeft = 60.;
*maxLeft = 95.;
}
// z = -30, -20, -10, 0, 10, 20
else if (zvalue < 30.) {
*minLeft = 55.;
*maxLeft = 95.;
}
// z = 30, 40, 50, 60
else if (zvalue < 70.) {
*minLeft = 55.;
*maxLeft = 90.;
}
// z = 70, 80
else {
*minLeft = 55.;
*maxLeft = 100.;
}
// Right.
parSCurves[0] = 14.;
parSCurves[1] = 125.;
parSCurves[2] = 10.;
parSCurves[3] = -35.;
parSCurves[4] = -999.;
parSCurves[5] = -999.;
parSCurves[6] = 110.;
parSCurves[7] = 140.;
parSCurves[8] = 3.;
parSCurves[9] = 20.;
parSCurves[10] = -999.;
parSCurves[11] = -999.;
// Left.
parSCurves[12] = 14.;
parSCurves[13] = 70.;
parSCurves[14] = 10.;
parSCurves[15] = -35.;
// fitLeft->SetParLimits(0,90.,170.);
parSCurves[18] = 50.;
parSCurves[19] = 80.;
parSCurves[20] = 3.;
parSCurves[21] = 20.;
// fitLeft->SetParLimits(3,150.,220.);
if (zvalue == 30.) {
parSCurves[12] = 20.;
parSCurves[13] = 65.;
parSCurves[14] = 10.;
parSCurves[15] = -30.;
}
}
return;
}
// assignParInfoSigma has the following parameters:
// - string filename, the filename of the measurement;
// - Float_t parSigma[], the array containing the parameters of the sigma function;
// - Float_t lengthParSigma, the length of the array containing the parameters of the sigma function.
void assignParInfoSigma(string filename, Float_t parSigma[], Int_t lengthParSigma)
{
if ((filename == "ProcessRawData-201501093-216") || (filename == "ProcessRawData-20150111-216") || (filename == "ProcessRawData-201501124-216") || (filename == "ProcessRawData-20150221-216")) {
// Sigma left.
parSigma[0] = 5.;
parSigma[1] = 600.;
parSigma[2] = 0.025;
parSigma[3] = 0.;
parSigma[4] = 10.;
parSigma[5] = 500.;
parSigma[6] = 700.;
// parSigma[7] = ?.;
// parSigma[8] = ?.;
// Sigma right.
parSigma[9] = 5.;
parSigma[10] = 600.;
parSigma[11] = 0.025;
parSigma[12] = 0.;
parSigma[13] = 10.;
parSigma[14] = 500.;
parSigma[15] = 700.;
// parSigma[16] = ?.;
// parSigma[17] = ?.;
}
else if (filename == "ProcessRawData-201502235-216") {
// Sigma left.
parSigma[0] = 5.;
parSigma[1] = 100.;
parSigma[2] = 6.;
parSigma[3] = 0.;
parSigma[4] = 10.;
parSigma[5] = 0.;
parSigma[6] = 200.;
// parSigma[7] = ?.;
// parSigma[8] = ?.;
// Sigma right.
parSigma[9] = 5.;
parSigma[10] = 100.;
parSigma[11] = 6.;
parSigma[12] = 0.;
parSigma[13] = 10.;
parSigma[14] = 0.;
parSigma[15] = 200.;
// parSigma[16] = ?.;
// parSigma[17] = ?.;
}
else if ((filename == "ProcessRawData-201502253-216") || (filename == "ProcessRawData-20150226-216") || (filename == "ProcessRawData-20150317-216")) {
// Sigma left.
parSigma[0] = 0.01;
parSigma[1] = 850.;
parSigma[2] = 5.;
parSigma[3] = 0.001;
parSigma[4] = 1.;
parSigma[5] = 800.;
parSigma[6] = 900.;
// parSigma[7] = ?.;
// parSigma[8] = ?.;
// Sigma right.
parSigma[9] = 0.01;
parSigma[10] = 850.;
parSigma[11] = 5.;
parSigma[12] = 0.001;
parSigma[13] = 1.;
parSigma[14] = 800.;
parSigma[15] = 900.;
// parSigma[16] = ?.;
// parSigma[17] = ?.;
}
else if ((filename == "ProcessRawData-201503232-216") || (filename == "ProcessRawData-201503273-216")) {
// Sigma left.
parSigma[0] = 0.01;
parSigma[1] = 550.;
parSigma[2] = 5.;
parSigma[3] = 0.001;
parSigma[4] = 1.;
parSigma[5] = 500.;
parSigma[6] = 700.;
// parSigma[7] = ?.;
// parSigma[8] = ?.;
// Sigma right.
parSigma[9] = 0.01;
parSigma[10] = 580.;
parSigma[11] = 5.;
parSigma[12] = 0.001;
parSigma[13] = 1.;
parSigma[14] = 500.;
parSigma[15] = 700.;
// parSigma[16] = ?.;
// parSigma[17] = ?.;
}
else if (filename == "ProcessRawData-201507248-216") {
// Sigma left.
parSigma[0] = 0.05;
parSigma[1] = 100.;
parSigma[2] = 5.;
parSigma[3] = 0.;
parSigma[4] = 10.;
parSigma[5] = 0.;
parSigma[6] = 200.;
// parSigma[7] = ?.;
// parSigma[8] = ?.;
// Sigma right.
parSigma[9] = 0.02;
parSigma[10] = 100.;
parSigma[11] = 3.;
parSigma[12] = 0.;
parSigma[13] = 10.;
parSigma[14] = 0.;
parSigma[15] = 200.;
// parSigma[16] = ?.;
// parSigma[17] = ?.;
}
else if (filename == "ProcessRawData-20150725-216") {
// Sigma left.
parSigma[0] = 0.03;
parSigma[1] = 100.;
parSigma[2] = 2.;
parSigma[3] = 0.;
parSigma[4] = 10.;
parSigma[5] = 0.;
parSigma[6] = 200.;
parSigma[7] = 1.;
parSigma[8] = 4.;
// Sigma right.
parSigma[9] = 0.02;
parSigma[10] = 100.;
parSigma[11] = 2.;
parSigma[12] = 0.;
parSigma[13] = 10.;
parSigma[14] = 0.;
parSigma[15] = 200.;
parSigma[16] = 1.;
parSigma[17] = 4.;
}
else if (filename == "ProcessRawData-201507273-216") {
// Sigma left.
parSigma[0] = 0.02;
parSigma[1] = 70.;
parSigma[2] = 3.;
parSigma[3] = 0.;
parSigma[4] = 2.;
parSigma[5] = 0.;
parSigma[6] = 200.;
parSigma[7] = 1.;
parSigma[8] = 5.;
// Sigma right.
parSigma[9] = 0.02;
parSigma[10] = 60.;
parSigma[11] = 3.2;
parSigma[12] = 0.;
parSigma[13] = 2.;
parSigma[14] = 0.;
parSigma[15] = 200.;
parSigma[16] = 1.;
parSigma[17] = 5.;
}
return;
}
// assignParInfoEta has the following parameters:
// - string filename, the filename of the measurement;
// - Float_t parEta[], the array containing the parameters of the erf function;
// - Float_t lengthParEta, the length of the array containing the parameters of the erf function.
void assignParInfoEta(string filename, Float_t parEta[], Int_t lengthParEta) {
if ((filename == "ProcessRawData-201502262-216") || (filename == "ProcessRawData-20150226-100V-216") || (filename == "ProcessRawData-20150226-120V-216") || (filename == "ProcessRawData-20150226-140V-216") || (filename == "ProcessRawData-20150226-160V-216") || (filename == "ProcessRawData-20150226-180V-216") || (filename == "ProcessRawData-20150226-200V-216") || (filename == "ProcessRawData-20150227-21.0dB-216") || (filename == "ProcessRawData-20150227-21.5dB-216") || (filename == "ProcessRawData-20150227-22.0dB-216") || (filename == "ProcessRawData-20150227-22.5dB-216") || (filename == "ProcessRawData-20150227-23.0dB-216") || (filename == "ProcessRawData-20150227-23.5dB-216") || (filename == "ProcessRawData-20150227-24.0dB-216") || (filename == "ProcessRawData-20150227-24.5dB-216") || (filename == "ProcessRawData-20150227-25.0dB-216") || (filename == "ProcessRawData-20150309-67ns-216") || (filename == "ProcessRawData-20150309-68ns-216") || (filename == "ProcessRawData-20150309-69ns-216") ||(filename == "ProcessRawData-20150309-70ns-216") || (filename == "ProcessRawData-20150309-71ns-216") || (filename == "ProcessRawData-20150309-72ns-216") || (filename == "ProcessRawData-20150309-73ns-216") || (filename == "ProcessRawData-20150309-74ns-216") || (filename == "ProcessRawData-20150309-75ns-216") || (filename == "ProcessRawData-20150315-67ns-216") || (filename == "ProcessRawData-20150315-68ns-216") || (filename == "ProcessRawData-20150315-69ns-216") ||(filename == "ProcessRawData-20150315-70ns-216") || (filename == "ProcessRawData-20150315-71ns-216") || (filename == "ProcessRawData-20150315-72ns-216") || (filename == "ProcessRawData-20150315-73ns-216") || (filename == "ProcessRawData-20150315-74ns-216") || (filename == "ProcessRawData-20150315-75ns-216")) {
parEta[0] = 1.5;
parEta[1] = 140.;
parEta[2] = 5.;
parEta[3] = 0.;
/*
parEta[4] = ;
parEta[5] = ;
parEta[6] = ;
parEta[7] = ;
parEta[8] = ;
parEta[9] = ;
parEta[10] = ;
parEta[11] = ;
*/
}
else if ((filename == "ProcessRawData-20150325-216") || (filename == "ProcessRawData-201503252-67ns-216") || (filename == "ProcessRawData-201503252-68ns-216") || (filename == "ProcessRawData-201503252-69ns-216") || (filename == "ProcessRawData-201503252-70ns-216") || (filename == "ProcessRawData-201503252-71ns-216") || (filename == "ProcessRawData-201503252-72ns-216") || (filename == "ProcessRawData-201503252-73ns-216") || (filename == "ProcessRawData-201503252-74ns-216") || (filename == "ProcessRawData-201503252-75ns-216")|| (filename == "ProcessRawData-201503253-150V-216") || (filename == "ProcessRawData-201503253-170V-216") || (filename == "ProcessRawData-201503253-190V-216") || (filename == "ProcessRawData-201503253-210V-216") || (filename == "ProcessRawData-201503253-230V-216") || (filename == "ProcessRawData-201503253-250V-216") || (filename == "ProcessRawData-20150326-24.5dB-216") || (filename == "ProcessRawData-20150326-25.0dB-216") || (filename == "ProcessRawData-20150326-25.5dB-216") || (filename == "ProcessRawData-20150326-26.0dB-216") || (filename == "ProcessRawData-20150326-26.5dB-216") || (filename == "ProcessRawData-20150326-27.0dB-216") || (filename == "ProcessRawData-20150326-27.5dB-216") || (filename == "ProcessRawData-20150326-28.0dB-216") || (filename == "ProcessRawData-20150326-28.5dB-216")) {
parEta[0] = 1.;
parEta[1] = 200.;
parEta[2] = 6.;
parEta[3] = 0.;
// parEta[4] = 0.5;
// parEta[5] = 2.5;
parEta[6] = 180.;
parEta[7] = 220.;
parEta[8] = 2.;
parEta[9] = 10.;
// parEta[10] = -0.2;
// parEta[11] = 0.2;
}
else if ((filename == "ProcessRawData-201507252-216") || (filename == "ProcessRawData-201507253-56ns-216") || (filename == "ProcessRawData-201507253-57ns-216") || (filename == "ProcessRawData-201507253-58ns-216") || (filename == "ProcessRawData-201507253-59ns-216") || (filename == "ProcessRawData-201507253-60ns-216") || (filename == "ProcessRawData-201507253-61ns-216") || (filename == "ProcessRawData-201507253-62ns-216") || (filename == "ProcessRawData-201507253-63ns-216") || (filename == "ProcessRawData-201507253-64ns-216") || (filename == "ProcessRawData-20150726--250V-216") || (filename == "ProcessRawData-20150726--230V-216") || (filename == "ProcessRawData-20150726--210V-216") || (filename == "ProcessRawData-20150726--190V-216") || (filename == "ProcessRawData-20150726--170V-216") || (filename == "ProcessRawData-20150726--150V-216") || (filename == "ProcessRawData-201507262-24.5dB-216") || (filename == "ProcessRawData-201507262-25.0dB-216") || (filename == "ProcessRawData-201507262-25.5dB-216") || (filename == "ProcessRawData-201507262-26.0dB-216") || (filename == "ProcessRawData-201507262-26.5dB-216") || (filename == "ProcessRawData-201507262-27.0dB-216") || (filename == "ProcessRawData-201507262-27.5dB-216") || (filename == "ProcessRawData-201507262-28.0dB-216") || (filename == "ProcessRawData-201507262-28.5dB-216")) {
parEta[0] = 1.;
parEta[1] = 20.;
parEta[2] = 6.;
parEta[3] = 0.;
// parEta[4] = 0.5;
// parEta[5] = 2.5;
parEta[6] = 0.;
parEta[7] = 30.;
parEta[8] = 2.;
parEta[9] = 10.;
// parEta[10] = -0.2;
// parEta[11] = 0.2;
}
else if ((filename == "ProcessRawData-201507263-56ns-216") || (filename == "ProcessRawData-201507263-57ns-216") || (filename == "ProcessRawData-201507263-58ns-216") || (filename == "ProcessRawData-201507263-59ns-216") || (filename == "ProcessRawData-201507263-60ns-216") || (filename == "ProcessRawData-201507263-61ns-216") || (filename == "ProcessRawData-201507263-62ns-216") || (filename == "ProcessRawData-201507263-63ns-216") || (filename == "ProcessRawData-201507263-64ns-216")) {
parEta[0] = 0.8;
parEta[1] = 19.;
parEta[2] = 6.;
parEta[3] = 0.;
parEta[4] = 0.5;
parEta[5] = 1.5;
parEta[6] = 10.;
parEta[7] = 25.;
parEta[8] = 2.;
parEta[9] = 20.;
parEta[10] = -0.2;
parEta[11] = 0.2;
}
return;
}
void fitSCurve(string edge, TGraph *graph, TF1 *f, Float_t parSCurves[], Int_t lengthParSCurves, Float_t *a, Float_t *mu, Float_t *sigma, Float_t *b, Float_t *ua, Float_t *umu, Float_t *usigma, Float_t *ub) {
// Set parameter names.
f->SetParName(0,"a");
f->SetParName(1,"#mu");
f->SetParName(2,"#sigma");
f->SetParName(3,"b");
// Set parameters and parameter limits.
if (edge == "right") {
f->SetParameter(0,parSCurves[0]);
f->SetParameter(1,parSCurves[1]);
f->SetParameter(2,parSCurves[2]);
f->SetParameter(3,parSCurves[3]);
if ((parSCurves[4] != -999) && (parSCurves[5] != -999))
f->SetParLimits(0,parSCurves[4],parSCurves[5]);
f->SetParLimits(1,parSCurves[6],parSCurves[7]);
f->SetParLimits(2,parSCurves[8],parSCurves[9]);
if ((parSCurves[10] != -999) && (parSCurves[11] != -999))
f->SetParLimits(3,parSCurves[10],parSCurves[11]);
}
else if ((edge == "left") && (lengthParSCurves>12)) {
f->SetParameter(0,parSCurves[12]);
f->SetParameter(1,parSCurves[13]);
f->SetParameter(2,parSCurves[14]);
f->SetParameter(3,parSCurves[15]);
// f->SetParLimits(0,parSCurves[16],parSCurves[17]);
f->SetParLimits(1,parSCurves[18],parSCurves[19]);
f->SetParLimits(2,parSCurves[20],parSCurves[21]);
// f->SetParLimits(3,parSCurves[22],parSCurves[23]);
}
// Fit.
graph->Fit(f,"BR");
// Get parameters and parameter uncertainties.
*a = f->GetParameter(0);
*mu = f->GetParameter(1);
*sigma = f->GetParameter(2);
*b = f->GetParameter(3);
*ua = f->GetParError(0);
*umu = f->GetParError(1);
*usigma = f->GetParError(2);
*ub = f->GetParError(3);
cout << edge << " side, a = " << *a << ", mu = " << *mu << ", sigma = " << *sigma << ", b = " << *b << endl;
return;
}
void fitMu(string edge, TGraph *graph, TF1 *f, Float_t *p0, Float_t *p1, Float_t *up0, Float_t *up1) {
// Set parameter names.
f->SetParName(0,"p_{0}");
f->SetParName(1,"p_{1}");
// Fit.
f->SetLineColor(kRed);
graph->Fit(f,"BR");
// Get parameters and parameter uncertainties.
*p0 = f->GetParameter(0);
*p1 = f->GetParameter(1);
*up0 = f->GetParError(0);
*up1 = f->GetParError(1);
return;
}
void fitSigma(string edge, TGraph *graph, TF1 *f, Float_t parSigma[], Int_t lengthParSigma, Float_t *q0, Float_t *q1, Float_t *q2, Float_t *uq0, Float_t *uq1, Float_t *uq2) {
// Set parameter names.
f->SetParName(0,"q_{0}");
f->SetParName(1,"q_{1}");
f->SetParName(2,"q_{2}");
// Set parameters and parameter limits.
if (edge == "left") {
f->SetParameter(0,parSigma[0]);
f->SetParameter(1,parSigma[1]);
f->SetParameter(2,parSigma[2]);
f->SetParLimits(0,parSigma[3],parSigma[4]);
f->SetParLimits(1,parSigma[5],parSigma[6]);
// f->SetParLimits(2,parSigma[7],parSigma[8]);
}
else if (edge == "right") {
f->SetParameter(0,parSigma[9]);
f->SetParameter(1,parSigma[10]);
f->SetParameter(2,parSigma[11]);
f->SetParLimits(0,parSigma[12],parSigma[13]);
f->SetParLimits(1,parSigma[14],parSigma[15]);
// f->SetParLimits(2,parSigma[16],parSigma[17]);
}
// Fit.
f->SetLineColor(kRed);
graph->Fit(f,"BR");
// Get parameters and parameter uncertainties.
*q0 = f->GetParameter(0);
*q1 = f->GetParameter(1);
*q2 = f->GetParameter(2);
*uq0 = f->GetParError(0);
*uq1 = f->GetParError(1);
*uq2 = f->GetParError(2);
return;
}
#endif
Federica Lionetto