PWGJE/EMCALJetTasks/macros/AddTaskFullpAJets.C

   1 // AddTaskFullpAJets.C 2013-02-03 16:58:39Z yaldo
   2
   3 void AddTaskFullpAJets()
   4 {
   5     const char *usedTracks="PicoTracks";
   6     const char *usedClusters="CaloClusters";
   7     const char *outClusName="CaloClustersCorr";
   8     const Double_t hadcorr=2.0;
   9     const Double_t minTrackPt=0.15;
  10     const Double_t minClusterPt=0.30;
  11     const Double_t minChargedJetPt=0.15;
  12     const Double_t minFullJetPt=0.15;
  13     const Double_t Eexcl=0.00;
  14     const Double_t phiMatch=0.03;
  15     const Double_t etaMatch=0.015;
  16
  17     // Some constants for the jet finders
  18     const Int_t cKT                 = 0;
  19     const Int_t cANTIKT             = 1;
  20     const Int_t cFULLJETS           = 0;
  21     const Int_t cCHARGEDJETS        = 1;
  22     const Int_t cNEUTRALJETS        = 2;
  23
  24     AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  25     if (!mgr)
  26     {
  27         Error("AddTaskJetCommon","No analysis manager found.");
  28         return 0;
  29     }
  30
  31     // Jet finders (RECONSTRUCTED DATA)
  32     TString tmpTaskName("");
  33     AliEmcalJetTask* jetFinderTask = NULL;
  34
  35     gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C");
  36     // ########## CHARGED JETS ##########
  37     // R=0.2
  38     jetFinderTask = AddTaskEmcalJet(usedTracks,"",cANTIKT,0.2,1,minTrackPt,minClusterPt);
  39     RequestMemory(jetFinderTask,250*1024);//more memory
  40
  41     // R=0.4
  42     jetFinderTask = AddTaskEmcalJet(usedTracks,"",cANTIKT,0.4,1,minTrackPt,minClusterPt);
  43     RequestMemory(jetFinderTask,250*1024);//more memory
  44
  45     // ########## FULL JETS ##########
  46     // last two settings are for min pt tracks/clusters
  47     // R=0.2, anti-kT
  48     jetFinderTask = AddTaskEmcalJet(usedTracks,outClusName,cANTIKT,0.2,cFULLJETS,minTrackPt,minClusterPt);
  49     RequestMemory(jetFinderTask,250*1024);//more memory
  50
  51     // R=0.2 kT
  52     jetFinderTask = AddTaskEmcalJet(usedTracks,outClusName,cKT,0.2,cFULLJETS,minTrackPt,minClusterPt);
  53     RequestMemory(jetFinderTask,250*1024);//more memory
  54
  55     // R=0.4
  56     jetFinderTask = AddTaskEmcalJet(usedTracks,outClusName,cANTIKT,0.4,cFULLJETS,minTrackPt,minClusterPt);
  57     RequestMemory(jetFinderTask,250*1024);//more memory
  58
  59     // R=0.4 kT
  60     jetFinderTask = AddTaskEmcalJet(usedTracks,outClusName,cKT,0.4,cFULLJETS,minTrackPt,minClusterPt);
  61     RequestMemory(jetFinderTask,250*1024);//more memory
  62
  63     // Add User Tasks'
  64     // Run with R=0.2
  65     AliAnalysisTaskFullpAJets *task1 = new AliAnalysisTaskFullpAJets("FileR2");
  66     mgr->AddTask(task1);
  67     task1->SetR_JET(2);
  68     AliAnalysisDataContainer *coutput1 =
  69     mgr->CreateContainer("R2List",
  70                          TList::Class(),
  71                          AliAnalysisManager::kOutputContainer,
  72                          "FullpAJetsR2.root");
  73     mgr->ConnectInput(task1,0,mgr->GetCommonInputContainer());
  74     mgr->ConnectOutput(task1,1,coutput1);
  75     RequestMemory(task1,250*1024);//more memory
  76
  77     // Run with R=0.4
  78     AliAnalysisTaskFullpAJets *task2 = new AliAnalysisTaskFullpAJets("FileR4");
  79     mgr->AddTask(task2);
  80     task2->SetR_JET(4);
  81     AliAnalysisDataContainer *coutput2 =
  82     mgr->CreateContainer("R4List",
  83                          TList::Class(),
  84                          AliAnalysisManager::kOutputContainer,
  85                          "FullpAJetsR4.root");
  86     mgr->ConnectInput(task2,0,mgr->GetCommonInputContainer());
  87     mgr->ConnectOutput(task2,1,coutput2);
  88     RequestMemory(task2,250*1024);//more memory
  89 }