#_______________________________________________________________________________ #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_Hannah #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/ConfigureEmcalEmbeddingHelperOnLEGOTrain.C #Module.MacroArgs #Module.Deps EmbeddingHelper #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); // Use to configure test of a specific pT Hard bin //__R_ADDTASK__->SetAutoConfigurePtHardBins(false); //__R_ADDTASK__->SetPtHardBin(11); __R_ADDTASK__->SetNPtHardBins(20); __R_ADDTASK__->SetFilePattern("alien:///alice/sim/2016/LHC16j5/%d/"); //__R_ADDTASK__->SetInputFilename("*/AOD/*/aod_archive.zip"); __R_ADDTASK__->SetInputFilename("*/AOD200/*/AliAOD.root"); //__R_ADDTASK__->SetInputFilename("*/AOD/*/AliAOD.root"); __R_ADDTASK__->SetTriggerMask(AliVEvent::kAny); // Equivalent to 0 __R_ADDTASK__->SetRandomFileAccess(kTRUE); __R_ADDTASK__->SetRandomEventNumberAccess(kTRUE); __R_ADDTASK__->SetMCRejectOutliers(); __R_ADDTASK__->SetPtHardJetPtRejectionFactor(4.); //__R_ADDTASK__->SetZVertexCut(10.); // Set YAML configuration, including good runlist of embedded events, and internal event selection feature __R_ADDTASK__->SetConfigurationPath("alien:///alice/cern.ch/user/r/rhaake/embeddingHelper_LHC15o_LHC16j5.yaml"); // Set the pt hard auto config identifier __R_ADDTASK__->SetAutoConfigureIdentifier("20190713bHannah"); // 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_Hannah_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 hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/h/hbossi/emcalCorrections_LHC15o_LHC16j5.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_Hannah_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 hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/h/hbossi/emcalCorrections_LHC15o_LHC16j5.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_Hannah_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 hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/h/hbossi/emcalCorrections_LHC15o_LHC16j5.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 JetExtractor_hybridLevel_full_02 #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskJetExtractor.C #Module.MacroArgs "tracks", "caloClustersCombined", "hybridLevelJets_AKTFullR020_tracks_pT0150_caloClustersCombined_E0300_pt_scheme", "Rho_Scaled", 0.2, "hybridLevelJets" #Module.Deps JetExtractor_Rho,JetExtractor_JetFinder_full_AKT_02_DetLevel,JetExtractor_JetFinder_full_AKT_02_HybridLevel,JetExtractor_JetFinder_full_AKT_02_PartLevel #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetNeedEmcalGeom(kTRUE); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(kTRUE); __R_ADDTASK__->SetSaveMCInformation(kTRUE); __R_ADDTASK__->SetSaveCaloClusters(kTRUE); __R_ADDTASK__->SetSaveConstituents(kTRUE); __R_ADDTASK__->SetSaveJetShapes(kTRUE); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(-20,0, 0.005); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(0,10, 0.007); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(10,20, 0.008); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(20,40, 0.05); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(40,60, 0.08); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(60,80, 0.15); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(80,100, 0.3); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(100,150, 0.4); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(150,200, 0.5); __R_ADDTASK__->GetClusterContainer(0)->SetIsEmbedding(kFALSE); __R_ADDTASK__->GetJetContainer(0)->SetJetAcceptanceType(AliEmcalJet::kEMCALfid); // add another cluster container to the event to handle the hybrid case __R_ADDTASK__->AddClusterContainer("caloClusters"); __R_ADDTASK__->GetClusterContainer(1)->SetIsEmbedding(kTRUE); __R_ADDTASK__->GetClusterContainer(1)->SetClusECut(0.); __R_ADDTASK__->GetClusterContainer(1)->SetClusPtCut(0.); __R_ADDTASK__->GetClusterContainer(1)->SetClusNonLinCorrEnergyCut(0.); __R_ADDTASK__->GetClusterContainer(1)->SetClusHadCorrEnergyCut(0.30); __R_ADDTASK__->GetClusterContainer(1)->SetDefaultClusterEnergy(AliVCluster::kHadCorr); __R_ADDTASK__->SetNeedEmbedClusterContainer(kTRUE); __R_ADDTASK__->ActivateJetMatching("detLevelJets_AKTFullR020_tracks_pT0150_caloClusters_E0300_pt_scheme", "partLevelJets_AKTFullR020_mcparticles_pT0150_pt_scheme", "", ""); // add pythia tracks AliParticleContainer * mcpartCont = new AliTrackContainer("tracks"); mcpartCont->SetIsEmbedding(kTRUE); __R_ADDTASK__->AdoptParticleContainer(mcpartCont); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetExtractor_hybridLevel_full_04 #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskJetExtractor.C #Module.MacroArgs "tracks", "caloClustersCombined", "hybridLevelJets_AKTFullR040_tracks_pT0150_caloClustersCombined_E0300_pt_scheme", "Rho_Scaled", 0.4, "hybridLevelJets" #Module.Deps JetExtractor_Rho,JetExtractor_JetFinder_full_AKT_04_DetLevel,JetExtractor_JetFinder_full_AKT_04_HybridLevel,JetExtractor_JetFinder_full_AKT_04_PartLevel,JetExtractor_Rho_JF #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetNeedEmcalGeom(kTRUE); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(kTRUE); __R_ADDTASK__->SetSaveMCInformation(kTRUE); __R_ADDTASK__->SetSaveCaloClusters(kTRUE); __R_ADDTASK__->SetSaveConstituents(kTRUE); __R_ADDTASK__->SetSaveJetShapes(kTRUE); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(-20,0, 0.004); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(0,10, 0.004); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(10,20, 0.04); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(20,40, 0.1); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(40,60, 0.2); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(60,80, 0.4); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(80,100, 0.5); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(100,150, 0.8); __R_ADDTASK__->GetJetTree()->AddExtractionPercentage(150,200, 1.0); __R_ADDTASK__->GetClusterContainer(0)->SetIsEmbedding(kFALSE); __R_ADDTASK__->GetJetContainer(0)->SetJetAcceptanceType(AliEmcalJet::kEMCALfid); // add another cluster container to the event to handle the hybrid case __R_ADDTASK__->AddClusterContainer("caloClusters"); __R_ADDTASK__->GetClusterContainer(1)->SetIsEmbedding(kTRUE); __R_ADDTASK__->GetClusterContainer(1)->SetClusECut(0.); __R_ADDTASK__->GetClusterContainer(1)->SetClusPtCut(0.); __R_ADDTASK__->GetClusterContainer(1)->SetClusNonLinCorrEnergyCut(0.); __R_ADDTASK__->GetClusterContainer(1)->SetClusHadCorrEnergyCut(0.30); __R_ADDTASK__->GetClusterContainer(1)->SetDefaultClusterEnergy(AliVCluster::kHadCorr); __R_ADDTASK__->SetNeedEmbedClusterContainer(kTRUE); __R_ADDTASK__->ActivateJetMatching("detLevelJets_AKTFullR040_tracks_pT0150_caloClusters_E0300_pt_scheme", "partLevelJets_AKTFullR040_mcparticles_pT0150_pt_scheme", "", ""); // add pythia tracks AliParticleContainer * mcpartCont = new AliTrackContainer("tracks"); mcpartCont->SetIsEmbedding(kTRUE); __R_ADDTASK__->AdoptParticleContainer(mcpartCont); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetExtractor_JetFinder_full_AKT_02_DetLevel #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, kTrackPtCut, kClusPtCut, kGhostArea, AliJetContainer::pt_scheme, "detLevelJets", 0., kFALSE, kFALSE #Module.Deps ConfigureEmcalCorrectionTask_Embedding_Hannah_combined,ConfigureEmcalCorrectionTask_Embedding_Hannah_data,ConfigureEmcalCorrectionTask_Embedding_Hannah_embed #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(kTRUE); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); AliTrackContainer * tracksDetLevel = __R_ADDTASK__->GetTrackContainer(0); // Get the det level tracks from the external event! tracksDetLevel->SetIsEmbedding(kTRUE); // Uses the name of the container passed into AliEmcalJetTask AliClusterContainer * clustersDetLevel =__R_ADDTASK__->GetClusterContainer(0); // Get the det level clusters from the external event! clustersDetLevel->SetIsEmbedding(kTRUE); // Additional configuration clustersDetLevel->SetClusTimeCut(-50e-9,100e-9); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetExtractor_JetFinder_full_AKT_02_HybridLevel #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "tracks", "caloClustersCombined", AliJetContainer::antikt_algorithm, 0.2, AliJetContainer::kFullJet, kTrackPtCut, kClusPtCut, kGhostArea, AliJetContainer::pt_scheme, "hybridLevelJets", -1000., kFALSE, kFALSE #Module.Deps ConfigureEmcalCorrectionTask_Embedding_Hannah_combined,ConfigureEmcalCorrectionTask_Embedding_Hannah_data,ConfigureEmcalCorrectionTask_Embedding_Hannah_embed #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(kTRUE); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); // artificial tracking efficiency __R_ADDTASK__->SetTrackEfficiency(.99); __R_ADDTASK__->SetTrackEfficiencyOnlyForEmbedding(kTRUE); __R_ADDTASK__->GetClusterContainer(0)->SetDefaultClusterEnergy(AliVCluster::kHadCorr); __R_ADDTASK__->GetClusterContainer(0)->SetClusTimeCut(-50e-9,100e-9); // Cells already combined, we don't need to take clusters from the external event /////// // External event (ie embedding) settings for PbPb jet finding (adds detector level PYTHIA) /////// // Add embedded tracks and clusters to jet finder // Tracks AliTrackContainer * tracksEmbedDetLevel = new AliTrackContainer("tracks"); // Get the det level tracks from the external event! tracksEmbedDetLevel->SetIsEmbedding(kTRUE); tracksEmbedDetLevel->SetParticlePtCut(0.15); __R_ADDTASK__->AdoptTrackContainer(tracksEmbedDetLevel); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetExtractor_JetFinder_full_AKT_02_PartLevel #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, kTrackPtCut, kClusPtCut, kGhostArea, AliJetContainer::pt_scheme, "partLevelJets", 0., kFALSE, kFALSE #Module.Deps ConfigureEmcalCorrectionTask_Embedding_Hannah_combined,ConfigureEmcalCorrectionTask_Embedding_Hannah_data,ConfigureEmcalCorrectionTask_Embedding_Hannah_embed #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(kTRUE); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); AliParticleContainer * partLevelTracks = __R_ADDTASK__->GetParticleContainer(0); //Called Embedded, but really just means get from an external event! partLevelTracks->SetIsEmbedding(kTRUE); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetExtractor_JetFinder_full_AKT_04_DetLevel #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "usedefault", "usedefault", AliJetContainer::antikt_algorithm, 0.4, AliJetContainer::kFullJet, kTrackPtCut, kClusPtCut, kGhostArea, AliJetContainer::pt_scheme, "detLevelJets", 0., kFALSE, kFALSE #Module.Deps ConfigureEmcalCorrectionTask_Embedding_Hannah_combined,ConfigureEmcalCorrectionTask_Embedding_Hannah_data,ConfigureEmcalCorrectionTask_Embedding_Hannah_embed #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(kTRUE); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); AliTrackContainer * tracksDetLevel = __R_ADDTASK__->GetTrackContainer(0); // Get the det level tracks from the external event! tracksDetLevel->SetIsEmbedding(kTRUE); // Uses the name of the container passed into AliEmcalJetTask AliClusterContainer * clustersDetLevel =__R_ADDTASK__->GetClusterContainer(0); // Get the det level clusters from the external event! clustersDetLevel->SetIsEmbedding(kTRUE); // Additional configuration clustersDetLevel->SetClusTimeCut(-50e-9,100e-9); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetExtractor_JetFinder_full_AKT_04_HybridLevel #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "tracks", "caloClustersCombined", AliJetContainer::antikt_algorithm, 0.4, AliJetContainer::kFullJet, kTrackPtCut, kClusPtCut, kGhostArea, AliJetContainer::pt_scheme, "hybridLevelJets", -1000., kFALSE, kFALSE #Module.Deps ConfigureEmcalCorrectionTask_Embedding_Hannah_combined,ConfigureEmcalCorrectionTask_Embedding_Hannah_data,ConfigureEmcalCorrectionTask_Embedding_Hannah_embed #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(kTRUE); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); // artificial tracking efficiency __R_ADDTASK__->SetTrackEfficiency(.99); __R_ADDTASK__->SetTrackEfficiencyOnlyForEmbedding(kTRUE); __R_ADDTASK__->GetClusterContainer(0)->SetDefaultClusterEnergy(AliVCluster::kHadCorr); __R_ADDTASK__->GetClusterContainer(0)->SetClusTimeCut(-50e-9,100e-9); // Cells already combined, we don't need to take clusters from the external event /////// // External event (ie embedding) settings for PbPb jet finding (adds detector level PYTHIA) /////// // Add embedded tracks and clusters to jet finder // Tracks AliTrackContainer * tracksEmbedDetLevel = new AliTrackContainer("tracks"); // Get the det level tracks from the external event! tracksEmbedDetLevel->SetIsEmbedding(kTRUE); tracksEmbedDetLevel->SetParticlePtCut(0.15); __R_ADDTASK__->AdoptTrackContainer(tracksEmbedDetLevel); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetExtractor_JetFinder_full_AKT_04_PartLevel #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "mcparticles", "", AliJetContainer::antikt_algorithm, 0.4, AliJetContainer::kFullJet, kTrackPtCut, kClusPtCut, kGhostArea, AliJetContainer::pt_scheme, "partLevelJets", 0., kFALSE, kFALSE #Module.Deps ConfigureEmcalCorrectionTask_Embedding_Hannah_combined,ConfigureEmcalCorrectionTask_Embedding_Hannah_data,ConfigureEmcalCorrectionTask_Embedding_Hannah_embed #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(kTRUE); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); AliParticleContainer * partLevelTracks = __R_ADDTASK__->GetParticleContainer(0); //Called Embedded, but really just means get from an external event! partLevelTracks->SetIsEmbedding(kTRUE); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetExtractor_Rho #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskRhoNew.C #Module.MacroArgs "tracks", "", "Rho", 0.2, AliJetContainer::kTPCfid, AliJetContainer::kChargedJet, kTRUE, AliJetContainer::E_scheme, "Rho_ExLJ" #Module.Deps JetExtractor_Rho_JF #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetVzRange(-10,10); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetExcludeLeadJets(2); __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(kTRUE); TString sFuncPath = "alien:///alice/cern.ch/user/j/jmulliga/scaleFactorEMCalLHC15o_PtDepTrackMatching.root"; TString sFuncName = "fScaleFactorEMCal"; __R_ADDTASK__->LoadRhoFunction(sFuncPath, sFuncName); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetExtractor_Rho_JF #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "tracks","", AliJetContainer::kt_algorithm, 0.2, AliJetContainer::kChargedJet, 0.15, 0.3, 0.005, AliJetContainer::E_scheme, "Jet", 0., kFALSE, kFALSE #Module.Deps ConfigureEmcalCorrectionTask_Embedding_Hannah_combined,ConfigureEmcalCorrectionTask_Embedding_Hannah_data,ConfigureEmcalCorrectionTask_Embedding_Hannah_embed #Module.Owner hbossi #Module.OutputFile AnalysisResults.root #Module.StartConfig AliTrackContainer* trackCont = new AliTrackContainer("tracks"); trackCont->SetIsEmbedding(kTRUE); __R_ADDTASK__->AdoptParticleContainer(trackCont); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); #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