[u/mrichter/AliRoot.git] / EMCAL / SMcalib / AliEMCALCalibAPD.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id: $ */

// Objects of this class read txt file with APD number data
//

#include <fstream>
#include <TString.h>

#include "AliEMCALCalibAPD.h"

const int kMaxLen = 1000; // maximum length of single line (# of characters)
// an OK line with complete info should have a certain number of characters
const int kMinLenAPDLine = 135;
const int kMaxLenAPDLine = 1000;

ClassImp(AliEMCALCalibAPD)

//____________________________________________________________________________
AliEMCALCalibAPD::AliEMCALCalibAPD() : 
  fNCalibAPD(0),
  fData(0)
{
  //Default constructor.
}

//____________________________________________________________________________
void AliEMCALCalibAPD::ReadCalibAPDInfo(Int_t nAPD, const TString &txtFileName)
{
  //Read data from txt file. ; coordinates given on SuperModule basis

  std::ifstream inputFile(txtFileName.Data());
  if (!inputFile) {
    printf("AliEMCALCalibAPD::ReadCalibAPDInfo - Cannot open the APD info file %s\n", txtFileName.Data());
    return;
  }

  fNCalibAPD = nAPD;
  if (fData) delete [] fData;
  fData = new AliEMCALCalibAPDData[fNCalibAPD];

  char line[kMaxLen];

  /* DS: header lines skipped when switched to white-spaced separeted dat files instead of csv:
     conversion from spreadsheet can be done a la
     sed 's/,/ /g' APD-database-Houston.csv | egrep ^2 | awk '{if (NF==19) {print $0}}' > APD-database-Houston.dat
     - meaning "replace , with whitespace, only get the columns that start with the number 2 (Houston APDs),
     and only get rows with the full 19 fields

  // get header lines:
  inputFile.getline(line, kMaxLen);
  //  printf(" 1st header line character count %d\n", inputFile.gcount());
  inputFile.getline(line, kMaxLen);
  //  printf(" 2nd header line character count %d\n", inputFile.gcount());
  */

  // variables for reading
  int i1,i2,i3,i4,i5;
  float f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12;
  char c1[10], c2[20];

  int j = 0; // index for the ones that were read OK
  for (Int_t i = 0; i < fNCalibAPD; i++) {
    AliEMCALCalibAPDData &t = fData[j];
    if (!inputFile) {
      printf("AliEMCALCalibAPD::ReadCalibAPDInfo - Error while reading input file; likely EOF..\n");
      fNCalibAPD = j; // that's how many we actually read succesfully
      printf("AliEMCALCalibAPD::ReadCalibAPDInfo - read %d OK\n", fNCalibAPD);
      return;
    }

    // use some temporary/local values to perhaps help with the order when 
    // trying to read all the many fields in a line..
    inputFile.getline(line, kMaxLen);
    int nchar = inputFile.gcount();
    //printf(" line %d ok %d - character count %d\n", i, j, nchar);

    if (nchar>kMinLenAPDLine && nchar<kMaxLenAPDLine) {
      // looks like the line has about the right number of characters, let's
      // try to decode it now..

      //      printf("input: %s\n",line);
      sscanf(line, "%d %u %s %s %d %d %f %f %f %f %f %f %f %f %f %d %f %f %f",
	     &i1, &i2, c1, c2, &i3, &i4, // header-type info
	     &f1, &f2, &f3, &f4, &f5, &f6, &f7, &f8, &f9, // measurements
	     &i5, &f10, &f11, &f12); // Hamamatsu

      //      printf("after scanf: %s\n",line);
      /*
      printf("%d,%u,%s,%s,%d,%d,%g,%g,%g,%g,%g,%g,%g,%g,%g,%d,%g,%13.11f,%g\n",
	     i1, i2, c1, c2, i3, i4, // header-type info
	     f1, f2, f3, f4, f5, f6, f7, f8, f9, // measurements
	     i5, f10, f11, f12); // Hamamatsu
      */
      // assign the variables:
      t.fAPDNum = i1; 
      t.fSerialNum = i2; 
      sprintf(t.fStatus,"%s",c1);
      sprintf(t.fLocation,"%s",c2);
      t.fRunNum = i3; 
      t.fTestPos = i4; 

      t.fV30 = f1; 
      t.fV50 = f2; 
      t.fVoltCoeff = f3;
      t.fPar[0] = f4;
      t.fPar[1] = f5;
      t.fPar[2] = f6;
      t.fParErr[0] = f7;
      t.fParErr[1] = f8;
      t.fParErr[2] = f9;

      t.fBreakDown = i5; 
      t.fHamV50 = f10;
      t.fDarkCurrent = f11;
      t.fTestTemp = f12;

      j++; // increment our 'OK' counter..
    } // line length appears OK      
  } // i, APD

  inputFile.close();

  return;
}

//____________________________________________________________________________
void AliEMCALCalibAPD::WriteCalibAPDInfo(const TString &txtFileName)
{
  // write data to txt file. ; coordinates given on SuperModule basis

  std::ofstream outputFile(txtFileName.Data());
  if (!outputFile) {
    printf("AliEMCALCalibAPD::WriteCalibAPDInfo - Cannot open the APD output file %s\n", txtFileName.Data());
    return;
  }

  char *comma = ",";

  for (Int_t i = 0; i < fNCalibAPD; i++) {
    AliEMCALCalibAPDData &t = fData[i];

    outputFile << t.fAPDNum << comma
	       << t.fSerialNum << comma // Serial Number; from Hamamatsu
	       << t.fStatus << comma //
	       << t.fLocation << comma //
	       << t.fRunNum << comma
	       << t.fTestPos << comma
	       << t.fV30 << comma // Catania/Houston Voltage V30 (V) at T = 25 deg C
	       << t.fV50 << comma 
	       << t.fVoltCoeff << comma // 1/M x dM/dV
	       << t.fPar[0] << comma // fit parameters, p0,p1,p2 - for ADC vs bias measurement
	       << t.fPar[1] << comma
	       << t.fPar[2] << comma
	       << t.fParErr[0] << comma // error on fit parameters	
	       << t.fParErr[1] << comma 
	       << t.fParErr[2] << comma 
	       << t.fBreakDown << comma // Hamamatsu Breakdown Voltage (V)	
	       << t.fHamV50 << comma; // Hamamatsu Voltage V50 (V)
    // I wasn't able to quite reproduce the values as they appeared on the 
    // original file: e.g. dark current is not always 11-field - if last digit is zero..
    // the other floats have 6 significant digits, but sometimes switch to 
    // scientific notation - don't know how to handle this for varying length of fields.. -leave it as is for now..
    outputFile << t.fDarkCurrent << comma; // Hamamatsu Dark Current (A)	
    /*
    // some tweaking for the dark-current field to get the output the same
    // as on the original CSV
    outputFile.precision(11);	
    outputFile << fixed << t.fDarkCurrent << comma; // Hamamatsu Dark Current (A)	
    // back to normal..
    outputFile.precision(6); outputFile.unsetf(ios_base::floatfield);	
    */

    outputFile << t.fTestTemp // Hamamatsu Testing Temperature (deg C)	
	       << endl; 

  } // i, APD

  outputFile.close();

  return;
}

//____________________________________________________________________________
AliEMCALCalibAPD::~AliEMCALCalibAPD()
{
  delete [] fData;
}

//____________________________________________________________________________
AliEMCALCalibAPD::AliEMCALCalibAPDData AliEMCALCalibAPD::GetCalibAPDDataId(Int_t apdIndex)const
{
  AliEMCALCalibAPDData t;  // just to maybe prevent a crash, but we are returning something not-initialized so maybe not better really..
  if (!fData)
    return t;

  return fData[apdIndex];
}

//____________________________________________________________________________
AliEMCALCalibAPD::AliEMCALCalibAPDData AliEMCALCalibAPD::GetCalibAPDDataNum(Int_t apdNum)const
{
  AliEMCALCalibAPDData t;  // just to maybe prevent a crash, but we are returning something not-initialized so maybe not better really..
  if (!fData)
    return t;

  for (int i=0; i<fNCalibAPD; i++) {
    if (fData[i].fAPDNum == apdNum) {
      return fData[i];
    }
  }

  // if we made it to here, then no match was found - just return then
  return t;
}
Commit	Line	Data
	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/* $Id: $ */
	17
	18	// Objects of this class read txt file with APD number data
	19	//
	20
	21	#include <fstream>
	22	#include <TString.h>
	23
	24	#include "AliEMCALCalibAPD.h"
	25
	26	const int kMaxLen = 1000; // maximum length of single line (# of characters)
	27	// an OK line with complete info should have a certain number of characters
	28	const int kMinLenAPDLine = 135;
	29	const int kMaxLenAPDLine = 1000;
	30
	31	ClassImp(AliEMCALCalibAPD)
	32
	33	//____________________________________________________________________________
	34	AliEMCALCalibAPD::AliEMCALCalibAPD() :
	35	fNCalibAPD(0),
	36	fData(0)
	37	{
	38	//Default constructor.
	39	}
	40
	41	//____________________________________________________________________________
	42	void AliEMCALCalibAPD::ReadCalibAPDInfo(Int_t nAPD, const TString &txtFileName)
	43	{
	44	//Read data from txt file. ; coordinates given on SuperModule basis
	45
	46	std::ifstream inputFile(txtFileName.Data());
	47	if (!inputFile) {
	48	printf("AliEMCALCalibAPD::ReadCalibAPDInfo - Cannot open the APD info file %s\n", txtFileName.Data());
	49	return;
	50	}
	51
	52	fNCalibAPD = nAPD;
	53	if (fData) delete [] fData;
	54	fData = new AliEMCALCalibAPDData[fNCalibAPD];
	55
	56	char line[kMaxLen];
	57
	58	/* DS: header lines skipped when switched to white-spaced separeted dat files instead of csv:
	59	conversion from spreadsheet can be done a la
	60	sed 's/,/ /g' APD-database-Houston.csv \| egrep ^2 \| awk '{if (NF==19) {print $0}}' > APD-database-Houston.dat
	61	- meaning "replace , with whitespace, only get the columns that start with the number 2 (Houston APDs),
	62	and only get rows with the full 19 fields
	63
	64	// get header lines:
	65	inputFile.getline(line, kMaxLen);
	66	// printf(" 1st header line character count %d\n", inputFile.gcount());
	67	inputFile.getline(line, kMaxLen);
	68	// printf(" 2nd header line character count %d\n", inputFile.gcount());
	69	*/
	70
	71	// variables for reading
	72	int i1,i2,i3,i4,i5;
	73	float f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12;
	74	char c1[10], c2[20];
	75
	76	int j = 0; // index for the ones that were read OK
	77	for (Int_t i = 0; i < fNCalibAPD; i++) {
	78	AliEMCALCalibAPDData &t = fData[j];
	79	if (!inputFile) {
	80	printf("AliEMCALCalibAPD::ReadCalibAPDInfo - Error while reading input file; likely EOF..\n");
	81	fNCalibAPD = j; // that's how many we actually read succesfully
	82	printf("AliEMCALCalibAPD::ReadCalibAPDInfo - read %d OK\n", fNCalibAPD);
	83	return;
	84	}
	85
	86	// use some temporary/local values to perhaps help with the order when
	87	// trying to read all the many fields in a line..
	88	inputFile.getline(line, kMaxLen);
	89	int nchar = inputFile.gcount();
	90	//printf(" line %d ok %d - character count %d\n", i, j, nchar);
	91
	92	if (nchar>kMinLenAPDLine && nchar<kMaxLenAPDLine) {
	93	// looks like the line has about the right number of characters, let's
	94	// try to decode it now..
	95
	96	// printf("input: %s\n",line);
	97	sscanf(line, "%d %u %s %s %d %d %f %f %f %f %f %f %f %f %f %d %f %f %f",
	98	&i1, &i2, c1, c2, &i3, &i4, // header-type info
	99	&f1, &f2, &f3, &f4, &f5, &f6, &f7, &f8, &f9, // measurements
	100	&i5, &f10, &f11, &f12); // Hamamatsu
	101
	102	// printf("after scanf: %s\n",line);
	103	/*
	104	printf("%d,%u,%s,%s,%d,%d,%g,%g,%g,%g,%g,%g,%g,%g,%g,%d,%g,%13.11f,%g\n",
	105	i1, i2, c1, c2, i3, i4, // header-type info
	106	f1, f2, f3, f4, f5, f6, f7, f8, f9, // measurements
	107	i5, f10, f11, f12); // Hamamatsu
	108	*/
	109	// assign the variables:
	110	t.fAPDNum = i1;
	111	t.fSerialNum = i2;
	112	sprintf(t.fStatus,"%s",c1);
	113	sprintf(t.fLocation,"%s",c2);
	114	t.fRunNum = i3;
	115	t.fTestPos = i4;
	116
	117	t.fV30 = f1;
	118	t.fV50 = f2;
	119	t.fVoltCoeff = f3;
	120	t.fPar[0] = f4;
	121	t.fPar[1] = f5;
	122	t.fPar[2] = f6;
	123	t.fParErr[0] = f7;
	124	t.fParErr[1] = f8;
	125	t.fParErr[2] = f9;
	126
	127	t.fBreakDown = i5;
	128	t.fHamV50 = f10;
	129	t.fDarkCurrent = f11;
	130	t.fTestTemp = f12;
	131
	132	j++; // increment our 'OK' counter..
	133	} // line length appears OK
	134	} // i, APD
	135
	136	inputFile.close();
	137
	138	return;
	139	}
	140
	141	//____________________________________________________________________________
	142	void AliEMCALCalibAPD::WriteCalibAPDInfo(const TString &txtFileName)
	143	{
	144	// write data to txt file. ; coordinates given on SuperModule basis
	145
	146	std::ofstream outputFile(txtFileName.Data());
	147	if (!outputFile) {
	148	printf("AliEMCALCalibAPD::WriteCalibAPDInfo - Cannot open the APD output file %s\n", txtFileName.Data());
	149	return;
	150	}
	151
	152	char *comma = ",";
	153
	154	for (Int_t i = 0; i < fNCalibAPD; i++) {
	155	AliEMCALCalibAPDData &t = fData[i];
	156
	157	outputFile << t.fAPDNum << comma
	158	<< t.fSerialNum << comma // Serial Number; from Hamamatsu
	159	<< t.fStatus << comma //
	160	<< t.fLocation << comma //
	161	<< t.fRunNum << comma
	162	<< t.fTestPos << comma
	163	<< t.fV30 << comma // Catania/Houston Voltage V30 (V) at T = 25 deg C
	164	<< t.fV50 << comma
	165	<< t.fVoltCoeff << comma // 1/M x dM/dV
	166	<< t.fPar[0] << comma // fit parameters, p0,p1,p2 - for ADC vs bias measurement
	167	<< t.fPar[1] << comma
	168	<< t.fPar[2] << comma
	169	<< t.fParErr[0] << comma // error on fit parameters
	170	<< t.fParErr[1] << comma
	171	<< t.fParErr[2] << comma
	172	<< t.fBreakDown << comma // Hamamatsu Breakdown Voltage (V)
	173	<< t.fHamV50 << comma; // Hamamatsu Voltage V50 (V)
	174	// I wasn't able to quite reproduce the values as they appeared on the
	175	// original file: e.g. dark current is not always 11-field - if last digit is zero..
	176	// the other floats have 6 significant digits, but sometimes switch to
	177	// scientific notation - don't know how to handle this for varying length of fields.. -leave it as is for now..
	178	outputFile << t.fDarkCurrent << comma; // Hamamatsu Dark Current (A)
	179	/*
	180	// some tweaking for the dark-current field to get the output the same
	181	// as on the original CSV
	182	outputFile.precision(11);
	183	outputFile << fixed << t.fDarkCurrent << comma; // Hamamatsu Dark Current (A)
	184	// back to normal..
	185	outputFile.precision(6); outputFile.unsetf(ios_base::floatfield);
	186	*/
	187
	188	outputFile << t.fTestTemp // Hamamatsu Testing Temperature (deg C)
	189	<< endl;
	190
	191	} // i, APD
	192
	193	outputFile.close();
	194
	195	return;
	196	}
	197
	198	//____________________________________________________________________________
	199	AliEMCALCalibAPD::~AliEMCALCalibAPD()
	200	{
	201	delete [] fData;
	202	}
	203
	204	//____________________________________________________________________________
	205	AliEMCALCalibAPD::AliEMCALCalibAPDData AliEMCALCalibAPD::GetCalibAPDDataId(Int_t apdIndex)const
	206	{
	207	AliEMCALCalibAPDData t; // just to maybe prevent a crash, but we are returning something not-initialized so maybe not better really..
	208	if (!fData)
	209	return t;
	210
	211	return fData[apdIndex];
	212	}
	213
	214	//____________________________________________________________________________
	215	AliEMCALCalibAPD::AliEMCALCalibAPDData AliEMCALCalibAPD::GetCalibAPDDataNum(Int_t apdNum)const
	216	{
	217	AliEMCALCalibAPDData t; // just to maybe prevent a crash, but we are returning something not-initialized so maybe not better really..
	218	if (!fData)
	219	return t;
	220
	221	for (int i=0; i<fNCalibAPD; i++) {
	222	if (fData[i].fAPDNum == apdNum) {
	223	return fData[i];
	224	}
	225	}
	226
	227	// if we made it to here, then no match was found - just return then
	228	return t;
	229	}
	230
	231