///*-- Author: Yves Schutz (SUBATECH)
// : Aleksei Pavlinov (WSU); Jun 30, 2006 - ALICE numbering scheme
// : Add decalibration and time calibration arrays: Jul 21, 2011 (GCB)
+// : adapted for DCAL by M.L. Wang CCNU & Subatech Oct-18-2012
///////////////////////////////////////////////////////////////////////////////
// //
// class for EMCAL calibration //
ClassImp(AliEMCALCalibData)
-//________________________________________________________________
-AliEMCALCalibData::AliEMCALCalibData()
+//__________________________________________________________
+AliEMCALCalibData::AliEMCALCalibData() :
+TNamed(), fADCchannelRef(0)
{
// Default constructor
Reset();
}
-//________________________________________________________________
-AliEMCALCalibData::AliEMCALCalibData(const char* name)
+//________________________________________________________________________
+AliEMCALCalibData::AliEMCALCalibData(const char* name) :
+TNamed(name,name),fADCchannelRef(0)
{
// Constructor
- TString namst = "Calib_";
- namst += name;
- SetName(namst.Data());
- SetTitle(namst.Data());
Reset();
}
-//________________________________________________________________
+//______________________________________________________________________
AliEMCALCalibData::AliEMCALCalibData(const AliEMCALCalibData& calibda) :
- TNamed(calibda)
+TNamed(calibda), fADCchannelRef(calibda.fADCchannelRef)
{
// copy constructor
- SetName(calibda.GetName());
+ SetName (calibda.GetName());
SetTitle(calibda.GetName());
Reset();
-
- Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
+
+ Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //22
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 supermodule = 0; supermodule < nSMod; supermodule++) {
+ nCol = AliEMCALGeoParams::fgkEMCALCols; //48
+ nRow = AliEMCALGeoParams::fgkEMCALRows; //24
+ // in reality they are 1/3 but leave them as 1/2
+
+ if(supermodule /2 == 5)
+ nRow = nRow/2;
+ if(supermodule > 11 && supermodule < 18)
+ nCol = nCol*2/3;
for(Int_t column = 0; column<nCol; column++) {
fADCpedestal[supermodule][column][row] =
calibda.GetADCpedestal(supermodule,column,row);
- fTimeChannel[supermodule][column][row] =
- calibda.GetTimeChannel(supermodule,column,row);
-
fTimeChannelDecal[supermodule][column][row] =
calibda.GetTimeChannelDecal(supermodule,column,row);
+ for(Int_t bc = 0; bc < 4; bc++)
+ fTimeChannel[supermodule][column][row][bc] =
+ calibda.GetTimeChannel(supermodule,column,row,bc);
}
}
}
}
-
//________________________________________________________________
AliEMCALCalibData &AliEMCALCalibData::operator =(const AliEMCALCalibData& calibda)
{
// assignment operator
- SetName(calibda.GetName());
+ SetName (calibda.GetName());
SetTitle(calibda.GetName());
Reset();
-
- Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
+
+ fADCchannelRef = calibda.GetADCchannelRef() ;
+
+ Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //22
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 supermodule = 0; supermodule < nSMod; supermodule++) {
+ nCol = AliEMCALGeoParams::fgkEMCALCols; //48
+ nRow = AliEMCALGeoParams::fgkEMCALRows; //24
+ // in reality they are 1/3 but leave them as 1/2
+
+ if(supermodule /2 == 5)
+ nRow = nRow/2;
+ if(supermodule > 11 && supermodule < 18)
+ nCol = nCol*2/3;
for(Int_t column = 0; column<nCol; column++) {
fADCpedestal[supermodule][column][row] =
calibda.GetADCpedestal(supermodule,column,row);
- fTimeChannel[supermodule][column][row] =
- calibda.GetTimeChannel(supermodule,column,row);
-
fTimeChannelDecal[supermodule][column][row] =
calibda.GetTimeChannelDecal(supermodule,column,row);
+ for(Int_t bc = 0; bc < 4; bc++)
+ fTimeChannel[supermodule][column][row][bc] =
+ calibda.GetTimeChannel(supermodule,column,row,bc);
+
}
}
}
+
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
+
+ fADCchannelRef = 0.0162;
+
+ Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //22
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 supermodule=0; supermodule<nSMod; supermodule++){
+ nCol = AliEMCALGeoParams::fgkEMCALCols; //48
+ nRow = AliEMCALGeoParams::fgkEMCALRows; //24
+ // in reality they are 1/3 but leave them as 1/2
+
+ if(supermodule /2 == 5)
+ nRow = nRow/2;
+ if(supermodule > 11 && supermodule < 18)
+ nCol = nCol*2/3;
+
for (Int_t column=0; column < nCol; column++){
for (Int_t row = 0; row < nRow; row++){
fADCchannelDecal [supermodule][column][row]=1.;
fADCchannel [supermodule][column][row]=1.;
- fTimeChannel [supermodule][column][row]=0.;
fTimeChannelDecal[supermodule][column][row]=0.;
-
+
+ for(Int_t bc = 0; bc < 4; bc++)
+ fTimeChannel[supermodule][column][row][bc]=0;
+
}
}
- }
+ }
}
-//________________________________________________________________
+//____________________________________________________
void AliEMCALCalibData::Print(Option_t *option) const
{
// Print tables of pedestals and ADC channels widths
// options are: "gain", "ped", "decal", "time", "all"
- Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
+ Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //22
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") || strstr(option,"all")) {
printf("\n ---- Pedestal values ----\n\n");
for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
- if(supermodule >= 10)
- nRow = nRow2;
+ nCol = AliEMCALGeoParams::fgkEMCALCols; //48
+ nRow = AliEMCALGeoParams::fgkEMCALRows; //24
+ // in reality they are 1/3 but leave them as 1/2
+ if(supermodule /2 == 5)
+ nRow = nRow/2;
+ if(supermodule > 11 && supermodule < 18)
+ nCol = nCol*2/3;
printf("============== Supermodule %d\n",supermodule+1);
for (Int_t column=0; column<nCol; column++){
for (Int_t row=0; row<nRow; row++){
if (strstr(option,"gain") || strstr(option,"all")) {
printf("\n ---- ADC channel values ----\n\n");
for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
- if(supermodule >= 10)
- nRow = nRow2;
+ nCol = AliEMCALGeoParams::fgkEMCALCols; //48
+ nRow = AliEMCALGeoParams::fgkEMCALRows; //24
+
+ // in reality they are 1/3 but leave them as 1/2
+ if(supermodule /2 == 5)
+ nRow = nRow/2;
+ if(supermodule > 11 && supermodule < 18)
+ nCol = nCol*2/3;
printf("============== Supermodule %d\n",supermodule+1);
for (Int_t column=0; column<nCol; column++){
for (Int_t row=0; row<nRow; row++){
if (strstr(option,"adcdecal") || strstr(option,"all")) {
printf("\n ---- ADC decalibration channel values ----\n\n");
for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
- if(supermodule >= 10)
- nRow = nRow2;
+ nCol = AliEMCALGeoParams::fgkEMCALCols; //48
+ nRow = AliEMCALGeoParams::fgkEMCALRows; //24
+ // in reality they are 1/3 but leave them as 1/2
+ if(supermodule /2 == 5)
+ nRow = nRow/2;
+ if(supermodule > 11 && supermodule < 18)
+ nCol = nCol*2/3;
printf("============== Supermodule %d\n",supermodule+1);
for (Int_t column=0; column<nCol; column++){
for (Int_t row=0; row<nRow; row++){
if (strstr(option,"time") || strstr(option,"all")) {
printf("\n ---- time channel values ----\n\n");
for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
- if(supermodule >= 10)
- nRow = nRow2;
+ nCol = AliEMCALGeoParams::fgkEMCALCols; //48
+ nRow = AliEMCALGeoParams::fgkEMCALRows; //24
+ // in reality they are 1/3 but leave them as 1/2
+ if(supermodule /2 == 5)
+ nRow = nRow/2;
+ if(supermodule > 11 && supermodule < 18)
+ nCol = nCol*2/3;
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 ",fTimeChannel[supermodule][column][row]);
+ for(Int_t bc = 0; bc < 4; bc++)
+ printf(" %2.4f ",fTimeChannel[supermodule][column][row][bc]);
}
printf("\n");
}
}
}
-
if (strstr(option,"time") || strstr(option,"all")) {
printf("\n ---- time decalibration channel values ----\n\n");
for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
- if(supermodule >= 10)
- nRow = nRow2;
+ nCol = AliEMCALGeoParams::fgkEMCALCols; //48
+ nRow = AliEMCALGeoParams::fgkEMCALRows; //24
+ // in reality they are 1/3 but leave them as 1/2
+ if(supermodule /2 == 5)
+ nRow = nRow/2;
+ if(supermodule > 11 && supermodule < 18)
+ nCol = nCol*2/3;
printf("============== Supermodule %d\n",supermodule+1);
for (Int_t column=0; column<nCol; column++){
for (Int_t row=0; row<nRow; row++){
}
}
}
-
-
}
//________________________________________________________________
// 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];
+ return fADCpedestal[supermodule][column][row];
}
//________________________________________________________________
-Float_t AliEMCALCalibData::GetTimeChannel(Int_t supermodule, Int_t column, Int_t row) const
+Float_t AliEMCALCalibData::GetTimeChannel(Int_t supermodule, Int_t column, Int_t row, Int_t bc) const
{
// Set channel time witdth values
- return fTimeChannel[supermodule][column][row];
+ return fTimeChannel[supermodule][column][row][bc];
}
//________________________________________________________________
}
//________________________________________________________________
-void AliEMCALCalibData::SetTimeChannel(Int_t supermodule, Int_t column, Int_t row, Float_t value)
+void AliEMCALCalibData::SetTimeChannel(Int_t supermodule, Int_t column, Int_t row, Int_t bc, Float_t value)
{
// Set ADC pedestal values
- fTimeChannel[supermodule][column][row] = value;
+ fTimeChannel[supermodule][column][row][bc] = value;
}
//________________________________________________________________
}
-//________________________________________________________________
-void AliEMCALCalibData::Fill(const AliEMCALCalibData *cd1, const AliEMCALCalibData *cd2, Bool_t print)
-{
- // Fill for (relative) recalibration to undo calibration
- // from 1 and apply calibration from 2.
-
- 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++){
- Double_t adc1 = cd1->GetADCchannel(supermodule, column, row);
- Double_t ped1 = cd1->GetADCpedestal(supermodule, column, row);
- Double_t adc2 = cd2->GetADCchannel(supermodule, column, row);
- Double_t ped2 = cd2->GetADCpedestal(supermodule, column, row);
- Double_t adc = adc2/adc1;
- Double_t ped = ped2-adc*ped1;
- if (print)
- printf("%d %d %d: %f %f %f %f -> %f %f\n",supermodule,column,row,adc1,ped1,adc2,ped2,adc,ped);
- fADCchannel [supermodule][column][row]=adc;
- fADCpedestal[supermodule][column][row]=ped;
- }
- }
- }
-}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff