L1TComparison Class Reference

Inheritance diagram for L1TComparison:

Public Member Functions
	L1TComparison (const ParameterSet &)
	~L1TComparison ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Static Public Member Functions
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)

Private Member Functions
virtual void	analyze (Event const &, EventSetup const &)
virtual void	beginJob ()
virtual void	beginLuminosityBlock (LuminosityBlock const &, EventSetup const &)
virtual void	beginRun (Run const &, EventSetup const &)
virtual void	endJob ()
virtual void	endLuminosityBlock (LuminosityBlock const &, EventSetup const &)
virtual void	endRun (Run const &, EventSetup const &)

Private Attributes
bool	algCheck_
int	algCount_
int	algFails_
edm::EDGetTokenT < GlobalAlgBlkBxCollection >	algTokenA_
edm::EDGetTokenT < GlobalAlgBlkBxCollection >	algTokenB_
bool	bxZeroOnly_
bool	egCheck_
int	egCount_
int	egFails_
edm::EDGetTokenT < EGammaBxCollection >	egTokenA_
edm::EDGetTokenT < EGammaBxCollection >	egTokenB_
bool	jetCheck_
int	jetCount_
int	jetFails_
edm::EDGetTokenT< JetBxCollection >	jetTokenA_
edm::EDGetTokenT< JetBxCollection >	jetTokenB_
bool	muonCheck_
int	muonCount_
int	muonFails_
edm::EDGetTokenT < MuonBxCollection >	muonTokenA_
edm::EDGetTokenT < MuonBxCollection >	muonTokenB_
bool	sumCheck_
int	sumCount_
int	sumFails_
edm::EDGetTokenT < EtSumBxCollection >	sumTokenA_
edm::EDGetTokenT < EtSumBxCollection >	sumTokenB_
string	tag_
bool	tauCheck_
int	tauCount_
int	tauFails_
edm::EDGetTokenT< TauBxCollection >	tauTokenA_
edm::EDGetTokenT< TauBxCollection >	tauTokenB_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 61 of file L1TComparison.cc.

Constructor & Destructor Documentation

L1TComparison::L1TComparison ( const ParameterSet &  iConfig )

explicit

Definition at line 122 of file L1TComparison.cc.

References gather_cfg::cout, and edm::ParameterSet::getParameter().

                                                        { 


  tag_        = iConfig.getParameter<string>("tag");
  egCheck_    = iConfig.getParameter<bool>("egCheck");
  tauCheck_   = iConfig.getParameter<bool>("tauCheck");
  jetCheck_   = iConfig.getParameter<bool>("jetCheck");
  sumCheck_   = iConfig.getParameter<bool>("sumCheck");
  muonCheck_  = iConfig.getParameter<bool>("muonCheck");
  algCheck_  = iConfig.getParameter<bool>("algCheck");
  bxZeroOnly_ = iConfig.getParameter<bool>("bxZeroOnly");

  cout << "L1T Summary for " << tag_ << "\n"; 
  cout << "DEBUG:  egCheck:    " << egCheck_ << "\n";
  cout << "DEBUG:  tauCheck:   " << tauCheck_ << "\n";
  cout << "DEBUG:  jetCheck:   " << jetCheck_ << "\n";
  cout << "DEBUG:  sumCheck:   " << sumCheck_ << "\n";
  cout << "DEBUG:  muonCheck:  " << muonCheck_ << "\n";
  cout << "DEBUG:  algCheck:   " << algCheck_ << "\n";

  if (egCheck_)   {egTokenA_   = consumes<EGammaBxCollection> (iConfig.getParameter<InputTag>("egTagA"));}
  if (tauCheck_)  {tauTokenA_  = consumes<TauBxCollection> (iConfig.getParameter<InputTag>("tauTagA"));}
  if (jetCheck_)  {jetTokenA_  = consumes<JetBxCollection>    (iConfig.getParameter<InputTag>("jetTagA"));}
  if (sumCheck_)  {sumTokenA_  = consumes<EtSumBxCollection>  (iConfig.getParameter<InputTag>("sumTagA"));}
  if (muonCheck_) {muonTokenA_ = consumes<MuonBxCollection>   (iConfig.getParameter<InputTag>("muonTagA"));}
  if (algCheck_)  {algTokenA_  = consumes<GlobalAlgBlkBxCollection>   (iConfig.getParameter<InputTag>("algTagA"));}

  if (egCheck_)   {egTokenB_   = consumes<EGammaBxCollection> (iConfig.getParameter<InputTag>("egTagB"));}
  if (tauCheck_)  {tauTokenB_   = consumes<TauBxCollection> (iConfig.getParameter<InputTag>("tauTagB"));}
  if (jetCheck_)  {jetTokenB_  = consumes<JetBxCollection>    (iConfig.getParameter<InputTag>("jetTagB"));}
  if (sumCheck_)  {sumTokenB_  = consumes<EtSumBxCollection>  (iConfig.getParameter<InputTag>("sumTagB"));}
  if (muonCheck_) {muonTokenB_ = consumes<MuonBxCollection>   (iConfig.getParameter<InputTag>("muonTagB"));}
  if (algCheck_)  {algTokenB_  = consumes<GlobalAlgBlkBxCollection>   (iConfig.getParameter<InputTag>("algTagB"));}

  egCount_   = 0;
  tauCount_  = 0;
  jetCount_  = 0;
  sumCount_  = 0;
  muonCount_ = 0;
  algCount_ = 0;

  egFails_   = 0;
  tauFails_  = 0;
  jetFails_  = 0;
  sumFails_  = 0;
  muonFails_ = 0;
  algFails_ = 0;

}

L1TComparison::~L1TComparison

(

)

Definition at line 172 of file L1TComparison.cc.

{}

Member Function Documentation

void L1TComparison::analyze ( Event const &  iEvent, EventSetup const &  iSetup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 176 of file L1TComparison.cc.

References compare_l1candidate(), gather_cfg::cout, cmsPerfPublish::fail(), edm::Event::getByToken(), and edm::HandleBase::isValid().

{

  cout << "L1TComparison Module output for " << tag_ << "\n";

  if (egCheck_){
    Handle<EGammaBxCollection> XTMPA;    
    iEvent.getByToken(egTokenA_, XTMPA);
    Handle<EGammaBxCollection> XTMPB;    
    iEvent.getByToken(egTokenB_, XTMPB);

    if (!(XTMPA.isValid() && XTMPB.isValid())){
      LogWarning("MissingProduct") << "L1Upgrade e-gamma's not found." << std::endl;
    } else {
      for (int ibx = XTMPA->getFirstBX(); ibx <= XTMPA->getLastBX(); ++ibx) {
    if (bxZeroOnly_ && (ibx != 0)) continue;
    if (ibx < XTMPB->getFirstBX()) continue;
    if (ibx > XTMPB->getLastBX()) continue;
    int sizeA = XTMPA->size(ibx);
    int sizeB = XTMPB->size(ibx);
    if (sizeA != sizeB){
      cout << "L1T COMPARISON FAILURE:  collections have different sizes for bx = " << ibx << "\n"; 
    } else {
      auto itB=XTMPB->begin(ibx);
      for (auto itA=XTMPA->begin(ibx); itA!=XTMPA->end(ibx); ++itA){        
        bool fail = compare_l1candidate(*itA, *itB);
        itB++;
        if (! fail) { egCount_++; }
        else        { egFails_++; }
      }

    }
      }
    }
  }


  if (tauCheck_){
    Handle<TauBxCollection> XTMPA;    
    iEvent.getByToken(tauTokenA_, XTMPA);
    Handle<TauBxCollection> XTMPB;    
    iEvent.getByToken(tauTokenB_, XTMPB);

    if (!(XTMPA.isValid() && XTMPB.isValid())){
      LogWarning("MissingProduct") << "L1Upgrade tau's not found." << std::endl;
    } else {
      for (int ibx = XTMPA->getFirstBX(); ibx <= XTMPA->getLastBX(); ++ibx) {
    if (bxZeroOnly_ && (ibx != 0)) continue;
    if (ibx < XTMPB->getFirstBX()) continue;
    if (ibx > XTMPB->getLastBX()) continue;
    int sizeA = XTMPA->size(ibx);
    int sizeB = XTMPB->size(ibx);
    if (sizeA != sizeB){
      cout << "L1T COMPARISON FAILURE:  collections have different sizes for bx = " << ibx << "\n"; 
    } else {
      auto itB=XTMPB->begin(ibx);
      for (auto itA=XTMPA->begin(ibx); itA!=XTMPA->end(ibx); ++itA){        
        bool fail = compare_l1candidate(*itA, *itB);
        itB++;
        if (! fail) { tauCount_++; }
        else        { tauFails_++; }
      }
    }
      }
    }
  }

  if (jetCheck_){
    Handle<JetBxCollection> XTMPA;    
    iEvent.getByToken(jetTokenA_, XTMPA);
    Handle<JetBxCollection> XTMPB;    
    iEvent.getByToken(jetTokenB_, XTMPB);

    if (!(XTMPA.isValid() && XTMPB.isValid())){
      LogWarning("MissingProduct") << "L1Upgrade jet's not found." << std::endl;
    } else {
      for (int ibx = XTMPA->getFirstBX(); ibx <= XTMPA->getLastBX(); ++ibx) {
    if (bxZeroOnly_ && (ibx != 0)) continue;
    if (ibx < XTMPB->getFirstBX()) continue;
    if (ibx > XTMPB->getLastBX()) continue;
    int sizeA = XTMPA->size(ibx);
    int sizeB = XTMPB->size(ibx);
    if (sizeA != sizeB){
      cout << "L1T COMPARISON FAILURE:  collections have different sizes for bx = " << ibx << "\n"; 
    } else {
      auto itB=XTMPB->begin(ibx);
      for (auto itA=XTMPA->begin(ibx); itA!=XTMPA->end(ibx); ++itA){        
        bool fail = compare_l1candidate(*itA, *itB);
        itB++;
        if (! fail) { jetCount_++; }
        else        { jetFails_++; }
      }
    }
      }
    }
  }

 if (sumCheck_){
    Handle<EtSumBxCollection> XTMPA;    
    iEvent.getByToken(sumTokenA_, XTMPA);
    Handle<EtSumBxCollection> XTMPB;    
    iEvent.getByToken(sumTokenB_, XTMPB);

    if (!(XTMPA.isValid() && XTMPB.isValid())){
      LogWarning("MissingProduct") << "L1Upgrade sum's not found." << std::endl;
    } else {
      for (int ibx = XTMPA->getFirstBX(); ibx <= XTMPA->getLastBX(); ++ibx) {
    if (bxZeroOnly_ && (ibx != 0)) continue;
    if (ibx < XTMPB->getFirstBX()) continue;
    if (ibx > XTMPB->getLastBX()) continue;
    int sizeA = XTMPA->size(ibx);
    int sizeB = XTMPB->size(ibx);

    if (sizeA != sizeB){
      cout << "L1T COMPARISON WARNING:  sums collections have different sizes for bx = " << ibx << "\n"; 
      cout << "L1T COMPARISON WARNING:  sums collections A size  = " << sizeA << "  sums collection B size = " << sizeB << "\n"; 
      cout << "L1T COMPARISON WARNING:  known issue because packer has not been udpated for Minbias\n";
    } 
    for (auto itA=XTMPA->begin(ibx); itA!=XTMPA->end(ibx); ++itA){      
          cout << "L1T COMPARISON :  EtSum type: A = " << itA->getType() <<  "\n";
        }
    for (auto itB=XTMPB->begin(ibx); itB!=XTMPB->end(ibx); ++itB){      
          cout << "L1T COMPARISON :  EtSum type: B = " << itB->getType() <<  "\n";
        }

    // temp workaround for sums not packed...
    if (sizeA > sizeB) sizeA = sizeB;
    if (sizeB > sizeA) sizeB = sizeA;
    
    if (sizeA != sizeB){
      cout << "L1T COMPARISON FAILURE:  collections have different sizes for bx = " << ibx << "\n"; 
    } else {
      auto itB=XTMPB->begin(ibx);
      for (auto itA=XTMPA->begin(ibx); itA!=XTMPA->end(ibx); ++itA){        
          cout << "L1T COMPARISON :  EtSum type: A = " << itA->getType() << " vs B = " << itB->getType() << "\n";
        if (itA->getType() != itB->getType()){
          cout << "L1T COMPARISON FAILURE:  Different types .... EtSum type:" << itA->getType() << " vs " << itB->getType() << "\n";
        }       
        if (itA->getType() == EtSum::kTotalEtEm) 
          cout << "L1T COMPARISON WARNING:  (known issue) sum of type " << itA->getType() << " when emulated has a dummy value (pending proper emulation)" << "\n"; 
        if (itA->getType() < EtSum::kMinBiasHFP0 || itA->getType() > EtSum::kMinBiasHFM1) {
          bool fail = compare_l1candidate(*itA, *itB);
          if (fail){ cout << "L1T COMPARISON FAILURE:  for type " << itA->getType() << "\n";}
          if (! fail) { sumCount_++; }
          else        { sumFails_++; }
        } else {
          cout << "L1T COMPARISON WARNING:  (known issue) not checking sum of type " << itA->getType() << "\n"; 
        }
        itB++;
      }
    }
      }
    }
  }


  if (muonCheck_){
    Handle<MuonBxCollection> XTMPA;    
    iEvent.getByToken(muonTokenA_, XTMPA);
    Handle<MuonBxCollection> XTMPB;    
    iEvent.getByToken(muonTokenB_, XTMPB);

    if (!(XTMPA.isValid() && XTMPB.isValid())){
      LogWarning("MissingProduct") << "L1Upgrade muon's not found." << std::endl;
    } else {
      for (int ibx = XTMPA->getFirstBX(); ibx <= XTMPA->getLastBX(); ++ibx) {
    if (bxZeroOnly_ && (ibx != 0)) continue;
    if (ibx < XTMPB->getFirstBX()) continue;
    if (ibx > XTMPB->getLastBX()) continue;
    int sizeA = XTMPA->size(ibx);
    int sizeB = XTMPB->size(ibx);
    if (sizeA != sizeB){
      cout << "L1T COMPARISON FAILURE:  collections have different sizes for bx = " << ibx << "\n"; 
    } else {
      auto itB=XTMPB->begin(ibx);
      for (auto itA=XTMPA->begin(ibx); itA!=XTMPA->end(ibx); ++itA){        
        bool fail = compare_l1candidate(*itA, *itB);
        itB++;
        if (! fail) { muonCount_++; }
        else        { muonFails_++; }
      }
    }
      }
    }
  }





}

void L1TComparison::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 369 of file L1TComparison.cc.

References gather_cfg::cout.

{
  cout << "INFO:  L1TComparison module beginJob called.\n";
}

void L1TComparison::beginLuminosityBlock ( LuminosityBlock const &  , EventSetup const &    )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 419 of file L1TComparison.cc.

{
}

void L1TComparison::beginRun ( Run const &  run, EventSetup const &  iSetup  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 409 of file L1TComparison.cc.

{
}

void L1TComparison::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 375 of file L1TComparison.cc.

References gather_cfg::cout, and cmsPerfPublish::fail().

                      {
  cout << "INFO:  L1T Comparison for " << tag_ << "\n";
  cout << "INFO: count of successful comparison for each type follows:\n";
  if (egCheck_)   cout << "eg:    " << egCount_ << "\n";
  if (tauCheck_)  cout << "tau:   " << tauCount_ << "\n";
  if (jetCheck_)  cout << "jet:   " << jetCount_ << "\n";
  if (sumCheck_)  cout << "sum:   " << sumCount_ << "\n";
  if (muonCheck_) cout << "muon:  " << muonCount_ << "\n";
  cout << "INFO: count of failed comparison for each type follows:\n";
  if (egCheck_)   cout << "eg:    " << egFails_ << "\n";
  if (tauCheck_)  cout << "tau:   " << tauFails_ << "\n";
  if (jetCheck_)  cout << "jet:   " << jetFails_ << "\n";
  if (sumCheck_)  cout << "sum:   " << sumFails_ << "\n";
  if (muonCheck_) cout << "muon:  " << muonFails_ << "\n";

  int fail = 0;
  if (egCheck_ &&   ((egFails_>0)   || (egCount_<=0))) fail = 1;
  if (tauCheck_ &&  ((tauFails_>0)  || (tauCount_<=0))) fail = 1;
  if (jetCheck_ &&  ((jetFails_>0)  || (jetCount_<=0))) fail = 1;
  if (sumCheck_ &&  ((sumFails_>0)  || (sumCount_<=0))) fail = 1;
  if (muonCheck_ && ((muonFails_>0) || (muonCount_<=0))) fail = 1;

  if (fail){
    cout << "SUMMARY:  L1T Comparison for " << tag_ << " was FAILURE\n";
  } else {
    cout << "SUMMARY:  L1T Comparison for " << tag_ << " was SUCCESS\n";
  }
      



}

void L1TComparison::endLuminosityBlock ( LuminosityBlock const &  , EventSetup const &    )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 424 of file L1TComparison.cc.

{
}

void L1TComparison::endRun ( Run const &  , EventSetup const &    )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 414 of file L1TComparison.cc.

{
}

void L1TComparison::fillDescriptions ( ConfigurationDescriptions &  descriptions )

static

Definition at line 429 of file L1TComparison.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                       {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

Member Data Documentation

bool L1TComparison::algCheck_

private

Definition at line 87 of file L1TComparison.cc.

int L1TComparison::algCount_

private

Definition at line 111 of file L1TComparison.cc.

int L1TComparison::algFails_

private

Definition at line 119 of file L1TComparison.cc.

edm::EDGetTokenT<GlobalAlgBlkBxCollection> L1TComparison::algTokenA_

private

Definition at line 96 of file L1TComparison.cc.

edm::EDGetTokenT<GlobalAlgBlkBxCollection> L1TComparison::algTokenB_

private

Definition at line 103 of file L1TComparison.cc.

bool L1TComparison::bxZeroOnly_

private

Definition at line 88 of file L1TComparison.cc.

bool L1TComparison::egCheck_

private

Definition at line 82 of file L1TComparison.cc.

int L1TComparison::egCount_

private

Definition at line 106 of file L1TComparison.cc.

int L1TComparison::egFails_

private

Definition at line 114 of file L1TComparison.cc.

edm::EDGetTokenT<EGammaBxCollection> L1TComparison::egTokenA_

private

Definition at line 91 of file L1TComparison.cc.

edm::EDGetTokenT<EGammaBxCollection> L1TComparison::egTokenB_

private

Definition at line 98 of file L1TComparison.cc.

bool L1TComparison::jetCheck_

private

Definition at line 84 of file L1TComparison.cc.

int L1TComparison::jetCount_

private

Definition at line 108 of file L1TComparison.cc.

int L1TComparison::jetFails_

private

Definition at line 116 of file L1TComparison.cc.

edm::EDGetTokenT<JetBxCollection> L1TComparison::jetTokenA_

private

Definition at line 93 of file L1TComparison.cc.

edm::EDGetTokenT<JetBxCollection> L1TComparison::jetTokenB_

private

Definition at line 100 of file L1TComparison.cc.

bool L1TComparison::muonCheck_

private

Definition at line 86 of file L1TComparison.cc.

int L1TComparison::muonCount_

private

Definition at line 110 of file L1TComparison.cc.

int L1TComparison::muonFails_

private

Definition at line 118 of file L1TComparison.cc.

edm::EDGetTokenT<MuonBxCollection> L1TComparison::muonTokenA_

private

Definition at line 95 of file L1TComparison.cc.

edm::EDGetTokenT<MuonBxCollection> L1TComparison::muonTokenB_

private

Definition at line 102 of file L1TComparison.cc.

bool L1TComparison::sumCheck_

private

Definition at line 85 of file L1TComparison.cc.

int L1TComparison::sumCount_

private

Definition at line 109 of file L1TComparison.cc.

int L1TComparison::sumFails_

private

Definition at line 117 of file L1TComparison.cc.

edm::EDGetTokenT<EtSumBxCollection> L1TComparison::sumTokenA_

private

Definition at line 94 of file L1TComparison.cc.

edm::EDGetTokenT<EtSumBxCollection> L1TComparison::sumTokenB_

private

Definition at line 101 of file L1TComparison.cc.

string L1TComparison::tag_

private

Definition at line 79 of file L1TComparison.cc.

bool L1TComparison::tauCheck_

private

Definition at line 83 of file L1TComparison.cc.

int L1TComparison::tauCount_

private

Definition at line 107 of file L1TComparison.cc.

int L1TComparison::tauFails_

private

Definition at line 115 of file L1TComparison.cc.

edm::EDGetTokenT<TauBxCollection> L1TComparison::tauTokenA_

private

Definition at line 92 of file L1TComparison.cc.

edm::EDGetTokenT<TauBxCollection> L1TComparison::tauTokenB_

private

Definition at line 99 of file L1TComparison.cc.