Coding Conventions
authorfca <fca@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 12 Apr 2010 14:24:11 +0000 (14:24 +0000)
committerfca <fca@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 12 Apr 2010 14:24:11 +0000 (14:24 +0000)
TPC/AliExternalComparison.cxx
TPC/AliExternalComparison.h
TPC/AliTPCCalROC.cxx
TPC/AliTPCCalROC.h

index dc90c5c..731ebe4 100644 (file)
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
-/*
-
-
-
-
-*/
-
-
+/* $Id:$ */
 
 #include "AliExternalComparison.h"
 #include "AliExternalTrackParam.h"
@@ -62,7 +55,7 @@ AliExternalComparison::AliExternalComparison(const Text_t *name, const Text_t *t
 }
 
 AliExternalComparison::AliExternalComparison(const AliExternalComparison& comp)
-  :TNamed(comp.fName,comp.fTitle),
+  :TNamed(comp),
    fResolHistos(new TObjArray(*(comp.fResolHistos))),
    fPullHistos(new TObjArray(*(comp.fPullHistos))),
    fRangeMatrix(new TMatrixD(*(comp.fRangeMatrix))),
@@ -73,17 +66,18 @@ AliExternalComparison::AliExternalComparison(const AliExternalComparison& comp)
   //
 }
 
-AliExternalComparison& AliExternalComparison::operator=(const AliExternalComparison&comp)
+AliExternalComparison& AliExternalComparison::operator=(const AliExternalComparison& comp)
 {
   //
   //
   //
-  SetName(comp.GetName());
-  SetTitle(comp.GetTitle());
-  fResolHistos=new TObjArray(*(comp.fResolHistos));
-  fPullHistos=new TObjArray(*(comp.fPullHistos));
-  fRangeMatrix=new TMatrixD(*(comp.fRangeMatrix));
-  fCutMatrix=new TMatrixD(*(comp.fCutMatrix));
+  if(this != &comp) {
+    TNamed::operator=(comp);
+    fResolHistos=new TObjArray(*(comp.fResolHistos));
+    fPullHistos=new TObjArray(*(comp.fPullHistos));
+    fRangeMatrix=new TMatrixD(*(comp.fRangeMatrix));
+    fCutMatrix=new TMatrixD(*(comp.fCutMatrix));
+  }
   return *this;
 }
 
@@ -108,7 +102,7 @@ void AliExternalComparison::Analyze() {
 }
 
 
-Long64_t AliExternalComparison::Merge(TCollection *li) {
+Long64_t AliExternalComparison::Merge(TCollection *const li) {
   //
   //
   //   
@@ -126,7 +120,7 @@ Long64_t AliExternalComparison::Merge(TCollection *li) {
   return 0;  
 }
 
-void  AliExternalComparison:: Add(AliExternalComparison*comp){
+void  AliExternalComparison::Add(AliExternalComparison *const comp){
   //
   // Add comparison object
   //
index bf232ca..4041df8 100644 (file)
@@ -6,7 +6,7 @@
 
 
 #include "TNamed.h"
-#include "TMatrixD.h"
+#include "TMatrixDfwd.h"
 class THnSparse;
 class TCollection;
 class AliExternalTrackParam;
@@ -18,12 +18,12 @@ class AliExternalComparison:public TNamed {
 public:
   AliExternalComparison(); 
   AliExternalComparison(const Text_t *name, const Text_t *title);
-  AliExternalComparison(const AliExternalComparison&);
-  AliExternalComparison& operator=(const AliExternalComparison&);
+  AliExternalComparison(const AliExternalComparison& comp);
+  AliExternalComparison& operator=(const AliExternalComparison& comp);
   //
   virtual ~AliExternalComparison();
-  virtual Long64_t       Merge(TCollection *li);
-  virtual void           Add(AliExternalComparison*comp);
+  virtual Long64_t       Merge(TCollection *const li);
+  virtual void           Add(AliExternalComparison *const comp);
   virtual void           Analyze();
   virtual void           Process(const AliExternalTrackParam *param0, const AliExternalTrackParam *param1);
   virtual void           Process(const AliExternalTrackParam *param0, TParticle *part);
@@ -45,7 +45,6 @@ public:
 
 protected:
   void    MakeHistos();
-public:
   TObjArray * fResolHistos;             // resolution histogram
   TObjArray * fPullHistos;              // pull       histogram
   TMatrixD  * fRangeMatrix;             // range matrix
index eb76054..5928d9d 100644 (file)
@@ -48,7 +48,7 @@ AliTPCCalROC::AliTPCCalROC()
              fSector(0),
              fNChannels(0),
              fNRows(0),
-             fIndexes(0),
+             fkIndexes(0),
              fData(0)
 {
   //
@@ -63,7 +63,7 @@ AliTPCCalROC::AliTPCCalROC(UInt_t  sector)
              fSector(0),
              fNChannels(0),
              fNRows(0),
-             fIndexes(0),
+             fkIndexes(0),
              fData(0)
 {
   //
@@ -72,7 +72,7 @@ AliTPCCalROC::AliTPCCalROC(UInt_t  sector)
   fSector = sector;
   fNChannels    =  AliTPCROC::Instance()->GetNChannels(fSector);
   fNRows        =  AliTPCROC::Instance()->GetNRows(fSector);
-  fIndexes      =  AliTPCROC::Instance()->GetRowIndexes(fSector);
+  fkIndexes      =  AliTPCROC::Instance()->GetRowIndexes(fSector);
   fData = new Float_t[fNChannels];
   for (UInt_t  idata = 0; idata< fNChannels; idata++) fData[idata] = 0.;
 }
@@ -83,7 +83,7 @@ AliTPCCalROC::AliTPCCalROC(const AliTPCCalROC &c)
              fSector(0),
              fNChannels(0),
              fNRows(0),
-             fIndexes(0),
+             fkIndexes(0),
              fData(0)
 {
   //
@@ -92,7 +92,7 @@ AliTPCCalROC::AliTPCCalROC(const AliTPCCalROC &c)
   fSector = c.fSector;
   fNChannels    =  AliTPCROC::Instance()->GetNChannels(fSector);
   fNRows        =  AliTPCROC::Instance()->GetNRows(fSector);
-  fIndexes      =  AliTPCROC::Instance()->GetRowIndexes(fSector);
+  fkIndexes      =  AliTPCROC::Instance()->GetRowIndexes(fSector);
   //
   fData   = new Float_t[fNChannels];
   for (UInt_t  idata = 0; idata< fNChannels; idata++) fData[idata] = c.fData[idata];
@@ -129,7 +129,7 @@ void AliTPCCalROC::Streamer(TBuffer &R__b)
    //
    if (R__b.IsReading()) {
       AliTPCCalROC::Class()->ReadBuffer(R__b, this);
-      fIndexes =  AliTPCROC::Instance()->GetRowIndexes(fSector);
+      fkIndexes =  AliTPCROC::Instance()->GetRowIndexes(fSector);
    } else {
       AliTPCCalROC::Class()->WriteBuffer(R__b,this);
    }
@@ -188,7 +188,7 @@ void AliTPCCalROC::Divide(const AliTPCCalROC*  roc) {
   }
 }
 
-Double_t AliTPCCalROC::GetMean(AliTPCCalROC* outlierROC) {
+Double_t AliTPCCalROC::GetMean(AliTPCCalROC *const outlierROC) const {
    //
    //  returns the mean value of the ROC
    //  pads with value != 0 in outlierROC are not used for the calculation
@@ -196,21 +196,21 @@ Double_t AliTPCCalROC::GetMean(AliTPCCalROC* outlierROC) {
    //
    if (!outlierROC) return TMath::Mean(fNChannels, fData);
    Double_t *ddata = new Double_t[fNChannels];
-   Int_t NPoints = 0;
+   Int_t nPoints = 0;
    for (UInt_t i=0;i<fNChannels;i++) {
       if (!(outlierROC->GetValue(i))) {
-         ddata[NPoints]= fData[i];
-         NPoints++;
+         ddata[nPoints]= fData[i];
+         nPoints++;
       }
    }
    Double_t mean = 0;
-   if(NPoints>0)
-     mean = TMath::Mean(NPoints, ddata);
+   if(nPoints>0)
+     mean = TMath::Mean(nPoints, ddata);
    delete [] ddata;
    return mean;
 }
 
-Double_t AliTPCCalROC::GetMedian(AliTPCCalROC* outlierROC) {
+Double_t AliTPCCalROC::GetMedian(AliTPCCalROC *const outlierROC) const {
    //
    //  returns the median value of the ROC
    //  pads with value != 0 in outlierROC are not used for the calculation
@@ -218,21 +218,21 @@ Double_t AliTPCCalROC::GetMedian(AliTPCCalROC* outlierROC) {
    //
    if (!outlierROC) return TMath::Median(fNChannels, fData);
    Double_t *ddata = new Double_t[fNChannels];
-   Int_t NPoints = 0;
+   Int_t nPoints = 0;
    for (UInt_t i=0;i<fNChannels;i++) {
       if (!(outlierROC->GetValue(i))) {
-         ddata[NPoints]= fData[i];
-         NPoints++;
+         ddata[nPoints]= fData[i];
+         nPoints++;
       }
    }
    Double_t median = 0;
-   if(NPoints>0)
-     median = TMath::Median(NPoints, ddata);
+   if(nPoints>0)
+     median = TMath::Median(nPoints, ddata);
    delete [] ddata;
    return median;
 }
 
-Double_t AliTPCCalROC::GetRMS(AliTPCCalROC* outlierROC) {
+Double_t AliTPCCalROC::GetRMS(AliTPCCalROC *const outlierROC) const {
    //
    //  returns the RMS value of the ROC
    //  pads with value != 0 in outlierROC are not used for the calculation
@@ -240,39 +240,39 @@ Double_t AliTPCCalROC::GetRMS(AliTPCCalROC* outlierROC) {
    //
    if (!outlierROC) return TMath::RMS(fNChannels, fData);
    Double_t *ddata = new Double_t[fNChannels];
-   Int_t NPoints = 0;
+   Int_t nPoints = 0;
    for (UInt_t i=0;i<fNChannels;i++) {
       if (!(outlierROC->GetValue(i))) {
-         ddata[NPoints]= fData[i];
-         NPoints++;
+         ddata[nPoints]= fData[i];
+         nPoints++;
       }
    }
    Double_t rms = 0;
-   if(NPoints>0)
-     rms = TMath::RMS(NPoints, ddata);
+   if(nPoints>0)
+     rms = TMath::RMS(nPoints, ddata);
    delete [] ddata;
    return rms;
 }
 
-Double_t AliTPCCalROC::GetLTM(Double_t *sigma, Double_t fraction, AliTPCCalROC* outlierROC){
+Double_t AliTPCCalROC::GetLTM(Double_t *const sigma, Double_t fraction, AliTPCCalROC *const outlierROC){
   //
   //  returns the LTM and sigma
   //  pads with value != 0 in outlierROC are not used for the calculation
    //  return 0 if no data is accepted by the outlier cuts 
   //
   Double_t *ddata = new Double_t[fNChannels];
-  UInt_t NPoints = 0;
+  UInt_t nPoints = 0;
   for (UInt_t i=0;i<fNChannels;i++) {
      if (!outlierROC || !(outlierROC->GetValue(i))) {
-        ddata[NPoints]= fData[i];
-        NPoints++;
+        ddata[nPoints]= fData[i];
+        nPoints++;
      }
   }
 
   Double_t ltm =0, lsigma=0;
-  if(NPoints>0) {
-    Int_t hh = TMath::Min(TMath::Nint(fraction *NPoints), Int_t(NPoints));
-    AliMathBase::EvaluateUni(NPoints,ddata, ltm, lsigma, hh);
+  if(nPoints>0) {
+    Int_t hh = TMath::Min(TMath::Nint(fraction *nPoints), Int_t(nPoints));
+    AliMathBase::EvaluateUni(nPoints,ddata, ltm, lsigma, hh);
     if (sigma) *sigma=lsigma;
   }
   
@@ -581,20 +581,20 @@ AliTPCCalROC * AliTPCCalROC::LocalFit(Int_t rowRadius, Int_t padRadius, AliTPCCa
   // chi2Threshold: Threshold for chi2 when EvalRobust is called
   // robustFraction: Fraction of data that will be used in EvalRobust
   //
-  AliTPCCalROC * ROCfitted = new AliTPCCalROC(fSector);
+  AliTPCCalROC * xROCfitted = new AliTPCCalROC(fSector);
   TLinearFitter fitterQ(6,"hyp5");
    // TLinearFitter fitterQ(6,"x0++x1++x2++x3++x4++x5");
   fitterQ.StoreData(kTRUE);
   for (UInt_t row=0; row < GetNrows(); row++) {
     //std::cout << "Entering row " << row << " of " << GetNrows() << " @ sector "<< fSector << " for local fitting... "<< std::endl;
     for (UInt_t pad=0; pad < GetNPads(row); pad++)
-      ROCfitted->SetValue(row, pad, GetNeighbourhoodValue(&fitterQ, row, pad, rowRadius, padRadius, ROCoutliers, robust, chi2Threshold, robustFraction));
+      xROCfitted->SetValue(row, pad, GetNeighbourhoodValue(&fitterQ, row, pad, rowRadius, padRadius, ROCoutliers, robust, chi2Threshold, robustFraction));
   }
-  return ROCfitted;
+  return xROCfitted;
 }
 
 
-Double_t AliTPCCalROC::GetNeighbourhoodValue(TLinearFitter* fitterQ, Int_t row, Int_t pad, Int_t rRadius, Int_t pRadius, AliTPCCalROC* ROCoutliers, Bool_t robust, Double_t chi2Threshold, Double_t robustFraction) {
+Double_t AliTPCCalROC::GetNeighbourhoodValue(TLinearFitter* fitterQ, Int_t row, Int_t pad, Int_t rRadius, Int_t pRadius, AliTPCCalROC *const ROCoutliers, Bool_t robust, Double_t chi2Threshold, Double_t robustFraction) {
   //
   // AliTPCCalROC::GetNeighbourhoodValue - smoothing - PRIVATE
   // in this function the fit for LocalFit is done
@@ -804,35 +804,35 @@ AliTPCCalROC* AliTPCCalROC::CreateGlobalFitCalROC(TVectorD &fitParam, Int_t sect
   Float_t dlx, dly;
   Float_t centerPad[3] = {0};
   Float_t localXY[3] = {0};
-  AliTPCCalROC * ROCfitted = new AliTPCCalROC(sector);
+  AliTPCCalROC * xROCfitted = new AliTPCCalROC(sector);
   AliTPCROC* tpcROCinstance = AliTPCROC::Instance();
-  tpcROCinstance->GetPositionLocal(sector, ROCfitted->GetNrows()/2, ROCfitted->GetNPads(ROCfitted->GetNrows()/2)/2, centerPad);  // calculate center of ROC 
+  tpcROCinstance->GetPositionLocal(sector, xROCfitted->GetNrows()/2, xROCfitted->GetNPads(xROCfitted->GetNrows()/2)/2, centerPad);  // calculate center of ROC 
   Int_t fitType = 1;
   if (fitParam.GetNoElements() == 6) fitType = 1;
   else fitType = 0;
   Double_t value = 0;
   if (fitType == 1) { // parabolic fit
-    for (UInt_t irow = 0; irow < ROCfitted->GetNrows(); irow++) {
-      for (UInt_t ipad = 0; ipad < ROCfitted->GetNPads(irow); ipad++) {
+    for (UInt_t irow = 0; irow < xROCfitted->GetNrows(); irow++) {
+      for (UInt_t ipad = 0; ipad < xROCfitted->GetNPads(irow); ipad++) {
        tpcROCinstance->GetPositionLocal(sector, irow, ipad, localXY);   // calculate position localXY by pad and row number
        dlx = localXY[0] - centerPad[0];
        dly = localXY[1] - centerPad[1];
        value = fitParam[0] + fitParam[1]*dlx + fitParam[2]*dly + fitParam[3]*dlx*dlx + fitParam[4]*dly*dly + fitParam[5]*dlx*dly;
-       ROCfitted->SetValue(irow, ipad, value);
+       xROCfitted->SetValue(irow, ipad, value);
       }
     }   
   }
   else {  // linear fit
-    for (UInt_t irow = 0; irow < ROCfitted->GetNrows(); irow++) {
-      for (UInt_t ipad = 0; ipad < ROCfitted->GetNPads(irow); ipad++) {
+    for (UInt_t irow = 0; irow < xROCfitted->GetNrows(); irow++) {
+      for (UInt_t ipad = 0; ipad < xROCfitted->GetNPads(irow); ipad++) {
        tpcROCinstance->GetPositionLocal(sector, irow, ipad, localXY);   // calculate position localXY by pad and row number
        dlx = localXY[0] - centerPad[0];
        dly = localXY[1] - centerPad[1];
        value = fitParam[0] + fitParam[1]*dlx + fitParam[2]*dly;
-       ROCfitted->SetValue(irow, ipad, value);
+       xROCfitted->SetValue(irow, ipad, value);
       }
     }   
   }
-  return ROCfitted;
+  return xROCfitted;
 }
 
index 111a2f5..7752f6d 100644 (file)
@@ -32,9 +32,9 @@ class AliTPCCalROC : public TNamed {
   UInt_t        GetNrows() const               { return fNRows;};
   UInt_t        GetNchannels()       const     { return fNChannels;};
   UInt_t        GetNPads(UInt_t row)  const     { return (row<fNRows)? AliTPCROC::Instance()->GetNPads(fSector,row):0;};
-  Float_t      GetValue(UInt_t row, UInt_t pad) const { return ( (row<fNRows) && (fIndexes[row]+pad)<fNChannels)? fData[fIndexes[row]+pad]: 0; };
+  Float_t      GetValue(UInt_t row, UInt_t pad) const { return ( (row<fNRows) && (fkIndexes[row]+pad)<fNChannels)? fData[fkIndexes[row]+pad]: 0; };
   Float_t      GetValue(UInt_t channel) const { return  fData[channel]; };
-  void         SetValue(UInt_t row, UInt_t pad, Float_t vd) { if ( row<fNRows && (fIndexes[row]+pad)<fNChannels)fData[fIndexes[row]+pad]= vd; };
+  void         SetValue(UInt_t row, UInt_t pad, Float_t vd) { if ( row<fNRows && (fkIndexes[row]+pad)<fNChannels)fData[fkIndexes[row]+pad]= vd; };
   void         SetValue(UInt_t channel, Float_t vd) {fData[channel]= vd; };
   virtual void Draw(Option_t* option = "");
   //
@@ -46,10 +46,10 @@ class AliTPCCalROC : public TNamed {
   void Divide(const AliTPCCalROC * roc);   // divide each channel of the ROC by the coresponding value of 'roc'
   // statistic
   //
-  Double_t GetMean(AliTPCCalROC* outlierROC = 0);
-  Double_t GetRMS(AliTPCCalROC* outlierROC = 0);
-  Double_t GetMedian(AliTPCCalROC* outlierROC = 0) ;
-  Double_t GetLTM(Double_t *sigma=0, Double_t fraction=0.9, AliTPCCalROC* outlierROC = 0);
+  Double_t GetMean(AliTPCCalROC *const outlierROC = 0) const;
+  Double_t GetRMS(AliTPCCalROC *const outlierROC = 0) const;
+  Double_t GetMedian(AliTPCCalROC *const outlierROC = 0) const;
+  Double_t GetLTM(Double_t *const sigma=0, Double_t fraction=0.9, AliTPCCalROC *const outlierROC = 0);
   TH1F * MakeHisto1D(Float_t min=4, Float_t max=-4, Int_t type=0);     
   TH2F * MakeHisto2D(Float_t min=4, Float_t max=-4, Int_t type=0);   
   TH2F * MakeHistoOutliers(Float_t delta=4, Float_t fraction=0.7, Int_t mode=0);
@@ -62,13 +62,13 @@ class AliTPCCalROC : public TNamed {
   static void Test();
  protected:
   
-  Double_t GetNeighbourhoodValue(TLinearFitter* fitterQ, Int_t row, Int_t pad, Int_t rRadius, Int_t pRadius, AliTPCCalROC* ROCoutliers, Bool_t robust, Double_t chi2Threshold, Double_t robustFraction);
+  Double_t GetNeighbourhoodValue(TLinearFitter* fitterQ, Int_t row, Int_t pad, Int_t rRadius, Int_t pRadius, AliTPCCalROC *const ROCoutliers, Bool_t robust, Double_t chi2Threshold, Double_t robustFraction);
   void GetNeighbourhood(TArrayI* &rowArray, TArrayI* &padArray, Int_t row, Int_t pad, Int_t rRadius, Int_t pRadius);
   
   UInt_t     fSector;          // sector number
   UInt_t     fNChannels;       // number of channels
   UInt_t     fNRows;           // number of rows
-  const UInt_t* fIndexes;      //!indexes
+  const UInt_t* fkIndexes;      //!indexes
   Float_t  *fData;            //[fNChannels] Data
   ClassDef(AliTPCCalROC,2)    //  TPC ROC calibration class