[u/mrichter/AliRoot.git] / EMCAL / macros / GeneratePedestalScript.C

// some global var/constants
const Int_t kNSM = 4; // for first LHC run
const Int_t kNRCU = 2;
AliCaloAltroMapping *fMapping[4]; // 1 for each side (A/C) and each RCU (0/1), i.e. 2*2 total
const Int_t kNBranch = 2;
const Int_t kNFEC = 10; // 0..9, when including LED Ref
const Int_t kNChip = 5; // really 0,2..4, i.e. skip #1
const Int_t kNChan = 16;
Float_t fMeanPed[kNSM][kNRCU][kNBranch][kNFEC][kNChip][kNChan];
Float_t fRmsPed[kNSM][kNRCU][kNBranch][kNFEC][kNChip][kNChan];
//
const int kNStrips = 24; // per SM
Int_t fHWAddrLEDRef[kNStrips][2]; // [2] is for Low/High gain

const Bool_t kDebug = kFALSE;
const Float_t kBadRMS = 20;

// help methods
void GetPedVal(const Int_t iSM, const Int_t igain, const TProfile2D *h2);
void GetPedValLEDRef(const Int_t iSM, const Int_t igain, const TProfile *h);
void PrintScript();
void Clear();
Int_t GetHWAddress(Int_t iside, Int_t icol, Int_t irow, Int_t igain);
Int_t GetHWAddressLEDRef(Int_t istrip, Int_t igain);
void DecodeHWAddress(Int_t hwAddr, Int_t & branch, Int_t & FEC, Int_t & chip, Int_t & chan);
void GetMapping();
void CreateMappingLEDRef();

// main method
void 
//GeneratePedestalScript(const char * filename = "alien/Run113790_113790_v1_s0.root") // 1st set
GeneratePedestalScript(const char * filename = "Run117756_117756_v1_s0.root") // 2nd set
{
  // get the DA info/object
  TFile *file = TFile::Open(filename); 
  //TMP  AliCaloCalibPedestal *emcCalibPedestal = AliCDBEntry->GetObject();
  if (kDebug) {
    file->ls();
  }

  // Get mapping file info, and clear arrays
  Clear();
  GetMapping();
  CreateMappingLEDRef();

  // Store the pedestal info
  for (Int_t iSM=0; iSM<kNSM; iSM++) {
    GetPedVal( iSM, 0, emcCalibPedestal->GetPedProfileLowGain(iSM) );
    GetPedVal( iSM, 1, emcCalibPedestal->GetPedProfileHighGain(iSM) );
    GetPedValLEDRef( iSM, 0, emcCalibPedestal->GetPedLEDRefProfileLowGain(iSM) );
    GetPedValLEDRef( iSM, 1, emcCalibPedestal->GetPedLEDRefProfileHighGain(iSM) );
  }

  // Generate the needed scripts
  PrintScript();

}

void 
GetPedVal(const Int_t iSM, const Int_t igain, const TProfile2D *h2)
{
  Int_t isect = iSM / 2; //
  Int_t iside = iSM % 2; // A or C side
  Int_t nCols = h2->GetNbinsX();
  Int_t nRows = h2->GetNbinsY();
  if (kDebug) {
    printf("GetPedVal: iSM %d isect %d iside %d igain %d nRows %d nCols %d\n",
	   iSM, isect, iside, igain, nRows, nCols);
  }
  Int_t hwAddress = 0;
  Int_t iRCU = 0;
  Int_t branch = 0;
  Int_t FEC = 0;
  Int_t chip = 0;
  Int_t chan = 0;

  Int_t icol = 0;
  Int_t irow = 0;
  Int_t bin = 0;

  for (icol=0; icol<nCols; icol++) {
    for (irow=0; irow<nRows; irow++) {

      hwAddress = GetHWAddress(iside, icol, irow, igain, iRCU);
      DecodeHWAddress(hwAddress, branch, FEC, chip, chan);  
      bin = h2->FindBin(icol, irow);
            
      // store the values
      fMeanPed[iSM][iRCU][branch][FEC][chip][chan] = h2->GetBinContent(bin);
      fRmsPed[iSM][iRCU][branch][FEC][chip][chan] = h2->GetBinError(bin);

      // report bad RMS channels:
      if (h2->GetBinError(bin) > kBadRMS) {
	printf(" bad pedestal RMS: iSM %d icol %d irow %d igain %d iRCU %d branch %d FEC %d chip %d chan %d - mean %4.1f rms %4.1f\n", 
	       iSM, icol, irow, igain, iRCU, branch, FEC, chip, chan,
	       h2->GetBinContent(bin), h2->GetBinError(bin));
      }
    }
  }

  return;
}


void 
GetPedValLEDRef(const Int_t iSM, const Int_t igain, const TProfile *h)
{
  Int_t isect = iSM / 2; //
  Int_t iside = iSM % 2; // A or C side
  Int_t nStrips = h->GetNbinsX();

  if (kDebug) {
    printf("GetPedValLEDRef: iSM %d isect %d iside %d igain %d nStrips %d\n",
	   iSM, isect, iside, igain, nStrips);
  }
  Int_t hwAddress = 0;
  Int_t iRCU = 0; // always true for LED Ref FEE
  Int_t branch = 0;
  Int_t FEC = 0;
  Int_t chip = 0;
  Int_t chan = 0;

  Int_t icol = 0;
  Int_t irow = 0;
  Int_t bin = 0;

  for (int istrip=0; istrip<nStrips; istrip++) {

    hwAddress = GetHWAddressLEDRef(istrip, igain);
    DecodeHWAddress(hwAddress, branch, FEC, chip, chan);  
    bin = h->FindBin(istrip);
            
    // store the values
    fMeanPed[iSM][iRCU][branch][FEC][chip][chan] = h->GetBinContent(bin);
    fRmsPed[iSM][iRCU][branch][FEC][chip][chan] = h->GetBinError(bin);
  }

  return;
}

void 
PrintScript()
{
  const char * sideStr[] = {"A","C"};
  const char * branchStr[] = {"A","B"};
  int VFPED = 0x06;
  int RCUWrite = 0x200000;

  char filename[100];
  char scriptLine[200];

  Int_t iSM = 0;
  Int_t iRCU = 0;
  Int_t ibranch = 0;
  Int_t iFEC = 0;
  Int_t ichip = 0;
  Int_t ichan = 0;
  Int_t Ped = 0;

  for (iSM=0; iSM<kNSM; iSM++) {
    int iside = iSM % 2;
    int isect = iSM / 2;
    for (iRCU=0; iRCU<kNRCU; iRCU++) {
      sprintf(filename, "setSM%1s%dRCU%d.scr", 
	      sideStr[iside], isect, iRCU);
      ofstream fout(filename);
      int nscriptLines = 0;

      for (ibranch=0; ibranch<kNBranch; ibranch++) {
	int firstFEC = 1; 
	if (ibranch==0 && iRCU==0) { // extra LED Ref FEE in this location
	  firstFEC = 0;
	}

	for (iFEC=firstFEC; iFEC<kNFEC; iFEC++) { // FEC 1..9 (or 0..9)
	  for (ichip=0; ichip<kNChip; ichip++) { // ALTRO 0,2,3,4
	    if (ichip!=1) {
	      for (ichan=0; ichan<kNChan; ichan++) {

		if (iFEC!=0 || (ichan<8 || ichan>11)) {

		  Ped = TMath::Nint(fMeanPed[iSM][iRCU][ibranch][iFEC][ichip][ichan]);
		  // raise Ped value to max for channels with exceptionally large RMS
		  if (fRmsPed[iSM][iRCU][ibranch][iFEC][ichip][ichan] > kBadRMS) {
		    printf(" bad pedestal RMS: iSM %d iRCU %d ibranch %d iFEC %d ichip %d ichan %d - raising from %d to 0x3ff\n", 
			   iSM, iRCU, ibranch, iFEC, ichip, ichan, Ped);
		    Ped = 0x3ff;
		  }
		  // 
		  int writeAddr = (ibranch << 16) | (iFEC << 12) | (ichip << 9) |
		    (ichan << 5) | VFPED | RCUWrite;  
		  sprintf(scriptLine, "w 0x%04x 0x%06x  # Branch %s, Card %d, Altro %d, Chan %d", 
			  nscriptLines, writeAddr, branchStr[ibranch],
			  iFEC, ichip, ichan);
		  fout << scriptLine << endl;
		  nscriptLines++;
		  
		  int writeVal = (Ped | RCUWrite);
		  sprintf(scriptLine, "w 0x%04x 0x%06x  # Pedestal 0x%x = %d", 
			  nscriptLines, writeVal, Ped, Ped);
		  fout << scriptLine << endl;
		  nscriptLines++;
		  
		}
	      }
	    }
	  }
	}
      }
      // ending, with execute and update step..	    
      sprintf(scriptLine, "w 0x%04x 0x380000  # End of the sequence", 
	      nscriptLines);
      fout << scriptLine << endl;
      nscriptLines++;

      sprintf(scriptLine, "w 0x%04x 0x3F0000  # End of the instruction memory", 
	      nscriptLines);
      fout << scriptLine << endl;
      nscriptLines++;

      fout << "wait 100 us" << endl;
      fout << "w 0x5304 0x0       \# execute and update registers" << endl; 
      fout.close();
    } // iRCU
  }// iSM

  return;
}

void 
Clear()
{
  for (Int_t iSM=0; iSM<kNSM; iSM++) {
    for (Int_t iRCU=0; iRCU<kNRCU; iRCU++) {
      for (Int_t ibranch=0; ibranch<kNBranch; ibranch++) {
	for (Int_t iFEC=0; iFEC<kNFEC; iFEC++) {
	  for (Int_t ichip=0; ichip<kNChip; ichip++) {
	    for (Int_t ichan=0; ichan<kNChan; ichan++) {
	      fMeanPed[iSM][iRCU][ibranch][iFEC][ichip][ichan] = 0;
	      fRmsPed[iSM][iRCU][ibranch][iFEC][ichip][ichan] = 0;
	    }
	  }
	}
      }
    }
  }

  for (int istrip=0; istrip<kNStrips; istrip++) {
    fHWAddrLEDRef[istrip][0] = 0; 
    fHWAddrLEDRef[istrip][1] = 0; 
  }

  return;
}

void 
DecodeHWAddress(Int_t hwAddr, Int_t & branch, Int_t & FEC, Int_t & chip, Int_t & chan)
{
  chan = hwAddr & 0xf;
  chip = (hwAddr >> 4) & 0x7;
  FEC = (hwAddr >> 7) & 0xf;
  branch = (hwAddr >> 11) & 0x1;
  return;
}

Int_t 
GetHWAddress(Int_t iside, Int_t icol, Int_t irow, Int_t igain, Int_t & iRCU)
{
  iRCU = -111;

  //RCU0
  if (0<=irow&&irow<8) iRCU=0; // first cable row
  else if (8<=irow&&irow<16 && 0<=icol&&icol<24) iRCU=0; // first half; 
  //second cable row
  //RCU1
  else if(8<=irow&&irow<16 && 24<=icol&&icol<48) iRCU=1; // second half; 
  //second cable row
  else if(16<=irow&&irow<24) iRCU=1; // third cable row

  // swap for odd=C side, to allow us to cable both sides the same
  Int_t iRCUSide = iRCU; 
 if (iside == 1) {
    iRCU = 1 - iRCU;
    iRCUSide = iRCU + 2; // to make it map file index
  }
  Int_t hwAddress = fMapping[iRCUSide]->GetHWAddress(irow, icol, igain);

  return hwAddress;
}

Int_t 
GetHWAddressLEDRef(Int_t istrip, Int_t igain)
{
  Int_t iRCU = 0; // for both sides; LED ref info is the same for both sides
  Int_t caloflag = 3; // AliCaloRawStreamV3::kLEDMonData;
 
  Int_t hwAddress = fHWAddrLEDRef[istrip][igain];

  return hwAddress;
}

void 
GetMapping()
{  
  TString sides[]={"A","C"};
  // Read mapping files from $ALICE_ROOT/CALO/mapping/*.data
  TString path = gSystem->Getenv("ALICE_ROOT");
  path += "/EMCAL/mapping/RCU";
  TString path2;
  for(Int_t j = 0; j < 2; j++){ // sides
    for(Int_t i = 0; i < 2; i++) { // RCU
      path2 = path;
      path2 += i;
      path2 += sides[j];
      path2 += ".data";
      if (kDebug) { printf("Mapping file: %s\n",path2.Data()); }
      fMapping[j*2 + i] = new AliCaloAltroMapping(path2.Data());
    }
  }
  return;
}

void 
CreateMappingLEDRef()
{ 
  Int_t iRCU = 0; // for both sides; LED ref info is the same for both sides
  Int_t caloflag = 3; // AliCaloRawStreamV3::kLEDMonData;

  Int_t maxAddr = 1 << 7; // LED Ref FEE is in FEC pos 0, i.e. addr space 0..127

  int nLEDRefFEEChan = 0;

  Int_t branch = 0;
  Int_t FEC = 0;
  Int_t chip = 0;
  Int_t chan = 0;
  for (int hwaddr = 0; hwaddr<maxAddr; hwaddr++) {

    DecodeHWAddress(hwaddr, branch, FEC, chip, chan);  
    if ( (chip!=1 && chip<kNChip) &&  // ALTROs 0,2,3,4
	 (chan<8 || chan>11) ) { // actual installed LED Ref FEE channels

      int istrip = fMapping[iRCU]->GetPad(hwaddr);
      int igain = fMapping[iRCU]->GetPadRow(hwaddr);
      int iflag = fMapping[iRCU]->GetSector(hwaddr);
      if (iflag == caloflag) {
	fHWAddrLEDRef[istrip][igain] = hwaddr; 
	nLEDRefFEEChan++;
      }
    }
  }

  if (kDebug) { cout << " nLEDRefFEEChan " << nLEDRefFEEChan << endl; }
}
Commit	Line	Data
29f94584	1	// some global var/constants
	2	const Int_t kNSM = 4; // for first LHC run
	3	const Int_t kNRCU = 2;
	4	AliCaloAltroMapping fMapping[4]; // 1 for each side (A/C) and each RCU (0/1), i.e. 22 total
	5	const Int_t kNBranch = 2;
	6	const Int_t kNFEC = 10; // 0..9, when including LED Ref
	7	const Int_t kNChip = 5; // really 0,2..4, i.e. skip #1
	8	const Int_t kNChan = 16;
	9	Float_t fMeanPed[kNSM][kNRCU][kNBranch][kNFEC][kNChip][kNChan];
	10	Float_t fRmsPed[kNSM][kNRCU][kNBranch][kNFEC][kNChip][kNChan];
	11	//
	12	const int kNStrips = 24; // per SM
	13	Int_t fHWAddrLEDRef[kNStrips][2]; // [2] is for Low/High gain
	14
75dbb105	15	const Bool_t kDebug = kFALSE;
	16	const Float_t kBadRMS = 20;
	17
29f94584	18	// help methods
	19	void GetPedVal(const Int_t iSM, const Int_t igain, const TProfile2D *h2);
	20	void GetPedValLEDRef(const Int_t iSM, const Int_t igain, const TProfile *h);
	21	void PrintScript();
	22	void Clear();
	23	Int_t GetHWAddress(Int_t iside, Int_t icol, Int_t irow, Int_t igain);
	24	Int_t GetHWAddressLEDRef(Int_t istrip, Int_t igain);
	25	void DecodeHWAddress(Int_t hwAddr, Int_t & branch, Int_t & FEC, Int_t & chip, Int_t & chan);
	26	void GetMapping();
	27	void CreateMappingLEDRef();
	28
	29	// main method
	30	void
75dbb105	31	//GeneratePedestalScript(const char * filename = "alien/Run113790_113790_v1_s0.root") // 1st set
75dbb105	32	GeneratePedestalScript(const char * filename = "Run117756_117756_v1_s0.root") // 2nd set
29f94584	33	{
	34	// get the DA info/object
	35	TFile *file = TFile::Open(filename);
75dbb105	36	//TMP AliCaloCalibPedestal *emcCalibPedestal = AliCDBEntry->GetObject();
	37	if (kDebug) {
	38	file->ls();
	39	}
29f94584	40
	41	// Get mapping file info, and clear arrays
	42	Clear();
	43	GetMapping();
	44	CreateMappingLEDRef();
	45
	46	// Store the pedestal info
	47	for (Int_t iSM=0; iSM<kNSM; iSM++) {
	48	GetPedVal( iSM, 0, emcCalibPedestal->GetPedProfileLowGain(iSM) );
	49	GetPedVal( iSM, 1, emcCalibPedestal->GetPedProfileHighGain(iSM) );
	50	GetPedValLEDRef( iSM, 0, emcCalibPedestal->GetPedLEDRefProfileLowGain(iSM) );
	51	GetPedValLEDRef( iSM, 1, emcCalibPedestal->GetPedLEDRefProfileHighGain(iSM) );
	52	}
	53
	54	// Generate the needed scripts
	55	PrintScript();
	56
	57	}
	58
	59	void
	60	GetPedVal(const Int_t iSM, const Int_t igain, const TProfile2D *h2)
	61	{
	62	Int_t isect = iSM / 2; //
	63	Int_t iside = iSM % 2; // A or C side
	64	Int_t nCols = h2->GetNbinsX();
	65	Int_t nRows = h2->GetNbinsY();
75dbb105	66	if (kDebug) {
	67	printf("GetPedVal: iSM %d isect %d iside %d igain %d nRows %d nCols %d\n",
	68	iSM, isect, iside, igain, nRows, nCols);
	69	}
29f94584	70	Int_t hwAddress = 0;
	71	Int_t iRCU = 0;
	72	Int_t branch = 0;
	73	Int_t FEC = 0;
	74	Int_t chip = 0;
	75	Int_t chan = 0;
	76
	77	Int_t icol = 0;
	78	Int_t irow = 0;
	79	Int_t bin = 0;
	80
	81	for (icol=0; icol<nCols; icol++) {
	82	for (irow=0; irow<nRows; irow++) {
	83
	84	hwAddress = GetHWAddress(iside, icol, irow, igain, iRCU);
	85	DecodeHWAddress(hwAddress, branch, FEC, chip, chan);
	86	bin = h2->FindBin(icol, irow);
	87
	88	// store the values
	89	fMeanPed[iSM][iRCU][branch][FEC][chip][chan] = h2->GetBinContent(bin);
	90	fRmsPed[iSM][iRCU][branch][FEC][chip][chan] = h2->GetBinError(bin);
75dbb105	91
	92	// report bad RMS channels:
	93	if (h2->GetBinError(bin) > kBadRMS) {
	94	printf(" bad pedestal RMS: iSM %d icol %d irow %d igain %d iRCU %d branch %d FEC %d chip %d chan %d - mean %4.1f rms %4.1f\n",
	95	iSM, icol, irow, igain, iRCU, branch, FEC, chip, chan,
	96	h2->GetBinContent(bin), h2->GetBinError(bin));
	97	}
29f94584	98	}
	99	}
	100
	101	return;
	102	}
	103
	104
	105	void
	106	GetPedValLEDRef(const Int_t iSM, const Int_t igain, const TProfile *h)
	107	{
	108	Int_t isect = iSM / 2; //
	109	Int_t iside = iSM % 2; // A or C side
	110	Int_t nStrips = h->GetNbinsX();
	111
75dbb105	112	if (kDebug) {
	113	printf("GetPedValLEDRef: iSM %d isect %d iside %d igain %d nStrips %d\n",
	114	iSM, isect, iside, igain, nStrips);
	115	}
29f94584	116	Int_t hwAddress = 0;
	117	Int_t iRCU = 0; // always true for LED Ref FEE
	118	Int_t branch = 0;
	119	Int_t FEC = 0;
	120	Int_t chip = 0;
	121	Int_t chan = 0;
	122
	123	Int_t icol = 0;
	124	Int_t irow = 0;
	125	Int_t bin = 0;
	126
	127	for (int istrip=0; istrip<nStrips; istrip++) {
	128
	129	hwAddress = GetHWAddressLEDRef(istrip, igain);
	130	DecodeHWAddress(hwAddress, branch, FEC, chip, chan);
	131	bin = h->FindBin(istrip);
	132
	133	// store the values
	134	fMeanPed[iSM][iRCU][branch][FEC][chip][chan] = h->GetBinContent(bin);
	135	fRmsPed[iSM][iRCU][branch][FEC][chip][chan] = h->GetBinError(bin);
	136	}
	137
	138	return;
	139	}
	140
	141	void
	142	PrintScript()
	143	{
	144	const char * sideStr[] = {"A","C"};
	145	const char * branchStr[] = {"A","B"};
	146	int VFPED = 0x06;
	147	int RCUWrite = 0x200000;
	148
	149	char filename[100];
	150	char scriptLine[200];
	151
	152	Int_t iSM = 0;
	153	Int_t iRCU = 0;
	154	Int_t ibranch = 0;
	155	Int_t iFEC = 0;
	156	Int_t ichip = 0;
	157	Int_t ichan = 0;
	158	Int_t Ped = 0;
	159
	160	for (iSM=0; iSM<kNSM; iSM++) {
	161	int iside = iSM % 2;
	162	int isect = iSM / 2;
	163	for (iRCU=0; iRCU<kNRCU; iRCU++) {
	164	sprintf(filename, "setSM%1s%dRCU%d.scr",
	165	sideStr[iside], isect, iRCU);
	166	ofstream fout(filename);
	167	int nscriptLines = 0;
	168
	169	for (ibranch=0; ibranch<kNBranch; ibranch++) {
	170	int firstFEC = 1;
	171	if (ibranch==0 && iRCU==0) { // extra LED Ref FEE in this location
	172	firstFEC = 0;
	173	}
	174
	175	for (iFEC=firstFEC; iFEC<kNFEC; iFEC++) { // FEC 1..9 (or 0..9)
	176	for (ichip=0; ichip<kNChip; ichip++) { // ALTRO 0,2,3,4
	177	if (ichip!=1) {
	178	for (ichan=0; ichan<kNChan; ichan++) {
	179
180	if (iFEC!=0 \|\| (ichan<8 \|\| ichan>11)) {
181
182	Ped = TMath::Nint(fMeanPed[iSM][iRCU][ibranch][iFEC][ichip][ichan]);
75dbb105	183	// raise Ped value to max for channels with exceptionally large RMS
	184	if (fRmsPed[iSM][iRCU][ibranch][iFEC][ichip][ichan] > kBadRMS) {
	185	printf(" bad pedestal RMS: iSM %d iRCU %d ibranch %d iFEC %d ichip %d ichan %d - raising from %d to 0x3ff\n",
	186	iSM, iRCU, ibranch, iFEC, ichip, ichan, Ped);
	187	Ped = 0x3ff;
	188	}
29f94584	189	//
	190	int writeAddr = (ibranch << 16) \| (iFEC << 12) \| (ichip << 9) \|
	191	(ichan << 5) \| VFPED \| RCUWrite;
	192	sprintf(scriptLine, "w 0x%04x 0x%06x # Branch %s, Card %d, Altro %d, Chan %d",
	193	nscriptLines, writeAddr, branchStr[ibranch],
	194	iFEC, ichip, ichan);
	195	fout << scriptLine << endl;
	196	nscriptLines++;
	197
	198	int writeVal = (Ped \| RCUWrite);
	199	sprintf(scriptLine, "w 0x%04x 0x%06x # Pedestal 0x%x = %d",
	200	nscriptLines, writeVal, Ped, Ped);
	201	fout << scriptLine << endl;
	202	nscriptLines++;
	203
	204	}
	205	}
	206	}
	207	}
	208	}
	209	}
	210	// ending, with execute and update step..
	211	sprintf(scriptLine, "w 0x%04x 0x380000 # End of the sequence",
	212	nscriptLines);
	213	fout << scriptLine << endl;
	214	nscriptLines++;
	215
	216	sprintf(scriptLine, "w 0x%04x 0x3F0000 # End of the instruction memory",
	217	nscriptLines);
	218	fout << scriptLine << endl;
	219	nscriptLines++;
	220
	221	fout << "wait 100 us" << endl;
	222	fout << "w 0x5304 0x0 \# execute and update registers" << endl;
	223	fout.close();
	224	} // iRCU
	225	}// iSM
	226
	227	return;
	228	}
	229
	230	void
	231	Clear()
	232	{
	233	for (Int_t iSM=0; iSM<kNSM; iSM++) {
	234	for (Int_t iRCU=0; iRCU<kNRCU; iRCU++) {
	235	for (Int_t ibranch=0; ibranch<kNBranch; ibranch++) {
	236	for (Int_t iFEC=0; iFEC<kNFEC; iFEC++) {
	237	for (Int_t ichip=0; ichip<kNChip; ichip++) {
	238	for (Int_t ichan=0; ichan<kNChan; ichan++) {
	239	fMeanPed[iSM][iRCU][ibranch][iFEC][ichip][ichan] = 0;
	240	fRmsPed[iSM][iRCU][ibranch][iFEC][ichip][ichan] = 0;
	241	}
	242	}
	243	}
	244	}
	245	}
	246	}
	247
	248	for (int istrip=0; istrip<kNStrips; istrip++) {
	249	fHWAddrLEDRef[istrip][0] = 0;
	250	fHWAddrLEDRef[istrip][1] = 0;
	251	}
	252
253	return;
254	}
255
256	void
257	DecodeHWAddress(Int_t hwAddr, Int_t & branch, Int_t & FEC, Int_t & chip, Int_t & chan)
258	{
259	chan = hwAddr & 0xf;
260	chip = (hwAddr >> 4) & 0x7;
261	FEC = (hwAddr >> 7) & 0xf;
262	branch = (hwAddr >> 11) & 0x1;
263	return;
264	}
265
266	Int_t
267	GetHWAddress(Int_t iside, Int_t icol, Int_t irow, Int_t igain, Int_t & iRCU)
268	{
269	iRCU = -111;
270
271	//RCU0
272	if (0<=irow&&irow<8) iRCU=0; // first cable row
273	else if (8<=irow&&irow<16 && 0<=icol&&icol<24) iRCU=0; // first half;
274	//second cable row
275	//RCU1
276	else if(8<=irow&&irow<16 && 24<=icol&&icol<48) iRCU=1; // second half;
277	//second cable row
278	else if(16<=irow&&irow<24) iRCU=1; // third cable row
279
280	// swap for odd=C side, to allow us to cable both sides the same
281	Int_t iRCUSide = iRCU;
282	if (iside == 1) {
283	iRCU = 1 - iRCU;
284	iRCUSide = iRCU + 2; // to make it map file index
285	}
286	Int_t hwAddress = fMapping[iRCUSide]->GetHWAddress(irow, icol, igain);
287
288	return hwAddress;
289	}
290
291	Int_t
292	GetHWAddressLEDRef(Int_t istrip, Int_t igain)
293	{
294	Int_t iRCU = 0; // for both sides; LED ref info is the same for both sides
295	Int_t caloflag = 3; // AliCaloRawStreamV3::kLEDMonData;
296
297	Int_t hwAddress = fHWAddrLEDRef[istrip][igain];
298
299	return hwAddress;
300	}
301
302	void
303	GetMapping()
304	{
305	TString sides[]={"A","C"};
306	// Read mapping files from $ALICE_ROOT/CALO/mapping/*.data
307	TString path = gSystem->Getenv("ALICE_ROOT");
308	path += "/EMCAL/mapping/RCU";
309	TString path2;
310	for(Int_t j = 0; j < 2; j++){ // sides
311	for(Int_t i = 0; i < 2; i++) { // RCU
312	path2 = path;
313	path2 += i;
314	path2 += sides[j];
315	path2 += ".data";
75dbb105	316	if (kDebug) { printf("Mapping file: %s\n",path2.Data()); }
29f94584	317	fMapping[j*2 + i] = new AliCaloAltroMapping(path2.Data());
	318	}
	319	}
	320	return;
	321	}
	322
	323	void
	324	CreateMappingLEDRef()
	325	{
	326	Int_t iRCU = 0; // for both sides; LED ref info is the same for both sides
	327	Int_t caloflag = 3; // AliCaloRawStreamV3::kLEDMonData;
	328
	329	Int_t maxAddr = 1 << 7; // LED Ref FEE is in FEC pos 0, i.e. addr space 0..127
	330
	331	int nLEDRefFEEChan = 0;
	332
	333	Int_t branch = 0;
	334	Int_t FEC = 0;
	335	Int_t chip = 0;
	336	Int_t chan = 0;
	337	for (int hwaddr = 0; hwaddr<maxAddr; hwaddr++) {
	338
	339	DecodeHWAddress(hwaddr, branch, FEC, chip, chan);
	340	if ( (chip!=1 && chip<kNChip) && // ALTROs 0,2,3,4
	341	(chan<8 \|\| chan>11) ) { // actual installed LED Ref FEE channels
	342
	343	int istrip = fMapping[iRCU]->GetPad(hwaddr);
	344	int igain = fMapping[iRCU]->GetPadRow(hwaddr);
	345	int iflag = fMapping[iRCU]->GetSector(hwaddr);
	346	if (iflag == caloflag) {
	347	fHWAddrLEDRef[istrip][igain] = hwaddr;
	348	nLEDRefFEEChan++;
	349	}
	350	}
	351	}
	352
75dbb105	353	if (kDebug) { cout << " nLEDRefFEEChan " << nLEDRefFEEChan << endl; }
29f94584	354	}