/////////////////////////////////////////////////////////////
//
//  CbmEmcReader
//
//  Reader of Emc geometry file (ASCII)
//
//  Created 14/08/06  by S.Spataro
//
/////////////////////////////////////////////////////////////// 

#include <iostream>
#include <fstream>
#include <string>
#include "TString.h"


#include "CbmEmcReader.h"

using namespace std;

CbmEmcReader::CbmEmcReader(TString name)
{ 
  sName = name;
  CbmEmcReader::fill_vector();
}

CbmEmcReader::~CbmEmcReader()
{
}

void CbmEmcReader::fill_vector()
{
    fstream f;
    string str;
    DataG4 data;
    int line_nr;
   
    f.open(sName.Data(), ios::in);
    cout << "CbmEmcReader: EMC geometry file == " << sName.Data() << endl;
    if(!f) {
	cout<<"CbmEmcReader Error: Data file cannot be open!"<<endl;
	exit(1);
    } else {
	line_nr = 1;
	while(line_nr<=3) {getline(f,str,'\n'); line_nr++;};
	while(f) {
	 
	    f>>data.module>>data.row>>data.crystal>>data.theta>>data.phi>>data.tau>>data.posX>>data.posY>>data.posZ>>data.pDz
	     >>data.pTheta>>data.pPhi>>data.pDy1>>data.pDx1>>data.pDx2>>data.pAlp1>>data.pDy2>>data.pDx3>>data.pDx4>>data.pAlp2;
	    CbmEmcReader::g4data.push_back(data); 
	};
	f.close();
    };
    f.clear();

    return;
}

int CbmEmcReader::GetMaxModules()
{
    int modules = 0;
    
    for(unsigned int i=0; i<CbmEmcReader::g4data.size(); i++) 
    if(modules<CbmEmcReader::g4data[i].module) 
    modules = CbmEmcReader::g4data[i].module;
    
    return modules;
}

int CbmEmcReader::GetMaxRows(int module)
{
    int rows = 0;
    
    for(unsigned int i=0; i<CbmEmcReader::g4data.size(); i++) 
    if(module==CbmEmcReader::g4data[i].module) 
    if(rows<CbmEmcReader::g4data[i].row) 
    rows = CbmEmcReader::g4data[i].row;
    
    return rows;
}

int CbmEmcReader::GetMaxCrystals(int module, int row)
{
    int crystals = 0;
    
    for(unsigned int i=0; i<CbmEmcReader::g4data.size(); i++) 
    if((module==CbmEmcReader::g4data[i].module)&&(row==CbmEmcReader::g4data[i].row)) 
    if(crystals<CbmEmcReader::g4data[i].crystal) 
    crystals = CbmEmcReader::g4data[i].crystal;
    
    return crystals;
}

int CbmEmcReader::GetMinModules()
{
    int modules = 1000000;
    
    for(unsigned int i=0; i<CbmEmcReader::g4data.size(); i++) 
    if(modules>CbmEmcReader::g4data[i].module) 
    modules = CbmEmcReader::g4data[i].module;
    
    return modules;
}

int CbmEmcReader::GetMinRows(int module)
{
    int rows = 1000000;
    
    for(unsigned int i=0; i<CbmEmcReader::g4data.size(); i++) 
    if(module==CbmEmcReader::g4data[i].module) 
    if(rows>CbmEmcReader::g4data[i].row) 
    rows = CbmEmcReader::g4data[i].row;
    
    return rows;
}

int CbmEmcReader::GetMinCrystals(int module, int row)
{
    int crystals = 1000000;
    
    for(unsigned int i=0; i<CbmEmcReader::g4data.size(); i++) 
    if((module==CbmEmcReader::g4data[i].module)&&(row==CbmEmcReader::g4data[i].row)) 
    if(crystals>CbmEmcReader::g4data[i].crystal) 
    crystals = CbmEmcReader::g4data[i].crystal;
    
    return crystals;
}

DataG4 CbmEmcReader::GetData(int module, int row, int crystal)
{
    DataG4 data;
    data.module = -1; // Module not present

    for(unsigned int i=0; i<CbmEmcReader::g4data.size(); i++) {
	if((module==CbmEmcReader::g4data[i].module)&&(row==CbmEmcReader::g4data[i].row)&&(crystal==CbmEmcReader::g4data[i].crystal)) {
	    data.module = CbmEmcReader::g4data[i].module;
	    data.row = CbmEmcReader::g4data[i].row;
	    data.crystal = CbmEmcReader::g4data[i].crystal;
	    data.theta = CbmEmcReader::g4data[i].theta;
	    data.phi = CbmEmcReader::g4data[i].phi;
	    data.tau = CbmEmcReader::g4data[i].tau;
	    data.posX = CbmEmcReader::g4data[i].posX;
	    data.posY = CbmEmcReader::g4data[i].posY;
	    data.posZ = CbmEmcReader::g4data[i].posZ;
	    data.pDz = CbmEmcReader::g4data[i].pDz;
	    data.pTheta = CbmEmcReader::g4data[i].pTheta;
	    data.pPhi = CbmEmcReader::g4data[i].pPhi;
	    data.pDy1 = CbmEmcReader::g4data[i].pDy1;
	    data.pDx1 = CbmEmcReader::g4data[i].pDx1;
	    data.pDx2 = CbmEmcReader::g4data[i].pDx2;
	    data.pAlp1 = CbmEmcReader::g4data[i].pAlp1;
	    data.pDy2 = CbmEmcReader::g4data[i].pDy2;
	    data.pDx3 = CbmEmcReader::g4data[i].pDx3;
	    data.pDx4 = CbmEmcReader::g4data[i].pDx4;
	    data.pAlp2 = CbmEmcReader::g4data[i].pAlp2;
	};
    };
    
    return data;
}