#_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_Embedder #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskParticleRandomizer.C #Module.MacroArgs "tracks_PbPb", "tracks_generated_embedded", kFALSE, kFALSE, kFALSE #Module.Deps ChargedJetsHadronCF_Embedding_Generator,ChargedJetsHadronCF_Embedding_JetFinderPbPb #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); __R_ADDTASK__->ActivateJetEmbedding("tracks_generated_pico"); if(kRHDoRandomize) { __R_ADDTASK__->SetRandomizeInPhi(kTRUE); __R_ADDTASK__->SetRandomizeInEta(kTRUE); __R_ADDTASK__->ActivateLeadingJetRemoval(kJetsAKTR030RH.Data(), 2, ""); 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")); __R_ADDTASK__->SetDistributionV2(flowV2); __R_ADDTASK__->SetDistributionV3(flowV3); } #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_Generator #Module.Libs pythia6_4_25 #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskJetEmbeddingFromPYTHIA.C #Module.MacroArgs "tracks_generated_pico", "", "", "", kPYTHIAPath, kNpTHardBins, kPtHardBinsScaling, "aodTree", "tracks", "", "", "", "lhc12a15e", kFALSE, -1, -1, 0, kMinPythiaJetPt, kFALSE, kTRUE #Module.Deps #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetRandomAccess(kTRUE); __R_ADDTASK__->SetTrackEfficiency(kRHTrackEfficiency); if(1) { __R_ADDTASK__->SetMarkMC(10000); } //######################## CONFIGURATION FOR DATA EMBEDDING // PWGJE/EMCALJetTasks/macros/AddTaskJetEmbeddingFromAOD.C // // "tracks_generated_pico", "", "", "", "alien:///alice/cern.ch/user/r/rhaake/filelist_10d.txt", "aodTree", "tracks", "", "", "", "lhc10d", kFALSE, -1, -1, 0, kFALSE, kTRUE //######################## CONFIGURATION FOR PYTHIA EMBEDDING // PWGJE/EMCALJetTasks/macros/AddTaskJetEmbeddingFromPYTHIA.C // //"tracks_generated_pico", "", "", "", kPYTHIAPath, kNpTHardBins, kPtHardBinsScaling, "aodTree", "tracks", "", "", "", "lhc12a15e", kFALSE, -1, -1, 0, kMinPythiaJetPt, kFALSE, kTRUE #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_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, "JetEmb", 0., kFALSE, kFALSE #Module.Deps ChargedJetsHadronCF_Embedding_Embedder #Module.Owner rhaake #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); AliTrackContainer* trackCont = new AliTrackContainer("tracks_generated_embedded"); trackCont->SetFilterHybridTracks(kTRUE); trackCont->SetParticlePtCut(kRHTrackPtCut); __R_ADDTASK__->AdoptParticleContainer(trackCont); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_JetFinderPbPb #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 EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); AliTrackContainer* trackCont = new AliTrackContainer("tracks_PbPb"); trackCont->SetFilterHybridTracks(kTRUE); trackCont->SetParticlePtCut(kRHTrackPtCut); __R_ADDTASK__->AdoptParticleContainer(trackCont); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_JetFinderPY #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, "JetPY", 0., kFALSE, kFALSE #Module.Deps ChargedJetsHadronCF_Embedding_Wrapper_PY #Module.Owner rhaake #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); AliTrackContainer* trackCont = new AliTrackContainer("tracks_generated_PY"); trackCont->SetFilterHybridTracks(kTRUE); trackCont->SetParticlePtCut(kRHTrackPtCut); __R_ADDTASK__->AdoptParticleContainer(trackCont); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_Rho #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskRhoNew.C #Module.MacroArgs "tracks_generated_embedded", "", kRHRhoName, 0.2, AliJetContainer::kTPCfid, AliJetContainer::kChargedJet, kTRUE, AliJetContainer::pt_scheme, "Rho_ExLJ" #Module.Deps ChargedJetsHadronCF_Embedding_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(kJetsKTR020RHEmbedded.Data()); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_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_Wrapper_PbPb #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); AliTrackContainer* trackCont = new AliTrackContainer("tracks_generated_embedded"); trackCont->SetFilterHybridTracks(kTRUE); trackCont->SetParticlePtCut(kRHTrackPtCut); __R_ADDTASK__->AdoptParticleContainer(trackCont); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_Support #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskChargedJetsHadronCF.C #Module.MacroArgs kPhysSelRH, "tracks_generated_embedded", kJetsAKTR030RHEmbedded.Data(), kRHRhoName, 0.3, 0.6, kRHMinLeadingHadronPt, 0.15, 0.15, 0.707, "" #Module.Deps ChargedJetsHadronCF_Embedding_JetFinderPY,ChargedJetsHadronCF_Embedding_JetFinder,ChargedJetsHadronCF_Embedding_JetFinderPbPb,ChargedJetsHadronCF_Embedding_Rho #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetVzRange(-10,10); __R_ADDTASK__->SelectCollisionCandidates(kPhysSelRH); __R_ADDTASK__->SetJetParticleArrayName("ChargedJetsHadronCF_Particles_Embedded_R030"); __R_ADDTASK__->SetTrackParticleArrayName("ChargedJetsHadronCF_TrackParticles_Embedded_R030"); __R_ADDTASK__->GetJetContainer()->SetMaxTrackPt(kRHMaxLeadingHadronPt); __R_ADDTASK__->SetNumberOfCentralityBins(10); __R_ADDTASK__->ActivateJetEmbedding(kJetsAKTR030RHPYTHIA.Data(), "tracks_generated_PY", 0.3, kRHNumberEmbeddedJets); __R_ADDTASK__->SetJetOutputMode(4); // 4: matching jets, 5: non-matching jets // Cuts for embedding Bool_t USE_LOW_CUT = kTRUE; Bool_t USE_NORMAL_CUT = kTRUE; Bool_t USE_HIGH_CUT = kTRUE; __R_ADDTASK__->SetCreateEmbeddingPtPlotPerCut(kFALSE); __R_ADDTASK__->SetJetVetoArrayName(kJetsAKTR030RHEmbeddedPbPb.Data()); __R_ADDTASK__->AddJetEmbeddingCut("NoCut", "ChargedJetsHadronCF_Particles_Embedded_R030_NoCut", 0., 999., 0, 100, 0.00, 0.0, 100.0, 999.); // ##################### MC percentage cut only, no jet vetos __R_ADDTASK__->AddJetEmbeddingCut(Form("lowerMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_lowerMCCut"), 0., 999., 0, 10, 0.15, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("lowerMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_lowerMCCut"), 0., 999., 10, 30, 0.20, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("lowerMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_lowerMCCut"), 0., 999., 30, 50, 0.20, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("lowerMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_lowerMCCut"), 0., 999., 50, 90, 0.20, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("normalMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_normalMCCut"), 0., 999., 0, 10, 0.20, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("normalMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_normalMCCut"), 0., 999., 10, 30, 0.25, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("normalMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_normalMCCut"), 0., 999., 30, 50, 0.25, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("normalMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_normalMCCut"), 0., 999., 50, 90, 0.25, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("higherMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_higherMCCut"), 0., 999., 0, 10, 0.25, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("higherMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_higherMCCut"), 0., 999., 10, 30, 0.30, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("higherMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_higherMCCut"), 0., 999., 30, 50, 0.30, 0.00, 100., 999.); __R_ADDTASK__->AddJetEmbeddingCut(Form("higherMCCut"), Form("ChargedJetsHadronCF_Particles_Embedded_R030_higherMCCut"), 0., 999., 50, 90, 0.30, 0.00, 100., 999.); /* Int_t vetos[2] = {15, 35}; for(Int_t i=0; i<2; i++) { Int_t jetVetoPt = vetos[i]; // output stream name, cut name, minPt, maxPt, minCent, maxCent, MC percentage low threshold, 2D cut offset, 2D cut slope if(USE_LOW_CUT) { __R_ADDTASK__->AddJetEmbeddingCut(Form("lowerCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_lowerCutVeto%02iGeV", jetVetoPt), 0., 999., 0, 10, 0.15, 0.05, 0.03, jetVetoPt); __R_ADDTASK__->AddJetEmbeddingCut(Form("lowerCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_lowerCutVeto%02iGeV", jetVetoPt), 0., 999., 10, 30, 0.20, 0.02, 0.04, jetVetoPt); __R_ADDTASK__->AddJetEmbeddingCut(Form("lowerCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_lowerCutVeto%02iGeV", jetVetoPt), 0., 999., 30, 50, 0.20, 0.02, 0.06, jetVetoPt); __R_ADDTASK__->AddJetEmbeddingCut(Form("lowerCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_lowerCutVeto%02iGeV", jetVetoPt), 0., 999., 50, 90, 0.20, 0.02, 0.06, jetVetoPt); } if(USE_NORMAL_CUT) { __R_ADDTASK__->AddJetEmbeddingCut(Form("normalCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_normalCutVeto%02iGeV", jetVetoPt), 0., 999., 0, 10, 0.20, 0.05, 0.03, jetVetoPt); __R_ADDTASK__->AddJetEmbeddingCut(Form("normalCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_normalCutVeto%02iGeV", jetVetoPt), 0., 999., 10, 30, 0.25, 0.02, 0.04, jetVetoPt); __R_ADDTASK__->AddJetEmbeddingCut(Form("normalCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_normalCutVeto%02iGeV", jetVetoPt), 0., 999., 30, 50, 0.25, 0.02, 0.06, jetVetoPt); __R_ADDTASK__->AddJetEmbeddingCut(Form("normalCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_normalCutVeto%02iGeV", jetVetoPt), 0., 999., 50, 90, 0.25, 0.02, 0.06, jetVetoPt); } if(USE_HIGH_CUT) { __R_ADDTASK__->AddJetEmbeddingCut(Form("higherCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_higherCutVeto%02iGeV", jetVetoPt), 0., 999., 0, 10, 0.25, 0.05, 0.03, jetVetoPt); __R_ADDTASK__->AddJetEmbeddingCut(Form("higherCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_higherCutVeto%02iGeV", jetVetoPt), 0., 999., 10, 30, 0.30, 0.02, 0.04, jetVetoPt); __R_ADDTASK__->AddJetEmbeddingCut(Form("higherCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_higherCutVeto%02iGeV", jetVetoPt), 0., 999., 30, 50, 0.30, 0.02, 0.06, jetVetoPt); __R_ADDTASK__->AddJetEmbeddingCut(Form("higherCutVeto%02iGeV", jetVetoPt), Form("ChargedJetsHadronCF_Particles_Embedded_R030_higherCutVeto%02iGeV", jetVetoPt), 0., 999., 50, 90, 0.30, 0.02, 0.06, jetVetoPt); } } */ #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_Support_PYTHIA #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskChargedJetsHadronCF.C #Module.MacroArgs kPhysSelRH, "tracks_generated_PY", kJetsAKTR030RHPYTHIA.Data(), "", 0.3, 0.6, kRHMinLeadingHadronPt, 0.15, 0.15, 0.707, "" #Module.Deps ChargedJetsHadronCF_Embedding_JetFinderPY,ChargedJetsHadronCF_Embedding_JetFinder,ChargedJetsHadronCF_Embedding_Rho #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetVzRange(-10,10); __R_ADDTASK__->SelectCollisionCandidates(kPhysSelRH); __R_ADDTASK__->SetJetParticleArrayName("ChargedJetsHadronCF_Particles_PY_R030"); __R_ADDTASK__->SetTrackParticleArrayName("ChargedJetsHadronCF_TrackParticles_PY_R030"); __R_ADDTASK__->GetJetContainer()->SetMaxTrackPt(kRHMaxLeadingHadronPt); __R_ADDTASK__->SetNumberOfCentralityBins(10); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_TaskCorrelations_HigherCut_NoVeto #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_Embedded_HigherCut_NoVeto", "ChargedJetsHadronCF" #Module.Deps ChargedJetsHadronCF_Embedding_Support #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetEventSelectionBit(kPhysSelRH); __R_ADDTASK__->SetCustomParticlesA("ChargedJetsHadronCF_Particles_Embedded_R030_higherMCCut"); __R_ADDTASK__->SetTriggersFromDetector(6); __R_ADDTASK__->SetCustomParticlesB("ChargedJetsHadronCF_TrackParticles_Embedded_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); __R_ADDTASK__->SetCheckEventNumberInMixedEvent(kFALSE); // 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); } __R_ADDTASK__->SetCentralityMethod("nano"); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_TaskCorrelations_LowerCut_NoVeto #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_Embedded_LowerCut_NoVeto", "ChargedJetsHadronCF" #Module.Deps ChargedJetsHadronCF_Embedding_Support #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetEventSelectionBit(kPhysSelRH); __R_ADDTASK__->SetCustomParticlesA("ChargedJetsHadronCF_Particles_Embedded_R030_lowerMCCut"); __R_ADDTASK__->SetTriggersFromDetector(6); __R_ADDTASK__->SetCustomParticlesB("ChargedJetsHadronCF_TrackParticles_Embedded_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); __R_ADDTASK__->SetCheckEventNumberInMixedEvent(kFALSE); // 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); } __R_ADDTASK__->SetCentralityMethod("nano"); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_TaskCorrelations_NoCut #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_Embedded_NoCut", "ChargedJetsHadronCF" #Module.Deps ChargedJetsHadronCF_Embedding_Support #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetEventSelectionBit(kPhysSelRH); __R_ADDTASK__->SetCustomParticlesA("ChargedJetsHadronCF_Particles_Embedded_R030_NoCut"); __R_ADDTASK__->SetTriggersFromDetector(6); __R_ADDTASK__->SetCustomParticlesB("ChargedJetsHadronCF_TrackParticles_Embedded_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); __R_ADDTASK__->SetCheckEventNumberInMixedEvent(kFALSE); // 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 = 60; __R_ADDTASK__->SetMixingTracks(15000); __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); } __R_ADDTASK__->SetCentralityMethod("nano"); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_TaskCorrelations_NormalCut_NoVeto #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_Embedded_NormalCut_NoVeto", "ChargedJetsHadronCF" #Module.Deps ChargedJetsHadronCF_Embedding_Support #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetEventSelectionBit(kPhysSelRH); __R_ADDTASK__->SetCustomParticlesA("ChargedJetsHadronCF_Particles_Embedded_R030_normalMCCut"); __R_ADDTASK__->SetTriggersFromDetector(6); __R_ADDTASK__->SetCustomParticlesB("ChargedJetsHadronCF_TrackParticles_Embedded_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); __R_ADDTASK__->SetCheckEventNumberInMixedEvent(kFALSE); // Load external pool manager if(kRHImportEventPool) { __R_ADDTASK__->SetMixingTracks(15000); TFile* fileInput = TFile::Open(Form("alien:///alice/cern.ch/user/r/rhaake/ExternalEventPool_%d.root", kPtHardBin)); 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); } __R_ADDTASK__->SetCentralityMethod("nano"); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ChargedJetsHadronCF_Embedding_TaskCorrelations_PYTHIA #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_PY", "ChargedJetsHadronCF" #Module.Deps ChargedJetsHadronCF_Embedding_Support_PYTHIA #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetEventSelectionBit(kPhysSelRH); __R_ADDTASK__->SetCustomParticlesA("ChargedJetsHadronCF_Particles_PY_R030"); __R_ADDTASK__->SetTriggersFromDetector(6); __R_ADDTASK__->SetCustomParticlesB("ChargedJetsHadronCF_TrackParticles_PY_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); __R_ADDTASK__->SetCheckEventNumberInMixedEvent(kTRUE); // 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 = 60; __R_ADDTASK__->SetMixingTracks(15000); __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(Form("alien:///alice/cern.ch/user/r/rhaake/ExternalEventPoolPY_%d.root", kPtHardBin)); 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_Embedding_Wrapper_PY #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskParticleRandomizer.C #Module.MacroArgs "", "tracks_generated_PY", kFALSE, kFALSE, kFALSE #Module.Deps ChargedJetsHadronCF_Embedding_Generator #Module.Owner rhaake #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); __R_ADDTASK__->ActivateJetEmbedding("tracks_generated_pico"); // No tracking eff. adaption here -> Was already done in the PYTHIA track "generation" #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(kComPhysSel); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); __R_ADDTASK__->SetTrackEfficiency(kRHTrackEfficiency); #Module.EndConfig