TMatrixD *mPyx= new TMatrixD();
TMatrixD *mY = new TMatrixD();
InitDetAnn(dEtSum,xData,vPx,vPy,mPyx,mY);
- if (fNin < fNclustMax) return;
-
+ if (fNin < fNclustMax){
+ delete [] xData[0], delete [] xData[1];
+ delete vPx;
+ delete vPy;
+ delete mPyx;
+ delete mY;
+ return;
+ }
Int_t nc=1, nk;
DoubleClusters(nc,nk,vPy,mY);
do{ //loop over beta
}
TRandom2 r;
r.SetSeed(0);
- for (Int_t iIn=fNin; iIn<fNeff; iIn++){
+ for (iIn=fNin; iIn<fNeff; iIn++){
xEta[iIn]=r.Uniform(-1*etaEff,etaEff);
xPhi[iIn]=r.Uniform(0.,2*TMath::Pi());
(*vPx)(iIn)=r.Uniform(0.01,0.02);