dEtSum+=(*vPx)(iIn);
iIn++;
}
- for (Int_t iIn=0; iIn<fNin; iIn++) (*vPx)(iIn)=(*vPx)(iIn)/dEtSum;
+ for (iIn=0; iIn<fNin; iIn++) (*vPx)(iIn)=(*vPx)(iIn)/dEtSum;
Int_t njdim=2*fNclustMax+1;
mPyx->ResizeTo(fNin,njdim);
(*vPy)(0)=1;
TMatrixDColumn(*mPyx,0)=1;
Double_t ypos=0,xpos=0;
- for (Int_t iIn=0; iIn<fNin; iIn++){
+ for (iIn=0; iIn<fNin; iIn++){
(*mY)(0,0)+=(*vPx)(iIn)*xEta[iIn];
ypos+=(*vPx)(iIn)*TMath::Sin(xPhi[iIn]);
xpos+=(*vPx)(iIn)*TMath::Cos(xPhi[iIn]);
}
//-----------------------------------------------------------------------------------
-void AliDAJetFinder::DoubleClusters(Int_t nc,Int_t &nk,TVectorD *vPy,TMatrixD *mY)
+void AliDAJetFinder::DoubleClusters(Int_t nc,Int_t &nk, TVectorD *vPy, TMatrixD *mY) const
{
for(Int_t iClust=0; iClust<nc; iClust++){
(*vPy)(iClust)=(*vPy)(iClust)/2;
}
//-----------------------------------------------------------------------------------
-void AliDAJetFinder::Annealing(Int_t nk,Double_t **xData,TVectorD *vPx,TVectorD *vPy,TMatrixD *mPyx,TMatrixD *mY)
+void AliDAJetFinder::Annealing(Int_t nk,Double_t **xData, TVectorD *vPx, TVectorD *vPy, TMatrixD *mPyx, TMatrixD *mY)
{
// Main part of the algorithm
const Double_t pi=TMath::Pi();
delete y;
delete y1;
delete ry;
- delete [] m;
+ delete [] m;
}
//-----------------------------------------------------------------------------------
-void AliDAJetFinder::NumCl(Int_t &nc,Int_t &nk,TVectorD *vPy,TMatrixD *mPyx,TMatrixD *mY)
+void AliDAJetFinder::NumCl(Int_t &nc,Int_t &nk,TVectorD *vPy, TMatrixD *mPyx,TMatrixD *mY)
{
static Bool_t growcl=true;
}
//-----------------------------------------------------------------------------------
-void AliDAJetFinder::ReduceClusters(Int_t **iSame,Int_t nc,Int_t &ncout,Int_t **cont,Int_t *nSameOut)
+void AliDAJetFinder::ReduceClusters(Int_t **iSame,Int_t nc,Int_t &ncout,Int_t **cont,Int_t *nSameOut) const
{
Int_t *nSame = new Int_t[nc];
Int_t *iperm = new Int_t[nc];
}
//-----------------------------------------------------------------------------------
-void AliDAJetFinder::StoreJets(Int_t nk,Double_t **xData,Int_t *xx,TMatrixD *mY)
+void AliDAJetFinder::StoreJets(Int_t nk,Double_t **xData, Int_t *xx, TMatrixD *mY)
{
//evaluate significant clusters properties
const Double_t pi=TMath::Pi();