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

/*
$Log$
Revision 1.1.2.2  2000/10/04 16:34:58  cblume
Replace include files by forward declarations

Revision 1.1.2.1  2000/09/22 14:47:52  cblume
Add the tracking code

*/                        
                                
#include "AliTRDcluster.h" 
#include "AliTRDconst.h"
#include "AliTRDtimeBin.h" 

ClassImp(AliTRDtimeBin)

//______________________________________________________

void AliTRDtimeBin::InsertCluster(AliTRDcluster* c, UInt_t index) {

// Insert cluster in TimeBin cluster array.
// Clusters are sorted according to Y coordinate.  

  if (fN==kMAX_CLUSTER_PER_TIME_BIN) {
    printf("AliTRDtimeBin::InsertCluster(): Too many clusters !\n"); 
    return;
  }
  if (fN==0) {fIndex[0]=index; fClusters[fN++]=c; return;}
  Int_t i=Find(c->GetY());
  memmove(fClusters+i+1 ,fClusters+i,(fN-i)*sizeof(AliTRDcluster*));
  memmove(fIndex   +i+1 ,fIndex   +i,(fN-i)*sizeof(UInt_t)); 
  fIndex[i]=index; fClusters[i]=c; fN++;
}  

//______________________________________________________

Int_t AliTRDtimeBin::Find(Double_t y) const {

// Returns index of the cluster nearest in Y    

  if (y <= fClusters[0]->GetY()) return 0;
  if (y > fClusters[fN-1]->GetY()) return fN;
  Int_t b=0, e=fN-1, m=(b+e)/2;
  for (; b<e; m=(b+e)/2) {
    if (y > fClusters[m]->GetY()) b=m+1;
    else e=m;
  }
  return m;
}    

//______________________________________________________
AliTRDcluster *AliTRDtimeBin::operator[](Int_t i)
{
  //
  // Index operator
  //
 
  return fClusters[i];

}
Commit	Line	Data
46d29e70	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	$Log$
	18	Revision 1.1.2.2 2000/10/04 16:34:58 cblume
	19	Replace include files by forward declarations
	20
	21	Revision 1.1.2.1 2000/09/22 14:47:52 cblume
	22	Add the tracking code
	23
	24	*/
	25
	26	#include "AliTRDcluster.h"
	27	#include "AliTRDconst.h"
	28	#include "AliTRDtimeBin.h"
	29
	30	ClassImp(AliTRDtimeBin)
	31
	32	//______________________________________________________
	33
	34	void AliTRDtimeBin::InsertCluster(AliTRDcluster* c, UInt_t index) {
	35
	36	// Insert cluster in TimeBin cluster array.
	37	// Clusters are sorted according to Y coordinate.
	38
	39	if (fN==kMAX_CLUSTER_PER_TIME_BIN) {
	40	printf("AliTRDtimeBin::InsertCluster(): Too many clusters !\n");
	41	return;
	42	}
	43	if (fN==0) {fIndex[0]=index; fClusters[fN++]=c; return;}
	44	Int_t i=Find(c->GetY());
	45	memmove(fClusters+i+1 ,fClusters+i,(fN-i)sizeof(AliTRDcluster));
	46	memmove(fIndex +i+1 ,fIndex +i,(fN-i)*sizeof(UInt_t));
	47	fIndex[i]=index; fClusters[i]=c; fN++;
	48	}
	49
	50	//______________________________________________________
	51
	52	Int_t AliTRDtimeBin::Find(Double_t y) const {
	53
	54	// Returns index of the cluster nearest in Y
	55
	56	if (y <= fClusters[0]->GetY()) return 0;
	57	if (y > fClusters[fN-1]->GetY()) return fN;
	58	Int_t b=0, e=fN-1, m=(b+e)/2;
	59	for (; b<e; m=(b+e)/2) {
	60	if (y > fClusters[m]->GetY()) b=m+1;
	61	else e=m;
	62	}
	63	return m;
	64	}
65
66	//______________________________________________________
67	AliTRDcluster *AliTRDtimeBin::operator[](Int_t i)
68	{
69	//
70	// Index operator
71	//
72
73	return fClusters[i];
74
75	}