//____________________________________________________________________
AliAnalysisTaskUE:: AliAnalysisTaskUE(const char* name):
AliAnalysisTask(name, ""),
-fDebug(kFALSE),
+fDebug(0),
fDeltaAOD(kFALSE),
fDeltaAODBranch(""),
- fArrayJets(0x0),
+fArrayJets(0x0),
fAOD(0x0),
fAODjets(0x0),
-//fArrayJets(0x0),
-//fDeltaAOD(kFALSE),
-//fDeltaAODBranch("jetsAOD"),
fListOfHistos(0x0),
fBinsPtInHist(30),
fMinJetPtInHist(0.),
Double_t xsec = fh1Xsec->GetBinContent(1);
Double_t ntrials = fh1Trials->GetBinContent(1);
Double_t normFactor = xsec/ntrials;
- Printf("xSec %f nTrials %f Norm %f \n",xsec,ntrials,normFactor);
+ if(fDebug > 1)Printf("xSec %f nTrials %f Norm %f \n",xsec,ntrials,normFactor);
TCanvas* c2 = new TCanvas("c2","Jet Pt dist",160,160,1200,800);
-AliAnalysisTaskJFSystematics *AddTaskJFSystematics(char *jf1 = "jetsMC",char *jf2 = "jets")\r
+AliAnalysisTaskJFSystematics *AddTaskJFSystematics(char *jf1 = "jets",char *jf2 = "jetsAODMC_UA104")\r
{\r
\r
// Get the pointer to the existing analysis manager via the static access method.\r
pwg4jfs->SetAODInput(kTRUE);\r
}\r
\r
- pwg4jfs->SetBranchGen(jf1); \r
- pwg4jfs->SetBranchRec(jf2); \r
+ pwg4jfs->SetBranchGen(jf2); \r
+ pwg4jfs->SetBranchRec(jf1); \r
mgr->AddTask(pwg4jfs);\r
\r
// Create ONLY the output containers for the data produced by the task.\r
\r
void ConfigTaskUE(AliAnalysisTaskUE * ueana); // common config, extend with different cases\r
\r
-AliAnalysisTaskUE *AddTaskUE()\r
+AliAnalysisTaskUE *AddTaskUE(char *jetBranch = "")\r
{\r
// Creates a jet fider task, configures it and adds it to the analysis manager.\r
\r
\r
AliAnalysisTaskUE* ueana = new AliAnalysisTaskUE("Underlying Event");\r
ConfigTaskUE(ueana);\r
+\r
+ TString jb(jetBranch);\r
+\r
+ if(jb.Length()>0){\r
+ ueana->ReadDeltaAOD(kTRUE);\r
+ ueana->SelectDeltaAODBranch(jb.Data());\r
+ }\r
+\r
mgr->AddTask(ueana);\r
\r
// Create ONLY the output containers for the data produced by the task.\r
// Get and connect other common input/output containers via the manager as below\r
//==============================================================================\r
- AliAnalysisDataContainer *coutput1_UE = mgr->CreateContainer("histosUE", TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:PWG4_UE",AliAnalysisManager::GetCommonFileName()));\r
+ AliAnalysisDataContainer *coutput1_UE = 0;\r
+ if(jb.Length()==0)coutput1_UE = mgr->CreateContainer("histosUE", TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:PWG4_UE",AliAnalysisManager::GetCommonFileName()));\r
+ else coutput1_UE = mgr->CreateContainer(Form("histosUE_%s",jb.Data()), TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:PWG4_UE_%s",AliAnalysisManager::GetCommonFileName(),jb.Data()));\r
\r
mgr->ConnectInput (ueana, 0, mgr->GetCommonInputContainer());\r
mgr->ConnectOutput (ueana, 0, coutput1_UE );\r
Double_t rad=0.7; \r
Double_t deltaPhiCut = 2.616;\r
\r
- ueana->SetDebugLevel(10); \r
+ ueana->SetDebugLevel(0); \r
ueana->SetPtRangeInHist(25, 0., 250.);\r
ueana->SetAnaTopology(anaType); \r
ueana->SetRegionType(regType); \r