[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRawDigit.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
/*
 
 
 Author: R. GUERNANE LPSC Grenoble CNRS/IN2P3
*/

// --- ROOT system ---
#include <Riostream.h>
#include <TMath.h>

// --- AliRoot header files ---
#include "AliEMCALRawDigit.h"
#include "AliLog.h"

ClassImp(AliEMCALRawDigit)

//____________________________________________________________________________
AliEMCALRawDigit::AliEMCALRawDigit() : TObject(),
fId(-1),
fNSamples(0),
fSamples(0x0),
fAmplitude(0),
fTime(0)
{
	// default ctor 
}

//____________________________________________________________________________
AliEMCALRawDigit::AliEMCALRawDigit(Int_t id, Int_t timeSamples[], Int_t nSamples) : TObject(),
fId(id),
fNSamples(nSamples),
fSamples(0x0),
fAmplitude(0),
fTime(0)
{
	//
	fSamples = new Int_t[fNSamples];
	for (Int_t i = 0; i < fNSamples; i++) fSamples[i] = timeSamples[i];
}

//____________________________________________________________________________
AliEMCALRawDigit::~AliEMCALRawDigit() 
{
  //dtor, delete array of time samples
  if(fSamples) delete [] fSamples;
}

//____________________________________________________________________________
void AliEMCALRawDigit::Clear(Option_t *) 
{
  // clear, delete array of time samples
  if(fSamples) delete [] fSamples;
}


//____________________________________________________________________________
Bool_t AliEMCALRawDigit::GetTimeSample(const Int_t iSample, Int_t& timeBin, Int_t& amp) const
{
	//
	if (iSample > fNSamples || iSample < 0) return kFALSE;
	
	amp     =  fSamples[iSample] & 0xFFF;
	timeBin = (fSamples[iSample] >> 12) & 0xFF;

	return kTRUE;
}

//____________________________________________________________________________
void AliEMCALRawDigit::SetTimeSamples(const Int_t timeSamples[], const Int_t nSamples) 
{
	//
	if (fSamples) 
	{
		AliDebug(1,"Samples already filled: delete first!");
		fNSamples = 0;
		delete [] fSamples;
	}
	
	fNSamples = nSamples;
	fSamples = new Int_t[fNSamples];
	for (Int_t i = 0; i < fNSamples; i++) fSamples[i] = timeSamples[i];
}

//____________________________________________________________________________
Bool_t AliEMCALRawDigit::GetMaximum(Int_t& amplitude, Int_t& time) const
{
	//
	if (!fNSamples)
	{
		AliDebug(1,"Digit has no time sample");
		return kFALSE;
	}
		
	amplitude = 0;
	for (Int_t i = 0; i < fNSamples; i++)
	{
		Int_t t, a;
		if (GetTimeSample(i, t, a))
		{
			if (a > amplitude)
			{
				amplitude = a;
				time      = t;
			}
		}
	}
	
	return kTRUE;
}

//____________________________________________________________________________
Int_t AliEMCALRawDigit::Compare(const TObject * obj) const
{
	// Compares two digits with respect to its Id
	// to sort according increasing Id
	
	Int_t rv=0;
	
	AliEMCALRawDigit* digit = (AliEMCALRawDigit*)obj; 
	
	Int_t iddiff = fId - digit->GetId();
	
	if (iddiff > 0) 
		rv =  1;
	else if (iddiff < 0)
		rv = -1; 
	else
		rv =  0;
	
	return rv; 
}

//____________________________________________________________________________
void AliEMCALRawDigit::Print(const Option_t* /*opt*/) const
{
	//
	printf("===\n| Digit id: %4d / %d Time Samples: \n",fId,fNSamples);
	for (Int_t i=0; i < fNSamples; i++) 
	{
		Int_t timeBin=-1, amp=0;
		GetTimeSample(i, timeBin, amp);
		printf("| (%d,%d) ",timeBin,amp);
	}
	
	printf("\n");
}
Commit	Line	Data
916f1e76	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
de39a0ff	15	/*
	16
	17
	18
	19
	20	Author: R. GUERNANE LPSC Grenoble CNRS/IN2P3
	21	*/
916f1e76	22
	23	// --- ROOT system ---
	24	#include <Riostream.h>
	25	#include <TMath.h>
	26
916f1e76	27	// --- AliRoot header files ---
916f1e76	28	#include "AliEMCALRawDigit.h"
de39a0ff	29	#include "AliLog.h"
916f1e76	30
	31	ClassImp(AliEMCALRawDigit)
	32
	33	//____________________________________________________________________________
	34	AliEMCALRawDigit::AliEMCALRawDigit() : TObject(),
	35	fId(-1),
	36	fNSamples(0),
de39a0ff	37	fSamples(0x0),
	38	fAmplitude(0),
	39	fTime(0)
916f1e76	40	{
	41	// default ctor
	42	}
	43
	44	//____________________________________________________________________________
	45	AliEMCALRawDigit::AliEMCALRawDigit(Int_t id, Int_t timeSamples[], Int_t nSamples) : TObject(),
	46	fId(id),
	47	fNSamples(nSamples),
de39a0ff	48	fSamples(0x0),
	49	fAmplitude(0),
	50	fTime(0)
916f1e76	51	{
916f1e76	52	//
de39a0ff	53	fSamples = new Int_t[fNSamples];
de39a0ff	54	for (Int_t i = 0; i < fNSamples; i++) fSamples[i] = timeSamples[i];
916f1e76	55	}
	56
	57	//____________________________________________________________________________
	58	AliEMCALRawDigit::~AliEMCALRawDigit()
	59	{
de39a0ff	60	//dtor, delete array of time samples
b57cc9b9	61	if(fSamples) delete [] fSamples;
916f1e76	62	}
916f1e76	63
e0dc3f7d	64	//____________________________________________________________________________
	65	void AliEMCALRawDigit::Clear(Option_t *)
	66	{
de39a0ff	67	// clear, delete array of time samples
b57cc9b9	68	if(fSamples) delete [] fSamples;
e0dc3f7d	69	}
	70
	71
916f1e76	72	//____________________________________________________________________________
	73	Bool_t AliEMCALRawDigit::GetTimeSample(const Int_t iSample, Int_t& timeBin, Int_t& amp) const
	74	{
de39a0ff	75	//
916f1e76	76	if (iSample > fNSamples \|\| iSample < 0) return kFALSE;
	77
	78	amp = fSamples[iSample] & 0xFFF;
	79	timeBin = (fSamples[iSample] >> 12) & 0xFF;
	80
	81	return kTRUE;
	82	}
	83
	84	//____________________________________________________________________________
de39a0ff	85	void AliEMCALRawDigit::SetTimeSamples(const Int_t timeSamples[], const Int_t nSamples)
916f1e76	86	{
de39a0ff	87	//
de39a0ff	88	if (fSamples)
916f1e76	89	{
aa878a7b	90	AliDebug(1,"Samples already filled: delete first!");
de39a0ff	91	fNSamples = 0;
de39a0ff	92	delete [] fSamples;
916f1e76	93	}
916f1e76	94
de39a0ff	95	fNSamples = nSamples;
	96	fSamples = new Int_t[fNSamples];
	97	for (Int_t i = 0; i < fNSamples; i++) fSamples[i] = timeSamples[i];
916f1e76	98	}
	99
	100	//____________________________________________________________________________
de39a0ff	101	Bool_t AliEMCALRawDigit::GetMaximum(Int_t& amplitude, Int_t& time) const
	102	{
	103	//
	104	if (!fNSamples)
	105	{
aa878a7b	106	AliDebug(1,"Digit has no time sample");
de39a0ff	107	return kFALSE;
	108	}
	109
	110	amplitude = 0;
	111	for (Int_t i = 0; i < fNSamples; i++)
	112	{
	113	Int_t t, a;
	114	if (GetTimeSample(i, t, a))
	115	{
	116	if (a > amplitude)
	117	{
	118	amplitude = a;
	119	time = t;
	120	}
	121	}
	122	}
	123
	124	return kTRUE;
	125	}
	126
	127	//____________________________________________________________________________
	128	Int_t AliEMCALRawDigit::Compare(const TObject * obj) const
916f1e76	129	{
	130	// Compares two digits with respect to its Id
	131	// to sort according increasing Id
de39a0ff	132
	133	Int_t rv=0;
	134
	135	AliEMCALRawDigit* digit = (AliEMCALRawDigit*)obj;
	136
	137	Int_t iddiff = fId - digit->GetId();
	138
	139	if (iddiff > 0)
	140	rv = 1;
	141	else if (iddiff < 0)
916f1e76	142	rv = -1;
916f1e76	143	else
de39a0ff	144	rv = 0;
de39a0ff	145
916f1e76	146	return rv;
	147	}
	148
	149	//____________________________________________________________________________
de39a0ff	150	void AliEMCALRawDigit::Print(const Option_t* /opt/) const
916f1e76	151	{
de39a0ff	152	//
	153	printf("===\n\| Digit id: %4d / %d Time Samples: \n",fId,fNSamples);
	154	for (Int_t i=0; i < fNSamples; i++)
	155	{
	156	Int_t timeBin=-1, amp=0;
	157	GetTimeSample(i, timeBin, amp);
	158	printf("\| (%d,%d) ",timeBin,amp);
	159	}
	160
	161	printf("\n");
916f1e76	162	}