|
|
|
#include <RctDigiToRctText.h>
Public Member Functions | |
| RctDigiToRctText (const edm::ParameterSet &) | |
| ~RctDigiToRctText () | |
Private Member Functions | |
| virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
Private Attributes | |
| std::ofstream | fdebug |
| handle for debug file | |
| std::ofstream | m_file [NUM_RCT_CRATES] |
| handles for output files | |
| bool | m_hexUpperCase |
| write upper case hex words | |
| edm::InputTag | m_rctInputLabel |
| label for RCT digis | |
| std::string | m_textFileName |
| basename for output files | |
Definition at line 32 of file RctDigiToRctText.h.
| RctDigiToRctText::RctDigiToRctText | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
open output text files
open info|debug file
Definition at line 12 of file RctDigiToRctText.cc.
References Exception, fdebug, convertXMLtoSQLite_cfg::fileName, i, m_file, m_textFileName, NUM_RCT_CRATES, and dbtoconf::out.
: m_rctInputLabel(iConfig.getParameter<edm::InputTag>("RctInputLabel")), m_textFileName(iConfig.getParameter<std::string>("TextFileName")), m_hexUpperCase(iConfig.getParameter<bool>("HexUpperCase")) { for (unsigned i=0; i<NUM_RCT_CRATES; i++){ std::stringstream fileStream; fileStream << m_textFileName << std::setw(2) << std::setfill('0') << i << ".txt"; std::string fileName(fileStream.str()); m_file[i].open(fileName.c_str(),std::ios::out); if(!m_file[i].good()) { throw cms::Exception("RctDigiToRctTextTextFileOpenError") << "RctDigiToRctText::RctDigiToRctText : " << " couldn't create the file " << fileName << std::endl; } } fdebug.open("rctdigitorcttext_debug.txt", std::ios::out); }
| RctDigiToRctText::~RctDigiToRctText | ( | ) |
close files
Definition at line 37 of file RctDigiToRctText.cc.
References fdebug, i, m_file, and NUM_RCT_CRATES.
| void RctDigiToRctText::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
count bunch crossing
get the RCT data
debug flags and stream
--- ELECTRONS ---
debug
--- REGIONS ---
print electrons
print regions
mip bits
quiet bits
region info
HF
flush data to files
Implements edm::EDAnalyzer.
Definition at line 46 of file RctDigiToRctText.cc.
References data, Exception, fdebug, edm::Event::getByLabel(), i, errorMatrix2Lands_multiChannel::id, j, LogDebug, m_file, m_hexUpperCase, m_rctInputLabel, nevt, NUM_RCT_CRATES, metsig::tau, and tmp.
{
static int nevt = -1;
nevt++;
edm::Handle<L1CaloEmCollection> em;
edm::Handle<L1CaloRegionCollection> rgn;
iEvent.getByLabel(m_rctInputLabel, em);
iEvent.getByLabel(m_rctInputLabel, rgn);
bool ldebug = false;
bool debug_NOTEMPTY[18]={0};
for(int i=0; i<18; i++) debug_NOTEMPTY[i]=false;
std::stringstream dstrm;
unsigned long int CAND[18][8];
int n_iso [18]={0};
int n_niso[18]={0};
bool iso;
int id;
for(L1CaloEmCollection::const_iterator iem=em->begin(); iem!=em->end(); iem++){
int crate = iem->rctCrate();
iso = iem->isolated();
unsigned data = iem->raw();
id = iso ? n_iso[crate]++ : 4 + n_niso[crate]++;
CAND[crate][id] = data;
if(crate>17 || id > 7)
throw cms::Exception("RctDigiToRctTextElectronIndexOutBounds")
<< "out of bounds indices crate:" << crate << "id:" << id << std::endl;
if(ldebug&&data!=0)
debug_NOTEMPTY[crate] = true;
dstrm.str("");
dstrm << "electron "
<< " bx:" << nevt
<< " crate:" << crate
<< " iso:" << iso
<< " raw:" << data
<< " \t cand:" << *iem;
if(debug_NOTEMPTY[crate]) fdebug << dstrm.str() << std::endl;
}
unsigned short MIPbits [18][7][2] = {{{0}}};
unsigned short QIEbits [18][7][2] = {{{0}}};
unsigned short RC [18][7][2] = {{{0}}};
unsigned short RCof [18][7][2] = {{{0}}};
unsigned short RCtau [18][7][2] = {{{0}}};
unsigned short HF [18][4][2] = {{{0}}};
for (L1CaloRegionCollection::const_iterator irgn=rgn->begin(); irgn!=rgn->end(); irgn++){
int crate = irgn->rctCrate();
int card = irgn->rctCard();
int rgnidx= irgn->rctRegionIndex();
dstrm.str("");
if (!irgn->id().isHf()){
RC [crate][card][rgnidx]=irgn->et();
RCof [crate][card][rgnidx]=irgn->overFlow();
RCtau [crate][card][rgnidx]=irgn->tauVeto();
MIPbits[crate][card][rgnidx]=irgn->mip();
QIEbits[crate][card][rgnidx]=irgn->quiet();
// debug info
dstrm << hex
<< "Et=" << irgn->et()
<< " OverFlow=" << irgn->overFlow()
<< " tauVeto=" << irgn->tauVeto()
<< " mip=" << irgn->mip()
<< " quiet=" << irgn->quiet()
<< " Card=" << irgn->rctCard()
<< " Region=" << irgn->rctRegionIndex()
<< " Crate=" << irgn->rctCrate()
<< dec;
if(ldebug) LogDebug("Regions") << dstrm.str();
} else {
HF [crate][irgn->id().rctEta()-7][irgn->id().rctPhi()]=irgn->et();
// debug info
dstrm << hex
<< "Et=" << irgn->et()
<< " FGrain=" << irgn->fineGrain()
<< " Eta=" << irgn->id().rctEta()
<< " Phi=" << irgn->id().rctPhi()
<< " Crate=" << irgn->rctCrate()
<< dec;
if(ldebug) LogDebug("HFRegions") << dstrm.str();
}
if(ldebug &&irgn->et()!=0) debug_NOTEMPTY[crate] = true; //debug
if(debug_NOTEMPTY[crate]) {
fdebug << "region"
<< " bx:" << nevt
<< " crate:" << crate
<< "\t";
fdebug << dstrm.str() << std::endl;
}
}
std::stringstream sstrm;
if(m_hexUpperCase)
sstrm << std::uppercase;
else
sstrm.unsetf(std::ios::uppercase);
for (unsigned crate=0; crate<NUM_RCT_CRATES; crate++){
sstrm.str("");
sstrm << "Crossing " << nevt << std::endl;
for(int j=0; j<8; j++) {
sstrm << setw(3) << setfill('0') << hex << (CAND[crate][j] & 0x3ff);
if(j<7) sstrm << " ";
}
sstrm << setfill(' ') << dec;
m_file[crate] << sstrm.str() << std::endl;
// debug
if(debug_NOTEMPTY[crate]) fdebug << sstrm.str() << std::endl;
if(ldebug) LogDebug("Electrons") << sstrm.str() << std::endl;
}
for (unsigned crate=0; crate<NUM_RCT_CRATES; crate++){
sstrm.str("");
for (int card=0; card<7; card++){
for (int j=0; j<2; j++){
sstrm << " " << MIPbits[crate][card][j];
}
}
m_file[crate] << sstrm.str() << std::endl;
if(debug_NOTEMPTY[crate]) fdebug << sstrm.str() << std::endl; //debug
sstrm.str("");
for (int card=0; card<7; card++){
for (int j=0; j<2; j++){
sstrm << " " << QIEbits[crate][card][j];
}
}
m_file[crate] << sstrm.str() << std::endl;
if(debug_NOTEMPTY[crate]) fdebug << sstrm.str() << std::endl; //debug
sstrm.str("");
for (int card=0; card<7; card++){
for (int j=0; j<2; j++){
unsigned long int tmp;
unsigned et = RC [crate][card][j];
unsigned ovf = RCof [crate][card][j];
unsigned tau = RCtau[crate][card][j];
//ovf = ovf || (et>=0x400);
tmp = ((tau & 0x1) << 11) | ((ovf &0x1)<<10) | ((et & 0x3ff));
sstrm << " " << setw(3) << setfill('0') << hex << tmp;
}
}
m_file[crate] << sstrm.str() << std::endl;
if(debug_NOTEMPTY[crate]) fdebug << sstrm.str() << std::endl << std::endl; //debug
sstrm.str("");
for (int ip=0; ip<2; ip++){
for (int ie=0; ie<4; ie++){
unsigned et = HF[crate][ie][ip] & 0xff;
sstrm << " " << setw(2) << setfill('0') << hex << et;
}
}
m_file[crate] << sstrm.str() << std::endl;
if(debug_NOTEMPTY[crate]) fdebug << sstrm.str() << std::endl; //debug
sstrm << setfill(' ') << dec;
}// end crate loop
for (unsigned i=0; i<NUM_RCT_CRATES; i++)
m_file[i] << std::flush;
fdebug << std::flush;
}
std::ofstream RctDigiToRctText::fdebug [private] |
handle for debug file
Definition at line 55 of file RctDigiToRctText.h.
Referenced by analyze(), RctDigiToRctText(), and ~RctDigiToRctText().
std::ofstream RctDigiToRctText::m_file[NUM_RCT_CRATES] [private] |
handles for output files
Definition at line 52 of file RctDigiToRctText.h.
Referenced by analyze(), RctDigiToRctText(), and ~RctDigiToRctText().
bool RctDigiToRctText::m_hexUpperCase [private] |
write upper case hex words
Definition at line 49 of file RctDigiToRctText.h.
Referenced by analyze().
std::string RctDigiToRctText::m_textFileName [private] |
basename for output files
Definition at line 46 of file RctDigiToRctText.h.
Referenced by RctDigiToRctText().
1.7.3