IOPrinter.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//   Class: IOPrinter
//
//   IOPrinter
//
//
//   Author :
//   G. Flouris               U Ioannina    Feb. 2015
//   modifications:  G Karathanasis    U Athens
//--------------------------------------------------
#include "L1Trigger/L1TTwinMux/interface/IOPrinter.h"
#include "L1Trigger/L1TTwinMux/interface/L1MuTMChambPhContainer.h"
#include "DataFormats/Common/interface/Handle.h"

using namespace std;

void IOPrinter::run(edm::Handle<L1MuDTChambPhContainer> inphiDigis,
                    const L1MuDTChambPhContainer& outphiDigis,
                    edm::Handle<RPCDigiCollection> rpcDigis,
                    const RPCGeometry& rpcGeometry) {
  cout << "======================================================" << endl;
  int bx = 0, wheel = 0, sector = 0, station = 1;

  cout << "DT Inputs/RPCDT Inputs" << endl;
  cout << "bx\twheel\tsector\tstation\tphi\tphib\tcode\tts2tag\tbxcnt\trpcbit" << endl;

  for (bx = -2; bx <= 2; bx++) {
    for (wheel = -2; wheel <= 2; wheel++) {
      for (sector = 0; sector < 12; sector++) {
        for (station = 1; station <= 4; station++) {
          auto dtts1 = inphiDigis->chPhiSegm1(wheel, station, sector, bx);
          auto dtts2 = inphiDigis->chPhiSegm2(wheel, station, sector, bx - 1);
          if (dtts1 && dtts1->code() != 7) {
            L1MuDTChambPhDigi dt_ts1 = *dtts1;
            cout << dtts1->bxNum() << "\t" << dtts1->whNum() << "\t" << dtts1->scNum() << "\t" << dtts1->stNum() << "\t"
                 << dtts1->phi() << "\t" << dtts1->phiB() << "\t" << dtts1->code() << "\t" << dtts1->Ts2Tag() << "\t"
                 << dtts1->BxCnt() << "\t0" << endl;
          }
          if (dtts2 && dtts2->code() != 7) {
            L1MuDTChambPhDigi dt_ts2 = *dtts2;
            cout << dtts2->bxNum() << "\t" << dtts2->whNum() << "\t" << dtts2->scNum() << "\t" << dtts2->stNum() << "\t"
                 << dtts2->phi() << "\t" << dtts2->phiB() << "\t" << dtts2->code() << "\t" << dtts2->Ts2Tag() << "\t"
                 << dtts2->BxCnt() << "\t0" << endl;
          }
        }
      }
    }
  }

  cout << endl;

  cout << "RPC Inputs" << endl;
  cout << "bx\tring\tsector\tstation\troll\tlayer\tstrip\tphi\tlocalX" << endl;
  //cout<<"RPCHitCleaner"<<endl;
  for (auto chamber = rpcDigis->begin(); chamber != rpcDigis->end(); ++chamber) {
    RPCDetId detid = (*chamber).first;
    for (auto digi = (*chamber).second.first; digi != (*chamber).second.second; ++digi) {
      RPCDigi digi_out(digi->strip(), digi->bx());
      //if(digi->bx()!=0) continue;
      int phi = RPCtoDTTranslator::radialAngle(detid, rpcGeometry, digi->strip()) << 2;
      float localx = RPCtoDTTranslator::localX(detid, rpcGeometry, digi->strip());
      cout << digi->bx() << "\t" << detid.ring() << "\t" << detid.sector() - 1 << "\t" << detid.station() << "\t"
           << detid.roll() << "\t" << detid.layer() << "\t" << digi->strip() << "\t" << phi << "\t" << localx << endl;
    }  
  }    
  cout << endl;

  cout << "TwinMux Output" << endl;
  cout << "bx\twheel\tsector\tstation\tphi\tphib\tcode\tts2tag\tbxcnt\trpcbit" << endl;

  for (bx = -2; bx <= 2; bx++) {
    for (wheel = -2; wheel <= 2; wheel++) {
      for (sector = 0; sector < 12; sector++) {
        for (station = 1; station <= 4; station++) {
          auto dtts1 = outphiDigis.chPhiSegm1(wheel, station, sector, bx);
          auto dtts2 = outphiDigis.chPhiSegm2(wheel, station, sector, bx - 1);
          if (dtts1 && dtts1->code() != 7) {
            L1MuDTChambPhDigi dt_ts1 = *dtts1;
            cout << dtts1->bxNum() << "\t" << dtts1->whNum() << "\t" << dtts1->scNum() << "\t" << dtts1->stNum() << "\t"
                 << dtts1->phi() << "\t" << dtts1->phiB() << "\t" << dtts1->code() << "\t" << dtts1->Ts2Tag() << "\t"
                 << dtts1->BxCnt() << "\t" << dtts1->RpcBit() << endl;
          }
          if (dtts2 && dtts2->code() != 7) {
            L1MuDTChambPhDigi dt_ts2 = *dtts2;
            cout << dtts2->bxNum() << "\t" << dtts2->whNum() << "\t" << dtts2->scNum() << "\t" << dtts2->stNum() << "\t"
                 << dtts2->phi() << "\t" << dtts2->phiB() << "\t" << dtts2->code() << "\t" << dtts2->Ts2Tag() << "\t"
                 << dtts2->BxCnt() << "\t" << dtts2->RpcBit() << endl;
          }
        }
      }
    }
  }

  cout << "======================================================" << endl;
}

void IOPrinter::run(L1MuDTChambPhContainer const* inphiDigis,
                    const L1MuDTChambPhContainer& outphiDigis,
                    RPCDigiCollection const* rpcDigis,
                    const RPCGeometry& rpcGeometry) {
  cout << "======================================================" << endl;
  int bx = 0, wheel = 0, sector = 0, station = 1;
  L1MuTMChambPhContainer inphiDigis_tm;
  const std::vector<L1MuDTChambPhDigi>* vInCon = inphiDigis->getContainer();
  inphiDigis_tm.setContainer(*vInCon);

  cout << "RPC->DT Inputs" << endl;
  cout << "bx\twheel\tsector\tstation\tphi\tphib\tcode\tts2tag\tbxcnt\trpcbit" << endl;

  for (bx = -2; bx <= 2; bx++) {
    for (wheel = -2; wheel <= 2; wheel++) {
      for (sector = 0; sector < 12; sector++) {
        for (station = 1; station <= 4; station++) {
          int nhits = DTRPCBxCorrection::nRPCHits(*inphiDigis, bx, wheel, sector, station);
          for (int hit = 0; hit < nhits; hit++) {
            auto dtts1 = inphiDigis_tm.chPhiSegm(wheel, station, sector, bx, hit);
            if (dtts1) {
              cout << dtts1->bxNum() << "\t" << dtts1->whNum() << "\t" << dtts1->scNum() << "\t" << dtts1->stNum()
                   << "\t" << dtts1->phi() << "\t" << dtts1->phiB() << "\t" << dtts1->code() << "\t" << dtts1->Ts2Tag()
                   << "\t" << dtts1->BxCnt() << "\t0" << endl;
            }
          }
        }
      }
    }
  }

  cout << endl;

  cout << "RPC Inputs" << endl;
  cout << "bx\tring\tsector\tstation\troll\tlayer\tstrip\tphi\tlocalX" << endl;
  //cout<<"RPCHitCleaner"<<endl;
  for (auto chamber = rpcDigis->begin(); chamber != rpcDigis->end(); ++chamber) {
    RPCDetId detid = (*chamber).first;
    for (auto digi = (*chamber).second.first; digi != (*chamber).second.second; ++digi) {
      RPCDigi digi_out(digi->strip(), digi->bx());
      //if(digi->bx()!=0) continue;
      int phi = RPCtoDTTranslator::radialAngle(detid, rpcGeometry, digi->strip()) << 2;
      float localx = RPCtoDTTranslator::localX(detid, rpcGeometry, digi->strip());
      cout << digi->bx() << "\t" << detid.ring() << "\t" << detid.sector() - 1 << "\t" << detid.station() << "\t"
           << detid.roll() << "\t" << detid.layer() << "\t" << digi->strip() << "\t" << phi << "\t" << localx << endl;
    }  
  }    
  cout << endl;

  cout << "TwinMux Output" << endl;
  cout << "bx\twheel\tsector\tstation\tphi\tphib\tcode\tts2tag\tbxcnt\trpcbit" << endl;

  for (bx = -2; bx <= 2; bx++) {
    for (wheel = -2; wheel <= 2; wheel++) {
      for (sector = 0; sector < 12; sector++) {
        for (station = 1; station <= 4; station++) {
          auto dtts1 = outphiDigis.chPhiSegm1(wheel, station, sector, bx);
          auto dtts2 = outphiDigis.chPhiSegm2(wheel, station, sector, bx - 1);
          if (dtts1 && dtts1->code() != 7) {
            L1MuDTChambPhDigi dt_ts1 = *dtts1;
            cout << dtts1->bxNum() << "\t" << dtts1->whNum() << "\t" << dtts1->scNum() << "\t" << dtts1->stNum() << "\t"
                 << dtts1->phi() << "\t" << dtts1->phiB() << "\t" << dtts1->code() << "\t" << dtts1->Ts2Tag() << "\t"
                 << dtts1->BxCnt() << "\t" << dtts1->RpcBit() << endl;
          }
          if (dtts2 && dtts2->code() != 7) {
            L1MuDTChambPhDigi dt_ts2 = *dtts2;
            cout << dtts2->bxNum() << "\t" << dtts2->whNum() << "\t" << dtts2->scNum() << "\t" << dtts2->stNum() << "\t"
                 << dtts2->phi() << "\t" << dtts2->phiB() << "\t" << dtts2->code() << "\t" << dtts2->Ts2Tag() << "\t"
                 << dtts2->BxCnt() << "\t" << dtts2->RpcBit() << endl;
          }
        }
      }
    }
  }

  cout << "======================================================" << endl;
}