Coding rule violations corrected
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 6 Mar 2011 14:24:54 +0000 (14:24 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 6 Mar 2011 14:24:54 +0000 (14:24 +0000)
33 files changed:
JETAN/AliAnalysisTaskJetCluster.cxx
JETAN/AliAnalysisTaskJets.cxx
JETAN/AliDAJetFinder.cxx
JETAN/AliDAJetFinder.h
JETAN/AliEventShape.cxx
JETAN/AliFastJetHeader.h
JETAN/AliFastJetHeaderV1.h
JETAN/AliFastJetInput.cxx
JETAN/AliFastJetInput.h
JETAN/AliJetAODFillUnitArrayEMCalDigits.cxx
JETAN/AliJetAODReaderHeader.h
JETAN/AliJetBkg.cxx
JETAN/AliJetBkg.h
JETAN/AliJetDummyGeo.cxx
JETAN/AliJetDummyGeo.h
JETAN/AliJetDummyShishKebabTrd1Module.h
JETAN/AliJetFinder.cxx
JETAN/AliJetFinder.h
JETAN/AliJetHadronCorrectionv0.cxx
JETAN/AliJetHadronCorrectionv0.h
JETAN/AliJetHadronCorrectionv1.cxx
JETAN/AliJetHadronCorrectionv1.h
JETAN/AliJetHistos.cxx
JETAN/AliJetKineReaderHeader.h
JETAN/AliJetProductionData.h
JETAN/AliJetReaderHeader.h
JETAN/AliJetSelector.h
JETAN/AliJetUnitArray.cxx
JETAN/AliJetUnitArray.h
JETAN/AliKMeansClustering.cxx
JETAN/AliKMeansClustering.h
JETAN/AliUA1JetFinderV2.cxx
JETAN/AliUA1JetFinderV2.h

index 8b35028..2265781 100644 (file)
@@ -45,6 +45,7 @@
 #include "AliESDEvent.h"
 #include "AliAODEvent.h"
 #include "AliAODHandler.h"
+#include "AliAODExtension.h"
 #include "AliAODTrack.h"
 #include "AliAODJet.h"
 #include "AliAODMCParticle.h"
index f15db37..bf4cbbc 100644 (file)
@@ -35,6 +35,7 @@
 #include "AliAODEvent.h"
 #include "AliAODJetEventBackground.h"
 #include "AliAODHandler.h"
+#include "AliAODExtension.h"
 #include "AliMCEventHandler.h"
 #include "AliESDInputHandler.h"
 #include "AliMCEvent.h"
index ee945e8..e27ee30 100644 (file)
@@ -156,8 +156,9 @@ void AliDAJetFinder::InitDetAnn(Double_t &dEtSum,Double_t **xData,TVectorD *vPx,
 }
 
 //-----------------------------------------------------------------------------------
-void AliDAJetFinder::DoubleClusters(Int_t nc,Int_t &nk,  TVectorD *vPy,  TMatrixD *mY) const
+void AliDAJetFinder::DoubleClusters(Int_t nc,Int_t &nk,  TVectorD *vPy, TMatrixD *mY) const
 {
+// Return double clusters
        for(Int_t iClust=0; iClust<nc; iClust++){
                (*vPy)(iClust)=(*vPy)(iClust)/2;
                (*vPy)(nc+iClust)=(*vPy)(iClust);
@@ -264,6 +265,7 @@ void AliDAJetFinder::Annealing(Int_t nk,Double_t **xData,  TVectorD *vPx,  TVect
 //-----------------------------------------------------------------------------------
 void AliDAJetFinder::NumCl(Int_t &nc,Int_t &nk,TVectorD *vPy,  TMatrixD *mPyx,TMatrixD *mY)
 {
+    // Number of clusters
        static Bool_t growcl=true;
        
        if (nk==2) growcl=true;
@@ -321,6 +323,7 @@ void AliDAJetFinder::NumCl(Int_t &nc,Int_t &nk,TVectorD *vPy,  TMatrixD *mPyx,TM
 //-----------------------------------------------------------------------------------
 void AliDAJetFinder::ReduceClusters(Int_t **iSame,Int_t nc,Int_t &ncout,Int_t **cont,Int_t *nSameOut) const
 {
+// Reduction step
        Int_t *nSame = new Int_t[nc];
        Int_t *iperm = new Int_t[nc];
        Int_t *go = new Int_t[nc];
index b9b1f9a..4e5a0c6 100644 (file)
@@ -26,7 +26,7 @@ public:
     void Annealing     (Int_t nk, Double_t **xData, TVectorD *vPx, TVectorD *vPy, TMatrixD *mPyx, TMatrixD *mY);
     void NumCl         (Int_t &nc, Int_t &nk, TVectorD *vPy, TMatrixD *mPyx, TMatrixD *mY);
     void ReduceClusters(Int_t **iSame, Int_t nc, Int_t &ncout, Int_t **cont, Int_t *nSameOut) const;
-    void DoubleClusters(Int_t nc, Int_t &nk, TVectorD *vPy, TMatrixD *mY) const;
+    void DoubleClusters(Int_t nc, Int_t &nk,  TVectorD *vPy,  TMatrixD *mY) const;
     void EndDetAnn     (Int_t &nk, Double_t **xData, Int_t *xx, Double_t etx, TVectorD *px, TVectorD *py, TMatrixD *pyx, TMatrixD *y);
     void StoreJets     (Int_t nk, Double_t **xData, Int_t *xx, TMatrixD *mY);
 
index 5864a47..e34640d 100644 (file)
@@ -259,7 +259,7 @@ Double_t AliEventShape::GetCircularityMC(AliMCEvent* mcEvent, Int_t  nstudymin,
 
 
   if (nmctracks < nstudymin) {
-      Printf("Too few particles, abort");
+      Printf("Too few particles, stopping");
       return -2;
   }
 
index 224e891..14f9c8b 100644 (file)
@@ -67,12 +67,12 @@ class AliFastJetHeader : public AliJetHeader
   fastjet::RecombinationScheme fRecombScheme; // = fastjet::BIpt_scheme;
   
   //fastjet::GhostedAreaSpec parameters
-  Double_t fGhostEtaMax;
-  Double_t fGhostArea;
-  Int_t    fActiveAreaRepeats;
+  Double_t fGhostEtaMax;       // max area of ghosts
+  Double_t fGhostArea;         // ghost area
+  Int_t    fActiveAreaRepeats; // number of repetitions of active area 
   
   //fastjet::AreaDefinition parameters
-  fastjet::AreaType fAreaType; 
+  fastjet::AreaType fAreaType; // the are type
   
   //fastjet::ClusterSequenceArea options parameters
   Double_t fPtMin; //jets with pt > ptmin
index dab3825..d8f7859 100644 (file)
@@ -69,19 +69,19 @@ class AliFastJetHeaderV1 : public AliJetHeader
  protected:
 
   //fastjet::JetDefinition parameters
-  Double_t fRparam;
+  Double_t fRparam;   // R param
   Double_t fRparamBkg;//R param for bkg calculation
   fastjet::JetAlgorithm fAlgorithm; //fastjet::kt_algorithm
   fastjet::Strategy fStrategy;  //= fastjet::Best;
   fastjet::RecombinationScheme fRecombScheme; // = fastjet::BIpt_scheme;
   
   //fastjet::GhostedAreaSpec parameters
-  Double_t fGhostEtaMax;
-  Double_t fGhostArea;
-  Int_t    fActiveAreaRepeats;
+  Double_t fGhostEtaMax;       // Max eta for ghosts
+  Double_t fGhostArea;         // Ghost area 
+  Int_t    fActiveAreaRepeats; // Active are repetitions
   
   //fastjet::AreaDefinition parameters
-  fastjet::AreaType fAreaType; 
+  fastjet::AreaType fAreaType; // area types
   
   //fastjet::ClusterSequenceArea options parameters
   Double_t fPtMin; //jets with pt > ptmin
index d075d11..6814f98 100644 (file)
@@ -189,7 +189,7 @@ Float_t  AliFastJetInput::EtaToTheta(Float_t arg)
 
 
 }
-Double_t AliFastJetInput::Thermalspectrum(Double_t *x, Double_t *par){
+Double_t AliFastJetInput::Thermalspectrum(const Double_t *x, const Double_t *par){
 
   return x[0]*TMath::Exp(-x[0]/par[0]);
 
index 17d70af..4049b0c 100644 (file)
@@ -42,9 +42,8 @@ class AliFastJetInput : public TObject
     void FillInput();
     vector<fastjet::PseudoJet> GetInputParticles()   const {return fInputParticles;}
     vector<fastjet::PseudoJet> GetInputParticlesCh() const {return fInputParticlesCh;}
-
     Float_t  EtaToTheta(Float_t arg);
-    static Double_t Thermalspectrum(Double_t *x, Double_t *par);
+    static Double_t Thermalspectrum(const Double_t *x, const Double_t *par);
 
  private:
    AliJetReader *fReader;  //! reader 
index 869d7c0..43fb930 100755 (executable)
@@ -386,7 +386,7 @@ void AliJetAODFillUnitArrayEMCalDigits::Exec(Option_t* const /*option*/)
 //____________________________________________________________________________
 void AliJetAODFillUnitArrayEMCalDigits::GetTracksPointingToCell(TArrayS* array,Double_t eta, Double_t phi, Double_t cut)
 {
+// Get all tracks pointing to cell 
   int size=0;
   
   for (Int_t itrk =  0; itrk <  fAOD->GetNumberOfTracks() ; itrk++) { //track loop
index 169d089..60cc637 100644 (file)
@@ -22,7 +22,7 @@ class AliJetAODReaderHeader : public AliJetReaderHeader
 
   // Getters
   Int_t   GetNaod()       const {return fNaod;}
-  UInt_t  GetTestFilterMask(){return fTestFilterMask;}     
+  UInt_t  GetTestFilterMask()   const {return fTestFilterMask;}            
   TString GetNonStdBranch()     const { return fNonStdBranch; }
   Bool_t  GetReadStdBranch()    const { return (fReadBranches & kReadStdBranch); }
   Bool_t  GetReadNonStdBranch() const { return (fReadBranches & kReadNonStdBranch); }
index 69185cb..ad3c944 100644 (file)
@@ -291,7 +291,7 @@ meanareab, false);
 
 //___________________________________________________________________
 Float_t AliJetBkg::BkgFastJet(){
-  
+// Return background  
   AliFastJetHeaderV1 *header = (AliFastJetHeaderV1*)fHeader;
   Bool_t debug  = header->GetDebug();     // debug option
 
@@ -313,7 +313,7 @@ Float_t AliJetBkg::BkgFastJet(){
 }
 //___________________________________________________________________
 Float_t  AliJetBkg::BkgChargedFastJet(){
-
+// Background for charged jets
   AliFastJetHeaderV1 *header = (AliFastJetHeaderV1*)fHeader;
   Bool_t debug  = header->GetDebug();     // debug option
 
@@ -593,8 +593,8 @@ Double_t  AliJetBkg::BkgFunction(Double_t */*x*/,Double_t */*par*/){
 }
 
 
-Bool_t AliJetBkg::EmcalAcceptance(const Float_t eta, const Float_t phi, const Float_t radius){
+Bool_t AliJetBkg::EmcalAcceptance(const Float_t eta, const Float_t phi, const Float_t radius) const{
+// Apply emcal acceptance cuts
   Float_t meanPhi=190./180.*TMath::Pi()-110./180.*TMath::Pi()/2;
   Float_t deltaphi=110./180.*TMath::Pi();
   Float_t phicut=deltaphi/2.-radius;
index b901e37..3764ace 100644 (file)
@@ -40,7 +40,7 @@ class AliJetBkg : public TObject
 //    Float_t BkgRemoveJetLeading(TClonesArray* fAODJets);
     Float_t BkgRemoveJetLeadingFromUArray(TClonesArray* fAODJets);
     Float_t EtaToTheta(Float_t arg);
-    Bool_t EmcalAcceptance(const Float_t eta, const Float_t phi, const Float_t radius);
+    Bool_t EmcalAcceptance(const Float_t eta, const Float_t phi, const Float_t radius) const;
     static Double_t BkgFunction(Double_t *x,Double_t *par);
     
  private:
index 052cb22..2892cbc 100644 (file)
@@ -280,7 +280,7 @@ AliJetDummyGeo::~AliJetDummyGeo()
 }
 
 //------------------------------------------------------------------------------------
-void AliJetDummyGeo::EtaPhiFromIndex(Int_t absId, Float_t& eta, Float_t& phi)
+void AliJetDummyGeo::EtaPhiFromIndex(Int_t absId, Float_t& eta, Float_t& phi) const
 {
   // Nov 16, 2006- float to double
   // version for TRD1 only
index ab5660d..418846a 100644 (file)
@@ -53,7 +53,7 @@ class AliJetDummyGeo : public TObject
   Float_t GetArm1PhiMin()           const  {return fArm1PhiMin;}
   Float_t GetArm1PhiMax()           const  {return fArm1PhiMax;}
   Float_t GetIPDistance()           const  {return fIPDistance;} 
-  void    EtaPhiFromIndex(Int_t id, Float_t& eta, Float_t& phi);
+  void    EtaPhiFromIndex(Int_t id, Float_t& eta, Float_t& phi) const;
   void    GetGlobal(const Double_t *loc, Double_t *glob, Int_t ind) const;
   void    GetGlobal(Int_t absId, Double_t glob[3]) const;
   void    GetGlobal(Int_t absId, TVector3 &vglob) const;
index 098065e..b7b42b4 100644 (file)
@@ -61,22 +61,8 @@ class AliJetDummyShishKebabTrd1Module : public TNamed {
     } else      {xr = fOK1.Y(); zr = fOK1.X();
     }
   }
-  void GetCenterOfCellInLocalCoordinateofSM3X3(Int_t ieta, Double_t &xr, Double_t &zr) const
-  {
-    // Center of cell   
-    // 3X3 case - Nov 9,2006
-    ieta = ieta<0? 0 : ieta; // check index
-    ieta = ieta>2? 2 : ieta;
-    xr   = fOK3X3[2-ieta].Y();
-    zr   = fOK3X3[2-ieta].X();
-  }
-  void GetCenterOfCellInLocalCoordinateofSM1X1(Double_t &xr, Double_t &zr) const
-  {
-    // 1X1 case - Nov 27,2006
-    // Center of cell is center of module
-    xr   = fOK.Y() - fgr;
-    zr   = fOK.X();
-  }
+  void GetCenterOfCellInLocalCoordinateofSM3X3(Int_t ieta, Double_t &xr, Double_t &zr) const;
+  void GetCenterOfCellInLocalCoordinateofSM1X1(Double_t &xr, Double_t &zr) const;
   // 15-may-06
   TVector2& GetCenterOfModuleFace() {return fOB;}  
   TVector2& GetCenterOfModuleFace(Int_t ieta) {
@@ -122,5 +108,23 @@ class AliJetDummyShishKebabTrd1Module : public TNamed {
   Int_t    fDebug;    // debug flag
   ClassDef(AliJetDummyShishKebabTrd1Module,1) // TRD1 Shish-Kebab module 
 };
+inline
+  void AliJetDummyShishKebabTrd1Module::GetCenterOfCellInLocalCoordinateofSM3X3(Int_t ieta, Double_t &xr, Double_t &zr) const
+  {
+    // Center of cell   
+    // 3X3 case - Nov 9,2006
+    ieta = ieta<0? 0 : ieta; // check index
+    ieta = ieta>2? 2 : ieta;
+    xr   = fOK3X3[2-ieta].Y();
+    zr   = fOK3X3[2-ieta].X();
+  }
+inline
+  void AliJetDummyShishKebabTrd1Module::GetCenterOfCellInLocalCoordinateofSM1X1(Double_t &xr, Double_t &zr) const
+  {
+    // 1X1 case - Nov 27,2006
+    // Center of cell is center of module
+    xr   = fOK.Y() - fgr;
+    zr   = fOK.X();
+  }
 
 #endif
index bdad604..07322fb 100644 (file)
@@ -157,7 +157,7 @@ void AliJetFinder::AddJet(AliAODJet p)
   new ((*fAODjets)[fNAODjets++]) AliAODJet(p);
 }
 
-void AliJetFinder::ConnectAOD(AliAODEvent* aod)
+void AliJetFinder::ConnectAOD(const AliAODEvent* aod)
 {
 // Connect to the AOD
     fAODjets = aod->GetJets();
index b8d70d4..199c74c 100755 (executable)
@@ -55,7 +55,7 @@ class AliJetFinder : public TObject
   virtual Bool_t        ProcessEvent();
   virtual Bool_t        ProcessEvent2();
   virtual void          ConnectTree(TTree* tree, TObject* data);
-  virtual void          ConnectAOD(AliAODEvent* aod);
+  virtual void          ConnectAOD(const AliAODEvent* aod);
   virtual void          ConnectAODNonStd(AliAODEvent* aod,const char* bname);
   virtual void          WriteHeaders();
 
index e9bf549..7c25beb 100644 (file)
@@ -61,6 +61,7 @@ AliJetHadronCorrectionv0::Instance()
 Double_t 
 AliJetHadronCorrectionv0::GetEnergy(Double_t pmom, Double_t eta, Int_t /*gid*/)
 {
+// Get energy for momentum pmom and pseudorapidity eta
   Int_t iEta=0; // index 
   Double_t etaw = TMath::Abs(eta);
   if(etaw > etaGrid[maxVariant-1]) etaw = etaGrid[maxVariant-1];
index d8a3001..5de3834 100644 (file)
 
 class AliJetHadronCorrectionv0: public AliJetHadronCorrection {
 
-  private:
-  static AliJetHadronCorrectionv0* fHadrCorr;
+  public:
+  static  AliJetHadronCorrectionv0* Instance();
+  virtual Double_t GetEnergy(Double_t pmom, Double_t eta, Int_t gid); 
+  Double_t GetEnergy(Double_t pmom, Double_t eta){return GetEnergy(pmom,eta,7);}
+  virtual ~AliJetHadronCorrectionv0() {}
 
   protected:
   AliJetHadronCorrectionv0() {;}
   AliJetHadronCorrectionv0(const char *name, const char *title);
 
-  public:
-  static  AliJetHadronCorrectionv0* Instance();
-  virtual Double_t GetEnergy(Double_t pmom, Double_t eta, Int_t gid); 
-  Double_t GetEnergy(Double_t pmom, Double_t eta){return GetEnergy(pmom,eta,7);}
 
-  virtual ~AliJetHadronCorrectionv0() {}
+  private:
+    static AliJetHadronCorrectionv0* fHadrCorr; // pointer to instance
 
   ClassDef(AliJetHadronCorrectionv0,1) // Hadron correction for EMC (version for MDC)
 };
index d72dd72..c0f8874 100644 (file)
@@ -63,7 +63,7 @@ AliJetHadronCorrectionv1::Instance()
   return fgHadrCorr;
 }
 
-void AliJetHadronCorrectionv1::SetGeometry2(AliJetDummyGeo *geometry)
+void AliJetHadronCorrectionv1::SetGeometry2(const AliJetDummyGeo *geometry)
 {
   // Initialise EMCAL geometry
     if (!geometry)
@@ -168,8 +168,9 @@ Double_t AliJetHadronCorrectionv1::GetEnergy(Double_t pmom, Double_t eta, Int_t
    
 }
 
-void AliJetHadronCorrectionv1::TrackPositionEMCal(AliAODTrack* track,Double_t &eta, Double_t &phi)
+void AliJetHadronCorrectionv1::TrackPositionEMCal(const AliAODTrack* track,Double_t &eta, Double_t &phi)
 {
+// Return track position on EMCal
   AliAODPid*    pid   = (AliAODPid*) track->GetDetPid();
     
   if(pid) {
index 643425a..b210685 100644 (file)
@@ -23,8 +23,8 @@ class AliJetHadronCorrectionv1 : public AliJetHadronCorrection
   Double_t GetEnergy(Double_t pmom, Double_t eta){return GetEnergy(pmom,eta,7);}
   
   void SetGeometry(TString name, Double_t fs = 1.); 
-  void SetGeometry2(AliJetDummyGeo *geometry);
-  void TrackPositionEMCal(AliAODTrack* track,Double_t &eta, Double_t &phi);
+  void SetGeometry2(const AliJetDummyGeo *geometry);
+  void TrackPositionEMCal(const AliAODTrack* track,Double_t &eta, Double_t &phi);
 
  protected:
   AliJetHadronCorrectionv1():fSamplingFraction(0) {for (Int_t i = 0; i < 8; i++) fPar[i] = 0.;}
@@ -35,7 +35,7 @@ class AliJetHadronCorrectionv1 : public AliJetHadronCorrection
 
   static AliJetHadronCorrectionv1* fgHadrCorr;  // Pointer to global instance (singleton)
   static Double_t fgParLookup[HCPARAMETERS][HCPARAMETERSETS]; // Global array with parameters for hadronic response
-  Double_t fPar[8];
+  Double_t fPar[8];            // Parameters
   Float_t  fSamplingFraction;  // Sampling fraction
     
   ClassDef(AliJetHadronCorrectionv1,2) // Hadron correction for EMC (version for MDC)
index 4e794ea..448b22e 100644 (file)
@@ -54,6 +54,7 @@ AliJetHistos::AliJetHistos():
 
 AliJetHistos::~AliJetHistos()
 {
+// Destructor
   delete fNJetsH;
   delete fPtH;
   delete fEtaH;
@@ -74,6 +75,7 @@ void AliJetHistos::SetProperties(TH1* h,const char* x, const char* y) const
 
 void AliJetHistos::AddHistosToList(TList *list)
 {
+// Add histos to the list
   list->Add(fNJetsH);
   list->Add(fPtH);
   list->Add(fEtaH);
@@ -84,7 +86,7 @@ void AliJetHistos::AddHistosToList(TList *list)
 
 void AliJetHistos::FillHistos(TClonesArray *jets)
 {
-
+// Fill histograms
   if(!jets)return;
   Int_t nj = jets->GetEntries();
   fNJetsH->Fill(nj,1);
index ff4c13c..8f60288 100644 (file)
@@ -30,9 +30,9 @@ class AliJetKineReaderHeader : public AliJetReaderHeader
          
  protected:
   //parameters set by user
-  Bool_t fFastSimTPC;
-  Bool_t fFastSimEMCAL;
-  Bool_t fChargedOnly;
+  Bool_t fFastSimTPC;   // TPC fast simulation flag
+  Bool_t fFastSimEMCAL; // EMCAL fast simulation flag
+  Bool_t fChargedOnly;  // Charged particle only flag
 
   ClassDef(AliJetKineReaderHeader,2);
 };
index 9bc9bc7..5a1ceae 100644 (file)
@@ -23,7 +23,7 @@ class AliJetProductionData : public TObject
     void    GetPtHardLimits(Int_t bin, Float_t& ptmin, Float_t& ptmax);
     TString GetRunTitle(Int_t bin);
     Float_t GetWeight(Int_t bin);
-
+ protected:
     Int_t     fNbins;         // Number of pt_hard bins used in the production
     Float_t*  fPtHardLimits;  //[fNbins+1]
     Float_t*  fWeights;       //[fNbins]
index c6093df..0d4a460 100755 (executable)
@@ -72,7 +72,7 @@ class AliJetReaderHeader : public TNamed
   Int_t   fCluster;        // cluster type
   Int_t   fDebug;          // debug option
   Bool_t  fDZ;             // include dead zones or not 
-  Int_t   fSignalPerBg;
+  Int_t   fSignalPerBg;    // Signal per Background
   Float_t fFiducialEtaMin; // Fiducial minimum eta
   Float_t fFiducialEtaMax; // Fiducial maximum eta
   Float_t fFiducialPhiMin; // Fiducial minimum phi
index 1009bf7..fca4fba 100644 (file)
@@ -35,7 +35,7 @@ class AliJetSelector : public TSelector
   virtual void   Terminate();
   
  protected:
-  AliJetFinder* fJetFinder;
+  AliJetFinder* fJetFinder; // Pointer to jet finder
   
   ClassDef(AliJetSelector, 1)
 };
index 77d28f3..cdad53b 100644 (file)
@@ -254,7 +254,7 @@ Float_t  AliJetUnitArray::EtaToTheta(Float_t arg) const
 }
 
 //------------------------------------------------------------------------
-Bool_t AliJetUnitArray::operator>(AliJetUnitArray* unit) const
+Bool_t AliJetUnitArray::operator>(const AliJetUnitArray* unit) const
 {
   // Greater than operator used by sort
   if( fUnitEnergy > unit->GetUnitEnergy())
@@ -264,7 +264,7 @@ Bool_t AliJetUnitArray::operator>(AliJetUnitArray* unit) const
 }
 
 //------------------------------------------------------------------------
-Bool_t AliJetUnitArray::operator<( AliJetUnitArray* unit) const
+Bool_t AliJetUnitArray::operator<(const AliJetUnitArray* unit) const
 {
   // Less than operator used by sort
   if( fUnitEnergy < unit->GetUnitEnergy())
@@ -274,7 +274,7 @@ Bool_t AliJetUnitArray::operator<( AliJetUnitArray* unit) const
 }
 
 //------------------------------------------------------------------------
-Bool_t AliJetUnitArray::operator==( AliJetUnitArray* unit) const
+Bool_t AliJetUnitArray::operator==(const AliJetUnitArray* unit) const
 {
   // equality operator used by sort
   if( fUnitEnergy == unit->GetUnitEnergy())
index 804e377..98b6ae5 100644 (file)
@@ -131,9 +131,9 @@ class AliJetUnitArray : public TObject
  protected:
   AliJetUnitArray(const AliJetUnitArray& rUnit);
   AliJetUnitArray& operator = (const AliJetUnitArray& rhs);
-  Bool_t operator>  ( AliJetUnitArray* unit1) const;
-  Bool_t operator<  ( AliJetUnitArray* unit1) const;
-  Bool_t operator== ( AliJetUnitArray* unit1) const;
+  Bool_t operator>  (const AliJetUnitArray* unit1) const;
+  Bool_t operator<  (const AliJetUnitArray* unit1) const;
+  Bool_t operator== (const AliJetUnitArray* unit1) const;
 
   Float_t                            fUnitEnergy;       // Energy (Pt,et) of the unit 
   Float_t                            fUnitEta;          // Eta of the unit
index f23ec64..b326433 100755 (executable)
@@ -32,7 +32,7 @@ ClassImp(AliKMeansClustering)
 Double_t AliKMeansClustering::fBeta = 10.;
 
  
-Int_t AliKMeansClustering::SoftKMeans(Int_t k, Int_t n, Double_t* x, Double_t* y, Double_t* mx, Double_t* my , Double_t* rk )
+Int_t AliKMeansClustering::SoftKMeans(Int_t k, Int_t n, const Double_t* x, const Double_t* y, Double_t* mx, Double_t* my , Double_t* rk )
 {
     //
     // The soft K-means algorithm
@@ -397,7 +397,7 @@ Double_t AliKMeansClustering::d(Double_t mx, Double_t my, Double_t x, Double_t y
 
 
 
-void AliKMeansClustering::OptimalInit(Int_t k, Int_t n, Double_t* x, Double_t* y, Double_t* mx, Double_t* my)
+void AliKMeansClustering::OptimalInit(Int_t k, Int_t n, const Double_t* x, const Double_t* y, Double_t* mx, Double_t* my)
 {
   //  
   // Optimal initialisation using the k-means++ algorithm
@@ -520,7 +520,7 @@ void AliKMeansResult::Sort()
   TMath::Sort(fK, fTarget, fInd);
 }
 
-void AliKMeansResult::Sort(Int_t n, Double_t* x, Double_t* y)
+void AliKMeansResult::Sort(Int_t n, const Double_t* x, const Double_t* y)
 {
   // Build target array and sort
   for (Int_t i = 0; i < fK; i++)
@@ -541,7 +541,7 @@ void AliKMeansResult::Sort(Int_t n, Double_t* x, Double_t* y)
   TMath::Sort(fK, fTarget, fInd);
 }
 
-void AliKMeansResult::CopyResults(AliKMeansResult* res)
+void AliKMeansResult::CopyResults(const AliKMeansResult* res)
 {
   fK = res->GetK();
   for (Int_t i = 0; i <fK; i++) {
index 7340e21..52edef5 100755 (executable)
@@ -20,16 +20,16 @@ class AliKMeansClustering : public TObject
   AliKMeansClustering()          {}
   virtual ~AliKMeansClustering() {}
   
-  static Int_t SoftKMeans (Int_t k, Int_t n, Double_t* x, Double_t* y, Double_t* mx, Double_t* my , Double_t* rk );
+  static Int_t SoftKMeans (Int_t k, Int_t n, const Double_t* x, const Double_t* y, Double_t* mx, Double_t* my , Double_t* rk );
   static Int_t SoftKMeans2(Int_t k, Int_t n, Double_t* x, Double_t* y, Double_t* mx, Double_t* my , Double_t* sigma2, 
                          Double_t* rk );
   static Int_t SoftKMeans3(Int_t k, Int_t n, Double_t* x, Double_t* y, Double_t* mx, Double_t* my , 
                           Double_t* sigmax2, Double_t* sigmay2, Double_t* rk );
-  static void  OptimalInit(Int_t k, Int_t n, Double_t* x, Double_t* y, Double_t* mx, Double_t* my);
+  static void  OptimalInit(Int_t k, Int_t n, const Double_t* x, const Double_t* y, Double_t* mx, Double_t* my);
   static void  SetBeta(Double_t beta) {fBeta = beta;}
   static Double_t d(Double_t mx, Double_t my, Double_t x, Double_t y);
 protected:
-  static Double_t fBeta;
+  static Double_t fBeta; // beta parameter
   
   ClassDef(AliKMeansClustering, 1)
 };
@@ -49,9 +49,9 @@ class AliKMeansResult : public TObject
   Double_t*  GetRk()     const  {return fRk;}
   Int_t*     GetInd()    const  {return fInd;}
   Double_t*  GetTarget() const  {return fTarget;}
-  void       CopyResults(AliKMeansResult* res);
+  void       CopyResults(const AliKMeansResult* res);
   void       Sort();
-  void       Sort(Int_t n, Double_t* x, Double_t* y);  
+  void       Sort(Int_t n, const Double_t* x, const Double_t* y);  
 protected:
   Int_t        fK;        //! Number of clusters
   Double_t*    fMx;       //! Position x
index df60158..0e73b96 100644 (file)
@@ -1484,7 +1484,7 @@ void AliUA1JetFinderV2::SubtractBackgStat(const Int_t& nIn, const Int_t&nJ,Float
 
 ////////////////////////////////////////////////////////////////////////
 void AliUA1JetFinderV2::SubtractBackgCone(const Int_t& nIn, const Int_t&nJ,Float_t& etbgTotalN,
-                      Float_t* ptT, Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
+                      const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
                       Float_t* const etJet, const Float_t* etaJet, const Float_t* phiJet,
                       Float_t* const etsigJet, Int_t* const multJet, Int_t* const injet)
 {
@@ -1598,7 +1598,7 @@ void AliUA1JetFinderV2::SubtractBackgCone(const Int_t& nIn, const Int_t&nJ,Float
 
 ////////////////////////////////////////////////////////////////////////
 void AliUA1JetFinderV2::SubtractBackgRatio(const Int_t& nIn, const Int_t&nJ,Float_t& etbgTotalN,
-                      Float_t* ptT, Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
+                      const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
                       Float_t* const etJet, const Float_t* etaJet, const Float_t* phiJet,
                       Float_t* const etsigJet, Int_t* const multJet, Int_t* const injet)
 {
index cc9d85f..df055fc 100644 (file)
@@ -49,12 +49,12 @@ class AliUA1JetFinderV2 : public AliJetFinder
                     Float_t* const etsigJet, Int_t* const multJet, Int_t* const injet);
 
   void SubtractBackgCone(const Int_t& nIn, const Int_t&nJ,Float_t& EtbgTotalN,
-                      Float_t* ptT, Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
+                      const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
                       Float_t* const etJet, const Float_t* etaJet, const Float_t* phiJet,
                       Float_t* const etsigJet, Int_t* const multJet, Int_t* const injet);
 
   void SubtractBackgRatio(const Int_t& nIn, const Int_t&nJ,Float_t& EtbgTotalN,
-                      Float_t* ptT, Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
+                      const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
                       Float_t* const etJet, const Float_t* etaJet, const Float_t* phiJet,
                       Float_t* const etsigJet, Int_t* const multJet, Int_t* const injet);