[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 = AliEMCALGeoParams::fgkEMCALModules; //12
  Int_t nCol  = AliEMCALGeoParams::fgkEMCALCols;    //48
  Int_t nRow  = AliEMCALGeoParams::fgkEMCALRows;    //24
  Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows;    //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 = AliEMCALGeoParams::fgkEMCALModules; //12
  Int_t nCol  = AliEMCALGeoParams::fgkEMCALCols;    //48
  Int_t nRow  = AliEMCALGeoParams::fgkEMCALRows;    //24
  Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows/2;  //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));
	Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
	Int_t nCol  = AliEMCALGeoParams::fgkEMCALCols;    //48
	Int_t nRow  = AliEMCALGeoParams::fgkEMCALRows;    //24
	Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows/2;  //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++){
					  fADCpedestal[supermodule][column][row]=0.;
					  fADCchannel [supermodule][column][row]=1.;
				}
		}
	} 
}

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

  Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
  Int_t nCol  = AliEMCALGeoParams::fgkEMCALCols;    //48
  Int_t nRow  = AliEMCALGeoParams::fgkEMCALRows;    //24
  Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows/2;  //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(" %2.4f ",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(" %2.4f ",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
0c5b726e	60	Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
	61	Int_t nCol = AliEMCALGeoParams::fgkEMCALCols; //48
	62	Int_t nRow = AliEMCALGeoParams::fgkEMCALRows; //24
	63	Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows; //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
0c5b726e	88	Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
	89	Int_t nCol = AliEMCALGeoParams::fgkEMCALCols; //48
	90	Int_t nRow = AliEMCALGeoParams::fgkEMCALRows; //24
	91	Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows/2; //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
0c5b726e	118	//memset(fADCchannel ,1,124824*sizeof(Float_t));
	119	//memset(fADCpedestal,0,124824*sizeof(Float_t));
	120	Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
	121	Int_t nCol = AliEMCALGeoParams::fgkEMCALCols; //48
	122	Int_t nRow = AliEMCALGeoParams::fgkEMCALRows; //24
	123	Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows/2; //12 - Modules 11 and 12 are half modules
	124	for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
	125	if(supermodule >= 10)
	126	nRow = nRow2;
	127	for (Int_t column=0; column<nCol; column++){
	128	for (Int_t row=0; row<nRow; row++){
	129	fADCpedestal[supermodule][column][row]=0.;
	130	fADCchannel [supermodule][column][row]=1.;
	131	}
	132	}
	133	}
f565d89d	134	}
	135
	136	//________________________________________________________________
	137	void AliEMCALCalibData::Print(Option_t *option) const
	138	{
	139	// Print tables of pedestals and ADC channels widths
	140
0c5b726e	141	Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
	142	Int_t nCol = AliEMCALGeoParams::fgkEMCALCols; //48
	143	Int_t nRow = AliEMCALGeoParams::fgkEMCALRows; //24
	144	Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows/2; //12 - Modules 11 and 12 are half modules
f565d89d	145	if (strstr(option,"ped")) {
	146	printf("\n ---- Pedestal values ----\n\n");
	147	for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
9219ceb6	148	if(supermodule >= 10)
beb28b95	149	nRow = nRow2;
f565d89d	150	printf("============== Supermodule %d\n",supermodule+1);
	151	for (Int_t column=0; column<nCol; column++){
	152	for (Int_t row=0; row<nRow; row++){
636af9b4	153	printf(" %2.4f ",fADCpedestal[supermodule][column][row]);
f565d89d	154	}
	155	printf("\n");
	156	}
	157	}
	158	}
	159
	160
	161	if (strstr(option,"gain")) {
	162	printf("\n ---- ADC channel values ----\n\n");
	163	for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
9219ceb6	164	if(supermodule >= 10)
beb28b95	165	nRow = nRow2;
f565d89d	166	printf("============== Supermodule %d\n",supermodule+1);
	167	for (Int_t column=0; column<nCol; column++){
	168	for (Int_t row=0; row<nRow; row++){
636af9b4	169	printf(" %2.4f ",fADCchannel[supermodule][column][row]);
f565d89d	170	}
	171	printf("\n");
	172	}
	173	}
	174	}
	175	}
	176
	177	//________________________________________________________________
	178	Float_t AliEMCALCalibData::GetADCchannel(Int_t supermodule, Int_t column, Int_t row) const
	179	{
	180	// Set ADC channel witdth values
63b315c5	181	// All indexes start from 0!
	182	// Supermodule, column,raw should follow the ALICE convention:
	183	// supermodule 0:11, column 0:47, row 0:23
f565d89d	184
63b315c5	185	return fADCchannel[supermodule][column][row];
f565d89d	186	}
	187
	188	//________________________________________________________________
	189	Float_t AliEMCALCalibData::GetADCpedestal(Int_t supermodule, Int_t column, Int_t row) const
	190	{
	191	// Get ADC pedestal values
63b315c5	192	return fADCpedestal[supermodule][column][row];
f565d89d	193	}
	194
	195	//________________________________________________________________
	196	void AliEMCALCalibData::SetADCchannel(Int_t supermodule, Int_t column, Int_t row, Float_t value)
	197	{
	198	// Set ADC channel width values
63b315c5	199	fADCchannel[supermodule][column][row] = value;
f565d89d	200	}
	201
	202	//________________________________________________________________
	203	void AliEMCALCalibData::SetADCpedestal(Int_t supermodule, Int_t column, Int_t row, Float_t value)
	204	{
	205	// Set ADC pedestal values
63b315c5	206	fADCpedestal[supermodule][column][row] = value;
f565d89d	207	}