function to create the standard jets added

diff --git a/PWG4/macros/AddTaskJets.C b/PWG4/macros/AddTaskJets.C
index 82c6a71396ca329ad4153d23b1ca9a21e10fc5eb..67163838d5120416b2f8b1ad23d865a4d941bfed 100644 (file)
--- a/PWG4/macros/AddTaskJets.C
+++ b/PWG4/macros/AddTaskJets.C
@@ -3,6 +3,24 @@ AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius = -1);
 \r
 AliAnalysisTaskJets *AddTaskJets(Char_t *jr, Char_t *jf,Float_t radius = -1); // for the new AF\r
 \r
 \r
 AliAnalysisTaskJets *AddTaskJets(Char_t *jr, Char_t *jf,Float_t radius = -1); // for the new AF\r
 \r

+AliAnalysisTaskJets *AddTaskJets(){\r
+  // fills the standard "jets" branch in the AOD\r
+  // need the ESDFilter to run before, to access the AODtracks\r
+  // Tracks selected by the first Filter (1<<0)\r
+  // needs to be adapted for different cuts\r
+  \r
+  // UA1 as standard chosen, since it is the most robust and simple JF\r
+  // R = 0.4 suffficient to provide accurate jet axis for correlation studies\r
+  // energy resolution suffers a little\r
+  // Acceptance of jets not limited by the Jet Finder but should be done\r
+  // by user to abs(eta) < 0.5 \r
+\r
+  return AddTaskJets("AOD","UA1",0.4);\r
+\r
+}\r
+\r
+\r
+\r

 AliAnalysisTaskJets *AddTaskJets(Char_t *jr, Char_t *jf, Float_t radius)\r
 {\r
   // Creates a jet finder task, configures it and adds it to the analysis manager.\r
 AliAnalysisTaskJets *AddTaskJets(Char_t *jr, Char_t *jf, Float_t radius)\r
 {\r
   // Creates a jet finder task, configures it and adds it to the analysis manager.\r


@@ -100,13 +118,6 @@ AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius){
     jh->SetRparam(0.4); // setup parameters                                  \r
     if(radius>0)jh->SetRparam(radius);\r
     jh->SetAlgorithm(2); // antikt from fastjet/JetDefinition.hh\r
     jh->SetRparam(0.4); // setup parameters                                  \r
     if(radius>0)jh->SetRparam(radius);\r
     jh->SetAlgorithm(2); // antikt from fastjet/JetDefinition.hh\r

-    // tmp CKB for checking with leticia\r
-    jh->SetRparam(0.2); // setup parameters \r
-    jh->SetPtMin(0.); \r
-    jh->SetRapRange(-0.9,0.9); \r
-    jh->SetGhostArea(0.01); \r
-    jh->SetGhostEtaMax(7.);\r
-    // BKC\r

     jetFinder = new AliFastJetFinder();\r
     if (jh) jetFinder->SetJetHeader(jh);\r
     break;\r
     jetFinder = new AliFastJetFinder();\r
     if (jh) jetFinder->SetJetHeader(jh);\r
     break;\r


@@ -124,7 +135,7 @@ AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius){
     jh->SetLegoNbinEta(274);\r
     jh->SetLegoEtaMin(-2);\r
     jh->SetLegoEtaMax(+2);\r
     jh->SetLegoNbinEta(274);\r
     jh->SetLegoEtaMin(-2);\r
     jh->SetLegoEtaMax(+2);\r

-    jh->SetMinJetEt(10.);\r
+    jh->SetMinJetEt(5.);\r

     jh->SetJetEtaMax(1.5);\r
     jh->SetJetEtaMin(-1.5);\r
 \r
     jh->SetJetEtaMax(1.5);\r
     jh->SetJetEtaMin(-1.5);\r
 \r


@@ -144,7 +155,7 @@ AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius){
     jh->SetLegoNbinEta(274);\r
     jh->SetLegoEtaMin(-2);\r
     jh->SetLegoEtaMax(+2);\r
     jh->SetLegoNbinEta(274);\r
     jh->SetLegoEtaMin(-2);\r
     jh->SetLegoEtaMax(+2);\r

-    jh->SetMinJetEt(10.);\r
+    jh->SetMinJetEt(5.);\r

     jh->SetJetEtaMax(1.5);\r
     jh->SetJetEtaMin(-1.5);\r
     jetFinder = new AliUA1JetFinderV1();\r
     jh->SetJetEtaMax(1.5);\r
     jh->SetJetEtaMin(-1.5);\r
     jetFinder = new AliUA1JetFinderV1();\r