@@ -124,10 +124,6 @@ class L1TrackObjectNtupleMaker : public edm::one::EDAnalyzer<edm::one::SharedRes
124124 explicit L1TrackObjectNtupleMaker (const edm::ParameterSet& iConfig);
125125 ~L1TrackObjectNtupleMaker () override ;
126126
127- template <typename T>
128- bool findHiggsToMuAncestor (T particle);
129- template <typename T>
130- bool findHiggsToBAncestor (T particle);
131127 template <typename T>
132128 bool isHard (T particle);
133129 double DoublePtFromBits (const L1Track&) const ;
@@ -207,7 +203,6 @@ class L1TrackObjectNtupleMaker : public edm::one::EDAnalyzer<edm::one::SharedRes
207203 edm::InputTag RecoVertexEmuInputTag;
208204 edm::InputTag GenParticleInputTag;
209205 edm::InputTag DisplacedVertexInputTag;
210- edm::InputTag DisplacedVertexEmulationInputTag;
211206
212207 edm::InputTag TrackFastJetsInputTag;
213208 edm::InputTag TrackJetsInputTag;
@@ -673,7 +668,6 @@ L1TrackObjectNtupleMaker::L1TrackObjectNtupleMaker(edm::ParameterSet const& iCon
673668 GenParticleInputTag = iConfig.getParameter <InputTag>(" GenParticleInputTag" );
674669 SimVertexInputTag = iConfig.getParameter <InputTag>(" SimVertexInputTag" );
675670 DisplacedVertexInputTag = iConfig.getParameter <InputTag>(" DisplacedVertexInputTag" );
676- DisplacedVertexEmulationInputTag = iConfig.getParameter <InputTag>(" DisplacedVertexEmulationInputTag" );
677671
678672 if (Displaced == " Prompt" || Displaced == " Both" ) {
679673 L1TrackInputTag = iConfig.getParameter <edm::InputTag>(" L1TrackInputTag" );
@@ -1133,63 +1127,6 @@ void L1TrackObjectNtupleMaker::endJob() {
11331127 delete m_trkfastjetExt_truetp_sumpt;
11341128}
11351129
1136- template <typename T>
1137- bool L1TrackObjectNtupleMaker::findHiggsToMuAncestor (T particle) {
1138- if ((particle->pdgId () == 13 || particle->pdgId () == -13 ) && particle->genParticles ().size () > 0 ) {
1139- reco::GenParticleRef genPart = particle->genParticles ()[0 ];
1140- reco::GenParticleRefVector parentParts = genPart->motherRefVector ();
1141- if (parentParts.size () > 0 ) {
1142- while (parentParts[0 ]->pdgId () == 13 || parentParts[0 ]->pdgId () == -13 )
1143- parentParts = parentParts[0 ]->motherRefVector ();
1144- reco::GenParticleRefVector daughters = parentParts[0 ]->daughterRefVector ();
1145- bool hasMuon = false ;
1146- bool hasAntiMuon = false ;
1147- for (auto daughter : daughters) {
1148- if (daughter->pdgId () == 13 )
1149- hasMuon = true ;
1150- if (daughter->pdgId () == -13 )
1151- hasAntiMuon = true ;
1152- }
1153- if (hasMuon && hasAntiMuon) {
1154- bool hAncestor = false ;
1155- while (parentParts.size () > 0 ) {
1156- if (parentParts[0 ]->pdgId () == 25 )
1157- hAncestor = true ;
1158- parentParts = parentParts[0 ]->motherRefVector ();
1159- }
1160- if (hAncestor) {
1161- return true ;
1162- }
1163- }
1164- }
1165- } // mu or anti mu
1166- return false ;
1167- }
1168-
1169- template <typename T>
1170- bool L1TrackObjectNtupleMaker::findHiggsToBAncestor (T particle) {
1171- TrackingVertexRef parentVert = particle->parentVertex ();
1172- TrackingParticle currentParticle = *particle;
1173- while (currentParticle.genParticles ().size () == 0 && parentVert->nSourceTracks () > 0 ) {
1174- TrackingParticleRefVector sourceTPs = parentVert->sourceTracks ();
1175- currentParticle = *sourceTPs[0 ];
1176- parentVert = currentParticle.parentVertex ();
1177- }
1178- bool bAncestor = false ;
1179- bool hAncestor = false ;
1180- reco::GenParticleRefVector genParticles = currentParticle.genParticles ();
1181- while (genParticles.size () > 0 ) {
1182- if (hAncestor == false && abs (genParticles[0 ]->pdgId ()) == 5 )
1183- bAncestor = true ;
1184- if (genParticles[0 ]->pdgId () == 25 )
1185- hAncestor = true ;
1186- genParticles = genParticles[0 ]->motherRefVector ();
1187- }
1188- if (bAncestor && hAncestor)
1189- return true ;
1190- return false ;
1191- }
1192-
11931130template <typename T>
11941131bool L1TrackObjectNtupleMaker::isHard (T particle) {
11951132 reco::GenParticleRefVector genParts = particle->genParticles ();
0 commit comments