AliOmegaDalitz::AliOmegaDalitz():
AliDecayer(),
fEPMass(0),
- fMPMass(0)
+ fMPMass(0),
+ fInit(0)
{
// Constructor
}
if ( kwHelp <= 0.0 )
{
AliFatal("Error in calculating Dalitz mass histogram binning!");
- if ( idecay == 1 )
- fEPMass = 0;
- else
- fMPMass = 0;
}
krollWada = (2.0 / mLL) * TMath::Exp(1.5 * TMath::Log(kwHelp))
* TMath::Sqrt(1.0 - 4.0 * epsilon / q)
//
//-----------------------------------------------------------------------------
+ if (!fInit) {
+ Init();
+ fInit=1;
+ }
+
Double_t pmass, lmass, omass, lpmass;
Double_t e1, p1, e3, p3;
Double_t betaSquare, lambda;
void AliOmegaDalitz::
Rot(Double_t pin[3], Double_t pout[3], Double_t costheta, Double_t sintheta,
- Double_t cosphi, Double_t sinphi)
+ Double_t cosphi, Double_t sinphi) const
{
+// Perform rotation
pout[0] = pin[0]*costheta*cosphi-pin[1]*sinphi+pin[2]*sintheta*cosphi;
pout[1] = pin[0]*costheta*sinphi+pin[1]*cosphi+pin[2]*sintheta*sinphi;
pout[2] = -1.0 * pin[0] * sintheta + pin[2] * costheta;
AliOmegaDalitz::AliOmegaDalitz(const AliOmegaDalitz &dalitz)
: AliDecayer(),
fEPMass(0),
- fMPMass(0)
+ fMPMass(0),
+ fInit(0)
{
// Copy constructor
dalitz.Copy(*this);