[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
\r
ClassImp(AliHLTCaloDigitMaker);\r
\r
//using namespace CaloHLTConst;\r
\r
AliHLTCaloDigitMaker::AliHLTCaloDigitMaker(TString det) :\r
  AliHLTCaloConstantsHandler(det),\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
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
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
            //printf("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
    //        printf("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
	44	\r
	45	ClassImp(AliHLTCaloDigitMaker);\r
	46	\r
	47	//using namespace CaloHLTConst;\r
	48	\r
	49	AliHLTCaloDigitMaker::AliHLTCaloDigitMaker(TString det) :\r
4f4b7ba4	50	AliHLTCaloConstantsHandler(det),\r
ef44ec64	51	fShmPtr(0),\r
	52	fDigitStructPtr(0),\r
	53	fDigitCount(0),\r
ef44ec64	54	fMapperPtr(0),\r
	55	fHighGainFactors(0),\r
	56	fLowGainFactors(0),\r
	57	fBadChannelMask(0),\r
f137c3c5	58	fChannelBook(0),\r
7173b37b	59	fMaxEnergy(900),\r
	60	fMinTime(0.0),\r
	61	fMaxTime(1008.0)\r
ef44ec64	62	{\r
ef44ec64	63	// See header file for documentation\r
ef44ec64	64	\r
4f4b7ba4	65	fShmPtr = new AliHLTCaloSharedMemoryInterfacev2(det);\r
ef44ec64	66	\r
	67	fHighGainFactors = new Float_t*[fCaloConstants->GetNXCOLUMNSMOD()];\r
	68	fLowGainFactors = new Float_t*[fCaloConstants->GetNXCOLUMNSMOD()];\r
4f4b7ba4	69	\r
c1e4a18c	70	fBadChannelMask = new Bool_t**[fCaloConstants->GetNXCOLUMNSMOD()];\r
4f4b7ba4	71	\r
ef44ec64	72	fChannelBook= new AliHLTCaloDigitDataStruct**[fCaloConstants->GetNXCOLUMNSMOD()];\r
4f4b7ba4	73	\r
ef44ec64	74	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	75	{\r
	76	fHighGainFactors[x] = new Float_t[fCaloConstants->GetNZROWSMOD()];\r
	77	fLowGainFactors[x] = new Float_t[fCaloConstants->GetNZROWSMOD()];\r
	78	\r
c1e4a18c	79	fBadChannelMask[x] = new Bool_t*[fCaloConstants->GetNZROWSMOD()];\r
ef44ec64	80	\r
	81	fChannelBook[x] = new AliHLTCaloDigitDataStruct*[fCaloConstants->GetNZROWSMOD()];\r
	82	\r
	83	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	84	{\r
	85	\r
7173b37b	86	fHighGainFactors[x][z] = 0.0153;\r
7173b37b	87	fLowGainFactors[x][z] = 0.245;\r
ef44ec64	88	\r
c1e4a18c	89	fBadChannelMask[x][z] = new Bool_t[fCaloConstants->GetNGAINS()];\r
	90	fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = false;\r
	91	fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false; \r
ef44ec64	92	\r
	93	fChannelBook[x][z] = 0;\r
	94	\r
	95	}\r
	96	} \r
c375e15d	97	}\r
c375e15d	98	\r
ef44ec64	99	AliHLTCaloDigitMaker::~AliHLTCaloDigitMaker() \r
	100	{\r
	101	//See header file for documentation\r
	102	}\r
	103	\r
	104	Int_t\r
	105	AliHLTCaloDigitMaker::MakeDigits(AliHLTCaloChannelDataHeaderStruct* channelDataHeader, AliHLTUInt32_t availableSize)\r
	106	{\r
	107	//See header file for documentation\r
	108	\r
37616445	109	Reset();\r
ad44d760	110	\r
ef44ec64	111	UInt_t totSize = sizeof(AliHLTCaloDigitDataStruct);\r
	112	\r
	113	// Int_t xMod = -1;\r
	114	// Int_t zMod = -1;\r
	115	\r
f137c3c5	116	\r
f137c3c5	117	AliHLTCaloCoordinate coord;\r
ef44ec64	118	\r
	119	\r
	120	AliHLTCaloChannelDataStruct* currentchannel = 0;\r
ef44ec64	121	\r
	122	fShmPtr->SetMemory(channelDataHeader);\r
	123	currentchannel = fShmPtr->NextChannel();\r
	124	\r
	125	while(currentchannel != 0)\r
	126	{\r
	127	if(availableSize < totSize) return -1;\r
	128	\r
f137c3c5	129	fMapperPtr->ChannelId2Coordinate(currentchannel->fChannelID, coord);\r
ef44ec64	130	\r
f137c3c5	131	// fMapperPtr->GetLocalCoord(currentchannel->fChannelID, locCoord);\r
f137c3c5	132	if(UseDigit(coord, currentchannel))\r
37616445	133	{\r
f137c3c5	134	AddDigit(currentchannel, coord);\r
ad44d760	135	// j++; \r
ad44d760	136	totSize += sizeof(AliHLTCaloDigitDataStruct);\r
37616445	137	}\r
37616445	138	currentchannel = fShmPtr->NextChannel(); // Get the next channel\r
ef44ec64	139	}\r
f137c3c5	140	// if(currentchannel)\r
	141	// {\r
	142	// fMapperPtr->GetLocalCoord(currentchannel->fChannelID, locCoord);\r
	143	// if(UseDigit(coord1, currentchannel))\r
	144	// {\r
	145	// AddDigit(currentchannel, coord1, locCoord);\r
	146	// j++; \r
	147	// totSize += sizeof(AliHLTCaloDigitDataStruct);\r
	148	// }\r
	149	// currentchannel = fShmPtr->NextChannel(); // Get the next channel\r
	150	// }\r
	151	// }\r
37616445	152	\r
ad44d760	153	// fDigitCount += j;\r
37616445	154	return fDigitCount; \r
ef44ec64	155	}\r
	156	\r
	157	void \r
	158	AliHLTCaloDigitMaker::SetGlobalHighGainFactor(Float_t factor)\r
	159	{\r
	160	//See header file for documentation\r
	161	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	162	{\r
	163	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	164	{\r
	165	fHighGainFactors[x][z] = factor;\r
	166	}\r
	167	}\r
	168	}\r
	169	\r
	170	void\r
	171	AliHLTCaloDigitMaker::SetGlobalLowGainFactor(Float_t factor)\r
	172	{\r
	173	//See header file for documentation\r
	174	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	175	{\r
	176	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	177	{\r
	178	fLowGainFactors[x][z] = factor;\r
	179	}\r
	180	}\r
	181	}\r
	182	\r
	183	void\r
	184	AliHLTCaloDigitMaker::SetBadChannelMask(TH2F* badChannelHGHist, TH2F* badChannelLGHist, Float_t qCut)\r
	185	{\r
	186	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	187	{\r
	188	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	189	{\r
	190	if(badChannelHGHist->GetBinContent(x, z) < qCut && badChannelHGHist->GetBinContent(x, z) > 0)\r
	191	{\r
c1e4a18c	192	fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = true;\r
ef44ec64	193	}\r
	194	else\r
	195	{\r
c1e4a18c	196	fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = false;\r
ef44ec64	197	}\r
	198	if(badChannelLGHist->GetBinContent(x, z) < qCut && badChannelLGHist->GetBinContent(x, z) > 0)\r
	199	{\r
c1e4a18c	200	fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false;\r
ef44ec64	201	}\r
	202	else\r
	203	{\r
c1e4a18c	204	fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false;\r
ef44ec64	205	}\r
	206	}\r
	207	}\r
	208	}\r
	209	\r
	210	void\r
	211	AliHLTCaloDigitMaker::Reset()\r
	212	{\r
	213	fDigitCount = 0;\r
	214	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	215	{\r
	216	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	217	{\r
	218	fChannelBook[x][z] = 0;\r
	219	}\r
	220	} \r
	221	\r
	222	}\r
	223	\r
	224	\r
f137c3c5	225	void AliHLTCaloDigitMaker::AddDigit(AliHLTCaloChannelDataStruct* channelData, AliHLTCaloCoordinate &coord)\r
ef44ec64	226	{\r
ef44ec64	227	\r
ad44d760	228	AliHLTCaloDigitDataStruct *tmpDigit = fDigitStructPtr + 1;\r
	229	\r
	230	if(fChannelBook[coord.fX][coord.fZ])\r
	231	{\r
	232	tmpDigit = fDigitStructPtr;\r
	233	fDigitStructPtr = fChannelBook[coord.fX][coord.fZ];\r
	234	fDigitCount--;\r
	235	// printf("Going to overwrite digit: x = %d, z = %d, gain = %d, energy = %f\n", fDigitStructPtr->fX, fDigitStructPtr->fZ, fDigitStructPtr->fGain, fDigitStructPtr->fEnergy);\r
	236	}\r
	237	\r
f137c3c5	238	fChannelBook[coord.fX][coord.fZ] = fDigitStructPtr;\r
ef44ec64	239	\r
7c80a370	240	\r
f137c3c5	241	fDigitStructPtr->fX = coord.fX;\r
f137c3c5	242	fDigitStructPtr->fZ = coord.fZ;\r
ad44d760	243	fDigitStructPtr->fGain = coord.fGain;\r
f137c3c5	244	fDigitStructPtr->fOverflow = false;\r
7c80a370	245	\r
7c80a370	246	fDigitStructPtr->fID = fDigitStructPtr->fZ * fCaloConstants->GetNXCOLUMNSMOD() + fDigitStructPtr->fX;\r
7173b37b	247	\r
f137c3c5	248	if(coord.fGain == fCaloConstants->GetHIGHGAIN() )\r
ef44ec64	249	{\r
f137c3c5	250	fDigitStructPtr->fEnergy = channelData->fEnergy*fHighGainFactors[coord.fX][coord.fZ];\r
f137c3c5	251	if(channelData->fEnergy >= fMaxEnergy)\r
ef44ec64	252	{\r
	253	fDigitStructPtr->fOverflow = true;\r
	254	}\r
cd9e2797	255	//printf("HG channel (x = %d, z = %d) with amplitude: %f --> Digit with energy: %f \n", coord.fX, coord.fZ, channelData->fEnergy, fDigitStructPtr->fEnergy);\r
ef44ec64	256	}\r
	257	else\r
	258	{\r
f137c3c5	259	fDigitStructPtr->fEnergy = channelData->fEnergy*fLowGainFactors[coord.fX][coord.fZ];\r
ac4d2ac8	260	if(channelData->fEnergy >= fMaxEnergy)\r
ef44ec64	261	{\r
	262	fDigitStructPtr->fOverflow = true;\r
	263	}\r
c1e4a18c	264	// printf("LG channel (x = %d, z = %d) with amplitude: %f --> Digit with energy: %f\n", coord.fX, coord.fZ, channelData->fEnergy, fDigitStructPtr->fEnergy); \r
ef44ec64	265	}\r
	266	fDigitStructPtr->fTime = channelData->fTime * 0.0000001; //TODO\r
	267	fDigitStructPtr->fCrazyness = channelData->fCrazyness;\r
f137c3c5	268	fDigitStructPtr->fModule = coord.fModuleId;\r
ad44d760	269	fDigitStructPtr = tmpDigit;\r
	270	// fDigitStructPtr++;\r
	271	fDigitCount++;\r
ef44ec64	272	}\r
ef44ec64	273	\r
f137c3c5	274	bool AliHLTCaloDigitMaker::UseDigit(AliHLTCaloCoordinate &channelCoordinates, AliHLTCaloChannelDataStruct *channel) \r
ef44ec64	275	{\r
7173b37b	276	\r
c1e4a18c	277	if(fBadChannelMask[channelCoordinates.fX][channelCoordinates.fZ][0] == true) return false;\r
7173b37b	278	if(channel->fTime < fMinTime \|\| channel->fTime > fMaxTime) return false;\r
7173b37b	279	\r
f137c3c5	280	AliHLTCaloDigitDataStruct *tmpDigit = fChannelBook[channelCoordinates.fX][channelCoordinates.fZ];\r
7c80a370	281	//printf("UseDigit: Got digit, x: %d, z: %d, gain: %d, amp: %f\n", channelCoordinates.fX, channelCoordinates.fZ, channelCoordinates.fGain, channel->fEnergy);\r
ef44ec64	282	if(tmpDigit)\r
ef44ec64	283	{\r
f137c3c5	284	if(channelCoordinates.fGain == fCaloConstants->GetLOWGAIN())\r
ef44ec64	285	{\r
cd9e2797	286	//printf("UseDigit: Already have digit with, x: %d, z: %d, with high gain \n", channelCoordinates.fX, channelCoordinates.fZ);\r
ef44ec64	287	if(tmpDigit->fOverflow)\r
ef44ec64	288	{\r
cd9e2797	289	// printf("But it was in overflow! Let's use this low gain!\n");\r
ef44ec64	290	return true;\r
	291	}\r
	292	return false;\r
	293	}\r
	294	else\r
	295	{\r
cd9e2797	296	//printf("UseDigit: Already have digit with, x: %d, z: %d, with low gain: %d\n", channelCoordinates.fX, channelCoordinates.fZ);\r
f137c3c5	297	if(channel->fEnergy > fMaxEnergy )\r
ef44ec64	298	{\r
	299	return false;\r
	300	}\r
	301	return true;\r
	302	}\r
	303	}\r
	304	return true;\r
	305	}\r
c1e4a18c	306	\r
	307	void AliHLTCaloDigitMaker::SetBadChannel(Int_t x, Int_t z, Bool_t bad)\r
	308	{\r
	309	// See header file for class documentation\r
c1e4a18c	310	fBadChannelMask[x][z][0] = bad;\r
	311	fBadChannelMask[x][z][1] = bad;\r
	312	}\r
	313	\r
	314	void AliHLTCaloDigitMaker::SetGain(Int_t x, Int_t z, Float_t ratio, Float_t gain)\r
	315	{\r
	316	// See header file for class documentation\r
	317	\r
	318	fHighGainFactors[x][z] = gain;\r
	319	fLowGainFactors[x][z] = gain * ratio;\r
	320	\r
	321	}\r