// -------------------------------------------------------------------------
// -----                     CbmMvdStationPar source file              -----
// -----                  Created 28/10/14  by P.Sitzmann              -----
// -------------------------------------------------------------------------

#include "CbmMvdStationPar.h"

#include <iostream>
#include <math.h>

using std::cout;
using std::endl;



// -----   Default constructor   -------------------------------------------
CbmMvdStationPar::CbmMvdStationPar() 
  : TNamed(),
    fStationCount(0),
    fZPositions(), 
    fThicknesses(), 
    fHeights(), 
    fWidths(), 
    fXResolutions(), 
    fYResolutions(), 
    fRadiationLength(), 
    fBeamHeights(), 
    fBeamWidths()
{
}
// -------------------------------------------------------------------------

// -----   Standart constructor   -------------------------------------------
CbmMvdStationPar::CbmMvdStationPar(Int_t count) 
  : TNamed(),
    fStationCount(count),
    fZPositions(), 
    fThicknesses(), 
    fHeights(), 
    fWidths(), 
    fXResolutions(), 
    fYResolutions(), 
    fRadiationLength(), 
    fBeamHeights(), 
    fBeamWidths()
{
}
// -------------------------------------------------------------------------


// -----   Destructor   ----------------------------------------------------
CbmMvdStationPar::~CbmMvdStationPar() {
}
// -------------------------------------------------------------------------

// -----   Public method Print   -------------------------------------------
void CbmMvdStationPar::Print(Option_t* opt) const {

cout 	<< endl << "--I--" << GetName() << "--I--" << endl
	<< endl << "Initialized ParameterFile with " << fStationCount << " Stations" << endl
	<< endl << "-------------------------------------------------------------------------" << endl
	<< endl << "fZPositions[0] " << fZPositions.at(0) << " z position" << endl
	<< endl << "fZPositions[1] " << fZPositions.at(1) << " z position" << endl
	<< endl << "fZPositions[2] " << fZPositions.at(2) << " z position" << endl
        << endl << "fZPositions[3] " << fZPositions.at(3) << " z position" << endl
	<< endl << "-------------------------------------------------------------------------" << endl
        << endl << "fBeamWidths[0] " << fBeamWidths.at(0) << " beam hole halfwidth" << endl
	<< endl << "fBeamWidths[1] " << fBeamWidths.at(1) << " beam hole halfwidth" << endl
	<< endl << "fBeamWidths[2] " << fBeamWidths.at(2) << " beam hole halfwidth" << endl
        << endl << "fBeamWidths[3] " << fBeamWidths.at(3) << " beam hole halfwidth" << endl
	<< endl << "-------------------------------------------------------------------------" << endl
        << endl << "fBeamHeights[0] " << fBeamHeights.at(0) << " beam hole halfheight" << endl
	<< endl << "fBeamHeights[1] " << fBeamHeights.at(1) << " beam hole halfheight" << endl
	<< endl << "fBeamHeights[2] " << fBeamHeights.at(2) << " beam hole halfheight" << endl
        << endl << "fBeamHeights[3] " << fBeamHeights.at(3) << " beam hole halfheight" << endl 
	<< endl << "-------------------------------------------------------------------------" << endl
	<< endl << "fThicknesses[0] " << fThicknesses.at(0) << " station thickness" << endl
	<< endl << "fThicknesses[1] " << fThicknesses.at(1) << " station thickness" << endl
	<< endl << "fThicknesses[2] " << fThicknesses.at(2) << " station thickness" << endl
        << endl << "fThicknesses[3] " << fThicknesses.at(3) << " station thickness" << endl   
	<< endl << "-------------------------------------------------------------------------" << endl
        << endl << "fWidths[0] " << fWidths.at(0) << " station halfwidth" << endl
	<< endl << "fWidths[1] " << fWidths.at(1) << " station halfwidth" << endl
	<< endl << "fWidths[2] " << fWidths.at(2) << " station halfwidth" << endl
        << endl << "fWidths[3] " << fWidths.at(3) << " station halfwidth" << endl
	<< endl << "-------------------------------------------------------------------------" << endl
        << endl << "fHeights[0] " << fHeights.at(0) << " station halfheight" << endl
	<< endl << "fHeights[1] " << fHeights.at(1) << " station halfheight" << endl
	<< endl << "fHeights[2] " << fHeights.at(2) << " station halfheight" << endl
        << endl << "fHeights[3] " << fHeights.at(3) << " station halfheight" << endl   
	<< endl << "--I--" << GetName() << "--I--" << endl;
}
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
Bool_t CbmMvdStationPar::Init()
{
for (Int_t i = 0; i < fStationCount; i++)
	{ fZPositions[i]=0.; 
    	  fThicknesses[i]=0.; 
          fHeights[i]=0.;
          fWidths[i]=0.;
          fXResolutions[i]=0.;
          fYResolutions[i]=0.;
          fRadiationLength[i]=0.; 
          fBeamHeights[i]=0.; 
          fBeamWidths[i]=0.; 
	}
}
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
Double_t CbmMvdStationPar::GetZPosition(Int_t stationNumber)
{
if(stationNumber <= fStationCount)
	{
	return fZPositions[stationNumber];
	}
else
	{
	cout <<  "Station number out of Range " << endl;
	}
return 0.;
}
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
Double_t CbmMvdStationPar::GetThickness(Int_t stationNumber)
{
if(stationNumber <= fStationCount)
	{
	return fThicknesses[stationNumber];
	}
else
	{
	cout <<  "Station number out of Range " << endl;
	}
return 0.;
}

// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
Double_t CbmMvdStationPar::GetHeight(Int_t stationNumber)
{
if(stationNumber <= fStationCount)
	{
	return fHeights[stationNumber];
	}
else
	{
	cout <<  "Station number out of Range " << endl;
	}
return 0.;
}

// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
Double_t CbmMvdStationPar::GetWidth(Int_t stationNumber)
{
if(stationNumber <= fStationCount)
	{
	return fWidths[stationNumber];
	
	}
else
	{
	cout <<  "Station number out of Range " << endl;
	}
return 0.;
}

// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
Double_t CbmMvdStationPar::GetXRes(Int_t stationNumber)
{
if(stationNumber <= fStationCount)
	{
	return fXResolutions[stationNumber]; 
	}
else
	{
	cout <<  "Station number out of Range " << endl;
	}
return 0.;
}

// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
Double_t CbmMvdStationPar::GetYRes(Int_t stationNumber)
{
if(stationNumber <= fStationCount)
	{
	return fYResolutions[stationNumber];
	}
else
	{
	cout <<  "Station number out of Range " << endl;
	}
return 0.;
} 

// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
Double_t CbmMvdStationPar::GetRadLength(Int_t stationNumber)
{
if(stationNumber <= fStationCount)
	{
	return fRadiationLength[stationNumber];
	}
else
	{
	cout <<  "Station number out of Range " << endl;
	}
return 0.;
}

// -------------------------------------------------------------------------

// -------------------------------------------------------------------------  
Double_t CbmMvdStationPar::GetBeamHeight(Int_t stationNumber)
{
if(stationNumber <= fStationCount)
	{
	return fBeamHeights[stationNumber];
	}
else
	{
	cout <<  "Station number out of Range " << endl;
	}
return 0.;
} 

// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
Double_t CbmMvdStationPar::GetBeamWidth(Int_t stationNumber)
{
if(stationNumber <= fStationCount)
	{
	return fBeamWidths[stationNumber];
	}
else
	{
	cout <<  "Station number out of Range " << endl;
	}
return 0.;
}

// -------------------------------------------------------------------------


// -------------------------------------------------------------------------
void CbmMvdStationPar::SetZPosition(Int_t stationNumber, Double_t z)
{
if(fZPositions[stationNumber]!=0)
	{
	fZPositions[stationNumber] = Int_t(((fZPositions[stationNumber] + z)/2)+0.5);
	}
else
	{
	fZPositions[stationNumber] = z;
	}

}
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
void CbmMvdStationPar::SetThickness(Int_t stationNumber, Double_t thickness)
{
Float_t Sensorthickness = fabs(thickness-fZPositions[stationNumber]);
if(Sensorthickness > fThicknesses[stationNumber])
	fThicknesses[stationNumber]=Sensorthickness;
}  
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
void CbmMvdStationPar::SetHeight(Int_t stationNumber, Double_t height)
{
if(fHeights[stationNumber] < height)
	{
	fHeights[stationNumber] = height;
	}

}  
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
void CbmMvdStationPar::SetWidth(Int_t stationNumber, Double_t width)
{
if(fWidths[stationNumber] < width)
	{
	fWidths[stationNumber] = width;
	}
}  
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
void CbmMvdStationPar::SetXRes(Int_t stationNumber, Double_t xres)
{
	fXResolutions[stationNumber] = xres;
}  
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
void CbmMvdStationPar::SetYRes(Int_t stationNumber, Double_t yres)
{
fYResolutions[stationNumber] = yres;
}  
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
void CbmMvdStationPar::SetRadLength(Int_t stationNumber, Double_t length)
{
fRadiationLength[stationNumber] = length;
}  
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
void CbmMvdStationPar::SetBeamHeight(Int_t stationNumber, Double_t beamheight)
{
if(fBeamHeights[stationNumber] != 0)
	{
	if(fBeamHeights[stationNumber] > beamheight)
		{
		fBeamHeights[stationNumber] = beamheight;
		}
	}
else
	{
	fBeamHeights[stationNumber] = beamheight;
	}
} 
// -------------------------------------------------------------------------

// -------------------------------------------------------------------------
void CbmMvdStationPar::SetBeamWidth(Int_t stationNumber, Double_t beamwidth)
{
if(fBeamWidths[stationNumber] != 0)
	{
	if(fBeamWidths[stationNumber] > beamwidth)
		{
		fBeamWidths[stationNumber] = beamwidth;
		}
	}
else
	{
	fBeamWidths[stationNumber] = beamwidth;
	}
}  
// -------------------------------------------------------------------------



ClassImp(CbmMvdStationPar)