#_______________________________________________________________________________ #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 ConfigureEmcalCorrectionTask_Embedding_James_combined #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/ConfigureEmcalCorrectionTaskOnLEGOTrain.C #Module.MacroArgs "combined" #Module.Deps EmcalCorrectionTask_Embedding_combined,EmcalCorrectionTask_Embedding_data,EmcalCorrectionTask_Embedding_embed #Module.Owner jmulliga #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/j/jmulliga/userJetEmbedding_50.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_James_data #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/ConfigureEmcalCorrectionTaskOnLEGOTrain.C #Module.MacroArgs "data" #Module.Deps EmcalCorrectionTask_Embedding_combined,EmcalCorrectionTask_Embedding_data,EmcalCorrectionTask_Embedding_embed #Module.Owner jmulliga #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/j/jmulliga/userJetEmbedding_50.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_James_embed #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/ConfigureEmcalCorrectionTaskOnLEGOTrain.C #Module.MacroArgs "embed" #Module.Deps EmcalCorrectionTask_Embedding_combined,EmcalCorrectionTask_Embedding_data,EmcalCorrectionTask_Embedding_embed #Module.Owner jmulliga #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/j/jmulliga/userJetEmbedding_50.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmbeddingHelper_Run2 #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__->SetPtHardBin(kPtHardBin); __R_ADDTASK__->SetNPtHardBins(20); __R_ADDTASK__->SetFilePattern("alien:///alice/sim/2016/LHC16j5/%d/"); //__R_ADDTASK__->SetFilePattern("alien:///alice/sim/2017/LHC17e2/"); //__R_ADDTASK__->SetInputFilename("*/AOD/*/aod_archive.zip"); __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 good runlist of embedded events __R_ADDTASK__->SetConfigurationPath("alien:///alice/cern.ch/user/j/jmulliga/LHC16j5_GoodRunlistEMCal.yaml"); __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/"); __R_ADDTASK__->SetAutoConfigureIdentifier("18Jan2018James"); __R_ADDTASK__->Initialize(); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmcalCorrectionTask_Embedding_combined #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/AddTaskEmcalCorrectionTask.C #Module.MacroArgs "combined" #Module.Deps EmbeddingHelper_Run2 #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_data #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/AddTaskEmcalCorrectionTask.C #Module.MacroArgs "data" #Module.Deps EmbeddingHelper_Run2 #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_embed #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/AddTaskEmcalCorrectionTask.C #Module.MacroArgs "embed" #Module.Deps EmbeddingHelper_Run2 #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 JetFinder_charged_KT_04_Embedding #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "usedefault", "", AliJetContainer::kt_algorithm, 0.4, AliJetContainer::kChargedJet, kTrackPtCut, kClusPtCut, kGhostArea, AliJetContainer::pt_scheme, "Jet", 0., kFALSE, kFALSE #Module.Deps MultSelectionTask #Module.Owner jmulliga #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetNCentBins(5); __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); AliTrackContainer* tracksDetLevel = new AliTrackContainer("tracks"); tracksDetLevel->SetIsEmbedding(kTRUE); tracksDetLevel->SetParticlePtCut(0.15); __R_ADDTASK__->AdoptTrackContainer(tracksDetLevel); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetFinder_full_AKT_02_Embedding #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, kTrackPtCut, kClusPtCut, kGhostArea, AliJetContainer::pt_scheme, "Jet", 0.1, kFALSE, kFALSE #Module.Deps EmcalCorrectionTask_Embedding_combined,EmcalCorrectionTask_Embedding_data,EmcalCorrectionTask_Embedding_embed #Module.Owner jmulliga #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetNCentBins(5); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); __R_ADDTASK__->GetClusterContainer(0)->SetClusNonLinCorrEnergyCut(0.); __R_ADDTASK__->GetClusterContainer(0)->SetClusHadCorrEnergyCut(0.3); __R_ADDTASK__->GetClusterContainer(0)->SetDefaultClusterEnergy(AliVCluster::kHadCorr); // Artifically reduce the embedded tracking efficiency __R_ADDTASK__->SetTrackEfficiency(0.98); __R_ADDTASK__->SetTrackEfficiencyOnlyForEmbedding(kTRUE); // Add track container for embedded event AliTrackContainer* tracksDetLevel = new AliTrackContainer("tracks"); tracksDetLevel->SetIsEmbedding(kTRUE); tracksDetLevel->SetParticlePtCut(0.15); __R_ADDTASK__->AdoptTrackContainer(tracksDetLevel); // EMCal clusters are already added to the jet finder in the combined cluster container #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetPerformance_Embedding #Module.Libs CGAL,fastjet,siscone,siscone_spherical,fastjetplugins,fastjettools,fastjetcontribfragile #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetPerformance.C #Module.MacroArgs "usedefault", "caloClustersCombined", "", 0.15, 0.3, "JetPerformance" #Module.Deps JetResponseMaker_Embedding_full_02 #Module.Owner jmulliga #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kINT7); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetNCentBins(5); // Set plotting options __R_ADDTASK__->SetPlotJetHistograms(kTRUE); __R_ADDTASK__->SetPlotClusterHistograms(kFALSE); __R_ADDTASK__->SetComputeBackground(kTRUE); __R_ADDTASK__->SetDoTriggerSimulation(kFALSE); __R_ADDTASK__->SetComputeMBDownscaling(kFALSE); __R_ADDTASK__->SetPlotMatchedJetHistograms(kTRUE); AliClusterContainer* clusCont = __R_ADDTASK__->GetClusterContainer("caloClustersCombined"); clusCont->SetClusHadCorrEnergyCut(0.3); clusCont->SetDefaultClusterEnergy(AliVCluster::kHadCorr); AliTrackContainer *trackContMC = __R_ADDTASK__->AddTrackContainer("tracks"); trackContMC->SetIsEmbedding(kTRUE); trackContMC->SetParticlePtCut(0.15); __R_ADDTASK__->SetUseAliEventCuts(kTRUE); __R_ADDTASK__->SetMakeGeneralHistograms(kFALSE); // Configure det-level jet container AliJetContainer* jetContFu02 = __R_ADDTASK__->AddJetContainer("Jet_AKTFullR020_tracks_pT0150_caloClustersCombined_E0300_pt_scheme", AliEmcalJet::kEMCALfid, 0.2); jetContFu02->SetRhoName("Rho_Scaled"); jetContFu02->SetPercAreaCut(0.6); jetContFu02->SetLeadingHadronType(0); jetContFu02->SetPtBiasJetTrack(5.); jetContFu02->SetPtBiasJetClus(5.); jetContFu02->SetJetPtCut(1.); jetContFu02->SetMaxTrackPt(100.); jetContFu02->SetIsParticleLevel(kFALSE); AliJetContainer* jetContGen = __R_ADDTASK__->AddJetContainer("Jet_AKTFullR020_mcparticles_pT0000_pt_scheme", AliEmcalJet::kEMCALfid, 0.2); jetContGen->SetPercAreaCut(0.6); jetContGen->SetLeadingHadronType(0); jetContGen->SetPtBiasJetTrack(5.); jetContGen->SetPtBiasJetClus(5.); jetContGen->SetJetPtCut(0.15); jetContGen->SetMaxTrackPt(1000.); jetContGen->SetIsParticleLevel(kTRUE); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetResponseMaker_Embedding_full_02 #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskJetResponseMaker.C #Module.MacroArgs "tracks", "caloClustersCombined", "Jet_AKTFullR020_tracks_pT0150_caloClustersCombined_E0300_pt_scheme", "Rho_Scaled", 0.2, "mcparticles", "", "Jet_AKTFullR020_mcparticles_pT0000_pt_scheme", "", 0.2, 1., 0., 5., 0, AliJetResponseMaker::kGeometrical, 0.25, 0.25, "EMCALfid",-999,0,10 #Module.Deps JetFinder_full_AKT_02_Embedding,MCJetFinder_full_AKT_02_Embedding,RhoTask_Embedding #Module.Owner jmulliga #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->GetJetContainer(0)->SetJetAcceptanceType(AliEmcalJet::kEMCALfid); __R_ADDTASK__->GetJetContainer(1)->SetJetAcceptanceType(AliEmcalJet::kEMCALfid); __R_ADDTASK__->GetJetContainer(1)->SetJetPtCut(1.); __R_ADDTASK__->GetJetContainer(1)->GetParticleContainer()->SetIsEmbedding(kTRUE); __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kINT7); __R_ADDTASK__->SetHistoType(0); __R_ADDTASK__->SetIsPythia(kFALSE); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetNCentBins(5); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin MCJetFinder_full_AKT_02_Embedding #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, 0., 0., kGhostArea, AliJetContainer::pt_scheme, "Jet", 0.1, kFALSE, kFALSE #Module.Deps EmbeddingHelper_Run2,MultSelectionTask #Module.Owner jmulliga #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetNCentBins(5); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); // Set mcparticles as embedded AliParticleContainer * truthTracks = __R_ADDTASK__->GetParticleContainer(0); truthTracks->SetIsEmbedding(kTRUE); #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 #_______________________________________________________________________________ #Module.Begin RhoTask_Embedding #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskRhoNew.C #Module.MacroArgs "usedefault", "usedefault", "Rho", 0.4, AliEmcalJet::kTPCfid, AliJetContainer::kChargedJet, kFALSE #Module.Deps JetFinder_charged_KT_04_Embedding #Module.Owner jmulliga #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SetExcludeLeadJets(2); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetNCentBins(5); __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->LoadRhoFunction("alien:///alice/cern.ch/user/j/jmulliga/scaleFactorEMCalLHC15o.root", "fScaleFactorEMCal"); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); #Module.EndConfig