IOPrinter Class Reference

#include <IOPrinter.h>

Public Member Functions
	IOPrinter ()
void	run (edm::Handle< L1MuDTChambPhContainer >, const L1MuDTChambPhContainer &, edm::Handle< RPCDigiCollection >, const RPCGeometry &)
void	run (L1MuDTChambPhContainer const *, const L1MuDTChambPhContainer &, RPCDigiCollection const *, const RPCGeometry &)
	~IOPrinter ()

Detailed Description

Definition at line 27 of file IOPrinter.h.

Constructor & Destructor Documentation

IOPrinter()

IOPrinter::IOPrinter ( )

inline

Definition at line 29 of file IOPrinter.h.

{};

~IOPrinter()

IOPrinter::~IOPrinter ( )

inline

Definition at line 30 of file IOPrinter.h.

{};

Member Function Documentation

run() [1/2]

void IOPrinter::run	(	edm::Handle< L1MuDTChambPhContainer >	inphiDigis,
		const L1MuDTChambPhContainer &	outphiDigis,
		edm::Handle< RPCDigiCollection >	rpcDigis,
		const RPCGeometry &	rpcGeometry
	)

Align track segments that are coming in bx-1.

for digicout

for chamber

Definition at line 18 of file IOPrinter.cc.

References l1GtPatternGenerator_cfi::bx, relativeConstraints::chamber, L1MuDTChambPhContainer::chPhiSegm1(), L1MuDTChambPhContainer::chPhiSegm2(), gather_cfg::cout, RPCDetId::layer(), RPCtoDTTranslator::localX(), RPCtoDTTranslator::radialAngle(), RPCDetId::ring(), RPCDetId::roll(), hgcalTBTopologyTester_cfi::sector, RPCDetId::sector(), relativeConstraints::station, RPCDetId::station(), and makeMuonMisalignmentScenario::wheel.

Referenced by L1TTwinMuxAlgorithm::run().

                                                    {
  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;
}

run() [2/2]

void IOPrinter::run	(	L1MuDTChambPhContainer const *	inphiDigis,
		const L1MuDTChambPhContainer &	outphiDigis,
		RPCDigiCollection const *	rpcDigis,
		const RPCGeometry &	rpcGeometry
	)

Align track segments that are coming in bx-1.

for digicout

for chamber

Definition at line 99 of file IOPrinter.cc.

References l1GtPatternGenerator_cfi::bx, relativeConstraints::chamber, L1MuTMChambPhContainer::chPhiSegm(), L1MuDTChambPhContainer::chPhiSegm1(), L1MuDTChambPhContainer::chPhiSegm2(), gather_cfg::cout, L1MuDTChambPhContainer::getContainer(), RPCDetId::layer(), RPCtoDTTranslator::localX(), nhits, DTRPCBxCorrection::nRPCHits(), RPCtoDTTranslator::radialAngle(), RPCDetId::ring(), RPCDetId::roll(), hgcalTBTopologyTester_cfi::sector, RPCDetId::sector(), L1MuTMChambPhContainer::setContainer(), relativeConstraints::station, RPCDetId::station(), and makeMuonMisalignmentScenario::wheel.

                                                    {
  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;
}