[u/mrichter/AliRoot.git] / HLT / CALO / AliHLTCaloDigitMaker.cxx

// $Id$\r
\r
/**************************************************************************\r
 * This file is property of and copyright by the ALICE HLT Project        * \r
 * All rights reserved.                                                   *\r
 *                                                                        *\r
 * Primary Authors: Oystein Djuvsland                                     *\r
 *                                                                        *\r
 * Permission to use, copy, modify and distribute this software and its   *\r
 * documentation strictly for non-commercial purposes is hereby granted   *\r
 * without fee, provided that the above copyright notice appears in all   *\r
 * copies and that both the copyright notice and this permission notice   *\r
 * appear in the supporting documentation. The authors make no claims     *\r
 * about the suitability of this software for any purpose. It is          * \r
 * provided "as is" without express or implied warranty.                  *\r
 **************************************************************************/\r
 /** \r
 * @file   AliHLTCALOClusterizer.cxx\r
 * @author Oystein Djuvsland\r
 * @date \r
 * @brief  Digit maker for CALO HLT  \r
 */\r
  \r
\r
    \r
\r
// see header file for class documentation\r
// or\r
// refer to README to build package\r
// or\r
// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt\r
\r
#include "AliHLTCaloDigitMaker.h"\r
#include "AliHLTCaloConstantsHandler.h"\r
#include "AliHLTCaloMapper.h"\r
#include "AliHLTCaloChannelDataStruct.h"\r
#include "AliHLTCaloChannelDataHeaderStruct.h"\r
#include "AliHLTCaloDigitDataStruct.h"\r
#include "AliHLTCaloCoordinate.h"\r
#include "AliHLTCaloSharedMemoryInterfacev2.h" // added by PTH\r
//#include "AliPHOSEMCAGeometry.h"\r
#include "TH2F.h"\r
#include "AliHLTCaloConstants.h"\r
#include "AliHLTLogging.h"\r
\r
ClassImp(AliHLTCaloDigitMaker);\r
\r
//using namespace CaloHLTConst;\r
\r
AliHLTCaloDigitMaker::AliHLTCaloDigitMaker(TString det) :\r
  AliHLTCaloConstantsHandler(det),\r
  AliHLTLogging(),\r
  fShmPtr(0),\r
  fDigitStructPtr(0),\r
  fDigitCount(0),\r
  fMapperPtr(0),\r
  fHighGainFactors(0),\r
  fLowGainFactors(0),\r
  fBadChannelMask(0),\r
  fChannelBook(0),\r
  fMaxEnergy(900),\r
  fMinTime(0.0),\r
  fMaxTime(1008.0)\r
{\r
  // See header file for documentation\r
\r
  fShmPtr = new AliHLTCaloSharedMemoryInterfacev2(det);\r
\r
  fHighGainFactors = new Float_t*[fCaloConstants->GetNXCOLUMNSMOD()];\r
  fLowGainFactors = new Float_t*[fCaloConstants->GetNXCOLUMNSMOD()];\r
  \r
  fBadChannelMask = new Bool_t**[fCaloConstants->GetNXCOLUMNSMOD()];\r
  \r
  fChannelBook= new AliHLTCaloDigitDataStruct**[fCaloConstants->GetNXCOLUMNSMOD()];\r
  \r
  for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
    {\r
      fHighGainFactors[x] = new Float_t[fCaloConstants->GetNZROWSMOD()];\r
      fLowGainFactors[x] = new Float_t[fCaloConstants->GetNZROWSMOD()];\r
\r
      fBadChannelMask[x] = new Bool_t*[fCaloConstants->GetNZROWSMOD()];\r
\r
      fChannelBook[x] = new AliHLTCaloDigitDataStruct*[fCaloConstants->GetNZROWSMOD()];\r
\r
      for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	{\r
\r
	  fHighGainFactors[x][z] = 0.0153;\r
	  fLowGainFactors[x][z] = 0.245;\r
	 \r
	  fBadChannelMask[x][z] = new Bool_t[fCaloConstants->GetNGAINS()];\r
	  fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = false;\r
	  fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false; \r
	  \r
	  fChannelBook[x][z] = 0;\r
	  \r
	}\r
    }	  \r
}\r
\r
AliHLTCaloDigitMaker::~AliHLTCaloDigitMaker() \r
{\r
  //See header file for documentation\r
}\r
\r
Int_t\r
AliHLTCaloDigitMaker::MakeDigits(AliHLTCaloChannelDataHeaderStruct* channelDataHeader, AliHLTUInt32_t availableSize)\r
{\r
  //See header file for documentation\r
  \r
  Reset();\r
\r
  UInt_t totSize = sizeof(AliHLTCaloDigitDataStruct);\r
  \r
//   Int_t xMod = -1;\r
//   Int_t zMod = -1;\r
  \r
\r
  AliHLTCaloCoordinate coord;\r
  \r
  \r
  AliHLTCaloChannelDataStruct* currentchannel = 0;\r
  \r
  fShmPtr->SetMemory(channelDataHeader);\r
  currentchannel = fShmPtr->NextChannel();\r
\r
  while(currentchannel != 0)\r
    {\r
      if(availableSize < totSize) return -1;\r
\r
      fMapperPtr->ChannelId2Coordinate(currentchannel->fChannelID, coord);\r
      \r
      //      fMapperPtr->GetLocalCoord(currentchannel->fChannelID, locCoord);\r
      if(UseDigit(coord, currentchannel))\r
      {\r
	AddDigit(currentchannel, coord);\r
	//	j++;	      \r
	totSize += sizeof(AliHLTCaloDigitDataStruct);\r
      }\r
      currentchannel = fShmPtr->NextChannel(); // Get the next channel\r
    }\r
//       if(currentchannel)\r
// 	{\r
// 	  fMapperPtr->GetLocalCoord(currentchannel->fChannelID, locCoord);\r
// 	  if(UseDigit(coord1, currentchannel))\r
// 	    {\r
// 	      AddDigit(currentchannel, coord1, locCoord);\r
// 	      j++;	      \r
// 	      totSize += sizeof(AliHLTCaloDigitDataStruct);\r
// 	    }\r
// 	  currentchannel = fShmPtr->NextChannel(); // Get the next channel\r
// 	}\r
//     }\r
   \r
//   fDigitCount += j;\r
   return fDigitCount; \r
}\r
\r
void \r
AliHLTCaloDigitMaker::SetGlobalHighGainFactor(Float_t factor)\r
{\r
  //See header file for documentation\r
  for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
    {\r
      for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	{\r
	  fHighGainFactors[x][z] = factor;\r
	}\r
    }\r
}\r
\r
\r
void\r
AliHLTCaloDigitMaker::SetGlobalLowGainFactor(Float_t factor)\r
{\r
  //See header file for documentation\r
  for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
    {\r
      for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	{\r
	  fLowGainFactors[x][z] = factor;\r
	}\r
    }\r
}\r
\r
\r
void\r
AliHLTCaloDigitMaker::SetBadChannelMask(TH2F* badChannelHGHist, TH2F* badChannelLGHist, Float_t qCut)\r
{\r
 for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
    {\r
      for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	{\r
	  if(badChannelHGHist->GetBinContent(x, z) < qCut && badChannelHGHist->GetBinContent(x, z) > 0)\r
	    {\r
	      fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = true;\r
	    }\r
	  else\r
	    {\r
	      fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = false;\r
	    }\r
	  if(badChannelLGHist->GetBinContent(x, z) < qCut && badChannelLGHist->GetBinContent(x, z) > 0)\r
	    {\r
	      fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false;\r
	    }\r
	  else\r
	    {\r
	      fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false;\r
	    }\r
	}\r
    }\r
}\r
\r
void\r
AliHLTCaloDigitMaker::Reset()\r
{\r
  fDigitCount = 0;\r
  for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
    {\r
      for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	{\r
	  fChannelBook[x][z] = 0;\r
	}\r
    }	  \r
\r
}\r
\r
\r
void AliHLTCaloDigitMaker::AddDigit(AliHLTCaloChannelDataStruct* channelData, AliHLTCaloCoordinate &coord)\r
{\r
\r
  AliHLTCaloDigitDataStruct *tmpDigit = fDigitStructPtr + 1;\r
\r
  if(fChannelBook[coord.fX][coord.fZ])\r
    {\r
      tmpDigit = fDigitStructPtr;\r
      fDigitStructPtr = fChannelBook[coord.fX][coord.fZ];\r
      fDigitCount--;\r
      //      printf("Going to overwrite digit: x = %d, z = %d, gain = %d, energy = %f\n", fDigitStructPtr->fX, fDigitStructPtr->fZ, fDigitStructPtr->fGain, fDigitStructPtr->fEnergy);\r
    }\r
  \r
  fChannelBook[coord.fX][coord.fZ] = fDigitStructPtr;\r
\r
 \r
  fDigitStructPtr->fX = coord.fX;\r
  fDigitStructPtr->fZ = coord.fZ;\r
  fDigitStructPtr->fGain = coord.fGain;\r
  fDigitStructPtr->fOverflow = false;\r
  \r
  fDigitStructPtr->fID = fDigitStructPtr->fZ * fCaloConstants->GetNXCOLUMNSMOD() + fDigitStructPtr->fX;\r
  \r
  if(coord.fGain == fCaloConstants->GetHIGHGAIN() )\r
    {\r
      fDigitStructPtr->fEnergy = channelData->fEnergy*fHighGainFactors[coord.fX][coord.fZ];\r
      if(channelData->fEnergy >= fMaxEnergy)\r
	{\r
	  fDigitStructPtr->fOverflow = true;\r
	}\r
            HLTDebug("HG channel (x = %d, z = %d) with amplitude: %f --> Digit with energy: %f \n", coord.fX, coord.fZ, channelData->fEnergy, fDigitStructPtr->fEnergy);\r
    }\r
  else\r
    {\r
      fDigitStructPtr->fEnergy = channelData->fEnergy*fLowGainFactors[coord.fX][coord.fZ];\r
      if(channelData->fEnergy >= fMaxEnergy)\r
	{\r
	  fDigitStructPtr->fOverflow = true;\r
	}\r
	    HLTDebug("LG channel (x = %d, z = %d) with amplitude: %f --> Digit with energy: %f\n", coord.fX, coord.fZ, channelData->fEnergy, fDigitStructPtr->fEnergy); \r
    }\r
  fDigitStructPtr->fTime = channelData->fTime * 0.0000001; //TODO\r
  fDigitStructPtr->fCrazyness = channelData->fCrazyness;\r
  fDigitStructPtr->fModule = coord.fModuleId;\r
  fDigitStructPtr = tmpDigit;\r
  //  fDigitStructPtr++;\r
  fDigitCount++;\r
}\r
\r
bool AliHLTCaloDigitMaker::UseDigit(AliHLTCaloCoordinate &channelCoordinates, AliHLTCaloChannelDataStruct *channel) \r
{\r
   \r
   if(fBadChannelMask[channelCoordinates.fX][channelCoordinates.fZ][0] == true) return false;\r
   if(channel->fTime < fMinTime || channel->fTime > fMaxTime) return false;\r
   \r
  AliHLTCaloDigitDataStruct *tmpDigit = fChannelBook[channelCoordinates.fX][channelCoordinates.fZ];\r
  //printf("UseDigit: Got digit, x: %d, z: %d, gain: %d, amp: %f\n", channelCoordinates.fX, channelCoordinates.fZ, channelCoordinates.fGain, channel->fEnergy);\r
  if(tmpDigit)\r
    {\r
      if(channelCoordinates.fGain == fCaloConstants->GetLOWGAIN())\r
	{\r
	  	  //printf("UseDigit: Already have digit with, x: %d, z: %d, with high gain \n", channelCoordinates.fX, channelCoordinates.fZ);\r
	  if(tmpDigit->fOverflow)\r
	    {\r
	      	     // printf("But it was in overflow! Let's use this low gain!\n");\r
	      return true;\r
	    }\r
	  return false;\r
	}\r
      else\r
	{\r
	  	  //printf("UseDigit: Already have digit with, x: %d, z: %d, with low gain: %d\n", channelCoordinates.fX, channelCoordinates.fZ);\r
	  if(channel->fEnergy > fMaxEnergy )\r
	    {\r
	      return false;\r
	    }\r
	  return true;\r
	}\r
    }\r
  return true;\r
}\r
\r
void AliHLTCaloDigitMaker::SetBadChannel(Int_t x, Int_t z, Bool_t bad)\r
{\r
   // See header file for class documentation\r
   fBadChannelMask[x][z][0] = bad;\r
   fBadChannelMask[x][z][1] = bad;\r
}\r
\r
void AliHLTCaloDigitMaker::SetGain(Int_t x, Int_t z, Float_t ratio, Float_t gain)\r
{\r
   // See header file for class documentation\r
   \r
   fHighGainFactors[x][z] = gain;\r
   fLowGainFactors[x][z] = gain * ratio;\r
   \r
}\r
Commit	Line	Data
ef44ec64	1	// $Id$\r
	2	\r
	3	/**************************************************************************\r
	4	* This file is property of and copyright by the ALICE HLT Project * \r
	5	* All rights reserved. *\r
	6	* *\r
	7	* Primary Authors: Oystein Djuvsland *\r
	8	* *\r
	9	* Permission to use, copy, modify and distribute this software and its *\r
	10	* documentation strictly for non-commercial purposes is hereby granted *\r
	11	* without fee, provided that the above copyright notice appears in all *\r
	12	* copies and that both the copyright notice and this permission notice *\r
	13	* appear in the supporting documentation. The authors make no claims *\r
	14	* about the suitability of this software for any purpose. It is * \r
	15	* provided "as is" without express or implied warranty. *\r
	16	**************************************************************************/\r
	17	/** \r
	18	* @file AliHLTCALOClusterizer.cxx\r
	19	* @author Oystein Djuvsland\r
	20	* @date \r
	21	* @brief Digit maker for CALO HLT \r
	22	*/\r
	23	\r
	24	\r
	25	\r
	26	\r
	27	// see header file for class documentation\r
	28	// or\r
	29	// refer to README to build package\r
	30	// or\r
	31	// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt\r
	32	\r
	33	#include "AliHLTCaloDigitMaker.h"\r
4f4b7ba4	34	#include "AliHLTCaloConstantsHandler.h"\r
ef44ec64	35	#include "AliHLTCaloMapper.h"\r
ef44ec64	36	#include "AliHLTCaloChannelDataStruct.h"\r
	37	#include "AliHLTCaloChannelDataHeaderStruct.h"\r
	38	#include "AliHLTCaloDigitDataStruct.h"\r
f137c3c5	39	#include "AliHLTCaloCoordinate.h"\r
ef44ec64	40	#include "AliHLTCaloSharedMemoryInterfacev2.h" // added by PTH\r
	41	//#include "AliPHOSEMCAGeometry.h"\r
	42	#include "TH2F.h"\r
ef44ec64	43	#include "AliHLTCaloConstants.h"\r
4655c95e	44	#include "AliHLTLogging.h"\r
ef44ec64	45	\r
	46	ClassImp(AliHLTCaloDigitMaker);\r
	47	\r
	48	//using namespace CaloHLTConst;\r
	49	\r
	50	AliHLTCaloDigitMaker::AliHLTCaloDigitMaker(TString det) :\r
4f4b7ba4	51	AliHLTCaloConstantsHandler(det),\r
4655c95e	52	AliHLTLogging(),\r
ef44ec64	53	fShmPtr(0),\r
	54	fDigitStructPtr(0),\r
	55	fDigitCount(0),\r
ef44ec64	56	fMapperPtr(0),\r
	57	fHighGainFactors(0),\r
	58	fLowGainFactors(0),\r
	59	fBadChannelMask(0),\r
f137c3c5	60	fChannelBook(0),\r
7173b37b	61	fMaxEnergy(900),\r
	62	fMinTime(0.0),\r
	63	fMaxTime(1008.0)\r
ef44ec64	64	{\r
ef44ec64	65	// See header file for documentation\r
ef44ec64	66	\r
4f4b7ba4	67	fShmPtr = new AliHLTCaloSharedMemoryInterfacev2(det);\r
ef44ec64	68	\r
	69	fHighGainFactors = new Float_t*[fCaloConstants->GetNXCOLUMNSMOD()];\r
	70	fLowGainFactors = new Float_t*[fCaloConstants->GetNXCOLUMNSMOD()];\r
4f4b7ba4	71	\r
c1e4a18c	72	fBadChannelMask = new Bool_t**[fCaloConstants->GetNXCOLUMNSMOD()];\r
4f4b7ba4	73	\r
ef44ec64	74	fChannelBook= new AliHLTCaloDigitDataStruct**[fCaloConstants->GetNXCOLUMNSMOD()];\r
4f4b7ba4	75	\r
ef44ec64	76	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	77	{\r
	78	fHighGainFactors[x] = new Float_t[fCaloConstants->GetNZROWSMOD()];\r
	79	fLowGainFactors[x] = new Float_t[fCaloConstants->GetNZROWSMOD()];\r
	80	\r
c1e4a18c	81	fBadChannelMask[x] = new Bool_t*[fCaloConstants->GetNZROWSMOD()];\r
ef44ec64	82	\r
	83	fChannelBook[x] = new AliHLTCaloDigitDataStruct*[fCaloConstants->GetNZROWSMOD()];\r
	84	\r
	85	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	86	{\r
	87	\r
7173b37b	88	fHighGainFactors[x][z] = 0.0153;\r
7173b37b	89	fLowGainFactors[x][z] = 0.245;\r
ef44ec64	90	\r
c1e4a18c	91	fBadChannelMask[x][z] = new Bool_t[fCaloConstants->GetNGAINS()];\r
	92	fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = false;\r
	93	fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false; \r
ef44ec64	94	\r
	95	fChannelBook[x][z] = 0;\r
	96	\r
	97	}\r
	98	} \r
c375e15d	99	}\r
c375e15d	100	\r
ef44ec64	101	AliHLTCaloDigitMaker::~AliHLTCaloDigitMaker() \r
	102	{\r
	103	//See header file for documentation\r
	104	}\r
	105	\r
	106	Int_t\r
	107	AliHLTCaloDigitMaker::MakeDigits(AliHLTCaloChannelDataHeaderStruct* channelDataHeader, AliHLTUInt32_t availableSize)\r
	108	{\r
	109	//See header file for documentation\r
	110	\r
37616445	111	Reset();\r
ad44d760	112	\r
ef44ec64	113	UInt_t totSize = sizeof(AliHLTCaloDigitDataStruct);\r
	114	\r
	115	// Int_t xMod = -1;\r
	116	// Int_t zMod = -1;\r
	117	\r
f137c3c5	118	\r
f137c3c5	119	AliHLTCaloCoordinate coord;\r
ef44ec64	120	\r
	121	\r
	122	AliHLTCaloChannelDataStruct* currentchannel = 0;\r
ef44ec64	123	\r
	124	fShmPtr->SetMemory(channelDataHeader);\r
	125	currentchannel = fShmPtr->NextChannel();\r
	126	\r
	127	while(currentchannel != 0)\r
	128	{\r
	129	if(availableSize < totSize) return -1;\r
	130	\r
f137c3c5	131	fMapperPtr->ChannelId2Coordinate(currentchannel->fChannelID, coord);\r
ef44ec64	132	\r
f137c3c5	133	// fMapperPtr->GetLocalCoord(currentchannel->fChannelID, locCoord);\r
f137c3c5	134	if(UseDigit(coord, currentchannel))\r
37616445	135	{\r
f137c3c5	136	AddDigit(currentchannel, coord);\r
ad44d760	137	// j++; \r
ad44d760	138	totSize += sizeof(AliHLTCaloDigitDataStruct);\r
37616445	139	}\r
37616445	140	currentchannel = fShmPtr->NextChannel(); // Get the next channel\r
ef44ec64	141	}\r
f137c3c5	142	// if(currentchannel)\r
	143	// {\r
	144	// fMapperPtr->GetLocalCoord(currentchannel->fChannelID, locCoord);\r
	145	// if(UseDigit(coord1, currentchannel))\r
	146	// {\r
	147	// AddDigit(currentchannel, coord1, locCoord);\r
	148	// j++; \r
	149	// totSize += sizeof(AliHLTCaloDigitDataStruct);\r
	150	// }\r
	151	// currentchannel = fShmPtr->NextChannel(); // Get the next channel\r
	152	// }\r
	153	// }\r
37616445	154	\r
ad44d760	155	// fDigitCount += j;\r
37616445	156	return fDigitCount; \r
ef44ec64	157	}\r
	158	\r
	159	void \r
	160	AliHLTCaloDigitMaker::SetGlobalHighGainFactor(Float_t factor)\r
	161	{\r
	162	//See header file for documentation\r
	163	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	164	{\r
	165	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	166	{\r
	167	fHighGainFactors[x][z] = factor;\r
	168	}\r
	169	}\r
	170	}\r
	171	\r
dc3d3428	172	\r
ef44ec64	173	void\r
	174	AliHLTCaloDigitMaker::SetGlobalLowGainFactor(Float_t factor)\r
	175	{\r
	176	//See header file for documentation\r
	177	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	178	{\r
	179	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	180	{\r
	181	fLowGainFactors[x][z] = factor;\r
	182	}\r
	183	}\r
	184	}\r
	185	\r
dc3d3428	186	\r
ef44ec64	187	void\r
	188	AliHLTCaloDigitMaker::SetBadChannelMask(TH2F* badChannelHGHist, TH2F* badChannelLGHist, Float_t qCut)\r
	189	{\r
	190	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	191	{\r
	192	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	193	{\r
	194	if(badChannelHGHist->GetBinContent(x, z) < qCut && badChannelHGHist->GetBinContent(x, z) > 0)\r
	195	{\r
c1e4a18c	196	fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = true;\r
ef44ec64	197	}\r
	198	else\r
	199	{\r
c1e4a18c	200	fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = false;\r
ef44ec64	201	}\r
	202	if(badChannelLGHist->GetBinContent(x, z) < qCut && badChannelLGHist->GetBinContent(x, z) > 0)\r
	203	{\r
c1e4a18c	204	fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false;\r
ef44ec64	205	}\r
	206	else\r
	207	{\r
c1e4a18c	208	fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false;\r
ef44ec64	209	}\r
	210	}\r
	211	}\r
	212	}\r
	213	\r
	214	void\r
	215	AliHLTCaloDigitMaker::Reset()\r
	216	{\r
	217	fDigitCount = 0;\r
	218	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	219	{\r
	220	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	221	{\r
	222	fChannelBook[x][z] = 0;\r
	223	}\r
	224	} \r
	225	\r
	226	}\r
	227	\r
	228	\r
f137c3c5	229	void AliHLTCaloDigitMaker::AddDigit(AliHLTCaloChannelDataStruct* channelData, AliHLTCaloCoordinate &coord)\r
ef44ec64	230	{\r
ef44ec64	231	\r
ad44d760	232	AliHLTCaloDigitDataStruct *tmpDigit = fDigitStructPtr + 1;\r
	233	\r
	234	if(fChannelBook[coord.fX][coord.fZ])\r
	235	{\r
	236	tmpDigit = fDigitStructPtr;\r
	237	fDigitStructPtr = fChannelBook[coord.fX][coord.fZ];\r
	238	fDigitCount--;\r
	239	// printf("Going to overwrite digit: x = %d, z = %d, gain = %d, energy = %f\n", fDigitStructPtr->fX, fDigitStructPtr->fZ, fDigitStructPtr->fGain, fDigitStructPtr->fEnergy);\r
	240	}\r
	241	\r
f137c3c5	242	fChannelBook[coord.fX][coord.fZ] = fDigitStructPtr;\r
ef44ec64	243	\r
7c80a370	244	\r
f137c3c5	245	fDigitStructPtr->fX = coord.fX;\r
f137c3c5	246	fDigitStructPtr->fZ = coord.fZ;\r
ad44d760	247	fDigitStructPtr->fGain = coord.fGain;\r
f137c3c5	248	fDigitStructPtr->fOverflow = false;\r
7c80a370	249	\r
7c80a370	250	fDigitStructPtr->fID = fDigitStructPtr->fZ * fCaloConstants->GetNXCOLUMNSMOD() + fDigitStructPtr->fX;\r
7173b37b	251	\r
f137c3c5	252	if(coord.fGain == fCaloConstants->GetHIGHGAIN() )\r
ef44ec64	253	{\r
f137c3c5	254	fDigitStructPtr->fEnergy = channelData->fEnergy*fHighGainFactors[coord.fX][coord.fZ];\r
f137c3c5	255	if(channelData->fEnergy >= fMaxEnergy)\r
ef44ec64	256	{\r
	257	fDigitStructPtr->fOverflow = true;\r
	258	}\r
4655c95e	259	HLTDebug("HG channel (x = %d, z = %d) with amplitude: %f --> Digit with energy: %f \n", coord.fX, coord.fZ, channelData->fEnergy, fDigitStructPtr->fEnergy);\r
ef44ec64	260	}\r
	261	else\r
	262	{\r
f137c3c5	263	fDigitStructPtr->fEnergy = channelData->fEnergy*fLowGainFactors[coord.fX][coord.fZ];\r
ac4d2ac8	264	if(channelData->fEnergy >= fMaxEnergy)\r
ef44ec64	265	{\r
	266	fDigitStructPtr->fOverflow = true;\r
	267	}\r
4655c95e	268	HLTDebug("LG channel (x = %d, z = %d) with amplitude: %f --> Digit with energy: %f\n", coord.fX, coord.fZ, channelData->fEnergy, fDigitStructPtr->fEnergy); \r
ef44ec64	269	}\r
	270	fDigitStructPtr->fTime = channelData->fTime * 0.0000001; //TODO\r
	271	fDigitStructPtr->fCrazyness = channelData->fCrazyness;\r
f137c3c5	272	fDigitStructPtr->fModule = coord.fModuleId;\r
ad44d760	273	fDigitStructPtr = tmpDigit;\r
	274	// fDigitStructPtr++;\r
	275	fDigitCount++;\r
ef44ec64	276	}\r
ef44ec64	277	\r
f137c3c5	278	bool AliHLTCaloDigitMaker::UseDigit(AliHLTCaloCoordinate &channelCoordinates, AliHLTCaloChannelDataStruct *channel) \r
ef44ec64	279	{\r
7173b37b	280	\r
c1e4a18c	281	if(fBadChannelMask[channelCoordinates.fX][channelCoordinates.fZ][0] == true) return false;\r
7173b37b	282	if(channel->fTime < fMinTime \|\| channel->fTime > fMaxTime) return false;\r
7173b37b	283	\r
f137c3c5	284	AliHLTCaloDigitDataStruct *tmpDigit = fChannelBook[channelCoordinates.fX][channelCoordinates.fZ];\r
7c80a370	285	//printf("UseDigit: Got digit, x: %d, z: %d, gain: %d, amp: %f\n", channelCoordinates.fX, channelCoordinates.fZ, channelCoordinates.fGain, channel->fEnergy);\r
ef44ec64	286	if(tmpDigit)\r
ef44ec64	287	{\r
f137c3c5	288	if(channelCoordinates.fGain == fCaloConstants->GetLOWGAIN())\r
ef44ec64	289	{\r
cd9e2797	290	//printf("UseDigit: Already have digit with, x: %d, z: %d, with high gain \n", channelCoordinates.fX, channelCoordinates.fZ);\r
ef44ec64	291	if(tmpDigit->fOverflow)\r
ef44ec64	292	{\r
cd9e2797	293	// printf("But it was in overflow! Let's use this low gain!\n");\r
ef44ec64	294	return true;\r
	295	}\r
	296	return false;\r
	297	}\r
	298	else\r
	299	{\r
cd9e2797	300	//printf("UseDigit: Already have digit with, x: %d, z: %d, with low gain: %d\n", channelCoordinates.fX, channelCoordinates.fZ);\r
f137c3c5	301	if(channel->fEnergy > fMaxEnergy )\r
ef44ec64	302	{\r
	303	return false;\r
	304	}\r
	305	return true;\r
	306	}\r
	307	}\r
	308	return true;\r
	309	}\r
c1e4a18c	310	\r
	311	void AliHLTCaloDigitMaker::SetBadChannel(Int_t x, Int_t z, Bool_t bad)\r
	312	{\r
	313	// See header file for class documentation\r
c1e4a18c	314	fBadChannelMask[x][z][0] = bad;\r
	315	fBadChannelMask[x][z][1] = bad;\r
	316	}\r
	317	\r
	318	void AliHLTCaloDigitMaker::SetGain(Int_t x, Int_t z, Float_t ratio, Float_t gain)\r
	319	{\r
	320	// See header file for class documentation\r
	321	\r
	322	fHighGainFactors[x][z] = gain;\r
	323	fLowGainFactors[x][z] = gain * ratio;\r
	324	\r
	325	}\r