Added debug flags
authorcholm <cholm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 29 Nov 2010 21:56:18 +0000 (21:56 +0000)
committercholm <cholm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 29 Nov 2010 21:56:18 +0000 (21:56 +0000)
PWG2/FORWARD/analysis2/AliFMDCorrections.cxx
PWG2/FORWARD/analysis2/AliFMDDensityCalculator.cxx
PWG2/FORWARD/analysis2/AliFMDDensityCalculator.h
PWG2/FORWARD/analysis2/AliFMDHistCollector.cxx
PWG2/FORWARD/analysis2/AliFMDHistCollector.h
PWG2/FORWARD/analysis2/AliFMDSharingFilter.cxx
PWG2/FORWARD/analysis2/AliFMDSharingFilter.h

index 3b3afb265067644d921f851e0c1eaab66edb6081..06dc4d1ce2ce39271227758649d07ff5c5e0c01a 100644 (file)
@@ -16,14 +16,16 @@ ClassImp(AliFMDCorrections)
 AliFMDCorrections::AliFMDCorrections()
   : TNamed(), 
     fRingHistos(),
-    fMultCut(0.3)
+    fMultCut(0.3),
+    fDebug(0)
 {}
 
 //____________________________________________________________________
 AliFMDCorrections::AliFMDCorrections(const char* title)
   : TNamed("fmdCorrections", title), 
     fRingHistos(), 
-    fMultCut(0.3)
+    fMultCut(0.3),
+    fDebug(0)
 {
   fRingHistos.SetName(GetName());
   fRingHistos.Add(new RingHistos(1, 'I'));
@@ -37,7 +39,8 @@ AliFMDCorrections::AliFMDCorrections(const char* title)
 AliFMDCorrections::AliFMDCorrections(const AliFMDCorrections& o)
   : TNamed(o), 
     fRingHistos(), 
-    fMultCut(o.fMultCut)
+    fMultCut(o.fMultCut),
+    fDebug(o.fDebug)
 {
   TIter    next(&o.fRingHistos);
   TObject* obj = 0;
@@ -54,11 +57,10 @@ AliFMDCorrections::~AliFMDCorrections()
 AliFMDCorrections&
 AliFMDCorrections::operator=(const AliFMDCorrections& o)
 {
-  SetName(o.GetName());
-  SetTitle(o.GetTitle());
+  TNamed::operator=(o);
 
   fMultCut = o.fMultCut;
-
+  fDebug   = o.fDebug;
   fRingHistos.Delete();
   TIter    next(&o.fRingHistos);
   TObject* obj = 0;
index 5da0a2f443e71ee89d9a96f55097ede69e6d3b86..a40651d853d23ce24ade5007abaf6e7d2f6502a4 100644 (file)
@@ -19,7 +19,8 @@ AliFMDDensityCalculator::AliFMDDensityCalculator()
     fMultCut(0.3),
     fSumOfWeights(0),
     fWeightedSum(0),
-    fCorrections(0)
+    fCorrections(0),
+    fDebug(0)
 {}
 
 //____________________________________________________________________
@@ -29,7 +30,8 @@ AliFMDDensityCalculator::AliFMDDensityCalculator(const char* title)
     fMultCut(0.3),
     fSumOfWeights(0),
     fWeightedSum(0),
-    fCorrections(0)
+    fCorrections(0),
+    fDebug(0)
 {
   fRingHistos.SetName(GetName());
   fRingHistos.Add(new RingHistos(1, 'I'));
@@ -54,7 +56,8 @@ AliFMDDensityCalculator::AliFMDDensityCalculator(const
     fMultCut(o.fMultCut),
     fSumOfWeights(o.fSumOfWeights),
     fWeightedSum(o.fWeightedSum),
-    fCorrections(o.fCorrections)
+    fCorrections(o.fCorrections),
+    fDebug(o.fDebug)
 {
   TIter    next(&o.fRingHistos);
   TObject* obj = 0;
@@ -71,10 +74,10 @@ AliFMDDensityCalculator::~AliFMDDensityCalculator()
 AliFMDDensityCalculator&
 AliFMDDensityCalculator::operator=(const AliFMDDensityCalculator& o)
 {
-  SetName(o.GetName());
-  SetTitle(o.GetTitle());
+  TNamed::operator=(o);
 
   fMultCut = o.fMultCut;
+  fDebug   = o.fDebug;
 
   fRingHistos.Delete();
   TIter    next(&o.fRingHistos);
index e2ec51265982b20763c8f67a9811774e00beb55e..4b761df484a6c14ef16be22262012e570467e3dc 100644 (file)
@@ -81,6 +81,12 @@ public:
    * @param dir List to write in
    */  
   void DefineOutput(TList* dir);
+  /** 
+   * Set the debug level.  The higher the value the more output 
+   * 
+   * @param dbg Debug level 
+   */
+  void SetDebug(Int_t dbg=1) { fDebug = dbg; }
 protected:
   /** 
    * Get the number of particles corresponding to the signal mult
@@ -184,11 +190,12 @@ protected:
    */
   RingHistos* GetRingHistos(UShort_t d, Char_t r) const;
 
-  TList    fRingHistos;    // List of histogram containers
-  Double_t fMultCut;       // Low cut on scaled energy loss
-  TH1D*    fSumOfWeights;
-  TH1D*    fWeightedSum;
-  TH1D*    fCorrections;
+  TList    fRingHistos;    //  List of histogram containers
+  Double_t fMultCut;       //  Low cut on scaled energy loss
+  TH1D*    fSumOfWeights;  //! Histogram
+  TH1D*    fWeightedSum;   //! Histogram
+  TH1D*    fCorrections;   //! Histogram
+  Int_t    fDebug;         //  Debug level 
 
   ClassDef(AliFMDDensityCalculator,1); // Calculate Nch density 
 };
index df71a5404977b4328f94c2a492bbcebeedba429a..13840f5b599828c831ea324996eb685501384d14 100644 (file)
@@ -22,14 +22,14 @@ ClassImp(AliFMDHistCollector)
 AliFMDHistCollector&
 AliFMDHistCollector::operator=(const AliFMDHistCollector& o)
 {
-  SetName(o.GetName());
-  SetTitle(o.GetTitle());
+  TNamed::operator=(o);
 
   fNCutBins       = o.fNCutBins;
   fCorrectionCut  = o.fCorrectionCut;
   fFirstBins      = o.fFirstBins;
   fLastBins       = o.fLastBins;
-  
+  fDebug          = o.fDebug;
+
   return *this;
 }
 
index 89441936321b21656572a5c710b6c3a00a1f4124..ad2c3b6b1131bcd6220ccf97297279aec0a19003 100644 (file)
@@ -30,7 +30,7 @@ public:
    * Constructor 
    */
   AliFMDHistCollector() 
-    : fNCutBins(0), fCorrectionCut(0), fFirstBins(), fLastBins()
+    : fNCutBins(0), fCorrectionCut(0), fFirstBins(), fLastBins(), fDebug(0)
   {}
   /** 
    * Constructor 
@@ -39,7 +39,8 @@ public:
    */
   AliFMDHistCollector(const char* title)
     : TNamed("fmdHistCollector", title), 
-      fNCutBins(1), fCorrectionCut(0.5), fFirstBins(1), fLastBins(1)
+      fNCutBins(1), fCorrectionCut(0.5), 
+      fFirstBins(1), fLastBins(1), fDebug(0)
   {}
   /** 
    * Copy constructor 
@@ -49,7 +50,7 @@ public:
   AliFMDHistCollector(const AliFMDHistCollector& o)
     : TNamed(o), 
       fNCutBins(o.fNCutBins), fCorrectionCut(o.fCorrectionCut),
-      fFirstBins(o.fFirstBins), fLastBins(o.fLastBins) 
+      fFirstBins(o.fFirstBins), fLastBins(o.fLastBins), fDebug(o.fDebug) 
   {}
 
   /** 
@@ -94,6 +95,12 @@ public:
    * @param cut Cut-off 
    */
   void SetCorrectionCut(Float_t cut=0.5) { fCorrectionCut = cut; }
+  /** 
+   * Set the debug level.  The higher the value the more output 
+   * 
+   * @param dbg Debug level 
+   */
+  void SetDebug(Int_t dbg=1) { fDebug = dbg; }
 protected:
   /** 
    * Get the first and last eta bin to use for a given ring and vertex 
@@ -223,6 +230,7 @@ protected:
   Float_t     fCorrectionCut;   // Cut-off on secondary corrections 
   TArrayI     fFirstBins;       // Array of first eta bins 
   TArrayI     fLastBins;        // Array of last eta bins 
+  Int_t       fDebug;           // Debug level 
 
   ClassDef(AliFMDHistCollector,1); // Calculate Nch density 
 };
index dcb0863be78fe2839fbea389a3d853dba6b56d9d..a3f07c40bc30b60d097974602a819d7f4077b424 100644 (file)
@@ -21,8 +21,8 @@ AliFMDSharingFilter::AliFMDSharingFilter()
   : TNamed(), 
     fRingHistos(),
     fLowCut(0.3),
-    fCorrectAngles(kFALSE),
-    fEtaCorr(0)
+    fCorrectAngles(kFALSE), 
+    fDebug(0)
 {}
 
 //____________________________________________________________________
@@ -30,18 +30,14 @@ AliFMDSharingFilter::AliFMDSharingFilter(const char* title)
   : TNamed("fmdSharingFilter", title), 
     fRingHistos(), 
     fLowCut(0.3),
-    fCorrectAngles(kFALSE),
-    fEtaCorr(0)
+    fCorrectAngles(kFALSE), 
+    fDebug(0)
 {
   fRingHistos.Add(new RingHistos(1, 'I'));
   fRingHistos.Add(new RingHistos(2, 'I'));
   fRingHistos.Add(new RingHistos(2, 'O'));
   fRingHistos.Add(new RingHistos(3, 'I'));
   fRingHistos.Add(new RingHistos(3, 'O'));
-  fEtaCorr = new TH2F("etaCorr", "Correction of #eta", 
-                     200, -4, 6,  200, -4, 6);
-  fEtaCorr->SetXTitle("#eta from ESD");
-  fEtaCorr->SetYTitle("#eta from AnaParameters");
 }
 
 //____________________________________________________________________
@@ -49,8 +45,8 @@ AliFMDSharingFilter::AliFMDSharingFilter(const AliFMDSharingFilter& o)
   : TNamed(o), 
     fRingHistos(), 
     fLowCut(o.fLowCut),
-    fCorrectAngles(o.fCorrectAngles),
-    fEtaCorr(o.fEtaCorr)
+    fCorrectAngles(o.fCorrectAngles), 
+    fDebug(o.fDebug)
 {
   TIter    next(&o.fRingHistos);
   TObject* obj = 0;
@@ -67,11 +63,11 @@ AliFMDSharingFilter::~AliFMDSharingFilter()
 AliFMDSharingFilter&
 AliFMDSharingFilter::operator=(const AliFMDSharingFilter& o)
 {
-  SetName(o.GetName());
-  SetTitle(o.GetTitle());
+  TNamed::operator=(o);
 
   fLowCut        = o.fLowCut;
   fCorrectAngles = o.fCorrectAngles;
+  fDebug         = o.fDebug;
 
   fRingHistos.Delete();
   TIter    next(&o.fRingHistos);
@@ -100,8 +96,7 @@ AliFMDSharingFilter::GetRingHistos(UShort_t d, Char_t r) const
 Bool_t
 AliFMDSharingFilter::Filter(const AliESDFMD& input, 
                            Bool_t           lowFlux,
-                           AliESDFMD&       output,
-                           Double_t         vz)
+                           AliESDFMD&       output)
 {
   output.Clear();
   TIter    next(&fRingHistos);
@@ -109,7 +104,7 @@ AliFMDSharingFilter::Filter(const AliESDFMD& input,
   while ((o = static_cast<RingHistos*>(next())))
     o->Clear();
 
-  AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();
+  // AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();
 
   for(UShort_t d = 1; d <= 3; d++) {
     Int_t nRings = (d == 1 ? 1 : 2);
@@ -136,10 +131,7 @@ AliFMDSharingFilter::Filter(const AliESDFMD& input,
          if (mult == AliESDFMD::kInvalidMult || mult == 0) continue;
 
          // Get the pseudo-rapidity 
-         Double_t eta1 = input.Eta(d,r,s,t);
-         Double_t eta2 = pars->GetEtaFromStrip(d,r,s,t,vz);
-         fEtaCorr->Fill(eta1, eta2);
-         Double_t eta = eta2;
+         Double_t eta = input.Eta(d,r,s,t);
          
          // Fill the diagnostics histogram 
          histos->fBefore->Fill(mult);
@@ -343,7 +335,6 @@ AliFMDSharingFilter::DefineOutput(TList* dir)
   TList* d = new TList;
   d->SetName(GetName());
   dir->Add(d);
-  d->Add(fEtaCorr);
   TIter    next(&fRingHistos);
   RingHistos* o = 0;
   while ((o = static_cast<RingHistos*>(next()))) {
index 71bf3cb77fdead72ea4bd4b253df796cca2e83b8..2c1bccf839f37819d40865c5b23980a0b3ec7c33 100644 (file)
@@ -71,6 +71,12 @@ public:
    * @param lowCut Low cut
    */
   void SetLowCut(Double_t lowCut=0.3) { fLowCut = lowCut; }
+  /** 
+   * Set the debug level.  The higher the value the more output 
+   * 
+   * @param dbg Debug level 
+   */
+  void SetDebug(Int_t dbg=1) { fDebug = dbg; }
 
   /** 
    * Enable use of angle corrected signals in the algorithm 
@@ -92,8 +98,7 @@ public:
    */
   Bool_t Filter(const AliESDFMD& input, 
                Bool_t           lowFlux, 
-               AliESDFMD&       output,
-               Double_t         vz);
+               AliESDFMD&       output);
   /** 
    * Scale the histograms to the total number of events 
    * 
@@ -244,7 +249,7 @@ protected:
   TList    fRingHistos;    // List of histogram containers
   Double_t fLowCut;        // Low cut on sharing
   Bool_t   fCorrectAngles; // Whether to work on angle corrected signals
-  TH2*     fEtaCorr;
+  Int_t    fDebug;         // Debug level 
 
   ClassDef(AliFMDSharingFilter,1); //
 };