|
|
|
#include <L1Trigger/TextToDigi/interface/RctDigiToSourceCardText.h>
Public Member Functions | |
| RctDigiToSourceCardText (const edm::ParameterSet &) | |
| ~RctDigiToSourceCardText () | |
Private Member Functions | |
| virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
Private Attributes | |
| std::ofstream | m_file |
| file handle | |
| unsigned short | m_nevt |
| event counter | |
| edm::InputTag | m_rctInputLabel |
| Label for RCT digis. | |
| SourceCardRouting | m_scRouting |
| source card router | |
| std::string | m_textFileName |
| Name out output file. | |
Description: Input RCT digis and output GCT text file to be loaded into the source cards for pattern tests.
Description: Input RCT digis and output text file to be loaded into the source cards for pattern tests.
Implementation: <Notes on="" implementation>="">
Definition at line 44 of file RctDigiToSourceCardText.h.
| RctDigiToSourceCardText::RctDigiToSourceCardText | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 28 of file RctDigiToSourceCardText.cc.
References Exception, m_file, m_scRouting, m_textFileName, and dbtoconf::out.
: m_rctInputLabel(iConfig.getParameter<edm::InputTag>("RctInputLabel")), m_textFileName(iConfig.getParameter<std::string>("TextFileName")), m_nevt(0) { // Open the output file m_file.open(m_textFileName.c_str(),std::ios::out); if(!m_file.good()) { throw cms::Exception("RctDigiToSourceCardTextFileOpenError") << "RctDigiToSourceCardText::RctDigiToSourceCardText : " << " couldn't open the file " << m_textFileName << " for writing" << std::endl; } // Make a SC routing object SourceCardRouting m_scRouting; }
| RctDigiToSourceCardText::~RctDigiToSourceCardText | ( | ) |
Definition at line 48 of file RctDigiToSourceCardText.cc.
References m_file.
{
// Close the output file
m_file.close();
}
| void RctDigiToSourceCardText::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 55 of file RctDigiToSourceCardText.cc.
References SourceCardRouting::EMUtoSTRING(), edm::Event::getByLabel(), LogDebug, m_file, m_nevt, m_rctInputLabel, m_scRouting, NUM_RCT_CRATES, SourceCardRouting::RC012toSTRING(), SourceCardRouting::RC234toSTRING(), SourceCardRouting::RC56HFtoSTRING(), and SourceCardRouting::RoutingModetoLogicalCardID().
{
// get the RCT data
Handle<L1CaloEmCollection> em;
Handle<L1CaloRegionCollection> rgn;
iEvent.getByLabel(m_rctInputLabel, em);
iEvent.getByLabel(m_rctInputLabel, rgn);
// General variables
int RoutingMode;
unsigned short logicalCardID;
std::string dataString;
// Arrays to hold electron variables
unsigned short eIsoRank[18][4]={{0}};
unsigned short eIsoCardId[18][4]={{0}};
unsigned short eIsoRegionId[18][4]={{0}};
unsigned short eNonIsoRank[18][4]={{0}};
unsigned short eNonIsoCardId[18][4]={{0}};
unsigned short eNonIsoRegionId[18][4]={{0}};
// Fill electrons
unsigned numIsoEM[18]={0};
unsigned numNonIsoEM[18]={0};
for (L1CaloEmCollection::const_iterator iem=em->begin(); iem!=em->end(); iem++){
if (iem->isolated()){
eIsoRank[iem->rctCrate()][numIsoEM[iem->rctCrate()]]=iem->rank();
eIsoCardId[iem->rctCrate()][numIsoEM[iem->rctCrate()]]=iem->rctCard();
eIsoRegionId[iem->rctCrate()][numIsoEM[iem->rctCrate()]]=iem->rctRegion();
numIsoEM[iem->rctCrate()]++;
} else {
eNonIsoRank[iem->rctCrate()][numNonIsoEM[iem->rctCrate()]]=iem->rank();
eNonIsoCardId[iem->rctCrate()][numNonIsoEM[iem->rctCrate()]]=iem->rctCard();
eNonIsoRegionId[iem->rctCrate()][numNonIsoEM[iem->rctCrate()]]=iem->rctRegion();
numNonIsoEM[iem->rctCrate()]++;
}
// Debug info
LogDebug("Electrons") << "Rank=" << iem->rank()
<< " Card=" << iem->rctCard()
<< " Region=" << iem->rctRegion()
<< " Crate=" << iem->rctCrate()
<< " Isolated=" << iem->isolated();
}
// Arrays to hold region variables
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}}};
unsigned short HFQ[18][4][2]={{{0}}};
unsigned short MIPbits[18][7][2]={{{0}}};
unsigned short Qbits[18][7][2]={{{0}}};
// Fill regions
for (L1CaloRegionCollection::const_iterator irgn=rgn->begin(); irgn!=rgn->end(); irgn++){
if (irgn->id().isHf()){
HF[irgn->rctCrate()][irgn->id().rctEta()-7][irgn->id().rctPhi()]=irgn->et();
HFQ[irgn->rctCrate()][irgn->id().rctEta()-7][irgn->id().rctPhi()]=irgn->fineGrain();
// Debug info
LogDebug("HFRegions") << "Et=" << irgn->et()
<< " FineGrain=" << irgn->fineGrain()
<< " Eta=" << irgn->id().rctEta()
<< " Phi=" << irgn->id().rctPhi()
<< " Crate=" << irgn->rctCrate();
} else {
RC[irgn->rctCrate()][irgn->rctCard()][irgn->rctRegionIndex()]=irgn->et();
RCof[irgn->rctCrate()][irgn->rctCard()][irgn->rctRegionIndex()]=irgn->overFlow();
RCtau[irgn->rctCrate()][irgn->rctCard()][irgn->rctRegionIndex()]=irgn->tauVeto();
MIPbits[irgn->rctCrate()][irgn->rctCard()][irgn->rctRegionIndex()]=irgn->mip();
Qbits[irgn->rctCrate()][irgn->rctCard()][irgn->rctRegionIndex()]=irgn->quiet();
// Debug info
LogDebug("Regions") << "Et=" << irgn->et()
<< " OverFlow=" << irgn->overFlow()
<< " tauVeto=" << irgn->tauVeto()
<< " mip=" << irgn->mip()
<< " quiet=" << irgn->quiet()
<< " Card=" << irgn->rctCard()
<< " Region=" << irgn->rctRegionIndex()
<< " Crate=" << irgn->rctCrate();
}
}
for (int crate=0; crate<NUM_RCT_CRATES; crate++){
// Logical Card ID = Source Card number
RoutingMode=0;
m_scRouting.RoutingModetoLogicalCardID(logicalCardID,RoutingMode,crate);
// Convert electrons to SC format
m_scRouting.EMUtoSTRING(logicalCardID,
m_nevt,
eIsoRank[crate],
eIsoCardId[crate],
eIsoRegionId[crate],
eNonIsoRank[crate],
eNonIsoCardId[crate],
eNonIsoRegionId[crate],
MIPbits[crate],
Qbits[crate],
dataString);
// Write electrons
m_file << dataString;
// Logical Card ID = Source Card number
RoutingMode=1;
m_scRouting.RoutingModetoLogicalCardID(logicalCardID,RoutingMode,crate);
// Convert regions to SC format
m_scRouting.RC56HFtoSTRING(logicalCardID,
m_nevt,
RC[crate],
RCof[crate],
RCtau[crate],
HF[crate],
HFQ[crate],
dataString);
// Write regions
m_file << dataString;
// Logical Card ID = Source Card number
RoutingMode=2;
m_scRouting.RoutingModetoLogicalCardID(logicalCardID,RoutingMode,crate);
// Convert regions to SC format
m_scRouting.RC012toSTRING(logicalCardID,
m_nevt,
RC[crate],
RCof[crate],
RCtau[crate],
dataString);
// Write regions
m_file << dataString;
// This is to 9 only as this is the shared source card
if (crate<9){
// Logical Card ID = Source Card number
RoutingMode=3;
m_scRouting.RoutingModetoLogicalCardID(logicalCardID,RoutingMode,crate);
// Convert regions to SC format
m_scRouting.RC234toSTRING(logicalCardID,
m_nevt,
RC[crate],
RCof[crate],
RCtau[crate],
RC[crate+9],
RCof[crate+9],
RCtau[crate+9],
dataString);
// Write regions
m_file << dataString;
}
}
// Force write to file
m_file << flush;
m_nevt++;
}
std::ofstream RctDigiToSourceCardText::m_file [private] |
file handle
Definition at line 59 of file RctDigiToSourceCardText.h.
Referenced by analyze(), RctDigiToSourceCardText(), and ~RctDigiToSourceCardText().
unsigned short RctDigiToSourceCardText::m_nevt [private] |
Label for RCT digis.
Definition at line 53 of file RctDigiToSourceCardText.h.
Referenced by analyze().
source card router
Definition at line 65 of file RctDigiToSourceCardText.h.
Referenced by analyze(), and RctDigiToSourceCardText().
std::string RctDigiToSourceCardText::m_textFileName [private] |
Name out output file.
Definition at line 56 of file RctDigiToSourceCardText.h.
Referenced by RctDigiToSourceCardText().
1.7.3