Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

New MET correction/uncertainty tool #9759

Merged
merged 16 commits into from Jun 27, 2015
Merged

New MET correction/uncertainty tool #9759

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
16 commits
Select commit Hold shift + click to select a range
10756c5
Fixing the few conflicts appeared during the development
Apr 9, 2015
fcf2d38
First commit of multiplicity based corrections
schoef May 24, 2014
6ac9bee
modification of MET std correction modules and XY correction module
schoef Feb 27, 2015
2c85c7c
modification of MET std correction modules and XY correction module, …
May 13, 2015
0151767
modification of MET std correction modules and XY correction module, …
May 13, 2015
c32841d
finish first migration
May 13, 2015
a470357
new MET uncertainty tool, part 2
May 13, 2015
e10dcb0
change of pat::MET format
May 13, 2015
7e5dce4
first working version, cleaning not yet done
May 13, 2015
ad399a1
last modifications, check for 74X compatibility, cleaning
May 16, 2015
c3bfb9a
cleaning
May 16, 2015
c596467
fixed JetMETCorrections/Type1MET/test/testMultShiftMETcorrInputProduc…
schoef May 17, 2015
e42570c
fix pfmet correction test, removed unphysical T0T2 MET flavors
May 18, 2015
b312f31
fix corrMET_pfMet_cfg test
May 19, 2015
8250fdd
added warning for negative JEC correction potentially noticed in the …
May 20, 2015
3607d86
fix multithread declaration, class version and TF1 pointers
Jun 26, 2015
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
91 changes: 76 additions & 15 deletions DataFormats/PatCandidates/interface/MET.h
View file
Open in desktop
Expand Up @@ -54,6 +54,8 @@ namespace pat {
MET(const edm::Ptr<reco::MET> & aMETRef);
/// copy constructor
MET( MET const&);
/// cosntructor for corrected METs (keeping specific informations from src MET)
MET(const reco::MET & corMET, const MET& srcMET );
/// destructor
virtual ~MET();

Expand All @@ -68,6 +70,13 @@ namespace pat {
/// set the associated GenMET
void setGenMET(const reco::GenMET & gm);

// ----- MET significance functions ----
// set the MET Significance
void setMETSignificance(const double& metSig);
// get the MET significance
double metSignificance() const;


// ---- methods for uncorrected MET ----
float uncorrectedPt() const;
float uncorrectedPhi() const;
Expand Down Expand Up @@ -105,7 +114,7 @@ namespace pat {
/// Returns the event electromagnetic energy extracted from HF
double emEtInHF() const {return caloSpecific().EmEtInHF;}
/// Returns the event MET Significance
double metSignificance() const {return caloSpecific().METSignificance;}
double caloMetSignificance() const {return caloSpecific().METSignificance;}
/// Returns the event SET in HF+
double CaloSETInpHF() const {return caloSpecific().CaloSETInpHF;}
/// Returns the event SET in HF-
Expand Down Expand Up @@ -140,28 +149,57 @@ namespace pat {

// ---- members for MET corrections ----
enum METUncertainty {
JetEnUp=0, JetEnDown=1, JetResUp=2, JetResDown=3,
JetResUp=0, JetResDown=1, JetEnUp=2, JetEnDown=3,
MuonEnUp=4, MuonEnDown=5, ElectronEnUp=6, ElectronEnDown=7, TauEnUp=8,TauEnDown=9,
UnclusteredEnUp=10,UnclusteredEnDown=11, NoShift=12, METUncertaintySize=13
UnclusteredEnUp=10,UnclusteredEnDown=11, NoShift=12, METUncertaintySize=13,
JetResUpSmear=14, JetResDownSmear=15, METFullUncertaintySize=16
};
enum METUncertaintyLevel {
Raw=0, Type1=1, Type1p2=2, Calo=3
enum METCorrectionLevel {
Raw=0, Type1=1, Type01=2, TypeXY=3, Type1XY=4, Type01XY=5,
Type1Smear=6, Type01Smear=7, Type1SmearXY=8,
Type01SmearXY=9, RawCalo=10, METCorrectionLevelSize=11
};
enum METCorrectionType {
None=0, T1=1, T0=2, TXY=3, TXYForRaw=4,
TXYForT01=5, TXYForT1Smear=6, TXYForT01Smear=7,
Smear=8, Calo=9, METCorrectionTypeSize=10
};

struct Vector2 {
double px, py;
double pt() const { return hypotf(px,py); }
double phi() const { return std::atan2(py,px); }
};
Vector2 shiftedP2(METUncertainty shift, METUncertaintyLevel level=Type1) const ;
Vector shiftedP3(METUncertainty shift, METUncertaintyLevel level=Type1) const ;
LorentzVector shiftedP4(METUncertainty shift, METUncertaintyLevel level=Type1) const ;
double shiftedPx(METUncertainty shift, METUncertaintyLevel level=Type1) const { return shiftedP2(shift,level).px; }
double shiftedPy(METUncertainty shift, METUncertaintyLevel level=Type1) const { return shiftedP2(shift,level).py; }
double shiftedPt(METUncertainty shift, METUncertaintyLevel level=Type1) const { return shiftedP2(shift,level).pt(); }
double shiftedPhi(METUncertainty shift, METUncertaintyLevel level=Type1) const { return shiftedP2(shift,level).phi(); }
double shiftedSumEt(METUncertainty shift, METUncertaintyLevel level=Type1) const ;

void setShift(double px, double py, double sumEt, METUncertainty shift, METUncertaintyLevel level=Type1) ;
Vector2 shiftedP2(METUncertainty shift, METCorrectionLevel level=Type1) const ;
Vector shiftedP3(METUncertainty shift, METCorrectionLevel level=Type1) const ;
LorentzVector shiftedP4(METUncertainty shift, METCorrectionLevel level=Type1) const ;
double shiftedPx(METUncertainty shift, METCorrectionLevel level=Type1) const { return shiftedP2(shift,level).px; };
double shiftedPy(METUncertainty shift, METCorrectionLevel level=Type1) const { return shiftedP2(shift,level).py; };
double shiftedPt(METUncertainty shift, METCorrectionLevel level=Type1) const { return shiftedP2(shift,level).pt(); };
double shiftedPhi(METUncertainty shift, METCorrectionLevel level=Type1) const { return shiftedP2(shift,level).phi(); };
double shiftedSumEt(METUncertainty shift, METCorrectionLevel level=Type1) const ;

Vector2 corP2(METCorrectionLevel level=Type1) const ;
Vector corP3(METCorrectionLevel level=Type1) const ;
LorentzVector corP4(METCorrectionLevel level=Type1) const ;
double corPx(METCorrectionLevel level=Type1) const { return corP2(level).px; };
double corPy(METCorrectionLevel level=Type1) const { return corP2(level).py; };
double corPt(METCorrectionLevel level=Type1) const { return corP2(level).pt(); };
double corPhi(METCorrectionLevel level=Type1) const { return corP2(level).phi(); };
double corSumEt(METCorrectionLevel level=Type1) const ;

Vector2 uncorP2() const;
Vector uncorP3() const ;
LorentzVector uncorP4() const ;
double uncorPx() const { return uncorP2().px; };
double uncorPy() const { return uncorP2().py; };
double uncorPt() const { return uncorP2().pt(); };
double uncorPhi() const { return uncorP2().phi(); };
double uncorSumEt() const;

void setUncShift(double px, double py, double sumEt, METUncertainty shift, bool isSmeared=false);
void setCorShift(double px, double py, double sumEt, METCorrectionType level);

// specific method to fill and retrieve the caloMET quickly from miniAODs,
//should be used by JetMET experts only for the beginning
Expand All @@ -172,6 +210,14 @@ namespace pat {
double caloMETPhi() const;
double caloMETSumEt() const;

//Functions for backward compatibility with 74X samples.
// allow access to 74X samples, transparent and not used in 75
Vector2 shiftedP2_74x(METUncertainty shift, METCorrectionLevel level) const ;
Vector shiftedP3_74x(METUncertainty shift, METCorrectionLevel level) const ;
LorentzVector shiftedP4_74x(METUncertainty shift, METCorrectionLevel level) const ;
double shiftedSumEt_74x(METUncertainty shift, METCorrectionLevel level) const ;


/// this below should be private but Reflex doesn't like it
class PackedMETUncertainty {
// defined as C++ class so that I can change the packing without having to touch the code elsewhere
Expand All @@ -183,6 +229,7 @@ namespace pat {
double dpy() const { if(!unpacked_) unpack(); return dpy_; }
double dsumEt() const { if(!unpacked_) unpack(); return dsumEt_; }
void set(float dpx, float dpy, float dsumEt) { dpx_ = dpx; dpy_ = dpy; dsumEt_ = dsumEt; pack(); unpack();}
void add(float dpx, float dpy, float dsumEt) { dpx_ += dpx; dpy_ += dpy; dsumEt_ += dsumEt; }
void unpack() const ;
void pack();

Expand All @@ -200,11 +247,25 @@ namespace pat {
// ---- holder for pfMET specific info ---
std::vector<SpecificPFMETData> pfMET_;

// MET significance
double metSig_;

const PackedMETUncertainty findMETTotalShift(MET::METCorrectionLevel cor, MET::METUncertainty shift) const;

std::map<MET::METCorrectionLevel, std::vector<MET::METCorrectionType> > corMap_;
void initCorMap();



protected:

// ---- non-public correction utilities ----
//kept for 74X backward-compatibility, not initialized and used in 75X
std::vector<PackedMETUncertainty> uncertaintiesRaw_, uncertaintiesType1_, uncertaintiesType1p2_;


std::vector<PackedMETUncertainty> uncertainties_;
std::vector<PackedMETUncertainty> corrections_;

PackedMETUncertainty caloPackedMet_;

};
Expand Down