How to access the corrected chunks after d0correct ------------------------------------------------------------------ EM chunk: --------------- see: analysis_example/src/D0ChunkAnalyze.cpp and analysis_example/rcp/D0ChunkAnalyze.rcp for an example. RCP emidAlgoRCP = rcp.getRCP("EMid_Algo"); // RCP for selecting the EMparticle chunk used vector emidNestedRCP = rcp.getVString("EMid_SearchRCPs"); EMparticleChunkSelector emAlgoSel(emidAlgoRCP,emidNestedRCP); edm::TKey emKey = TKey(emAlgoSel); THandle emChunk = emKey.find(event); with, in your RCP: RCP EMid_Algo = < emreco EMReco-scone-id > string EMid_SearchRCPs = ("clusterer","HMReco",) Muon chunk: ----------------- see: MuoCandidate documentation for details (http://www-d0.fnal.gov/phys_id/muon_id/d0_private/certif/p14/index.html) RCP muonidAlgoRCP = rcp.getRCP("Muonid_Algo"); MuonParticleChunkSelector muAlgoSel(muonidAlgoRCP); TKey candKey = TKey(muAlgoSel); THandle candchunk= candKey.find(event); with, in your RCP: RCP Muonid_Algo = Jet chunk: -------------- create a JetAlgoInfo for either JCCA or JCCB and use the method: JetAlgoInfo JetAlgoInfo::CorrectedJetAlgoInfo(JetAlgoInfo uncorrectedJetAlgoInfo); to access the corresponding corrected jet chunk with good jets corrected for JES without muons and JetAlgoInfo JetAlgoInfo::CorrectedMuJetAlgoInfo(JetAlgoInfo uncorrectedJetAlgoInfo); to access the corresponding corrected jet chunk with good jets corrected for JES with muons string jetType = rcp.getString("JetAlgo_type"); std::vector jetValues = rcp.getVFloat("JetAlgo_values"); std::vector jetNames = rcp.getVString("JetAlgo_names"); if ( jetNames.size() > 0 ) while ( *jetNames.rbegin() == "" ) jetNames.pop_back(); if ( jetValues.size() < jetNames.size() ) { cout << "D0ChunkAnalyze: inconsistent RCP values for the jet selection, disable it" << endl; jetNames.clear(); jetValues.clear(); } JetAlgoInfo uncorrected_jetAlgoInfo(jetType,jetValues,jetNames); JetAlgoInfo corrected_jetAlgoInfo=JetAlgoInfo::CorrectedJetAlgoInfo(uncorrected_jetAlgoInfo); JetChunkSelector jsel=JetChunkSelector(corrected_jetAlgoInfo); const TKey jetKey(jsel); THandle hjet = jetKey.find(event); with for JCCA, in your RCP: string JetAlgo_type = "PreSCilcone" string JetAlgo_names = ( "towers" "coneSize" "Radius_of_Cone" "Min_Jet_ET") float JetAlgo_values = ( 0. 0.3 0.7 8. ) Missing ET chunk: ------------------------ you can access the corrected chunk as the uncorrected one: const TKey metKey; THandle metchunk = metKey.find(event);