[u/mrichter/AliRoot.git] / HLT / MUON / AliHLTMUONTriggerRecord.cxx

/**************************************************************************
 * This file is property of and copyright by the ALICE HLT Project        * 
 * All rights reserved.                                                   *
 *                                                                        *
 * Primary Authors:                                                       *
 *   Artur Szostak <artursz@iafrica.com>                                  *
 *                                                                        *
 * 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.                  *
 **************************************************************************/

/* $Id$ */

/**
 * @file   AliHLTMUONTriggerRecord.cxx
 * @author Artur Szostak <artursz@iafrica.com>
 * @date   
 * @brief  Implementation of the AliHLTMUONTriggerRecord class.
 */

#include "AliHLTMUONTriggerRecord.h"
#include "AliLog.h"
#include "mapping/AliMpDEManager.h"
#include <cstring>
#include <iostream>
#include <iomanip>
using namespace std;

ClassImp(AliHLTMUONTriggerRecord);


std::ostream& operator << (
		std::ostream& stream,
		const AliHLTMUONTriggerRecord& trigrec
	)
{
// Stream operator.

	stream	<< "ID: " << trigrec.fId
		<< "; sign: " << trigrec.fSign
		<< "; p = (" << trigrec.Px()
		<< ", " << trigrec.Py()
		<< ", " << trigrec.Pz()
		<< ")";
	return stream;
}


AliHLTMUONTriggerRecord::AliHLTMUONTriggerRecord(
		Int_t id, Int_t sign,
		Float_t px, Float_t py, Float_t pz,
		Int_t sourceDDL, Float_t zf, Float_t qbl
	) :
	fId(id), fSign(sign), fMomentum(px, py, pz), fSourceDDL(sourceDDL),
	fZmiddle(zf), fQBL(qbl)
{
	if (sign < -1 or 1 < sign)
	{
		AliError(Form("Trying to set the sign to %d. This is outside the"
			" valid range of [-1..1]", sign
		));
		fSign = 0;
	}

	// Fill the debugging information to invalid values by default.
	for (int i = 0; i < 4; i++)
	{
		fDetElemId[i] = -1;
		fPatternX[i] = -1;
		fPatternY[i] = -1;
	}
}


const TVector3& AliHLTMUONTriggerRecord::Hit(Int_t chamber) const
{
// Returns the hit on the specified chamber.

	if (11 <= chamber and chamber <= 14) return fHit[chamber - 11];
	
	AliError(Form(
		"Chamber number %d is not in the valid range [11..14].",
		int(chamber)
	));
	return fHit[0];
}


Int_t AliHLTMUONTriggerRecord::DetElemId(Int_t chamber) const
{
// Returns the detector element ID for the specified chamber associated
// to the hit on that chamber.

	if (11 <= chamber and chamber <= 14) return fDetElemId[chamber - 11];
	
	AliError(Form(
		"Chamber number %d is not in the valid range [11..14].",
		int(chamber)
	));
	return fDetElemId[0];
}


Int_t AliHLTMUONTriggerRecord::PatternX(Int_t chamber) const
{
// Returns the raw data X pattern of the hit on the specified chamber.

	if (11 <= chamber and chamber <= 14) return fPatternX[chamber - 11];
	
	AliError(Form(
		"Chamber number %d is not in the valid range [11..14].",
		int(chamber)
	));
	return fPatternX[0];
}


Int_t AliHLTMUONTriggerRecord::PatternY(Int_t chamber) const
{
// Returns the raw data Y pattern of the hit on the specified chamber.

	if (11 <= chamber and chamber <= 14) return fPatternY[chamber - 11];
	
	AliError(Form(
		"Chamber number %d is not in the valid range [11..14].",
		int(chamber)
	));
	return fPatternY[0];
}


void AliHLTMUONTriggerRecord::SetHit(Int_t chamber, Float_t x, Float_t y, Float_t z)
{
// Fills the hit coordinate on the specified chamber.

	if (11 <= chamber and chamber <= 14)
	{
		fHit[chamber - 11].SetXYZ(x, y, z);
	}
	else
	{
		AliError(Form(
			"Chamber number %d is not in the valid range [11..14].",
			int(chamber)
		));
	}
}


void AliHLTMUONTriggerRecord::SetHitDebugInfo(
		Int_t chamber,
		Int_t detElemId, UShort_t patternX, UShort_t patternY
	)
{
// Fills the debugging information corresponding to the hit on the specified chamber.

	if (11 <= chamber and chamber <= 14)
	{
		fDetElemId[chamber - 11] = detElemId;
		fPatternX[chamber - 11] = patternX;
		fPatternY[chamber - 11] = patternY;
	}
	else
	{
		AliError(Form(
			"Chamber number %d is not in the valid range [11..14].",
			int(chamber)
		));
	}
}


void AliHLTMUONTriggerRecord::Print(Option_t* option) const
{
// Prints the trigger record to standard output (screen).

	if (	option == NULL or strcmp(option, "") == 0 or
		strcmp(option, "compact") == 0
	   )
	{
		cout << *this << endl;
	}
	else if (strcmp(option, "detail") == 0)
	{
		cout << "Trigger record ID = " << fId << "; sign = ";
		if (fSign != 0)
			cout << fSign;
		else
			cout << "unknown";
		cout << "; momentum: (px = " << Px()
			<< " GeV/c, py = " << Py()
			<< " GeV/c, pz = " << Pz()
			<< " GeV/c)" << endl;
		cout << "Used Zmiddle = " << fZmiddle << " cm and QBL = "
			<< fQBL << " T.m for the momentum calculation." << endl;
	}
	else if (strcmp(option, "all") == 0)
	{
		cout << "Trigger record ID = " << fId << "; sign = ";
		if (fSign != 0)
			cout << fSign;
		else
			cout << "unknown";
		cout << "; momentum: (px = " << Px()
			<< " GeV/c, py = " << Py()
			<< " GeV/c, pz = " << Pz()
			<< " GeV/c)" << endl;
		cout << "Used Zmiddle = " << fZmiddle << " cm and QBL = "
			<< fQBL << " T.m for the momentum calculation." << endl;
		
		streamsize w = cout.width();
		ios::fmtflags f = cout.flags();
		cout << setw(9) << "Chamber" << setw(0) << "  Hit:"
			<< setw(8) << "X (cm)"
			<< setw(12) << "Y (cm)"
			<< setw(12) << "Z (cm)"
			<< setw(12) << "DetElemID"
			<< setw(18) << "X bit pattern"
			<< setw(18) << "Y bit pattern" << endl;
		for (int i = 0; i < 4; i++)
		{
			cout << setw(9) << i+11
				<< setw(14) << fHit[i].X()
				<< setw(12) << fHit[i].Y()
				<< setw(12) << fHit[i].Z()
				<< setw(12) << fDetElemId[i];
			if (fPatternX[i] != -1)
			{
				// Print the X pattern as a bit pattern.
				cout << "  ";
				for (Int_t j = 15; j >= 0; j--)
					cout << (((fPatternX[i] & (1 << j)) > 0) ? "1" : "0");
			}
			else
			{
				cout << "  ----------------";
			}
			if (fPatternY[i] != -1)
			{
				// Print the Y pattern as a bit pattern.
				cout << "  ";
				for (Int_t j = 15; j >= 0; j--)
					cout << (((fPatternY[i] & (1 << j)) > 0) ? "1" : "0");
			}
			else
			{
				cout << "  ----------------";
			}
			cout << endl;
		}
		cout.width(w); // reset the field width to previous value.
		cout.flags(f); // reset the flags to previous values.
	}
	else
	{
		AliError("Unknown option specified. Can only be one of 'compact',"
			" 'detail' or 'all'."
		);
	}
}


Int_t AliHLTMUONTriggerRecord::Compare(const TObject* obj) const
{
// We compare this object with 'obj' first by trigger record ID, then
// by sign and finally by momentum.

	if (obj->IsA() == AliHLTMUONTriggerRecord::Class())
	{
		const AliHLTMUONTriggerRecord* tr =
			static_cast<const AliHLTMUONTriggerRecord*>(obj);
		if (fId < tr->fId) return -1;
		if (fId > tr->fId) return 1;
		if (fSign < tr->fSign) return -1;
		if (fSign > tr->fSign) return 1;
		if (Px() < tr->Px()) return -1;
		if (Px() > tr->Px()) return 1;
		if (Py() < tr->Py()) return -1;
		if (Py() > tr->Py()) return 1;
		if (Pz() < tr->Pz()) return -1;
		if (Pz() > tr->Pz()) return 1;
		return 0;
	}
	else
	{
		AliError(Form("Do not know how to compare %s to %s.",
			this->ClassName(),
			obj->ClassName()
		));
		return -999;
	}
}


bool AliHLTMUONTriggerRecord::operator == (const AliHLTMUONTriggerRecord& trigrec) const
{
// Compares the trigger record 'trigrec' to this one.

	return	fId == trigrec.fId and fSign == trigrec.fSign
		and fMomentum == trigrec.fMomentum
		and fHit[0] == trigrec.fHit[0] and fHit[1] == trigrec.fHit[1]
		and fHit[2] == trigrec.fHit[2] and fHit[3] == trigrec.fHit[3];
}
Commit	Line	Data
753896cd	1	/**************************************************************************
	2	* This file is property of and copyright by the ALICE HLT Project *
	3	* All rights reserved. *
	4	* *
	5	* Primary Authors: *
	6	* Artur Szostak <artursz@iafrica.com> *
	7	* *
	8	* Permission to use, copy, modify and distribute this software and its *
	9	* documentation strictly for non-commercial purposes is hereby granted *
	10	* without fee, provided that the above copyright notice appears in all *
	11	* copies and that both the copyright notice and this permission notice *
	12	* appear in the supporting documentation. The authors make no claims *
	13	* about the suitability of this software for any purpose. It is *
	14	* provided "as is" without express or implied warranty. *
	15	**************************************************************************/
	16
	17	/* $Id$ */
	18
	19	/**
	20	* @file AliHLTMUONTriggerRecord.cxx
	21	* @author Artur Szostak <artursz@iafrica.com>
	22	* @date
	23	* @brief Implementation of the AliHLTMUONTriggerRecord class.
	24	*/
	25
	26	#include "AliHLTMUONTriggerRecord.h"
	27	#include "AliLog.h"
	28	#include "mapping/AliMpDEManager.h"
	29	#include <cstring>
	30	#include <iostream>
	31	#include <iomanip>
	32	using namespace std;
	33
	34	ClassImp(AliHLTMUONTriggerRecord);
	35
	36
	37	std::ostream& operator << (
	38	std::ostream& stream,
	39	const AliHLTMUONTriggerRecord& trigrec
	40	)
	41	{
	42	// Stream operator.
	43
	44	stream << "ID: " << trigrec.fId
	45	<< "; sign: " << trigrec.fSign
	46	<< "; p = (" << trigrec.Px()
	47	<< ", " << trigrec.Py()
	48	<< ", " << trigrec.Pz()
	49	<< ")";
	50	return stream;
	51	}
	52
	53
	54	AliHLTMUONTriggerRecord::AliHLTMUONTriggerRecord(
	55	Int_t id, Int_t sign,
	56	Float_t px, Float_t py, Float_t pz,
	57	Int_t sourceDDL, Float_t zf, Float_t qbl
	58	) :
	59	fId(id), fSign(sign), fMomentum(px, py, pz), fSourceDDL(sourceDDL),
	60	fZmiddle(zf), fQBL(qbl)
	61	{
	62	if (sign < -1 or 1 < sign)
	63	{
	64	AliError(Form("Trying to set the sign to %d. This is outside the"
65	" valid range of [-1..1]", sign
66	));
67	fSign = 0;
68	}
69
70	// Fill the debugging information to invalid values by default.
71	for (int i = 0; i < 4; i++)
72	{
73	fDetElemId[i] = -1;
74	fPatternX[i] = -1;
75	fPatternY[i] = -1;
76	}
77	}
78
79
80	const TVector3& AliHLTMUONTriggerRecord::Hit(Int_t chamber) const
81	{
82	// Returns the hit on the specified chamber.
83
84	if (11 <= chamber and chamber <= 14) return fHit[chamber - 11];
85
86	AliError(Form(
87	"Chamber number %d is not in the valid range [11..14].",
88	int(chamber)
89	));
90	return fHit[0];
91	}
92
93
94	Int_t AliHLTMUONTriggerRecord::DetElemId(Int_t chamber) const
95	{
96	// Returns the detector element ID for the specified chamber associated
97	// to the hit on that chamber.
98
99	if (11 <= chamber and chamber <= 14) return fDetElemId[chamber - 11];
100
101	AliError(Form(
102	"Chamber number %d is not in the valid range [11..14].",
103	int(chamber)
104	));
105	return fDetElemId[0];
106	}
107
108
109	Int_t AliHLTMUONTriggerRecord::PatternX(Int_t chamber) const
110	{
111	// Returns the raw data X pattern of the hit on the specified chamber.
112
113	if (11 <= chamber and chamber <= 14) return fPatternX[chamber - 11];
114
115	AliError(Form(
116	"Chamber number %d is not in the valid range [11..14].",
117	int(chamber)
118	));
119	return fPatternX[0];
120	}
121
122
123	Int_t AliHLTMUONTriggerRecord::PatternY(Int_t chamber) const
124	{
125	// Returns the raw data Y pattern of the hit on the specified chamber.
126
127	if (11 <= chamber and chamber <= 14) return fPatternY[chamber - 11];
128
129	AliError(Form(
130	"Chamber number %d is not in the valid range [11..14].",
131	int(chamber)
132	));
133	return fPatternY[0];
134	}
135
136
137	void AliHLTMUONTriggerRecord::SetHit(Int_t chamber, Float_t x, Float_t y, Float_t z)
138	{
139	// Fills the hit coordinate on the specified chamber.
140
141	if (11 <= chamber and chamber <= 14)
142	{
143	fHit[chamber - 11].SetXYZ(x, y, z);
144	}
145	else
146	{
147	AliError(Form(
148	"Chamber number %d is not in the valid range [11..14].",
149	int(chamber)
150	));
151	}
152	}
153
154
155	void AliHLTMUONTriggerRecord::SetHitDebugInfo(
156	Int_t chamber,
157	Int_t detElemId, UShort_t patternX, UShort_t patternY
158	)
159	{
160	// Fills the debugging information corresponding to the hit on the specified chamber.
161
162	if (11 <= chamber and chamber <= 14)
163	{
164	fDetElemId[chamber - 11] = detElemId;
165	fPatternX[chamber - 11] = patternX;
166	fPatternY[chamber - 11] = patternY;
167	}
168	else
169	{
170	AliError(Form(
171	"Chamber number %d is not in the valid range [11..14].",
172	int(chamber)
173	));
174	}
175	}
176
177
178	void AliHLTMUONTriggerRecord::Print(Option_t* option) const
179	{
180	// Prints the trigger record to standard output (screen).
181
182	if ( option == NULL or strcmp(option, "") == 0 or
183	strcmp(option, "compact") == 0
184	)
185	{
186	cout << *this << endl;
187	}
188	else if (strcmp(option, "detail") == 0)
189	{
190	cout << "Trigger record ID = " << fId << "; sign = ";
191	if (fSign != 0)
192	cout << fSign;
193	else
194	cout << "unknown";
195	cout << "; momentum: (px = " << Px()
196	<< " GeV/c, py = " << Py()
197	<< " GeV/c, pz = " << Pz()
198	<< " GeV/c)" << endl;
199	cout << "Used Zmiddle = " << fZmiddle << " cm and QBL = "
200	<< fQBL << " T.m for the momentum calculation." << endl;
201	}
202	else if (strcmp(option, "all") == 0)
203	{
204	cout << "Trigger record ID = " << fId << "; sign = ";
205	if (fSign != 0)
206	cout << fSign;
207	else
208	cout << "unknown";
209	cout << "; momentum: (px = " << Px()
210	<< " GeV/c, py = " << Py()
211	<< " GeV/c, pz = " << Pz()
212	<< " GeV/c)" << endl;
213	cout << "Used Zmiddle = " << fZmiddle << " cm and QBL = "
214	<< fQBL << " T.m for the momentum calculation." << endl;
215
216	streamsize w = cout.width();
217	ios::fmtflags f = cout.flags();
218	cout << setw(9) << "Chamber" << setw(0) << " Hit:"
219	<< setw(8) << "X (cm)"
220	<< setw(12) << "Y (cm)"
221	<< setw(12) << "Z (cm)"
222	<< setw(12) << "DetElemID"
223	<< setw(18) << "X bit pattern"
224	<< setw(18) << "Y bit pattern" << endl;
225	for (int i = 0; i < 4; i++)
226	{
227	cout << setw(9) << i+11
228	<< setw(14) << fHit[i].X()
229	<< setw(12) << fHit[i].Y()
230	<< setw(12) << fHit[i].Z()
231	<< setw(12) << fDetElemId[i];
232	if (fPatternX[i] != -1)
233	{
234	// Print the X pattern as a bit pattern.
235	cout << " ";
236	for (Int_t j = 15; j >= 0; j--)
237	cout << (((fPatternX[i] & (1 << j)) > 0) ? "1" : "0");
238	}
239	else
240	{
241	cout << " ----------------";
242	}
243	if (fPatternY[i] != -1)
244	{
245	// Print the Y pattern as a bit pattern.
246	cout << " ";
247	for (Int_t j = 15; j >= 0; j--)
248	cout << (((fPatternY[i] & (1 << j)) > 0) ? "1" : "0");
249	}
250	else
251	{
252	cout << " ----------------";
253	}
254	cout << endl;
255	}
256	cout.width(w); // reset the field width to previous value.
257	cout.flags(f); // reset the flags to previous values.
258	}
259	else
260	{
261	AliError("Unknown option specified. Can only be one of 'compact',"
262	" 'detail' or 'all'."
263	);
264	}
265	}
266
267
268	Int_t AliHLTMUONTriggerRecord::Compare(const TObject* obj) const
269	{
270	// We compare this object with 'obj' first by trigger record ID, then
271	// by sign and finally by momentum.
272
273	if (obj->IsA() == AliHLTMUONTriggerRecord::Class())
274	{
275	const AliHLTMUONTriggerRecord* tr =
276	static_cast<const AliHLTMUONTriggerRecord*>(obj);
277	if (fId < tr->fId) return -1;
278	if (fId > tr->fId) return 1;
279	if (fSign < tr->fSign) return -1;
280	if (fSign > tr->fSign) return 1;
281	if (Px() < tr->Px()) return -1;
282	if (Px() > tr->Px()) return 1;
283	if (Py() < tr->Py()) return -1;
284	if (Py() > tr->Py()) return 1;
285	if (Pz() < tr->Pz()) return -1;
286	if (Pz() > tr->Pz()) return 1;
287	return 0;
288	}
289	else
290	{
291	AliError(Form("Do not know how to compare %s to %s.",
292	this->ClassName(),
293	obj->ClassName()
294	));
295	return -999;
296	}
297	}
298
299
300	bool AliHLTMUONTriggerRecord::operator == (const AliHLTMUONTriggerRecord& trigrec) const
301	{
302	// Compares the trigger record 'trigrec' to this one.
303
304	return fId == trigrec.fId and fSign == trigrec.fSign
305	and fMomentum == trigrec.fMomentum
306	and fHit[0] == trigrec.fHit[0] and fHit[1] == trigrec.fHit[1]
307	and fHit[2] == trigrec.fHit[2] and fHit[3] == trigrec.fHit[3];
308	}