function to create the standard jets added
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 29 Oct 2009 09:10:38 +0000 (09:10 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 29 Oct 2009 09:10:38 +0000 (09:10 +0000)
PWG4/macros/AddTaskJets.C

index 82c6a71..6716383 100644 (file)
@@ -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
+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
@@ -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
-    // 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
@@ -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->SetMinJetEt(10.);\r
+    jh->SetMinJetEt(5.);\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->SetMinJetEt(10.);\r
+    jh->SetMinJetEt(5.);\r
     jh->SetJetEtaMax(1.5);\r
     jh->SetJetEtaMin(-1.5);\r
     jetFinder = new AliUA1JetFinderV1();\r