RegionalCaloTrigger/test/testL1RCTJetSummaryCard.cpp

0001 #include "L1Trigger/RegionalCaloTrigger/interface/L1RCTJetSummaryCard.h"
0002 #include "L1Trigger/RegionalCaloTrigger/interface/L1RCTLookupTables.h"
0003 #include "CondFormats/L1TObjects/interface/L1RCTParameters.h"
0004
0005 #include <vector>
0006 using std::vector;
0007
0008 int main() {
0009   // For testing use 1:1 LUT
0010   std::vector<double> eGammaECalScaleFactors(32, 1.0);
0011   std::vector<double> eGammaHCalScaleFactors(32, 1.0);
0012   std::vector<double> jetMETECalScaleFactors(32, 1.0);
0013   std::vector<double> jetMETHCalScaleFactors(32, 1.0);
0014   std::vector<double> c,d,e,f,g,h;
0015   L1RCTParameters* rctParameters =
0016     new L1RCTParameters(1.0,                       // eGammaLSB
0017             1.0,                       // jetMETLSB
0018             3.0,                       // eMinForFGCut
0019             40.0,                      // eMaxForFGCut
0020             0.5,                       // hOeCut
0021             1.0,                       // eMinForHoECut
0022             50.0,                      // eMaxForHoECut
0023             1.0,                       // hMinForHoECut
0024             2.0,                       // eActivityCut
0025             3.0,                       // hActivityCut
0026             3,                         // eicIsolationThreshold
0027                         3,                         // jscQuietThresholdBarrel
0028                         3,                         // jscQuietThresholdEndcap
0029             false,                     // noiseVetoHB
0030             false,                     // noiseVetoHEplus
0031             false,                     // noiseVetoHEminus
0032             false,                     // use Lindsey
0033             eGammaECalScaleFactors,
0034             eGammaHCalScaleFactors,
0035             jetMETECalScaleFactors,
0036             jetMETHCalScaleFactors,
0037             c,
0038             d,
0039             e,
0040             f,
0041             g,
0042             h
0043             );
0044   L1RCTLookupTables* lut = new L1RCTLookupTables();
0045   lut->setRCTParameters(rctParameters);  // transcoder and etScale are not used
0046   L1RCTJetSummaryCard jsc(0,lut);
0047   std::vector<unsigned short> hfregions(8);
0048   std::vector<unsigned short> bregions(14);
0049   std::vector<unsigned short> tauBits(14);
0050   std::vector<unsigned short> mipBits(14);
0051   std::vector<unsigned short> isoElectrons(14);
0052   std::vector<unsigned short> nonIsoElectrons(14);
0053   isoElectrons.at(0) = 10;
0054   isoElectrons.at(1) = 20;
0055   isoElectrons.at(2) = 30;
0056   isoElectrons.at(3) = 40;
0057   isoElectrons.at(4) = 50;
0058   nonIsoElectrons.at(0) = 80;
0059   nonIsoElectrons.at(1) = 35;
0060   nonIsoElectrons.at(2) = 92;
0061   nonIsoElectrons.at(3) = 50;
0062   nonIsoElectrons.at(4) = 49;
0063   nonIsoElectrons.at(5) = 34;
0064   mipBits.at(0) = 1;
0065   mipBits.at(1) = 1;
0066   mipBits.at(10) = 1;
0067   bregions.at(0) = 100;
0068   bregions.at(2) = 50;
0069   bregions.at(12) = 50;
0070   jsc.fillMIPBits(mipBits);
0071   jsc.fillTauBits(tauBits);
0072   jsc.fillNonIsolatedEGObjects(nonIsoElectrons);
0073   jsc.fillIsolatedEGObjects(isoElectrons);
0074   jsc.fillRegionSums(bregions);
0075   jsc.fillHFRegionSums(hfregions);
0076   jsc.fillQuietBits();
0077   jsc.fillJetRegions();
0078   jsc.print();
0079 }