added storage of 1D hits histograms for diagnostics

diff --git a/FMD/analysis/AliFMDAnalysisTaskBackgroundCorrection.cxx b/FMD/analysis/AliFMDAnalysisTaskBackgroundCorrection.cxx
index ee9f95646b879eb803ebccfe21d36f0ce717b898..b9fb17632de1d6dd4366a2f8619f49ff14595ea7 100644 (file)
--- a/FMD/analysis/AliFMDAnalysisTaskBackgroundCorrection.cxx
+++ b/FMD/analysis/AliFMDAnalysisTaskBackgroundCorrection.cxx
@@ -189,8 +189,22 @@ void AliFMDAnalysisTaskBackgroundCorrection::Exec(Option_t */*option*/)
 }
 //_____________________________________________________________________
 void AliFMDAnalysisTaskBackgroundCorrection::Terminate(Option_t */*option*/) {
+  AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();
   
+  Int_t nVtxbins = pars->GetNvtxBins();
   
+  for(UShort_t det=1;det<=3;det++) {
+    Int_t nRings = (det==1 ? 1 : 2);
+    for (UShort_t ir = 0; ir < nRings; ir++) {
+      Char_t ringChar = (ir == 0 ? 'I' : 'O');
+      for(Int_t i =0; i<nVtxbins; i++) {
+       TH2F* hHits      = (TH2F*)fOutputList->FindObject(Form("hits_FMD%d%c_vtxbin%d",det,ringChar,i));
+       TH1D* hHitsproj  = hHits->ProjectionX(Form("hits_FMD%d%c_vtxbin%d_proj",det,ringChar,i),1,hHits->GetNbinsY());
+       fHitList->Add(hHitsproj);
+       
+      }
+    }
+  }
 }
 //_____________________________________________________________________
 //

diff --git a/FMD/analysis/AliFMDAnalysisTaskDndeta.cxx b/FMD/analysis/AliFMDAnalysisTaskDndeta.cxx
index 934da9b7aa785324d930d45671424809e615ee6b..668d705db9e25bcb5c511aada32f1fcd7de01775 100644 (file)
--- a/FMD/analysis/AliFMDAnalysisTaskDndeta.cxx
+++ b/FMD/analysis/AliFMDAnalysisTaskDndeta.cxx
@@ -105,7 +105,7 @@ void AliFMDAnalysisTaskDndeta::CreateOutputObjects()
                              hBg->GetXaxis()->GetXmax(),
                              nSec, 0, 2*TMath::Pi());
            
-           hHits  = new TH1F(Form("hHits_FMD%d%c_vtxbin%d",det,ringChar,i),Form("hHits_FMD%d%c_vtxbin%d",det,ringChar,i),
+           hHits  = new TH1F(Form("hMCHits_FMD%d%c_vtxbin%d",det,ringChar,i),Form("hMCHits_FMD%d%c_vtxbin%d",det,ringChar,i),
                              hBg->GetNbinsX(),
                              hBg->GetXaxis()->GetXmin(),
                              hBg->GetXaxis()->GetXmax());
@@ -274,7 +274,7 @@ void AliFMDAnalysisTaskDndeta::ProcessPrimary() {
          //Double_t x,y,z;
          
          Float_t   eta   = pars->GetEtaFromStrip(det,ring,sec,strip,vertex.At(2));//-1*TMath::Log(TMath::Tan(0.5*theta));
-         TH1F* hHits = (TH1F*)fOutputList->FindObject(Form("hHits_FMD%d%c_vtxbin%d",det,ring,vertexBin));
+         TH1F* hHits = (TH1F*)fOutputList->FindObject(Form("hMCHits_FMD%d%c_vtxbin%d",det,ring,vertexBin));
          hHits->Fill(eta);
          Float_t nstrips = (ring =='O' ? 256 : 512);
          

diff --git a/FMD/analysis/AliFMDAnalysisTaskSE.cxx b/FMD/analysis/AliFMDAnalysisTaskSE.cxx
index 7f5d414914c8d4736340ed488f1c0eb4b87019f0..2fc8e55d0d8ad54f2dfa2edb57ae512a59538490 100644 (file)
--- a/FMD/analysis/AliFMDAnalysisTaskSE.cxx
+++ b/FMD/analysis/AliFMDAnalysisTaskSE.cxx
@@ -101,6 +101,7 @@ void AliFMDAnalysisTaskSE::UserExec(Option_t */*option*/)
 //_____________________________________________________________________
 void AliFMDAnalysisTaskSE::Terminate(Option_t */*option*/)
 {
+  fBackground.Terminate("");
   fDndeta.Terminate("");
 
 }