#include <CondFormats/L1TObjects/interface/L1RCTParameters.h>

Public Member Functions
unsigned short	calcCard (unsigned short rct_iphi, unsigned short absIeta) const
unsigned short	calcCrate (unsigned short rct_iphi, short ieta) const
unsigned short	calcIAbsEta (unsigned short iCrate, unsigned short iCard, unsigned short iTower) const
short	calcIEta (unsigned short iCrate, unsigned short iCard, unsigned short iTower) const
unsigned short	calcIPhi (unsigned short iCrate, unsigned short iCard, unsigned short iTower) const
unsigned short	calcTower (unsigned short rct_iphi, unsigned short absIeta) const
double	eActivityCut () const
const std::vector< double > &	eGammaECalScaleFactors () const
const std::vector< double > &	eGammaHCalScaleFactors () const
double	eGammaLSB () const
float	EGammaTPGSum (const float &ecal, const float &hcal, const unsigned &iAbsEta) const
unsigned	eicIsolationThreshold () const
double	eMaxForFGCut () const
double	eMaxForHoECut () const
double	eMinForFGCut () const
double	eMinForHoECut () const
double	hActivityCut () const
double	hMinForHoECut () const
double	hOeCut () const
const std::vector< double > &	jetMETECalScaleFactors () const
const std::vector< double > &	jetMETHCalScaleFactors () const
double	jetMETLSB () const
float	JetMETTPGSum (const float &ecal, const float &hcal, const unsigned &iAbsEta) const
unsigned	jscQuietThresholdBarrel () const
unsigned	jscQuietThresholdEndcap () const
	L1RCTParameters (double eGammaLSB, double jetMETLSB, double eMinForFGCut, double eMaxForFGCut, double hOeCut, double eMinForHoECut, double eMaxForHoECut, double hMinForHoECut, double eActivityCut, double hActivityCut, unsigned eicIsolationThreshold, unsigned jscQuietThresholdBarrel, unsigned jscQuietThresholdEndcap, bool noiseVetoHB, bool noiseVetoHEplus, bool noiseVetoHEminus, bool useLindsey, const std::vector< double > &eGammaECalScaleFactors, const std::vector< double > &eGammaHCalScaleFactors, const std::vector< double > &jetMETECalScaleFactors, const std::vector< double > &jetMETHCalScaleFactors, const std::vector< double > &ecal_calib, const std::vector< double > &hcal_calib, const std::vector< double > &hcal_high_calib, const std::vector< double > &cross_terms, const std::vector< double > &lowHoverE_smear, const std::vector< double > &highHoverE_smear)
	L1RCTParameters ()
bool	noiseVetoHB () const
bool	noiseVetoHEminus () const
bool	noiseVetoHEplus () const
void	print (std::ostream &s) const
	~L1RCTParameters ()
Private Member Functions
float	correctedTPGSum (const float &ecal, const float &hcal, const unsigned &index) const
Private Attributes
std::vector< std::vector < double > >	cross_terms_
double	eActivityCut_
std::vector< std::vector < double > >	ecal_calib_
std::vector< double >	eGammaECalScaleFactors_
std::vector< double >	eGammaHCalScaleFactors_
double	eGammaLSB_
unsigned	eicIsolationThreshold_
double	eMaxForFGCut_
double	eMaxForHoECut_
double	eMinForFGCut_
double	eMinForHoECut_
double	hActivityCut_
std::vector< std::vector < double > >	hcal_calib_
std::vector< std::vector < double > >	hcal_high_calib_
double	hMinForHoECut_
double	hOeCut_
std::vector< double >	HoverE_smear_high_
std::vector< double >	HoverE_smear_low_
std::vector< double >	jetMETECalScaleFactors_
std::vector< double >	jetMETHCalScaleFactors_
double	jetMETLSB_
unsigned	jscQuietThresholdBarrel_
unsigned	jscQuietThresholdEndcap_
bool	noiseVetoHB_
bool	noiseVetoHEminus_
bool	noiseVetoHEplus_
bool	useCorrections_

Detailed Description

Description: Class to contain parameters which define RCT Lookup Tables

Usage: <usage>

Definition at line 26 of file L1RCTParameters.h.

Constructor & Destructor Documentation

L1RCTParameters::L1RCTParameters ( ) [inline]

Definition at line 31 of file L1RCTParameters.h.

{}

L1RCTParameters::L1RCTParameters	(	double	eGammaLSB,
		double	jetMETLSB,
		double	eMinForFGCut,
		double	eMaxForFGCut,
		double	hOeCut,
		double	eMinForHoECut,
		double	eMaxForHoECut,
		double	hMinForHoECut,
		double	eActivityCut,
		double	hActivityCut,
		unsigned	eicIsolationThreshold,
		unsigned	jscQuietThresholdBarrel,
		unsigned	jscQuietThresholdEndcap,
		bool	noiseVetoHB,
		bool	noiseVetoHEplus,
		bool	noiseVetoHEminus,
		bool	useLindsey,
		const std::vector< double > &	eGammaECalScaleFactors,
		const std::vector< double > &	eGammaHCalScaleFactors,
		const std::vector< double > &	jetMETECalScaleFactors,
		const std::vector< double > &	jetMETHCalScaleFactors,
		const std::vector< double > &	ecal_calib,
		const std::vector< double > &	hcal_calib,
		const std::vector< double > &	hcal_high_calib,
		const std::vector< double > &	cross_terms,
		const std::vector< double > &	lowHoverE_smear,
		const std::vector< double > &	highHoverE_smear
	)

Definition at line 17 of file L1RCTParameters.cc.

References cross_terms_, ecal_calib_, hcal_calib_, hcal_high_calib_, and i.

                                   :
  eGammaLSB_(eGammaLSB),
  jetMETLSB_(jetMETLSB),
  eMinForFGCut_(eMinForFGCut),
  eMaxForFGCut_(eMaxForFGCut),
  hOeCut_(hOeCut),
  eMinForHoECut_(eMinForHoECut),
  eMaxForHoECut_(eMaxForHoECut),
  hMinForHoECut_(hMinForHoECut),
  eActivityCut_(eActivityCut),
  hActivityCut_(hActivityCut),
  eicIsolationThreshold_(eicIsolationThreshold),
  jscQuietThresholdBarrel_(jscQuietThresholdBarrel),
  jscQuietThresholdEndcap_(jscQuietThresholdEndcap),
  noiseVetoHB_(noiseVetoHB),
  noiseVetoHEplus_(noiseVetoHEplus),
  noiseVetoHEminus_(noiseVetoHEminus),
  useCorrections_(useCorrections),
  eGammaECalScaleFactors_(eGammaECalScaleFactors),
  eGammaHCalScaleFactors_(eGammaHCalScaleFactors),
  jetMETECalScaleFactors_(jetMETECalScaleFactors),
  jetMETHCalScaleFactors_(jetMETHCalScaleFactors),
  HoverE_smear_low_(lowHoverE_smear),
  HoverE_smear_high_(highHoverE_smear)
{
  ecal_calib_.resize(28);
  hcal_calib_.resize(28);
  hcal_high_calib_.resize(28);
  cross_terms_.resize(28);

  for(unsigned i = 0; i < ecal_calib.size(); ++i)
    ecal_calib_[i/3].push_back(ecal_calib[i]);
  for(unsigned i = 0; i < hcal_calib.size(); ++i)
    hcal_calib_[i/3].push_back(hcal_calib[i]);
  for(unsigned i = 0; i < hcal_high_calib.size(); ++i)
    hcal_high_calib_[i/3].push_back(hcal_high_calib[i]);
  for(unsigned i = 0; i < cross_terms.size(); ++i)
    cross_terms_[i/6].push_back(cross_terms[i]);
}

L1RCTParameters::~L1RCTParameters ( ) [inline]

Definition at line 63 of file L1RCTParameters.h.

{;}

Member Function Documentation

unsigned short L1RCTParameters::calcCard	(	unsigned short	rct_iphi,
		unsigned short	absIeta
	)		const

Definition at line 96 of file L1RCTParameters.cc.

Referenced by L1RCT::digiInput().

{
  unsigned short card = 999;
  // Note absIeta counts from 1-32 (not 0-31)
  if (absIeta <= 24){
    card =  ((absIeta-1)/8)*2 + (rct_iphi%8)/4;
  }
  // 25 <= absIeta <= 28 (card 6)
  else if ((absIeta >= 25) && (absIeta <= 28)){
    card = 6;
  }
  else{}
  return card;
}

unsigned short L1RCTParameters::calcCrate	(	unsigned short	rct_iphi,
		short	ieta
	)		const

Definition at line 85 of file L1RCTParameters.cc.

References abs.

Referenced by L1RCT::digiInput(), and L1RCTLutWriter::writeRcLutFile().

{
  unsigned short crate = rct_iphi/8;
  if(abs(ieta) > 28) crate = rct_iphi / 2;
  if (ieta > 0){
    crate = crate + 9;
  }
  return crate;
}

unsigned short L1RCTParameters::calcIAbsEta	(	unsigned short	iCrate,
		unsigned short	iCard,
		unsigned short	iTower
	)		const

Definition at line 181 of file L1RCTParameters.cc.

Referenced by calcIEta().

{
  unsigned short absIEta;
  if(iCard < 6) 
    absIEta = (iCard / 2) * 8 + (iTower / 4) + 1;  // rm -1 JLL 7Nov07
  else if(iCard == 6) {
    // card 6, region 0
    if(iTower < 16) // 17->16
      absIEta = 25 + iTower / 4;  // rm -1
    // card 6, region 1
    else
      absIEta = 28 - ((iTower - 16) / 4);  // 17->16
  }
  // HF regions
  else
    absIEta = 29 + iTower % 4;
  return absIEta;
}

short L1RCTParameters::calcIEta	(	unsigned short	iCrate,
		unsigned short	iCard,
		unsigned short	iTower
	)		const

Definition at line 148 of file L1RCTParameters.cc.

References calcIAbsEta().

Referenced by L1RCTLookupTables::lookup(), and RctInputTextToDigi::produce().

{
  unsigned short absIEta = calcIAbsEta(iCrate, iCard, iTower);
  short iEta;
  if(iCrate < 9) iEta = -absIEta;
  else iEta = absIEta;
  return iEta;
}

unsigned short L1RCTParameters::calcIPhi	(	unsigned short	iCrate,
		unsigned short	iCard,
		unsigned short	iTower
	)		const

Definition at line 159 of file L1RCTParameters.cc.

Referenced by L1RCTLookupTables::lookup(), and RctInputTextToDigi::produce().

{
  short iPhi;
  if(iCard < 6)
    iPhi = (iCrate % 9) * 8 + (iCard % 2) * 4 + (iTower % 4); // rm -1 7Nov07
  else if(iCard == 6){
    // region 0
    if(iTower < 16)  // 17->16
      iPhi = (iCrate % 9) * 8 + (iTower % 4);  // rm -1
    // region 1
    else
      iPhi = (iCrate % 9) * 8 + ((iTower - 16) % 4) + 4; // 17 -> 16
  }
  // HF regions
  else
    iPhi = (iCrate % 9) * 2 + iTower / 4;
  return iPhi;
}

unsigned short L1RCTParameters::calcTower	(	unsigned short	rct_iphi,
		unsigned short	absIeta
	)		const

Definition at line 113 of file L1RCTParameters.cc.

Referenced by L1RCT::digiInput(), and L1RCTLutWriter::writeRcLutFile().

{
  unsigned short tower = 999;
  unsigned short iphi = rct_iphi;
  unsigned short regionPhi = (iphi % 8)/4;

  // Note absIeta counts from 1-32 (not 0-31)
  if (absIeta <= 24){
    // assume iphi between 0 and 71; makes towers from 0-31, mod. 7Nov07
    tower = ((absIeta-1)%8)*4 + (iphi%4);  // REMOVED +1
  }
  // 25 <= absIeta <= 28 (card 6)
  else if ((absIeta >= 25) && (absIeta <= 28)){
    if (regionPhi == 0){
      // towers from 0-31, modified 7Nov07 Jessica Leonard
      tower = (absIeta-25)*4 + (iphi%4);  // REMOVED +1
    }
    else {
      tower = 28 + iphi % 4 + (25 - absIeta) * 4;  // mod. 7Nov07 JLL
    }
  }
  // absIeta >= 29 (HF regions)
  else if ((absIeta >= 29) && (absIeta <= 32)){
    // SPECIAL DEFINITION OF REGIONPHI FOR HF SINCE HF IPHI IS 0-17 
    // Sept. 19 J. Leonard
    regionPhi = iphi % 2;
    // HF MAPPING, just regions now, don't need to worry about towers
    // just calling it "tower" for convenience
    tower = (regionPhi) * 4 + absIeta - 29;
  }
  return tower;
}

float L1RCTParameters::correctedTPGSum	(	const float &	ecal,
		const float &	hcal,
		const unsigned &	index
	)		const `[private]`

Definition at line 239 of file L1RCTParameters.cc.

References trackerHits::c, cross_terms_, alignCSCRings::e, patCandidatesForDimuonsSequences_cff::ecal, ecal_calib_, h, patCandidatesForDimuonsSequences_cff::hcal, hcal_calib_, hcal_high_calib_, HoverE_smear_high_, HoverE_smear_low_, getHLTprescales::index, and funct::pow().

Referenced by EGammaTPGSum(), and JetMETTPGSum().

{
  if(index >= 28 && ecal > 120 && hcal > 120) return (ecal + hcal); // return plain sum if outside of calibration range or index is too high

  // let's make sure we're asking for items that are there.
  if(ecal_calib_.at(index).size() != 3 || hcal_calib_.at(index).size() != 3 ||
     hcal_high_calib_.at(index).size() != 3 || cross_terms_.at(index).size() != 6 ||
     HoverE_smear_high_.size() <= index || HoverE_smear_low_.size() <= index)
    return (ecal+hcal);

  double e = ecal, h = hcal;
  double ec = 0.0, hc = 0.0, c = 0.0;
  
  ec = (ecal_calib_.at(index).at(0)*std::pow(e,3.) +
        ecal_calib_.at(index).at(1)*std::pow(e,2.) +
        ecal_calib_.at(index).at(2)*e);
  
  if(e + h < 23)
    {
      hc = (hcal_calib_.at(index).at(0)*std::pow(h,3.) + 
            hcal_calib_.at(index).at(1)*std::pow(h,2.) + 
            hcal_calib_.at(index).at(2)*h);
      
      c = (cross_terms_.at(index).at(0)*std::pow(e, 2.)*h +
           cross_terms_.at(index).at(1)*std::pow(h, 2.)*e +
           cross_terms_.at(index).at(2)*e*h +
           cross_terms_.at(index).at(3)*std::pow(e, 3.)*h +
           cross_terms_.at(index).at(4)*std::pow(h, 3.)*e +
           cross_terms_.at(index).at(5)*std::pow(h, 2.)*std::pow(e, 2.));
    }
  else
    {
      hc = (hcal_high_calib_.at(index).at(0)*std::pow(h,3.) + 
            hcal_high_calib_.at(index).at(1)*std::pow(h,2.) + 
            hcal_high_calib_.at(index).at(2)*h);
    }
    
  if(h/(e+h) >= 0.05)
    {
      ec *= HoverE_smear_high_.at(index);
      hc *= HoverE_smear_high_.at(index);
      c *= HoverE_smear_high_.at(index);
    }
  else
    {
      ec *= HoverE_smear_low_.at(index);
    }
  return ec + hc + c;
}

double L1RCTParameters::eActivityCut ( ) const [inline]

Definition at line 75 of file L1RCTParameters.h.

References eActivityCut_.

Referenced by L1RCTLookupTables::activityBit().

{return eActivityCut_;}

const std::vector<double>& L1RCTParameters::eGammaECalScaleFactors ( ) const [inline]

Definition at line 83 of file L1RCTParameters.h.

References eGammaECalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

{return eGammaECalScaleFactors_;}

const std::vector<double>& L1RCTParameters::eGammaHCalScaleFactors ( ) const [inline]

Definition at line 84 of file L1RCTParameters.h.

References eGammaHCalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

{return eGammaHCalScaleFactors_;}

double L1RCTParameters::eGammaLSB ( ) const [inline]

Definition at line 67 of file L1RCTParameters.h.

References eGammaLSB_.

Referenced by L1RCTLookupTables::convertEcal(), and L1RCTLookupTables::eGammaETCode().

{return eGammaLSB_;}

float L1RCTParameters::EGammaTPGSum	(	const float &	ecal,
		const float &	hcal,
		const unsigned &	iAbsEta
	)		const

Definition at line 221 of file L1RCTParameters.cc.

References correctedTPGSum(), eGammaECalScaleFactors_, eGammaHCalScaleFactors_, query::result, and useCorrections_.

Referenced by L1RCTLookupTables::eGammaETCode().

{
  float ecal_c = ecal*eGammaECalScaleFactors_.at(iAbsEta-1);
  float hcal_c = hcal*eGammaHCalScaleFactors_.at(iAbsEta-1);
  float result = ecal_c + hcal_c;

  if(useCorrections_)
    {
      if(eGammaHCalScaleFactors_.at(iAbsEta-1) != 0)
        hcal_c = hcal;
      
      result = correctedTPGSum(ecal_c,hcal_c,iAbsEta-1);
    }

  return result;
}

unsigned L1RCTParameters::eicIsolationThreshold ( ) const [inline]

Definition at line 77 of file L1RCTParameters.h.

References eicIsolationThreshold_.

Referenced by L1RCTLutWriter::analyze(), and L1RCTElectronIsolationCard::calcElectronCandidates().

{return eicIsolationThreshold_;}

double L1RCTParameters::eMaxForFGCut ( ) const [inline]

Definition at line 70 of file L1RCTParameters.h.

References eMaxForFGCut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

{return eMaxForFGCut_;}

double L1RCTParameters::eMaxForHoECut ( ) const [inline]

Definition at line 73 of file L1RCTParameters.h.

References eMaxForHoECut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

{return eMaxForHoECut_;}

double L1RCTParameters::eMinForFGCut ( ) const [inline]

Definition at line 69 of file L1RCTParameters.h.

References eMinForFGCut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

{return eMinForFGCut_;}

double L1RCTParameters::eMinForHoECut ( ) const [inline]

Definition at line 72 of file L1RCTParameters.h.

References eMinForHoECut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

{return eMinForHoECut_;}

double L1RCTParameters::hActivityCut ( ) const [inline]

Definition at line 76 of file L1RCTParameters.h.

References hActivityCut_.

Referenced by L1RCTLookupTables::activityBit().

{return hActivityCut_;}

double L1RCTParameters::hMinForHoECut ( ) const [inline]

Definition at line 74 of file L1RCTParameters.h.

References hMinForHoECut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

{return hMinForHoECut_;}

double L1RCTParameters::hOeCut ( ) const [inline]

Definition at line 71 of file L1RCTParameters.h.

References hOeCut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

{return hOeCut_;}

const std::vector<double>& L1RCTParameters::jetMETECalScaleFactors ( ) const [inline]

Definition at line 85 of file L1RCTParameters.h.

References jetMETECalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

{return jetMETECalScaleFactors_;}

const std::vector<double>& L1RCTParameters::jetMETHCalScaleFactors ( ) const [inline]

Definition at line 86 of file L1RCTParameters.h.

References jetMETHCalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

{return jetMETHCalScaleFactors_;}

double L1RCTParameters::jetMETLSB ( ) const [inline]

Definition at line 68 of file L1RCTParameters.h.

References jetMETLSB_.

Referenced by L1RCTLookupTables::convertHcal(), L1RCTLookupTables::jetMETETCode(), and L1RCTLookupTables::lookup().

{return jetMETLSB_;}

float L1RCTParameters::JetMETTPGSum	(	const float &	ecal,
		const float &	hcal,
		const unsigned &	iAbsEta
	)		const

Definition at line 201 of file L1RCTParameters.cc.

References correctedTPGSum(), eGammaECalScaleFactors_, jetMETECalScaleFactors_, jetMETHCalScaleFactors_, query::result, and useCorrections_.

Referenced by L1RCTLookupTables::jetMETETCode().

{
  float ecal_c = ecal*jetMETECalScaleFactors_.at(iAbsEta-1);
  float hcal_c = hcal*jetMETHCalScaleFactors_.at(iAbsEta-1);
  float result = ecal_c + hcal_c;

  if(useCorrections_)
    {
      if(jetMETHCalScaleFactors_.at(iAbsEta-1) != 0)
        hcal_c = hcal;

      if(jetMETECalScaleFactors_.at(iAbsEta-1) != 0)
        ecal_c = ecal*eGammaECalScaleFactors_.at(iAbsEta-1); // Use eGamma Corrections

      result = correctedTPGSum(ecal_c,hcal_c,iAbsEta-1);
    }

  return result;
}

unsigned L1RCTParameters::jscQuietThresholdBarrel ( ) const [inline]

Definition at line 78 of file L1RCTParameters.h.

References jscQuietThresholdBarrel_.

Referenced by L1RCTLutWriter::analyze(), and L1RCTJetSummaryCard::fillQuietBits().

{return jscQuietThresholdBarrel_;}

unsigned L1RCTParameters::jscQuietThresholdEndcap ( ) const [inline]

Definition at line 79 of file L1RCTParameters.h.

References jscQuietThresholdEndcap_.

Referenced by L1RCTLutWriter::analyze(), and L1RCTJetSummaryCard::fillQuietBits().

{return jscQuietThresholdEndcap_;}

bool L1RCTParameters::noiseVetoHB ( ) const [inline]

Definition at line 80 of file L1RCTParameters.h.

References noiseVetoHB_.

Referenced by L1RCTLookupTables::lookup().

{return noiseVetoHB_;}

bool L1RCTParameters::noiseVetoHEminus ( ) const [inline]

Definition at line 82 of file L1RCTParameters.h.

References noiseVetoHEminus_.

Referenced by L1RCTLookupTables::lookup().

{return noiseVetoHEminus_;}

bool L1RCTParameters::noiseVetoHEplus ( ) const [inline]

Definition at line 81 of file L1RCTParameters.h.

References noiseVetoHEplus_.

Referenced by L1RCTLookupTables::lookup().

{return noiseVetoHEplus_;}

void L1RCTParameters::print ( std::ostream & s ) const

Definition at line 290 of file L1RCTParameters.cc.

References cross_terms_, eActivityCut_, ecal_calib_, eGammaECalScaleFactors_, eGammaHCalScaleFactors_, eGammaLSB_, eicIsolationThreshold_, eMaxForFGCut_, eMaxForHoECut_, eMinForFGCut_, eMinForHoECut_, end, hActivityCut_, hcal_calib_, hcal_high_calib_, hMinForHoECut_, hOeCut_, HoverE_smear_high_, HoverE_smear_low_, i, j, jetMETECalScaleFactors_, jetMETHCalScaleFactors_, jetMETLSB_, jscQuietThresholdBarrel_, jscQuietThresholdEndcap_, noiseVetoHB_, noiseVetoHEminus_, noiseVetoHEplus_, and useCorrections_.

                                           {

    s << "\nPrinting record L1RCTParametersRcd" <<endl;
    s << "\n\"Parameter description\" \n  \"Parameter name\"  \"Value\" "
            << endl;
    s << "\ne/gamma least significant bit energy transmitted from receiver cards to EIC cards. \n  "
            << "eGammaLSB = " << eGammaLSB_ << endl ;
    s << "\nLSB of region Et scale from RCT to GCT (GeV) \n  "
            << "jetMETLSB = " << jetMETLSB_ << endl ;
    s << "\nminimum ECAL Et for which fine-grain veto is applied (GeV) \n  "
            << " eMinForFGCut = " << eMinForFGCut_ << endl ;
    s << "\nmaximum ECAL Et for which fine-grain veto is applied (GeV) \n  "
            << "eMaxForFGCut = " << eMaxForFGCut_ << endl ;
    s << "\nmaximum value of (HCAL Et / ECAL Et) \n  "
            << "hOeCut = " << hOeCut_ << endl ;
    s << "\nminimum ECAL Et for which H/E veto is applied (GeV) \n  "
            << "eMinForHoECut = " << eMinForHoECut_ << endl ;
    s << "\nmaximum ECAL Et for which H/E veto is applied (GeV) \n  "
            << "eMaxForHoECut = " << eMaxForHoECut_ << endl ;
    s << "\nminimum HCAL Et for which H/E veto is applied (GeV)  \n  "
            << "hMinForHoECut = " << hMinForHoECut_ << endl ;
    s << "\nECAL Et threshold above which tau activity bit is set (GeV)  \n  "
            << "eActivityCut = " << eActivityCut_ << endl ;
    s << "\nHCAL Et threshold above which tau activity bit is set (GeV)  \n  "
            << "hActivityCut = " << hActivityCut_ << endl ;
    s << "\nNeighboring trigger tower energy minimum threshold that marks candidate as non-isolated. (LSB bits) \n  "
            << "eicIsolationThreshold = " << eicIsolationThreshold_ << endl ;
    s << "\nIf jetMet energy in RCT Barrel Region is below this value, a quiet bit is set. (LSB bits)\n  "
            << "jscQuietThreshBarrel = " << jscQuietThresholdBarrel_ << endl ;
    s << "\nIf jetMet energy in RCT Endcap Region is below this value, a quiet bit is set. (LSB bits) \n  "
            << "jscQuietThreshEndcap = " << jscQuietThresholdEndcap_ << endl ;
    s << "\nWhen set to TRUE, HCAL energy is ignored if no ECAL energy is present in corresponding trigger tower for RCT Barrel \n  "
            << "noiseVetoHB = " << noiseVetoHB_ << endl ;
    s << "\nWhen set to TRUE, HCAL energy is ignored if no ECAL energy is present in corresponding trigger tower for RCT Encap+ \n  "
            << "noiseVetoHEplus = " << noiseVetoHEplus_ << endl ;
    s << "\nWhen set to TRUE, HCAL energy is ignored if no ECAL energy is present in corresponding trigger tower for RCT Endcap- \n  "
            << "noiseVetoHEminus = " << noiseVetoHEminus_ << endl ;

    s << "\n\neta-dependent multiplicative factors for ECAL Et before summation \n  "
            << "eGammaECal Scale Factors " << endl;
    s << "ieta  ScaleFactor" <<endl;
    for(int i = 0 ; i<28; i++)
        s << setw(4) << i+1 << "  " << eGammaECalScaleFactors_.at(i) <<endl;

    s << "\n\neta-dependent multiplicative factors for HCAL Et before summation \n  "
            <<"eGammaHCal Scale Factors "<<endl;
    s << "ieta  ScaleFactor" <<endl;
    for(int i = 0 ; i<28; i++)
        s << setw(4) << i+1 << "  " << eGammaHCalScaleFactors_.at(i) <<endl;
     
    s << "\n\neta-dependent multiplicative factors for ECAL Et before summation \n  "
            <<"jetMETECal Scale Factors "<<endl;
    s << "ieta  ScaleFactor" <<endl;
    for(int i = 0 ; i<28; i++)
        s<< setw(4) << i+1 << "  " << jetMETECalScaleFactors_.at(i) <<endl;
     
    s << "\n\neta-dependent multiplicative factors for HCAL Et before summation \n"
            <<"jetMETHCal Scale Factors "<<endl;
    s << "ieta  ScaleFactor" <<endl;
    for(int i = 0 ; i<28; i++)
        s << setw(4) <<i+1 << "  " << jetMETHCalScaleFactors_.at(i) <<endl;


    if(useCorrections_) {
        s<< "\n\nUSING calibration variables " <<endl;


        s << "\n\nH over E smear low Correction Factors "<<endl;
        s << "ieta  Correction Factor" <<endl;
        for(int i = 0 ; i<28; i++)
            s << setw(4) <<i+1 << "  " << HoverE_smear_low_.at(i) <<endl;


        s << "\n\nH over E smear high Correction Factors "<<endl;
        s <<"ieta  Correction Factor" <<endl;
        for(int i = 0 ; i<28; i++)
            s << setw(4) <<i+1 << "  " << HoverE_smear_high_.at(i) <<endl;

        s << "\n\necal calibrations "<<endl;
        s << "ieta  CorrFactor1  CorrFactor2  CorrFactor3" <<endl;
        int end =  ecal_calib_[0].size();
        for(int i = 0 ; i<28; i++) {
            s << setw(4) << i;
            for(int j = 0; j< end ; j++)
                s << setw(11) << setprecision(8) << ecal_calib_[i][j] ;
         
            s << endl;

        }

        s <<"\n\nhcal calibrations "<<endl;
        s <<"ieta  CorrFactor1  CorrFactor2  CorrFactor3" <<endl;
        end =  hcal_calib_[0].size();
        for(int i = 0 ; i<28; i++) {
            s << setw(4) << i;
            for(int j = 0; j< end ; j++)
                s << setw(11) << setprecision(8) << hcal_calib_[i][j] ;
         
            s << endl;

        }
        s <<"\n\nhcal_high calibrations "<<endl;
        s <<"ieta  CorrFactor1  CorrFactor2  CorrFactor3" <<endl;
        end =  hcal_high_calib_[0].size();
        for(int i = 0 ; i<28; i++) {
            s << setw(4) << i;
            for(int j = 0; j< end ; j++)
                s << setw(11) << setprecision(8) << hcal_high_calib_[i][j] ;
         
            s << endl;

        }
        end = cross_terms_[0].size();
        s <<"\n\ncross terms calibrations "<<endl;
        s <<"ieta  CorrFactor1  CorrFactor2  CorrFactor3  CorrFactor4  CorrFactor5  CorrFactor6" <<endl;
        for(int i = 0 ; i<28; i++) {
            s << setw(4) << i;
            for(int j = 0; j< end ; j++)
                s << setw(11) << setprecision(8) << cross_terms_[i][j] ;
         
            s << endl;

        }
 
    } else
        s<< "\n\nNOT USING calibration variables " <<endl;

    s << "\n\n" <<endl;

}