File indexing completed on 2024-04-06 12:19:58
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016 #include "L1Trigger/GlobalMuonTrigger/test/L1MuGMTPattern.h"
0017
0018
0019
0020
0021
0022 #include <fstream>
0023 #include <iomanip>
0024 #include <cmath>
0025
0026
0027
0028
0029 #include "FWCore/Framework/interface/Frameworkfwd.h"
0030 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0031
0032 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTReadoutCollection.h"
0033
0034 using namespace std;
0035
0036
0037
0038
0039 L1MuGMTPattern::L1MuGMTPattern(const edm::ParameterSet& ps) {
0040 m_inputTag = ps.getUntrackedParameter<edm::InputTag>("GMTInputTag", edm::InputTag("gmt"));
0041 m_inputCaloTag = ps.getUntrackedParameter<edm::InputTag>("GMTInputCaloTag", edm::InputTag("source"));
0042 m_outfilename = ps.getUntrackedParameter<string>("OutputFile", "gmt_testfile.dat");
0043 m_outputType = ps.getUntrackedParameter<int>("OutputType", 1);
0044 }
0045
0046
0047
0048
0049 L1MuGMTPattern::~L1MuGMTPattern() {}
0050
0051 void L1MuGMTPattern::beginJob() {
0052 ofstream of(m_outfilename.c_str(), ios::trunc);
0053 if (m_outputType == 1) {
0054 of << "-- GMT HW Test file generated by L1MuGMTPattern" << endl;
0055 of << "--" << endl;
0056 }
0057 of.close();
0058 if (m_outputType == 4) {
0059 of.open((m_outfilename + string("DT")).c_str(), ios::trunc);
0060 of.close();
0061 of.open((m_outfilename + string("CSC")).c_str(), ios::trunc);
0062 of.close();
0063 of.open((m_outfilename + string("BRPC")).c_str(), ios::trunc);
0064 of.close();
0065 of.open((m_outfilename + string("FRPC")).c_str(), ios::trunc);
0066 of.close();
0067 }
0068 }
0069
0070 void L1MuGMTPattern::endJob() {}
0071
0072
0073
0074
0075
0076 void L1MuGMTPattern::analyze(const edm::Event& e, const edm::EventSetup& es) {
0077
0078 if (m_outputType == 1) {
0079 ofstream of(m_outfilename.c_str(), ios::app);
0080 int runn = e.id().run();
0081 int eventn = e.id().event();
0082 of << "--" << endl;
0083 of << "RUN " << runn << endl;
0084 of << "EVT " << eventn << endl;
0085 of.close();
0086 }
0087
0088
0089
0090 edm::Handle<L1MuGMTReadoutCollection> gmtrc_handle;
0091 e.getByLabel(m_inputTag, gmtrc_handle);
0092 L1MuGMTReadoutCollection const* gmtrc = gmtrc_handle.product();
0093
0094 L1MuGMTReadoutRecord gmt_record = gmtrc->getRecord(0);
0095
0096
0097 printRegional("DT", gmt_record.getDTBXCands());
0098 printRegional("CSC", gmt_record.getCSCCands());
0099 printRegional("BRPC", gmt_record.getBrlRPCCands());
0100 printRegional("FRPC", gmt_record.getFwdRPCCands());
0101
0102
0103
0104
0105 edm::Handle<L1CaloRegionCollection> calocoll_h;
0106 e.getByLabel(m_inputCaloTag, calocoll_h);
0107 if (calocoll_h.isValid()) {
0108 printMipIso(calocoll_h.product());
0109 } else {
0110 edm::LogVerbatim("GMT_PSB_info") << " Calorimeter MIP/QUIET bits not found in the Event ";
0111 }
0112
0113
0114 printGMT("GMT", gmt_record.getGMTCands());
0115 printGMT("GMTB", gmt_record.getGMTBrlCands());
0116 printGMT("GMTF", gmt_record.getGMTFwdCands());
0117
0118 edm::Handle<vector<unsigned> > mi_h;
0119 e.getByLabel(m_inputTag, mi_h);
0120 if (mi_h.isValid()) {
0121 printMI(mi_h.product());
0122 } else {
0123 edm::LogVerbatim("GMT_PSB_info") << " MI product invalid ";
0124 }
0125
0126 printCANC();
0127 }
0128
0129 void L1MuGMTPattern::printRegional(string tag, const vector<L1MuRegionalCand>& rmc) {
0130 ofstream of;
0131 if (m_outputType == 4) {
0132 of.open((m_outfilename + tag).c_str(), ios::app);
0133 } else {
0134 of.open(m_outfilename.c_str(), ios::app);
0135 }
0136
0137 int i = 0;
0138 vector<L1MuRegionalCand>::const_iterator iter1;
0139 edm::LogVerbatim("GMTPattern") << "L1 " << tag.data() << " TRIGGER muons (bx0): " << endl;
0140 for (iter1 = rmc.begin(); iter1 != rmc.end(); iter1++) {
0141 if (i > 3)
0142 break;
0143 edm::LogVerbatim("GMTPattern") << " " << ++i << " : ";
0144 (*iter1).print();
0145
0146 if (m_outputType == 1 && !((*iter1).empty())) {
0147 of << left << setw(4) << tag.data() << right << " " << setw(1) << (*iter1).charge_valid_packed() << " " << setw(1)
0148 << (*iter1).charge_packed() << " " << setw(1) << (*iter1).finehalo_packed() << " " << setw(2)
0149 << (*iter1).eta_packed() << " " << setw(1) << (*iter1).quality_packed() << " " << setw(2)
0150 << (*iter1).pt_packed() << " " << setw(3) << (*iter1).phi_packed();
0151
0152 of << " -- " << left << setw(4) << tag.data() << right << " vch=" << (*iter1).charge_valid_packed()
0153 << " ch=" << (*iter1).chargeValue() << " hf=" << (*iter1).finehalo_packed() << " eta=" << (*iter1).etaValue()
0154 << " q=" << (*iter1).quality() << " pt=" << (*iter1).ptValue() << " phi=" << (*iter1).phiValue() << endl;
0155 }
0156 if (m_outputType == 2 || m_outputType == 4) {
0157 of << " " << hex << setfill('0') << setw(8) << invertQPt((*iter1).getDataWord()) << dec;
0158 }
0159 }
0160 if (m_outputType == 2 || m_outputType == 4) {
0161 for (; i < 4; i++) {
0162 of << " 0000ffff";
0163 }
0164 }
0165 if (m_outputType == 2) {
0166 of << " |";
0167 }
0168 if (m_outputType == 4) {
0169 of << endl;
0170 }
0171 }
0172
0173 void L1MuGMTPattern::printGMT(string tag, const vector<L1MuGMTExtendedCand>& exc) {
0174 ofstream of(m_outfilename.c_str(), ios::app);
0175 int i = 0;
0176 vector<L1MuGMTExtendedCand>::const_iterator gmt_iter;
0177 const char* sysign_txt[4] = {"POS ", "NEG ", "UNDEF", "SYNC "};
0178 edm::LogVerbatim("GMTPattern") << "L1 " << tag.data() << " muons (bx0): ";
0179 for (gmt_iter = exc.begin(); gmt_iter != exc.end(); gmt_iter++) {
0180 if (i > 3)
0181 break;
0182 edm::LogVerbatim("GMTPattern") << " " << ++i << " : ";
0183 (*gmt_iter).print();
0184
0185 if (m_outputType == 1 && !((*gmt_iter).empty())) {
0186 of << left << setw(4) << tag.data() << right << " " << setw(1) << (*gmt_iter).sysign() << " " << setw(1)
0187 << (*gmt_iter).mip() << " " << setw(1) << (*gmt_iter).isol() << " " << setw(2) << (*gmt_iter).etaIndex() << " "
0188 << setw(1) << (*gmt_iter).quality() << " " << setw(2) << (*gmt_iter).ptIndex() << " " << setw(3)
0189 << (*gmt_iter).phiIndex() << " " << setw(3) << (*gmt_iter).rank();
0190
0191 of << " -- " << left << setw(4) << tag.data() << right << " sy_sign/charge=" << sysign_txt[(*gmt_iter).sysign()]
0192 << " MIP=" << (*gmt_iter).mip() << " ISO=" << (*gmt_iter).isol() << " eta=" << (*gmt_iter).etaValue()
0193 << " q=" << (*gmt_iter).quality() << " pt=" << (*gmt_iter).ptValue() << " phi=" << (*gmt_iter).phiValue()
0194 << " rank=" << (*gmt_iter).rank() << endl;
0195 }
0196 if (m_outputType > 1 && tag == "GMT") {
0197 of << " " << hex << setfill('0') << setw(8) << (*gmt_iter).getDataWord() << dec;
0198 }
0199 }
0200 if (m_outputType > 1 && tag == "GMT") {
0201 for (; i < 4; i++) {
0202 of << " 00000000";
0203 }
0204 of << endl;
0205 }
0206 }
0207
0208 void L1MuGMTPattern::printMipIso(L1CaloRegionCollection const* regions) {
0209 int nmip = 0;
0210 int nnq = 0;
0211
0212 if (m_outputType == 1) {
0213 ofstream of(m_outfilename.c_str(), ios::app);
0214 L1CaloRegionCollection::const_iterator iter;
0215
0216 edm::LogVerbatim("GMT_PSB_info") << "MIP/QUIET bits rceived by the GMT :";
0217
0218 for (iter = regions->begin(); iter != regions->end(); iter++) {
0219 if ((*iter).id().ieta() < 4 || (*iter).id().ieta() > 17 || (*iter).id().iphi() > 17)
0220 continue;
0221 if ((*iter).mip())
0222 nmip++;
0223 if (!(*iter).quiet())
0224 nnq++;
0225
0226
0227
0228
0229
0230
0231 }
0232
0233 of << "MIP " << nmip << " ";
0234 for (iter = regions->begin(); iter != regions->end(); iter++) {
0235 if ((*iter).id().ieta() < 4 || (*iter).id().ieta() > 17 || (*iter).id().iphi() > 17)
0236 continue;
0237 if ((*iter).mip())
0238 of << (*iter).id().ieta() - 4 << " " << (((*iter).id().iphi() + 9) % 18) << " ";
0239 }
0240 of << endl;
0241
0242 of << "NQ " << nnq << " ";
0243 for (iter = regions->begin(); iter != regions->end(); iter++) {
0244 if ((*iter).id().ieta() < 4 || (*iter).id().ieta() > 17 || (*iter).id().iphi() > 17)
0245 continue;
0246 if (!(*iter).quiet())
0247 of << (*iter).id().ieta() - 4 << " " << (((*iter).id().iphi() + 9) % 18) << " ";
0248 }
0249 of << endl;
0250 }
0251 }
0252
0253 void L1MuGMTPattern::printMI(const vector<unsigned>* mi) {
0254 if (m_outputType == 1) {
0255 ofstream of(m_outfilename.c_str(), ios::app);
0256 of << "MI ";
0257 vector<unsigned>::const_iterator imi;
0258 for (imi = mi->begin(); imi != mi->end(); imi++) {
0259 of << " " << *imi;
0260 }
0261 of << endl;
0262 }
0263 }
0264
0265 void L1MuGMTPattern::printCANC() {
0266 if (m_outputType == 1) {
0267 ofstream of(m_outfilename.c_str(), ios::app);
0268 of << "CANC 0 0 0 0" << endl;
0269 }
0270 }
0271
0272 unsigned L1MuGMTPattern::invertQPt(unsigned w) {
0273 unsigned qpt = (w >> 8) & 0xff;
0274 w &= 0xffff00ff;
0275 w |= ((~qpt) << 8) & 0x0000ff00;
0276 return w;
0277 }
0278
0279
0280 DEFINE_FWK_MODULE(L1MuGMTPattern);
