L1TLSBlock Class Reference

#include <L1TLSBlock.h>

Public Types
enum	BLOCKBY { STATISTICS, N_BLOCKINGBY }
typedef std::pair< int, int >	LumiRange
typedef std::vector< LumiRange >	LumiRangeList
typedef std::vector< std::pair< int, double > >	LumiTestDoubleList
typedef std::vector< std::pair< int, double > >	LumiTestIntList

Public Member Functions
LumiRangeList	doBlocking (const LumiTestDoubleList &, double, BLOCKBY)
LumiRangeList	doBlocking (const LumiTestIntList &, int, BLOCKBY)
	L1TLSBlock ()
virtual	~L1TLSBlock ()

Private Member Functions
void	blockByStatistics ()
double	computeErrorFromRange (LumiRange &)
void	initializeIO (bool)
void	orderTestDoubleList ()
void	orderTestIntList ()

Private Attributes
LumiTestDoubleList	inputDoubleList_
LumiTestIntList	inputIntList_
LumiRangeList	outputList_
double	thresholdD_
int	thresholdI_

Detailed Description

Description: offline DQM class for LS blocking

Implementation: <TODO: enter implementation details>

Author

: Pietro Vischia - LIP Lisbon pietr.nosp@m.o.vi.nosp@m.schia.nosp@m.@gma.nosp@m.il.co.nosp@m.m

Changelog: 2012/10/23 12:01:01: Creation, infrastructure and blocking by statistics

Todo:

• Activate the "method" string usage instead of base one
• improve the switch method (instead of calling functions one should define a template - it is probably an uberpain, but would be more neat)
• in doBlocking, substitute the switch default case with a cms::Exception
• Cleanup includes
• Add other blocking methods

Description: offline DQM class for LS blocking

Implementation: <TODO: enter implementation details>

Author

: Pietro Vischia - LIP Lisbon pietr.nosp@m.o.vi.nosp@m.schia.nosp@m.@gma.nosp@m.il.co.nosp@m.m

Changelog: 2012/10/23 12:01:01: Class

Todo: see header file

Definition at line 72 of file L1TLSBlock.h.

Member Typedef Documentation

typedef std::pair<int,int> L1TLSBlock::LumiRange

Definition at line 78 of file L1TLSBlock.h.

typedef std::vector<LumiRange> L1TLSBlock::LumiRangeList

Definition at line 79 of file L1TLSBlock.h.

typedef std::vector<std::pair<int,double> > L1TLSBlock::LumiTestDoubleList

Definition at line 76 of file L1TLSBlock.h.

typedef std::vector<std::pair<int,double> > L1TLSBlock::LumiTestIntList

Definition at line 77 of file L1TLSBlock.h.

Member Enumeration Documentation

enum L1TLSBlock::BLOCKBY

Enumerator
STATISTICS
N_BLOCKINGBY

Definition at line 81 of file L1TLSBlock.h.

{STATISTICS, N_BLOCKINGBY};

Constructor & Destructor Documentation

L1TLSBlock::L1TLSBlock

(

)

Definition at line 62 of file L1TLSBlock.cc.

                      {

  // Initialize LumiRangeList object. This will be redone at every doBlocking call. Perhaps rethink this
  initializeIO(false);
  
}

L1TLSBlock::~L1TLSBlock ( )

virtual

Definition at line 72 of file L1TLSBlock.cc.

{}

Member Function Documentation

void L1TLSBlock::blockByStatistics ( )

private

Definition at line 150 of file L1TLSBlock.cc.

References mps_fire::i.

                                  {
  LumiRange currentRange;
  double currentError(0);
  bool resetFlag(true);

  for(LumiTestDoubleList::iterator i=inputDoubleList_.begin(); i!=inputDoubleList_.end(); i++){
    if(resetFlag){
      currentRange = std::make_pair(i->first,i->first);
      resetFlag = false;
    }
    else
      currentRange = std::make_pair(currentRange.first,i->first);
    currentError = computeErrorFromRange(currentRange);
    if(currentError < thresholdD_){
      outputList_.push_back(currentRange);
      resetFlag = true;
    }
    
  }
}

double L1TLSBlock::computeErrorFromRange ( LumiRange &  lumiRange )

private

Definition at line 171 of file L1TLSBlock.cc.

References relativeConstraints::error, plotBeamSpotDB::first, mps_fire::i, and edm::second().

                                                            {
  std::vector<double> errorList;
  errorList.clear();

  for(size_t i=0; i < inputDoubleList_.size(); i++){
    if(inputDoubleList_[i].first>lumiRange.first && inputDoubleList_[i].first<lumiRange.second)
      errorList.push_back(inputDoubleList_[i].second);
  }
  
  double error(-1);
  for(size_t i=0; i<errorList.size(); i++)
    error += 1 / (errorList[i] * errorList[i] );
  return error;

}

L1TLSBlock::LumiRangeList L1TLSBlock::doBlocking	(	const LumiTestDoubleList &	inputList,
		double	threshold,
		BLOCKBY	blockingMethod
	)

Definition at line 77 of file L1TLSBlock.cc.

References gather_cfg::cout, and electronIdCutBased_cfi::threshold.

{
  inputDoubleList_ = inputList;
  thresholdD_ = threshold;
  
  initializeIO(true);
  
  orderTestDoubleList();
  
  switch(blockingMethod){
  case STATISTICS:
    blockByStatistics();
    break;
  default:
    cout << "[L1TLSBlock::doBlocking()]: Blocking method does not exist or is not implemented" << endl;
  }
  
  
  return outputList_;
}

L1TLSBlock::LumiRangeList L1TLSBlock::doBlocking	(	const LumiTestIntList &	inputList,
		int	threshold,
		BLOCKBY	blockingMethod
	)

Definition at line 101 of file L1TLSBlock.cc.

References gather_cfg::cout, and electronIdCutBased_cfi::threshold.

{
  inputIntList_ = inputList;
  thresholdI_ = threshold;
  
  initializeIO(true);
  
  orderTestIntList();
  
  switch(blockingMethod){
  case STATISTICS:
    cout << "[L1TLSBlock::doBlocking()]: Blocking by statistics require doubles as inputs for test variable and threshold" << endl;
    break;
  default:
    cout << "[L1TLSBlock::doBlocking()]: Blocking method does not exist or is not implemented" << endl;
  }
  
  return outputList_;
}

void L1TLSBlock::initializeIO ( bool  outputOnly )

private

Definition at line 124 of file L1TLSBlock.cc.

                                            {
  if(!outputOnly){
    inputIntList_.clear();
    inputDoubleList_.clear();
  }
  outputList_.clear();
}

void L1TLSBlock::orderTestDoubleList ( )

private

Definition at line 135 of file L1TLSBlock.cc.

References jetUpdater_cfi::sort.

                                    {
  std::sort(inputDoubleList_.begin(), inputDoubleList_.end(), sort_pair_first<int, double>());
  // std::sort(v.begin(), v.end(), sort_pair_second<int, std::greater<int> >());
}

void L1TLSBlock::orderTestIntList ( )

private

Definition at line 143 of file L1TLSBlock.cc.

References jetUpdater_cfi::sort.

                                 {
  std::sort(inputIntList_.begin(), inputIntList_.end(), sort_pair_first<int, int>());
}

Member Data Documentation

LumiTestDoubleList L1TLSBlock::inputDoubleList_

private

Definition at line 102 of file L1TLSBlock.h.

LumiTestIntList L1TLSBlock::inputIntList_

private

Definition at line 101 of file L1TLSBlock.h.

LumiRangeList L1TLSBlock::outputList_

private

Definition at line 103 of file L1TLSBlock.h.

double L1TLSBlock::thresholdD_

private

Definition at line 104 of file L1TLSBlock.h.

int L1TLSBlock::thresholdI_

private

Definition at line 105 of file L1TLSBlock.h.