added unsubtracted jet container

diff --git a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskEmcalQGTagging.cxx b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskEmcalQGTagging.cxx
index dbf75909b8a0b96a88b38225cfede92b47c8567a..f3c405137a6b653b4c006bdb143e291dd99042ff 100644 (file)
--- a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskEmcalQGTagging.cxx
+++ b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskEmcalQGTagging.cxx
@@ -198,20 +198,39 @@ Bool_t AliAnalysisTaskEmcalQGTagging::FillHistograms()
       if (!jet1) continue;
       AliEmcalJet* jet2 = 0x0;
       fPtJet->Fill(jet1->Pt());
+      AliEmcalJet *jetUS = NULL;
+      Int_t ifound=0;
+      Int_t ilab=-1;
 
-      if (!(fJetShapeType == kData)) {
+        if (!(fJetShapeType == kData)) {
        AliPythiaInfo *partonsInfo = 0x0;
+
        if((fJetShapeType == kTrueDet) || (fJetShapeType == kDetEmb)){
          AliJetContainer *jetContTrue = GetJetContainer(1);
-         jet2 = jet1->ClosestJet();
+          AliJetContainer *jetContUS = GetJetContainer(2);
+          if(fJetShapeSub==kConstSub){
+         for(Int_t i = 0; i<jetContUS->GetNJets(); i++) {
+         jetUS = jetContUS->GetJet(i);
+         if(jetUS->GetLabel()==jet1->GetLabel()) {
+           ifound++;
+           if(ifound==1) ilab = i;
+         }
+       }   
+          if(ilab==-1) continue;
+         jetUS=jetContUS->GetJet(ilab);
+
+           jet2=jetUS->ClosestJet();}
+         if(!(fJetShapeSub==kConstSub)) jet2 = jet1->ClosestJet();
          if (!jet2) {
            Printf("jet2 not exists, returning");
            continue;
          }
-          fh2ResponseUW->Fill(jet1->Pt(),jet2->Pt());    
 
+          fh2ResponseUW->Fill(jet1->Pt(),jet2->Pt());    
 
-         Double_t fraction = jetCont->GetFractionSharedPt(jet1);
+          Double_t fraction=0;
+         if(!(fJetShapeSub==kConstSub))  fraction = jetCont->GetFractionSharedPt(jet1);
+          if(fJetShapeSub==kConstSub) fraction = jetContUS->GetFractionSharedPt(jetUS);
           cout<<"hey a jet"<<fraction<<" "<<jet1->Pt()<<" "<<jet2->Pt()<<" "<<fCent<<endl;
           
          if(fraction<fMinFractionShared) continue;

diff --git a/PWGJE/EMCALJetTasks/macros/AddTaskEmcalQGTagging.C b/PWGJE/EMCALJetTasks/macros/AddTaskEmcalQGTagging.C
index 7f0fd6f9d829b4eb85eeab2e9487fa2b579043ae..4d5b358661cee4fe42efabf1469bd51cb776fb4a 100644 (file)
--- a/PWGJE/EMCALJetTasks/macros/AddTaskEmcalQGTagging.C
+++ b/PWGJE/EMCALJetTasks/macros/AddTaskEmcalQGTagging.C
@@ -1,8 +1,10 @@
 AliAnalysisTaskEmcalQGTagging* AddTaskEmcalQGTagging(const char * njetsBase,
+                                                     const char * njetsUS,
                                                     const char * njetsTrue,
                                                     const Double_t R,
                                                     const char * nrhoBase, 
                                                     const char * ntracks, 
+                                                     const char * natracksUS,
                                                     const char * nclusters,
                                                     const char * ntracksTrue,
                                                     const char *type,                                
@@ -44,10 +46,12 @@ AliAnalysisTaskEmcalQGTagging* AddTaskEmcalQGTagging(const char * njetsBase,
   //task->SetVzRange(-10.,10.);
 
   AliParticleContainer *trackCont  = task->AddParticleContainer(ntracks);
+   AliParticleContainer *trackContUS  = task->AddParticleContainer(ntracksUS);
   AliParticleContainer *trackContTrue  = task->AddParticleContainer(ntracksTrue);
   AliClusterContainer *clusterCont = task->AddClusterContainer(nclusters);
 
   AliJetContainer *jetContBase=0x0;
+   AliJetContainer *jetContUS=0x0;
   AliJetContainer *jetContTrue=0x0;
 
   TString strType(type);
@@ -110,6 +114,16 @@ AliAnalysisTaskEmcalQGTagging* AddTaskEmcalQGTagging(const char * njetsBase,
       jetContTrue->SetPercAreaCut(0.6); 
       jetContTrue->SetPythiaInfoName("PythiaInfo");
     }
+    if(jetShapeSub==AliAnalysisTaskEmcalQGTagging::kConstSub){
+  if(jetContUS) {
+      jetContUS->SetRhoName(nrhoBase);
+      jetContUS->ConnectParticleContainer(trackContUS);
+      jetContUS->SetPercAreaCut(0.6); 
+      jetContUS->SetPythiaInfoName("PythiaInfo");
+  }}
+
+
+
     
   }