|
|
|
#include <CSCTFUnpacker.h>
Public Member Functions | |
| CSCTFUnpacker (const edm::ParameterSet &pset) | |
| void | produce (edm::Event &e, const edm::EventSetup &c) |
| ~CSCTFUnpacker (void) | |
Private Attributes | |
| int | m_maxBX |
| int | m_minBX |
| CSCTriggerMapping * | mapping |
| edm::InputTag | producer |
| std::vector< int > | slot2sector |
| bool | swapME1strips |
| CSCTFEvent | tfEvent |
Definition at line 16 of file CSCTFUnpacker.h.
| CSCTFUnpacker::CSCTFUnpacker | ( | const edm::ParameterSet & | pset | ) |
Definition at line 39 of file CSCTFUnpacker.cc.
References CSCTriggerMapping::addRecord(), Reference_intrackfit_cff::endcap, Exception, edm::ParameterSet::getParameter(), LogDebug, m_maxBX, m_minBX, mapping, producer, slot2sector, relativeConstraints::station, and swapME1strips.
:edm::EDProducer(),mapping(0){ LogDebug("CSCTFUnpacker|ctor")<<"Started ..."; // Edges of the time window, which LCTs are put into (unlike tracks, which are always centred around 0): m_minBX = pset.getParameter<int>("MinBX"); //3 m_maxBX = pset.getParameter<int>("MaxBX"); //9 // Swap: if(swapME1strips && me1b && !zplus) strip = 65 - strip; // 1-64 -> 64-1 : swapME1strips = pset.getParameter<bool>("swapME1strips"); // Initialize slot<->sector assignment slot2sector = pset.getParameter< std::vector<int> >("slot2sector"); LogDebug("CSCTFUnpacker|ctor")<<"Verifying slot<->sector map from 'vint32 slot2sector'"; for(int slot=0; slot<22; slot++) if( slot2sector[slot]<0 || slot2sector[slot]>12 ) throw cms::Exception("Invalid configuration")<<"CSCTFUnpacker: sector index is set out of range (slot2sector["<<slot<<"]="<<slot2sector[slot]<<", should be [0-12])"; // Just for safety (in case of bad data): slot2sector.resize(32); // As we use standard CSC digi containers, we have to initialize mapping: std::string mappingFile = pset.getParameter<std::string>("mappingFile"); if( mappingFile.length() ){ LogDebug("CSCTFUnpacker|ctor") << "Define ``mapping'' only if you want to screw up real geometry"; mapping = new CSCTriggerMappingFromFile(mappingFile); } else { LogDebug("CSCTFUnpacker|ctor") << "Generating default hw<->geometry mapping"; class M: public CSCTriggerSimpleMapping{ void fill(void){} }; mapping = new M(); for(int endcap=1; endcap<=2; endcap++) for(int station=1; station<=4; station++) for(int sector=1; sector<=6; sector++) for(int csc=1; csc<=9; csc++){ if( station==1 ){ mapping->addRecord(endcap,station,sector,1,csc,endcap,station,sector,1,csc); mapping->addRecord(endcap,station,sector,2,csc,endcap,station,sector,2,csc); } else mapping->addRecord(endcap,station,sector,0,csc,endcap,station,sector,0,csc); } } producer = pset.getParameter<edm::InputTag>("producer"); produces<CSCCorrelatedLCTDigiCollection>(); produces<L1CSCTrackCollection>(); produces<L1CSCStatusDigiCollection>(); produces<CSCTriggerContainer<csctf::TrackStub> >("DT"); }
| CSCTFUnpacker::~CSCTFUnpacker | ( | void | ) |
Definition at line 87 of file CSCTFUnpacker.cc.
References mapping.
| void CSCTFUnpacker::produce | ( | edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [virtual] |
Implements edm::EDProducer.
Definition at line 91 of file CSCTFUnpacker.cc.
References L1CSCSPStatusDigi::af, CSCSPRecord::AFs(), newFWLiteAna::base, L1CSCSPStatusDigi::bx, CSCSPHeader::BXN(), CSCSPRecord::BXs(), CSCSPEvent::counters(), FEDRawData::data(), CSCTriggerMapping::detId(), CSCSPHeader::endcap(), Reference_intrackfit_cff::endcap, L1CSCSPStatusDigi::fmm_status, edm::Event::getByLabel(), CSCSPEvent::header(), edm::InputTag::instance(), CSCSPHeader::L1A(), L1CSCSPStatusDigi::l1a_bxn, edm::InputTag::label(), CSCSPRecord::LCT(), link(), L1CSCSPStatusDigi::link_status, m_maxBX, m_minBX, mapping, FEDNumbering::MAXCSCTFFEDID, CSCSPRecord::mbStubs(), FEDNumbering::MINCSCTFFEDID, L1CSCSPStatusDigi::mpc_link_id, CSCSPHeader::nTBINs(), L1CSCSPStatusDigi::orbit_cnt, CSCSPCounters::orbit_counter(), listBenchmarks::pattern, producer, CSCSPRecord::ptSpy(), edm::Event::put(), CSCSPEvent::record(), L1CSCSPStatusDigi::se, CSCSPHeader::sector(), CSCSPRecord::SEs(), FEDRawData::size(), CSCSPHeader::slot(), slot2sector, L1CSCSPStatusDigi::sm, CSCSPRecord::SMs(), L1CSCSPStatusDigi::sp_slot, CSCTFEvent::SPs_fast(), relativeConstraints::station, CSCSPHeader::status(), ntuplemaker::status, strip(), swapME1strips, tfEvent, L1CSCSPStatusDigi::track_cnt, CSCSPCounters::track_counter(), CSCSPRecord::tracks(), testEve_cfg::tracks, CSCTFEvent::unpack(), L1CSCSPStatusDigi::vp, CSCSPRecord::VPs(), and cms::Exception::what().
{
// Get a handle to the FED data collection
edm::Handle<FEDRawDataCollection> rawdata;
e.getByLabel(producer.label(),producer.instance(),rawdata);
// create the collection of CSC wire and strip digis as well as of DT stubs, which we receive from DTTF
std::auto_ptr<CSCCorrelatedLCTDigiCollection> LCTProduct(new CSCCorrelatedLCTDigiCollection);
std::auto_ptr<L1CSCTrackCollection> trackProduct(new L1CSCTrackCollection);
std::auto_ptr<L1CSCStatusDigiCollection> statusProduct(new L1CSCStatusDigiCollection);
std::auto_ptr<CSCTriggerContainer<csctf::TrackStub> > dtProduct(new CSCTriggerContainer<csctf::TrackStub>);
for(int fedid=FEDNumbering::MINCSCTFFEDID; fedid<=FEDNumbering::MAXCSCTFFEDID; fedid++){
const FEDRawData& fedData = rawdata->FEDData(fedid);
if( fedData.size()==0 ) continue;
//LogDebug("CSCTFUnpacker|produce");
//if( monitor ) monitor->process((unsigned short*)fedData.data());
unsigned int unpacking_status = tfEvent.unpack((unsigned short*)fedData.data(),fedData.size()/2);
if( unpacking_status==0 ){
// There may be several SPs in event
std::vector<const CSCSPEvent*> SPs = tfEvent.SPs_fast();
// Cycle over all of them
for(std::vector<const CSCSPEvent *>::const_iterator spItr=SPs.begin(); spItr!=SPs.end(); spItr++){
const CSCSPEvent *sp = *spItr;
L1CSCSPStatusDigi status;
status.sp_slot = sp->header().slot();
status.l1a_bxn = sp->header().BXN();
status.fmm_status = sp->header().status();
status.track_cnt = sp->counters().track_counter();
status.orbit_cnt = sp->counters().orbit_counter();
// Finds central LCT BX
// assumes window is odd number of bins
int central_lct_bx = (m_maxBX + m_minBX)/2;
// Find central SP BX
// assumes window is odd number of bins
int central_sp_bx = int(sp->header().nTBINs()/2);
for(unsigned int tbin=0; tbin<sp->header().nTBINs(); tbin++){
status.se |= sp->record(tbin).SEs();
status.sm |= sp->record(tbin).SMs();
status.bx |= sp->record(tbin).BXs();
status.af |= sp->record(tbin).AFs();
status.vp |= sp->record(tbin).VPs();
for(unsigned int FPGA=0; FPGA<5; FPGA++)
for(unsigned int MPClink=0; MPClink<3; ++MPClink){
std::vector<CSCSP_MEblock> lct = sp->record(tbin).LCT(FPGA,MPClink);
if( lct.size()==0 ) continue;
status.link_status[lct[0].spInput()] |=
(1<<lct[0].receiver_status_frame1())|
(1<<lct[0].receiver_status_frame2())|
((lct[0].aligment_fifo()?1:0)<<4);
status.mpc_link_id |= (lct[0].link()<<2)|lct[0].mpc();
int station = ( FPGA ? FPGA : 1 );
int endcap=0, sector=0;
if( slot2sector[sp->header().slot()] ){
endcap = slot2sector[sp->header().slot()]/7 + 1;
sector = slot2sector[sp->header().slot()];
if( sector>6 ) sector -= 6;
} else {
endcap = (sp->header().endcap()?1:2);
sector = sp->header().sector();
}
int subsector = ( FPGA>1 ? 0 : FPGA+1 );
int cscid = lct[0].csc() ;
try{
CSCDetId id = mapping->detId(endcap,station,sector,subsector,cscid,0);
// corrlcts now have no layer associated with them
LCTProduct->insertDigi(id,
CSCCorrelatedLCTDigi(
0,lct[0].vp(),lct[0].quality(),lct[0].wireGroup(),
(swapME1strips && cscid<=3 && station==1 && endcap==2 && lct[0].strip()<65 ? 65 - lct[0].strip() : lct[0].strip() ),
lct[0].pattern(),lct[0].l_r(),
(lct[0].tbin()+(central_lct_bx-central_sp_bx)),
lct[0].link(), lct[0].BXN(), 0, cscid )
);
} catch(cms::Exception &e) {
edm::LogInfo("CSCTFUnpacker|produce") << e.what() << "Not adding digi to collection in event "
<<sp->header().L1A()<<" (endcap="<<endcap<<",station="<<station<<",sector="<<sector<<",subsector="<<subsector<<",cscid="<<cscid<<",spSlot="<<sp->header().slot()<<")";
}
}
std::vector<CSCSP_MBblock> mbStubs = sp->record(tbin).mbStubs();
for(std::vector<CSCSP_MBblock>::const_iterator iter=mbStubs.begin(); iter!=mbStubs.end(); iter++){
int endcap, sector;
if( slot2sector[sp->header().slot()] ){
endcap = slot2sector[sp->header().slot()]/7 + 1;
sector = slot2sector[sp->header().slot()];
if( sector>6 ) sector -= 6;
} else {
endcap = (sp->header().endcap()?1:2);
sector = sp->header().sector();
}
const unsigned int csc2dt[6][2] = {{2,3},{4,5},{6,7},{8,9},{10,11},{12,1}};
DTChamberId id((endcap==1?2:-2),1, csc2dt[sector-1][iter->id()-1]);
CSCCorrelatedLCTDigi base(0,iter->vq(),iter->quality(),iter->cal(),iter->flag(),iter->bc0(),iter->phi_bend(),tbin+(central_lct_bx-central_sp_bx),iter->id(),iter->bxn(),iter->timingError(),iter->BXN());
csctf::TrackStub dtStub(base,id,iter->phi(),0);
dtProduct->push_back(dtStub);
}
std::vector<CSCSP_SPblock> tracks = sp->record(tbin).tracks();
unsigned int trkNumber=0;
for(std::vector<CSCSP_SPblock>::const_iterator iter=tracks.begin(); iter!=tracks.end(); iter++,trkNumber++){
L1CSCTrack track;
if( slot2sector[sp->header().slot()] ){
track.first.m_endcap = slot2sector[sp->header().slot()]/7 + 1;
track.first.m_sector = slot2sector[sp->header().slot()];
if( track.first.m_sector>6 ) track.first.m_sector -= 6;
} else {
track.first.m_endcap = (sp->header().endcap()?1:2);
track.first.m_sector = sp->header().sector();
}
track.first.m_lphi = iter->phi();
track.first.m_ptAddress = iter->ptLUTaddress();
track.first.m_fr = iter->f_r();
track.first.m_ptAddress|=(iter->f_r() << 21);
track.first.setStationIds(iter->ME1_id(),iter->ME2_id(),iter->ME3_id(),iter->ME4_id(),iter->MB_id());
track.first.setTbins(iter->ME1_tbin(), iter->ME2_tbin(), iter->ME3_tbin(), iter->ME4_tbin(), iter->MB_tbin() );
track.first.setBx(iter->tbin()-central_sp_bx);
track.first.setBits(iter->syncErr(), iter->bx0(), iter->bc0());
track.first.setLocalPhi(iter->phi());
track.first.setEtaPacked(iter->eta());
track.first.setChargePacked(iter->charge());
track.first.m_output_link = iter->id();
if( track.first.m_output_link ){
track.first.m_rank = (iter->f_r()?sp->record(tbin).ptSpy()&0x7F:(sp->record(tbin).ptSpy()&0x7F00)>>8);
track.first.setChargeValidPacked((iter->f_r()?(sp->record(tbin).ptSpy()&0x80)>>8:(sp->record(tbin).ptSpy()&0x8000)>>15));
} else {
track.first.m_rank = 0;
track.first.setChargeValidPacked(0);
}
track.first.setFineHaloPacked(iter->halo());
track.first.m_winner = iter->MS_id()&(1<<trkNumber);
std::vector<CSCSP_MEblock> lcts = iter->LCTs();
for(std::vector<CSCSP_MEblock>::const_iterator lct=lcts.begin(); lct!=lcts.end(); lct++){
int station = ( lct->spInput()>6 ? (lct->spInput()-1)/3 : 1 );
int subsector = ( lct->spInput()>6 ? 0 : (lct->spInput()-1)/3 + 1 );
try{
CSCDetId id = mapping->detId(track.first.m_endcap,station,track.first.m_sector,subsector,lct->csc(),0);
track.second.insertDigi(id,
CSCCorrelatedLCTDigi(
0,lct->vp(),lct->quality(),lct->wireGroup(),
(swapME1strips && lct->csc()<=3 && station==1 && track.first.m_endcap==2 && lct[0].strip()<65 ? 65 - lct[0].strip() : lct[0].strip() ),
lct->pattern(),lct->l_r(),
(lct->tbin()+(central_lct_bx-central_sp_bx)),
lct->link(), lct->BXN(), 0, lct->csc() )
);
} catch(cms::Exception &e) {
edm::LogInfo("CSCTFUnpacker|produce") << e.what() << "Not adding track digi to collection in event"
<<sp->header().L1A()<<" (endcap="<<track.first.m_endcap<<",station="<<station<<",sector="<<track.first.m_sector<<",subsector="<<subsector<<",cscid="<<lct->csc()<<",spSlot="<<sp->header().slot()<<")";
}
}
std::vector<CSCSP_MBblock> mbStubs = iter->dtStub();
for(std::vector<CSCSP_MBblock>::const_iterator iter=mbStubs.begin(); iter!=mbStubs.end(); iter++){
CSCDetId id = mapping->detId(track.first.m_endcap,1,track.first.m_sector,iter->id(),1,0);
track.second.insertDigi(id,
CSCCorrelatedLCTDigi(iter->phi(),iter->vq(),iter->quality()+100,iter->cal(),iter->flag(),iter->bc0(),iter->phi_bend(),tbin+(central_lct_bx-central_sp_bx),iter->id(),iter->bxn(),iter->timingError(),iter->BXN())
);
}
trackProduct->push_back( track );
}
}
statusProduct->second.push_back( status );
}
} else {
edm::LogError("CSCTFUnpacker|produce")<<" problem of unpacking TF event: 0x"<<std::hex<<unpacking_status<<std::dec<<" code";
}
statusProduct->first = unpacking_status;
} //end of fed cycle
e.put(dtProduct,"DT");
e.put(LCTProduct); // put processed lcts into the event.
e.put(trackProduct);
e.put(statusProduct);
}
int CSCTFUnpacker::m_maxBX [private] |
Definition at line 18 of file CSCTFUnpacker.h.
Referenced by CSCTFUnpacker(), and produce().
int CSCTFUnpacker::m_minBX [private] |
Definition at line 18 of file CSCTFUnpacker.h.
Referenced by CSCTFUnpacker(), and produce().
CSCTriggerMapping* CSCTFUnpacker::mapping [private] |
Definition at line 21 of file CSCTFUnpacker.h.
Referenced by CSCTFUnpacker(), produce(), and ~CSCTFUnpacker().
edm::InputTag CSCTFUnpacker::producer [private] |
Definition at line 30 of file CSCTFUnpacker.h.
Referenced by CSCTFUnpacker(), and produce().
std::vector<int> CSCTFUnpacker::slot2sector [private] |
Definition at line 27 of file CSCTFUnpacker.h.
Referenced by CSCTFUnpacker(), and produce().
bool CSCTFUnpacker::swapME1strips [private] |
Definition at line 19 of file CSCTFUnpacker.h.
Referenced by CSCTFUnpacker(), and produce().
CSCTFEvent CSCTFUnpacker::tfEvent [private] |
Definition at line 23 of file CSCTFUnpacker.h.
Referenced by produce().
1.7.3