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

/* $Id: AliTRDtrackletGTU.cxx 28397 2008-09-02 09:33:00Z cblume $ */

////////////////////////////////////////////////////////////////////////////
//                                                                        //
//  GTU tracklet                                                          //
//                                                                        //
//  Author: J. Klein (Jochen.Klein@cern.ch)                               //
//                                                                        //
////////////////////////////////////////////////////////////////////////////

#include "TMath.h"
#include "TClass.h"

#include "AliTRDtrackletGTU.h"
#include "AliTRDtrackletWord.h"
#include "AliTRDmcmTracklet.h"
#include "AliTRDtrackletMCM.h"
#include "AliLog.h"
#include "AliTRDgtuParam.h"
#include "AliTRDgeometry.h"
#include "AliTRDpadPlane.h"

ClassImp(AliTRDtrackletGTU)

AliTRDtrackletBase* AliTRDtrackletGTU::fgkDummyTracklet = new AliTRDtrackletWord(0);

AliTRDtrackletGTU::AliTRDtrackletGTU() :
  AliTRDtrackletBase(),
  fGtuParam(AliTRDgtuParam::Instance()),
  fTracklet(fgkDummyTracklet), 
  fSubChannel(0x0),
  fAssignedZ(kFALSE),
  fAlpha(0),
  fYProj(0),
  fYPrime(0),
  fIndex(0)
{
  // ctor for any tracklet deriving from AliTRDtrackletBase

  fSubChannel = new Int_t[fGtuParam->GetNZChannels()];
  for (Int_t zch = 0; zch < fGtuParam->GetNZChannels(); zch++) 
    fSubChannel[zch] = 0;
}

AliTRDtrackletGTU::AliTRDtrackletGTU(AliTRDtrackletBase *tracklet) :
  AliTRDtrackletBase(*tracklet),
  fGtuParam(AliTRDgtuParam::Instance()),
  fTracklet(0x0), 
  fSubChannel(0x0),
  fAssignedZ(kFALSE),
  fAlpha(0),
  fYProj(0),
  fYPrime(0),
  fIndex(0)
{
  // ctor for any tracklet deriving from AliTRDtrackletBase

  fSubChannel = new Int_t[fGtuParam->GetNZChannels()];
  for (Int_t zch = 0; zch < fGtuParam->GetNZChannels(); zch++) 
    fSubChannel[zch] = 0;
  fTracklet = tracklet;
  if ( fTracklet->IsA() == TClass::GetClass("AliTRDtrackletMCM")) {
      AliInfo(Form("label from mcm tracklet: %i", ((AliTRDtrackletMCM*) fTracklet)->GetLabel()));
  }
}

AliTRDtrackletGTU::AliTRDtrackletGTU(const AliTRDtrackletGTU& tracklet) :
  AliTRDtrackletBase(tracklet),
  fGtuParam(AliTRDgtuParam::Instance()),
  fTracklet(tracklet.fTracklet), 
  fSubChannel(0x0),
  fAssignedZ(tracklet.fAssignedZ),
  fAlpha(tracklet.fAlpha),
  fYProj(tracklet.fYProj),
  fYPrime(tracklet.fYPrime), 
  fIndex(tracklet.fIndex)
{
  // copy ctor

  fSubChannel = new Int_t[fGtuParam->GetNZChannels()];
  for (Int_t zch = 0; zch < fGtuParam->GetNZChannels(); zch++) 
    fSubChannel[zch] = tracklet.fSubChannel[zch];
}

AliTRDtrackletGTU& AliTRDtrackletGTU::operator=(const AliTRDtrackletGTU &rhs)
{
  if (&rhs != this) {
    fTracklet = rhs.fTracklet;
    for (Int_t zch = 0; zch < fGtuParam->GetNZChannels(); zch++) 
      fSubChannel[zch] = rhs.fSubChannel[zch];
    fIndex = rhs.fIndex;
    fYPrime = rhs.fYPrime;
    fYProj = rhs.fYProj;
    fAlpha = rhs.fAlpha;
    fAssignedZ = rhs.fAssignedZ;
  }
  
  return *this;
}

AliTRDtrackletGTU::~AliTRDtrackletGTU() 
{
  // dtor
  if (fSubChannel)
    delete [] fSubChannel; 
}

Int_t AliTRDtrackletGTU::Compare(const TObject *o) const {
  // sorting w. r. t. Z, Y if no z-channel assigned (as needed in input unit)
  // otherwise w. r. t. Z, Y_proj
  // must be changed to Z-channel, Y_proj
  // will be changed

  if (!o) 
    return 0;

  if (!o->InheritsFrom("AliTRDtrackletGTU")) {
      AliError("Cannot compare to object not deriving from AliTRDtrackletGTU");
      return 0;
  }

  if (!fAssignedZ) {
    if ( GetZbin() < ((AliTRDtrackletGTU*) o)->GetZbin()) 
      return -1;
    else if (GetZbin() > ((AliTRDtrackletGTU*) o)->GetZbin()) 
      return 1;
    else 
      if (GetYbin() < ((AliTRDtrackletGTU*) o)->GetYbin())
	return -1;
      else if (GetYbin() > ((AliTRDtrackletGTU*) o)->GetYbin())
	return 1;
      else 
	return 0;
  }
  else {
    // sorting should be according to zsubindex, not to Z !!!
    // therefore this depends on the zch
    if (GetZbin() < ((AliTRDtrackletGTU*) o)->GetZbin()) 
      return -1;
    else if (GetZbin() > ((AliTRDtrackletGTU*) o)->GetZbin()) 
      return 1;
    else 
      if (GetYProj() < ((AliTRDtrackletGTU*) o)->GetYProj())
	return -1;
      else if (GetYProj() > ((AliTRDtrackletGTU*) o)->GetYProj())
	return 1;
      else 
	return 0;
  }
}

void AliTRDtrackletGTU::SetSubChannel(Int_t zch, Int_t subch) 
{
  // set the subchannel in the given z-channel
  fAssignedZ = kTRUE;
  fSubChannel[zch] = subch;
}

Int_t AliTRDtrackletGTU::GetSubChannel(Int_t zch) 
{
  // get the subchannel in the given z-channel
  return fSubChannel[zch];
}

Int_t AliTRDtrackletGTU::GetLabel() const
{
    if ( fTracklet->IsA() == TClass::GetClass("AliTRDtrackletMCM"))
	return ((AliTRDtrackletMCM*) fTracklet)->GetLabel();
    else
	return -1;
}

/*
Float_t AliTRDtrackletGTU::GetPhysX(Int_t layer) 
{
  // get the x-position (in the local system) assuming the tracklet is in the given layer
  return fGtuParam->GetGeo()->GetTime0(layer);
}

Float_t AliTRDtrackletGTU::GetPhysY() 
{
  // 
  return GetYbin() * 0.0160; 
}

Float_t AliTRDtrackletGTU::GetPhysAlpha() 
{
  return GetAlpha() * 0.01; // wrong factor!
}

Float_t AliTRDtrackletGTU::GetPhysZ(Int_t stack, Int_t layer) 
{
  return fGtuParam->GetGeo()->GetPadPlane(layer, stack)->GetRowPos(GetZbin()); // not the middle of a pad!
}
*/
Commit	Line	Data
52c19022	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	/* $Id: AliTRDtrackletGTU.cxx 28397 2008-09-02 09:33:00Z cblume $ */
	17
	18	////////////////////////////////////////////////////////////////////////////
	19	// //
	20	// GTU tracklet //
	21	// //
	22	// Author: J. Klein (Jochen.Klein@cern.ch) //
	23	// //
	24	////////////////////////////////////////////////////////////////////////////
	25
	26	#include "TMath.h"
4cc89512	27	#include "TClass.h"
52c19022	28
	29	#include "AliTRDtrackletGTU.h"
	30	#include "AliTRDtrackletWord.h"
	31	#include "AliTRDmcmTracklet.h"
4cc89512	32	#include "AliTRDtrackletMCM.h"
52c19022	33	#include "AliLog.h"
	34	#include "AliTRDgtuParam.h"
	35	#include "AliTRDgeometry.h"
	36	#include "AliTRDpadPlane.h"
	37
	38	ClassImp(AliTRDtrackletGTU)
	39
	40	AliTRDtrackletBase* AliTRDtrackletGTU::fgkDummyTracklet = new AliTRDtrackletWord(0);
	41
	42	AliTRDtrackletGTU::AliTRDtrackletGTU() :
	43	AliTRDtrackletBase(),
	44	fGtuParam(AliTRDgtuParam::Instance()),
	45	fTracklet(fgkDummyTracklet),
	46	fSubChannel(0x0),
	47	fAssignedZ(kFALSE),
	48	fAlpha(0),
	49	fYProj(0),
	50	fYPrime(0),
	51	fIndex(0)
	52	{
	53	// ctor for any tracklet deriving from AliTRDtrackletBase
	54
	55	fSubChannel = new Int_t[fGtuParam->GetNZChannels()];
	56	for (Int_t zch = 0; zch < fGtuParam->GetNZChannels(); zch++)
	57	fSubChannel[zch] = 0;
	58	}
	59
	60	AliTRDtrackletGTU::AliTRDtrackletGTU(AliTRDtrackletBase *tracklet) :
	61	AliTRDtrackletBase(*tracklet),
	62	fGtuParam(AliTRDgtuParam::Instance()),
	63	fTracklet(0x0),
	64	fSubChannel(0x0),
	65	fAssignedZ(kFALSE),
	66	fAlpha(0),
	67	fYProj(0),
	68	fYPrime(0),
	69	fIndex(0)
	70	{
	71	// ctor for any tracklet deriving from AliTRDtrackletBase
	72
	73	fSubChannel = new Int_t[fGtuParam->GetNZChannels()];
	74	for (Int_t zch = 0; zch < fGtuParam->GetNZChannels(); zch++)
	75	fSubChannel[zch] = 0;
	76	fTracklet = tracklet;
4cc89512	77	if ( fTracklet->IsA() == TClass::GetClass("AliTRDtrackletMCM")) {
	78	AliInfo(Form("label from mcm tracklet: %i", ((AliTRDtrackletMCM*) fTracklet)->GetLabel()));
	79	}
52c19022	80	}
	81
	82	AliTRDtrackletGTU::AliTRDtrackletGTU(const AliTRDtrackletGTU& tracklet) :
	83	AliTRDtrackletBase(tracklet),
	84	fGtuParam(AliTRDgtuParam::Instance()),
	85	fTracklet(tracklet.fTracklet),
	86	fSubChannel(0x0),
	87	fAssignedZ(tracklet.fAssignedZ),
	88	fAlpha(tracklet.fAlpha),
	89	fYProj(tracklet.fYProj),
	90	fYPrime(tracklet.fYPrime),
	91	fIndex(tracklet.fIndex)
	92	{
	93	// copy ctor
	94
	95	fSubChannel = new Int_t[fGtuParam->GetNZChannels()];
	96	for (Int_t zch = 0; zch < fGtuParam->GetNZChannels(); zch++)
	97	fSubChannel[zch] = tracklet.fSubChannel[zch];
	98	}
	99
	100	AliTRDtrackletGTU& AliTRDtrackletGTU::operator=(const AliTRDtrackletGTU &rhs)
	101	{
	102	if (&rhs != this) {
	103	fTracklet = rhs.fTracklet;
	104	for (Int_t zch = 0; zch < fGtuParam->GetNZChannels(); zch++)
	105	fSubChannel[zch] = rhs.fSubChannel[zch];
	106	fIndex = rhs.fIndex;
	107	fYPrime = rhs.fYPrime;
	108	fYProj = rhs.fYProj;
	109	fAlpha = rhs.fAlpha;
	110	fAssignedZ = rhs.fAssignedZ;
	111	}
	112
	113	return *this;
	114	}
	115
	116	AliTRDtrackletGTU::~AliTRDtrackletGTU()
	117	{
	118	// dtor
	119	if (fSubChannel)
	120	delete [] fSubChannel;
	121	}
	122
	123	Int_t AliTRDtrackletGTU::Compare(const TObject *o) const {
	124	// sorting w. r. t. Z, Y if no z-channel assigned (as needed in input unit)
	125	// otherwise w. r. t. Z, Y_proj
	126	// must be changed to Z-channel, Y_proj
	127	// will be changed
	128
	129	if (!o)
	130	return 0;
	131
	132	if (!o->InheritsFrom("AliTRDtrackletGTU")) {
	133	AliError("Cannot compare to object not deriving from AliTRDtrackletGTU");
	134	return 0;
	135	}
	136
	137	if (!fAssignedZ) {
	138	if ( GetZbin() < ((AliTRDtrackletGTU*) o)->GetZbin())
	139	return -1;
	140	else if (GetZbin() > ((AliTRDtrackletGTU*) o)->GetZbin())
	141	return 1;
	142	else
	143	if (GetYbin() < ((AliTRDtrackletGTU*) o)->GetYbin())
144	return -1;
145	else if (GetYbin() > ((AliTRDtrackletGTU*) o)->GetYbin())
146	return 1;
147	else
148	return 0;
149	}
150	else {
151	// sorting should be according to zsubindex, not to Z !!!
152	// therefore this depends on the zch
153	if (GetZbin() < ((AliTRDtrackletGTU*) o)->GetZbin())
154	return -1;
155	else if (GetZbin() > ((AliTRDtrackletGTU*) o)->GetZbin())
156	return 1;
157	else
158	if (GetYProj() < ((AliTRDtrackletGTU*) o)->GetYProj())
159	return -1;
160	else if (GetYProj() > ((AliTRDtrackletGTU*) o)->GetYProj())
161	return 1;
162	else
163	return 0;
164	}
165	}
166
167	void AliTRDtrackletGTU::SetSubChannel(Int_t zch, Int_t subch)
168	{
169	// set the subchannel in the given z-channel
170	fAssignedZ = kTRUE;
171	fSubChannel[zch] = subch;
172	}
173
174	Int_t AliTRDtrackletGTU::GetSubChannel(Int_t zch)
175	{
176	// get the subchannel in the given z-channel
177	return fSubChannel[zch];
178	}
179
4cc89512	180	Int_t AliTRDtrackletGTU::GetLabel() const
	181	{
	182	if ( fTracklet->IsA() == TClass::GetClass("AliTRDtrackletMCM"))
	183	return ((AliTRDtrackletMCM*) fTracklet)->GetLabel();
	184	else
	185	return -1;
	186	}
	187
52c19022	188	/*
	189	Float_t AliTRDtrackletGTU::GetPhysX(Int_t layer)
	190	{
	191	// get the x-position (in the local system) assuming the tracklet is in the given layer
	192	return fGtuParam->GetGeo()->GetTime0(layer);
	193	}
	194
	195	Float_t AliTRDtrackletGTU::GetPhysY()
	196	{
	197	//
	198	return GetYbin() * 0.0160;
	199	}
	200
	201	Float_t AliTRDtrackletGTU::GetPhysAlpha()
	202	{
	203	return GetAlpha() * 0.01; // wrong factor!
	204	}
	205
	206	Float_t AliTRDtrackletGTU::GetPhysZ(Int_t stack, Int_t layer)
	207	{
	208	return fGtuParam->GetGeo()->GetPadPlane(layer, stack)->GetRowPos(GetZbin()); // not the middle of a pad!
	209	}
	210	*/