-AliAnalysisTaskJetCluster *AddTaskJetCluster(char* bRec = "AOD",char* bGen = "",UInt_t filterMask = 16, Int_t iPhysicsSelection = 1,Char_t *jf = "KT", Float_t radius = 0.4,Int_t nSkip = 0);\r
+AliAnalysisTaskJetCluster *AddTaskJetCluster(char* bRec = "AOD",char* bGen = "",UInt_t filterMask = 16, Int_t iPhysicsSelection = 1,Char_t *jf = "KT", Float_t radius = 0.4,Int_t nSkip = 0,Int_t kWriteAOD);\r
\r
\r
AliAnalysisTaskJetCluster *AddTaskJetClusterDelta(UInt_t filterMask = 16,Bool_t kUseAODMC = kFALSE,Int_t iPhysicsSelection = 1,Char_t *jf = "KT", UInt_t iFlag){\r
if(iFlag&(1<<0))js = AddTaskJetCluster("AOD","",filterMask,iPhysicsSelection,jf,0.00001);\r
if(iFlag&(1<<1))js = AddTaskJetCluster("AOD","",filterMask,iPhysicsSelection,jf,0.1);\r
if(iFlag&(1<<2))js = AddTaskJetCluster("AOD","",filterMask,iPhysicsSelection,jf,0.2);\r
- if(iFlag&(1<<4))js = AddTaskJetCluster("AOD","",filterMask,iPhysicsSelection,jf,0.4);\r
+ if(iFlag&(1<<4))js = AddTaskJetCluster("AOD","",filterMask,iPhysicsSelection,jf,0.4,0,1);\r
if(iFlag&(1<<6))js = AddTaskJetCluster("AOD","",filterMask,iPhysicsSelection,jf,0.6);\r
if(iFlag&(1<<8))js = AddTaskJetCluster("AOD","",filterMask,iPhysicsSelection,jf,0.8);\r
if(iFlag&(1<<10))js = AddTaskJetCluster("AOD","",filterMask,iPhysicsSelection,jf,1.0);\r
}\r
\r
\r
-AliAnalysisTaskJetCluster *AddTaskJetCluster(char* bRec,char* bGen ,UInt_t filterMask,Int_t iPhysicsSelection,Char_t *jf,Float_t radius,Int_t nSkip)\r
+AliAnalysisTaskJetCluster *AddTaskJetCluster(char* bRec,char* bGen ,UInt_t filterMask,Int_t iPhysicsSelection,Char_t *jf,Float_t radius,Int_t nSkip,Int_t kWriteAOD)\r
{\r
// Creates a jet fider task, configures it and adds it to the analysis manager.\r
\r
}\r
\r
\r
- if(TMath::Abs(radius-0.4)<0.01){\r
+ if(kWriteAOD){\r
pwg4spec->SetJetOutputBranch(Form("clusters%s_%s%s",bRec,jf,cRadius));\r
pwg4spec->SetJetOutputMinPt(0); // store only jets / clusters above a certain threshold\r
}\r
-AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2(const char* bRec = "jets",const char* bGen = "jetsAODMC_UA104",const char* bBkg="jeteventbackground_jetsAOD_FASTKT04",UInt_t filterMask = 16, Int_t iPhysicsSelection = 1,UInt_t iEventSelectionMask = 0,Bool_t kBackground=kFALSE);\r
+AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2(const char* bRec = "jets",const char* bGen = "jetsAODMC_UA104",const char* bBkg="jeteventbackground_clustersAOD_KT06",UInt_t filterMask = 16, Int_t iPhysicsSelection = 1,UInt_t iEventSelectionMask = 0,Bool_t kBackground=kTRUE);\r
\r
\r
AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2Delta(UInt_t filterMask = 16,Bool_t kUseAODMC = kFALSE,Int_t iPhysicsSelection = 1,UInt_t iFlag = 0xfffffff, UInt_t iEventSelectionMask = 0,Bool_t kBackground = kFALSE){\r
AliAnalysisTaskJetSpectrum2 *js = 0;\r
\r
- if(!kBackground){\r
- if(kUseAODMC){\r
- if(iFlag&(1<<0))js = AddTaskJetSpectrum2("jets","jetsAODMC_UA104","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- if(iFlag&(1<<1))js = AddTaskJetSpectrum2("jets","jetsAODMC2_UA104","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
-\r
- if(iFlag&(1<<2))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","jetsAODMC_FASTJET04","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- if(iFlag&(1<<3))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","jetsAODMC2_FASTJET04","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
-\r
- if(iFlag&(1<<4)){\r
- js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","jetsAODMC_FASTKT04","",filterMask,iPhysicsSelection,iEventSelectionMask);\r
- }\r
- if(iFlag&(1<<5))js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","jetsAODMC2_FASTKT04","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- if(iFlag&(1<<6))js = AddTaskJetSpectrum2("jetsAOD_UA107","jetsAODMC_UA107","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- }\r
-\r
- if(iFlag&(1<<7))js = AddTaskJetSpectrum2("jets","jetsAOD_FASTJET04","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
-\r
- if(iFlag&(1<<8))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- if(iFlag&(1<<9))js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- if(iFlag&(1<<10))js = AddTaskJetSpectrum2("jetsAOD_SISCONE04","","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- \r
- if(iFlag&(1<<11)){\r
- js = AddTaskJetSpectrum2("jetsAOD_UA107","","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- js->SetRecEtaWindow(0.2);\r
- }\r
- if(iFlag&(1<<12)){\r
- js = AddTaskJetSpectrum2("jetsAOD_FASTJET07","","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- js->SetRecEtaWindow(0.2);\r
- }\r
- if(iFlag&(1<<13)){\r
- js = AddTaskJetSpectrum2("jetsAOD_FASTKT07","","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- js->SetRecEtaWindow(0.2);\r
- }\r
- if(iFlag&(1<<14)){\r
- js = AddTaskJetSpectrum2("jetsAOD_SISCONE07","","",filterMask,iPhysicsSelection, iEventSelectionMask);\r
- js->SetRecEtaWindow(0.2);\r
- }}\r
+ TString cBack = "";\r
+ if(kBackground)cBack = "jeteventbackground_clustersAOD_KT06";\r
\r
-\r
- if(kBackground){\r
if(kUseAODMC){\r
- if(iFlag&(1<<0))js = AddTaskJetSpectrum2("jets","jetsAODMC_UA104","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
- if(iFlag&(1<<1))js = AddTaskJetSpectrum2("jets","jetsAODMC2_UA104","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
-\r
- if(iFlag&(1<<2))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","jetsAODMC_FASTJET04","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
- if(iFlag&(1<<3))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","jetsAODMC2_FASTJET04","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
-\r
+ if(iFlag&(1<<0))js = AddTaskJetSpectrum2("jets","jetsAODMC_UA104",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ if(iFlag&(1<<1))js = AddTaskJetSpectrum2("jets","jetsAODMC2_UA104",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ \r
+ if(iFlag&(1<<2))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","jetsAODMC_FASTJET04",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ if(iFlag&(1<<3))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","jetsAODMC2_FASTJET04",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ \r
if(iFlag&(1<<4)){\r
- js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","jetsAODMC_FASTKT04","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection,iEventSelectionMask,kBackground);\r
+ js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","jetsAODMC_FASTKT04",cBack.Data(),filterMask,iPhysicsSelection,iEventSelectionMask,kBackground);\r
}\r
- if(iFlag&(1<<5))js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","jetsAODMC2_FASTKT04","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
- if(iFlag&(1<<6))js = AddTaskJetSpectrum2("jetsAOD_UA107","jetsAODMC_UA107","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ if(iFlag&(1<<5))js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","jetsAODMC2_FASTKT04",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ if(iFlag&(1<<6))js = AddTaskJetSpectrum2("jetsAOD_UA107","jetsAODMC_UA107",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
}\r
-\r
- if(iFlag&(1<<7))js = AddTaskJetSpectrum2("jets","jetsAOD_FASTJET04","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
-\r
- if(iFlag&(1<<8))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
- if(iFlag&(1<<9))js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
- if(iFlag&(1<<10))js = AddTaskJetSpectrum2("jetsAOD_SISCONE04","","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ \r
+ if(iFlag&(1<<7))js = AddTaskJetSpectrum2("jets","jetsAOD_FASTJET04",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ \r
+ if(iFlag&(1<<8))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ if(iFlag&(1<<9))js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ if(iFlag&(1<<10))js = AddTaskJetSpectrum2("jetsAOD_SISCONE04","",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
\r
if(iFlag&(1<<11)){\r
- js = AddTaskJetSpectrum2("jetsAOD_UA107","","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ js = AddTaskJetSpectrum2("jetsAOD_UA107","",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
js->SetRecEtaWindow(0.2);\r
}\r
if(iFlag&(1<<12)){\r
- js = AddTaskJetSpectrum2("jetsAOD_FASTJET07","","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ js = AddTaskJetSpectrum2("jetsAOD_FASTJET07","",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
js->SetRecEtaWindow(0.2);\r
}\r
if(iFlag&(1<<13)){\r
- js = AddTaskJetSpectrum2("jetsAOD_FASTKT07","","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ js = AddTaskJetSpectrum2("jetsAOD_FASTKT07","",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
js->SetRecEtaWindow(0.2);\r
}\r
if(iFlag&(1<<14)){\r
- js = AddTaskJetSpectrum2("jetsAOD_SISCONE07","","jeteventbackground_jetsAOD_FASTKT04",filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
+ js = AddTaskJetSpectrum2("jetsAOD_SISCONE07","",cBack.Data(),filterMask,iPhysicsSelection, iEventSelectionMask,kBackground);\r
js->SetRecEtaWindow(0.2);\r
- }}\r
-\r
+ }\r
\r
return js;\r
}\r
\r
pwg4spec->SetBranchGen(bGen); \r
pwg4spec->SetBranchRec(bRec); \r
- pwg4spec->SetBranchBkg(bBkg);\r
+ if(strlen(bBkg)>0)pwg4spec->SetBranchBkg(bBkg);\r
\r
pwg4spec->SetFilterMask(filterMask); \r
pwg4spec->SetUseGlobalSelection(kTRUE); \r
pwg4spec->SetMinJetPt(5.);\r
- if(kBackground) pwg4spec->SetBkgSubtraction(kTRUE);\r
- \r
+ if(kBackground&&strlen(bBkg)>0){\r
+ pwg4spec->SetBkgSubtraction(kTRUE);\r
+ } \r
+\r
if(type == "AOD"){\r
// Assume all jet are produced already\r
pwg4spec->SetAODJetInput(kTRUE);\r
TString kLocalDataList = "local_deltaaod.txt"; // Change local xml dataset for local interactive analysis
// == local process variables
+TString kPluginMode;
// Temporaries.
// Main analysis train macro. If a configuration file is provided, all parameters
// are taken from there but may be altered by CheckModuleFlags.
- if (strlen(config_file) && !LoadConfig(config_file)) return;
+ // these flag may be needed by the config file
+ kPluginMode = plugin_mode;
+
+ if (strlen(config_file) && !LoadConfig(config_file)) return;
if(iOffset)kProofOffset = iOffset;
TString smode(analysis_mode);
smode.ToUpper();
if (kSaveTrain)WriteConfig();
- // Check compatibility of selected modules
+ // Check compatibility of selected modules
CheckModuleFlags(smode);
// gROOT->ProcessLine(".trace");
gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/AddTaskJetCluster.C");
AliAnalysisTaskJetCluster *taskCl = 0;
if(iPWG4Cluster&1){
- taskCl = AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelection,"KT");
+ taskCl = AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelection,"KT",0.4,0,1); // this one is for the background jets
+ taskCl = AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelection,"KT",0.5,0,1); // this one is for the background jets
+ taskCl = AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelection,"KT",0.6,0,1); // this one is for the background jets
+ taskCl = AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelection,"KT",0.7,0,1); // this one is for the background jets
+
if (!taskCl) ::Warning("AnalysisTrainPWG4Jets", "AliAnalysisTaskCluster cannot run for this train conditions - EXCLUDED");
}
if(iPWG4Cluster&2){
}
else selection = 1<<0|1<<1|1<<2|1<<3|1<<4|1<<5|1<<7|1<<8|1<<9;
AddTaskJetClusterDelta(kHighPtFilterMask,kUseAODMC,iPhysicsSelection,"ANTIKT",selection);
- AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelection,"ANTIKT",0.4,1);
+ // AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelection,"ANTIKT",0.4,0,1);
}
}