Adding protection against taking Sqrt of negative values...

[u/mrichter/AliRoot.git] / JETAN / AliDAJetHeader.cxx

//  **************************************************************************
//  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
//  *                                                                        *
//  * Author: The ALICE Off-line Project.                                    *
//  * Contributors are mentioned in the code where appropriate.              *
//  *                                                                        *
//  * Permission to use, copy, modify and distribute this software and its   *
//  * documentation strictly for non-commercial purposes is hereby granted   *
//  * without fee, provided that the above copyright notice appears in all   *
//  * copies and that both the copyright notice and this permission notice   *
//  * appear in the supporting documentation. The authors make no claims     *
//  * about the suitability of this software for any purpose. It is          *
//  * provided "as is" without express or implied warranty.                  *
//  **************************************************************************

//----------------------------------------------------------------------------
// Deterministic Annealing Jet header class
// Stores parameters of DA jet algorithm
// Author: Davide Perrino (davide.perrino@ba.infn.it, davide.perrino@cern.ch)
//----------------------------------------------------------------------------

#include "AliDAJetHeader.h"

ClassImp(AliDAJetHeader)


//---------------------------------------------------------------------
AliDAJetHeader::AliDAJetHeader():
	AliJetHeader("AliDAJetHeader"),
	fSelectJets(kTRUE),
	fRadius(0.7),
	fNclustMax(10),
	fFixedCl(kFALSE),
	fEtMin(10.)
{
    // Constructor
}

//---------------------------------------------------------------------
void AliDAJetHeader::SetRadius(Float_t radius)
{
    // The radius requested is used to estimate the number of clusters
    // to be found, in order to obtain jets with the expected area.
    // It must not be intended as a sharp limit on the cluster shape
    
  Float_t fEtaLim = 0.9;
  Int_t nclust = (Int_t) (4.*fEtaLim/(radius*radius)) + 1;
  SetNclust(nclust);
}