added jet finders for lower p_T (lo seed UA1, seed less cone, kT
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 16 Dec 2009 18:55:01 +0000 (18:55 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 16 Dec 2009 18:55:01 +0000 (18:55 +0000)
PWG4/macros/AddTaskJets.C

index f5befdd..7fca31d 100644 (file)
@@ -58,16 +58,16 @@ Int_t AddTaskJetsDelta(char *nonStdFile,UInt_t filterMask,Bool_t kUseAODMC){
   Int_t iCount = 0;\r
 \r
 \r
-  const char *cJF[7]        = {"UA1","UA1","UA1","CDF","DA","SISCONE","FASTJET"};\r
-  const Float_t radius[7]   = {  0.4,  0.7,  1.0,  0.7, 0.7,      0.4,      0.4};\r
-  UInt_t  flag[7]     = {    6,    7,    7,    7,   7,        7,        7};\r
-  flag[5] = 0; // set siscone to 0 for proof mode...\r
+  const char *cJF[9]        = {"UA1","UA1","UA1","CDF","DA","SISCONE","FASTJET","FASTKT","UA1LO"};\r
+  const Float_t radius[9]   = {  0.4,  0.7,  1.0,  0.7, 0.7,      0.4,      0.4,     0.4,    0.4};\r
+  UInt_t  flag[9]           = {    6,    7,    7,    7,   7,        7,        7,       7,       7};\r
+  // flag[5] = 0; // set siscone to 0 for proof mode...\r
   // flag first bit AOD, second bit AODMC2 third bit AODMC2\r
   // i.e. 7 all, 6 only MC2 and MC\r
   // this stay at three\r
   const char *cReader[3] = {"AOD","AODMC","AODMC2"};  \r
 \r
-  for(int i = 0; i< 7;i++){\r
+  for(int i = 0; i< 9;i++){\r
     if(!kUseAODMC)flag[i]&=1;\r
     for(int ib = 0;ib<3;ib++){      \r
       if(flag[i]&(1<<ib)){\r
@@ -196,6 +196,7 @@ AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius){
     if(radius>0)jh->SetRparam(radius);\r
     jh->SetAlgorithm(2); // antikt from fastjet/JetDefinition.hh\r
     jetFinder = new AliFastJetFinder();\r
+    jh->SetPtMin(1);\r
     if (jh) jetFinder->SetJetHeader(jh);\r
     break;\r
 \r
@@ -204,6 +205,7 @@ AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius){
     jh->SetRparam(0.4); // setup parameters                                  \r
     if(radius>0)jh->SetRparam(radius);\r
     jh->SetAlgorithm(0); // kt from fastjet/JetDefinition.hh\r
+    jh->SetPtMin(1);\r
     jetFinder = new AliFastJetFinder();\r
     if (jh) jetFinder->SetJetHeader(jh);\r
     break;\r
@@ -246,7 +248,6 @@ AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius){
     jh->SetRadius(0.4);\r
     if(radius>0)jh->SetRadius(radius);\r
     jh->SetEtSeed(4.);\r
-    jh->SetEtSeed(4.);\r
     jh->SetNAcceptJets(6);\r
     jh->SetLegoNbinPhi(432);\r
     jh->SetLegoNbinEta(274);\r
@@ -259,6 +260,27 @@ AliJetFinder *CreateJetFinder(Char_t *jf,Float_t radius){
     jetFinder = new AliUA1JetFinderV1();\r
     if (jh) jetFinder->SetJetHeader(jh);\r
     break;\r
+  case "UA1LO":\r
+    AliUA1JetHeaderV1 *jh=new AliUA1JetHeaderV1();\r
+    jh->SetComment("UA1 jet code with Lo Pt settings parameters");\r
+    jh->BackgMode(0);\r
+    jh->SetRadius(0.4);\r
+    if(radius>0)jh->SetRadius(radius);\r
+    jh->SetEtSeed(1.);\r
+    jh->SetNAcceptJets(6);\r
+    jh->SetLegoNbinPhi(432);\r
+    jh->SetLegoNbinEta(274);\r
+    jh->SetLegoEtaMin(-2);\r
+    jh->SetLegoEtaMax(+2);\r
+    jh->SetMinJetEt(1.);\r
+    jh->SetJetEtaMax(1.5);\r
+    jh->SetJetEtaMin(-1.5);\r
+\r
+    jetFinder = new AliUA1JetFinderV1();\r
+    if (jh) jetFinder->SetJetHeader(jh);\r
+    break;\r
+\r
+\r
 \r
   case "UA1MC":\r
     AliUA1JetHeaderV1 *jh=new AliUA1JetHeaderV1();\r