.. _program_listing_file_Root_JetContainer.cxx: Program Listing for File JetContainer.cxx ========================================= |exhale_lsh| :ref:`Return to documentation for file ` (``Root/JetContainer.cxx``) .. |exhale_lsh| unicode:: U+021B0 .. UPWARDS ARROW WITH TIP LEFTWARDS .. code-block:: cpp #include "xAODAnaHelpers/JetContainer.h" #include #include #include // std::domain_error #include "xAODTruth/TruthEventContainer.h" #include "xAODBTagging/BTaggingUtilities.h" using namespace xAH; JetContainer::JetContainer(const std::string& name, const std::string& detailStr, float units, bool mc) : ParticleContainer(name,detailStr,units,mc), m_trkSelTool(nullptr) { // rapidity if(m_infoSwitch.m_rapidity) { m_rapidity =new std::vector(); } // trigger if ( m_infoSwitch.m_trigger ) { m_isTrigMatched = new std::vector (); m_isTrigMatchedToChain = new std::vector > (); m_listTrigChains = new std::vector (); } // clean if(m_infoSwitch.m_cleanTrig && ! (m_infoSwitch.m_clean || m_infoSwitch.m_cleanLight || m_infoSwitch.m_cleanLLP)) { std::cout << "JetContainer WARNING You asked for cleanTrig for " << name << "but didn't specify clean, cleanLight or cleanLLP. Going to assume you wanted clean." << std::endl; m_infoSwitch.m_clean = true; } if(m_infoSwitch.m_clean || m_infoSwitch.m_cleanLight || m_infoSwitch.m_cleanLLP) { if(m_infoSwitch.m_clean){ m_Timing =new std::vector(); m_LArQuality =new std::vector(); m_HECQuality =new std::vector(); m_NegativeE =new std::vector(); m_AverageLArQF =new std::vector(); m_BchCorrCell =new std::vector(); m_N90Constituents =new std::vector(); m_LArBadHVEnergyFrac =new std::vector(); m_LArBadHVNCell =new std::vector(); m_ChargedFraction =new std::vector(); m_OotFracClusters5 =new std::vector(); m_OotFracClusters10 =new std::vector(); m_LeadingClusterPt =new std::vector(); m_LeadingClusterSecondLambda=new std::vector(); m_LeadingClusterCenterLambda=new std::vector(); m_LeadingClusterSecondR =new std::vector(); if(m_infoSwitch.m_cleanTrig) { m_clean_passLooseBadTriggerUgly=new std::vector(); } else { m_clean_passLooseBadUgly =new std::vector(); m_clean_passTightBadUgly =new std::vector(); } } if(m_infoSwitch.m_cleanTrig) { m_clean_passLooseBadTrigger =new std::vector(); } else { m_clean_passLooseBad =new std::vector(); m_clean_passTightBad =new std::vector(); } if(m_infoSwitch.m_cleanLLP) { m_clean_passLooseBadLLP =new std::vector(); } } if(m_infoSwitch.m_timing && !m_infoSwitch.m_clean){ m_Timing = new std::vector(); } // energy if ( m_infoSwitch.m_energy || m_infoSwitch.m_energyLight ) { if ( m_infoSwitch.m_energy ) { m_HECFrac = new std::vector(); m_CentroidR = new std::vector(); m_LowEtConstituentsFrac = new std::vector(); } m_EMFrac = new std::vector(); m_FracSamplingMax = new std::vector(); m_FracSamplingMaxIndex = new std::vector(); m_GhostMuonSegmentCount = new std::vector(); m_Width = new std::vector(); } // scales if ( m_infoSwitch.m_scales ) { m_emScalePt = new std::vector(); m_constScalePt = new std::vector(); m_pileupScalePt = new std::vector(); m_originConstitScalePt = new std::vector(); m_etaJESScalePt = new std::vector(); m_gscScalePt = new std::vector(); m_jmsScalePt = new std::vector(); m_insituScalePt = new std::vector(); m_emScaleM = new std::vector(); m_constScaleM = new std::vector(); m_pileupScaleM = new std::vector(); m_originConstitScaleM = new std::vector(); m_etaJESScaleM = new std::vector(); m_gscScaleM = new std::vector(); m_jmsScaleM = new std::vector(); m_insituScaleM = new std::vector(); } // constscale eta if ( m_infoSwitch.m_constscaleEta ) { m_constScaleEta = new std::vector(); } // detector eta if ( m_infoSwitch.m_detectorEta ) { m_detectorEta = new std::vector(); } // layer if ( m_infoSwitch.m_layer ) { m_EnergyPerSampling = new std::vector< std::vector >(); } // tracksAll if ( m_infoSwitch.m_trackAll ) { m_NumTrkPt1000 = new std::vector< std::vector > (); m_SumPtTrkPt1000 = new std::vector< std::vector > (); m_TrackWidthPt1000 = new std::vector< std::vector > (); m_NumTrkPt500 = new std::vector< std::vector > (); m_SumPtTrkPt500 = new std::vector< std::vector > (); m_TrackWidthPt500 = new std::vector< std::vector > (); m_JVF = new std::vector< std::vector > (); } // trackPV if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_jvt ) { if (m_infoSwitch.m_trackPV){ m_NumTrkPt1000PV = new std::vector(); m_SumPtTrkPt1000PV = new std::vector(); m_TrackWidthPt1000PV = new std::vector(); m_NumTrkPt500PV = new std::vector(); m_SumPtTrkPt500PV = new std::vector(); m_TrackWidthPt500PV = new std::vector(); m_JVFPV = new std::vector(); m_JvtJvfcorr = new std::vector(); m_JvtRpt = new std::vector(); } m_Jvt = new std::vector(); } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sfJVTName == "Loose" ) { m_JvtPass_Loose = new std::vector(); if ( m_mc ) { m_JvtEff_SF_Loose = new std::vector< std::vector > (); } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sfJVTName == "Medium" ) { m_JvtPass_Medium = new std::vector(); if ( m_mc ) { m_JvtEff_SF_Medium = new std::vector< std::vector > (); } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sfJVTName == "Tight" ) { m_JvtPass_Tight = new std::vector(); if ( m_mc ) { m_JvtEff_SF_Tight = new std::vector< std::vector > (); } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sfJVTName == "FixedEffPt" ) { m_JvtPass_FixedEffPt = new std::vector(); if ( m_mc ) { m_JvtEff_SF_FixedEffPt = new std::vector< std::vector > (); } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sfJVTName == "TightFwd" ) { m_JvtPass_TightFwd = new std::vector(); if ( m_mc ) { m_JvtEff_SF_TightFwd = new std::vector< std::vector > (); } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sffJVTName == "Loose" ) { m_fJvtPass_Loose = new std::vector(); if ( m_mc ) { m_fJvtEff_SF_Loose = new std::vector< std::vector > (); } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sffJVTName == "Medium" ) { m_fJvtPass_Medium = new std::vector(); if ( m_mc ) { m_fJvtEff_SF_Medium = new std::vector< std::vector > (); } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sffJVTName == "Tight" ) { m_fJvtPass_Tight = new std::vector(); if ( m_mc ) { m_fJvtEff_SF_Tight = new std::vector< std::vector > (); } } if (m_infoSwitch.m_fJvt) { m_fJvt = new std::vector(); } if (m_infoSwitch.m_NNJvt) { m_NNJvt = new std::vector(); m_NNJvtPass = new std::vector(); } // chargedPFOPV if ( m_infoSwitch.m_chargedPFOPV ) { m_SumPtChargedPFOPt500PV = new std::vector (); m_fCharged = new std::vector (); } // allTrack // trackAll or trackPV if ( m_infoSwitch.m_allTrack ) { m_GhostTrackCount = new std::vector (); m_GhostTrackPt = new std::vector (); m_GhostTrack_pt = new std::vector< std::vector > (); m_GhostTrack_qOverP = new std::vector< std::vector > (); m_GhostTrack_eta = new std::vector< std::vector > (); m_GhostTrack_phi = new std::vector< std::vector > (); m_GhostTrack_e = new std::vector< std::vector > (); m_GhostTrack_d0 = new std::vector< std::vector > (); m_GhostTrack_z0 = new std::vector< std::vector > (); // allTrackDetail if(m_infoSwitch.m_allTrackDetail){ m_GhostTrack_nPixelHits = new std::vector< std::vector >(); m_GhostTrack_nSCTHits = new std::vector< std::vector >(); m_GhostTrack_nTRTHits = new std::vector< std::vector >(); m_GhostTrack_nPixelSharedHits = new std::vector< std::vector >(); m_GhostTrack_nPixelSplitHits = new std::vector< std::vector >(); m_GhostTrack_nInnermostPixelLayerHits = new std::vector< std::vector >(); m_GhostTrack_nInnermostPixelLayerSharedHits = new std::vector< std::vector >(); m_GhostTrack_nInnermostPixelLayerSplitHits = new std::vector< std::vector >(); m_GhostTrack_nNextToInnermostPixelLayerHits = new std::vector< std::vector >(); m_GhostTrack_nNextToInnermostPixelLayerSharedHits = new std::vector< std::vector >(); m_GhostTrack_nNextToInnermostPixelLayerSplitHits = new std::vector< std::vector >(); } } // constituent if ( m_infoSwitch.m_constituent ) { m_numConstituents = new std::vector< int > (); } if ( m_infoSwitch.m_constituentAll ) { m_constituentWeights = new std::vector< std::vector > (); m_constituent_pt = new std::vector< std::vector > (); m_constituent_eta = new std::vector< std::vector > (); m_constituent_phi = new std::vector< std::vector > (); m_constituent_e = new std::vector< std::vector > (); } // flavorTag if( m_infoSwitch.m_flavorTag || m_infoSwitch.m_flavorTagHLT || m_infoSwitch.m_flavorTagTLA ) { m_HadronConeExclTruthLabelID =new std::vector(); m_HadronConeExclExtendedTruthLabelID=new std::vector(); } if( m_infoSwitch.m_flavorTag || m_infoSwitch.m_flavorTagHLT ) { m_DL1r =new std::vector(); m_DL1r_pu =new std::vector(); m_DL1r_pc =new std::vector(); m_DL1r_pb =new std::vector(); m_DL1dv00 =new std::vector(); m_DL1dv00_pu =new std::vector(); m_DL1dv00_pc =new std::vector(); m_DL1dv00_pb =new std::vector(); m_DL1dv01 =new std::vector(); m_DL1dv01_pu =new std::vector(); m_DL1dv01_pc =new std::vector(); m_DL1dv01_pb =new std::vector(); m_GN1 =new std::vector(); m_GN1_pu =new std::vector(); m_GN1_pc =new std::vector(); m_GN1_pb =new std::vector(); m_GN2v00LegacyWP =new std::vector(); m_GN2v00LegacyWP_pu =new std::vector(); m_GN2v00LegacyWP_pc =new std::vector(); m_GN2v00LegacyWP_pb =new std::vector(); m_GN2v00NewAliasWP =new std::vector(); m_GN2v00NewAliasWP_pu =new std::vector(); m_GN2v00NewAliasWP_pc =new std::vector(); m_GN2v00NewAliasWP_pb =new std::vector(); m_GN2v01 =new std::vector(); m_GN2v01_pu =new std::vector(); m_GN2v01_pc =new std::vector(); m_GN2v01_pb =new std::vector(); m_GN2v01_ptau =new std::vector(); // Jet Fitter if( m_infoSwitch.m_jetFitterDetails){ m_JetFitter_nVTX = new std::vector(); m_JetFitter_nSingleTracks = new std::vector(); m_JetFitter_nTracksAtVtx = new std::vector(); m_JetFitter_mass = new std::vector(); m_JetFitter_energyFraction = new std::vector(); m_JetFitter_significance3d = new std::vector(); m_JetFitter_deltaeta = new std::vector(); m_JetFitter_deltaphi = new std::vector(); m_JetFitter_N2Tpar = new std::vector(); } // SV Details if( m_infoSwitch.m_svDetails){ m_SV0 = new std::vector(); m_sv0_NGTinSvx = new std::vector(); m_sv0_N2Tpair = new std::vector(); m_sv0_massvx = new std::vector(); m_sv0_efracsvx = new std::vector(); m_sv0_normdist = new std::vector(); m_SV1 = new std::vector(); m_SV1IP3D = new std::vector(); m_COMBx = new std::vector(); m_sv1_pu = new std::vector(); m_sv1_pb = new std::vector(); m_sv1_pc = new std::vector(); m_sv1_c = new std::vector(); m_sv1_cu = new std::vector(); m_sv1_NGTinSvx = new std::vector(); m_sv1_N2Tpair = new std::vector(); m_sv1_massvx = new std::vector(); m_sv1_efracsvx = new std::vector(); m_sv1_normdist = new std::vector(); m_sv1_Lxy = new std::vector(); m_sv1_sig3d = new std::vector(); m_sv1_L3d = new std::vector(); m_sv1_distmatlay = new std::vector(); m_sv1_dR = new std::vector(); } // IP3D if( m_infoSwitch.m_ipDetails){ m_IP2D_pu = new std::vector(); m_IP2D_pb = new std::vector(); m_IP2D_pc = new std::vector(); m_IP2D = new std::vector(); m_IP2D_c = new std::vector(); m_IP2D_cu = new std::vector(); m_nIP2DTracks = new std::vector(); m_IP2D_gradeOfTracks = new std::vector >(); m_IP2D_flagFromV0ofTracks = new std::vector >(); m_IP2D_valD0wrtPVofTracks = new std::vector >(); m_IP2D_sigD0wrtPVofTracks = new std::vector >(); m_IP2D_weightBofTracks = new std::vector >(); m_IP2D_weightCofTracks = new std::vector >(); m_IP2D_weightUofTracks = new std::vector >(); m_IP3D = new std::vector(); m_IP3D_pu = new std::vector(); m_IP3D_pb = new std::vector(); m_IP3D_pc = new std::vector(); m_IP3D_c = new std::vector(); m_IP3D_cu = new std::vector(); m_nIP3DTracks = new std::vector(); m_IP3D_gradeOfTracks = new std::vector >(); m_IP3D_flagFromV0ofTracks = new std::vector >(); m_IP3D_valD0wrtPVofTracks = new std::vector >(); m_IP3D_sigD0wrtPVofTracks = new std::vector >(); m_IP3D_valZ0wrtPVofTracks = new std::vector >(); m_IP3D_sigZ0wrtPVofTracks = new std::vector >(); m_IP3D_weightBofTracks = new std::vector >(); m_IP3D_weightCofTracks = new std::vector >(); m_IP3D_weightUofTracks = new std::vector >(); } if( m_infoSwitch.m_JVC ) { m_JetVertexCharge_discriminant = new std::vector(); } } // flavorTagHLT if( m_infoSwitch.m_flavorTagHLT ) { m_vtxOnlineValid = new std::vector(); m_vtxHadDummy = new std::vector(); m_bs_online_vx = new std::vector(); m_bs_online_vy = new std::vector(); m_bs_online_vz = new std::vector(); m_vtx_offline_x0 = new std::vector(); m_vtx_offline_y0 = new std::vector(); m_vtx_offline_z0 = new std::vector(); m_vtx_online_x0 = new std::vector(); m_vtx_online_y0 = new std::vector(); m_vtx_online_z0 = new std::vector(); m_vtx_online_bkg_x0 = new std::vector(); m_vtx_online_bkg_y0 = new std::vector(); m_vtx_online_bkg_z0 = new std::vector(); } if( m_infoSwitch.m_flavorTagTLA ) { m_fastDIPS =new std::vector(); m_fastDIPS_pu =new std::vector(); m_fastDIPS_pc =new std::vector(); m_fastDIPS_pb =new std::vector(); } if( !m_infoSwitch.m_jetBTag.empty() ) { std::stringstream ss_wpName; for(const auto& btaginfo : m_infoSwitch.m_jetBTag) { for(auto wp : btaginfo.second) { ss_wpName.str(""); ss_wpName << wp.first << "_" << std::setfill('0') << std::setw(2) << wp.second; m_btags.push_back(new btagOpPoint(m_mc, btaginfo.first, ss_wpName.str())); } } } if ( !m_infoSwitch.m_jetBTagCts.empty() ){ for(const auto& tagger: m_infoSwitch.m_jetBTagCts) { m_btags.push_back(new btagOpPoint(m_mc,tagger,"Continuous")); } } // area if( m_infoSwitch.m_area ) { m_GhostArea = new std::vector(); m_ActiveArea = new std::vector(); m_VoronoiArea = new std::vector(); m_ActiveArea4vec_pt = new std::vector(); m_ActiveArea4vec_eta = new std::vector(); m_ActiveArea4vec_phi = new std::vector(); m_ActiveArea4vec_m = new std::vector(); } // truth if ( m_infoSwitch.m_truth && m_mc ) { m_ConeTruthLabelID =new std::vector(); m_TruthCount =new std::vector(); m_TruthLabelDeltaR_B=new std::vector; m_TruthLabelDeltaR_C=new std::vector; m_TruthLabelDeltaR_T=new std::vector; m_PartonTruthLabelID=new std::vector(); m_GhostTruthAssociationFraction=new std::vector; m_truth_E =new std::vector; m_truth_pt =new std::vector; m_truth_phi=new std::vector; m_truth_eta=new std::vector; } // truth detail if ( m_infoSwitch.m_truthDetails ) { m_GhostBHadronsFinalCount = new std::vector (); m_GhostBHadronsInitialCount = new std::vector (); m_GhostBQuarksFinalCount = new std::vector (); m_GhostBHadronsFinalPt = new std::vector (); m_GhostBHadronsInitialPt = new std::vector (); m_GhostBQuarksFinalPt = new std::vector (); m_GhostCHadronsFinalCount = new std::vector (); m_GhostCHadronsInitialCount = new std::vector (); m_GhostCQuarksFinalCount = new std::vector (); m_GhostCHadronsFinalPt = new std::vector (); m_GhostCHadronsInitialPt = new std::vector (); m_GhostCQuarksFinalPt = new std::vector (); m_GhostTausFinalCount = new std::vector (); m_GhostTausFinalPt = new std::vector (); m_truth_pdgId = new std::vector (); m_truth_partonPt = new std::vector (); m_truth_partonDR = new std::vector (); } // charge if ( m_infoSwitch.m_charge ) { m_charge =new std::vector(); } // passSel if ( m_infoSwitch.m_passSel ) { m_passSel =new std::vector(); } // passOR if ( m_infoSwitch.m_passOR ) { m_passOR =new std::vector(); } } JetContainer::~JetContainer() { if(m_debug) std::cout << " Deleting JetContainer " << std::endl; if(m_infoSwitch.m_rapidity) { delete m_rapidity; } // trigger if ( m_infoSwitch.m_trigger ) { delete m_isTrigMatched ; delete m_isTrigMatchedToChain ; delete m_listTrigChains ; } // clean if(m_infoSwitch.m_clean || m_infoSwitch.m_cleanLight || m_infoSwitch.m_cleanLLP) { if(m_infoSwitch.m_clean){ delete m_Timing; delete m_LArQuality; delete m_HECQuality; delete m_NegativeE; delete m_AverageLArQF; delete m_BchCorrCell; delete m_N90Constituents; delete m_LArBadHVEnergyFrac; delete m_LArBadHVNCell; delete m_ChargedFraction; delete m_OotFracClusters5; delete m_OotFracClusters10; delete m_LeadingClusterPt; delete m_LeadingClusterSecondLambda; delete m_LeadingClusterCenterLambda; delete m_LeadingClusterSecondR; if(m_infoSwitch.m_cleanTrig) { delete m_clean_passLooseBadTriggerUgly; } else { delete m_clean_passLooseBadUgly; delete m_clean_passTightBadUgly; } } if(m_infoSwitch.m_cleanTrig) { delete m_clean_passLooseBadTrigger; } else { delete m_clean_passLooseBad; delete m_clean_passTightBad; } if(m_infoSwitch.m_cleanLLP) { delete m_clean_passLooseBadLLP; } } if(m_infoSwitch.m_timing && !m_infoSwitch.m_clean){ delete m_Timing; } // energy if ( m_infoSwitch.m_energy || m_infoSwitch.m_energyLight ) { if ( m_infoSwitch.m_energy ){ delete m_HECFrac; delete m_CentroidR; delete m_LowEtConstituentsFrac; } delete m_EMFrac; delete m_FracSamplingMax; delete m_FracSamplingMaxIndex; delete m_GhostMuonSegmentCount; delete m_Width; } // scales if ( m_infoSwitch.m_scales ) { delete m_emScalePt ; delete m_constScalePt ; delete m_pileupScalePt ; delete m_originConstitScalePt ; delete m_etaJESScalePt ; delete m_gscScalePt ; delete m_jmsScalePt ; delete m_insituScalePt ; delete m_emScaleM ; delete m_constScaleM ; delete m_pileupScaleM ; delete m_originConstitScaleM ; delete m_etaJESScaleM ; delete m_gscScaleM ; delete m_jmsScaleM ; delete m_insituScaleM ; } // constscale eta if ( m_infoSwitch.m_constscaleEta ) { delete m_constScaleEta ; } // detector eta if ( m_infoSwitch.m_detectorEta ) { delete m_detectorEta; } // layer if ( m_infoSwitch.m_layer ) { delete m_EnergyPerSampling; } // tracksAll if ( m_infoSwitch.m_trackAll ) { delete m_NumTrkPt1000 ; delete m_SumPtTrkPt1000 ; delete m_TrackWidthPt1000 ; delete m_NumTrkPt500 ; delete m_SumPtTrkPt500 ; delete m_TrackWidthPt500 ; delete m_JVF ; } // trackPV if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_jvt ) { if ( m_infoSwitch.m_trackPV ){ delete m_NumTrkPt1000PV; delete m_SumPtTrkPt1000PV; delete m_TrackWidthPt1000PV; delete m_NumTrkPt500PV; delete m_SumPtTrkPt500PV; delete m_TrackWidthPt500PV; delete m_JVFPV; delete m_JvtJvfcorr; delete m_JvtRpt; } delete m_Jvt; } if (m_infoSwitch.m_fJvt){ delete m_fJvt; } if (m_infoSwitch.m_NNJvt){ delete m_NNJvt; delete m_NNJvtPass; } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sfJVTName == "Loose" ) { delete m_JvtPass_Loose; if ( m_mc ) { delete m_JvtEff_SF_Loose; } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sfJVTName == "Medium" ) { delete m_JvtPass_Medium; if ( m_mc ) { delete m_JvtEff_SF_Medium; } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sfJVTName == "Tight" ) { delete m_JvtPass_Tight; if ( m_mc ) { delete m_JvtEff_SF_Tight; } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sfJVTName == "FixedEffPt" ) { delete m_JvtPass_FixedEffPt; if ( m_mc ) { delete m_JvtEff_SF_FixedEffPt; } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sfJVTName == "TightFwd" ) { delete m_JvtPass_TightFwd; if ( m_mc ) { delete m_JvtEff_SF_TightFwd; } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sffJVTName == "Loose" ) { delete m_fJvtPass_Loose; if ( m_mc ) { delete m_fJvtEff_SF_Loose; } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sffJVTName == "Medium" ) { delete m_fJvtPass_Medium; if ( m_mc ) { delete m_fJvtEff_SF_Medium; } } if ( m_infoSwitch.m_trackPV || m_infoSwitch.m_sffJVTName == "Tight" ) { delete m_fJvtPass_Tight; if ( m_mc ) { delete m_fJvtEff_SF_Tight; } } // chargedPFOPV if ( m_infoSwitch.m_chargedPFOPV ) { delete m_SumPtChargedPFOPt500PV; delete m_fCharged; } // allTrack // trackAll or trackPV if ( m_infoSwitch.m_allTrack ) { delete m_GhostTrackCount ; delete m_GhostTrackPt ; delete m_GhostTrack_pt ; delete m_GhostTrack_qOverP ; delete m_GhostTrack_eta ; delete m_GhostTrack_phi ; delete m_GhostTrack_e ; delete m_GhostTrack_d0 ; delete m_GhostTrack_z0 ; // allTrackDetail if(m_infoSwitch.m_allTrackDetail){ delete m_GhostTrack_nPixelHits ; delete m_GhostTrack_nSCTHits ; delete m_GhostTrack_nTRTHits ; delete m_GhostTrack_nPixelSharedHits ; delete m_GhostTrack_nPixelSplitHits ; delete m_GhostTrack_nInnermostPixelLayerHits ; delete m_GhostTrack_nInnermostPixelLayerSharedHits ; delete m_GhostTrack_nInnermostPixelLayerSplitHits ; delete m_GhostTrack_nNextToInnermostPixelLayerHits ; delete m_GhostTrack_nNextToInnermostPixelLayerSharedHits ; delete m_GhostTrack_nNextToInnermostPixelLayerSplitHits ; } } // constituent if ( m_infoSwitch.m_constituent ) { delete m_numConstituents; } if ( m_infoSwitch.m_constituentAll ) { delete m_constituentWeights; delete m_constituent_pt ; delete m_constituent_eta ; delete m_constituent_phi ; delete m_constituent_e ; } // flavorTag if( m_infoSwitch.m_flavorTag || m_infoSwitch.m_flavorTagHLT || m_infoSwitch.m_flavorTagTLA ) { delete m_HadronConeExclTruthLabelID; delete m_HadronConeExclExtendedTruthLabelID; } if( m_infoSwitch.m_flavorTag || m_infoSwitch.m_flavorTagHLT ) { // flavorTag delete m_DL1r; delete m_DL1r_pu; delete m_DL1r_pc; delete m_DL1r_pb; delete m_DL1dv00; delete m_DL1dv00_pu; delete m_DL1dv00_pc; delete m_DL1dv00_pb; delete m_DL1dv01; delete m_DL1dv01_pu; delete m_DL1dv01_pc; delete m_DL1dv01_pb; delete m_GN1; delete m_GN1_pu; delete m_GN1_pc; delete m_GN1_pb; delete m_GN2v00LegacyWP; delete m_GN2v00LegacyWP_pu; delete m_GN2v00LegacyWP_pc; delete m_GN2v00LegacyWP_pb; delete m_GN2v00NewAliasWP; delete m_GN2v00NewAliasWP_pu; delete m_GN2v00NewAliasWP_pc; delete m_GN2v00NewAliasWP_pb; delete m_GN2v01; delete m_GN2v01_pu; delete m_GN2v01_pc; delete m_GN2v01_pb; delete m_GN2v01_ptau; // Jet Fitter if( m_infoSwitch.m_jetFitterDetails){ delete m_JetFitter_nVTX ; delete m_JetFitter_nSingleTracks ; delete m_JetFitter_nTracksAtVtx ; delete m_JetFitter_mass ; delete m_JetFitter_energyFraction; delete m_JetFitter_significance3d; delete m_JetFitter_deltaeta ; delete m_JetFitter_deltaphi ; delete m_JetFitter_N2Tpar ; } // SV Details if( m_infoSwitch.m_svDetails){ delete m_SV0; delete m_sv0_NGTinSvx ; delete m_sv0_N2Tpair ; delete m_sv0_massvx ; delete m_sv0_efracsvx ; delete m_sv0_normdist ; delete m_SV1; delete m_SV1IP3D; delete m_COMBx; delete m_sv1_pb ; delete m_sv1_pc ; delete m_sv1_c ; delete m_sv1_cu ; delete m_sv1_NGTinSvx ; delete m_sv1_N2Tpair ; delete m_sv1_massvx ; delete m_sv1_efracsvx ; delete m_sv1_normdist ; delete m_sv1_Lxy ; delete m_sv1_sig3d ; delete m_sv1_L3d ; delete m_sv1_distmatlay; delete m_sv1_dR ; } // IP3D if( m_infoSwitch.m_ipDetails){ delete m_IP2D_pu ; delete m_IP2D_pb ; delete m_IP2D_pc ; delete m_IP2D ; delete m_IP2D_c ; delete m_IP2D_cu ; delete m_nIP2DTracks ; delete m_IP2D_gradeOfTracks ; delete m_IP2D_flagFromV0ofTracks; delete m_IP2D_valD0wrtPVofTracks; delete m_IP2D_sigD0wrtPVofTracks; delete m_IP2D_weightBofTracks ; delete m_IP2D_weightCofTracks ; delete m_IP2D_weightUofTracks ; delete m_IP3D; delete m_IP3D_pu ; delete m_IP3D_pb ; delete m_IP3D_pc ; delete m_IP3D_c ; delete m_IP3D_cu ; delete m_nIP3DTracks ; delete m_IP3D_gradeOfTracks ; delete m_IP3D_flagFromV0ofTracks ; delete m_IP3D_valD0wrtPVofTracks ; delete m_IP3D_sigD0wrtPVofTracks ; delete m_IP3D_valZ0wrtPVofTracks ; delete m_IP3D_sigZ0wrtPVofTracks ; delete m_IP3D_weightBofTracks ; delete m_IP3D_weightCofTracks ; delete m_IP3D_weightUofTracks ; } if( m_infoSwitch.m_JVC ) { delete m_JetVertexCharge_discriminant ; } } // flavorTagHLT if( m_infoSwitch.m_flavorTagHLT ) { delete m_vtxOnlineValid ; delete m_vtxHadDummy ; delete m_bs_online_vx ; delete m_bs_online_vy ; delete m_bs_online_vz ; delete m_vtx_offline_x0 ; delete m_vtx_offline_y0 ; delete m_vtx_offline_z0 ; delete m_vtx_online_x0 ; delete m_vtx_online_y0 ; delete m_vtx_online_z0 ; delete m_vtx_online_bkg_x0 ; delete m_vtx_online_bkg_y0 ; delete m_vtx_online_bkg_z0 ; } if( m_infoSwitch.m_flavorTagTLA ) { delete m_fastDIPS; delete m_fastDIPS_pu; delete m_fastDIPS_pc; delete m_fastDIPS_pb; } for(auto btag : m_btags) delete btag; // area if( m_infoSwitch.m_area ) { delete m_GhostArea ; delete m_ActiveArea ; delete m_VoronoiArea ; delete m_ActiveArea4vec_pt ; delete m_ActiveArea4vec_eta ; delete m_ActiveArea4vec_phi ; delete m_ActiveArea4vec_m ; } // truth if ( m_infoSwitch.m_truth && m_mc ) { delete m_ConeTruthLabelID; delete m_TruthCount; delete m_TruthLabelDeltaR_B; delete m_TruthLabelDeltaR_C; delete m_TruthLabelDeltaR_T; delete m_PartonTruthLabelID; delete m_GhostTruthAssociationFraction; delete m_truth_E; delete m_truth_pt; delete m_truth_phi; delete m_truth_eta; } // truth detail if ( m_infoSwitch.m_truthDetails ) { delete m_GhostBHadronsFinalCount ; delete m_GhostBHadronsInitialCount ; delete m_GhostBQuarksFinalCount ; delete m_GhostBHadronsFinalPt ; delete m_GhostBHadronsInitialPt ; delete m_GhostBQuarksFinalPt ; delete m_GhostCHadronsFinalCount ; delete m_GhostCHadronsInitialCount ; delete m_GhostCQuarksFinalCount ; delete m_GhostCHadronsFinalPt ; delete m_GhostCHadronsInitialPt ; delete m_GhostCQuarksFinalPt ; delete m_GhostTausFinalCount ; delete m_GhostTausFinalPt ; delete m_truth_pdgId ; delete m_truth_partonPt ; delete m_truth_partonDR ; } // charge if ( m_infoSwitch.m_charge ) { delete m_charge; } // passSel if ( m_infoSwitch.m_passSel ) { delete m_passSel; } // passOR if ( m_infoSwitch.m_passOR ) { delete m_passOR; } } void JetContainer::setTree(TTree *tree) { // // Connect branches ParticleContainer::setTree(tree); if(m_infoSwitch.m_rapidity) { connectBranch(tree,"rapidity", &m_rapidity); } if ( m_infoSwitch.m_trigger ){ connectBranch (tree, "isTrigMatched", &m_isTrigMatched); connectBranch >(tree, "isTrigMatchedToChain", &m_isTrigMatchedToChain ); connectBranch (tree, "listTrigChains", &m_listTrigChains ); } if(m_infoSwitch.m_clean || m_infoSwitch.m_cleanLight || m_infoSwitch.m_cleanLLP) { if(m_infoSwitch.m_clean){ connectBranch(tree, "Timing", &m_Timing); connectBranch(tree, "LArQuality", &m_LArQuality); connectBranch (tree, "LArBadHVNCell", &m_LArBadHVNCell); connectBranch(tree, "LArQuality", &m_LArQuality); connectBranch(tree, "HECQuality", &m_HECQuality); connectBranch(tree, "NegativeE", &m_NegativeE); connectBranch(tree, "AverageLArQF", &m_AverageLArQF); connectBranch(tree, "BchCorrCell", &m_BchCorrCell); connectBranch(tree, "N90Constituents", &m_N90Constituents); connectBranch(tree, "LArBadHVEnergyFrac", &m_LArBadHVEnergyFrac); connectBranch(tree, "ChargedFraction", &m_ChargedFraction); connectBranch(tree, "OotFracClusters5", &m_OotFracClusters5); connectBranch(tree, "OotFracClusters10", &m_OotFracClusters10); connectBranch(tree, "LeadingClusterPt", &m_LeadingClusterPt); connectBranch(tree, "LeadingClusterSecondLambda", &m_LeadingClusterSecondLambda); connectBranch(tree, "LeadingClusterCenterLambda", &m_LeadingClusterCenterLambda); connectBranch(tree, "LeadingClusterSecondR", &m_LeadingClusterSecondR); if(m_infoSwitch.m_cleanTrig) { connectBranch (tree, "clean_passLooseBadTriggerUgly",&m_clean_passLooseBadTriggerUgly); } else { connectBranch (tree, "clean_passLooseBadUgly", &m_clean_passLooseBadUgly); connectBranch (tree, "clean_passTightBadUgly", &m_clean_passTightBadUgly); } } if(m_infoSwitch.m_cleanTrig) { connectBranch (tree, "clean_passLooseBadTrigger", &m_clean_passLooseBadTrigger); } else { connectBranch (tree, "clean_passLooseBad", &m_clean_passLooseBad); connectBranch (tree, "clean_passTightBad", &m_clean_passTightBad); } if(m_infoSwitch.m_cleanLLP) { connectBranch (tree, "clean_passLooseBadLLP", &m_clean_passLooseBadLLP); } } if(m_infoSwitch.m_timing && !m_infoSwitch.m_clean){ connectBranch(tree, "Timing", &m_Timing); } if(m_infoSwitch.m_energy || m_infoSwitch.m_energyLight ) { if ( m_infoSwitch.m_energy ){ connectBranch(tree, "HECFrac", &m_HECFrac); connectBranch(tree, "CentroidR", &m_CentroidR); connectBranch(tree, "LowEtConstituentsFrac", &m_LowEtConstituentsFrac); } connectBranch(tree, "EMFrac", &m_EMFrac); connectBranch(tree, "FracSamplingMax", &m_FracSamplingMax); connectBranch(tree, "FracSamplingMaxIndex", &m_FracSamplingMaxIndex); connectBranch(tree, "GhostMuonSegmentCount", &m_GhostMuonSegmentCount); connectBranch(tree, "Width", &m_Width); } if(m_infoSwitch.m_trackPV) { connectBranch(tree, "NumTrkPt1000PV", &m_NumTrkPt1000PV); connectBranch(tree, "SumPtTrkPt1000PV", &m_SumPtTrkPt1000PV); connectBranch(tree, "TrackWidthPt1000PV", &m_TrackWidthPt1000PV); connectBranch(tree, "NumTrkPt500PV", &m_NumTrkPt500PV); connectBranch(tree, "SumPtTrkPt500PV", &m_SumPtTrkPt500PV); connectBranch(tree, "TrackWidthPt500PV", &m_TrackWidthPt500PV); connectBranch(tree, "JVFPV", &m_JVFPV); } if(m_infoSwitch.m_trackAll || m_infoSwitch.m_trackPV || m_infoSwitch.m_jvt) { if(m_infoSwitch.m_trackAll || m_infoSwitch.m_trackPV){ connectBranch(tree, "JvtJvfcorr", &m_JvtJvfcorr); connectBranch(tree, "JvtRpt", &m_JvtRpt); } connectBranch(tree, "Jvt", &m_Jvt); } if ( m_infoSwitch.m_chargedPFOPV ) { connectBranch(tree, "SumPtChargedPFOPt500PV", &m_SumPtChargedPFOPt500PV); connectBranch(tree, "fCharged", &m_fCharged); } if(m_infoSwitch.m_JVC) { connectBranch(tree,"JetVertexCharge_discriminant", &m_JetVertexCharge_discriminant); } if(m_infoSwitch.m_flavorTag || m_infoSwitch.m_flavorTagHLT || m_infoSwitch.m_flavorTagTLA) { connectBranch (tree,"HadronConeExclTruthLabelID" ,&m_HadronConeExclTruthLabelID); connectBranch (tree,"HadronConeExclExtendedTruthLabelID",&m_HadronConeExclExtendedTruthLabelID); } if(m_infoSwitch.m_flavorTag || m_infoSwitch.m_flavorTagHLT) { connectBranch(tree,"DL1r" ,&m_DL1r ); connectBranch(tree,"DL1r_pu" ,&m_DL1r_pu ); connectBranch(tree,"DL1r_pc" ,&m_DL1r_pc ); connectBranch(tree,"DL1r_pb" ,&m_DL1r_pb ); connectBranch(tree,"DL1dv00" ,&m_DL1dv00 ); connectBranch(tree,"DL1dv00_pu" ,&m_DL1dv00_pu ); connectBranch(tree,"DL1dv00_pc" ,&m_DL1dv00_pc ); connectBranch(tree,"DL1dv00_pb" ,&m_DL1dv00_pb ); connectBranch(tree,"DL1dv01" ,&m_DL1dv01 ); connectBranch(tree,"DL1dv01_pu" ,&m_DL1dv01_pu ); connectBranch(tree,"DL1dv01_pc" ,&m_DL1dv01_pc ); connectBranch(tree,"DL1dv01_pb" ,&m_DL1dv01_pb ); connectBranch(tree,"GN1" ,&m_GN1 ); connectBranch(tree,"GN1_pu" ,&m_GN1_pu ); connectBranch(tree,"GN1_pc" ,&m_GN1_pc ); connectBranch(tree,"GN1_pb" ,&m_GN1_pb ); connectBranch(tree,"GN2v00LegacyWP" ,&m_GN2v00LegacyWP ); connectBranch(tree,"GN2v00LegacyWP_pu" ,&m_GN2v00LegacyWP_pu ); connectBranch(tree,"GN2v00LegacyWP_pc" ,&m_GN2v00LegacyWP_pc ); connectBranch(tree,"GN2v00LegacyWP_pb" ,&m_GN2v00LegacyWP_pb ); connectBranch(tree,"GN2v00NewAliasWP" ,&m_GN2v00NewAliasWP ); connectBranch(tree,"GN2v00NewAliasWP_pu" ,&m_GN2v00NewAliasWP_pu ); connectBranch(tree,"GN2v00NewAliasWP_pc" ,&m_GN2v00NewAliasWP_pc ); connectBranch(tree,"GN2v00NewAliasWP_pb" ,&m_GN2v00NewAliasWP_pb ); connectBranch(tree,"GN2v01" ,&m_GN2v01 ); connectBranch(tree,"GN2v01_pu" ,&m_GN2v01_pu ); connectBranch(tree,"GN2v01_pc" ,&m_GN2v01_pc ); connectBranch(tree,"GN2v01_pb" ,&m_GN2v01_pb ); connectBranch(tree,"GN2v01_ptau" ,&m_GN2v01_ptau); } if(m_infoSwitch.m_flavorTagHLT) { connectBranch(tree,"vtxHadDummy", &m_vtxHadDummy); connectBranch(tree,"bs_online_vx", &m_bs_online_vx); connectBranch(tree,"bs_online_vy", &m_bs_online_vy); connectBranch(tree,"bs_online_vz", &m_bs_online_vz); connectBranch(tree,"vtx_offline_x0", &m_vtx_offline_x0); connectBranch(tree,"vtx_offline_y0", &m_vtx_offline_y0); connectBranch(tree,"vtx_offline_z0", &m_vtx_offline_z0); connectBranch(tree,"vtx_online_x0", &m_vtx_online_x0); connectBranch(tree,"vtx_online_y0", &m_vtx_online_y0); connectBranch(tree,"vtx_online_z0", &m_vtx_online_z0); connectBranch(tree,"vtx_online_bkg_x0", &m_vtx_online_bkg_x0); connectBranch(tree,"vtx_online_bkg_y0", &m_vtx_online_bkg_y0); connectBranch(tree,"vtx_online_bkg_z0", &m_vtx_online_bkg_z0); } if(m_infoSwitch.m_flavorTagTLA) { connectBranch(tree,"fastDIPS", &m_fastDIPS); connectBranch(tree,"fastDIPS_pu", &m_fastDIPS_pu); connectBranch(tree,"fastDIPS_pc", &m_fastDIPS_pc); connectBranch(tree,"fastDIPS_pb", &m_fastDIPS_pb); } if(m_infoSwitch.m_jetFitterDetails) { connectBranch(tree,"JetFitter_nVTX" , &m_JetFitter_nVTX ); connectBranch(tree,"JetFitter_nSingleTracks" , &m_JetFitter_nSingleTracks ); connectBranch(tree,"JetFitter_nTracksAtVtx" , &m_JetFitter_nTracksAtVtx ); connectBranch(tree,"JetFitter_mass" , &m_JetFitter_mass ); connectBranch(tree,"JetFitter_energyFraction", &m_JetFitter_energyFraction ); connectBranch(tree,"JetFitter_significance3d", &m_JetFitter_significance3d ); connectBranch(tree,"JetFitter_deltaeta" , &m_JetFitter_deltaeta ); connectBranch(tree,"JetFitter_deltaphi" , &m_JetFitter_deltaphi ); connectBranch(tree,"JetFitter_N2Tpair" , &m_JetFitter_N2Tpar ); } if( m_infoSwitch.m_svDetails){ connectBranch(tree, "SV0", &m_SV0); connectBranch(tree, "sv0_NGTinSvx", &m_sv0_NGTinSvx ); connectBranch(tree, "sv0_N2Tpair", &m_sv0_N2Tpair ); connectBranch(tree, "sv0_massvx", &m_sv0_massvx ); connectBranch(tree, "sv0_efracsvx", &m_sv0_efracsvx ); connectBranch(tree, "sv0_normdist", &m_sv0_normdist ); connectBranch(tree, "SV1", &m_SV1); connectBranch(tree, "SV1IP3D", &m_SV1IP3D); connectBranch(tree, "COMBx", &m_COMBx); connectBranch(tree, "sv1_pu", &m_sv1_pu ); connectBranch(tree, "sv1_pb", &m_sv1_pb ); connectBranch(tree, "sv1_pc", &m_sv1_pc ); connectBranch(tree, "sv1_c", &m_sv1_c ); connectBranch(tree, "sv1_cu", &m_sv1_cu ); connectBranch(tree, "sv1_NGTinSvx", &m_sv1_NGTinSvx ); connectBranch(tree, "sv1_N2Tpair", &m_sv1_N2Tpair ); connectBranch(tree, "sv1_massvx", &m_sv1_massvx ); connectBranch(tree, "sv1_efracsvx", &m_sv1_efracsvx ); connectBranch(tree, "sv1_normdist", &m_sv1_normdist ); connectBranch(tree, "sv1_Lxy", &m_sv1_Lxy ); connectBranch