#_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_PbPb_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_Wrapper_PbPb #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kPhysSelRH); //__R_ADDTASK__->SetJetEtaRange(-0.4,0.4); AliTrackContainer* trackCont = new AliTrackContainer("tracks_PbPb"); trackCont->SetFilterHybridTracks(kTRUE); trackCont->SetParticlePtCut(kRHTrackPtCut); __R_ADDTASK__->AdoptParticleContainer(trackCont); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_PbPb_Rho #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskRhoNew.C #Module.MacroArgs "tracks_PbPb", "", "RhoR020KT", 0.2, AliJetContainer::kTPCfid, AliJetContainer::kChargedJet, kTRUE, AliJetContainer::pt_scheme, "Rho_ExLJ" #Module.Deps ChargedJetsHadronCF_PbPb_Rho_JetFinderKT #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetVzRange(-10,10); __R_ADDTASK__->SelectCollisionCandidates(kPhysSelRH); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); __R_ADDTASK__->SetExcludeLeadJets(2); // default AliJetContainer* cont = __R_ADDTASK__->GetJetContainer(0); cont->SetParticlePtCut(kRHTrackPtCut); cont->SetArrayName(Form("Jet_KTChargedR020_tracks_PbPb_pT%04d_pt_scheme",(Int_t)(kRHTrackPtCut*1000))); /* AliTrackContainer* trackCont = __R_ADDTASK__->GetTrackContainer(0); trackCont->SetFilterHybridTracks(kTRUE); trackCont->SetParticlePtCut(kRHTrackPtCut); __R_ADDTASK__->AdoptParticleContainer(trackCont); */ #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_PbPb_Rho_JetFinderKT #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_Wrapper_PbPb #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kPhysSelRH); AliTrackContainer* trackCont = new AliTrackContainer("tracks_PbPb"); trackCont->SetFilterHybridTracks(kTRUE); trackCont->SetParticlePtCut(kRHTrackPtCut); __R_ADDTASK__->AdoptParticleContainer(trackCont); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_PbPb_Support #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskChargedJetsHadronCF.C #Module.MacroArgs kPhysSelRH, "tracks_PbPb", Form("Jet_AKTChargedR030_tracks_PbPb_pT%04d_pt_scheme",(Int_t)(kRHTrackPtCut*1000)), "RhoR020KT", 0.3, 0.6, kRHMinLeadingHadronPt, 0.15, kRHTrackPtCut, 0.707, "" #Module.Deps ChargedJetsHadronCF_PbPb_JetFinder,ChargedJetsHadronCF_PbPb_Rho #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetVzRange(-10,10); __R_ADDTASK__->SetNumberOfCentralityBins(10); __R_ADDTASK__->SelectCollisionCandidates(kPhysSelRH); __R_ADDTASK__->SetJetParticleArrayName("ChargedJetsHadronCF_Jets"); __R_ADDTASK__->SetTrackParticleArrayName("ChargedJetsHadronCF_Tracks"); __R_ADDTASK__->GetJetContainer()->SetMaxTrackPt(kRHMaxLeadingHadronPt); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_PbPb_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_PbPb_Support #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetEventSelectionBit(kPhysSelRH); __R_ADDTASK__->SetCustomParticlesA("ChargedJetsHadronCF_Jets"); __R_ADDTASK__->SetTriggersFromDetector(6); __R_ADDTASK__->SetCustomParticlesB("ChargedJetsHadronCF_Tracks"); __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); __R_ADDTASK__->SetCheckEventNumberInMixedEvent(kTRUE); __R_ADDTASK__->SetUseVtxAxis(2); // Efficiency map corrections if(kRHUseEfficiencyMaps) { cout << "### Using efficiency maps ###\n"; // PbPb, AOD149, 2011, hybrid tracks file = TFile::Open("alien:///alice/cern.ch/user/j/jgrosseo/pbpb/correction_2011_aod149_hybrid.root"); //file = TFile::Open("alien:///alice/cern.ch/user/j/jgrosseo/pp/correction_276_global.root"); __R_ADDTASK__->SetEfficiencyCorrectionAssociated((THnF*) file->Get("correction")); } // Event pool saving if(kRHExportEventPool) { Double_t minCent = 0; Double_t maxCent = 50; Double_t minZvtx = -10; Double_t maxZvtx = 10; Double_t minPt = 0.5; Double_t maxPt = 40; __R_ADDTASK__->SetMixingTracks(15000); __R_ADDTASK__->SetUsePtBinnedEventPool(kTRUE); __R_ADDTASK__->AddEventPoolsToOutput(minCent, maxCent, minZvtx, maxZvtx, minPt, maxPt); } else { __R_ADDTASK__->SetMixingTracks(15000); TFile* fileInput = 0; if(kRHMinLeadingHadronPt > 4.0) { std::cout << "Using 5 GeV/c biased event pool.\n"; fileInput = TFile::Open("alien:///alice/cern.ch/user/r/rhaake/ExternalPoolPbPb_Large_5GeVBias.root"); } else { if(kRHTrackPtCut > 2.99) { std::cout << "alien:///alice/cern.ch/user/r/rhaake/ExternalPoolPbPb3GeV_Large.root" << std::endl; fileInput = TFile::Open("alien:///alice/cern.ch/user/r/rhaake/ExternalPoolPbPb3GeV_Large.root"); std::cout << fileInput << std::endl; } else if(kRHTrackPtCut > 1.99) { std::cout << "alien:///alice/cern.ch/user/r/rhaake/ExternalPoolPbPb2GeV_Large.root" << std::endl; fileInput = TFile::Open("alien:///alice/cern.ch/user/r/rhaake/ExternalPoolPbPb2GeV_Large.root"); std::cout << fileInput << std::endl; } else if(kRHTrackPtCut > 0.99) { std::cout << "alien:///alice/cern.ch/user/r/rhaake/ExternalPoolPbPb1GeV_Large.root" << std::endl; fileInput = TFile::Open("alien:///alice/cern.ch/user/r/rhaake/ExternalPoolPbPb1GeV_Large.root"); std::cout << fileInput << std::endl; } else fileInput = TFile::Open("alien:///alice/cern.ch/user/r/rhaake/ExternalPoolPbPb_Large.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_Wrapper_PbPb #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskParticleRandomizer.C #Module.MacroArgs "tracks", "tracks_PbPb", kFALSE, kFALSE, kFALSE #Module.Deps #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kPhysSelRH); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); __R_ADDTASK__->SetTrackEfficiency(kRHTrackEfficiency); #Module.EndConfig