#_______________________________________________________________________________ #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_Leticia #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/ConfigureEmcalEmbeddingHelperOnLEGOTrain.C #Module.MacroArgs #Module.Deps EmbeddingHelper_Leticia #Module.Owner lcunquei #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kCentral); // 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/2019/LHC19f4_2/%d/"); //__R_ADDTASK__->SetInputFilename("*/AOD/*/aod_archive.zip"); __R_ADDTASK__->SetInputFilename("*/AOD/*/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/l/lhavener/embeddingHelper_LHC18_LHC19f4_kCentral.yaml"); // Set the pt hard auto config identifier __R_ADDTASK__->SetAutoConfigureIdentifier("20190810Leticia"); // 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 EmbeddingHelper_Leticia #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/AddTaskEmcalEmbeddingHelper.C #Module.MacroArgs #Module.Deps CDBConnectTask,MultSelectionTask_Leticia #Module.Owner lcunquei #Module.OutputFile EventStat_temp.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 JetFinderAKTCharged_R04_EschemeDetLevelNew #Module.Libs CGAL,fastjet,siscone,siscone_spherical,fastjetplugins,fastjettools,fastjetcontribfragile,JETAN,FASTJETAN,ESDfilter,PWGGAEMCALTasks,PWGJEEMCALJetTasks #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "usedefault","",AliJetContainer::antikt_algorithm,0.4,AliJetContainer::kChargedJet,0.15,0.3,kGhostArea,AliJetContainer::E_scheme,"detLevelJets",0,kFALSE,kFALSE #Module.Deps CDBConnectTask,MultSelectionTask_Leticia,EmbeddingHelper_Leticia #Module.Owner lcunquei #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kCentral); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kpp); // Setup the tracks properly to be retrieved from the external event AliParticleContainer * tracks = __R_ADDTASK__->GetParticleContainer(0); tracks->SetIsEmbedding(kTRUE); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetFinderAKTCharged_R04_EschemeNew #Module.Libs CGAL,fastjet,siscone,siscone_spherical,fastjetplugins,fastjettools,fastjetcontribfragile,JETAN,FASTJETAN,ESDfilter,PWGGAEMCALTasks,PWGJEEMCALJetTasks #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "tracks","",AliJetContainer::antikt_algorithm,0.4,AliJetContainer::kChargedJet,0.15,0.3,kGhostArea,AliJetContainer::E_scheme,"hybridLevelJets",0,kFALSE,kFALSE #Module.Deps EmbeddingHelper_Leticia,JetFinderKtCharged_R02_EschemeNew,MultSelectionTask_Leticia,RhoESchemeNewSubstructure,RhoMassEschemeNewSubstructure #Module.Owner lcunquei #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kCentral); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); // 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); AliEmcalJetUtilityConstSubtractor* constUtil = (AliEmcalJetUtilityConstSubtractor*)__R_ADDTASK__->AddUtility(new AliEmcalJetUtilityConstSubtractor("ConstSubtractor")); AliEmcalJetUtilityGenSubtractor* genSub =(AliEmcalJetUtilityGenSubtractor*) __R_ADDTASK__->AddUtility(new AliEmcalJetUtilityGenSubtractor("GenSubtractor")); //genSub->SetGenericSubtractionJetMass(kTRUE); //genSub->SetGenericSubtractionExtraJetShapes(kTRUE); genSub->SetUseExternalBkg(kTRUE); genSub->SetRhoName("Rho"); genSub->SetRhomName("Rhomass"); constUtil->SetJetsSubName(Form("%sConstSub",__R_ADDTASK__->GetName())); constUtil->SetParticlesSubName("tracksSubR04"); constUtil->SetUseExternalBkg(kTRUE); constUtil->SetRhoName("Rho"); constUtil->SetRhomName("Rhomass"); // Particle containers AliTrackContainer * tracksDetLevel = new AliTrackContainer("tracks"); // Get the det level tracks from the external event! tracksDetLevel->SetIsEmbedding(kTRUE); __R_ADDTASK__->AdoptTrackContainer(tracksDetLevel); // Recycle unused embedded events __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetFinderAKTCharged_R04_EschemePartLevelNew #Module.Libs CGAL,fastjet,siscone,siscone_spherical,fastjetplugins,fastjettools,fastjetcontribfragile,JETAN,FASTJETAN,ESDfilter,PWGGAEMCALTasks,PWGJEEMCALJetTasks #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "mcparticles","",AliJetContainer::antikt_algorithm,0.4,AliJetContainer::kChargedJet,0,kClusPtCut,kGhostArea,AliJetContainer::E_scheme,"partLevelJets",0,kFALSE,kFALSE #Module.Deps CDBConnectTask,EmbeddingHelper_Leticia,MultSelectionTask_Leticia #Module.Owner lcunquei #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kCentral); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kpp); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); // 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 JetFinderKtCharged_R02_EschemeNew #Module.Libs CGAL,fastjet,siscone,siscone_spherical,fastjetplugins,fastjettools,fastjetcontribfragile,JETAN,FASTJETAN,ESDfilter,PWGGAEMCALTasks,PWGJEEMCALJetTasks #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C #Module.MacroArgs "usedefault","",AliJetContainer::kt_algorithm,0.2,AliJetContainer::kChargedJet,0.15,0,kGhostArea,AliJetContainer::E_scheme,"Jet",0,kFALSE,kFALSE #Module.Deps CDBConnectTask,EmbeddingHelper_Leticia,MultSelectionTask_Leticia #Module.Owner lcunquei #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kCentral); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetTaggerR04DetHybridNew #Module.Libs CGAL,fastjet,siscone,siscone_spherical,fastjetplugins,fastjettools,fastjetcontribfragile,PWGJEEMCALJetTasks #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetTaggerFast.C #Module.MacroArgs "hybridLevelJets_AKTChargedR040_tracks_pT0150_E_scheme","detLevelJets_AKTChargedR040_tracks_pT0150_E_scheme",0.4,"Rho","","tracks","","TPCfid","V0M",kComPhysSel,"" #Module.Deps JetFinderAKTCharged_R04_EschemeDetLevelNew,JetFinderAKTCharged_R04_EschemePartLevelNew,JetFinderAKTCharged_R04_EschemeNew #Module.Owner lcunquei #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(AliVEvent::kCentral); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); // 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.EndConfig #_______________________________________________________________________________ #Module.Begin JetTaggerR04DetPartNew #Module.Libs CGAL,fastjet,siscone,siscone_spherical,fastjetplugins,fastjettools,fastjetcontribfragile,PWGJEEMCALJetTasks #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetTaggerFast.C #Module.MacroArgs "detLevelJets_AKTChargedR040_tracks_pT0150_E_scheme","partLevelJets_AKTChargedR040_mcparticles_pT0000_E_scheme",0.4,"","","tracks","","TPCfid","V0M",kComPhysSel,"" #Module.Deps JetFinderAKTCharged_R04_EschemeDetLevelNew,JetFinderAKTCharged_R04_EschemePartLevelNew,JetTaggerR04DetHybridNew #Module.Owner lcunquei #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(AliVEvent::kCentral); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); // Reapply the max track pt cut off to maintain energy resolution and avoid fake tracks auto detLevelJetCont = __R_ADDTASK__->GetJetContainer(0); detLevelJetCont->SetMaxTrackPt(100); auto partLevelJetCont = __R_ADDTASK__->GetJetContainer(1); partLevelJetCont->SetMaxTrackPt(1000); #Module.EndConfig #Module.EndConfig #_______________________________________________________________________________ #Module.Begin MultSelectionTask_Leticia #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 lcunquei #Module.OutputFile EventStat_temp.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kCentral); /* 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 RhoESchemeNewSubstructure #Module.Libs CGAL,fastjet,siscone,siscone_spherical,fastjetplugins,fastjettools,fastjetcontribfragile,PWGJEEMCALJetTasks #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskRhoNew.C #Module.MacroArgs "usedefault", "","Rho", 0.2,AliJetContainer::kTPCfid,AliJetContainer::kChargedJet,kTRUE,AliJetContainer::E_scheme #Module.Deps CDBConnectTask,EmbeddingHelper_Leticia,JetFinderKtCharged_R02_EschemeNew,MultSelectionTask_Leticia #Module.Owner lcunquei #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetExcludeLeadJets(2); __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kCentral); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetNCentBins(5); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); __R_ADDTASK__->SetZvertexDiffValue(0.1); AliJetContainer *cont = __R_ADDTASK__->GetJetContainer(0); cont->SetJetRadius(0.2); cont->SetJetAcceptanceType(AliJetContainer::kTPCfid); cont->SetMaxTrackPt(100); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin RhoMassEschemeNewSubstructure #Module.Libs CGAL,fastjet,siscone,siscone_spherical,fastjetplugins,fastjettools,fastjetcontribfragile,PWGJEEMCALJetTasks #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskRhoMass.C #Module.MacroArgs "Jet_KTChargedR020_tracks_pT0150_E_scheme","tracks","","Rhomass",0.2,"TPC",0.01,0,0,2,kTRUE,"RhoMass" #Module.Deps CDBConnectTask,EmbeddingHelper_Leticia,JetFinderKtCharged_R02_EschemeNew,MultSelectionTask_Leticia #Module.Owner lcunquei #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kCentral); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetHistoBins(250,0,250); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); AliJetContainer *cont = __R_ADDTASK__->GetJetContainer(0); cont->SetJetRadius(0.2); cont->SetJetAcceptanceType(AliJetContainer::kTPCfid); cont->SetMaxTrackPt(100); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); __R_ADDTASK__->SetZvertexDiffValue(0.1); __R_ADDTASK__->SetVzRange(-10,10); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin SoftDropResponse_ChargedJets_R04_1 #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalSoftdropResponse.C #Module.MacroArgs 0.4, AliJetContainer::kChargedJet, AliJetContainer::E_scheme, kTRUE, "mcparticles", "1" #Module.Deps CDBConnectTask,JetTaggerR04DetHybridNew,JetTaggerR04DetPartNew #Module.Owner lhavener #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetMakeGeneralHistograms(kTRUE); __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kCentral); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetForceBeamType(AliAnalysisTaskEmcal::kAA); __R_ADDTASK__->SetRecycleUnusedEmbeddedEventsMode(kTRUE); AliJetContainer *contConstSub =__R_ADDTASK__->GetJetContainer(1); contConstSub->SetRhoName("Rho"); contConstSub->SetRhoMassName("RhoMass"); AliJetContainer* jetContCh04 = __R_ADDTASK__->AddJetContainer("hybridLevelJets_AKTChargedR040_tracks_pT0150_E_schemeConstSub"); __R_ADDTASK__->SetNameDetLevelJetContainer(jetContCh04->GetName()); jetContCh04->SetJetPtCut(0.); jetContCh04->SetMaxTrackPt(1000.); AliParticleContainer * partCont = __R_ADDTASK__->GetMCParticleContainer(0); partCont->SetIsEmbedding(kTRUE); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin SubstructureEmbedding_ConstSub #Module.Libs CGAL,fastjet,siscone,siscone_spherical,fastjetplugins,fastjettools,fastjetcontribfragile,PWGJEEMCALJetTasks #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskNewJetSubstructure.C #Module.MacroArgs "hybridLevelJets_AKTChargedR040_tracks_pT0150_E_schemeConstSub","hybridLevelJets_AKTChargedR040_tracks_pT0150_E_scheme","detLevelJets_AKTChargedR040_tracks_pT0150_E_scheme","partLevelJets_AKTChargedR040_mcparticles_pT0000_E_scheme",0.4,"Rho","tracksSubR04","tracks","mcparticles","","mcparticles","TPC","V0M",1<<31,"","","Raw",AliAnalysisTaskNewJetSubstructure::kDetEmbPartPythia, AliAnalysisTaskNewJetSubstructure::kConstSub,AliAnalysisTaskNewJetSubstructure::kInclusive #Module.Deps RhoESchemeNewSubstructure,RhoMassEschemeNewSubstructure,JetTaggerR04DetPartNew,JetTaggerR04DetHybridNew #Module.Owner lcunquei #Module.OutputFile AnalysisResults.root #Module.StartConfig AliJetContainer *cont = __R_ADDTASK__->GetJetContainer(0); cont->SetRhoName("Rho"); cont->SetRhoMassName("RhoMass"); cont->SetJetRadius(0.4); cont->SetJetAcceptanceType(AliJetContainer::kTPCfid); cont->SetMaxTrackPt(100); auto detLevelJetCont = __R_ADDTASK__->GetJetContainer(2); detLevelJetCont->SetMaxTrackPt(100); auto partLevelJetCont = __R_ADDTASK__->GetJetContainer(3); partLevelJetCont->SetMaxTrackPt(1000); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); __R_ADDTASK__->SetJetPtThreshold(20); __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kCentral); __R_ADDTASK__->SetNeedEmcalGeom(kFALSE); //__R_ADDTASK__->SetZvertexDiffValue(0.1); //__R_ADDTASK__->SetVzRange(-10,10); //__R_ADDTASK__->SetDoTwoTrack(kTRUE); __R_ADDTASK__->SetHardCutoff(0.1); __R_ADDTASK__->SetMinCentrality(0); __R_ADDTASK__->SetMaxCentrality(10); __R_ADDTASK__->SetPowerAlgorithm(0); __R_ADDTASK__->SetMinFractionShared(0.5); __R_ADDTASK__->SetUseNewCentralityEstimation(kTRUE); #Module.EndConfig