Background subraction disabled since not correct. (Elena Bruna)

diff --git a/JETAN/AliFastJetFinder.cxx b/JETAN/AliFastJetFinder.cxx
index 6db9adac5b5897b6cd26d60acc7de2ac845d9e92..a98f162fec53add4b45c0100c3c369d6593beeb1 100644 (file)
--- a/JETAN/AliFastJetFinder.cxx
+++ b/JETAN/AliFastJetFinder.cxx
@@ -469,15 +469,15 @@ Bool_t AliFastJetFinder::ProcessEvent2()
   fJetBkg->SetHeader(fHeader);
   fJetBkg->SetReader(fReader);
   fJetBkg->SetFastJetInput(fInputFJ);
-  Double_t bkg1=fJetBkg->BkgFastJet();
-  Double_t bkg2=fJetBkg->BkgChargedFastJet();
-  Double_t bkg3=fJetBkg->BkgFastJetCone(fAODjets);
-  Double_t bkg4=fJetBkg->BkgRemoveJetLeading(fAODjets);
+//  Double_t bkg1=fJetBkg->BkgFastJet();
+//  Double_t bkg2=fJetBkg->BkgChargedFastJet();
+//  Double_t bkg3=fJetBkg->BkgFastJetCone(fAODjets);
+//  Double_t bkg4=fJetBkg->BkgRemoveJetLeading(fAODjets);
   
-  fAODEvBkg->SetBackground(0,bkg1);
-  fAODEvBkg->SetBackground(1,bkg2);
-  fAODEvBkg->SetBackground(2,bkg3);
-  fAODEvBkg->SetBackground(3,bkg4);
+//  fAODEvBkg->SetBackground(0,bkg1);
+//  fAODEvBkg->SetBackground(1,bkg2);
+//  fAODEvBkg->SetBackground(2,bkg3);
+//  fAODEvBkg->SetBackground(3,bkg4);
   
   Int_t nEntRef    = ref->GetEntries();
 

diff --git a/JETAN/AliJetBkg.cxx b/JETAN/AliJetBkg.cxx
index f160ade077f5ce980c5a34e1a81b9ef6e5bbb8d1..27d5b894bc0d9d12dacd0133932f7703ba8a92bf 100644 (file)
--- a/JETAN/AliJetBkg.cxx
+++ b/JETAN/AliJetBkg.cxx
@@ -149,160 +149,6 @@ Float_t AliJetBkg::BkgStat()
   return enTot/accEMCal;
 
 }
-/////////////////////////////////
-Float_t AliJetBkg::BkgRemoveJetLeading(TClonesArray* fAODJets)
-{
-  // Remove the particles of the
-  // two largest jets using the track references  stored in the AODJet from the estimation of new rho. 
-
-  AliFastJetHeaderV1 *header = (AliFastJetHeaderV1*)fHeader;
-  Bool_t debug  = header->GetDebug();     // debug option
-
-  if(debug)cout<<"==============AliJetBkg::BkgRemoveJetLeading()============="<<endl;
-
- 
-  // check if we are reading AOD jets
-  TRefArray *refs = 0;
-  Bool_t fromAod = !strcmp(fReader->ClassName(),"AliJetESDReader");
-  if (fromAod) { refs = fReader->GetReferences(); }
-  
-  //Hard wired Calorimeter area (get it later from the AliJetReaderHeader.h)
-  Double_t accEMCal=2*0.7*110./180*TMath::Pi();//2.68 area of EMCal
-
-  Int_t nJ=fAODJets->GetEntries(); //this must be the # of jets... 
-  if(debug)cout<<"nJets:  "<<nJ<<endl;
-  
-
-  //begin unit array 
-  TClonesArray* fUnit = fReader->GetUnitArray();
-  if(fUnit == 0) { cout << "Could not get the momentum array" << endl; return -99; }
-
-  Int_t nIn = fUnit->GetEntries();
-  if(nIn == 0) { cout << "entries = 0 ; Event empty !!!" << endl ; return -99; }
-  
-  Float_t rhoback=0.0;
-  Float_t jetarea1=0.0,jetarea2=0.0;
-    
-  Int_t particlejet1=-99;
-  Int_t particlejet2=-99;
-  TRefArray *refarray1 = 0;
-  TRefArray *refarray2 = 0;
-  Int_t nJettracks1 = 0, nJettracks2 = 0;
-  Int_t acc=0,acc1=0;
-  AliAODJet *jet1;
-  AliAODJet *jet2;
-
-  if(nJ==1){
-    jet1 = dynamic_cast<AliAODJet*>(fAODJets->At(0));
-    jetarea1=jet1->EffectiveAreaCharged();
-    Float_t jetPhi=jet1->Phi();
-    Float_t jetEta=jet1->Eta();
-    if(jetPhi>1.396 && jetPhi<3.316 && jetEta>-0.7 && jetEta<0.7)acc=1;
-    refarray1=jet1->GetRefTracks();
-    nJettracks1=refarray1->GetEntries();
-    if(debug)cout<<"nJ = 1, acc="<<acc<<"  jetarea1="<<jetarea1<<endl;
-
-  }
-
-  if(nJ>=2){
-    jet1 = dynamic_cast<AliAODJet*>(fAODJets->At(0));
-    jetarea1=jet1->EffectiveAreaCharged();
-    Float_t jetPhi1=jet1->Phi();
-    Float_t jetEta1=jet1->Eta();
-    if(jetPhi1>1.396 && jetPhi1<3.316 && jetEta1>-0.7 && jetEta1<0.7)acc=1;
-    refarray1=jet1->GetRefTracks();
-    nJettracks1=refarray1->GetEntries();
-    if(debug)cout<<"npart = "<<nJettracks1<<endl;
-   
-    jet2 = dynamic_cast<AliAODJet*>(fAODJets->At(1));
-    jetarea2=jet2->EffectiveAreaCharged();
-    Float_t jetPhi2=jet2->Phi();
-    Float_t jetEta2=jet2->Eta();
-    if(jetPhi2>1.396 && jetPhi2<3.316 && jetEta2>-0.7 && jetEta2<0.7)acc1=1;
-    refarray2=jet2->GetRefTracks();
-    nJettracks2=refarray2->GetEntries();
-    if(debug)cout<<"nJ = "<<nJ<<", acc="<<acc<<"  acc1="<<acc1<<"  jetarea1="<<jetarea1<<"  jetarea2="<<jetarea2<<endl;
-  }
-
-
-  
-  // cout<<" nIn = "<<nIn<<endl;
-  Float_t sumPt=0;
-  Float_t eta,phi,pt;
-  Int_t ipart=0;
-
-  for(Int_t i=0; i<nIn; i++) 
-    { //Unit Array Loop
-      AliJetUnitArray *uArray = (AliJetUnitArray*)fUnit->At(i);
-
-      if(uArray->GetUnitEnergy()>0.){
-       eta   = uArray->GetUnitEta();
-       phi   = uArray->GetUnitPhi();
-       pt = uArray->GetUnitEnergy();
-       //      cout<<"ipart = "<<ipart<<" eta="<<eta<<"  phi="<<phi<<endl;
-       if(phi>1.396 && phi<3.316 && eta>-0.7 && eta<0.7){
-         //cout<<sumPt<<endl;
-           sumPt+=pt;
-
-           
-           if(nJ==1 && acc==1){
-             for(Int_t ii=0; ii<nJettracks1;ii++){    
-               
-               particlejet1  = ((AliJetUnitArray*)refarray1->At(ii))->GetUnitTrackID();
-        
-               if(ipart==particlejet1) {
-                 sumPt-=pt;
-               }
-             }
-           }
-         
-         
-           if(nJ>=2){
-             
-             //first jet
-             if(acc==1){
-               for(Int_t ii=0; ii<nJettracks1;ii++){    
-                 particlejet1  = ((AliJetUnitArray*)refarray1->At(ii))->GetUnitTrackID();
-                 
-               //cout<<"uArr loop = "<<i<<"  ipart in uArr (1/2)="<<ipart<<"  part in jet="<<ii<<"  partID="<<particlejet1<<" sumPt="<<sumPt<<endl; 
-                 if(ipart==particlejet1) {
-                   sumPt-=pt;
-                 }
-               }
-             }
-             if(acc1==1){
-               //second jet
-               for(Int_t ii=0; ii<nJettracks2;ii++){ 
-                 particlejet2  = ((AliJetUnitArray*)refarray2->At(ii))->GetUnitTrackID();
-                 //cout<<"uArr loop = "<<i<<"  ipart in uArr (2/2)="<<ipart<<"  part in jet="<<ii<<"  partID="<<particlejet2<<" sumPt="<<sumPt<<endl; 
-                 if(ipart==particlejet2) {
-                   sumPt-=pt;
-                 }
-               }
-             }
-             
-             
-           }
-       
-       }//if phi,eta
-       ipart++;
-      }//end if energy
-    }// end unit array loop          
-  
-
-  Float_t areasum = accEMCal-acc*jetarea1-acc1*jetarea2;
-  if(debug)cout<<"pt sum   "<<sumPt<<" area  "<<areasum<<endl;
-   
-  if(nJ>0) rhoback=sumPt/areasum;
-  else rhoback=0.;
-  if(debug)cout<<" rho from leading jet paricle array removed   "<<rhoback<<endl;
- 
-  return rhoback;
- 
-}
-
-
-
 
 ////////////////////////////////////////////////////////////////////////
 Float_t AliJetBkg::BkgFastJetCone(TClonesArray* fAODJets)