[u/mrichter/AliRoot.git] / HLTANALYSIS / JET / AliHLTJETBase.cxx

//-*- Mode: C++ -*-
//**************************************************************************
//* This file is property of and copyright by the ALICE HLT Project        * 
//* ALICE Experiment at CERN, All rights reserved.                         *
//*                                                                        *
//* Primary Authors: Jochen Thaeder <thaeder@kip.uni-heidelberg.de>        *
//*                  for The ALICE HLT Project.                            *
//*                                                                        *
//* 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.                  *
//**************************************************************************

/** @file   AliHLTJETBase.cxx
    @author Jochen Thaeder
    @date   
    @brief  Base functionality for HLT JET package
*/

// see header file for class documentation
// or
// refer to README to build package
// or
// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt

#include "AliHLTJETBase.h"

using namespace std;

/** ROOT macro for the implementation of ROOT specific class methods */
ClassImp(AliHLTJETBase)

/*
 * ---------------------------------------------------------------------------------
 *                            Constructor / Destructor
 * ---------------------------------------------------------------------------------
 */

//##################################################################################
  AliHLTJETBase::AliHLTJETBase() {
  // see header file for class documentation
  // or
  // refer to README to build package
  // or
  // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
}


//##################################################################################
AliHLTJETBase::~AliHLTJETBase() {
  // see header file for class documentation
}

/*
 * ---------------------------------------------------------------------------------
 *                              Initialize static const
 * ---------------------------------------------------------------------------------
 */

//##################################################################################
const Char_t *AliHLTJETBase::fgkJetAlgorithmType[] = { 
  "Anti Kt", 
  "Kt",
  "FFSC SquareCell",
  "FFSC RadiusCell"
};

#if 0

//################################################################################## 
Float_t AliHLTJETBase::GetDistance2( const Float_t eta1, const Float_t phi1, 
				     const Float_t eta2, const Float_t phi2) {
  // see header file for class documentation
  
  return ( (eta1-eta2)*(eta1-eta2) ) + ( (phi1-phi2)*(phi1-phi2) );
}

//##################################################################################
void AliHLTJETBase::XYZtoRPhiEta( const Double_t *xyz, Double_t *rpe ) {
  // see header file for class documentation

  // Spherical [r,phi,eta] - rpe
  // Cartesian [x,y,z] - xyz
  
  rpe[0] = TMath::Sqrt( xyz[0]*xyz[0] + xyz[1]*xyz[1] + xyz[2]*xyz[2] );
  rpe[1] = AliHLTJETBase::GetPhiFromXYZ( xyz );
  rpe[2] = 0.5 * TMath::Log( (rpe[0]+xyz[2])/(rpe[0]-xyz[2]) );
}

//##################################################################################
void AliHLTJETBase::XYZEtoRPhiEtaPt( const Double_t *xyze, Double_t *rpep ) {
  // see header file for class documentation

  // Spherical [r,phi,eta,pt] - rpep
  // Cartesian [x,y,z,E] - xyze
  
  rpep[0] = TMath::Sqrt( xyze[0]*xyze[0] + xyze[1]*xyze[1] + xyze[2]*xyze[2] );
  rpep[1] = AliHLTJETBase::GetPhiFromXYZ( xyze );
  rpep[2] = 0.5 * TMath::Log( (rpep[0]+xyze[2])/(rpep[0]-xyze[2]) );
  rpep[3] = AliHLTJETBase::GetPtFromXYZ( xyze );
}

//##################################################################################
void AliHLTJETBase::XYZEtoRPhiEtaPt( const Float_t *xyze, Double_t *rpep ) {
  // see header file for class documentation

  // Spherical [r,phi,eta,pt] - rpep
  // Cartesian [x,y,z,E] - xyze
  Double_t x = (Double_t) xyze[0];
  Double_t y = (Double_t) xyze[1];
  Double_t z = (Double_t) xyze[2];
  
  rpep[0] = TMath::Sqrt( x*x + y*y + z*z );
  rpep[1] = AliHLTJETBase::GetPhiFromXYZ( xyze );
  rpep[2] = 0.5 * TMath::Log( (rpep[0]+z)/(rpep[0]-z) );
  rpep[3] = AliHLTJETBase::GetPtFromXYZ( xyze );
}

//##################################################################################
void AliHLTJETBase::XYZEtoEPhiEtaPt( const Double_t *xyze, Double_t *epep ) {
  // see header file for class documentation

  // Spherical [r,phi,eta,pt] - rpep
  // Cartesian [x,y,z,E] - xyze

  Double_t r = TMath::Sqrt( xyze[0]*xyze[0] + xyze[1]*xyze[1] + xyze[2]*xyze[2] );
  
  epep[0] = xyze[3];
  epep[1] = AliHLTJETBase::GetPhiFromXYZ( xyze );
  epep[2] = 0.5 * TMath::Log( (r+xyze[2])/(r-xyze[2]) );
  epep[3] = AliHLTJETBase::GetPtFromXYZ( xyze );
}

//##################################################################################
void AliHLTJETBase::XYZEtoEPhiEtaPt( const Float_t *xyze, Double_t *epep ) {
  // see header file for class documentation

  // Spherical [r,phi,eta,pt] - rpep
  // Cartesian [x,y,z,E] - xyze
  Double_t x = (Double_t) xyze[0];
  Double_t y = (Double_t) xyze[1];
  Double_t z = (Double_t) xyze[2];
  Double_t e = (Double_t) xyze[3];
  Double_t r = TMath::Sqrt( x*x + y*y + z*z );
  
  epep[0] = e;
  epep[1] = AliHLTJETBase::GetPhiFromXYZ( xyze );
  epep[2] = 0.5 * TMath::Log( (r+z)/(r-z) );
  epep[3] = AliHLTJETBase::GetPtFromXYZ( xyze );
}

//##################################################################################
Double_t AliHLTJETBase::GetPtFromXYZ( const Double_t *pxpypz ) {
  // see header file for class documentation
  // Cartesian [px,py,pz] 
  
  return TMath::Sqrt( pxpypz[0]*pxpypz[0] + pxpypz[1]*pxpypz[1] );
}

//##################################################################################
Double_t AliHLTJETBase::GetPtFromXYZ( const Float_t *pxpypz ) {
  // see header file for class documentation
  // Cartesian [px,py,pz] 
  
  Double_t px = (Double_t) pxpypz[0];
  Double_t py = (Double_t) pxpypz[1];

  return TMath::Sqrt( px*px + py*py );
}

//##################################################################################
Double_t AliHLTJETBase::GetPhiFromXYZ( const Double_t *xyz ) {
  // see header file for class documentation
  // Cartesian [x,y,z,E] - xyz
 
  Double_t phi = TMath::ATan2( xyz[1], xyz[0] );  
  if ( phi < 0.0 ) phi += TMath::TwoPi();

  return phi;
}

//##################################################################################
Double_t AliHLTJETBase::GetPhiFromXYZ( const Float_t *xyz ) {
  // see header file for class documentation
  // Cartesian [x,y,z,E] - xyz
 
  Double_t x = (Double_t) xyz[0];
  Double_t y = (Double_t) xyz[1];

  Double_t phi = TMath::ATan2( y, x );  
  if ( phi < 0.0 ) phi += TMath::TwoPi();

  return phi;
}

//##################################################################################
Double_t AliHLTJETBase::GetEtaFromXYZ( const Double_t *xyz ) {
  // see header file for class documentation
  // Cartesian [x,y,z,E] - xyz
  
  Double_t r = TMath::Sqrt( xyz[0]*xyz[0] + xyz[1]*xyz[1] + xyz[2]*xyz[2] );
  return 0.5 * TMath::Log( (r+xyz[2])/(r-xyz[2]) );
}
//##################################################################################
Double_t AliHLTJETBase::GetEtaFromXYZ( const Float_t *xyz ) {
  // see header file for class documentation
  // Cartesian [x,y,z,E] - xyz
  
  Double_t x = (Double_t) xyz[0];
  Double_t y = (Double_t) xyz[1];
  Double_t z = (Double_t) xyz[2];
  
  Double_t r = TMath::Sqrt( x*x + y*y + z*z );
  return 0.5 * TMath::Log( (r+z)/(r-z) );
}


#endif
Commit	Line	Data
7c3c85cd	1	//-- Mode: C++ --
	2	//**************************************************************************
	3	//* This file is property of and copyright by the ALICE HLT Project *
	4	//* ALICE Experiment at CERN, All rights reserved. *
	5	//* *
	6	//* Primary Authors: Jochen Thaeder <thaeder@kip.uni-heidelberg.de> *
	7	//* for The ALICE HLT Project. *
	8	//* *
	9	//* Permission to use, copy, modify and distribute this software and its *
	10	//* documentation strictly for non-commercial purposes is hereby granted *
	11	//* without fee, provided that the above copyright notice appears in all *
	12	//* copies and that both the copyright notice and this permission notice *
	13	//* appear in the supporting documentation. The authors make no claims *
	14	//* about the suitability of this software for any purpose. It is *
	15	//* provided "as is" without express or implied warranty. *
	16	//**************************************************************************
	17
	18	/** @file AliHLTJETBase.cxx
	19	@author Jochen Thaeder
	20	@date
	21	@brief Base functionality for HLT JET package
	22	*/
	23
	24	// see header file for class documentation
	25	// or
	26	// refer to README to build package
	27	// or
	28	// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
	29
7c3c85cd	30	#include "AliHLTJETBase.h"
7c3c85cd	31
a7f38894	32	using namespace std;
a7f38894	33
7c3c85cd	34	/** ROOT macro for the implementation of ROOT specific class methods */
	35	ClassImp(AliHLTJETBase)
	36
	37	/*
	38	* ---------------------------------------------------------------------------------
	39	* Constructor / Destructor
	40	* ---------------------------------------------------------------------------------
	41	*/
	42
	43	//##################################################################################
	44	AliHLTJETBase::AliHLTJETBase() {
	45	// see header file for class documentation
	46	// or
	47	// refer to README to build package
	48	// or
	49	// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
	50	}
	51
	52
	53	//##################################################################################
	54	AliHLTJETBase::~AliHLTJETBase() {
	55	// see header file for class documentation
	56	}
1f9fec4a	57
fdc39952	58	/*
1f9fec4a	59	* ---------------------------------------------------------------------------------
	60	* Initialize static const
	61	* ---------------------------------------------------------------------------------
	62	*/
	63
	64	//##################################################################################
	65	const Char_t *AliHLTJETBase::fgkJetAlgorithmType[] = {
	66	"Anti Kt",
	67	"Kt",
	68	"FFSC SquareCell",
	69	"FFSC RadiusCell"
	70	};
	71
	72	#if 0
	73
fdc39952	74	//##################################################################################
	75	Float_t AliHLTJETBase::GetDistance2( const Float_t eta1, const Float_t phi1,
	76	const Float_t eta2, const Float_t phi2) {
	77	// see header file for class documentation
	78
	79	return ( (eta1-eta2)(eta1-eta2) ) + ( (phi1-phi2)(phi1-phi2) );
	80	}
7c3c85cd	81
	82	//##################################################################################
	83	void AliHLTJETBase::XYZtoRPhiEta( const Double_t xyz, Double_t rpe ) {
	84	// see header file for class documentation
	85
	86	// Spherical [r,phi,eta] - rpe
	87	// Cartesian [x,y,z] - xyz
	88
	89	rpe[0] = TMath::Sqrt( xyz[0]xyz[0] + xyz[1]xyz[1] + xyz[2]*xyz[2] );
	90	rpe[1] = AliHLTJETBase::GetPhiFromXYZ( xyz );
	91	rpe[2] = 0.5 * TMath::Log( (rpe[0]+xyz[2])/(rpe[0]-xyz[2]) );
	92	}
	93
	94	//##################################################################################
	95	void AliHLTJETBase::XYZEtoRPhiEtaPt( const Double_t xyze, Double_t rpep ) {
	96	// see header file for class documentation
	97
	98	// Spherical [r,phi,eta,pt] - rpep
	99	// Cartesian [x,y,z,E] - xyze
	100
	101	rpep[0] = TMath::Sqrt( xyze[0]xyze[0] + xyze[1]xyze[1] + xyze[2]*xyze[2] );
	102	rpep[1] = AliHLTJETBase::GetPhiFromXYZ( xyze );
	103	rpep[2] = 0.5 * TMath::Log( (rpep[0]+xyze[2])/(rpep[0]-xyze[2]) );
	104	rpep[3] = AliHLTJETBase::GetPtFromXYZ( xyze );
	105	}
	106
	107	//##################################################################################
	108	void AliHLTJETBase::XYZEtoRPhiEtaPt( const Float_t xyze, Double_t rpep ) {
	109	// see header file for class documentation
	110
	111	// Spherical [r,phi,eta,pt] - rpep
	112	// Cartesian [x,y,z,E] - xyze
	113	Double_t x = (Double_t) xyze[0];
	114	Double_t y = (Double_t) xyze[1];
	115	Double_t z = (Double_t) xyze[2];
	116
	117	rpep[0] = TMath::Sqrt( xx + yy + z*z );
	118	rpep[1] = AliHLTJETBase::GetPhiFromXYZ( xyze );
	119	rpep[2] = 0.5 * TMath::Log( (rpep[0]+z)/(rpep[0]-z) );
	120	rpep[3] = AliHLTJETBase::GetPtFromXYZ( xyze );
	121	}
	122
	123	//##################################################################################
	124	void AliHLTJETBase::XYZEtoEPhiEtaPt( const Double_t xyze, Double_t epep ) {
	125	// see header file for class documentation
	126
	127	// Spherical [r,phi,eta,pt] - rpep
	128	// Cartesian [x,y,z,E] - xyze
	129
	130	Double_t r = TMath::Sqrt( xyze[0]xyze[0] + xyze[1]xyze[1] + xyze[2]*xyze[2] );
	131
	132	epep[0] = xyze[3];
	133	epep[1] = AliHLTJETBase::GetPhiFromXYZ( xyze );
	134	epep[2] = 0.5 * TMath::Log( (r+xyze[2])/(r-xyze[2]) );
	135	epep[3] = AliHLTJETBase::GetPtFromXYZ( xyze );
	136	}
	137
	138	//##################################################################################
	139	void AliHLTJETBase::XYZEtoEPhiEtaPt( const Float_t xyze, Double_t epep ) {
	140	// see header file for class documentation
	141
	142	// Spherical [r,phi,eta,pt] - rpep
	143	// Cartesian [x,y,z,E] - xyze
	144	Double_t x = (Double_t) xyze[0];
145	Double_t y = (Double_t) xyze[1];
146	Double_t z = (Double_t) xyze[2];
147	Double_t e = (Double_t) xyze[3];
148	Double_t r = TMath::Sqrt( xx + yy + z*z );
149
150	epep[0] = e;
151	epep[1] = AliHLTJETBase::GetPhiFromXYZ( xyze );
152	epep[2] = 0.5 * TMath::Log( (r+z)/(r-z) );
153	epep[3] = AliHLTJETBase::GetPtFromXYZ( xyze );
154	}
155
156	//##################################################################################
157	Double_t AliHLTJETBase::GetPtFromXYZ( const Double_t *pxpypz ) {
158	// see header file for class documentation
159	// Cartesian [px,py,pz]
160
161	return TMath::Sqrt( pxpypz[0]pxpypz[0] + pxpypz[1]pxpypz[1] );
162	}
163
164	//##################################################################################
165	Double_t AliHLTJETBase::GetPtFromXYZ( const Float_t *pxpypz ) {
166	// see header file for class documentation
167	// Cartesian [px,py,pz]
168
169	Double_t px = (Double_t) pxpypz[0];
170	Double_t py = (Double_t) pxpypz[1];
171
172	return TMath::Sqrt( pxpx + pypy );
173	}
174
175	//##################################################################################
176	Double_t AliHLTJETBase::GetPhiFromXYZ( const Double_t *xyz ) {
177	// see header file for class documentation
178	// Cartesian [x,y,z,E] - xyz
179
180	Double_t phi = TMath::ATan2( xyz[1], xyz[0] );
181	if ( phi < 0.0 ) phi += TMath::TwoPi();
182
183	return phi;
184	}
185
186	//##################################################################################
187	Double_t AliHLTJETBase::GetPhiFromXYZ( const Float_t *xyz ) {
188	// see header file for class documentation
189	// Cartesian [x,y,z,E] - xyz
190
191	Double_t x = (Double_t) xyz[0];
192	Double_t y = (Double_t) xyz[1];
193
194	Double_t phi = TMath::ATan2( y, x );
195	if ( phi < 0.0 ) phi += TMath::TwoPi();
196
197	return phi;
198	}
199
200	//##################################################################################
201	Double_t AliHLTJETBase::GetEtaFromXYZ( const Double_t *xyz ) {
202	// see header file for class documentation
203	// Cartesian [x,y,z,E] - xyz
204
205	Double_t r = TMath::Sqrt( xyz[0]xyz[0] + xyz[1]xyz[1] + xyz[2]*xyz[2] );
206	return 0.5 * TMath::Log( (r+xyz[2])/(r-xyz[2]) );
207	}
208	//##################################################################################
209	Double_t AliHLTJETBase::GetEtaFromXYZ( const Float_t *xyz ) {
210	// see header file for class documentation
211	// Cartesian [x,y,z,E] - xyz
212
213	Double_t x = (Double_t) xyz[0];
214	Double_t y = (Double_t) xyz[1];
215	Double_t z = (Double_t) xyz[2];
216
217	Double_t r = TMath::Sqrt( xx + yy + z*z );
218	return 0.5 * TMath::Log( (r+z)/(r-z) );
219	}
220
7c3c85cd	221
	222
	223	#endif