#_______________________________________________________________________________ #Module.Begin CDBConnectTask #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGPP/PilotTrain/AddTaskCDBconnect.C #Module.MacroArgs #Module.Deps #Module.Owner mfasel #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SetFallBackToRaw(kTRUE); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmbeddingHelper_MC #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/ConfigureEmcalEmbeddingHelperOnLEGOTrain.C #Module.MacroArgs #Module.Deps EmbeddingHelper #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig //__R_ADDTASK__->SetPtHardBin(kPtHardBin); __R_ADDTASK__->SetNPtHardBins(11); __R_ADDTASK__->SetFilePattern("alien:///alice/sim/2012/LHC12a15e_fix/169838/%d/AOD149/"); __R_ADDTASK__->SetTriggerMask(0); __R_ADDTASK__->SetMCRejectOutliers(); // Set the pt hard auto config identifier __R_ADDTASK__->SetAutoConfigureIdentifier("13Apr2018Raymond"); // Set additional configuration options in the embedding helper __R_ADDTASK__->SetConfigurationPath("alien:///alice/cern.ch/user/r/rehlersi/embedding/embeddingConfigurationLHC12a15e.yaml"); // Ensure that the internal and external event acceptances are similar (within 4 cm) // Do _not_ use this option - it wastes too many embedded events and never finishes! //__R_ADDTASK__->SetMaxVertexDistance(4); // It is very important to call Initialize() with true! __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_MC_combined #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/ConfigureEmcalCorrectionTaskOnLEGOTrain.C #Module.MacroArgs "combined" #Module.Deps EmcalCorrectionTask_Embedding_CellLevel_data,EmcalCorrectionTask_Embedding_CellLevel_embed,EmcalCorrectionTask_Embedding_ClusterLevel_combined #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/r/rehlersi/embedding/embedMCUserConfiguration.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_MC_data #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/ConfigureEmcalCorrectionTaskOnLEGOTrain.C #Module.MacroArgs "data" #Module.Deps EmcalCorrectionTask_Embedding_CellLevel_data,EmcalCorrectionTask_Embedding_CellLevel_embed,EmcalCorrectionTask_Embedding_ClusterLevel_combined #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/r/rehlersi/embedding/embedMCUserConfiguration.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_MC_embed #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/ConfigureEmcalCorrectionTaskOnLEGOTrain.C #Module.MacroArgs "embed" #Module.Deps EmcalCorrectionTask_Embedding_CellLevel_data,EmcalCorrectionTask_Embedding_CellLevel_embed,EmcalCorrectionTask_Embedding_ClusterLevel_combined #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/r/rehlersi/embedding/embedMCUserConfiguration.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmbeddingHelper #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/AddTaskEmcalEmbeddingHelper.C #Module.MacroArgs #Module.Deps CDBConnectTask,MultSelectionTask #Module.Owner jmulliga #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetRandomFileAccess(kTRUE); // Setup pt hard autoconfiguration // The configure wagon should still set the identifier! __R_ADDTASK__->SetAutoConfigurePtHardBins(true); //__R_ADDTASK__->SetAutoConfigureBasePath("alien:///alice/cern.ch/user/a/alitrain/"); __R_ADDTASK__->SetAutoConfigureBasePath("/home/alitrain/train-workdir/"); __R_ADDTASK__->SetAutoConfigureTrainTypePath("PWGJE/Jets_EMC_PbPb/"); // DO NOT call Initialize() here! #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmcalCorrectionTask_Embedding_CellLevel_data #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/AddTaskEmcalCorrectionTask.C #Module.MacroArgs "data" #Module.Deps EmbeddingHelper #Module.Owner jmulliga #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetUseNewCentralityEstimation(kIsRun2); __R_ADDTASK__->SetNCentBins(5); __R_ADDTASK__->SetForceBeamType(AliEmcalCorrectionTask::kAA); //Do not call Initialize() here! Instead, call it in the configure macro. See the wagon "ConfigureEmcalCorrectionTask_Embedding" //__R_ADDTASK__->Initialize(); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmcalCorrectionTask_Embedding_CellLevel_embed #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/AddTaskEmcalCorrectionTask.C #Module.MacroArgs "embed" #Module.Deps EmbeddingHelper #Module.Owner jmulliga #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetUseNewCentralityEstimation(kIsRun2); __R_ADDTASK__->SetNCentBins(5); __R_ADDTASK__->SetForceBeamType(AliEmcalCorrectionTask::kAA); //Do not call Initialize() here! Instead, call it in the configure macro. See the wagon "ConfigureEmcalCorrectionTask_Embedding" //__R_ADDTASK__->Initialize(); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmcalCorrectionTask_Embedding_ClusterLevel_combined #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/AddTaskEmcalCorrectionTask.C #Module.MacroArgs "combined" #Module.Deps EmcalCorrectionTask_Embedding_CellLevel_data,EmcalCorrectionTask_Embedding_CellLevel_embed #Module.Owner jmulliga #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetUseNewCentralityEstimation(kIsRun2); __R_ADDTASK__->SetNCentBins(5); __R_ADDTASK__->SetForceBeamType(AliEmcalCorrectionTask::kAA); //Do not call Initialize() here! Instead, call it in the configure macro. See the wagon "ConfigureEmcalCorrectionTask_Embedding" //__R_ADDTASK__->Initialize(); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmcalJetSampleTask_detLevel_Embedding_MC #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetSample.C #Module.MacroArgs "usedefault", "usedefault", "", "detLevel" #Module.Deps JetFinder_full_AKT_02_3GeVConstituent_detLevelJets_Embedding_MC #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetHistoBins(600, 0, 300); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); // Embedded tracks AliTrackContainer * tracksDetLevel = __R_ADDTASK__->GetTrackContainer(0); tracksDetLevel->SetIsEmbedding(kTRUE); tracksDetLevel->SetParticlePtCut(0.15); // Cluster containers __R_ADDTASK__->GetClusterContainer(0)->SetIsEmbedding(kTRUE); __R_ADDTASK__->GetClusterContainer(0)->SetClusECut(0.); __R_ADDTASK__->GetClusterContainer(0)->SetClusPtCut(0.); __R_ADDTASK__->GetClusterContainer(0)->SetClusNonLinCorrEnergyCut(0.); __R_ADDTASK__->GetClusterContainer(0)->SetClusHadCorrEnergyCut(0.30); __R_ADDTASK__->GetClusterContainer(0)->SetDefaultClusterEnergy(AliVCluster::kHadCorr); __R_ADDTASK__->SetHistoBins(600, 0, 300); // Jets __R_ADDTASK__->AddJetContainer("detLevelJets_AKTFullR020_tracks_pT3000_caloClusters_E3000_pt_scheme", AliEmcalJet::kEMCALfid, 0.2); //__R_ADDTASK__->AddJetContainer(AliJetContainer::kFullJet, AliJetContainer::antikt_algorithm, AliJetContainer::pt_scheme, 0.2, AliEmcalJet::kEMCALfid, "detLevelJets"); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmcalJetSampleTask_hybridLevel_Embedding_MC #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetSample.C #Module.MacroArgs "usedefault", "caloClustersCombined", "emcalCellsCombined", "hybridLevel" #Module.Deps JetFinder_full_AKT_02_3GeVConstituent_hybridLevelJets_Embedding_MC #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetHistoBins(600, 0, 300); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); // Particle containers __R_ADDTASK__->GetParticleContainer(0)->SetParticlePtCut(0.15); // Embedded tracks AliTrackContainer * tracksDetLevel = new AliTrackContainer("tracks"); tracksDetLevel->SetIsEmbedding(kTRUE); tracksDetLevel->SetParticlePtCut(0.15); __R_ADDTASK__->AdoptTrackContainer(tracksDetLevel); // Cluster containers __R_ADDTASK__->GetClusterContainer(0)->SetClusECut(0.); __R_ADDTASK__->GetClusterContainer(0)->SetClusPtCut(0.); __R_ADDTASK__->GetClusterContainer(0)->SetClusNonLinCorrEnergyCut(0.); __R_ADDTASK__->GetClusterContainer(0)->SetClusHadCorrEnergyCut(0.30); __R_ADDTASK__->GetClusterContainer(0)->SetDefaultClusterEnergy(AliVCluster::kHadCorr); __R_ADDTASK__->SetHistoBins(600, 0, 300); // Jets __R_ADDTASK__->AddJetContainer("hybridLevelJets_AKTFullR020_tracks_pT3000_caloClustersCombined_E3000_pt_scheme", AliEmcalJet::kEMCALfid, 0.2); //__R_ADDTASK__->AddJetContainer(AliJetContainer::kFullJet, AliJetContainer::antikt_algorithm, AliJetContainer::pt_scheme, 0.2, AliEmcalJet::kEMCALfid, "hybridLevelJets"); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmcalJetSampleTask_partLevel_Embedding_MC #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetSample.C #Module.MacroArgs "mcparticles", "", "", "partLevel" #Module.Deps JetFinder_full_AKT_02_3GeVConstituent_partLevelJets_Embedding_MC #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetHistoBins(600, 0, 300); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); // Particle containers AliMCParticleContainer * partCont = dynamic_cast(__R_ADDTASK__->GetParticleContainer(0)); partCont->SetIsEmbedding(kTRUE); partCont->SetParticlePtCut(0.15); // Jets __R_ADDTASK__->AddJetContainer("partLevelJets_AKTFullR020_mcparticles_pT3000_pt_scheme", AliEmcalJet::kEMCALfid, 0.2); //__R_ADDTASK__->AddJetContainer(AliJetContainer::kFullJet, AliJetContainer::antikt_algorithm, AliJetContainer::pt_scheme, 0.2, AliEmcalJet::kEMCALfid, "partLevelJets"); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetFinder_full_AKT_02_3GeVConstituent_detLevelJets_Embedding_MC #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "usedefault", "usedefault", AliJetContainer::antikt_algorithm, 0.2, AliJetContainer::kFullJet, 3.0, 3.0, kGhostArea, AliJetContainer::pt_scheme, "detLevelJets", 0., kFALSE, kFALSE #Module.Deps EmcalCorrectionTask_Embedding_ClusterLevel_combined #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); // Setup the tracks properly to be retrieved from the external event AliParticleContainer * tracks = __R_ADDTASK__->GetParticleContainer(0); tracks->SetIsEmbedding(kTRUE); AliClusterContainer * clusters = __R_ADDTASK__->GetClusterContainer(0); clusters->SetIsEmbedding(kTRUE); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetFinder_full_AKT_02_3GeVConstituent_hybridLevelJets_Embedding_MC #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "usedefault", "caloClustersCombined", AliJetContainer::antikt_algorithm, 0.2, AliJetContainer::kFullJet, 3.0, 3.0, kGhostArea, AliJetContainer::pt_scheme, "hybridLevelJets", 0., kFALSE, kFALSE #Module.Deps EmcalCorrectionTask_Embedding_ClusterLevel_combined #Module.Owner rehlersi #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); // Artificial tracking efficiency derived from the full jet R_AA paper. // 10-30 \approx .99 // 0-10% \approx .98 // Apply only to embedded particles __R_ADDTASK__->SetTrackEfficiency(.98); __R_ADDTASK__->SetTrackEfficiencyOnlyForEmbedding(kTRUE); // Particle containers AliTrackContainer * tracksDetLevel = new AliTrackContainer("tracks"); // Get the det level tracks from the external event! tracksDetLevel->SetIsEmbedding(kTRUE); tracksDetLevel->SetParticlePtCut(3.); __R_ADDTASK__->AdoptTrackContainer(tracksDetLevel); // Cluster containers __R_ADDTASK__->GetClusterContainer(0)->SetDefaultClusterEnergy(AliVCluster::kHadCorr); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetFinder_full_AKT_02_3GeVConstituent_partLevelJets_Embedding_MC #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "mcparticles", "", AliJetContainer::antikt_algorithm, 0.2, AliJetContainer::kFullJet, 3.0, 3.0, kGhostArea, AliJetContainer::pt_scheme, "partLevelJets", 0., kFALSE, kFALSE #Module.Deps EmbeddingHelper #Module.Owner rehlersi #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); // Setup the tracks properly to be retrieved from the external event AliMCParticleContainer * truthTracks = __R_ADDTASK__->GetMCParticleContainer(0); truthTracks->SetIsEmbedding(kTRUE); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetResponseMaker_Embedding_MC_Cent010_trackBias #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskJetResponseMaker.C #Module.MacroArgs "tracks", "caloClustersCombined", "hybridLevelJets_AKTFullR020_tracks_pT3000_caloClustersCombined_E3000_pt_scheme", "", 0.2, "mcparticles", "", "partLevelJets_AKTFullR020_mcparticles_pT3000_pt_scheme", "", 0.2, 0., 0., 0., 0, AliJetResponseMaker::kGeometrical, 1.2, 1.2, "EMCALfid", -999, 0, 10 #Module.Deps JetFinder_full_AKT_02_3GeVConstituent_partLevelJets_Embedding_MC,JetFinder_full_AKT_02_3GeVConstituent_hybridLevelJets_Embedding_MC #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig ///////////// // Note: kPtHardBin in the arguments should _not_ be set in the arguments!! ///////////// __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); // Configure __R_ADDTASK__->SetHistoType(1); // THnSparse __R_ADDTASK__->SetJetRelativeEPAngleAxis(kTRUE); // Jets // Jet-Hadron jet area cut (det level) // 60% of R=0.2 jets __R_ADDTASK__->GetJetContainer(0)->SetJetAreaCut(0.075); // This will crash otherwise, as it will attempt to access this particle containers (although it will not use it for anything) __R_ADDTASK__->GetJetContainer(1)->GetParticleContainer()->SetIsEmbedding(kTRUE); // Remove a large number of trivial minimal pt jets which will pollute the matching __R_ADDTASK__->GetJetContainer(1)->SetJetPtCut(0.15); // Disable pythia support in AliAnalysisTaskEmcal since it interferes with AliEmcalEmbeddingQA __R_ADDTASK__->SetIsPythia(kFALSE); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetResponseMaker_Embedding_MC_Cent3050_trackBias #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskJetResponseMaker.C #Module.MacroArgs "tracks", "caloClustersCombined", "hybridLevelJets_AKTFullR020_tracks_pT3000_caloClustersCombined_E3000_pt_scheme", "", 0.2, "mcparticles", "", "partLevelJets_AKTFullR020_mcparticles_pT3000_pt_scheme", "", 0.2, 0., 0., 0., 0, AliJetResponseMaker::kGeometrical, 1.2, 1.2, "EMCALfid", -999, 30, 50 #Module.Deps JetFinder_full_AKT_02_3GeVConstituent_partLevelJets_Embedding_MC,JetFinder_full_AKT_02_3GeVConstituent_hybridLevelJets_Embedding_MC #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig ///////////// // Note: kPtHardBin in the arguments should _not_ be set in the arguments!! ///////////// __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); // Configure __R_ADDTASK__->SetHistoType(1); // THnSparse __R_ADDTASK__->SetJetRelativeEPAngleAxis(kTRUE); // Jets // Jet-Hadron jet area cut (det level) // 60% of R=0.2 jets __R_ADDTASK__->GetJetContainer(0)->SetJetAreaCut(0.075); // This will crash otherwise, as it will attempt to access this particle containers (although it will not use it for anything) __R_ADDTASK__->GetJetContainer(1)->GetParticleContainer()->SetIsEmbedding(kTRUE); // Remove a large number of trivial minimal pt jets which will pollute the matching __R_ADDTASK__->GetJetContainer(1)->SetJetPtCut(0.15); // Disable pythia support in AliAnalysisTaskEmcal since it interferes with AliEmcalEmbeddingQA __R_ADDTASK__->SetIsPythia(kFALSE); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetResponseMaker_Embedding_MC_Cent3050_trackBias_detLevel #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskJetResponseMaker.C #Module.MacroArgs "tracks", "caloClusters", "detLevelJets_AKTFullR020_tracks_pT3000_caloClusters_E3000_pt_scheme", "", 0.2, "mcparticles", "", "partLevelJets_AKTFullR020_mcparticles_pT3000_pt_scheme", "", 0.2, 0., 0., 0., 0, AliJetResponseMaker::kGeometrical, 1.2, 1.2, "EMCALfid", -999, 30, 50 #Module.Deps JetFinder_full_AKT_02_3GeVConstituent_partLevelJets_Embedding_MC,JetFinder_full_AKT_02_3GeVConstituent_detLevelJets_Embedding_MC #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig ///////////// // Note: kPtHardBin in the arguments should _not_ be set in the arguments!! ///////////// __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); // Configure __R_ADDTASK__->SetHistoType(1); // THnSparse __R_ADDTASK__->SetJetRelativeEPAngleAxis(kTRUE); // Jets // Jet-Hadron jet area cut (det level) // 60% of R=0.2 jets __R_ADDTASK__->GetJetContainer(0)->SetJetAreaCut(0.075); // This will crash otherwise, as it will attempt to access this particle containers (although it will not use it for anything) __R_ADDTASK__->GetJetContainer(0)->GetParticleContainer()->SetIsEmbedding(kTRUE); __R_ADDTASK__->GetJetContainer(0)->GetClusterContainer()->SetIsEmbedding(kTRUE); __R_ADDTASK__->GetJetContainer(1)->GetParticleContainer()->SetIsEmbedding(kTRUE); // Remove a large number of trivial minimal pt jets which will pollute the matching __R_ADDTASK__->GetJetContainer(1)->SetJetPtCut(0.15); // Disable pythia support in AliAnalysisTaskEmcal since it interferes with AliEmcalEmbeddingQA __R_ADDTASK__->SetIsPythia(kFALSE); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin MultSelectionTask #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/OADB/COMMON/MULTIPLICITY/macros/AddTaskMultSelection.C #Module.MacroArgs kFALSE #Module.Deps CDBConnectTask #Module.Owner jmazer #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); /* Actually the multiplicity task does need to run on AODs as well if (strcmp(kDatatype,"AOD")==0) { // AOD data sets don't need centrality task cout << "*** Removing centrality task because AOD dataset" << endl; AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); TObjArray *array = mgr->GetTasks(); array->Remove(__R_ADDTASK__); } */ #Module.EndConfig