[u/mrichter/AliRoot.git] / EMCAL / AliEMCALCalibData.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$ */

//_________________________________________________________________________
///*-- Author: Yves Schutz (SUBATECH)
//           : Aleksei Pavlinov (WSU); Jun 30, 2006 - ALICE numbering scheme
///////////////////////////////////////////////////////////////////////////////
//                                                                           //
// class for EMCAL calibration                                               //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include <TMath.h>

#include "AliEMCALCalibData.h"

ClassImp(AliEMCALCalibData)

//________________________________________________________________
AliEMCALCalibData::AliEMCALCalibData()
{
  // Default constructor
  Reset();
}

//________________________________________________________________
AliEMCALCalibData::AliEMCALCalibData(const char* name)
{
  // Constructor
  TString namst = "Calib_";
  namst += name;
  SetName(namst.Data());
  SetTitle(namst.Data());
  Reset();
}

//________________________________________________________________
AliEMCALCalibData::AliEMCALCalibData(const AliEMCALCalibData& calibda) :
  TNamed(calibda)
{
  // copy constructor
  SetName(calibda.GetName());
  SetTitle(calibda.GetName());
  Reset();

  Int_t nSMod = 12;
  Int_t nCol  = 48;
  Int_t nRow  = 24;
  Int_t nRow2 = 12; //Modules 11 and 12 are half modules

  for(Int_t supermodule=0; supermodule<nSMod; supermodule++) {
    if(supermodule >= 10)
      nRow = nRow2;
    for(Int_t column=0; column<nCol; column++) {
      for(Int_t row=0; row<nRow; row++) {
	fADCchannel[supermodule][column][row] = 
	  calibda.GetADCchannel(supermodule,column,row);
	fADCpedestal[supermodule][column][row] = 
	  calibda.GetADCpedestal(supermodule,column,row);
      }
    }
  }
}


//________________________________________________________________
AliEMCALCalibData &AliEMCALCalibData::operator =(const AliEMCALCalibData& calibda)
{
  // assignment operator
  SetName(calibda.GetName());
  SetTitle(calibda.GetName());
  Reset();

  Int_t nSMod = 12;
  Int_t nCol  = 48;
  Int_t nRow  = 24;
  Int_t nRow2 = 12; //Modules 11 and 12 are half modules

  for(Int_t supermodule=0; supermodule<nSMod; supermodule++) {
    if(supermodule >= 10)
      nRow = nRow2;
    for(Int_t column=0; column<nCol; column++) {
      for(Int_t row=0; row<nRow; row++) {
	fADCchannel[supermodule][column][row] = 
	  calibda.GetADCchannel(supermodule,column,row);
	fADCpedestal[supermodule][column][row] = 
	  calibda.GetADCpedestal(supermodule,column,row);
      }
    }
  }
  return *this;
}

//________________________________________________________________
AliEMCALCalibData::~AliEMCALCalibData()
{
  // Destructor
}

//________________________________________________________________
void AliEMCALCalibData::Reset()
{
  // Set all pedestals to 0 and all ADC channels widths to 1
  memset(fADCchannel ,1,12*48*24*sizeof(Float_t));
  memset(fADCpedestal,0,12*48*24*sizeof(Float_t));
}

//________________________________________________________________
void  AliEMCALCalibData::Print(Option_t *option) const
{
  // Print tables of pedestals and ADC channels widths

  Int_t nSMod = 12;
  Int_t nCol  = 48;
  Int_t nRow  = 24;
  Int_t nRow2 = 12; //Modules 11 and 12 are half modules
  if (strstr(option,"ped")) {
    printf("\n	----	Pedestal values	----\n\n");
    for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
      if(supermodule >= 10)
	nRow = nRow2;
      printf("============== Supermodule %d\n",supermodule+1);
      for (Int_t column=0; column<nCol; column++){
	for (Int_t row=0; row<nRow; row++){
	  printf("%4.1f",fADCpedestal[supermodule][column][row]);
	}
	printf("\n");
      }
    } 
  }


  if (strstr(option,"gain")) {
    printf("\n	----	ADC channel values	----\n\n");
    for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
      if(supermodule >= 10) 
	nRow = nRow2;
      printf("============== Supermodule %d\n",supermodule+1);
      for (Int_t column=0; column<nCol; column++){
	for (Int_t row=0; row<nRow; row++){
	  printf("%4.1f",fADCchannel[supermodule][column][row]);
	}
	printf("\n");
      }
    }   
  }
}

//________________________________________________________________
Float_t AliEMCALCalibData::GetADCchannel(Int_t supermodule, Int_t column, Int_t row) const
{
  // Set ADC channel witdth values
  // All indexes start from 0!
  // Supermodule, column,raw should follow the ALICE convention:
  // supermodule 0:11, column 0:47, row 0:23

  return fADCchannel[supermodule][column][row];
}

//________________________________________________________________
Float_t AliEMCALCalibData::GetADCpedestal(Int_t supermodule, Int_t column, Int_t row) const
{
  // Get ADC pedestal values
 return fADCpedestal[supermodule][column][row];
}

//________________________________________________________________
void AliEMCALCalibData::SetADCchannel(Int_t supermodule, Int_t column, Int_t row, Float_t value)
{ 
  // Set ADC channel width values
  fADCchannel[supermodule][column][row] = value;
}

//________________________________________________________________
void AliEMCALCalibData::SetADCpedestal(Int_t supermodule, Int_t column, Int_t row, Float_t value)
{
  // Set ADC pedestal values
  fADCpedestal[supermodule][column][row] = value;
}
Commit	Line	Data
f565d89d	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
63b315c5	18	//_________________________________________________________________________
	19	///*-- Author: Yves Schutz (SUBATECH)
	20	// : Aleksei Pavlinov (WSU); Jun 30, 2006 - ALICE numbering scheme
f565d89d	21	///////////////////////////////////////////////////////////////////////////////
	22	// //
	23	// class for EMCAL calibration //
	24	// //
	25	///////////////////////////////////////////////////////////////////////////////
	26
090026bf	27	#include <TMath.h>
090026bf	28
f565d89d	29	#include "AliEMCALCalibData.h"
	30
	31	ClassImp(AliEMCALCalibData)
	32
	33	//________________________________________________________________
	34	AliEMCALCalibData::AliEMCALCalibData()
	35	{
	36	// Default constructor
	37	Reset();
	38	}
	39
	40	//________________________________________________________________
	41	AliEMCALCalibData::AliEMCALCalibData(const char* name)
	42	{
	43	// Constructor
	44	TString namst = "Calib_";
	45	namst += name;
	46	SetName(namst.Data());
	47	SetTitle(namst.Data());
	48	Reset();
	49	}
	50
	51	//________________________________________________________________
	52	AliEMCALCalibData::AliEMCALCalibData(const AliEMCALCalibData& calibda) :
	53	TNamed(calibda)
	54	{
	55	// copy constructor
	56	SetName(calibda.GetName());
	57	SetTitle(calibda.GetName());
	58	Reset();
	59
	60	Int_t nSMod = 12;
	61	Int_t nCol = 48;
	62	Int_t nRow = 24;
beb28b95	63	Int_t nRow2 = 12; //Modules 11 and 12 are half modules
f565d89d	64
f565d89d	65	for(Int_t supermodule=0; supermodule<nSMod; supermodule++) {
9219ceb6	66	if(supermodule >= 10)
beb28b95	67	nRow = nRow2;
f565d89d	68	for(Int_t column=0; column<nCol; column++) {
	69	for(Int_t row=0; row<nRow; row++) {
	70	fADCchannel[supermodule][column][row] =
	71	calibda.GetADCchannel(supermodule,column,row);
	72	fADCpedestal[supermodule][column][row] =
	73	calibda.GetADCpedestal(supermodule,column,row);
	74	}
	75	}
	76	}
	77	}
	78
	79
	80	//________________________________________________________________
	81	AliEMCALCalibData &AliEMCALCalibData::operator =(const AliEMCALCalibData& calibda)
	82	{
	83	// assignment operator
	84	SetName(calibda.GetName());
	85	SetTitle(calibda.GetName());
	86	Reset();
	87
	88	Int_t nSMod = 12;
	89	Int_t nCol = 48;
	90	Int_t nRow = 24;
beb28b95	91	Int_t nRow2 = 12; //Modules 11 and 12 are half modules
f565d89d	92
f565d89d	93	for(Int_t supermodule=0; supermodule<nSMod; supermodule++) {
9219ceb6	94	if(supermodule >= 10)
beb28b95	95	nRow = nRow2;
f565d89d	96	for(Int_t column=0; column<nCol; column++) {
	97	for(Int_t row=0; row<nRow; row++) {
	98	fADCchannel[supermodule][column][row] =
	99	calibda.GetADCchannel(supermodule,column,row);
	100	fADCpedestal[supermodule][column][row] =
	101	calibda.GetADCpedestal(supermodule,column,row);
	102	}
	103	}
	104	}
	105	return *this;
	106	}
	107
	108	//________________________________________________________________
	109	AliEMCALCalibData::~AliEMCALCalibData()
	110	{
	111	// Destructor
	112	}
	113
	114	//________________________________________________________________
	115	void AliEMCALCalibData::Reset()
	116	{
	117	// Set all pedestals to 0 and all ADC channels widths to 1
	118	memset(fADCchannel ,1,124824*sizeof(Float_t));
	119	memset(fADCpedestal,0,124824*sizeof(Float_t));
	120	}
	121
	122	//________________________________________________________________
	123	void AliEMCALCalibData::Print(Option_t *option) const
	124	{
	125	// Print tables of pedestals and ADC channels widths
	126
	127	Int_t nSMod = 12;
	128	Int_t nCol = 48;
	129	Int_t nRow = 24;
beb28b95	130	Int_t nRow2 = 12; //Modules 11 and 12 are half modules
f565d89d	131	if (strstr(option,"ped")) {
	132	printf("\n ---- Pedestal values ----\n\n");
	133	for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
9219ceb6	134	if(supermodule >= 10)
beb28b95	135	nRow = nRow2;
f565d89d	136	printf("============== Supermodule %d\n",supermodule+1);
	137	for (Int_t column=0; column<nCol; column++){
	138	for (Int_t row=0; row<nRow; row++){
	139	printf("%4.1f",fADCpedestal[supermodule][column][row]);
	140	}
	141	printf("\n");
	142	}
	143	}
	144	}
	145
	146
	147	if (strstr(option,"gain")) {
	148	printf("\n ---- ADC channel values ----\n\n");
	149	for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
9219ceb6	150	if(supermodule >= 10)
beb28b95	151	nRow = nRow2;
f565d89d	152	printf("============== Supermodule %d\n",supermodule+1);
	153	for (Int_t column=0; column<nCol; column++){
	154	for (Int_t row=0; row<nRow; row++){
	155	printf("%4.1f",fADCchannel[supermodule][column][row]);
	156	}
	157	printf("\n");
	158	}
	159	}
	160	}
	161	}
	162
	163	//________________________________________________________________
	164	Float_t AliEMCALCalibData::GetADCchannel(Int_t supermodule, Int_t column, Int_t row) const
	165	{
	166	// Set ADC channel witdth values
63b315c5	167	// All indexes start from 0!
	168	// Supermodule, column,raw should follow the ALICE convention:
	169	// supermodule 0:11, column 0:47, row 0:23
f565d89d	170
63b315c5	171	return fADCchannel[supermodule][column][row];
f565d89d	172	}
	173
	174	//________________________________________________________________
	175	Float_t AliEMCALCalibData::GetADCpedestal(Int_t supermodule, Int_t column, Int_t row) const
	176	{
	177	// Get ADC pedestal values
63b315c5	178	return fADCpedestal[supermodule][column][row];
f565d89d	179	}
	180
	181	//________________________________________________________________
	182	void AliEMCALCalibData::SetADCchannel(Int_t supermodule, Int_t column, Int_t row, Float_t value)
	183	{
	184	// Set ADC channel width values
63b315c5	185	fADCchannel[supermodule][column][row] = value;
f565d89d	186	}
	187
	188	//________________________________________________________________
	189	void AliEMCALCalibData::SetADCpedestal(Int_t supermodule, Int_t column, Int_t row, Float_t value)
	190	{
	191	// Set ADC pedestal values
63b315c5	192	fADCpedestal[supermodule][column][row] = value;
f565d89d	193	}