From c020dafce74feb1bc555d5d65a7fdb0ffd53c789 Mon Sep 17 00:00:00 2001 From: Frank Winklmeier <fwinkl@cern> Date: Tue, 24 Nov 2020 12:58:42 +0100 Subject: [PATCH] SpecialUtils: delete obsolete package (ATLASG-1572) --- .../SpecialUtils/CMakeLists.txt | 21 ---- .../SpecialUtils/ElectronTypeDefs.h | 15 --- .../SpecialUtils/SpecialUtils/JetTypeDefs.h | 15 --- .../SpecialUtils/SpecialUtils/MuonTypeDefs.h | 15 --- .../SpecialUtils/NeutrinoTypeDefs.h | 15 --- .../SpecialUtils/SpecialUtils/NeutrinoUtils.h | 40 ------- .../SpecialUtils/NeutrinoUtils.icc | 103 ------------------ .../SpecialUtils/ParticleBaseTypeDefs.h | 15 --- .../SpecialUtils/PhotonTypeDefs.h | 15 --- .../SpecialUtils/TauJetTypeDefs.h | 15 --- .../SpecialUtils/doc/packagedoc.h | 25 ----- 11 files changed, 294 deletions(-) delete mode 100644 PhysicsAnalysis/AnalysisCommon/SpecialUtils/CMakeLists.txt delete mode 100755 PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/ElectronTypeDefs.h delete mode 100755 PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/JetTypeDefs.h delete mode 100755 PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/MuonTypeDefs.h delete mode 100755 PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoTypeDefs.h delete mode 100755 PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoUtils.h delete mode 100755 PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoUtils.icc delete mode 100755 PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/ParticleBaseTypeDefs.h delete mode 100755 PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/PhotonTypeDefs.h delete mode 100755 PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/TauJetTypeDefs.h delete mode 100644 PhysicsAnalysis/AnalysisCommon/SpecialUtils/doc/packagedoc.h diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/CMakeLists.txt b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/CMakeLists.txt deleted file mode 100644 index 586cc7f5d367..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/CMakeLists.txt +++ /dev/null @@ -1,21 +0,0 @@ -################################################################################ -# Package: SpecialUtils -################################################################################ - -# Declare the package name: -atlas_subdir( SpecialUtils ) - -# Declare the package's dependencies: -atlas_depends_on_subdirs( PUBLIC - PhysicsAnalysis/AnalysisCommon/ParticleEvent - Reconstruction/Jet/JetEvent - Reconstruction/MuonIdentification/muonEvent - Reconstruction/egamma/egammaEvent - Reconstruction/tauEvent ) - -# External dependencies: -find_package( CLHEP ) - -# Install files from the package: -atlas_install_headers( SpecialUtils ) - diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/ElectronTypeDefs.h b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/ElectronTypeDefs.h deleted file mode 100755 index 9484c2b1236e..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/ElectronTypeDefs.h +++ /dev/null @@ -1,15 +0,0 @@ -/* - Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration -*/ - -#ifndef ELECTRON_TYPE_DEFS_H -#define ELECTRON_TYPE_DEFS_H - -#include "egammaEvent/ElectronContainer.h" - -typedef std::vector<const ElectronContainer::base_value_type*> ElectronVect; - -typedef std::vector<ElectronContainer::base_value_type*> ListOfElectrons; - -#endif /// ELECTRON_TYPE_DEFS_H - diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/JetTypeDefs.h b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/JetTypeDefs.h deleted file mode 100755 index e399b6e4e7c9..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/JetTypeDefs.h +++ /dev/null @@ -1,15 +0,0 @@ -/* - Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration -*/ - -#ifndef PARTICLEJET_TYPE_DEFS_H -#define PARTICLEJET_TYPE_DEFS_H - -#include "JetEvent/JetCollection.h" - -typedef std::vector<const JetCollection::base_value_type*> JetVect; - -typedef std::vector<JetCollection::base_value_type*> ListOfJets; - -#endif /// PARTICLEJET_TYPE_DEFS_H - diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/MuonTypeDefs.h b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/MuonTypeDefs.h deleted file mode 100755 index daf3a1a349fc..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/MuonTypeDefs.h +++ /dev/null @@ -1,15 +0,0 @@ -/* - Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration -*/ - -#ifndef MUON_TYPE_DEFS_H -#define MUON_TYPE_DEFS_H - -#include "muonEvent/MuonContainer.h" - -typedef std::vector<const Analysis::MuonContainer::base_value_type*> MuonVect; - -typedef std::vector<Analysis::MuonContainer::base_value_type*> ListOfMuons; - -#endif /// MUON_TYPE_DEFS_H - diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoTypeDefs.h b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoTypeDefs.h deleted file mode 100755 index 113f4cd86a73..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoTypeDefs.h +++ /dev/null @@ -1,15 +0,0 @@ -/* - Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration -*/ - -#ifndef NEUTRINO_TYPE_DEFS_H -#define NEUTRINO_TYPE_DEFS_H - -#include "ParticleEvent/NeutrinoContainer.h" - -typedef std::vector<const NeutrinoContainer::base_value_type*> NeutrinoVect; - -typedef std::vector<NeutrinoContainer::base_value_type*> ListOfNeutrinos; - -#endif /// NEUTRINO_TYPE_DEFS_H - diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoUtils.h b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoUtils.h deleted file mode 100755 index ced88ad6296f..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoUtils.h +++ /dev/null @@ -1,40 +0,0 @@ -/* - Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration -*/ - -/** - Utilities for Neutrino - - @author Tadashi Maeno - */ - -#ifndef SPECIALUTILS_NEUTRINOUTILS_H -#define SPECIALUTILS_NEUTRINOUTILS_H - -namespace NeutrinoUtils { - - /** - construct neutrino 4-momontum from lepton, missing pt and the W-mass. - The neutrino mass is zero. This method returns the container of neutrino candidates. - */ - template <class LEPTON, class COLL> - bool candidatesFromWMass(LEPTON* lepton, const double pxMiss, const double pyMiss, - COLL & neutrinoContainer, bool DealWithNegative = false); - - /** - tau tau -> a+b+ptmiss utility. This method returns two neutrinos - given two particles and ptmiss. Here the constraint is that a and b are - collinear with the taus (either leptons or a TauJet. Thus 'a' carries a - fraction x_a of one tau's momentum and that tau's neutrinos carry (1-x_a) - of that tau's momentum. - */ - template <class LEPTON1, class LEPTON2, class COLL> - void neutrinosFromColinearApproximation(LEPTON1* particleA, LEPTON2* particleB, - const double pxMiss, const double pyMiss, - COLL & neutrinoContainer, - double &x_a, double &x_b); -} - -#include "SpecialUtils/NeutrinoUtils.icc" - -#endif diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoUtils.icc b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoUtils.icc deleted file mode 100755 index 2194271d5570..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/NeutrinoUtils.icc +++ /dev/null @@ -1,103 +0,0 @@ -/* - Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration -*/ - -#include "ParticleEvent/Neutrino.h" -#include "CLHEP/Units/SystemOfUnits.h" - -#include <cmath> - -template <class LEPTON, class COLL> inline -bool NeutrinoUtils::candidatesFromWMass(LEPTON* lepton, const double pxMiss, const double pyMiss, - COLL & neutrinoContainer, bool DealWithNegative) -{ - // W mass - const double WMass = 80.4*CLHEP::GeV; - - // clear output - neutrinoContainer.clear(); - - // solve the quadratic equation - - double ptMiss = sqrt (pxMiss*pxMiss + pyMiss*pyMiss); - - double alpha = pow(WMass,2)+pow((pxMiss+lepton->px()),2)+pow((pyMiss+lepton->py()),2) - -pow(lepton->e(),2); - - double beta = 0.5 * ( alpha-pow(ptMiss,2)+pow(lepton->pz(),2) ); - - double gamma = -( beta*beta - ( pow(lepton->e(),2)*pow(ptMiss,2) ) ) - / ( pow(lepton->e(),2)-pow(lepton->pz(),2) ); - - double lambda = 2*beta*lepton->pz() / (pow(lepton->e(),2)-pow(lepton->pz(),2)); - - double delta = pow(lambda,2)-4*gamma; - - // if no solution - if ( delta < 0 ){ - if ( DealWithNegative ){ - delta = 0; - } - else{ - return false; - } - } - - delta = sqrt(delta); - - // instantiate Neutrino - - double pz = (lambda-delta)/2.0; - double e = sqrt(pxMiss*pxMiss+pyMiss*pyMiss+pz*pz); - - Neutrino *nu1 = new Neutrino; - nu1->set4Mom(CLHEP::HepLorentzVector(pxMiss,pyMiss,pz,e)); - - neutrinoContainer.push_back(nu1); - if ( delta == 0 ) return true; - - pz = (lambda+delta)/2.0; - e = sqrt(pxMiss*pxMiss+pyMiss*pyMiss+pz*pz); - - Neutrino *nu2 = new Neutrino; - nu2->set4Mom(CLHEP::HepLorentzVector(pxMiss,pyMiss,pz,e)); - - neutrinoContainer.push_back(nu2); - - return true; -} - - -template <class LEPTON1, class LEPTON2, class COLL> inline -void NeutrinoUtils::neutrinosFromColinearApproximation(LEPTON1* particleA, LEPTON2* particleB, - const double pxMiss, const double pyMiss, - COLL & neutrinoContainer, - double &x_a, double &x_b) -{ - // numerator - double numerator = particleA->px()*particleB->py() - particleA->py()*particleB->px(); - - // fraction - x_a = numerator/(particleB->py()*(particleA->px()+pxMiss)-particleB->px()*(particleA->py()+pyMiss)); - x_b = numerator/(particleA->px()*(particleB->py()+pyMiss)-particleA->py()*(particleB->px()+pxMiss)); - - // instantiate Neutrino - double px = (1./x_a-1.)*particleA->px(); - double py = (1./x_a-1.)*particleA->py(); - double pz = (1./x_a-1.)*particleA->pz(); - double e = sqrt(px*px+py*py+pz*pz); - - Neutrino *nu1 = new Neutrino; - nu1->set4Mom(CLHEP::HepLorentzVector(px,py,pz,e)); - - px = (1./x_b-1.)*particleB->px(); - py = (1./x_b-1.)*particleB->py(); - pz = (1./x_b-1.)*particleB->pz(); - e = sqrt(px*px+py*py+pz*pz); - - Neutrino *nu2 = new Neutrino; - nu2->set4Mom(CLHEP::HepLorentzVector(px,py,pz,e)); - - neutrinoContainer.push_back(nu1); - neutrinoContainer.push_back(nu2); -} diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/ParticleBaseTypeDefs.h b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/ParticleBaseTypeDefs.h deleted file mode 100755 index d0794f48799e..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/ParticleBaseTypeDefs.h +++ /dev/null @@ -1,15 +0,0 @@ -/* - Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration -*/ - -#ifndef PARTICLEBASE_TYPE_DEFS_H -#define PARTICLEBASE_TYPE_DEFS_H - -#include "ParticleEvent/ParticleBaseContainer.h" - -typedef std::vector<const ParticleBaseContainer::base_value_type*> ParticleVect; - -typedef std::vector<ParticleBaseContainer::base_value_type*> ListOfParticles; - -#endif /// PARTICLEBASE_TYPE_DEFS_H - diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/PhotonTypeDefs.h b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/PhotonTypeDefs.h deleted file mode 100755 index cb63eaef0430..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/PhotonTypeDefs.h +++ /dev/null @@ -1,15 +0,0 @@ -/* - Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration -*/ - -#ifndef PHOTON_TYPE_DEFS_H -#define PHOTON_TYPE_DEFS_H - -#include "egammaEvent/PhotonContainer.h" - -typedef std::vector<const PhotonContainer::base_value_type*> PhotonVect; - -typedef std::vector<PhotonContainer::base_value_type*> ListOfPhotons; - -#endif /// PHOTON_TYPE_DEFS_H - diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/TauJetTypeDefs.h b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/TauJetTypeDefs.h deleted file mode 100755 index 4a8c9ae33f70..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/SpecialUtils/TauJetTypeDefs.h +++ /dev/null @@ -1,15 +0,0 @@ -/* - Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration -*/ - -#ifndef TAUJET_TYPE_DEFS_H -#define TAUJET_TYPE_DEFS_H - -#include "tauEvent/TauJetContainer.h" - -typedef std::vector<const Analysis::TauJetContainer::base_value_type*> TauJetVect; - -typedef std::vector<Analysis::TauJetContainer::base_value_type*> ListOfTauJets; - -#endif /// TAUJET_TYPE_DEFS_H - diff --git a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/doc/packagedoc.h b/PhysicsAnalysis/AnalysisCommon/SpecialUtils/doc/packagedoc.h deleted file mode 100644 index 5d76aced634b..000000000000 --- a/PhysicsAnalysis/AnalysisCommon/SpecialUtils/doc/packagedoc.h +++ /dev/null @@ -1,25 +0,0 @@ -/* - Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration -*/ - -/** -@page SpecialUtils_page SpecialUtils - -@section SpecialUtils_introductionSpecialUtils Introduction -This package contains non-general tools, i.e., tools for Neutrino. General tools should be put in AnalysisTools/Utils. - -@section SpecialUtils_packagecontentSpecialUtils Package Contents -SpecialUtils contains the following classes: - -- NeutrinoUtils : tools for Neutrino. There 2 utilities here. One where the W mass constraint is used to find -the longitudinal component of the neutrino momentum, for example in W --> e nu, assuming that the transverse -missing energy is carried away by the neutrino. Up to 2 solutions can be obtained. The other is where one uses -the collinear approximation to solve for the neutrinos in the decay X --> a + b + Missing Energy, for example Z ---> tautau -> tauJet + tauJet + Missing Energy; The constraint is that a and b are collinear with the taus -(either leptons or a TauJet. Thus, 'a' carries a fraction x_a of one tau's momentum and that tau's neutrinos -carry (1-x_a) of that tau's momentum. Both of these 2 utilities return a container of Neutrino Objects. There is -a class for the neutrino, Netrino.h, in the package PhyscisAnalysis/AnalysisCommon/ParticleEvent. - - - -*/ -- GitLab