Updates for jet core task
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 8 Nov 2012 12:03:32 +0000 (12:03 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 8 Nov 2012 12:03:32 +0000 (12:03 +0000)
JETAN/AliAnalysisTaskJetCluster.cxx
JETAN/DEV/AliAnalysisTaskJetCluster.cxx
PWGJE/AliAnalysisTaskJetCore.cxx

index e83e0f7..c994580 100644 (file)
@@ -1507,11 +1507,19 @@ void AliAnalysisTaskJetCluster::UserExec(Option_t */*option*/)
      Double_t sigma2=0.;
      Double_t meanarea2=0.;
 
+     Double_t bkg4=0;
+     Double_t sigma4=0.;
+     Double_t meanarea4=0.;
+
      clustSeq.get_median_rho_and_sigma(jets2, range, true, bkg1, sigma1, meanarea1, true);
      fAODJetBackgroundOut->SetBackground(0,bkg1,sigma1,meanarea1);
 
      //     fh1BiARandomCones[0]->Fill(omCone-(bkg1*areaRandomCone));    
-     //  fh1BiARandomConesRan[0]->Fill(ptRandomConeRan-(bkg1*areaRandomCone));    
+     //  fh1BiARandomConesRan[0]->Fill(ptRandomConeRan-(bkg1*areaRandomCone));      
+       clustSeq.get_median_rho_and_sigma(sortedJets, range, true, bkg4, sigma4, meanarea4, true);
+     fAODJetBackgroundOut->SetBackground(3,bkg4,sigma4,meanarea4);
+
+     //////////////////////
      
      clustSeq.get_median_rho_and_sigma(jets2, range, false, bkg2, sigma2, meanarea2, true);
      fAODJetBackgroundOut->SetBackground(1,bkg2,sigma2,meanarea2);
index e83e0f7..c994580 100644 (file)
@@ -1507,11 +1507,19 @@ void AliAnalysisTaskJetCluster::UserExec(Option_t */*option*/)
      Double_t sigma2=0.;
      Double_t meanarea2=0.;
 
+     Double_t bkg4=0;
+     Double_t sigma4=0.;
+     Double_t meanarea4=0.;
+
      clustSeq.get_median_rho_and_sigma(jets2, range, true, bkg1, sigma1, meanarea1, true);
      fAODJetBackgroundOut->SetBackground(0,bkg1,sigma1,meanarea1);
 
      //     fh1BiARandomCones[0]->Fill(omCone-(bkg1*areaRandomCone));    
-     //  fh1BiARandomConesRan[0]->Fill(ptRandomConeRan-(bkg1*areaRandomCone));    
+     //  fh1BiARandomConesRan[0]->Fill(ptRandomConeRan-(bkg1*areaRandomCone));      
+       clustSeq.get_median_rho_and_sigma(sortedJets, range, true, bkg4, sigma4, meanarea4, true);
+     fAODJetBackgroundOut->SetBackground(3,bkg4,sigma4,meanarea4);
+
+     //////////////////////
      
      clustSeq.get_median_rho_and_sigma(jets2, range, false, bkg2, sigma2, meanarea2, true);
      fAODJetBackgroundOut->SetBackground(1,bkg2,sigma2,meanarea2);
index 22f15f2..7fa5c5d 100644 (file)
@@ -537,12 +537,12 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
    AliInputEventHandler* inputHandler = (AliInputEventHandler*)
    ((AliAnalysisManager::GetAnalysisManager())->GetInputEventHandler());
         std::cout<<inputHandler->IsEventSelected()<<" "<<fOfflineTrgMask<<std::endl;
-   if(!(inputHandler->IsEventSelected() & fOfflineTrgMask)){
-      if(fDebug) Printf(" Trigger Selection: event REJECTED ... ");
-      fHistEvtSelection->Fill(2);
-      PostData(1, fOutputList);
-      return;
-   }
+          if(!(inputHandler->IsEventSelected() & fOfflineTrgMask)){
+        if(fDebug) Printf(" Trigger Selection: event REJECTED ... ");
+        fHistEvtSelection->Fill(2);
+          PostData(1, fOutputList);
+          return;
+           }
 
    // vertex selection
    if(!aod){
@@ -620,6 +620,8 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
      if(externalBackground)rho = externalBackground->GetBackground(0);}
    if(fFlagRandom==1){
       if(externalBackground)rho = externalBackground->GetBackground(2);}
+   if(fFlagRandom==2){
+      if(externalBackground)rho = externalBackground->GetBackground(3);}
 
    // fetch jets
    TClonesArray *aodJets[2];