[u/mrichter/AliRoot.git] / TEvtGen / EvtGenBase / EvtStringParticle.cxx

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of the EvtGen package developed jointly
//      for the BaBar and CLEO collaborations.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
//      Copyright (C) 1998      Caltech, UCSB
//
// Module: EvtStringParticle.cc
//
// Description: Class to describe the partons that are produced in JetSet.
//
// Modification history:
//
//    RYD       Febuary 27,1998       Module created
//
//------------------------------------------------------------------------
// 
#include "EvtGenBase/EvtPatches.hh"
#include <iostream>
#include <math.h>
#include <stdlib.h>
#include "EvtGenBase/EvtStringParticle.hh"
#include "EvtGenBase/EvtVector4R.hh"
#include "EvtGenBase/EvtReport.hh"


EvtStringParticle::~EvtStringParticle(){

  if (_npartons!=0){
    
    delete [] _p4partons;
    delete [] _idpartons;
    
  }

}


EvtStringParticle::EvtStringParticle(){

  _p4partons=0;
  _idpartons=0;
  _npartons=0;

  return;
  
}

void EvtStringParticle::init(EvtId id, const EvtVector4R& p4){

  _validP4=true;
  setp(p4);
  setpart_num(id);

}


void EvtStringParticle::initPartons(int npartons,
			      EvtVector4R* p4partons,EvtId* idpartons){

  _p4partons = new EvtVector4R[npartons];
  _idpartons = new EvtId[npartons];

  int i;

  _npartons=npartons;

  for(i=0;i<npartons;i++){

    _p4partons[i]=p4partons[i];
    _idpartons[i]=idpartons[i];

  }

}

int EvtStringParticle::getNPartons(){

  return _npartons;

}

EvtId EvtStringParticle::getIdParton(int i){

  return _idpartons[i];

}

EvtVector4R EvtStringParticle::getP4Parton(int i){

  return _p4partons[i];

}
  
EvtSpinDensity EvtStringParticle::rotateToHelicityBasis() const{

  report(ERROR,"EvtGen") << "rotateToHelicityBasis not implemented for strin particle."; 
  report(ERROR,"EvtGen") << "Will terminate execution."; 

  ::abort();

  EvtSpinDensity rho;
  return rho;

}

EvtSpinDensity EvtStringParticle::rotateToHelicityBasis(double,
							double,
							double) const{

  report(ERROR,"EvtGen") << "rotateToHelicityBasis(alpha,beta,gamma) not implemented for string particle."; 
  report(ERROR,"EvtGen") << "Will terminate execution."; 

  ::abort();

  EvtSpinDensity rho;
  return rho;

}
Commit	Line	Data
da0e9ce3	1	//--------------------------------------------------------------------------
	2	//
	3	// Environment:
	4	// This software is part of the EvtGen package developed jointly
	5	// for the BaBar and CLEO collaborations. If you use all or part
	6	// of it, please give an appropriate acknowledgement.
	7	//
	8	// Copyright Information: See EvtGen/COPYRIGHT
	9	// Copyright (C) 1998 Caltech, UCSB
	10	//
	11	// Module: EvtStringParticle.cc
	12	//
	13	// Description: Class to describe the partons that are produced in JetSet.
	14	//
	15	// Modification history:
	16	//
	17	// RYD Febuary 27,1998 Module created
	18	//
	19	//------------------------------------------------------------------------
	20	//
	21	#include "EvtGenBase/EvtPatches.hh"
	22	#include <iostream>
	23	#include <math.h>
	24	#include <stdlib.h>
	25	#include "EvtGenBase/EvtStringParticle.hh"
	26	#include "EvtGenBase/EvtVector4R.hh"
	27	#include "EvtGenBase/EvtReport.hh"
	28
	29
	30	EvtStringParticle::~EvtStringParticle(){
	31
	32	if (_npartons!=0){
	33
	34	delete [] _p4partons;
	35	delete [] _idpartons;
	36
	37	}
	38
	39	}
	40
	41
	42	EvtStringParticle::EvtStringParticle(){
	43
	44	_p4partons=0;
	45	_idpartons=0;
	46	_npartons=0;
	47
	48	return;
	49
	50	}
	51
	52	void EvtStringParticle::init(EvtId id, const EvtVector4R& p4){
	53
	54	_validP4=true;
	55	setp(p4);
	56	setpart_num(id);
	57
	58	}
	59
	60
	61	void EvtStringParticle::initPartons(int npartons,
	62	EvtVector4R* p4partons,EvtId* idpartons){
	63
	64	_p4partons = new EvtVector4R[npartons];
65	_idpartons = new EvtId[npartons];
66
67	int i;
68
69	_npartons=npartons;
70
71	for(i=0;i<npartons;i++){
72
73	_p4partons[i]=p4partons[i];
74	_idpartons[i]=idpartons[i];
75
76	}
77
78	}
79
80	int EvtStringParticle::getNPartons(){
81
82	return _npartons;
83
84	}
85
86	EvtId EvtStringParticle::getIdParton(int i){
87
88	return _idpartons[i];
89
90	}
91
92	EvtVector4R EvtStringParticle::getP4Parton(int i){
93
94	return _p4partons[i];
95
96	}
97
98	EvtSpinDensity EvtStringParticle::rotateToHelicityBasis() const{
99
100	report(ERROR,"EvtGen") << "rotateToHelicityBasis not implemented for strin particle.";
101	report(ERROR,"EvtGen") << "Will terminate execution.";
102
103	::abort();
104
105	EvtSpinDensity rho;
106	return rho;
107
108	}
109
110	EvtSpinDensity EvtStringParticle::rotateToHelicityBasis(double,
111	double,
112	double) const{
113
114	report(ERROR,"EvtGen") << "rotateToHelicityBasis(alpha,beta,gamma) not implemented for string particle.";
115	report(ERROR,"EvtGen") << "Will terminate execution.";
116
117	::abort();
118
119	EvtSpinDensity rho;
120	return rho;
121
122	}
123
124