[u/mrichter/AliRoot.git] / HLT / MUON / AliHLTMUONRecHit.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   AliHLTMUONRecHit.cxx
 * @author Artur Szostak <artursz@iafrica.com>
 * @date   
 * @brief  Implementation of the AliHLTMUONRecHit class.
 */

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

ClassImp(AliHLTMUONRecHit);
ClassImp(AliHLTMUONRecHit::Channel);


std::ostream& operator << (std::ostream& stream, const AliHLTMUONRecHit& hit)
{
// Stream operator for std::ostream classes.

	stream << "(" << hit.X() << ", " << hit.Y() << ", " << hit.Z() << ")";
	return stream;
}


void AliHLTMUONRecHit::SetDebugInfo(
		Int_t detElemId, Int_t clusterId, UInt_t nChExp, Int_t sourceDDL
	)
{
// Sets the debugging information.

	fSourceDDL = sourceDDL;
	fDetElemId = detElemId;
	fClusterId = clusterId;
	fNchExp = nChExp;
}
	

Int_t AliHLTMUONRecHit::Chamber(bool warn) const
{
// Returns the chamber ID for this hit.

	if (fSourceDDL != -1)
	{
		if (fSourceDDL < 1 or fSourceDDL > 20)
		{
			return ((fSourceDDL-1) / 2) + 1;
		}
		else if (warn)
		{
			AliError(Form("The DDL source number: %d is out of range."
				" Valid values are [1..20]", fSourceDDL
			));
		}
	}
	if (fDetElemId != -1) return AliMpDEManager::GetChamberId(fDetElemId, warn);
	
	if (warn)
	{
		AliWarning("Neither the DDL source nor the detector element ID was not set,"
			" so we do not know on which chamber this hit was reconstructed."
		);
	}
	return -1;
}


void AliHLTMUONRecHit::AddChannel(
		Short_t manu, Short_t channel, Short_t signal,
		UInt_t rawDataWord
	)
{
// Adds a new channel to the channels list forming this hit's cluster.

	Int_t index = fChannels.GetEntriesFast();
	new (fChannels[index]) Channel(manu, channel, signal, rawDataWord);
}


void AliHLTMUONRecHit::Print(Option_t* option) const
{
// Prints the coordinates of this hit 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 << "(x = " << X() << " cm, y = " << Y()
			<< " cm, z = " << Z()
			<< " cm); source DDL = " << fSourceDDL
			<< "; DetElemID = " << fDetElemId
			<< "; cluster ID = " << fClusterId
			<< "; expected #ch = " << fNchExp << endl;
	}
	else if (strcmp(option, "all") == 0)
	{
		cout << "(x = " << X() << " cm, y = " << Y()
			<< " cm, z = " << Z()
			<< " cm); source DDL = " << fSourceDDL
			<< "; DetElemID = " << fDetElemId
			<< "; cluster ID = " << fClusterId
			<< "; expected #ch = " << fNchExp << endl;
		if (fChannels.GetEntriesFast() == 0)
		{
			cout << "No channels found for this hit." << endl;
		}
		else
		{
			streamsize w = cout.width();
			ios::fmtflags f = cout.flags();
			cout << setw(12) << "MANU"
				<< setw(12) << "Channel"
				<< setw(12) << "Signal"
				<< setw(15) << "Raw data word" << endl;
			cout << showbase;
			for (Int_t i = 0; i < fChannels.GetEntriesFast(); i++)
			{
				const AliHLTMUONRecHit::Channel* c =
					static_cast<const AliHLTMUONRecHit::Channel*>(fChannels[i]);
				cout << dec << setw(12) << c->Manu()
					<< setw(12) << c->Address()
					<< setw(12) << c->Signal()
					<< "     " << hex << setw(10) << internal;
				ios::char_type fc = cout.fill('0');
				cout << c->RawDataWord() << right << endl;
				cout.fill(fc);
			}
			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 AliHLTMUONRecHit::Compare(const TObject* obj) const
{
// We compare this object with 'obj' first by X, then Y, then Z.

	if (obj->IsA() == AliHLTMUONRecHit::Class())
	{
		const AliHLTMUONRecHit* h = static_cast<const AliHLTMUONRecHit*>(obj);
		if (X() < h->X()) return -1;
		if (X() > h->X()) return 1;
		if (Y() < h->Y()) return -1;
		if (Y() > h->Y()) return 1;
		if (Z() < h->Z()) return -1;
		if (Z() > h->Z()) return 1;
		return 0;
	}
	else
	{
		AliError(Form("Do not know how to compare %s to %s.",
			this->ClassName(),
			obj->ClassName()
		));
		return -999;
	}
}


std::ostream& operator << (std::ostream& stream, const AliHLTMUONRecHit::Channel& c)
{
// Stream operator for std::ostream classes.

	stream << "Channel: " << c.fManu << " , " << c.fAddress
		<< "; ADC: " << c.fSignal;
	return stream;
}


void AliHLTMUONRecHit::Channel::Print(Option_t* option) const
{
// Prints the details of this channel 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)
	{
		streamsize w = cout.width();
		ios::fmtflags f = cout.flags();
		cout << "MANU = " << fManu << ", Channel address = " << fAddress
			<< ", Signal = " << fSignal
			<< "; Raw data word = " << hex << setw(10) << internal;
		ios::char_type fc = cout.fill('0');
		cout << fRawDataWord << endl;
		cout.fill(fc); // reset fill character
		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' or 'detail'."
		);
	}
}


Int_t AliHLTMUONRecHit::Channel::Compare(const TObject* obj) const
{
// We compare this object with 'obj' first by MANU number, then by MANU channel
// address, then ADC signal.

	if (obj->IsA() == Channel::Class())
	{
		const Channel* c = static_cast<const Channel*>(obj);
		if (fManu < c->Manu()) return -1;
		if (fManu > c->Manu()) return 1;
		if (fAddress < c->Address()) return -1;
		if (fAddress > c->Address()) return 1;
		if (fSignal < c->Signal()) return -1;
		if (fSignal > c->Signal()) return 1;
		return 0;
	}
	else
	{
		AliError(Form("Do not know how to compare %s to %s.",
			this->ClassName(),
			obj->ClassName()
		));
		return -999;
	}
}
Commit	Line	Data
	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 AliHLTMUONRecHit.cxx
	21	* @author Artur Szostak <artursz@iafrica.com>
	22	* @date
	23	* @brief Implementation of the AliHLTMUONRecHit class.
	24	*/
	25
	26	#include "AliHLTMUONRecHit.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(AliHLTMUONRecHit);
	35	ClassImp(AliHLTMUONRecHit::Channel);
	36
	37
	38	std::ostream& operator << (std::ostream& stream, const AliHLTMUONRecHit& hit)
	39	{
	40	// Stream operator for std::ostream classes.
	41
	42	stream << "(" << hit.X() << ", " << hit.Y() << ", " << hit.Z() << ")";
	43	return stream;
	44	}
	45
	46
	47	void AliHLTMUONRecHit::SetDebugInfo(
	48	Int_t detElemId, Int_t clusterId, UInt_t nChExp, Int_t sourceDDL
	49	)
	50	{
	51	// Sets the debugging information.
	52
	53	fSourceDDL = sourceDDL;
	54	fDetElemId = detElemId;
	55	fClusterId = clusterId;
	56	fNchExp = nChExp;
	57	}
	58
	59
	60	Int_t AliHLTMUONRecHit::Chamber(bool warn) const
	61	{
	62	// Returns the chamber ID for this hit.
	63
	64	if (fSourceDDL != -1)
	65	{
	66	if (fSourceDDL < 1 or fSourceDDL > 20)
	67	{
	68	return ((fSourceDDL-1) / 2) + 1;
	69	}
	70	else if (warn)
	71	{
	72	AliError(Form("The DDL source number: %d is out of range."
	73	" Valid values are [1..20]", fSourceDDL
	74	));
	75	}
	76	}
	77	if (fDetElemId != -1) return AliMpDEManager::GetChamberId(fDetElemId, warn);
	78
	79	if (warn)
	80	{
	81	AliWarning("Neither the DDL source nor the detector element ID was not set,"
	82	" so we do not know on which chamber this hit was reconstructed."
	83	);
	84	}
	85	return -1;
	86	}
	87
	88
	89	void AliHLTMUONRecHit::AddChannel(
	90	Short_t manu, Short_t channel, Short_t signal,
	91	UInt_t rawDataWord
	92	)
	93	{
	94	// Adds a new channel to the channels list forming this hit's cluster.
	95
	96	Int_t index = fChannels.GetEntriesFast();
	97	new (fChannels[index]) Channel(manu, channel, signal, rawDataWord);
	98	}
	99
	100
	101	void AliHLTMUONRecHit::Print(Option_t* option) const
	102	{
	103	// Prints the coordinates of this hit to standard output (screen).
	104
	105	if ( option == NULL or strcmp(option, "") == 0 or
	106	strcmp(option, "compact") == 0
	107	)
	108	{
	109	cout << *this << endl;
	110	}
	111	else if (strcmp(option, "detail") == 0)
	112	{
	113	cout << "(x = " << X() << " cm, y = " << Y()
	114	<< " cm, z = " << Z()
	115	<< " cm); source DDL = " << fSourceDDL
	116	<< "; DetElemID = " << fDetElemId
	117	<< "; cluster ID = " << fClusterId
	118	<< "; expected #ch = " << fNchExp << endl;
	119	}
	120	else if (strcmp(option, "all") == 0)
	121	{
	122	cout << "(x = " << X() << " cm, y = " << Y()
	123	<< " cm, z = " << Z()
	124	<< " cm); source DDL = " << fSourceDDL
	125	<< "; DetElemID = " << fDetElemId
	126	<< "; cluster ID = " << fClusterId
	127	<< "; expected #ch = " << fNchExp << endl;
	128	if (fChannels.GetEntriesFast() == 0)
	129	{
	130	cout << "No channels found for this hit." << endl;
	131	}
	132	else
	133	{
	134	streamsize w = cout.width();
	135	ios::fmtflags f = cout.flags();
	136	cout << setw(12) << "MANU"
	137	<< setw(12) << "Channel"
	138	<< setw(12) << "Signal"
	139	<< setw(15) << "Raw data word" << endl;
	140	cout << showbase;
	141	for (Int_t i = 0; i < fChannels.GetEntriesFast(); i++)
	142	{
	143	const AliHLTMUONRecHit::Channel* c =
	144	static_cast<const AliHLTMUONRecHit::Channel*>(fChannels[i]);
	145	cout << dec << setw(12) << c->Manu()
	146	<< setw(12) << c->Address()
	147	<< setw(12) << c->Signal()
	148	<< " " << hex << setw(10) << internal;
	149	ios::char_type fc = cout.fill('0');
	150	cout << c->RawDataWord() << right << endl;
	151	cout.fill(fc);
	152	}
	153	cout.width(w); // reset the field width to previous value.
	154	cout.flags(f); // reset the flags to previous values.
	155	}
	156	}
	157	else
	158	{
	159	AliError("Unknown option specified. Can only be one of 'compact',"
	160	" 'detail' or 'all'."
	161	);
	162	}
	163	}
	164
	165
	166	Int_t AliHLTMUONRecHit::Compare(const TObject* obj) const
	167	{
	168	// We compare this object with 'obj' first by X, then Y, then Z.
	169
	170	if (obj->IsA() == AliHLTMUONRecHit::Class())
	171	{
	172	const AliHLTMUONRecHit* h = static_cast<const AliHLTMUONRecHit*>(obj);
	173	if (X() < h->X()) return -1;
	174	if (X() > h->X()) return 1;
	175	if (Y() < h->Y()) return -1;
	176	if (Y() > h->Y()) return 1;
	177	if (Z() < h->Z()) return -1;
	178	if (Z() > h->Z()) return 1;
	179	return 0;
	180	}
	181	else
	182	{
	183	AliError(Form("Do not know how to compare %s to %s.",
	184	this->ClassName(),
	185	obj->ClassName()
	186	));
	187	return -999;
	188	}
	189	}
	190
	191
	192	std::ostream& operator << (std::ostream& stream, const AliHLTMUONRecHit::Channel& c)
	193	{
	194	// Stream operator for std::ostream classes.
	195
	196	stream << "Channel: " << c.fManu << " , " << c.fAddress
	197	<< "; ADC: " << c.fSignal;
	198	return stream;
	199	}
	200
	201
	202	void AliHLTMUONRecHit::Channel::Print(Option_t* option) const
	203	{
	204	// Prints the details of this channel to standard output (screen).
	205
	206	if ( option == NULL or strcmp(option, "") == 0 or
	207	strcmp(option, "compact") == 0
	208	)
	209	{
	210	cout << *this << endl;
	211	}
	212	else if (strcmp(option, "detail") == 0)
	213	{
	214	streamsize w = cout.width();
	215	ios::fmtflags f = cout.flags();
	216	cout << "MANU = " << fManu << ", Channel address = " << fAddress
	217	<< ", Signal = " << fSignal
	218	<< "; Raw data word = " << hex << setw(10) << internal;
	219	ios::char_type fc = cout.fill('0');
	220	cout << fRawDataWord << endl;
	221	cout.fill(fc); // reset fill character
	222	cout.width(w); // reset the field width to previous value.
	223	cout.flags(f); // reset the flags to previous values.
	224	}
	225	else
	226	{
	227	AliError("Unknown option specified. Can only be one of"
	228	" 'compact' or 'detail'."
	229	);
	230	}
	231	}
	232
	233
	234	Int_t AliHLTMUONRecHit::Channel::Compare(const TObject* obj) const
	235	{
	236	// We compare this object with 'obj' first by MANU number, then by MANU channel
	237	// address, then ADC signal.
	238
	239	if (obj->IsA() == Channel::Class())
	240	{
	241	const Channel* c = static_cast<const Channel*>(obj);
	242	if (fManu < c->Manu()) return -1;
	243	if (fManu > c->Manu()) return 1;
	244	if (fAddress < c->Address()) return -1;
	245	if (fAddress > c->Address()) return 1;
	246	if (fSignal < c->Signal()) return -1;
	247	if (fSignal > c->Signal()) return 1;
	248	return 0;
	249	}
	250	else
	251	{
	252	AliError(Form("Do not know how to compare %s to %s.",
	253	this->ClassName(),
	254	obj->ClassName()
	255	));
	256	return -999;
	257	}
	258	}