More debug guard stuff for tracing execution.
authorcholm <cholm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 22 May 2012 13:57:03 +0000 (13:57 +0000)
committercholm <cholm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 22 May 2012 13:57:03 +0000 (13:57 +0000)
16 files changed:
PWGLF/FORWARD/analysis2/AliCentralMCMultiplicityTask.cxx
PWGLF/FORWARD/analysis2/AliCentralMultiplicityTask.cxx
PWGLF/FORWARD/analysis2/AliDisplacedVertexSelection.cxx
PWGLF/FORWARD/analysis2/AliFMDCorrector.cxx
PWGLF/FORWARD/analysis2/AliFMDDensityCalculator.cxx
PWGLF/FORWARD/analysis2/AliFMDEnergyFitter.cxx
PWGLF/FORWARD/analysis2/AliFMDEnergyFitterTask.cxx
PWGLF/FORWARD/analysis2/AliFMDEventInspector.cxx
PWGLF/FORWARD/analysis2/AliFMDEventPlaneFinder.cxx
PWGLF/FORWARD/analysis2/AliFMDEventPlaneTask.cxx
PWGLF/FORWARD/analysis2/AliFMDHistCollector.cxx
PWGLF/FORWARD/analysis2/AliFMDMCTrackDensity.cxx
PWGLF/FORWARD/analysis2/AliFMDSharingFilter.cxx
PWGLF/FORWARD/analysis2/AliForwardFlowTaskQC.cxx
PWGLF/FORWARD/analysis2/AliForwarddNdetaTask.cxx
PWGLF/FORWARD/analysis2/trains/MakeAODTrain.C

index 999cd01..612296c 100644 (file)
@@ -44,6 +44,7 @@ AliCentralMCMultiplicityTask::AliCentralMCMultiplicityTask(const char* name)
   // 
   // Constructor 
   //   
+  DGUARD(fDebug,0,"Named CTOR of AliCentralMCMultiplicityTask: %s", name);
 }
 //____________________________________________________________________
 AliCentralMCMultiplicityTask::AliCentralMCMultiplicityTask() 
@@ -54,6 +55,7 @@ AliCentralMCMultiplicityTask::AliCentralMCMultiplicityTask()
   // 
   // Constructor 
   // 
+  DGUARD(fDebug,0,"Default CTOR of AliCentralMCMultiplicityTask");
 }
 //____________________________________________________________________
 AliCentralMCMultiplicityTask::AliCentralMCMultiplicityTask(const AliCentralMCMultiplicityTask& o)
@@ -64,6 +66,7 @@ AliCentralMCMultiplicityTask::AliCentralMCMultiplicityTask(const AliCentralMCMul
   //
   // Copy constructor 
   // 
+  DGUARD(fDebug,0,"COPY CTOR of AliCentralMCMultiplicityTask");
 }
 //____________________________________________________________________
 AliCentralMCMultiplicityTask&
@@ -72,6 +75,8 @@ AliCentralMCMultiplicityTask::operator=(const AliCentralMCMultiplicityTask& o)
   // 
   // Assignment operator 
   //
+  DGUARD(fDebug,3,"Assignment of AliCentralMCMultiplicityTask");
+  if (&o == this) return *this; 
   AliCentralMultiplicityTask::operator=(o);
   fAODMCCentral     = o.fAODMCCentral;
   fTrackDensity     = o.fTrackDensity;
@@ -84,6 +89,7 @@ void AliCentralMCMultiplicityTask::UserCreateOutputObjects()
   // Create output objects 
   // 
   //
+  DGUARD(fDebug,1,"Create user output in AliCentralMCMultiplicityTask");
   AliCentralMultiplicityTask::UserCreateOutputObjects();
 
   AliAnalysisManager* am = AliAnalysisManager::GetAnalysisManager();
@@ -108,6 +114,7 @@ void AliCentralMCMultiplicityTask::UserExec(Option_t* option)
   // Parameters:
   //    option Not used
   //  
+  DGUARD(fDebug,1,"Process event in AliCentralMCMultiplicityTask");
   fAODMCCentral.Clear("");
   
   // Call base class 
@@ -143,6 +150,7 @@ void AliCentralMCMultiplicityTask::Terminate(Option_t* option)
   // Parameters:
   //    option Not used 
   //
+  DGUARD(fDebug,1,"Final analysis of merge in AliCentralMCMultiplicityTask");
   AliCentralMultiplicityTask::Terminate(option);
 }
 //____________________________________________________________________
index 14415a1..64fe763 100644 (file)
@@ -49,6 +49,7 @@ AliCentralMultiplicityTask::AliCentralMultiplicityTask(const char* name)
   // 
   // Constructor 
   //   
+  DGUARD(fDebug,0,"Named CTOR of AliCentralMultiplicityTask: %s", name);
   DefineOutput(1, TList::Class());
   fBranchNames = 
     "ESD:AliESDRun.,AliESDHeader.,AliMultiplicity.,"
@@ -76,6 +77,7 @@ AliCentralMultiplicityTask::AliCentralMultiplicityTask()
   // 
   // Constructor 
   // 
+  DGUARD(fDebug,0,"Default CTOR of AliCentralMultiplicityTask");
 }
 //____________________________________________________________________
 AliCentralMultiplicityTask::AliCentralMultiplicityTask(const AliCentralMultiplicityTask& o)
@@ -99,6 +101,8 @@ AliCentralMultiplicityTask::AliCentralMultiplicityTask(const AliCentralMultiplic
   //
   // Copy constructor 
   // 
+  DGUARD(fDebug,0,"COPY CTOR of AliCentralMultiplicityTask");
+
 }
 //____________________________________________________________________
 AliCentralMultiplicityTask&
@@ -107,6 +111,7 @@ AliCentralMultiplicityTask::operator=(const AliCentralMultiplicityTask& o)
   // 
   // Assignment operator 
   //
+  DGUARD(fDebug,3,"Assignment of AliCentralMultiplicityTask");
   if (&o == this) return *this; 
   fInspector         = o.fInspector;
   fData              = o.fData;
@@ -132,6 +137,7 @@ void AliCentralMultiplicityTask::UserCreateOutputObjects()
   // Create output objects 
   // 
   //
+  DGUARD(fDebug,1,"Create user output in AliCentralMultiplicityTask");
 
   AliAnalysisManager* am = AliAnalysisManager::GetAnalysisManager();
   AliAODHandler*      ah = 
@@ -159,6 +165,7 @@ AliCentralMultiplicityTask::GetESDEvent()
   //
   // Get the ESD event. IF this is the first event, initialise
   //
+  DGUARD(fDebug,1,"Get ESD event in AliCentralMultiplicityTask");
   AliESDEvent* esd = dynamic_cast<AliESDEvent*>(InputEvent());
   if (!esd) {
     AliWarning("No ESD event found for input event");
@@ -233,6 +240,7 @@ void
 AliCentralMultiplicityTask::MarkEventForStore() const
 {
   // Make sure the AOD tree is filled 
+  DGUARD(fDebug,1,"Mark AOD event for store in AliCentralMultiplicityTask");
   AliAnalysisManager* am = AliAnalysisManager::GetAnalysisManager();
   AliAODHandler*      ah = 
     dynamic_cast<AliAODHandler*>(am->GetOutputEventHandler());
@@ -245,6 +253,10 @@ AliCentralMultiplicityTask::MarkEventForStore() const
 //____________________________________________________________________
 void AliCentralMultiplicityTask::FindEtaLimits()
 {
+  // Find our pseudo-rapidity limits 
+  // 
+  // Uses the secondary map to do so.
+  DGUARD(fDebug,1,"Find eta limits in AliCentralMultiplicityTask");
   AliCentralCorrSecondaryMap* secMap = GetManager().GetSecMap();
 
   const TAxis& vaxis = secMap->GetVertexAxis();
@@ -261,10 +273,11 @@ void AliCentralMultiplicityTask::FindEtaLimits()
   secs->SetOwner();
   secs->SetName("secondaryMaps");
   fList->Add(secs);
+  unsigned short s = 1;
   TH2D* hCoverage = new TH2D("coverage", "#eta coverage per v_{z}", 
-                            secMap->GetCorrection(UShort_t(1))->GetXaxis()->GetNbins(),
-                            secMap->GetCorrection(UShort_t(1))->GetXaxis()->GetXmin(),
-                            secMap->GetCorrection(UShort_t(1))->GetXaxis()->GetXmax(),
+                            secMap->GetCorrection(s)->GetXaxis()->GetNbins(),
+                            secMap->GetCorrection(s)->GetXaxis()->GetXmin(),
+                            secMap->GetCorrection(s)->GetXaxis()->GetXmax(),
                             vaxis.GetNbins(),vaxis.GetXmin(),vaxis.GetXmax());
   hCoverage->SetDirectory(0);
   hCoverage->SetXTitle("#eta");
@@ -352,6 +365,7 @@ void AliCentralMultiplicityTask::UserExec(Option_t* /*option*/)
   // Parameters:
   //    option Not used
   //  
+  DGUARD(fDebug,1,"Process event in AliCentralMultiplicityTask");
   fAODCentral.Clear("");
   fIvz = 0;
 
@@ -401,6 +415,7 @@ void
 AliCentralMultiplicityTask::ProcessESD(TH2D& aodHist, 
                                       const AliMultiplicity* spdmult) const
 {
+  DGUARD(fDebug,1,"Process the ESD in AliCentralMultiplicityTask");
   fNTracklet->Reset();
   fNCluster->Reset();
 
@@ -432,6 +447,7 @@ void
 AliCentralMultiplicityTask::CorrectData(TH2D& aodHist, UShort_t vtxbin) const
 {  
   // Corrections
+  DGUARD(fDebug,1,"Correct data in AliCentralMultiplicityTask");
   TH1D* hAcceptance = fManager.GetAcceptanceCorrection(vtxbin);
   TH2D* hSecMap     = fManager.GetSecMapCorrection(vtxbin);
   
@@ -498,6 +514,7 @@ void AliCentralMultiplicityTask::Terminate(Option_t* /*option*/)
   // Parameters:
   //    option Not used 
   //
+  DGUARD(fDebug,1,"Process merged output in AliCentralMultiplicityTask");
 }
 //____________________________________________________________________
 void
index 03a14d3..e38a73f 100644 (file)
@@ -36,11 +36,13 @@ AliDisplacedVertexSelection::Output(TList* /*l*/, const char* /* name*/) const
 void
 AliDisplacedVertexSelection::Print(Option_t*) const
 {
+#if 0
   char ind[gROOT->GetDirLevel()+1];
   for (Int_t i = 0; i < gROOT->GetDirLevel(); i++) ind[i] = ' ';
   ind[gROOT->GetDirLevel()] = '\0';
   std::cout << std::boolalpha 
            << std::noboolalpha << std::endl;
+#endif
 }
 //____________________________________________________________________
 
@@ -107,11 +109,11 @@ for(Int_t i = 0; i < 4; ++i)
 //Global Event Filter//
 ///////////////////////
  Bool_t zdcAccSat[21];
- Bool_t zdcAccSatRunClass[21];
-for(Int_t k = -10; k <= 10; k++) 
+ // Bool_t zdcAccSatRunClass[21];
+ for(Int_t k = -10; k <= 10; k++) 
    {
       zdcAccSat[k+10]         = kFALSE;
-      zdcAccSatRunClass[k+10] = kFALSE;
+      // zdcAccSatRunClass[k+10] = kFALSE;
    }
 
 
@@ -220,11 +222,11 @@ for(Int_t i = 0; i < 4; ++i)
 //Global Event Filter//
 ///////////////////////
  Bool_t zdcAccSat[21];
- Bool_t zdcAccSatRunClass[21];
+ // Bool_t zdcAccSatRunClass[21];
 for(Int_t k = -10; k <= 10; k++) 
    {
       zdcAccSat[k+10]         = kFALSE;
-      zdcAccSatRunClass[k+10] = kFALSE;
+      // zdcAccSatRunClass[k+10] = kFALSE;
    }
 
 
index 637599b..10e4368 100644 (file)
@@ -35,6 +35,7 @@ AliFMDCorrector::AliFMDCorrector()
     fDebug(0)
 {
   // Constructor
+  DGUARD(fDebug, 0, "Default CTOR of AliFMDCorrector");
 }
 
 //____________________________________________________________________
@@ -51,6 +52,7 @@ AliFMDCorrector::AliFMDCorrector(const char* title)
   // 
   // Parameters: 
   //   title   Title
+  DGUARD(fDebug, 0, "Named CTOR of AliFMDCorrector: %s", title);
   fRingHistos.SetName(GetName());
   fRingHistos.Add(new RingHistos(1, 'I'));
   fRingHistos.Add(new RingHistos(2, 'I'));
@@ -73,6 +75,7 @@ AliFMDCorrector::AliFMDCorrector(const AliFMDCorrector& o)
   // 
   // Parameters: 
   //  o  Object to copy from 
+  DGUARD(fDebug, 0, "Copy CTOR of AliFMDCorrector");
   TIter    next(&o.fRingHistos);
   TObject* obj = 0;
   while ((obj = next())) fRingHistos.Add(obj);
@@ -84,6 +87,7 @@ AliFMDCorrector::~AliFMDCorrector()
   // Destructor 
   // 
   // 
+  DGUARD(fDebug, 3, "DTOR of AliFMDCorrector");
   fRingHistos.Delete();
 }
 
@@ -95,6 +99,7 @@ AliFMDCorrector::operator=(const AliFMDCorrector& o)
   // 
   // Parameters:
   //   o   Object to assign from 
+  DGUARD(fDebug, 3, "Assignment of AliFMDCorrector");
   if (&o == this) return *this; 
   TNamed::operator=(o);
 
@@ -121,6 +126,7 @@ AliFMDCorrector::Init(const TAxis&)
   // Parameters:
   //    etaAxis Eta axis to use
   //
+  DGUARD(fDebug, 1, "Initialization of AliFMDCorrector");
   if (!fUseSecondaryMap)
     AliWarning("Secondary maps not used - BE CAREFUL");
   if (!fUseVertexBias)
@@ -167,6 +173,7 @@ AliFMDCorrector::Correct(AliForwardUtil::Histos& hists,
   // Return:
   //    true on successs 
   //
+  DGUARD(fDebug, 1, "Correct histograms of AliFMDCorrector");
   AliForwardCorrectionManager& fcm = AliForwardCorrectionManager::Instance();
 
   UShort_t uvb = vtxbin;
@@ -299,6 +306,7 @@ AliFMDCorrector::ScaleHistograms(const TList* dir, Int_t nEvents)
   //    dir     Where the output is stored
   //    nEvents Number of events 
   //
+  DGUARD(fDebug, 1, "Scale histograms of AliFMDCorrector");
   if (nEvents <= 0) return;
   TList* d = static_cast<TList*>(dir->FindObject(GetName()));
   if (!d) return;
@@ -329,6 +337,7 @@ AliFMDCorrector::DefineOutput(TList* dir)
   // Parameters:
   //    dir List to write in
   //  
+  DGUARD(fDebug, 1, "Define output of AliFMDCorrector");
   TList* d = new TList;
   d->SetName(GetName());
   dir->Add(d);
index a1a77fb..690766b 100644 (file)
@@ -51,7 +51,7 @@ AliFMDDensityCalculator::AliFMDDensityCalculator()
   // 
   // Constructor 
   //
-   
+  DGUARD(fDebug, 0, "Default CTOR of FMD density calculator");
 }
 
 //____________________________________________________________________
@@ -85,6 +85,7 @@ AliFMDDensityCalculator::AliFMDDensityCalculator(const char* title)
   // Parameters:
   //    name Name of object
   //
+  DGUARD(fDebug, 0, "Named CTOR of FMD density calculator: %s", title);
   fRingHistos.SetName(GetName());
   fRingHistos.SetOwner();
   fRingHistos.Add(new RingHistos(1, 'I'));
@@ -151,6 +152,7 @@ AliFMDDensityCalculator::AliFMDDensityCalculator(const
   // Parameters:
   //    o Object to copy from 
   //
+  DGUARD(fDebug, 0, "Copy CTOR of FMD density calculator");
   TIter    next(&o.fRingHistos);
   TObject* obj = 0;
   while ((obj = next())) fRingHistos.Add(obj);
@@ -162,6 +164,7 @@ AliFMDDensityCalculator::~AliFMDDensityCalculator()
   // 
   // Destructor 
   //
+  DGUARD(fDebug, 3, "DTOR of FMD density calculator");
   fRingHistos.Delete();
 }
 
@@ -178,6 +181,7 @@ AliFMDDensityCalculator::operator=(const AliFMDDensityCalculator& o)
   // Return:
   //    Reference to this object
   //
+  DGUARD(fDebug, 3, "Assignment of FMD density calculator");
   if (&o == this) return *this; 
   TNamed::operator=(o);
 
@@ -215,6 +219,7 @@ AliFMDDensityCalculator::Init(const TAxis& axis)
   //
   // Parameters:
   //   etaAxis   Not used
+  DGUARD(fDebug, 1, "Initialize FMD density calculator");
   CacheMaxWeights();
 
   TIter    next(&fRingHistos);
@@ -306,8 +311,7 @@ AliFMDDensityCalculator::Calculate(const AliESDFMD&        fmd,
   // 
   // Return:
   //    true on successs 
-
-  
+  DGUARD(fDebug, 1, "Calculate density in FMD density calculator");
   
   for (UShort_t d=1; d<=3; d++) { 
     UShort_t nr = (d == 1 ? 1 : 2);
@@ -431,6 +435,7 @@ AliFMDDensityCalculator::FindMaxWeight(const AliFMDCorrELossFit* cor,
   //    r     Ring 
   //    iEta  Eta bin 
   //
+  DGUARD(fDebug, 3, "Find maximum weight in FMD density calculator");
   AliFMDCorrELossFit::ELossFit* fit = cor->GetFit(d,r,iEta);
   if (!fit) { 
     // AliWarning(Form("No energy loss fit for FMD%d%c at eta=%f", d, r, eta));
@@ -446,6 +451,7 @@ AliFMDDensityCalculator::CacheMaxWeights()
   // 
   // Find the max weights and cache them 
   // 
+  DGUARD(fDebug, 2, "Cache maximum weights in FMD density calculator");
   AliForwardCorrectionManager&  fcm = AliForwardCorrectionManager::Instance();
   AliFMDCorrELossFit*           cor = fcm.GetELossFit();
   const TAxis&                  eta = cor->GetEtaAxis();
@@ -583,6 +589,7 @@ AliFMDDensityCalculator::NParticles(Float_t  mult,
   //    The number of particles 
   //
   // if (mult <= GetMultCut()) return 0;
+  DGUARD(fDebug, 3, "Calculate Nch in FMD density calculator");
   if (lowFlux) return 1;
   
   AliForwardCorrectionManager&  fcm = AliForwardCorrectionManager::Instance();
@@ -640,6 +647,7 @@ AliFMDDensityCalculator::Correction(UShort_t d,
   // Return:
   //    
   //
+  DGUARD(fDebug, 3, "Apply correction in FMD density calculator");
   AliForwardCorrectionManager&  fcm = AliForwardCorrectionManager::Instance();
 
   Float_t correction = 1; 
@@ -674,6 +682,7 @@ AliFMDDensityCalculator::GenerateAcceptanceCorrection(Char_t r) const
   // Return:
   //    Newly allocated histogram of acceptance corrections
   //
+  DGUARD(fDebug, 3, "Make acceptance correction in FMD density calculator");
   const Double_t ic1[] = { 4.9895, 15.3560 };
   const Double_t ic2[] = { 1.8007, 17.2000 };
   const Double_t oc1[] = { 4.2231, 26.6638 };
@@ -768,6 +777,7 @@ AliFMDDensityCalculator::ScaleHistograms(const TList* dir, Int_t nEvents)
   //    dir     where to put the output
   //    nEvents Number of events 
   //
+  DGUARD(fDebug, 1, "Scale histograms in FMD density calculator");
   if (nEvents <= 0) return;
   TList* d = static_cast<TList*>(dir->FindObject(GetName()));
   if (!d) return;
@@ -798,6 +808,7 @@ AliFMDDensityCalculator::DefineOutput(TList* dir)
   // Parameters:
   //    dir List to write in
   //  
+  DGUARD(fDebug, 1, "Define output FMD density calculator");
   TList* d = new TList;
   d->SetOwner();
   d->SetName(GetName());
index 152784c..1a3c1ea 100644 (file)
@@ -49,6 +49,7 @@ AliFMDEnergyFitter::AliFMDEnergyFitter()
   // 
   // Default Constructor - do not use 
   //
+  DGUARD(fDebug, 0, "Default CTOR of AliFMDEnergyFitter");
 }
 
 //____________________________________________________________________
@@ -77,6 +78,7 @@ AliFMDEnergyFitter::AliFMDEnergyFitter(const char* title)
   // Parameters:
   //    title Title of object  - not significant 
   //
+  DGUARD(fDebug, 0, "Named CTOR of AliFMDEnergyFitter: %s", title);
   fEtaAxis.SetName("etaAxis");
   fEtaAxis.SetTitle("#eta");
   fCentralityAxis.SetName("centralityAxis");
@@ -114,6 +116,7 @@ AliFMDEnergyFitter::AliFMDEnergyFitter(const AliFMDEnergyFitter& o)
   // Parameters:
   //    o Object to copy from 
   //
+  DGUARD(fDebug, 0, "Copy CTOR of AliFMDEnergyFitter");
   TIter    next(&o.fRingHistos);
   TObject* obj = 0;
   while ((obj = next())) fRingHistos.Add(obj);
@@ -141,6 +144,7 @@ AliFMDEnergyFitter::operator=(const AliFMDEnergyFitter& o)
   // Return:
   //    Reference to this 
   //
+  DGUARD(fDebug, 3, "Assignment of AliFMDEnergyFitter");
   if (&o == this) return *this; 
   TNamed::operator=(o);
 
@@ -211,6 +215,7 @@ AliFMDEnergyFitter::Init(const TAxis& eAxis)
   // has already been set (using SetEtaAxis), then this parameter will be 
   // ignored
   //
+  DGUARD(fDebug, 1, "Initialize of AliFMDEnergyFitter");
   if (fEtaAxis.GetNbins() == 0 || 
       TMath::Abs(fEtaAxis.GetXmax() - fEtaAxis.GetXmin()) < 1e-7) 
     SetEtaAxis(eAxis);
@@ -284,6 +289,7 @@ AliFMDEnergyFitter::Accumulate(const AliESDFMD& input,
   // Return:
   //    True on success, false otherwise 
   //
+  DGUARD(fDebug, 1, "Accumulate statistics in AliFMDEnergyFitter");
   Int_t icent = fCentralityAxis.FindBin(cent);
   if (icent < 1 || icent > fCentralityAxis.GetNbins()) icent = 0;
 
@@ -329,6 +335,7 @@ AliFMDEnergyFitter::Fit(const TList* dir)
   // Parameters:
   //    dir Where the histograms are  
   //
+  DGUARD(fDebug, 1, "Fit distributions in AliFMDEnergyFitter");
   if (!fDoFits) return;
 
   TList* d = static_cast<TList*>(dir->FindObject(GetName()));
@@ -379,6 +386,7 @@ AliFMDEnergyFitter::MakeCorrectionsObject(TList* d)
   // Parameters:
   //    dir List to analyse 
   //
+  DGUARD(fDebug, 1, "Make the correction objec in AliFMDEnergyFitter");
   AliForwardCorrectionManager& mgr = AliForwardCorrectionManager::Instance();
     
   AliFMDCorrELossFit* obj = new AliFMDCorrELossFit;
@@ -412,6 +420,7 @@ AliFMDEnergyFitter::DefineOutput(TList* dir)
   // Parameters:
   //    dir Directory to add to 
   //
+  DGUARD(fDebug, 1, "Define output in AliFMDEnergyFitter");
   TList* d = new TList;
   d->SetName(GetName());
   dir->Add(d);
@@ -484,6 +493,7 @@ AliFMDEnergyFitter::RingHistos::RingHistos()
   // 
   // Default CTOR
   //
+  DGUARD(fDebug, 0, "Default CTOR AliFMDEnergyFitter::RingHistos");
 }
 
 //____________________________________________________________________
@@ -504,6 +514,8 @@ AliFMDEnergyFitter::RingHistos::RingHistos(UShort_t d, Char_t r)
   //    r ring 
   //
   fEtaEDists.SetName("EDists");
+  DGUARD(fDebug, 0, "Named CTOR AliFMDEnergyFitter::RingHistos: FMD%d%c",
+        d, r);
 }
 //____________________________________________________________________
 AliFMDEnergyFitter::RingHistos::RingHistos(const RingHistos& o)
@@ -521,6 +533,7 @@ AliFMDEnergyFitter::RingHistos::RingHistos(const RingHistos& o)
   // Parameters:
   //    o Object to copy from 
   //
+  DGUARD(fDebug, 0, "Copy CTOR AliFMDEnergyFitter::RingHistos");
   fFits.Clear();
   TIter next(&o.fEtaEDists);
   TObject* obj = 0;
@@ -546,6 +559,7 @@ AliFMDEnergyFitter::RingHistos::operator=(const RingHistos& o)
   // Return:
   //    Reference to this 
   //
+  DGUARD(fDebug, 3, "Assignment of AliFMDEnergyFitter::RingHistos");
   if (&o == this) return *this; 
   AliForwardUtil::RingHistos::operator=(o);
   
@@ -577,6 +591,7 @@ AliFMDEnergyFitter::RingHistos::~RingHistos()
   // 
   // Destructor 
   //
+  DGUARD(fDebug, 3, "DTOR of AliFMDEnergyFitter::RingHistos");
   if (fEDist) delete fEDist;
   fEtaEDists.Delete();
 }
@@ -595,6 +610,7 @@ AliFMDEnergyFitter::RingHistos::Fill(Bool_t empty, Int_t ieta,
   //    icent  Centrality bin (1-based)
   //    mult   Signal 
   //
+  DGUARD(fDebug, 2, "Filling in AliFMDEnergyFitter::RingHistos");
   if (empty) { 
     fEmpty->Fill(mult);
     return;
@@ -775,6 +791,7 @@ AliFMDEnergyFitter::RingHistos::Init(const TAxis& eAxis,
   //    nDEbins    Number of bins 
   //    useIncrBin Whether to use an increasing bin size 
   //
+  DGUARD(fDebug, 2, "Initialize in AliFMDEnergyFitter::RingHistos");
   fEDist = new TH1D(Form("%s_edist", fName.Data()), 
                    Form("#DeltaE/#DeltaE_{mip} for %s", fName.Data()), 
                    200, 0, 6);
@@ -873,6 +890,7 @@ AliFMDEnergyFitter::RingHistos::Fit(TList*           dir,
   //    chi2nuCut   Cut on @f$ \chi^2/\nu@f$ - 
   //                    the reduced @f$\chi^2@f$ 
   //
+  DGUARD(fDebug, 2, "Fit in AliFMDEnergyFitter::RingHistos");
 
   // Get our ring list from the output container 
   TList* l = GetOutputList(dir);
@@ -1112,6 +1130,8 @@ AliFMDEnergyFitter::RingHistos::FitHist(TH1*     dist,
   // Return:
   //    The best fit function 
   //
+  DGUARD(fDebug, 3, "Fit histogram in AliFMDEnergyFitter::RingHistos: %s",
+        dist->GetName());
   Double_t maxRange = 10;
 
   // Create a fitter object 
@@ -1357,6 +1377,7 @@ AliFMDEnergyFitter::RingHistos::Output(TList* dir)
   // Parameters:
   //    dir 
   //
+  DGUARD(fDebug, 2, "Define output in AliFMDEnergyFitter::RingHistos");
   fList = DefineOutputList(dir);
 }
 
index 5178a5d..aac5410 100644 (file)
@@ -44,6 +44,7 @@ AliFMDEnergyFitterTask::AliFMDEnergyFitterTask()
   // 
   // Constructor
   //
+  DGUARD(fDebug,0,"Default CTOR of AliFMDEnergyFitterTask");
 }
 
 //____________________________________________________________________
@@ -62,6 +63,7 @@ AliFMDEnergyFitterTask::AliFMDEnergyFitterTask(const char* name)
   // Parameters:
   //    name Name of task 
   //
+  DGUARD(fDebug,0,"Named CTOR of AliFMDEnergyFitterTask: %s", name);
   DefineOutput(1, TList::Class());
   DefineOutput(2, TList::Class());
 }
@@ -82,6 +84,7 @@ AliFMDEnergyFitterTask::AliFMDEnergyFitterTask(const AliFMDEnergyFitterTask& o)
   // Parameters:
   //    o Object to copy from 
   //
+  DGUARD(fDebug,0,"COPY CTOR of AliFMDEnergyFitterTask");
   DefineOutput(1, TList::Class());
   DefineOutput(2, TList::Class());
 }
@@ -99,6 +102,7 @@ AliFMDEnergyFitterTask::operator=(const AliFMDEnergyFitterTask& o)
   // Return:
   //    Reference to this object 
   //
+  DGUARD(fDebug,3,"Assignment of AliFMDEnergyFitterTask");
   if (&o == this) return *this; 
   AliAnalysisTaskSE::operator=(o);
 
@@ -134,6 +138,7 @@ AliFMDEnergyFitterTask::Init()
   // Initialize the task 
   // 
   //
+  DGUARD(fDebug,1,"Initialize of AliFMDEnergyFitterTask");
   fFirstEvent = true;
 }
 
@@ -145,6 +150,7 @@ AliFMDEnergyFitterTask::InitializeSubs()
   // Initialise the sub objects and stuff.  Called on first event 
   // 
   //
+  DGUARD(fDebug,1,"Initialize subs of AliFMDEnergyFitterTask");
   AliForwardCorrectionManager& fcm = AliForwardCorrectionManager::Instance();
   UShort_t sys = fEventInspector.GetCollisionSystem();
   UShort_t sNN = fEventInspector.GetEnergy();
@@ -165,6 +171,7 @@ AliFMDEnergyFitterTask::UserCreateOutputObjects()
   // Create output objects 
   // 
   //
+  DGUARD(fDebug,1,"Create output objects of AliFMDEnergyFitterTask");
   fList = new TList;
   fList->SetOwner();
 
@@ -187,6 +194,7 @@ AliFMDEnergyFitterTask::UserExec(Option_t*)
   // static Int_t cnt = 0;
   // cnt++;
   // Get the input data 
+  DGUARD(fDebug,1,"Analyse event of AliFMDEnergyFitterTask");
   
   AliMCEvent* mcevent = MCEvent();
   if(mcevent) {
@@ -275,7 +283,8 @@ AliFMDEnergyFitterTask::Terminate(Option_t*)
   // Parameters:
   //    option Not used 
   //
-  AliInfo(Form("Running terminate of %s", GetName()));
+  DGUARD(fDebug,1,"Processing merged output of AliFMDEnergyFitterTask");
+
   TList* list = dynamic_cast<TList*>(GetOutputData(1));
   if (!list) {
     AliError(Form("No output list defined (%p)", GetOutputData(1)));
index 2fb78fa..a544d70 100644 (file)
@@ -70,6 +70,7 @@ AliFMDEventInspector::AliFMDEventInspector()
   // 
   // Constructor 
   //
+  DGUARD(fDebug,0,"Default CTOR of AliFMDEventInspector");
 }
 
 //____________________________________________________________________
@@ -106,6 +107,7 @@ AliFMDEventInspector::AliFMDEventInspector(const char* name)
   // Parameters:
   //   name Name of object
   //
+  DGUARD(fDebug,0,"Named CTOR of AliFMDEventInspector: %s", name);
 }
 
 //____________________________________________________________________
@@ -142,6 +144,7 @@ AliFMDEventInspector::AliFMDEventInspector(const AliFMDEventInspector& o)
   // Parameters:
   //   o Object to copy from 
   //
+  DGUARD(fDebug,0,"Copy CTOR of AliFMDEventInspector");
 }
 
 //____________________________________________________________________
@@ -165,6 +168,7 @@ AliFMDEventInspector::operator=(const AliFMDEventInspector& o)
   // Return:
   //    Reference to this object
   //
+  DGUARD(fDebug,3,"Assignment of AliFMDEventInspector");
   if (&o == this) return *this; 
   TNamed::operator=(o);
   fHEventsTr         = o.fHEventsTr;
@@ -224,6 +228,7 @@ AliFMDEventInspector::FetchHistograms(const TList* d,
   // Return:
   //    true on success, false otherwise 
   //
+  DGUARD(fDebug,3,"Fetch histograms in AliFMDEventInspector");
   hEventsTr    = 0;
   hEventsTrVtx = 0;
   hTriggers    = 0;
@@ -247,7 +252,7 @@ AliFMDEventInspector::Init(const TAxis& vtxAxis)
   // Parameters:
   //   vtxAxis Vertex axis in use 
   //
-  AliInfo("Initialising the event inspector");
+  DGUARD(fDebug,1,"Initialize in AliFMDEventInspector");
   
   // -1.5 -0.5 0.5 1.5 ... 89.5 ... 100.5
   // ----- 92 number --------- ---- 1 ---
@@ -369,6 +374,7 @@ AliFMDEventInspector::StoreInformation(Int_t runNo)
 {
   // Write TNamed objects to output list containing information about
   // the running conditions 
+  DGUARD(fDebug,2,"Store information from AliFMDEventInspector");
   if (!fList) return;
 
   TNamed* sys = new TNamed("sys", "");
@@ -416,6 +422,7 @@ AliFMDEventInspector::DefineOutput(TList* dir)
   // 
   //   dir Directory to add to 
   //
+  DGUARD(fDebug,1,"Define output from AliFMDEventInspector");
   fList = new TList;
   fList->SetName(GetName());
   dir->Add(fList);
@@ -447,6 +454,7 @@ AliFMDEventInspector::Process(const AliESDEvent* event,
   // Return:
   //    0 (or kOk) on success, otherwise a bit mask of error codes 
   //
+  DGUARD(fDebug,1,"Process event in AliFMDEventInspector");
 
   // --- Check that we have an event ---------------------------------
   if (!event) { 
@@ -542,6 +550,8 @@ AliFMDEventInspector::ReadCentrality(const AliESDEvent* esd,
   // Return:
   //    False on error, true otherwise 
   //
+  DGUARD(fDebug,2,"Read the centrality in AliFMDEventInspector");
+  
   cent = -1;
   qual = 0;
   AliCentrality* centObj = const_cast<AliESDEvent*>(esd)->GetCentrality();
@@ -613,6 +623,7 @@ AliFMDEventInspector::ReadTriggers(const AliESDEvent* esd, UInt_t& triggers,
   // Return:
   //    @c true on success, @c false otherwise 
   //
+  DGUARD(fDebug,2,"Read the triggers in AliFMDEventInspector");
   triggers = 0;
 
   // Get the analysis manager - should always be there 
@@ -925,6 +936,7 @@ AliFMDEventInspector::ReadVertex(const AliESDEvent* esd,
   // Return:
   //    @c true on success, @c false otherwise 
   //
+  DGUARD(fDebug,2,"Read the vertex in AliFMDEventInspector");
   vz = 0;
   vx = 1024;
   vy = 1024;
@@ -1057,6 +1069,7 @@ AliFMDEventInspector::ReadRunDetails(const AliESDEvent* esd)
   // AliInfo(Form("Parameters from 1st ESD event: cms=%s, sNN=%f, field=%f",
   //          esd->GetBeamType(), 2*esd->GetBeamEnergy(), 
   //          esd->GetMagneticField()));
+  DGUARD(fDebug,2,"Read the run details in AliFMDEventInspector");
   fCollisionSystem = 
     AliForwardUtil::ParseCollisionSystem(esd->GetBeamType());
   fEnergy          = 
index 84f3e07..bb0eed9 100644 (file)
@@ -57,6 +57,7 @@ AliFMDEventPlaneFinder::AliFMDEventPlaneFinder()
   // 
   // Constructor 
   //
+  DGUARD(fDebug,0,"Default CTOR of AliFMDEventPlaneFinder");
    
 }
 
@@ -94,6 +95,7 @@ AliFMDEventPlaneFinder::AliFMDEventPlaneFinder(const char* title)
   // Parameters:
   //    name Name of object
   //
+  DGUARD(fDebug,0,"Named CTOR of AliFMDEventPlaneFinder: %s", title);
 }
 
 //____________________________________________________________________
@@ -131,6 +133,7 @@ AliFMDEventPlaneFinder::AliFMDEventPlaneFinder(const
   // Parameters:
   //    o Object to copy from 
   //
+  DGUARD(fDebug,0,"Copy CTOR of AliFMDEventPlaneFinder");
 }
 
 //____________________________________________________________________
@@ -154,6 +157,7 @@ AliFMDEventPlaneFinder::operator=(const AliFMDEventPlaneFinder& o)
   // Return:
   //    Reference to this object
   //
+  DGUARD(fDebug,3,"Assignment of AliFMDEventPlaneFinder");
   if (&o == this) return *this; 
   TNamed::operator=(o);
 
@@ -194,6 +198,7 @@ AliFMDEventPlaneFinder::Init(const TAxis& etaAxis)
   // Parameters:
   //   etaAxis   fmd eta axis binning
   //
+  DGUARD(fDebug,1,"Initalization of AliFMDEventPlaneFinder");
   fHistFMDEventplane = new TH1D("hFMDEventplane", "FMD eventplane", 
                                 100, 0., TMath::Pi());
   fHistFMDEventplane->Sumw2();
@@ -298,6 +303,7 @@ AliFMDEventPlaneFinder::DefineOutput(TList* dir)
   // Parameters:
   //    dir List to write in
   //  
+  DGUARD(fDebug,1,"Define output of AliFMDEventPlaneFinder");
   fList = new TList();
   fList->SetOwner();
   fList->SetName(GetName());
@@ -322,6 +328,7 @@ AliFMDEventPlaneFinder::FindEventplane(AliVEvent* ev,
   // 
   // Return:
   //    true on successs 
+  DGUARD(fDebug,1,"Find the event plane in AliFMDEventPlaneFinder");
 
   fQt.Set(0., 0.);
   fQa.Set(0., 0.);
@@ -363,6 +370,7 @@ AliFMDEventPlaneFinder::CalcQVectors(TH2D* h, TH1D* eHist)
   //
   // Calculate the Q vectors
   //
+  DGUARD(fDebug,2,"Calculate Q-vectors in AliFMDEventPlaneFinder");
   Double_t phi = 0, eta = 0, weight = 0;
   for (Int_t e = 1; e <= h->GetNbinsX(); e++) {
     Double_t qx = 0, qy = 0;
@@ -406,6 +414,7 @@ AliFMDEventPlaneFinder::CalcEventplane(TVector2 v) const
   //
   // Calculate the eventplane
   //
+  DGUARD(fDebug,2,"Calculate Event plane in AliFMDEventPlaneFinder");
   Double_t ep = -1;
  
   if (v.Mod() == 0.) return ep;
@@ -426,6 +435,7 @@ AliFMDEventPlaneFinder::FillHists(AliAODForwardEP* fmdEP)
   //
   // Fill diagnostics histograms
   //
+  DGUARD(fDebug,2,"Fill histograms in AliFMDEventPlaneFinder");
   Double_t fmdEPt = fmdEP->GetEventplane();
   Double_t fmdEPa = fmdEP->GetEventplaneA();
   Double_t fmdEPc = fmdEP->GetEventplaneC();
index a1a503a..d07beb3 100644 (file)
@@ -41,6 +41,7 @@ AliFMDEventPlaneTask::AliFMDEventPlaneTask()
   // 
   // Default constructor
   //
+  DGUARD(fDebug,0,"Default CTOR of AliFMDEventPlaneTask");
 }
 //_____________________________________________________________________
 AliFMDEventPlaneTask::AliFMDEventPlaneTask(const char* name) 
@@ -63,6 +64,7 @@ AliFMDEventPlaneTask::AliFMDEventPlaneTask(const char* name)
   // Parameters:
   //  name: Name of task
   //
+  DGUARD(fDebug,0,"Named CTOR of AliFMDEventPlaneTask: %s", name);
 
   DefineOutput(1, TList::Class());
   DefineOutput(2, TList::Class());
@@ -88,6 +90,7 @@ AliFMDEventPlaneTask::AliFMDEventPlaneTask(const AliFMDEventPlaneTask& o)
   // Parameters:
   //    o Object to copy from 
   //
+  DGUARD(fDebug,0,"Copy CTOR of AliFMDEventPlaneTask");
 }
 //_____________________________________________________________________
 AliFMDEventPlaneTask&
@@ -96,6 +99,7 @@ AliFMDEventPlaneTask::operator=(const AliFMDEventPlaneTask& o)
   // 
   // Assignment operator 
   //
+  DGUARD(fDebug,3,"Assignment of AliFMDEventPlaneTask");
   if (&o == this) return *this;
   fSumList           = o.fSumList;
   fOutputList        = o.fOutputList;
@@ -116,6 +120,7 @@ void AliFMDEventPlaneTask::UserCreateOutputObjects()
   //
   // Create output objects
   //
+  DGUARD(fDebug,1,"Create user objects of AliFMDEventPlaneTask");
   if (!fSumList)
     fSumList = new TList();
   fSumList->SetName("Sums");
@@ -147,6 +152,7 @@ void AliFMDEventPlaneTask::UserExec(Option_t */*option*/)
   // Parameters:
   //  option: Not used
   //
+  DGUARD(fDebug,1,"Process an event in AliFMDEventPlaneTask");
 
   // Reset data members
   fCent = -1;
@@ -181,6 +187,7 @@ void AliFMDEventPlaneTask::Terminate(Option_t */*option*/)
   // Parameters:
   //  option: Not used
   //
+  DGUARD(fDebug,1,"Process merged output of AliFMDEventPlaneTask");
 
   // Reinitiate lists if Terminate is called separately!
   fSumList = dynamic_cast<TList*> (GetOutputData(1));
@@ -208,6 +215,7 @@ Bool_t AliFMDEventPlaneTask::AODCheck(const AliAODForwardMult* aodfm)
   // Parameters: 
   //  AliAODForwardMult: forward mult object with trigger and vertex info
   //
+  DGUARD(fDebug,2,"Check AOD in AliFMDEventPlaneTask");
 
   if (!aodfm->IsTriggerBits(AliAODForwardMult::kOffline)) return kFALSE;
 
index 9ab4e50..a2ac1a7 100644 (file)
@@ -47,7 +47,9 @@ AliFMDHistCollector::AliFMDHistCollector()
     fFiducialMethod(kByCut),
     fSkipFMDRings(0),
     fBgAndHitMaps(false)
-{}
+{
+  DGUARD(fDebug, 0, "Default CTOR of AliFMDHistCollector");
+}
 
 //____________________________________________________________________
 AliFMDHistCollector::AliFMDHistCollector(const char* title)
@@ -65,6 +67,7 @@ AliFMDHistCollector::AliFMDHistCollector(const char* title)
     fSkipFMDRings(0),
     fBgAndHitMaps(false)
 {
+  DGUARD(fDebug, 0, "Named CTOR of AliFMDHistCollector: %s", title);
 }
 //____________________________________________________________________
 AliFMDHistCollector::AliFMDHistCollector(const AliFMDHistCollector& o)
@@ -81,11 +84,14 @@ AliFMDHistCollector::AliFMDHistCollector(const AliFMDHistCollector& o)
     fFiducialMethod(o.fFiducialMethod),
     fSkipFMDRings(o.fSkipFMDRings),
     fBgAndHitMaps(o.fBgAndHitMaps)
-{}
+{
+  DGUARD(fDebug, 0, "Copy CTOR of AliFMDHistCollector");
+}
 
 //____________________________________________________________________
 AliFMDHistCollector::~AliFMDHistCollector()
 { 
+  DGUARD(fDebug, 3, "DTOR of AliFMDHistCollector");
   if (fList) delete fList;
 }
 //____________________________________________________________________
@@ -101,6 +107,7 @@ AliFMDHistCollector::operator=(const AliFMDHistCollector& o)
   // Return:
   //    Reference to this object
   //
+  DGUARD(fDebug, 3, "Assignment of AliFMDHistCollector");
   if (&o == this) return *this; 
   TNamed::operator=(o);
 
@@ -131,6 +138,7 @@ AliFMDHistCollector::Init(const TAxis& vtxAxis,
   // Parameters:
   //    vtxAxis  Vertex axis 
   //  
+  DGUARD(fDebug, 1, "Initialization of AliFMDHistCollector");
 
   AliForwardCorrectionManager& fcm = AliForwardCorrectionManager::Instance();
 
@@ -283,6 +291,7 @@ AliFMDHistCollector::DefineOutput(TList* dir)
   // Parameters:
   //    dir List to write in
   //  
+  DGUARD(fDebug, 1, "Define output of AliFMDHistCollector");
   fList = new TList;
   fList->SetOwner();
   fList->SetName(GetName());
@@ -570,6 +579,7 @@ AliFMDHistCollector::Collect(const AliForwardUtil::Histos& hists,
   // Return:
   //    true on successs 
   //
+  DGUARD(fDebug, 1, "Collect final histogram of AliFMDHistCollector");
   AliForwardCorrectionManager& fcm = AliForwardCorrectionManager::Instance();
   const TAxis* vtxAxis = fcm.GetVertexAxis();
   Double_t vMin    = vtxAxis->GetBinLowEdge(vtxbin);
@@ -666,7 +676,8 @@ AliFMDHistCollector::Collect(const AliForwardUtil::Histos& hists,
       // Remove temporary histogram 
       if(fBgAndHitMaps) {
        TH2D* hRingSumVtx 
-         = static_cast<TH2D*>(vtxList->FindObject(Form("hitMapFMD%d%c", d, r)));
+         = static_cast<TH2D*>(vtxList->FindObject(Form("hitMapFMD%d%c", 
+                                                       d, r)));
        hRingSumVtx->Add(t);
       }
       delete t;
index a5cb9b5..4418616 100644 (file)
@@ -30,6 +30,7 @@ AliFMDMCTrackDensity::State::Clear(Bool_t alsoCount)
 AliFMDMCTrackDensity::State&
 AliFMDMCTrackDensity::State::operator=(const State& o)
 {
+  if (&o == this) return *this;
   angle          = o.angle;
   oldDetector    = o.oldDetector;
   oldRing        = o.oldRing;
index 2e451d3..d30be27 100644 (file)
@@ -69,6 +69,7 @@ AliFMDSharingFilter::AliFMDSharingFilter()
   // 
   // Default Constructor - do not use 
   //
+  DGUARD(fDebug,0, "Default CTOR for AliFMDSharingFilter");
 }
 
 //____________________________________________________________________
@@ -94,6 +95,7 @@ AliFMDSharingFilter::AliFMDSharingFilter(const char* title)
   // Parameters:
   //    title Title of object  - not significant 
   //
+  DGUARD(fDebug,0, "Named CTOR for AliFMDSharingFilter: %s", title);
   fRingHistos.SetName(GetName());
   fRingHistos.SetOwner();
   
@@ -131,6 +133,7 @@ AliFMDSharingFilter::AliFMDSharingFilter(const AliFMDSharingFilter& o)
   // Parameters:
   //    o Object to copy from 
   //
+  DGUARD(fDebug,0, "Copy CTOR for AliFMDSharingFilter");
   TIter    next(&o.fRingHistos);
   TObject* obj = 0;
   while ((obj = next())) fRingHistos.Add(obj);
@@ -142,6 +145,7 @@ AliFMDSharingFilter::~AliFMDSharingFilter()
   // 
   // Destructor
   //
+  DGUARD(fDebug,3, "DTOR for AliFMDSharingFilter");
   fRingHistos.Delete();
 }
 
@@ -158,6 +162,7 @@ AliFMDSharingFilter::operator=(const AliFMDSharingFilter& o)
   // Return:
   //    Reference to this 
   //
+  DGUARD(fDebug,3, "Assigment for AliFMDSharingFilter");
   if (&o == this) return *this;
   TNamed::operator=(o);
 
@@ -212,6 +217,7 @@ void
 AliFMDSharingFilter::Init()
 {
   // Initialise 
+  DGUARD(fDebug,1, "Initialize for AliFMDSharingFilter");
   AliForwardCorrectionManager&  fcm = AliForwardCorrectionManager::Instance();
  
   // Get the high cut.  The high cut is defined as the 
@@ -270,6 +276,7 @@ AliFMDSharingFilter::Filter(const AliESDFMD& input,
   // Return:
   //    True on success, false otherwise 
   //
+  DGUARD(fDebug,1, "Filter event in AliFMDSharingFilter");
   output.Clear();
   TIter    next(&fRingHistos);
   RingHistos* o      = 0;
@@ -493,7 +500,7 @@ AliFMDSharingFilter::Filter(const AliESDFMD& input,
   fSummed->Fill(kMergedWithOther, nMerged);
   fSummed->Fill(kMergedInto,      nSummed);
 
-  DBGL(1, Form("none=%9d, candidate=%9d, merged=%9d, summed=%9d", 
+  DBGL(5, Form("none=%9d, candidate=%9d, merged=%9d, summed=%9d", 
               nNone, nCandidate, nMerged, nSummed));
   next.Reset();
   while ((o = static_cast<RingHistos*>(next())))
@@ -620,14 +627,14 @@ AliFMDSharingFilter::MultiplicityOfStrip(Double_t thisE,
   // strip is considered a candidate.
   if (thisE < lowCut || thisE > 20) { 
     thisStatus = kNone;
-    DBGL(3,Form(" %9f<%9f || %9f>20, 0'ed", thisE, lowCut, thisE));
+    DBGL(5,Form(" %9f<%9f || %9f>20, 0'ed", thisE, lowCut, thisE));
     if (prevStatus == kCandidate) prevStatus = kNone;
     return 0;  
   }
   // It this strip was merged with the previous strip, then 
   // make the next strip a candidate and zero the value in this strip. 
   if (thisStatus == kMergedWithOther) { 
-    DBGL(3,Form(" Merged with other, 0'ed"));
+    DBGL(5,Form(" Merged with other, 0'ed"));
     return 0;
   }
 
@@ -675,7 +682,7 @@ AliFMDSharingFilter::MultiplicityOfStrip(Double_t thisE,
       thisStatus = kMergedInto;
     else 
       thisStatus = kNone;
-    DBGL(3, Form(" %9f>%f9 (was %9f) -> %9f %s", 
+    DBGL(5, Form(" %9f>%f9 (was %9f) -> %9f %s", 
                 totalE, highCut, thisE, totalE,status2String(thisStatus)));
     return totalE;
   }
@@ -687,14 +694,14 @@ AliFMDSharingFilter::MultiplicityOfStrip(Double_t thisE,
   
   // If the total signal is smaller than low cut then zero this and kill this 
   if (totalE < lowCut)  {
-    DBGL(3, Form(" %9f<%9f (was %9f), zeroed", totalE, lowCut, thisE));
+    DBGL(5, Form(" %9f<%9f (was %9f), zeroed", totalE, lowCut, thisE));
     thisStatus = kNone;
     return 0;
   }
 
   // If total signal not above high cut or lower than low cut, 
   // mark this as a candidate for merging into the next, and zero signal 
-  DBGL(3, Form(" %9f<%9f<%9f (was %9f), zeroed, candidate", 
+  DBGL(5, Form(" %9f<%9f<%9f (was %9f), zeroed, candidate", 
                   lowCut, totalE, highCut, thisE));
   thisStatus = kCandidate;
   return (fZeroSharedHitsBelowThreshold ? 0 : totalE); 
@@ -868,6 +875,7 @@ AliFMDSharingFilter::ScaleHistograms(const TList* dir, Int_t nEvents)
   //    dir     Where the output is 
   //    nEvents Number of events 
   //
+  DGUARD(fDebug,1, "Scale histograms in AliFMDSharingFilter");
   if (nEvents <= 0) return;
   TList* d = static_cast<TList*>(dir->FindObject(GetName()));
   if (!d) return;
@@ -899,6 +907,7 @@ AliFMDSharingFilter::DefineOutput(TList* dir)
   // Parameters:
   //    dir Directory to add to 
   //
+  DGUARD(fDebug,1, "Define output in AliFMDSharingFilter");
   TList* d = new TList;
   d->SetName(GetName());
   dir->Add(d);
index 489d57c..42daec4 100644 (file)
@@ -526,13 +526,13 @@ void AliForwardFlowTaskQC::VertexBin::AddOutput(TList* outputlist)
     Int_t nBins = accMap->GetNbinsX();
     Double_t eta[48] = { 0. };
     Int_t n = 0;
-    Double_t newOcc[48] = { 0. };
+    // Double_t newOcc[48] = { 0. };
     Double_t prev = -1;
     for (Int_t i = 0; i < nBins; i++) {
       Double_t occ = accMap->GetBinContent(i+1);
       if (prev != occ && (((occ > 0.6 || occ == 0) && i*0.25-6 < 4) || ((occ == 0) && i*0.25-6 >= 4))) {
         eta[n] = i*0.25-6.;
-        newOcc[n] = occ;
+        // newOcc[n] = occ;
         n++;
         if (fDebug > 5) AliInfo(Form("eta: %f \t occ: %f \t Vertex: %d \n", eta[n-1], occ, fVzMin));
       }
@@ -896,7 +896,7 @@ void AliForwardFlowTaskQC::VertexBin::CumulantsTerminate(TList* inlist, TList* o
   }
 
   // For flow calculations
-  Double_t two = 0, qc2 = 0, vnTwo = 0, four = 0, qc4 = 0, vnFour = 0; 
+  Double_t two = 0, qc2 = 0, /*vnTwo = 0,*/ four = 0, qc4 = 0/*, vnFour = 0*/; 
   Double_t twoPrime = 0, qc2Prime = 0, vnTwoDiff = 0, fourPrime = 0, qc4Prime = 0, vnFourDiff = 0;
   Double_t w2 = 0, w4 = 0, w2p = 0, w4p = 0;
   Double_t w2Two = 0, w2pTwoPrime = 0, w4Four = 0, w4pFourPrime = 0;
@@ -931,7 +931,7 @@ void AliForwardFlowTaskQC::VertexBin::CumulantsTerminate(TList* inlist, TList* o
        if (fDebug > 0) AliInfo(Form("%s: QC_%d{2} = %1.3f for eta = %1.2f and centrality %3.1f - skipping", fType.Data(), fMoment, qc2, eta, cent));
        continue;
       }
-      vnTwo = TMath::Sqrt(qc2);
+      // vnTwo = TMath::Sqrt(qc2);
  //     if (!TMath::IsNaN(vnTwo*mult)) 
  //       cumu2->Fill(eta, cent, vnTwo, fCumuHist->GetBinContent(0,cBin,0)); 
 
@@ -979,7 +979,7 @@ void AliForwardFlowTaskQC::VertexBin::CumulantsTerminate(TList* inlist, TList* o
        if (fDebug > 0) AliInfo(Form("%s: QC_%d{4} = %1.3f for eta = %1.2f and centrality %3.1f - skipping", fType.Data(), fMoment, qc2, eta, cent));
        continue;
       }
-      vnFour = TMath::Power(-qc4, 0.25);
+      // vnFour = TMath::Power(-qc4, 0.25);
  //     if (!TMath::IsNaN(vnFour*mult)) 
  //         cumu4->Fill(eta, cent, vnFour, fCumuHist->GetBinContent(0,cBin,0));
 
index d2cc33f..3dcb3bf 100644 (file)
@@ -20,6 +20,7 @@ AliForwarddNdetaTask::AliForwarddNdetaTask()
   //
   // Constructor 
   // 
+  DGUARD(fDebug, 0, "Default CTOR of AliForwarddNdetaTask");
 }
 
 //____________________________________________________________________
@@ -32,6 +33,7 @@ AliForwarddNdetaTask::AliForwarddNdetaTask(const char* /* name */)
   // Paramters
   //   name    Name of task 
   SetTitle("FMD");
+  DGUARD(fDebug, 0, "Named CTOR of AliForwarddNdetaTask");
 }
 
 //____________________________________________________________________
@@ -41,6 +43,7 @@ AliForwarddNdetaTask::AliForwarddNdetaTask(const AliForwarddNdetaTask& o)
   // 
   // Copy constructor
   // 
+  DGUARD(fDebug, 0, "Copy CTOR of AliForwarddNdetaTask");
 }
 
 //____________________________________________________________________
@@ -59,6 +62,8 @@ AliForwarddNdetaTask::MakeCentralityBin(const char* name, Short_t l, Short_t h)
   // Return:
   //    Newly allocated object (of our type)
   //
+  DGUARD(fDebug, 3,"Make a centrality bin for AliForwarddNdetaTask: %s [%d,%d]",
+        name, l, h);
   return new AliForwarddNdetaTask::CentralityBin(name, l, h);
 }
 
@@ -108,13 +113,14 @@ AliForwarddNdetaTask::CentralityBin::End(TList*      sums,
                                         TList*      mclist,
                                         TList*      truthlist )
 {
-  AliInfo(Form("In End of %s with corrEmpty=%d, cutEdges=%d, rootProj=%d", 
-              GetName(), corrEmpty, cutEdges, rootProj));
+  DGUARD(fDebug, 1, "In End of %s with corrEmpty=%d, cutEdges=%d, rootProj=%d", 
+        GetName(), corrEmpty, cutEdges, rootProj);
   AliBasedNdetaTask::CentralityBin::End(sums, results, scheme, 
                                        shapeCorr, trigEff, 
                                        symmetrice, rebin, 
                                        rootProj, corrEmpty, cutEdges,
-                                       triggerMask, marker, color, mclist, truthlist);
+                                       triggerMask, marker, color, mclist, 
+                                       truthlist);
 
   if (!IsAllBin()) return;
   TFile* file = TFile::Open("forward.root", "READ");
index 9f0c7d0..268fb59 100644 (file)
@@ -75,13 +75,16 @@ public:
    * @param mc       If true, assume simulated events 
    * @param usePar   If true, use PARs 
    */
-  void Run(const char* mode, const char* oper, 
-          Int_t nEvents=-1, Bool_t mc=false,
-          Bool_t usePar=false)
+  void Run(const char* mode, 
+          const char* oper, 
+          Int_t       nEvents=-1, 
+          Bool_t      mc=false,
+          Bool_t      usePar=false, 
+          Int_t       dbg=0)
   {
     Info("Run", "Running in mode=%s, oper=%s, events=%d, MC=%d, Par=%d", 
         mode, oper, nEvents, mc, usePar);
-    Exec("ESD", mode, oper, nEvents, mc, usePar);
+    Exec("ESD", mode, oper, nEvents, mc, usePar, dbg);
   }
   /** 
    * Run this analysis 
@@ -92,12 +95,16 @@ public:
    * @param mc       If true, assume simulated events 
    * @param usePar   If true, use PARs 
    */
-  void Run(EMode mode, EOper oper, Int_t nEvents=-1, Bool_t mc=false, 
-          Bool_t usePar = false)
+  void Run(EMode  mode, 
+          EOper  oper, 
+          Int_t  nEvents=-1, 
+          Bool_t mc=false, 
+          Bool_t usePar = false,
+          Int_t  dbg = 0)
   {
     Info("Run", "Running in mode=%d, oper=%d, events=%d, MC=%d, Par=%d", 
         mode, oper, nEvents, mc, usePar);
-    Exec(kESD, mode, oper, nEvents, mc, usePar);
+    Exec(kESD, mode, oper, nEvents, mc, usePar, dbg);
   }
 protected:
   /**