#_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Randomizer #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskParticleRandomizer.C #Module.MacroArgs "tracks", "tracks_randomized", kTRUE, kTRUE, kFALSE #Module.Deps #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); // Load flow TFile* fileInput = TFile::Open("alien:///alice/cern.ch/user/r/rhaake/FlowV2V3V4.root"); TH2D* flowV2 = static_cast(fileInput->Get("v2_EP_PbPb")); TH2D* flowV3 = static_cast(fileInput->Get("v3_EP_PbPb")); TH2D* flowV4 = static_cast(fileInput->Get("v4_EP_PbPb")); __R_ADDTASK__->SetDistributionV2(flowV2); __R_ADDTASK__->SetDistributionV3(flowV3); __R_ADDTASK__->SetDistributionV4(flowV4); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_R030_JetFinder #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "","", 1, 0.3, 1, kRHTrackPtCut, 0.0, kGhostArea, 1, "Jet", 0., kFALSE, kFALSE #Module.Deps ChargedJetsHadronCF_Randomizer #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetJetEtaRange(-0.4,0.4); AliTrackContainer* trackCont = new AliTrackContainer(kTracksRH); trackCont->SetFilterHybridTracks(kTRUE); trackCont->SetParticlePtCut(kRHTrackPtCut); __R_ADDTASK__->AdoptParticleContainer(trackCont); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Rho_JetFinderKT_R020 #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "","", 0, 0.2, 1, kRHTrackPtCut, 0.0, kGhostArea, 1, "Jet", 0., kFALSE, kFALSE #Module.Deps ChargedJetsHadronCF_Randomizer #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); AliTrackContainer* trackCont = new AliTrackContainer(kTracksRH); trackCont->SetFilterHybridTracks(kTRUE); trackCont->SetParticlePtCut(kRHTrackPtCut); __R_ADDTASK__->AdoptParticleContainer(trackCont); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_R030_TaskCorrelations #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGCF/Correlations/macros/dphicorrelations/AddTaskPhiCorrelations.C #Module.MacroArgs 0, 0, 0, kTRUE, 0, "histosPhiCorrelations_R030", "ChargedJetsHadronCF" #Module.Deps ChargedJetsHadronCF_R030_Support #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetEventSelectionBit(kPhysSelRH); __R_ADDTASK__->SetCustomParticlesA("ChargedJetsHadronCF_Particles_R030"); __R_ADDTASK__->SetTriggersFromDetector(6); __R_ADDTASK__->SetCustomParticlesB("ChargedJetsHadronCF_TrackParticles_R030"); __R_ADDTASK__->SetAssociatedFromDetector(7); __R_ADDTASK__->SetFilterBit(272); // hybrid tracks (same as in the jet task) __R_ADDTASK__->SetCustomBinning(kRHCustomBinning.Data()); __R_ADDTASK__->SetMixingTracks(50000); // Event pool saving if(kRHExportEventPool) { Double_t minCent = 0; Double_t maxCent = 100; Double_t minZvtx = -10; Double_t maxZvtx = 10; Double_t minPt = 4; Double_t maxPt = 20; __R_ADDTASK__->SetMixingTracks(20000); __R_ADDTASK__->SetUsePtBinnedEventPool(kTRUE); __R_ADDTASK__->AddEventPoolsToOutput(minCent, maxCent, minZvtx, maxZvtx, minPt, maxPt); } // Load external pool manager if(kRHImportEventPool) { __R_ADDTASK__->SetMixingTracks(20000); TFile* fileInput = TFile::Open("alien:///alice/cern.ch/user/r/rhaake/ExternalEventPool.root"); AliEventPoolManager* poolMgr = static_cast(fileInput->Get("AliEventPoolManager")); if(!poolMgr) { std::cout << "Pool manager not successfully loaded.\n"; return -1; } __R_ADDTASK__->SetUsePtBinnedEventPool(kTRUE); __R_ADDTASK__->SetExternalEventPoolManager(poolMgr); } #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Rho #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskRhoNew.C #Module.MacroArgs kTracksRH, "", "RhoR020KT", 0.2, AliJetContainer::kTPCfid, AliJetContainer::kChargedJet, kTRUE, AliJetContainer::pt_scheme, "Rho_ExLJ" #Module.Deps ChargedJetsHadronCF_Rho_JetFinderKT_R020 #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetVzRange(-10,10); __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); __R_ADDTASK__->SetExcludeLeadJets(2); // default AliJetContainer* cont = __R_ADDTASK__->GetJetContainer(0); cont->SetParticlePtCut(kRHTrackPtCut); cont->SetArrayName(kJetsKTR020RH.Data()); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_R030_Support #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskChargedJetsHadronCF.C #Module.MacroArgs kPhysSelRH, kTracksRH, kJetsAKTR030RH, kRHRhoName, 0.3, 0.6, 0.15 #Module.Deps ChargedJetsHadronCF_R030_JetFinder,ChargedJetsHadronCF_Rho #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetVzRange(-10,10); __R_ADDTASK__->SetNumberOfCentralityBins(100); __R_ADDTASK__->SelectCollisionCandidates(kPhysSelRH); __R_ADDTASK__->SetJetParticleArrayName("ChargedJetsHadronCF_Particles_R030"); __R_ADDTASK__->SetTrackParticleArrayName("ChargedJetsHadronCF_TrackParticles_R030"); __R_ADDTASK__->GetJetContainer()->SetMaxTrackPt(kRHMaxLeadingHadronPt); if(kRHUseDijetTrigger) { __R_ADDTASK__->SetJetOutputMode(kRHDijetTriggerMode); __R_ADDTASK__->SetEventCriteriumSelection(3); __R_ADDTASK__->SetEventCriteriumLeadingJets(kRHDijetTriggerLeadingPt,kRHDijetTriggerSubleadingPt, 5./6 * TMath::Pi()); } if(kRHUseFakefactorRejection) { TFile* inputFile = TFile::Open("alien:///alice/cern.ch/user/r/rhaake/fakefactorsR030.root"); TProfile* cutprofile = static_cast(inputFile->Get("hFakeFactorMean")); __R_ADDTASK__->ActivateFakejetRejection(0., 1., *cutprofile); } #Module.EndConfig