In AliMUONESDInterface:

[u/mrichter/AliRoot.git] / JETAN / AliFastJetFinder.cxx

diff --git a/JETAN/AliFastJetFinder.cxx b/JETAN/AliFastJetFinder.cxx
index 896c457e8e2026b9c85acf83a0ed502d44a9cae0..f324a21d91e1d793c3ce8192dedb1edfc799217c 100644 (file)
--- a/JETAN/AliFastJetFinder.cxx
+++ b/JETAN/AliFastJetFinder.cxx
@@ -384,7 +384,7 @@ void AliFastJetFinder::RunTest(const char* datafile)
 
 //____________________________________________________________________________
 
-void AliFastJetFinder::WriteJHeaderToFile()
+void AliFastJetFinder::WriteJHeaderToFile() const
 {
   fHeader->Write();
 }
@@ -425,8 +425,23 @@ Bool_t AliFastJetFinder::ProcessEvent()
   // Jets
   FindJets(); 
 
-  fJetBkg->SetHeader(fHeader);
-  fJetBkg->SetReader(fReader);
+
+  if( fAODEvBkg){
+    fJetBkg->SetHeader(fHeader);
+    fJetBkg->SetReader(fReader);
+    Double_t sigma1 = 0,meanarea1= 0,sigma2 = 0,meanarea2 = 0;
+    Double_t bkg1 = 0,bkg2 = 0;
+    
+    fJetBkg->SetFastJetInput(fInputFJ);
+    fJetBkg->BkgFastJetb(bkg1,sigma1,meanarea1);
+    fJetBkg->BkgFastJetWoHardest(bkg2,sigma2,meanarea2);
+    fAODEvBkg->SetBackground(0,bkg1,sigma1,meanarea1);
+    fAODEvBkg->SetBackground(1,bkg2,sigma2,meanarea2);
+  }
+
+
+
+
   /*
   fJetBkg->SetFastJetInput(fInputFJ);
   Double_t bkg1=fJetBkg->BkgFastJet();
@@ -466,9 +481,19 @@ Bool_t AliFastJetFinder::ProcessEvent2()
   // Jets
   FindJets();
   
-  fJetBkg->SetHeader(fHeader);
-  fJetBkg->SetReader(fReader);
-  fJetBkg->SetFastJetInput(fInputFJ);
+  if( fAODEvBkg){
+    fJetBkg->SetHeader(fHeader);
+    fJetBkg->SetReader(fReader);
+    fJetBkg->SetFastJetInput(fInputFJ);
+    Double_t sigma1,meanarea1,sigma2,meanarea2;
+    Double_t bkg1,bkg2;
+    fJetBkg->BkgFastJetb(bkg1,sigma1,meanarea1);
+    fJetBkg->BkgFastJetWoHardest(bkg2,sigma2,meanarea2);
+    fAODEvBkg->SetBackground(0,bkg1,sigma1,meanarea1);
+    fAODEvBkg->SetBackground(1,bkg2,sigma2,meanarea2);
+  }
+
+
 //  Double_t bkg1=fJetBkg->BkgFastJet();
 //  Double_t bkg2=fJetBkg->BkgChargedFastJet();
 //  Double_t bkg3=fJetBkg->BkgFastJetCone(fAODjets);