|
|
|
#include <CSCTFDTReceiver.h>
Public Member Functions | |
| CSCTFDTReceiver () | |
| CSCTriggerContainer < csctf::TrackStub > | process (const L1MuDTChambPhContainer *) |
| ~CSCTFDTReceiver () | |
Private Attributes | |
| CSCTriggerContainer < csctf::TrackStub > | dtstubs |
Definition at line 12 of file CSCTFDTReceiver.h.
| CSCTFDTReceiver::CSCTFDTReceiver | ( | ) | [inline] |
Definition at line 16 of file CSCTFDTReceiver.h.
References CSCTriggerContainer< T >::clear(), and dtstubs.
| CSCTFDTReceiver::~CSCTFDTReceiver | ( | ) | [inline] |
Definition at line 17 of file CSCTFDTReceiver.h.
{}
| CSCTriggerContainer< csctf::TrackStub > CSCTFDTReceiver::process | ( | const L1MuDTChambPhContainer * | dttrig | ) |
shift all by one and take mod 8, since DT quality of 7 is a null stub
Definition at line 12 of file CSCTFDTReceiver.cc.
References L1MuDTChambPhContainer::chPhiSegm1(), L1MuDTChambPhContainer::chPhiSegm2(), CSCTriggerContainer< T >::clear(), L1MuDTChambPhDigi::code(), gather_cfg::cout, dtstubs, ExpressReco_HICollisions_FallBack::e, L1MuDTTFConfig::getBxMax(), L1MuDTTFConfig::getBxMin(), CSCBitWidths::kGlobalPhiDataBitWidth, CSCTFDTReceiverLUT::lut, CSCDetId::maxEndcapId(), CSCTriggerNumbering::maxTriggerSectorId(), CSCDetId::minEndcapId(), CSCTriggerNumbering::minTriggerSectorId(), phi, L1MuDTChambPhDigi::phiB(), CSCTriggerContainer< T >::push_back(), DTConfigTraco::RESOLPSI, and asciidump::s.
Referenced by CSCTFTrackProducer::produce().
{
dtstubs.clear();
if( !dttrig ) return dtstubs;
const int dt_minBX = L1MuDTTFConfig::getBxMin();
const int dt_maxBX = L1MuDTTFConfig::getBxMax();
const int dt_toffs = 0;// changed since DT tpg now centers around zero //abs(dt_maxBX - dt_minBX);
// consider all BX
for(int bx = dt_minBX + dt_toffs; bx <= dt_maxBX + dt_toffs; ++bx)
for(int e = CSCDetId::minEndcapId(); e <= CSCDetId::maxEndcapId(); ++e)
for(int s = CSCTriggerNumbering::minTriggerSectorId(); s <= CSCTriggerNumbering::maxTriggerSectorId(); ++s)
{
int wheel = (e == 1) ? 2 : -2;
int sector = 2*s - 1;
int csc_bx = bx + 6;//Delay DT stubs by 6 bx.
// combine two 30 degree DT sectors into a 60 degree CSC
// sector.
for(int is = sector; is <= sector+1; ++is)
{
int iss = (is == 12) ? 0 : is;
L1MuDTChambPhDigi* dtts[2];
for(int stub = 0; stub < 2; ++stub)
{
dtts[stub] = (stub == 0) ? dttrig->chPhiSegm1(wheel,1,iss,bx) :
dttrig->chPhiSegm2(wheel,1,iss,bx);
if(dtts[stub])
{
// --------------------------------------------------------------
// IKF: this code has been reformulated ...
// --------------------------------------------------------------
// // Convert stubs to CSC format (signed -> unsigned)
// // phi was 12 bits (signed) for pi radians = 57.3 deg
// // relative to center of 30 degree DT sector
// double tmp = static_cast<const double> (dtts[stub]->phi()) /
// DTConfigTraco::RESOLPSIR * 180./M_PI + 15.;
// int phi = static_cast<int> (tmp/62. * (1<<(CSCBitWidths::kGlobalPhiDataBitWidth)));
// --------------------------------------------------------------
// IKF ...and is now this line, actually works a tiny bit better.
// --------------------------------------------------------------
// float tmp = dtts[stub] -> phi() * 1.0;
//
// tmp *= 90.0;
// tmp /= 31.0;
// // tmp /= M_PI;
// tmp /= 3.1416;
// tmp += 1057.0;
//
// int phi = static_cast<int> (tmp);
// --------------------------------------------------------------
// IKF ...and is now this line, actually works a tiny bit better.
// --------------------------------------------------------------
int phi = dtts[stub] -> phi();
if (phi < 0) phi += 4096;
if (phi > 4096)
{ std::cout << "AAAAAAAAAAGH TOO BIG PHI:" << phi << std::endl;
continue;
}
if (phi < 0){
std::cout << "AAAAAAAAH NEG PHI" << phi << std::endl;
continue;
}
phi = CSCTFDTReceiverLUT::lut[phi];
// --------------------------------------------------------------
// DT chambers may lie outside CSC sector boundary
// Eventually we need to extend CSC phi definition
// --------------------------------------------------------------
// IKF: this is a protection, physically can't happen in data (bus too narrow) -
// - what really happens in data?
// --------------------------------------------------------------
phi = (phi>0) ? phi : 0;
phi = (phi<(1<<(CSCBitWidths::kGlobalPhiDataBitWidth))) ? phi :
(1<<(CSCBitWidths::kGlobalPhiDataBitWidth))-1;
// change phib from 10 bits to 6
int phib = (dtts[stub]->phiB() + DTConfigTraco::RESOLPSI) / 16;
int qual = dtts[stub]->code();
// barrel allows quality=0!
qual = (qual + 1)%8;
CSCCorrelatedLCTDigi dtinfo(stub+1,1, qual, 0, stub, 0, phib, csc_bx+stub, 1+(is+1)%2);
DTChamberId dtid(wheel,1,iss+1);
csctf::TrackStub tsCSC(dtinfo,dtid, phi, 0);
dtstubs.push_back(tsCSC);
}
}
}
}
return dtstubs;
}
Definition at line 24 of file CSCTFDTReceiver.h.
Referenced by CSCTFDTReceiver(), and process().
1.7.3