new settings, larger binning for QA, p_T threshold for spectrum, number of processed...
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 28 May 2010 05:40:00 +0000 (05:40 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 28 May 2010 05:40:00 +0000 (05:40 +0000)
PWG4/macros/AddTaskJetSpectrum2.C
PWG4/macros/AddTaskPWG4HighPtQATPConly.C
PWG4/macros/AnalysisTrainPWG4Jets.C

index a4bc4d4..af9e364 100644 (file)
@@ -70,6 +70,8 @@ AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2(const char* bRec,const char* bG
 \r
    pwg4spec->SetFilterMask(filterMask); \r
    pwg4spec->SetUseGlobalSelection(kTRUE); \r
+   pwg4spec->SetMinJetPt(5.);\r
+\r
 \r
    if(type == "AOD"){\r
      // Assume all jet are produced already\r
index 7f5006f..eb24796 100644 (file)
@@ -74,14 +74,14 @@ AliPWG4HighPtQATPConly* AddTaskPWG4HighPtQATPConly(int cuts=1)//1: Standard Cuts
   taskPWG4QA->SetCuts(trackCuts);
   taskPWG4QA->SetCutsITS(trackCutsITS);
   taskPWG4QA->SetMaxCosmicAngle(0.008);
-  
+  taskPWG4QA->SetCutType(cuts);
+
  
   // E. Create ONLY the output containers for the data produced by the task.
   // Get and connect other common input/output containers via the manager as below
   //==============================================================================
 
   //------ input data ------
-  //  AliAnalysisDataContainer *cinput0  = mgr->GetCommonInputContainer();
   TString outputfile = "";
   outputfile = AliAnalysisManager::GetCommonFileName();
   outputfile += Form(":PWG4_HighPtQATPConly%d",cuts);
@@ -90,27 +90,29 @@ AliPWG4HighPtQATPConly* AddTaskPWG4HighPtQATPConly(int cuts=1)//1: Standard Cuts
   AliAnalysisDataContainer *cout_hist1;
   AliAnalysisDataContainer *cout_hist2;
   AliAnalysisDataContainer *cout_hist3;
+  AliAnalysisDataContainer *cout_cuts0;
+  AliAnalysisDataContainer *cout_cuts1;
 
-  if(cuts==1) {
-    cout_hist0 = mgr->CreateContainer("qa_histsCuts1", TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);
-    cout_hist1 = mgr->CreateContainer("qa_histsTPCCuts1", TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);
-    cout_hist2 = mgr->CreateContainer("qa_histsITSCuts1", TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);  
-    cout_hist3 = mgr->CreateContainer("qa_histsCosmicsCuts1", TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);  
-  }
-  else if(cuts==2) {
-    cout_hist0 = mgr->CreateContainer("qa_histsCuts2", TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);
-    cout_hist1 = mgr->CreateContainer("qa_histsTPCCuts2", TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);
-    cout_hist2 = mgr->CreateContainer("qa_histsITSCuts2", TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);  
-    cout_hist3 = mgr->CreateContainer("qa_histsCosmicsCuts2", TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);  
-  }
+  cout_hist0 = mgr->CreateContainer(Form("qa_histsCuts%d",cuts), TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);
+  cout_hist1 = mgr->CreateContainer(Form("qa_histsTPCCuts%d",cuts), TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);
+  cout_hist2 = mgr->CreateContainer(Form("qa_histsITSCuts%d",cuts), TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);  
+  cout_hist3 = mgr->CreateContainer(Form("qa_histsCosmicsCuts%d",cuts), TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);  
+  cout_cuts0 = mgr->CreateContainer(Form("qa_trackCuts%d",cuts), AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
+  cout_cuts1 = mgr->CreateContainer(Form("qa_trackCutsITS%d",cuts), AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
+
+  //Add task to manager
   mgr->AddTask(taskPWG4QA);
 
+  //Connect input containter to manager
   mgr->ConnectInput(taskPWG4QA,0,mgr->GetCommonInputContainer());
+
+  //Connect output containers to manager
   mgr->ConnectOutput(taskPWG4QA,0,cout_hist0);
   mgr->ConnectOutput(taskPWG4QA,1,cout_hist1);
   mgr->ConnectOutput(taskPWG4QA,2,cout_hist2);
   mgr->ConnectOutput(taskPWG4QA,3,cout_hist3);
+  mgr->ConnectOutput(taskPWG4QA,4,cout_cuts0);
+  mgr->ConnectOutput(taskPWG4QA,5,cout_cuts1);
 
   // Return task pointer at the end
   return taskPWG4QA;
index 4f2308b..0856713 100644 (file)
@@ -52,6 +52,7 @@ Bool_t      kSkipTerminate      = kTRUE; // Do not call Teminate
 Bool_t      kUseDate            = kFALSE; // use date in train name
 Bool_t      kUseDebug           = kTRUE; // activate debugging
 Int_t       kUseSysInfo         = 0; // activate debugging
+Long_t      kNumberOfEvents     = 1234567890; // number of events to process from the chain
 Bool_t      kUseMC              = kTRUE;  // use MC info
 Bool_t      kIsMC               = kTRUE;  // is MC info, if false it overwrites Use(AOD)MC
 Bool_t      kUseAODMC           = kTRUE;  // use MC infA
@@ -116,8 +117,8 @@ Int_t       kProofOffset = 0;
 //== grid plugin setup variables
 Bool_t      kPluginUse         = kTRUE;   // do not change
 Bool_t      kPluginUseProductionMode  = kFALSE;   // use the plugin in production mode
-TString     kPluginRootVersion       = "v5-26-00b-2";  // *CHANGE ME IF MORE RECENT IN GRID*
-TString     kPluginAliRootVersion    = "v4-19-04-AN";  // *CHANGE ME IF MORE RECENT IN GRID*                                          
+TString     kPluginRootVersion       = "v5-26-00b-6";  // *CHANGE ME IF MORE RECENT IN GRID*
+TString     kPluginAliRootVersion    = "v4-19-13-AN";  // *CHANGE ME IF MORE RECENT IN GRID*                                          
 Bool_t      kPluginMergeViaJDL       = kTRUE;  // merge via JDL
 Bool_t      kPluginFastReadOption   = kFALSE;  // use xrootd tweaks
 Bool_t      kPluginOverwriteMode    = kTRUE;  // overwrite existing collections
@@ -350,8 +351,16 @@ void AnalysisTrainPWG4Jets(const char *analysis_mode="local",
    if (iDIJETAN) {
       gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/AddTaskDiJets.C");
       AliAnalysisTaskDiJets *taskdijets = 0;
-      taskdijets = AddTaskDiJets(); 
+      if(iDIJETAN&1)taskdijets = AddTaskDiJets(); 
       if (!taskdijets) ::Warning("AnalysisTrainPWG4Jets", "AliAnalysisTaskJets cannot run for this train conditions - EXCLUDED");
+      if(iDIJETAN&2){
+       taskdijets = AddTaskDiJets("jetsAOD_CDF07"); 
+       taskdijets = AddTaskDiJets("jetsAOD_DA07"); 
+       taskdijets = AddTaskDiJets("jetsAOD_FASTJET07"); 
+       taskdijets = AddTaskDiJets("jetsAOD_FASTKT07"); 
+       taskdijets = AddTaskDiJets("jetsAOD_SISCONE07"); 
+       taskdijets = AddTaskDiJets("jetsAOD_UA107");
+      }
    }
 
    if(iPWG1QASym){
@@ -539,7 +548,9 @@ void AnalysisTrainPWG4Jets(const char *analysis_mode="local",
       //      TString gcArguments = "-run-on-train -run-jet -run-chic -run-neutralmeson -run-cf";
       //      TString gcArguments = "-run-on-train -run-jet -run-neutralmeson -run-cf -use-own-xyz";
       //      TString gcArguments = "-run-on-train -run-jet -run-neutralmeson -run-cf -use-own-xyz";
-      TString gcArguments = "-run-on-train -run-jet -run-neutralmeson -run-cf -use-own-xyz -bg-off";
+      TString gcArguments = "-run-on-train -run-jet -run-omega-meson -run-neutralmeson";
+      TString kGCAnalysisCutSelectionId="9002111000";
+      gcArguments.Append(Form("-set-cut-selection  %s ",kGCAnalysisCutSelectionId.Data()));
       if(!kIsMC)gcArguments += " -mc-off";
       AliAnalysisTaskGammaConversion * taskGammaConversion = AddTaskGammaConversion(gcArguments,mgr->GetCommonInputContainer());
       gSystem->ChangeDirectory(cdir);
@@ -637,7 +648,7 @@ void StartAnalysis(const char *mode, TChain *chain) {
             ::Error("AnalysisTrainPWG4Jets.C::StartAnalysis", "Cannot create the chain");
             return;
          }   
-         mgr->StartAnalysis(mode, chain);
+         mgr->StartAnalysis(mode, chain,kNumberOfEvents);
          return;
       case 1:
          if (!kProofDataSet.Length()) {