[u/mrichter/AliRoot.git] / ITS / AliITSRawClusterSPD.cxx

/**************************************************************************
 * Copyright(c) 2000-2004, 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.                  *
 **************************************************************************/
#include <Riostream.h>
 
#include "AliITSRawClusterSPD.h"

////////////////////////////////////////////////////
//  Cluster classes for set:ITS                   //
//  Raw Clusters for SPD                          //
//                                                //
////////////////////////////////////////////////////

ClassImp(AliITSRawClusterSPD)

//______________________________________________________________________
AliITSRawClusterSPD::AliITSRawClusterSPD(){
  // Default constructor
  fX=fZ=fQ=0.0;
  fZStart=fZStop=0.0;
  fNClZ=fNClX=fXStart=fXStop=fZend=fNTracks=0;
  fTracks[0]=fTracks[1]=fTracks[2]=-3,fModule=0;
}

//______________________________________________________________________
AliITSRawClusterSPD::AliITSRawClusterSPD(Float_t clz,Float_t clx,
					 Float_t Charge,Int_t ClusterSizeZ,
					 Int_t ClusterSizeX,Int_t xstart,
					 Int_t xstop,
					 Float_t zstart,Float_t zstop,
					 Int_t zend,Int_t module) {
    // constructor

    fZ       = clz;
    fX       = clx;
    fQ       = Charge;
    fNClZ    = ClusterSizeZ;
    fNClX    = ClusterSizeX;
    fXStart  = xstart;
    fXStop   = xstop;
    fZStart  = zstart;
    fZStop   = zstop;
    fZend    = zend;
    fModule  = module;
}
//______________________________________________________________________
void AliITSRawClusterSPD::Add(AliITSRawClusterSPD* clJ) {
    // Recolculate the new center of gravity coordinate and cluster sizes
    // in both directions after grouping of clusters

    if(this->fZStop < clJ->ZStop()) this->fZStop = clJ->ZStop();
    this->fZ      = this->fZ + clJ->Z();
    this->fX      = (this->fX + clJ->X())/2.;
    this->fQ      = this->fQ + clJ->Q();
    this->fXStart = clJ->XStart(); // for a comparison with the next
    this->fXStop  = clJ->XStop();  // z column
    if(this->fZend < clJ->Zend())       this->fZend    = clJ->Zend();
    this->fNClX   = this->fXStop - this->fXStart + 1; 
    (this->fNClZ)++;

    return;
}
//______________________________________________________________________
Bool_t AliITSRawClusterSPD::Brother(AliITSRawClusterSPD* cluster,
				    Float_t dz,Float_t dx) const {
    // fXStart, fXstop and fZend information is used now instead of dz and dx
    // to check an absent (or a present) of the gap between two pixels in 
    // both x and z directions. The increasing order of fZend is used.
    Bool_t brother = kFALSE;  
    Bool_t test2 = kFALSE;  
    Bool_t test3 = kFALSE;

    dx = dz = 0; // to remove unused variable warning.
    // Diagonal clusters are included:
    if(fXStop >= (cluster->XStart() -1) && 
       fXStart <= (cluster->XStop()+1)) test2 = kTRUE;

    // Diagonal clusters are excluded:   
    // if(fXStop >= cluster->XStart() &&
    //    fXStart <= cluster->XStop()) test2 = kTRUE;
    if(cluster->Zend() == (fZend + 1)) test3 = kTRUE; 
    if(test2 && test3) {
	// cout<<"test 2,3 0k, brother = true "<<endl;
	return brother = kTRUE;
    } // end if
    return brother;
}
//______________________________________________________________________
void AliITSRawClusterSPD::PrintInfo() const{
    // print

    cout << ", Z: " << fZ << ", X: " << fX << ", Charge: " << fQ<<endl;
    cout << " Z cluster size: " << fNClZ <<", X cluster size "<< fNClX <<endl;
    cout <<" XStart, XStop,Zend, Module ="<<fXStart<<","
	 <<fXStop<<","<<fZend << "," << fModule<<endl;
}

//______________________________________________________________________
void AliITSRawClusterSPD::SetTracks(Int_t track0, Int_t track1, Int_t track2) {
  // set tracks in cluster (no more than three ones)
  fTracks[0]=track0;
  fTracks[1]=track1;
  fTracks[2]=track2;
}
Commit	Line	Data
41b19549	1	/**************************************************************************
	2	* Copyright(c) 2000-2004, 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	#include <Riostream.h>
	16
	17	#include "AliITSRawClusterSPD.h"
	18
	19	////////////////////////////////////////////////////
	20	// Cluster classes for set:ITS //
	21	// Raw Clusters for SPD //
	22	// //
	23	////////////////////////////////////////////////////
	24
	25	ClassImp(AliITSRawClusterSPD)
	26
	27	//______________________________________________________________________
	28	AliITSRawClusterSPD::AliITSRawClusterSPD(){
	29	// Default constructor
	30	fX=fZ=fQ=0.0;
	31	fZStart=fZStop=0.0;
	32	fNClZ=fNClX=fXStart=fXStop=fZend=fNTracks=0;
	33	fTracks[0]=fTracks[1]=fTracks[2]=-3,fModule=0;
	34	}
	35
	36	//______________________________________________________________________
	37	AliITSRawClusterSPD::AliITSRawClusterSPD(Float_t clz,Float_t clx,
	38	Float_t Charge,Int_t ClusterSizeZ,
	39	Int_t ClusterSizeX,Int_t xstart,
	40	Int_t xstop,
	41	Float_t zstart,Float_t zstop,
	42	Int_t zend,Int_t module) {
	43	// constructor
	44
	45	fZ = clz;
	46	fX = clx;
	47	fQ = Charge;
	48	fNClZ = ClusterSizeZ;
	49	fNClX = ClusterSizeX;
	50	fXStart = xstart;
	51	fXStop = xstop;
	52	fZStart = zstart;
	53	fZStop = zstop;
	54	fZend = zend;
	55	fModule = module;
	56	}
	57	//______________________________________________________________________
	58	void AliITSRawClusterSPD::Add(AliITSRawClusterSPD* clJ) {
	59	// Recolculate the new center of gravity coordinate and cluster sizes
	60	// in both directions after grouping of clusters
	61
	62	if(this->fZStop < clJ->ZStop()) this->fZStop = clJ->ZStop();
	63	this->fZ = this->fZ + clJ->Z();
	64	this->fX = (this->fX + clJ->X())/2.;
65	this->fQ = this->fQ + clJ->Q();
66	this->fXStart = clJ->XStart(); // for a comparison with the next
67	this->fXStop = clJ->XStop(); // z column
68	if(this->fZend < clJ->Zend()) this->fZend = clJ->Zend();
69	this->fNClX = this->fXStop - this->fXStart + 1;
70	(this->fNClZ)++;
71
72	return;
73	}
74	//______________________________________________________________________
75	Bool_t AliITSRawClusterSPD::Brother(AliITSRawClusterSPD* cluster,
76	Float_t dz,Float_t dx) const {
77	// fXStart, fXstop and fZend information is used now instead of dz and dx
78	// to check an absent (or a present) of the gap between two pixels in
79	// both x and z directions. The increasing order of fZend is used.
80	Bool_t brother = kFALSE;
81	Bool_t test2 = kFALSE;
82	Bool_t test3 = kFALSE;
83
84	dx = dz = 0; // to remove unused variable warning.
85	// Diagonal clusters are included:
86	if(fXStop >= (cluster->XStart() -1) &&
87	fXStart <= (cluster->XStop()+1)) test2 = kTRUE;
88
89	// Diagonal clusters are excluded:
90	// if(fXStop >= cluster->XStart() &&
91	// fXStart <= cluster->XStop()) test2 = kTRUE;
92	if(cluster->Zend() == (fZend + 1)) test3 = kTRUE;
93	if(test2 && test3) {
94	// cout<<"test 2,3 0k, brother = true "<<endl;
95	return brother = kTRUE;
96	} // end if
97	return brother;
98	}
99	//______________________________________________________________________
100	void AliITSRawClusterSPD::PrintInfo() const{
101	// print
102
103	cout << ", Z: " << fZ << ", X: " << fX << ", Charge: " << fQ<<endl;
104	cout << " Z cluster size: " << fNClZ <<", X cluster size "<< fNClX <<endl;
105	cout <<" XStart, XStop,Zend, Module ="<<fXStart<<","
106	<<fXStop<<","<<fZend << "," << fModule<<endl;
107	}
108
109	//______________________________________________________________________
110	void AliITSRawClusterSPD::SetTracks(Int_t track0, Int_t track1, Int_t track2) {
111	// set tracks in cluster (no more than three ones)
112	fTracks[0]=track0;
113	fTracks[1]=track1;
114	fTracks[2]=track2;
115	}