extra k_T jet finders configs and AODMC2b reading (charged particles in accepatance...
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 31 Jan 2010 10:02:39 +0000 (10:02 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 31 Jan 2010 10:02:39 +0000 (10:02 +0000)
PWG4/macros/AddTaskJets.C

index 1bfbb71..bddccba 100644 (file)
@@ -58,21 +58,21 @@ Int_t AddTaskJetsDelta(char *nonStdFile,UInt_t filterMask,Bool_t kUseAODMC,UInt_
   Int_t iCount = 0;\r
 \r
   // Jet Fidners Selected by run flag first bit 2^0 second by 2^1 etc\r
-  const char *cJF[13]        = {"UA1","UA1","UA1","CDF","DA","SISCONE","FASTJET","FASTKT","UA1LO","FASTKT","FASTKT","FASTKT","FASTKT"};\r
-  const Float_t radius[13]   = {  0.4,  0.7,  1.0,  0.7, 0.7,      0.4,      0.4,     0.4,    0.7,    0.1 ,    0.2 ,    0.6 ,    0.8 };\r
-  UInt_t  flag[13]           = {    6,    7,    7,    7,   7,        7,        7,       7,      7,      7 ,      7 ,      7 ,      7 };\r
+  const char *cJF[14]        = {"UA1","UA1","UA1","CDF","DA","SISCONE","FASTJET","FASTKT","UA1LO","FASTKT","FASTKT","FASTKT","FASTKT","FASTKT"};\r
+  const Float_t radius[14]   = {  0.4,  0.7,  1.0,  0.7, 0.7,      0.4,      0.4,     0.4,    0.7,    0.1 ,    0.2 ,    0.6 ,    0.8 ,   0.001};\r
+  UInt_t  flag[14]           = {    6,    7,    7,    7,   7,        7,       7,       15,      7,     9 ,     9,     9 ,     9 ,      9};\r
   // flag[5] = 0; // set siscone to 0 for proof mode...\r
-  // flag first bit AOD, second bit AODMC2 third bit AODMC2\r
+  // flag first bit AOD, second bit AODMC2 third bit AODMC2 third (8) bit AOODMC2b (limited acceptance)\r
   // i.e. 7 all, 6 only MC2 and MC\r
   // this stay at three\r
-  const char *cReader[3] = {"AOD","AODMC","AODMC2"};  \r
+  const char *cReader[4] = {"AOD","AODMC","AODMC2","AODMC2b"};  \r
 \r
   \r
 \r
-  for(int i = 0; i< 13;i++){\r
+  for(int i = 0; i< 14;i++){\r
     if(!(runFlag&(1<<i)))continue;\r
     if(!kUseAODMC)flag[i]&=1; // switch OFF MC if we do not have it\r
-    for(int ib = 0;ib<3;ib++){      \r
+    for(int ib = 0;ib<4;ib++){      \r
       if(flag[i]&(1<<ib)){\r
        jetana = AddTaskJets(cReader[ib],cJF[i],radius[i],filterMask);\r
        if(jetana){\r
@@ -381,6 +381,16 @@ AliJetReader *CreateJetReader(Char_t *jr,UInt_t filterMask){
     er = new AliJetAODReader();\r
     er->SetReaderHeader(jrh);\r
     break;\r
+  case "AODMC2b":\r
+    AliJetAODReaderHeader *jrh = new AliJetAODReaderHeader();\r
+    jrh->SetComment("AOD MC Reader");\r
+    jrh->SetPtCut(0.);\r
+    jrh->SetFiducialEta(-0.9,0.9); // to take all MC particles default is 0.9\r
+    jrh->SetReadAODMC(2);// 1 all primary MC , 2 all primary charged\r
+    // Define reader and set its header\r
+    er = new AliJetAODReader();\r
+    er->SetReaderHeader(jrh);\r
+    break;\r
 \r
   default:\r
     ::Error("AddTaskJets", "Wrong jet reader selected\n");\r