#_______________________________________________________________________________ #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_Data_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 rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/r/rehlersi/embedding/embedDataUserConfiguration.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_Data_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 rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/r/rehlersi/embedding/embedDataUserConfiguration.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin ConfigureEmcalCorrectionTask_Embedding_Data_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 rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetUserConfigurationFilename("alien:///alice/cern.ch/user/r/rehlersi/embedding/embedDataUserConfiguration.yaml"); __R_ADDTASK__->Initialize(true); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmbeddingHelper_data #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWG/EMCAL/macros/AddTaskEmcalEmbeddingHelper.C #Module.MacroArgs #Module.Deps CDBConnectTask #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SetConfigurationPath("alien:///alice/cern.ch/user/r/rehlersi/embedding/embeddingConfigurationLHC11a.yaml"); __R_ADDTASK__->SetFilePattern("alien:///alice/data/2011/LHC11a/"); __R_ADDTASK__->SetInputFilename("*/ESDs/pass4_with_SDD/AOD113/*/AliAOD.root"); // We lose lots of statistics by selecting events here. Better to select at the task level in exchange for longer processing time //__R_ADDTASK__->SetTriggerMask(AliVEvent::kEMC1 | AliVEvent::kAnyINT); __R_ADDTASK__->SetTriggerMask(AliVEvent::kAny); __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetRandomFileAccess(kTRUE); __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_data #Module.Owner rehlersi #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_data #Module.Owner rehlersi #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_data #Module.Owner rehlersi #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_hybridJets_Embedding_Data #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetSample.C #Module.MacroArgs "usedefault", "caloClustersCombined", "emcalCellsCombined", "hybrid" #Module.Deps JetFinder_full_AKT_02_3GeVConstituent_hybridJets_Embedding_Data #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetHistoBins(600, 0, 300); // Particle containers __R_ADDTASK__->GetParticleContainer(0)->SetParticlePtCut(0.15); // Embedded tracks AliTrackContainer * tracksDetLevel = new AliTrackContainer("tracks"); tracksDetLevel->SetIsEmbedding(kTRUE); tracksDetLevel->SetParticlePtCut(0.15); // Track cuts are different for LHC11a and LHC11h tracksDetLevel->SetTrackCutsPeriod("lhc11a"); __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("hybridJets_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, "hybridJets"); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin EmcalJetSampleTask_truthJets_Embedding_Data #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetSample.C #Module.MacroArgs "usedefault", "usedefault", "", "truth" #Module.Deps JetFinder_full_AKT_02_3GeVConstituent_truthJets_Embedding_Data #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(kComPhysSel); __R_ADDTASK__->SetHistoBins(600, 0, 300); // Particle containers AliTrackContainer * tracks = __R_ADDTASK__->GetTrackContainer(0); tracks->SetIsEmbedding(kTRUE); tracks->SetParticlePtCut(0.15); // Track cuts are different for LHC11a and LHC11h tracks->SetTrackCutsPeriod("lhc11a"); // Cluster containers AliClusterContainer * clusterCont = __R_ADDTASK__->GetClusterContainer(0); clusterCont->SetIsEmbedding(kTRUE); // Jets __R_ADDTASK__->AddJetContainer("truthJets_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, "hybridJets"); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetFinder_full_AKT_02_3GeVConstituent_hybridJets_Embedding_Data #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, "hybridJets", 0., kFALSE, kFALSE #Module.Deps EmcalCorrectionTask_Embedding_combined,EmcalCorrectionTask_Embedding_data,EmcalCorrectionTask_Embedding_embed #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.); // Track cuts are different between LHC11a and LHC11h tracksDetLevel->SetTrackCutsPeriod("lhc11a"); __R_ADDTASK__->AdoptTrackContainer(tracksDetLevel); // Cluster containers __R_ADDTASK__->GetClusterContainer(0)->SetDefaultClusterEnergy(AliVCluster::kHadCorr); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetFinder_full_AKT_02_3GeVConstituent_truthJets_Embedding_Data #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, "truthJets", 0., kFALSE, kFALSE #Module.Deps EmcalCorrectionTask_Embedding_combined,EmcalCorrectionTask_Embedding_data,EmcalCorrectionTask_Embedding_embed #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 AliParticleContainer * tracks = __R_ADDTASK__->GetParticleContainer(0); tracks->SetIsEmbedding(kTRUE); AliTrackContainer * trackCont = dynamic_cast(tracks); if (trackCont) { trackCont->SetTrackCutsPeriod("lhc11a"); } AliClusterContainer * clusters = __R_ADDTASK__->GetClusterContainer(0); clusters->SetIsEmbedding(kTRUE); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetHadron_10GeVBias_embed_pp #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetHCorrelations.C #Module.MacroArgs "usedefault", "usedefault", AliAnalysisTaskEmcalJetHCorrelations::kDisableBias,10.0,50000,5000,1,10, AliVEvent::kEMC1 | AliVEvent::kAnyINT , AliVEvent::kEMC1 | AliVEvent::kAnyINT , kFALSE,kFALSE,0,kFALSE, "", "", "clusbias10R2EmbedNoTrackEff" #Module.Deps JetResponseMaker_Embedding_Data #Module.Owner rehlersi #Module.OutputFile AnalysisResults.root #Module.StartConfig __R_ADDTASK__->SelectCollisionCandidates(AliVEvent::kEMCEGA | AliVEvent::kMB | AliVEvent::kCentral | AliVEvent::kSemiCentral); __R_ADDTASK__->ConfigureForEmbeddingAnalysis(); #Module.EndConfig #_______________________________________________________________________________ #Module.Begin JetResponseMaker_Embedding_Data #Module.Libs #Module.DataTypes ESD, AOD, MC #Module.MacroName $ALICE_PHYSICS/PWGJE/EMCALJetTasks/macros/AddTaskJetResponseMaker.C #Module.MacroArgs "tracks", "caloClustersCombined", "hybridJets_AKTFullR020_tracks_pT3000_caloClustersCombined_E3000_pt_scheme", "", 0.2, "tracks", "caloClusters", "truthJets_AKTFullR020_tracks_pT3000_caloClusters_E3000_pt_scheme", "", 0.2, 0., 0., 0., 1, AliJetResponseMaker::kGeometrical, 1.2, 1.2, "EMCALfid", -999, -999, -999 #Module.Deps JetFinder_full_AKT_02_3GeVConstituent_hybridJets_Embedding_Data,JetFinder_full_AKT_02_3GeVConstituent_truthJets_Embedding_Data #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 the EMCal containers (although it will not use it for anything) __R_ADDTASK__->GetJetContainer(1)->GetParticleContainer()->SetIsEmbedding(kTRUE); __R_ADDTASK__->GetJetContainer(1)->GetClusterContainer()->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); #Module.EndConfig