Added fastjet

diff --git a/PWG4/macros/AddTaskJets.C b/PWG4/macros/AddTaskJets.C
index 6dfde860d3c9745d1edac99a8acd9bfee408d2e6..a4eb1ceae4c36be1475e765dc30940903bea4c82 100644 (file)
--- a/PWG4/macros/AddTaskJets.C
+++ b/PWG4/macros/AddTaskJets.C
@@ -61,6 +61,7 @@ AliAnalysisTaskJets *AddTaskJets(Char_t *jr, Char_t *jf, Float_t radius)
 }\r
 \r
 AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius){\r
+  AliJetFinder *jetFinder = 0;\r
 \r
   switch (jf) {\r
   case "CDF":\r
@@ -83,12 +84,12 @@ AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius){
     break;\r
 \r
 \r
-  case "Fastjet":\r
-    AliFastJetHeader *jh = new AliFastJetHeader();\r
-    jh->SetRparam(0.7); // setup parameters                                  \r
-\r
+  case "FASTJET":\r
+    AliFastJetHeaderV1 *jh = new AliFastJetHeaderV1();\r
+    jh->SetRparam(0.4); // setup parameters                                  \r
+    if(radius>0)jh->SetRparam(radius);\r
+    jh->SetAlgorithm(2); // antikt from fastjet/JetDefinition.hh\r
     jetFinder = new AliFastJetFinder();\r
-    jetFinder->SetOutputFile("jets.root");\r
     if (jh) jetFinder->SetJetHeader(jh);\r
     break;\r
 \r
@@ -125,13 +126,12 @@ AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius){
     jh->SetMinJetEt(10.);\r
     jh->SetJetEtaMax(1.5);\r
     jh->SetJetEtaMin(-1.5);\r
-\r
     jetFinder = new AliUA1JetFinderV1();\r
     if (jh) jetFinder->SetJetHeader(jh);\r
     break;\r
-  case default:\r
-    ::Error("AddTaskJets", "Wrong jet finder selected\n");\r
-    return 0;\r
+  default:\r
+    Printf("\n >>>>>>> AddTaskJets Error Wrong jet finder selected\n");\r
+    break;\r
   }\r
 \r
   return jetFinder;\r