[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
// 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
      \r
      HLTDebug("HG channel (x = %d, z = %d) with amplitude: %f --> Digit with energy: %f \n", \r
	       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
bdaaafe2	23	\r
ef44ec64	24	\r
	25	// see header file for class documentation\r
	26	// or\r
	27	// refer to README to build package\r
	28	// or\r
	29	// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt\r
	30	\r
	31	#include "AliHLTCaloDigitMaker.h"\r
4f4b7ba4	32	#include "AliHLTCaloConstantsHandler.h"\r
ef44ec64	33	#include "AliHLTCaloMapper.h"\r
ef44ec64	34	#include "AliHLTCaloChannelDataStruct.h"\r
	35	#include "AliHLTCaloChannelDataHeaderStruct.h"\r
	36	#include "AliHLTCaloDigitDataStruct.h"\r
f137c3c5	37	#include "AliHLTCaloCoordinate.h"\r
ef44ec64	38	#include "AliHLTCaloSharedMemoryInterfacev2.h" // added by PTH\r
	39	//#include "AliPHOSEMCAGeometry.h"\r
	40	#include "TH2F.h"\r
ef44ec64	41	#include "AliHLTCaloConstants.h"\r
4655c95e	42	#include "AliHLTLogging.h"\r
ef44ec64	43	\r
	44	ClassImp(AliHLTCaloDigitMaker);\r
	45	\r
	46	//using namespace CaloHLTConst;\r
	47	\r
	48	AliHLTCaloDigitMaker::AliHLTCaloDigitMaker(TString det) :\r
4f4b7ba4	49	AliHLTCaloConstantsHandler(det),\r
4655c95e	50	AliHLTLogging(),\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
dc3d3428	170	\r
ef44ec64	171	void\r
	172	AliHLTCaloDigitMaker::SetGlobalLowGainFactor(Float_t factor)\r
	173	{\r
	174	//See header file for documentation\r
	175	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	176	{\r
	177	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	178	{\r
	179	fLowGainFactors[x][z] = factor;\r
	180	}\r
	181	}\r
	182	}\r
	183	\r
dc3d3428	184	\r
ef44ec64	185	void\r
	186	AliHLTCaloDigitMaker::SetBadChannelMask(TH2F* badChannelHGHist, TH2F* badChannelLGHist, Float_t qCut)\r
	187	{\r
	188	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	189	{\r
	190	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	191	{\r
	192	if(badChannelHGHist->GetBinContent(x, z) < qCut && badChannelHGHist->GetBinContent(x, z) > 0)\r
	193	{\r
c1e4a18c	194	fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = true;\r
ef44ec64	195	}\r
	196	else\r
	197	{\r
c1e4a18c	198	fBadChannelMask[x][z][fCaloConstants->GetHIGHGAIN()] = false;\r
ef44ec64	199	}\r
	200	if(badChannelLGHist->GetBinContent(x, z) < qCut && badChannelLGHist->GetBinContent(x, z) > 0)\r
	201	{\r
c1e4a18c	202	fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false;\r
ef44ec64	203	}\r
	204	else\r
	205	{\r
c1e4a18c	206	fBadChannelMask[x][z][fCaloConstants->GetLOWGAIN()] = false;\r
ef44ec64	207	}\r
	208	}\r
	209	}\r
	210	}\r
	211	\r
	212	void\r
	213	AliHLTCaloDigitMaker::Reset()\r
	214	{\r
	215	fDigitCount = 0;\r
	216	for(int x = 0; x < fCaloConstants->GetNXCOLUMNSMOD(); x++)\r
	217	{\r
	218	for(int z = 0; z < fCaloConstants->GetNZROWSMOD(); z++)\r
	219	{\r
	220	fChannelBook[x][z] = 0;\r
	221	}\r
	222	} \r
	223	\r
	224	}\r
	225	\r
	226	\r
f137c3c5	227	void AliHLTCaloDigitMaker::AddDigit(AliHLTCaloChannelDataStruct* channelData, AliHLTCaloCoordinate &coord)\r
ef44ec64	228	{\r
ef44ec64	229	\r
ad44d760	230	AliHLTCaloDigitDataStruct *tmpDigit = fDigitStructPtr + 1;\r
	231	\r
	232	if(fChannelBook[coord.fX][coord.fZ])\r
	233	{\r
	234	tmpDigit = fDigitStructPtr;\r
	235	fDigitStructPtr = fChannelBook[coord.fX][coord.fZ];\r
	236	fDigitCount--;\r
	237	// printf("Going to overwrite digit: x = %d, z = %d, gain = %d, energy = %f\n", fDigitStructPtr->fX, fDigitStructPtr->fZ, fDigitStructPtr->fGain, fDigitStructPtr->fEnergy);\r
	238	}\r
	239	\r
f137c3c5	240	fChannelBook[coord.fX][coord.fZ] = fDigitStructPtr;\r
ef44ec64	241	\r
7c80a370	242	\r
f137c3c5	243	fDigitStructPtr->fX = coord.fX;\r
f137c3c5	244	fDigitStructPtr->fZ = coord.fZ;\r
ad44d760	245	fDigitStructPtr->fGain = coord.fGain;\r
f137c3c5	246	fDigitStructPtr->fOverflow = false;\r
7c80a370	247	\r
7c80a370	248	fDigitStructPtr->fID = fDigitStructPtr->fZ * fCaloConstants->GetNXCOLUMNSMOD() + fDigitStructPtr->fX;\r
7173b37b	249	\r
f137c3c5	250	if(coord.fGain == fCaloConstants->GetHIGHGAIN() )\r
ef44ec64	251	{\r
f137c3c5	252	fDigitStructPtr->fEnergy = channelData->fEnergy*fHighGainFactors[coord.fX][coord.fZ];\r
f137c3c5	253	if(channelData->fEnergy >= fMaxEnergy)\r
ef44ec64	254	{\r
	255	fDigitStructPtr->fOverflow = true;\r
	256	}\r
bdaaafe2	257	\r
bdaaafe2	258	HLTDebug("HG channel (x = %d, z = %d) with amplitude: %f --> Digit with energy: %f \n", \r
bdaaafe2	259	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
bdaaafe2	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