IOPrinter Class Reference

#include <IOPrinter.h>

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

Detailed Description

Definition at line 28 of file IOPrinter.h.

Constructor & Destructor Documentation

IOPrinter::IOPrinter ( )

inline

Definition at line 30 of file IOPrinter.h.

{};

IOPrinter::~IOPrinter ( )

inline

Definition at line 31 of file IOPrinter.h.

References run().

{};

Member Function Documentation

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

Align track segments that are coming in bx-1.

for digicout

for chamber

Definition at line 17 of file IOPrinter.cc.

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

Referenced by L1TTwinMuxAlgorithm::run(), and ~IOPrinter().

                                                                                                                                                                          {

    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, c, digi->strip()) << 2 ;
        float localx = RPCtoDTTranslator::localX(detid, c, 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 edm::EventSetup &	c
	)

Align track segments that are coming in bx-1.

for digicout

for chamber

Definition at line 95 of file IOPrinter.cc.

References 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(), 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, c, digi->strip()) <<2 ;
       float localx = RPCtoDTTranslator::localX(detid, c, 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;
}