File indexing completed on 2023-03-17 10:47:08
0001 void TestCorrections(double rawPt, double rawEta, double rawPhi, double rawE, double JPTE, int NPV) {
0002 gROOT->ProcessLine("#include <vector>");
0003 gSystem->Load("libFWCoreFWLite.so");
0004 FWLiteEnabler::enable();
0005
0006 string L1Tag = "../data/Jec11V0_L1Offset_AK5JPT.txt";
0007 string L1JPTTag = "../data/Jec11V0_L1JPTOffset_AK5JPT.txt";
0008 JetCorrectorParameters *L1Par = new JetCorrectorParameters(L1Tag);
0009 JetCorrectorParameters *L1JPTPar = new JetCorrectorParameters(L1JPTTag);
0010 vector<JetCorrectorParameters> vPar;
0011 vPar.push_back(*L1Par);
0012 vPar.push_back(*L1JPTPar);
0013
0014 FactorizedJetCorrector *JetCorrector = new FactorizedJetCorrector(vPar);
0015
0016 TLorentzVector rawJetP4;
0017 rawJetP4.SetPtEtaPhiE(rawPt, rawEta, rawPhi, rawE);
0018 JetCorrector->setJPTrawP4(rawJetP4);
0019 JetCorrector->setJetE(JPTE);
0020 JetCorrector->setNPV(NPV);
0021 vector<float> vcor;
0022 vcor = JetCorrector->getSubCorrections();
0023 cout << "Correction applied to JPT jet after L1Offset = " << vcor[0] << endl;
0024 cout << "Correction applied to JPT jet after L1JPTOffset = " << vcor[1] << endl;
0025 }