1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
21 #include "AliRICHRawCluster.h"
24 ClassImp(AliRICHRawCluster)
25 Int_t AliRICHRawCluster::Compare(TObject *obj)
28 // Compare two clusters
30 AliRICHRawCluster *raw=(AliRICHRawCluster *)obj;
34 else if (y<yo) return -1;
39 Int_t AliRICHRawCluster::
40 BinarySearch(Float_t y, TArrayF coord, Int_t from, Int_t upto)
42 // Find object using a binary search. Array must first have been sorted.
43 // Search can be limited by setting upto to desired index.
45 Int_t low=from, high=upto-1, half;
48 if(y>coord[half]) low=half;
54 void AliRICHRawCluster::SortMin(Int_t *idx,Float_t *xdarray,Float_t *xarray,Float_t *yarray,Float_t *qarray, Int_t ntr)
57 // Get the 3 closest points(cog) one can find on the second cathode
58 // starting from a given cog on first cathode
62 // Loop over deltax, only 3 times
67 Int_t id[3] = {-2,-2,-2};
68 Float_t jx[3] = {0.,0.,0.};
69 Float_t jy[3] = {0.,0.,0.};
70 Float_t jq[3] = {0.,0.,0.};
71 Int_t jid[3] = {-2,-2,-2};
81 if ((i == 1 && j == id[i-1])
82 ||(i == 2 && (j == id[i-1] || j == id[i-2]))) continue;
83 if (TMath::Abs(xdarray[j]) < xmin) {
84 xmin = TMath::Abs(xdarray[j]);
115 Int_t AliRICHRawCluster::PhysicsContribution()
118 // Type of physics processes
123 for (Int_t i=0; i<fMultiplicity; i++) {
124 if (fPhysicsMap[i]==2) iPhys++;
125 if (fPhysicsMap[i]==1) iMixed++;
126 if (fPhysicsMap[i]==0) iBg++;
128 if (iMixed==0 && iBg==0) {
130 } else if ((iPhys != 0 && iBg !=0) || iMixed != 0) {