[u/mrichter/AliRoot.git] / HLT / PHOS / AliHLTPHOSEmcCalibData.cxx

// $Id$

/**************************************************************************
 * 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.                  * 
 *                                                                        *
 * 7/1- 2007                                                              *
 * Author: Per Thomas Hille                                               *
 * code modified to include the two differnt gains of each detection      *
 * channel and also to include the Pedestals loaded into the ALTRO        * 
 * registers for basline subtraction and zero supression                  * 
 *                                                                        *
 *                                                                        * 
 **************************************************************************/

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
// class for PHOS EmCal calibration                                          //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include "AliHLTPHOSEmcCalibData.h"

ClassImp(AliHLTPHOSEmcCalibData)

//________________________________________________________________
  AliHLTPHOSEmcCalibData::AliHLTPHOSEmcCalibData(): TNamed()
{
  // Default constructor
  Reset();
}

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

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


  //  for(Int_t gain = 0; gain < N_GAINS; gain ++){
  for(Int_t module=0; module<NMODULES; module++) {
    for(Int_t column=0; column< NXCOLUMNSMOD; column++) {
      for(Int_t row=0; row<NZROWSMOD; row++) {
	for(Int_t gain = 0; gain < NGAINS; gain ++){
	  fADCchannelEnergy[module][column][row][gain] =  calibda.fADCchannelEnergy[module][column][row][gain];
	  fADCpedestalEmcMeasured[module][column][row][gain] = calibda.fADCpedestalEmcMeasured[module][column][row][gain];
	}
      }  
    }
  }
}

//________________________________________________________________
AliHLTPHOSEmcCalibData &AliHLTPHOSEmcCalibData::operator =(const AliHLTPHOSEmcCalibData& calibda)
{
  // assignment operator

  if(this != &calibda) { 

    SetName(calibda.GetName());
    SetTitle(calibda.GetName());
    //    for(Int_t gain = 0; gain < N_GAINS; gain ++){
    for(Int_t module=0; module<NMODULES; module++) {
      for(Int_t column=0; column< NXCOLUMNSMOD; column++) {
	for(Int_t row=0; row<NZROWSMOD; row++) {
	  for(Int_t gain = 0; gain < NGAINS; gain ++){
	    fADCchannelEnergy[module][column][row][gain] = calibda.fADCchannelEnergy[module][column][row][gain];
	    fADCpedestalEmcMeasured[module][column][row][gain] = calibda.fADCpedestalEmcMeasured[module][column][row][gain];
	  }
	}
      }
    }
  }
  return *this;
}

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

//________________________________________________________________
void AliHLTPHOSEmcCalibData::Reset()
{
  // Set all pedestals and all ADC channels to its ideal values = 1.

  //  for(Int_t gain = 0; gain < PHOS_GAINS; gain ++){
    for (Int_t module=0; module<NMODULES; module++){
      for (Int_t column=0; column< NXCOLUMNSMOD; column++){
	for (Int_t row=0; row<NZROWSMOD; row++){
	  for(Int_t gain = 0; gain < NGAINS; gain ++){
	  fADCpedestalEmcMeasured[module][column][row][gain] = 0.;
	  fADCchannelEnergy[module][column][row][gain]  = 1.;
	}
      }
    }
  }
}

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

  if (strstr(option,"ped")) {
    printf("\n	----	EMC Pedestal values	----\n\n");
    //    for(Int_t gain = 0; gain < N_GAINS; gain ++){
      for (Int_t module=0; module<NMODULES; module++){
	printf("============== Module %d\n",module+1);
	for (Int_t column=0; column< NXCOLUMNSMOD; column++){
	  for (Int_t row=0; row<NZROWSMOD; row++){
	    for(Int_t gain = 0; gain < NGAINS; gain ++){
	    printf("%4.1f",fADCpedestalEmcMeasured[module][column][row][gain]);
	  }
	  printf("\n");
	}
      }
    }
  }

  if (strstr(option,"gain")) {
    printf("\n	----	EMC ADC channel values	----\n\n");
    //   for(Int_t gain = 0; gain < NGAINS; gain ++){ 
    for (Int_t module=0; module<NMODULES; module++){
      printf("============== Module %d\n",module+1);
      for (Int_t column=0; column< NXCOLUMNSMOD; column++){
	for (Int_t row=0; row<NZROWSMOD; row++){
	  for(Int_t gain = 0; gain < NGAINS; gain ++){ 
	    printf("%4.1f",fADCchannelEnergy[module][column][row][gain]);
	  }
	  printf("\n");
	}
      }
    }  
  }
}

//________________________________________________________________
void  
AliHLTPHOSEmcCalibData::MakeADCpedestalCorrectionTable()
{
  for (Int_t module=0; module<NMODULES; module++){
    printf("============== Module %d\n",module+1);
    for (Int_t column=0; column< NXCOLUMNSMOD; column++){
      for (Int_t row=0; row<NZROWSMOD; row++){
	for(Int_t gain = 0; gain < NGAINS; gain ++){ 
	  fADCpedestalCorrectionTable[module][column][row][gain] = fADCpedestalEmcMeasured[module][column][row][gain] - fADCpedestalAltroReg[module][column][row][gain];
	  //	  printf("%4.1f",fADCchannelEnergy[module][column][row][gain]);
	  //
	}
	printf("\n");
      }
    }
  }    
}


//________________________________________________________________
Float_t AliHLTPHOSEmcCalibData::GetADCchannelEnergy(Int_t module, Int_t column, Int_t row, Int_t gain) const
{
  //Return EMC calibration coefficient
  //module, column,raw should follow the internal PHOS convention:
  //module 1:5, column 1:56, row 1:64

  return fADCchannelEnergy[module-1][column-1][row-1][gain];
}

//________________________________________________________________
Float_t AliHLTPHOSEmcCalibData::GetADCpedestalEmcMeasured(Int_t module, Int_t column, Int_t row, Int_t gain) const
{
  //Return EMC pedestal
  //module, column,raw should follow the internal PHOS convention:
  //module 1:5, column 1:56, row 1:64


  return fADCpedestalEmcMeasured[module-1][column-1][row-1][gain];
}

//________________________________________________________________
void AliHLTPHOSEmcCalibData::SetADCchannelEnergy(Int_t module, Int_t column, Int_t row, Int_t gain, Float_t value)
{
  //Set EMC calibration coefficient
  //module, column,raw should follow the internal PHOS convention:
  //module 1:5, column 1:56, row 1:64

  fADCchannelEnergy[module-1][column-1][row-1][gain] = value;
}

//________________________________________________________________
void AliHLTPHOSEmcCalibData::SetADCpedestalEmcMeasured(Int_t module, Int_t column, Int_t row, Int_t gain, Float_t value)
{
  //Set EMC pedestal
  //module, column,raw should follow the internal PHOS convention:
  //module 1:5, column 1:56, row 1:64
  fADCpedestalEmcMeasured[module-1][column-1][row-1][gain] = value;
}
Commit	Line	Data
1b41ab20	1	// $Id$
1b41ab20	2
228cb536	3	/**************************************************************************
	4	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* *
	6	* Author: The ALICE Off-line Project. *
	7	* Contributors are mentioned in the code where appropriate. *
	8	* *
	9	* Permission to use, copy, modify and distribute this software and its *
	10	* documentation strictly for non-commercial purposes is hereby granted *
	11	* without fee, provided that the above copyright notice appears in all *
	12	* copies and that both the copyright notice and this permission notice *
	13	* appear in the supporting documentation. The authors make no claims *
	14	* about the suitability of this software for any purpose. It is *
	15	* provided "as is" without express or implied warranty. *
	16	* *
	17	* 7/1- 2007 *
	18	* Author: Per Thomas Hille *
	19	* code modified to include the two differnt gains of each detection *
	20	* channel and also to include the Pedestals loaded into the ALTRO *
	21	* registers for basline subtraction and zero supression *
	22	* *
	23	* *
	24	**************************************************************************/
	25
	26	///////////////////////////////////////////////////////////////////////////////
	27	// //
	28	// class for PHOS EmCal calibration //
	29	// //
	30	///////////////////////////////////////////////////////////////////////////////
	31
	32	#include "AliHLTPHOSEmcCalibData.h"
	33
	34	ClassImp(AliHLTPHOSEmcCalibData)
	35
	36	//________________________________________________________________
	37	AliHLTPHOSEmcCalibData::AliHLTPHOSEmcCalibData(): TNamed()
	38	{
	39	// Default constructor
	40	Reset();
	41	}
	42
	43	//________________________________________________________________
	44	AliHLTPHOSEmcCalibData::AliHLTPHOSEmcCalibData(const char* name)
	45	{
	46	// Constructor
	47	TString namst = "Calib_";
	48	namst += name;
	49	SetName(namst.Data());
	50	SetTitle(namst.Data());
	51	Reset();
	52	}
	53
	54	//________________________________________________________________
	55	AliHLTPHOSEmcCalibData::AliHLTPHOSEmcCalibData(const AliHLTPHOSEmcCalibData& calibda) :
	56	TNamed(calibda)
	57	{
	58	// copy constructor
	59	SetName(calibda.GetName());
	60	SetTitle(calibda.GetName());
	61
	62
bde48b84	63	// for(Int_t gain = 0; gain < N_GAINS; gain ++){
27029341	64	for(Int_t module=0; module<NMODULES; module++) {
	65	for(Int_t column=0; column< NXCOLUMNSMOD; column++) {
	66	for(Int_t row=0; row<NZROWSMOD; row++) {
	67	for(Int_t gain = 0; gain < NGAINS; gain ++){
228cb536	68	fADCchannelEnergy[module][column][row][gain] = calibda.fADCchannelEnergy[module][column][row][gain];
	69	fADCpedestalEmcMeasured[module][column][row][gain] = calibda.fADCpedestalEmcMeasured[module][column][row][gain];
	70	}
	71	}
	72	}
	73	}
	74	}
	75
	76	//________________________________________________________________
	77	AliHLTPHOSEmcCalibData &AliHLTPHOSEmcCalibData::operator =(const AliHLTPHOSEmcCalibData& calibda)
	78	{
	79	// assignment operator
	80
	81	if(this != &calibda) {
	82
	83	SetName(calibda.GetName());
	84	SetTitle(calibda.GetName());
bde48b84	85	// for(Int_t gain = 0; gain < N_GAINS; gain ++){
27029341	86	for(Int_t module=0; module<NMODULES; module++) {
	87	for(Int_t column=0; column< NXCOLUMNSMOD; column++) {
	88	for(Int_t row=0; row<NZROWSMOD; row++) {
	89	for(Int_t gain = 0; gain < NGAINS; gain ++){
228cb536	90	fADCchannelEnergy[module][column][row][gain] = calibda.fADCchannelEnergy[module][column][row][gain];
	91	fADCpedestalEmcMeasured[module][column][row][gain] = calibda.fADCpedestalEmcMeasured[module][column][row][gain];
	92	}
	93	}
	94	}
	95	}
	96	}
	97	return *this;
	98	}
	99
	100	//________________________________________________________________
	101	AliHLTPHOSEmcCalibData::~AliHLTPHOSEmcCalibData()
	102	{
	103	// Destructor
	104	}
	105
	106	//________________________________________________________________
	107	void AliHLTPHOSEmcCalibData::Reset()
	108	{
	109	// Set all pedestals and all ADC channels to its ideal values = 1.
	110
	111	// for(Int_t gain = 0; gain < PHOS_GAINS; gain ++){
27029341	112	for (Int_t module=0; module<NMODULES; module++){
	113	for (Int_t column=0; column< NXCOLUMNSMOD; column++){
	114	for (Int_t row=0; row<NZROWSMOD; row++){
	115	for(Int_t gain = 0; gain < NGAINS; gain ++){
228cb536	116	fADCpedestalEmcMeasured[module][column][row][gain] = 0.;
	117	fADCchannelEnergy[module][column][row][gain] = 1.;
	118	}
	119	}
	120	}
	121	}
	122	}
	123
	124	//________________________________________________________________
	125	void AliHLTPHOSEmcCalibData::Print(Option_t *option) const
	126	{
	127	// Print tables of pedestals and ADC channels
	128
	129	if (strstr(option,"ped")) {
	130	printf("\n ---- EMC Pedestal values ----\n\n");
bde48b84	131	// for(Int_t gain = 0; gain < N_GAINS; gain ++){
27029341	132	for (Int_t module=0; module<NMODULES; module++){
228cb536	133	printf("============== Module %d\n",module+1);
27029341	134	for (Int_t column=0; column< NXCOLUMNSMOD; column++){
	135	for (Int_t row=0; row<NZROWSMOD; row++){
	136	for(Int_t gain = 0; gain < NGAINS; gain ++){
228cb536	137	printf("%4.1f",fADCpedestalEmcMeasured[module][column][row][gain]);
	138	}
	139	printf("\n");
	140	}
	141	}
	142	}
	143	}
	144
	145	if (strstr(option,"gain")) {
	146	printf("\n ---- EMC ADC channel values ----\n\n");
27029341	147	// for(Int_t gain = 0; gain < NGAINS; gain ++){
27029341	148	for (Int_t module=0; module<NMODULES; module++){
228cb536	149	printf("============== Module %d\n",module+1);
27029341	150	for (Int_t column=0; column< NXCOLUMNSMOD; column++){
	151	for (Int_t row=0; row<NZROWSMOD; row++){
	152	for(Int_t gain = 0; gain < NGAINS; gain ++){
228cb536	153	printf("%4.1f",fADCchannelEnergy[module][column][row][gain]);
	154	}
	155	printf("\n");
	156	}
	157	}
	158	}
	159	}
	160	}
	161
d504c864	162	//________________________________________________________________
228cb536	163	void
	164	AliHLTPHOSEmcCalibData::MakeADCpedestalCorrectionTable()
	165	{
27029341	166	for (Int_t module=0; module<NMODULES; module++){
228cb536	167	printf("============== Module %d\n",module+1);
27029341	168	for (Int_t column=0; column< NXCOLUMNSMOD; column++){
	169	for (Int_t row=0; row<NZROWSMOD; row++){
	170	for(Int_t gain = 0; gain < NGAINS; gain ++){
228cb536	171	fADCpedestalCorrectionTable[module][column][row][gain] = fADCpedestalEmcMeasured[module][column][row][gain] - fADCpedestalAltroReg[module][column][row][gain];
	172	// printf("%4.1f",fADCchannelEnergy[module][column][row][gain]);
	173	//
	174	}
	175	printf("\n");
	176	}
	177	}
	178	}
	179	}
	180
	181
	182	//________________________________________________________________
	183	Float_t AliHLTPHOSEmcCalibData::GetADCchannelEnergy(Int_t module, Int_t column, Int_t row, Int_t gain) const
	184	{
	185	//Return EMC calibration coefficient
	186	//module, column,raw should follow the internal PHOS convention:
	187	//module 1:5, column 1:56, row 1:64
	188
	189	return fADCchannelEnergy[module-1][column-1][row-1][gain];
	190	}
	191
	192	//________________________________________________________________
	193	Float_t AliHLTPHOSEmcCalibData::GetADCpedestalEmcMeasured(Int_t module, Int_t column, Int_t row, Int_t gain) const
	194	{
	195	//Return EMC pedestal
	196	//module, column,raw should follow the internal PHOS convention:
	197	//module 1:5, column 1:56, row 1:64
	198
	199
	200	return fADCpedestalEmcMeasured[module-1][column-1][row-1][gain];
	201	}
	202
	203	//________________________________________________________________
	204	void AliHLTPHOSEmcCalibData::SetADCchannelEnergy(Int_t module, Int_t column, Int_t row, Int_t gain, Float_t value)
	205	{
	206	//Set EMC calibration coefficient
	207	//module, column,raw should follow the internal PHOS convention:
	208	//module 1:5, column 1:56, row 1:64
	209
	210	fADCchannelEnergy[module-1][column-1][row-1][gain] = value;
	211	}
	212
	213	//________________________________________________________________
	214	void AliHLTPHOSEmcCalibData::SetADCpedestalEmcMeasured(Int_t module, Int_t column, Int_t row, Int_t gain, Float_t value)
	215	{
	216	//Set EMC pedestal
	217	//module, column,raw should follow the internal PHOS convention:
	218	//module 1:5, column 1:56, row 1:64
	219	fADCpedestalEmcMeasured[module-1][column-1][row-1][gain] = value;
	220	}