#_______________________________________________________________________________ #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_jetHadron_response #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__->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("26Feb2019Raymond"); // Set additional configuration options in the embedding helper __R_ADDTASK__->SetConfigurationPath("alien:///alice/cern.ch/user/r/rehlersi/config/embedding/embeddingHelper_LHC11h_LHC12a15e_semiCentralResponse.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_jetHadron_response_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/config/embedding/emcalCorrections_LHC11h_LHC12a15e.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_MC_jetHadron_response_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/config/embedding/emcalCorrections_LHC11h_LHC12a15e.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_MC_jetHadron_response_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/config/embedding/emcalCorrections_LHC11h_LHC12a15e.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 EmcalJetTagger_Embedding_MC_DetToPartLevel #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetTaggerFast.C #Module.MacroArgs "detLevelJets_AKTFullR020_tracks_pT3000_caloClusters_E3000_pt_scheme", "partLevelJets_AKTFullR020_mcparticles_pT3000_pt_scheme", 0.2, "", "", "", "", "EMCALfid", "V0M", kComPhysSel, "" #Module.Deps EmcalJetTagger_Embedding_MC_hybridToDetLevel,JetFinder_full_AKT_02_3GeVConstituent_detLevelJets_Embedding_MC,JetFinder_full_AKT_02_3GeVConstituent_partLevelJets_Embedding_MC #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig // Task level settings __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); // Tag via geometrical matching __R_ADDTASK__->SetJetTaggingMethod(PWGJE::EMCALJetTasks::AliEmcalJetTaggerTaskFast::kGeo); // Tag the closest jet __R_ADDTASK__->SetJetTaggingType(PWGJE::EMCALJetTasks::AliEmcalJetTaggerTaskFast::kClosest); // Don't impose any additional acceptance cuts beyond the jet containers __R_ADDTASK__->SetTypeAcceptance(PWGJE::EMCALJetTasks::AliEmcalJetTaggerTaskFast::kNoLimit); // Use default matching distance __R_ADDTASK__->SetMaxDistance(0.3); // Redundant, but done for completeness __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); // Reapply the max track pt cut off to maintain energy resolution and avoid fake tracks // However, don't apply to the particle level jets which don't suffer this effect auto detLevelJetCont = __R_ADDTASK__->GetJetContainer(0); detLevelJetCont->SetMaxTrackPt(100); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmcalJetTagger_Embedding_MC_hybridToDetLevel #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetTaggerFast.C #Module.MacroArgs "hybridLevelJets_AKTFullR020_tracks_pT3000_caloClustersCombined_E3000_pt_scheme", "detLevelJets_AKTFullR020_tracks_pT3000_caloClusters_E3000_pt_scheme", 0.2, "", "", "", "", "EMCALfid", "V0M", kComPhysSel, "" #Module.Deps JetFinder_full_AKT_02_3GeVConstituent_detLevelJets_Embedding_MC,JetFinder_full_AKT_02_3GeVConstituent_hybridLevelJets_Embedding_MC #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig // Task level settings __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); // Tag via geometrical matching __R_ADDTASK__->SetJetTaggingMethod(PWGJE::EMCALJetTasks::AliEmcalJetTaggerTaskFast::kGeo); // Tag the closest jet __R_ADDTASK__->SetJetTaggingType(PWGJE::EMCALJetTasks::AliEmcalJetTaggerTaskFast::kClosest); // Don't impose any additional acceptance cuts beyond the jet containers __R_ADDTASK__->SetTypeAcceptance(PWGJE::EMCALJetTasks::AliEmcalJetTaggerTaskFast::kNoLimit); // Use default matching distance __R_ADDTASK__->SetMaxDistance(0.3); // Redundant, but done for completeness __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); // Reapply the max track pt cut off to maintain energy resolution and avoid fake tracks auto hybridJetCont = __R_ADDTASK__->GetJetContainer(0); hybridJetCont->SetMaxTrackPt(100); auto detLevelJetCont = __R_ADDTASK__->GetJetContainer(1); detLevelJetCont->SetMaxTrackPt(100); #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); clusters->SetClusTimeCut(-50e-9, 100e-9); // 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(.99); __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); __R_ADDTASK__->GetClusterContainer(0)->SetClusTimeCut(-50e-9, 100e-9); // 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 JetHPerformance_Embedding_MC_response #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetHPerformance.C #Module.MacroArgs #Module.Deps EmcalJetTagger_Embedding_MC_DetToPartLevel,EmcalJetTagger_Embedding_MC_hybridToDetLevel #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->AddConfigurationFile("alien:///alice/cern.ch/user/r/rehlersi/config/embedding/jetHPerformance_LHC11h_base.yaml", "base"); __R_ADDTASK__->AddConfigurationFile("alien:///alice/cern.ch/user/r/rehlersi/config/embedding/jetHPerformance_LHC11h_LHC12a15e.yaml", "periodSpecific"); __R_ADDTASK__->AddConfigurationFile("alien:///alice/cern.ch/user/r/rehlersi/config/embedding/jetHPerformance_LHC11h_override_semiCentral.yaml", "semiCentralOverride"); __R_ADDTASK__->Initialize(); // Just to be certain __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