if ( q <= 4.0 * epsilon )
{
AliFatal("Error in calculating Dalitz mass histogram binning!");
- if ( idecay == 1 )
- fEPMass = 0;
- else
- fMPMass = 0;
}
kwHelp = (1.0 + q / (1.0 - delta)) * (1.0 + q / (1.0 - delta))
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)
Rot(Double_t pin[3], Double_t pout[3], Double_t costheta, Double_t sintheta,
Double_t cosphi, Double_t sinphi)
{
+// 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;