[u/mrichter/AliRoot.git] / EMCAL / AliEMCALTriggerElectronics.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.                  *
 **************************************************************************/

/*
 
 
EMCal trigger electronics manager L0/L1
can handle both simulated digits and raw data
Author: R. GUERNANE LPSC Grenoble CNRS/IN2P3
*/

#include "AliEMCALTriggerElectronics.h"
#include "AliEMCALTriggerTRU.h"
#include "AliEMCALTriggerSTU.h"
#include "AliEMCALGeometry.h"
#include "AliRunLoader.h"
#include "AliEMCAL.h" 
#include "AliRun.h" 
#include "AliEMCALTriggerDCSConfig.h"
#include "AliEMCALTriggerData.h"
#include "AliEMCALDigit.h"
#include "AliCaloRawStreamV3.h"
#include "AliEMCALTriggerSTURawStream.h"
#include "AliEMCALDigit.h"
#include "AliEMCALTriggerRawDigit.h"
#include "AliEMCALTriggerPatch.h"

#include <TVector2.h>
#include <TClonesArray.h>

namespace
{
	const Int_t kNTRU = 30;
}

ClassImp(AliEMCALTriggerElectronics)

//__________________
AliEMCALTriggerElectronics::AliEMCALTriggerElectronics(const AliEMCALTriggerDCSConfig *dcsConf) : TObject(),
fTRU(new TClonesArray("AliEMCALTriggerTRU",32)),
fSTU(0x0)
{
	TVector2 rSize;
	
	rSize.Set( 24.,  4. );

	// 32 TRUs
	for (Int_t i=0;i<kNTRU;i++) 
	{
		AliEMCALTriggerTRUDCSConfig* truConf = dcsConf->GetTRUDCSConfig(i);
		new ((*fTRU)[i]) AliEMCALTriggerTRU(truConf, rSize, i % 2);
	}
	
	rSize.Set( 48., 64. );
	
	// 1 STU
	AliEMCALTriggerSTUDCSConfig* stuConf = dcsConf->GetSTUDCSConfig();
	fSTU = new AliEMCALTriggerSTU(stuConf, rSize);
	
	TString str = "map";
	for (Int_t i=0;i<kNTRU;i++) fSTU->Build(str,
											i,
											(static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->Map(),
											(static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->RegionSize() 
										    );
}

//________________
AliEMCALTriggerElectronics::~AliEMCALTriggerElectronics()
{
	//
	fTRU->Delete();
 	delete fSTU;
}

//__________________
void AliEMCALTriggerElectronics::Digits2Trigger(TClonesArray* digits, const Int_t V0M[], AliEMCALTriggerData* data)
{
	//
	AliEMCALGeometry* geom = 0x0;
	
	AliRunLoader *rl = AliRunLoader::Instance();
	if (rl->GetAliRun() && rl->GetAliRun()->GetDetector("EMCAL")){
	  AliEMCAL* emcal = dynamic_cast<AliEMCAL*>(rl->GetAliRun()->GetDetector("EMCAL"));
	  if(emcal)geom = emcal->GetGeometry();
	}
	
	if(!geom) geom =  AliEMCALGeometry::GetInstance(AliEMCALGeometry::GetDefaultGeometryName());
	
	if(!geom) AliError("Cannot access geometry!");
	
	//	digits->Sort();
	
	Int_t region[48][64], posMap[48][64];
	for (Int_t i = 0; i < 48; i++) for (Int_t j = 0; j < 64; j++) 
	{
		region[i][j] =  0;
		posMap[i][j] = -1;
	}
	
	for (Int_t i = 0; i < digits->GetEntriesFast(); i++)
	{
		AliEMCALTriggerRawDigit* digit = (AliEMCALTriggerRawDigit*)digits->At(i);
		
		Int_t id = digit->GetId();
		
		Int_t iTRU, iADC;
		
		Bool_t isOK1 = geom->GetTRUFromAbsFastORIndex(id, iTRU, iADC);
		
		if ((isOK1 && iTRU >= kNTRU) || !isOK1) continue;

		for (Int_t j = 0; j < digit->GetNSamples(); j++)
		{
			Int_t time, amp;
			Bool_t isOK2 = digit->GetTimeSample(j, time, amp);
			
			if (isOK1 && isOK2 && amp) (static_cast<AliEMCALTriggerTRU*>(fTRU->At(iTRU)))->SetADC(iADC, time, amp);
		}
		
		Int_t px, py;
		if (geom->GetPositionInEMCALFromAbsFastORIndex(id, px, py))
		{
			posMap[px][py] = i;
			
			if (fSTU->GetRawData() && digit->GetL1TimeSum() >= 0) 
			{
				region[px][py] = digit->GetL1TimeSum();
			}
		}
	}

	Int_t iL0 = 0;

	for (Int_t i=0; i<kNTRU; i++) 
	{
		AliDebug(999, Form("===========< TRU %2d >============\n", i));
		
		AliEMCALTriggerTRU* iTRU = static_cast<AliEMCALTriggerTRU*>(fTRU->At(i));

		// L0 is always computed from F-ALTRO
		if (iTRU->L0()) 
		{
			iL0 += iTRU->L0();
			
			Int_t sizeX = (Int_t) ((iTRU->PatchSize())->X() * (iTRU->SubRegionSize())->X());
			
			Int_t sizeY = (Int_t) ((iTRU->PatchSize())->Y() * (iTRU->SubRegionSize())->Y());
			
			// transform local to global 
			TIter Next(&iTRU->Patches());
			while (AliEMCALTriggerPatch* p = (AliEMCALTriggerPatch*)Next())
			{
				Int_t px, py, id; p->Position(px, py);
				
				if (geom->GetAbsFastORIndexFromPositionInTRU(i, px, py, id) 
					&& 
					geom->GetPositionInEMCALFromAbsFastORIndex(id, px, py)) p->SetPosition(px, py);
				
				if (!data->GetMode()) // Simulation
				{
					Int_t peaks = p->Peaks();
					
					Int_t pos;
					AliEMCALTriggerRawDigit* dig = 0x0;
					
					for (Int_t j = 0; j < sizeX * sizeY; j++)
					{
						if (peaks & (1 << j))
						{
							pos = posMap[px + j % sizeX][py + j / sizeX];
							
							if (pos == -1)
							{
								// Add a new digit
								new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
								
								dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
							}
							else
							{
								dig = (AliEMCALTriggerRawDigit*)digits->At(pos);
							}
							
							dig->SetL0Time(p->Time());
						}
					}
					
					pos = posMap[px][py];
					
					if (pos == -1)
					{
						// Add a new digit
						new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
						
						dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
					}
					else
					{
						dig = (AliEMCALTriggerRawDigit*)digits->At(pos);
					}
					
					dig->SetTriggerBit(kL0,0);
				}
			}
			
			data->SetL0Trigger(0, i, 1);
		}
		else
			data->SetL0Trigger(0, i, 0);
	}

	// A L0 has been issued, run L1
	// Depending on raw data enabled or not in STU data: L1 computation 
	// should be done from F-ALTRO or directly on TRU time sums in STU raw data
	if (iL0) 
	{
		// Use L1 threshold from raw data when reconstructing raw data
		if (data->GetMode())
		{
			fSTU->SetThreshold(kL1Gamma, data->GetL1GammaThreshold());
			fSTU->SetThreshold(kL1Jet,   data->GetL1JetThreshold()  );			
		}
		else
		{
			fSTU->ComputeThFromV0(V0M); // C/A
			data->SetL1GammaThreshold(fSTU->GetThreshold(kL1Gamma));
			data->SetL1JetThreshold(  fSTU->GetThreshold(kL1Jet)  );
		}
		
		if (fSTU->GetRawData())
		{
			// Compute L1 from STU raw data
			fSTU->SetRegion(region);
		}
		else
		{
			// Build STU raw data from F-ALTRO
			TString str = "region";
			for (Int_t i = 0; i < kNTRU; i++) fSTU->Build(str,
														  i, 
														  (static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->Region(), 
														  (static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->RegionSize());
		}

		fSTU->L1(kL1Gamma);
		
		Int_t id, px, py;
		AliEMCALTriggerRawDigit* dig = 0x0;
		
		TIterator* nP = 0x0;
		
		nP = (fSTU->Patches()).MakeIterator();
		
		while (AliEMCALTriggerPatch* p = (AliEMCALTriggerPatch*)nP->Next()) 
		{			
			p->Position(px, py);
			
			if (geom->GetAbsFastORIndexFromPositionInEMCAL(px, py, id))
			{
				if (posMap[px][py] == -1)
				{
					// Add a new digit
					new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
					
					dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
				} 
				else
				{
					dig = (AliEMCALTriggerRawDigit*)digits->At(posMap[px][py]);								
				}
				
				dig->SetTriggerBit(kL1Gamma,0);
			}
		}

		fSTU->Reset();

		fSTU->L1(kL1Jet);
		
		nP = (fSTU->Patches()).MakeIterator();
		
		while (AliEMCALTriggerPatch* p = (AliEMCALTriggerPatch*)nP->Next()) 
		{			
			p->Position(px, py);

			px *= (Int_t)((fSTU->SubRegionSize())->X());

			py *= (Int_t)((fSTU->SubRegionSize())->Y());
			
			if (geom->GetAbsFastORIndexFromPositionInEMCAL(px, py, id))
			{
				if (posMap[px][py] == -1)
				{
					// Add a new digit
					new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
					
					dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
				} 
				else
				{
					dig = (AliEMCALTriggerRawDigit*)digits->At(posMap[px][py]);								
				}
				
				dig->SetTriggerBit(kL1Jet,0);
			}
		}
		
		Int_t** reg = fSTU->Region();
		
		if (!fSTU->GetRawData())
		{
			// Update digits w/ L1 time sum
			// Done in raw digit maker when raw data enabled
			for (Int_t i = 0; i < 48; i++)
			{
				for (Int_t j = 0; j < 64; j++)
				{
					if (reg[i][j])
					{
						if (geom->GetAbsFastORIndexFromPositionInEMCAL(i, j, id))
						{
							if (posMap[i][j] == -1)
							{
								// Add a new digit with L1 time sum
								new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
								
								dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
							} 
							else
							{
								dig = (AliEMCALTriggerRawDigit*)digits->At(posMap[i][j]);								
							}
							
							dig->SetL1TimeSum(reg[i][j]);
						}
					}
				}
			}
		}
	}

	if (AliDebugLevel() >= 999) data->Scan();
	
	// Now reset the electronics for a fresh start with next event
	Reset();
}

//__________________
void AliEMCALTriggerElectronics::Reset()
{
	//
	TIter NextTRU(fTRU);
	while ( AliEMCALTriggerTRU *TRU = (AliEMCALTriggerTRU*)NextTRU() ) TRU->Reset();
	
	fSTU->Reset();
}
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	**************************************************************************/
	15
	16	/*
	17
	18
	19	EMCal trigger electronics manager L0/L1
	20	can handle both simulated digits and raw data
	21	Author: R. GUERNANE LPSC Grenoble CNRS/IN2P3
	22	*/
	23
	24	#include "AliEMCALTriggerElectronics.h"
	25	#include "AliEMCALTriggerTRU.h"
	26	#include "AliEMCALTriggerSTU.h"
	27	#include "AliEMCALGeometry.h"
	28	#include "AliRunLoader.h"
	29	#include "AliEMCAL.h"
	30	#include "AliRun.h"
fff39dd1	31	#include "AliEMCALTriggerDCSConfig.h"
916f1e76	32	#include "AliEMCALTriggerData.h"
	33	#include "AliEMCALDigit.h"
	34	#include "AliCaloRawStreamV3.h"
916f1e76	35	#include "AliEMCALTriggerSTURawStream.h"
916f1e76	36	#include "AliEMCALDigit.h"
de39a0ff	37	#include "AliEMCALTriggerRawDigit.h"
de39a0ff	38	#include "AliEMCALTriggerPatch.h"
916f1e76	39
	40	#include <TVector2.h>
	41	#include <TClonesArray.h>
916f1e76	42
	43	namespace
	44	{
8cf986c8	45	const Int_t kNTRU = 30;
916f1e76	46	}
	47
	48	ClassImp(AliEMCALTriggerElectronics)
	49
	50	//__________________
fff39dd1	51	AliEMCALTriggerElectronics::AliEMCALTriggerElectronics(const AliEMCALTriggerDCSConfig *dcsConf) : TObject(),
916f1e76	52	fTRU(new TClonesArray("AliEMCALTriggerTRU",32)),
	53	fSTU(0x0)
	54	{
	55	TVector2 rSize;
	56
	57	rSize.Set( 24., 4. );
	58
	59	// 32 TRUs
fff39dd1	60	for (Int_t i=0;i<kNTRU;i++)
	61	{
	62	AliEMCALTriggerTRUDCSConfig* truConf = dcsConf->GetTRUDCSConfig(i);
de39a0ff	63	new ((*fTRU)[i]) AliEMCALTriggerTRU(truConf, rSize, i % 2);
fff39dd1	64	}
fff39dd1	65
916f1e76	66	rSize.Set( 48., 64. );
	67
	68	// 1 STU
fff39dd1	69	AliEMCALTriggerSTUDCSConfig* stuConf = dcsConf->GetSTUDCSConfig();
fff39dd1	70	fSTU = new AliEMCALTriggerSTU(stuConf, rSize);
916f1e76	71
de39a0ff	72	TString str = "map";
	73	for (Int_t i=0;i<kNTRU;i++) fSTU->Build(str,
	74	i,
	75	(static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->Map(),
	76	(static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->RegionSize()
	77	);
916f1e76	78	}
	79
	80	//________________
	81	AliEMCALTriggerElectronics::~AliEMCALTriggerElectronics()
	82	{
	83	//
	84	fTRU->Delete();
	85	delete fSTU;
	86	}
	87
	88	//__________________
de39a0ff	89	void AliEMCALTriggerElectronics::Digits2Trigger(TClonesArray* digits, const Int_t V0M[], AliEMCALTriggerData* data)
916f1e76	90	{
	91	//
	92	AliEMCALGeometry* geom = 0x0;
	93
	94	AliRunLoader *rl = AliRunLoader::Instance();
a51e676d	95	if (rl->GetAliRun() && rl->GetAliRun()->GetDetector("EMCAL")){
	96	AliEMCAL* emcal = dynamic_cast<AliEMCAL*>(rl->GetAliRun()->GetDetector("EMCAL"));
	97	if(emcal)geom = emcal->GetGeometry();
	98	}
	99
	100	if(!geom) geom = AliEMCALGeometry::GetInstance(AliEMCALGeometry::GetDefaultGeometryName());
	101
	102	if(!geom) AliError("Cannot access geometry!");
916f1e76	103
de39a0ff	104	// digits->Sort();
	105
	106	Int_t region[48][64], posMap[48][64];
	107	for (Int_t i = 0; i < 48; i++) for (Int_t j = 0; j < 64; j++)
	108	{
	109	region[i][j] = 0;
	110	posMap[i][j] = -1;
	111	}
	112
	113	for (Int_t i = 0; i < digits->GetEntriesFast(); i++)
916f1e76	114	{
de39a0ff	115	AliEMCALTriggerRawDigit* digit = (AliEMCALTriggerRawDigit*)digits->At(i);
	116
	117	Int_t id = digit->GetId();
	118
	119	Int_t iTRU, iADC;
	120
	121	Bool_t isOK1 = geom->GetTRUFromAbsFastORIndex(id, iTRU, iADC);
	122
c33e980f	123	if ((isOK1 && iTRU >= kNTRU) \|\| !isOK1) continue;
c33e980f	124
de39a0ff	125	for (Int_t j = 0; j < digit->GetNSamples(); j++)
916f1e76	126	{
de39a0ff	127	Int_t time, amp;
de39a0ff	128	Bool_t isOK2 = digit->GetTimeSample(j, time, amp);
916f1e76	129
de39a0ff	130	if (isOK1 && isOK2 && amp) (static_cast<AliEMCALTriggerTRU*>(fTRU->At(iTRU)))->SetADC(iADC, time, amp);
	131	}
	132
	133	Int_t px, py;
	134	if (geom->GetPositionInEMCALFromAbsFastORIndex(id, px, py))
	135	{
	136	posMap[px][py] = i;
916f1e76	137
de39a0ff	138	if (fSTU->GetRawData() && digit->GetL1TimeSum() >= 0)
916f1e76	139	{
de39a0ff	140	region[px][py] = digit->GetL1TimeSum();
916f1e76	141	}
	142	}
	143	}
916f1e76	144
de39a0ff	145	Int_t iL0 = 0;
916f1e76	146
	147	for (Int_t i=0; i<kNTRU; i++)
	148	{
de39a0ff	149	AliDebug(999, Form("===========< TRU %2d >============\n", i));
916f1e76	150
de39a0ff	151	AliEMCALTriggerTRU* iTRU = static_cast<AliEMCALTriggerTRU*>(fTRU->At(i));
916f1e76	152
de39a0ff	153	// L0 is always computed from F-ALTRO
de39a0ff	154	if (iTRU->L0())
916f1e76	155	{
de39a0ff	156	iL0 += iTRU->L0();
de39a0ff	157
ddd131cc	158	Int_t sizeX = (Int_t) ((iTRU->PatchSize())->X() * (iTRU->SubRegionSize())->X());
de39a0ff	159
ddd131cc	160	Int_t sizeY = (Int_t) ((iTRU->PatchSize())->Y() * (iTRU->SubRegionSize())->Y());
de39a0ff	161
	162	// transform local to global
	163	TIter Next(&iTRU->Patches());
	164	while (AliEMCALTriggerPatch* p = (AliEMCALTriggerPatch*)Next())
	165	{
	166	Int_t px, py, id; p->Position(px, py);
	167
	168	if (geom->GetAbsFastORIndexFromPositionInTRU(i, px, py, id)
	169	&&
	170	geom->GetPositionInEMCALFromAbsFastORIndex(id, px, py)) p->SetPosition(px, py);
	171
a61738e1	172	if (!data->GetMode()) // Simulation
de39a0ff	173	{
	174	Int_t peaks = p->Peaks();
	175
a61738e1	176	Int_t pos;
	177	AliEMCALTriggerRawDigit* dig = 0x0;
	178
de39a0ff	179	for (Int_t j = 0; j < sizeX * sizeY; j++)
	180	{
	181	if (peaks & (1 << j))
	182	{
a61738e1	183	pos = posMap[px + j % sizeX][py + j / sizeX];
de39a0ff	184
	185	if (pos == -1)
	186	{
a61738e1	187	// Add a new digit
de39a0ff	188	new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
	189
	190	dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
	191	}
	192	else
	193	{
	194	dig = (AliEMCALTriggerRawDigit*)digits->At(pos);
	195	}
	196
	197	dig->SetL0Time(p->Time());
	198	}
a61738e1	199	}
	200
	201	pos = posMap[px][py];
	202
	203	if (pos == -1)
	204	{
	205	// Add a new digit
	206	new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
	207
	208	dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
	209	}
	210	else
	211	{
	212	dig = (AliEMCALTriggerRawDigit*)digits->At(pos);
	213	}
	214
	215	dig->SetTriggerBit(kL0,0);
de39a0ff	216	}
	217	}
	218
	219	data->SetL0Trigger(0, i, 1);
de39a0ff	220	}
	221	else
	222	data->SetL0Trigger(0, i, 0);
916f1e76	223	}
	224
	225	// A L0 has been issued, run L1
de39a0ff	226	// Depending on raw data enabled or not in STU data: L1 computation
	227	// should be done from F-ALTRO or directly on TRU time sums in STU raw data
	228	if (iL0)
916f1e76	229	{
de39a0ff	230	// Use L1 threshold from raw data when reconstructing raw data
	231	if (data->GetMode())
	232	{
	233	fSTU->SetThreshold(kL1Gamma, data->GetL1GammaThreshold());
	234	fSTU->SetThreshold(kL1Jet, data->GetL1JetThreshold() );
	235	}
	236	else
	237	{
	238	fSTU->ComputeThFromV0(V0M); // C/A
	239	data->SetL1GammaThreshold(fSTU->GetThreshold(kL1Gamma));
	240	data->SetL1JetThreshold( fSTU->GetThreshold(kL1Jet) );
	241	}
916f1e76	242
de39a0ff	243	if (fSTU->GetRawData())
	244	{
	245	// Compute L1 from STU raw data
	246	fSTU->SetRegion(region);
	247	}
	248	else
	249	{
	250	// Build STU raw data from F-ALTRO
	251	TString str = "region";
	252	for (Int_t i = 0; i < kNTRU; i++) fSTU->Build(str,
	253	i,
	254	(static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->Region(),
	255	(static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->RegionSize());
	256	}
916f1e76	257
de39a0ff	258	fSTU->L1(kL1Gamma);
916f1e76	259
a61738e1	260	Int_t id, px, py;
	261	AliEMCALTriggerRawDigit* dig = 0x0;
	262
	263	TIterator* nP = 0x0;
	264
	265	nP = (fSTU->Patches()).MakeIterator();
	266
	267	while (AliEMCALTriggerPatch* p = (AliEMCALTriggerPatch*)nP->Next())
	268	{
	269	p->Position(px, py);
	270
	271	if (geom->GetAbsFastORIndexFromPositionInEMCAL(px, py, id))
	272	{
	273	if (posMap[px][py] == -1)
	274	{
	275	// Add a new digit
	276	new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
	277
	278	dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
	279	}
	280	else
	281	{
	282	dig = (AliEMCALTriggerRawDigit*)digits->At(posMap[px][py]);
	283	}
	284
	285	dig->SetTriggerBit(kL1Gamma,0);
	286	}
	287	}
916f1e76	288
	289	fSTU->Reset();
	290
de39a0ff	291	fSTU->L1(kL1Jet);
916f1e76	292
a61738e1	293	nP = (fSTU->Patches()).MakeIterator();
	294
	295	while (AliEMCALTriggerPatch* p = (AliEMCALTriggerPatch*)nP->Next())
	296	{
	297	p->Position(px, py);
	298
	299	px *= (Int_t)((fSTU->SubRegionSize())->X());
	300
	301	py *= (Int_t)((fSTU->SubRegionSize())->Y());
	302
	303	if (geom->GetAbsFastORIndexFromPositionInEMCAL(px, py, id))
	304	{
	305	if (posMap[px][py] == -1)
	306	{
	307	// Add a new digit
	308	new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
	309
	310	dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
	311	}
	312	else
	313	{
	314	dig = (AliEMCALTriggerRawDigit*)digits->At(posMap[px][py]);
	315	}
	316
	317	dig->SetTriggerBit(kL1Jet,0);
	318	}
	319	}
de39a0ff	320
	321	Int_t** reg = fSTU->Region();
	322
	323	if (!fSTU->GetRawData())
	324	{
	325	// Update digits w/ L1 time sum
	326	// Done in raw digit maker when raw data enabled
	327	for (Int_t i = 0; i < 48; i++)
	328	{
	329	for (Int_t j = 0; j < 64; j++)
	330	{
	331	if (reg[i][j])
	332	{
de39a0ff	333	if (geom->GetAbsFastORIndexFromPositionInEMCAL(i, j, id))
de39a0ff	334	{
de39a0ff	335	if (posMap[i][j] == -1)
	336	{
	337	// Add a new digit with L1 time sum
	338	new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
	339
	340	dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
	341	}
	342	else
	343	{
	344	dig = (AliEMCALTriggerRawDigit*)digits->At(posMap[i][j]);
	345	}
	346
	347	dig->SetL1TimeSum(reg[i][j]);
	348	}
	349	}
	350	}
	351	}
	352	}
916f1e76	353	}
916f1e76	354
de39a0ff	355	if (AliDebugLevel() >= 999) data->Scan();
39c05eac	356
916f1e76	357	// Now reset the electronics for a fresh start with next event
	358	Reset();
	359	}
	360
	361	//__________________
	362	void AliEMCALTriggerElectronics::Reset()
	363	{
	364	//
	365	TIter NextTRU(fTRU);
	366	while ( AliEMCALTriggerTRU TRU = (AliEMCALTriggerTRU)NextTRU() ) TRU->Reset();
	367
	368	fSTU->Reset();
	369	}