EcalTPGParamBuilder Class Reference

#include <EcalTPGParamBuilder.h>

Inheritance diagram for EcalTPGParamBuilder:

Public Member Functions
virtual void	analyze (const edm::Event &evt, const edm::EventSetup &evtSetup)
virtual void	beginJob ()
bool	checkIfOK (const EcalPedestals::Item &item)
	EcalTPGParamBuilder (edm::ParameterSet const &pSet)
	~EcalTPGParamBuilder ()
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 (const std::string &iProcessName, std::vector< const char * > &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 ()

Private Member Functions
void	computeFineGrainEBParameters (uint &lowRatio, uint &highRatio, uint &lowThreshold, uint &highThreshold, uint &lut)
void	computeFineGrainEEParameters (uint &threshold, uint &lut_strip, uint &lut_tower)
bool	computeLinearizerParam (double theta, double gainRatio, double calibCoeff, std::string subdet, int &mult, int &shift)
void	computeLUT (int *lut, std::string det="EB")
std::vector< unsigned int >	computeWeights (EcalShapeBase &shape, TH1F *histo)
void	create_header ()
void	getCoeff (coeffStruc &coeff, const EcalIntercalibConstantMap &calibMap, uint rawId)
void	getCoeff (coeffStruc &coeff, const EcalGainRatioMap &gainMap, uint rawId)
void	getCoeff (coeffStruc &coeff, const EcalPedestalsMap &pedMap, uint rawId)
void	getCoeff (coeffStruc &coeff, const std::map< EcalLogicID, MonPedestalsDat > &pedMap, const EcalLogicID &logicId)
std::pair< std::string, int >	getCrate (int tcc)
std::string	getDet (int tcc)
int	getEtaSlice (int tccId, int towerInTCC)
int	getGCTRegionEta (int tteta)
int	getGCTRegionPhi (int ttphi)
bool	realignBaseline (linStruc &lin, float forceBase12)
int	uncodeWeight (double weight, int complement2=7)
double	uncodeWeight (int iweight, int complement2=7)

Private Attributes
int	bst_conf_id_
int	btt_conf_id_
int	bxt_conf_id_
unsigned int	complement2_
EcalTPGDBApp *	db_
bool	DBEE_
int	del_conf_id_
std::map< int, std::vector< int > >	delays_EB_
std::map< int, std::vector< int > >	delays_EE_
double	Et_sat_EB_
double	Et_sat_EE_
edm::ESHandle < EcalTrigTowerConstituentsMap >	eTTmap_
double	FG_highRatio_EB_
double	FG_highThreshold_EB_
double	FG_lowRatio_EB_
double	FG_lowThreshold_EB_
unsigned int	FG_lut_EB_
unsigned int	FG_lut_strip_EE_
unsigned int	FG_lut_tower_EE_
double	FG_Threshold_EE_
int	fgr_conf_id_
int	forcedPedestalValue_
bool	forceEtaSlice_
std::ofstream *	geomFile_
bool	H2_
int	lin_conf_id_
int	lut_conf_id_
double	LUT_constant_EB_
double	LUT_constant_EE_
double	LUT_noise_EB_
double	LUT_noise_EE_
std::string	LUT_option_
double	LUT_stochastic_EB_
double	LUT_stochastic_EE_
double	LUT_threshold_EB_
double	LUT_threshold_EE_
int	m_write_bst
int	m_write_btt
int	m_write_bxt
int	m_write_del
int	m_write_fgr
int	m_write_lin
int	m_write_lut
int	m_write_ped
int	m_write_sli
int	m_write_spi
int	m_write_wei
unsigned int	nSample_
Char_t	ntupleCrate_ [10]
Char_t	ntupleDet_ [10]
Int_t *	ntupleInts_
std::ofstream *	out_file_
int	ped_conf_id_
unsigned int	pedestal_offset_
std::map< int, std::vector< int > >	phases_EB_
std::map< int, std::vector< int > >	phases_EE_
unsigned int	sampleMax_
unsigned int	SFGVB_lut_
int	SFGVB_SpikeKillingThreshold_
unsigned int	SFGVB_Threshold_
int	sli_conf_id_
unsigned int	sliding_
int	spi_conf_id_
std::string	tag_
const CaloSubdetectorGeometry *	theBarrelGeometry_
const CaloSubdetectorGeometry *	theEndcapGeometry_
const EcalElectronicsMapping *	theMapping_
std::string	TimingDelays_EB_
std::string	TimingDelays_EE_
std::string	TimingPhases_EB_
std::string	TimingPhases_EE_
double	TTF_highThreshold_EB_
double	TTF_highThreshold_EE_
double	TTF_lowThreshold_EB_
double	TTF_lowThreshold_EE_
bool	useInterCalibration_
bool	useTransverseEnergy_
int	version_
int	wei_conf_id_
double	weight_timeShift_
bool	weight_unbias_recovery_
bool	writeToDB_
bool	writeToFiles_
double	xtal_LSB_EB_
double	xtal_LSB_EE_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
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 54 of file EcalTPGParamBuilder.h.

Constructor & Destructor Documentation

EcalTPGParamBuilder::EcalTPGParamBuilder ( edm::ParameterSet const &  pSet )

explicit

Definition at line 52 of file EcalTPGParamBuilder.cc.

References bst_conf_id_, btt_conf_id_, bxt_conf_id_, gather_cfg::cout, db_, DBEE_, del_conf_id_, delays_EB_, delays_EE_, alignCSCRings::e, Et_sat_EB_, Et_sat_EE_, cppFunctionSkipper::exception, FG_highRatio_EB_, FG_highThreshold_EB_, FG_lowRatio_EB_, FG_lowThreshold_EB_, FG_lut_EB_, FG_lut_strip_EE_, FG_lut_tower_EE_, FG_Threshold_EE_, fgr_conf_id_, forcedPedestalValue_, forceEtaSlice_, geomFile_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), H2_, lin_conf_id_, lut_conf_id_, LUT_constant_EB_, LUT_constant_EE_, LUT_noise_EB_, LUT_noise_EE_, LUT_option_, LUT_stochastic_EB_, LUT_stochastic_EE_, LUT_threshold_EB_, LUT_threshold_EE_, m_write_bst, m_write_btt, m_write_bxt, m_write_del, m_write_fgr, m_write_lin, m_write_lut, m_write_ped, m_write_sli, m_write_spi, m_write_wei, dbtoconf::out, out_file_, GetRecoTauVFromDQM_MC_cff::outFile, ped_conf_id_, pedestal_offset_, phases_EB_, phases_EE_, sampleMax_, SFGVB_lut_, SFGVB_SpikeKillingThreshold_, SFGVB_Threshold_, funct::sin(), sli_conf_id_, sliding_, spi_conf_id_, AlCaHLTBitMon_QueryRunRegistry::string, tag_, TimingDelays_EB_, TimingDelays_EE_, TimingPhases_EB_, TimingPhases_EE_, TTF_highThreshold_EB_, TTF_highThreshold_EE_, TTF_lowThreshold_EB_, TTF_lowThreshold_EE_, useInterCalibration_, useTransverseEnergy_, version_, wei_conf_id_, weight_timeShift_, weight_unbias_recovery_, writeToDB_, and writeToFiles_.

  : xtal_LSB_EB_(0), xtal_LSB_EE_(0), nSample_(5), complement2_(7)
{
  ped_conf_id_=0;
  lin_conf_id_=0;
  lut_conf_id_=0;
  wei_conf_id_=0;
  fgr_conf_id_=0;
  sli_conf_id_=0;
  spi_conf_id_=0; //modif-alex 21/01/11
  del_conf_id_=0; //modif-alex 21/01/11
  bxt_conf_id_=0;
  btt_conf_id_=0;
  bst_conf_id_=0;
  tag_="";
  version_=0;

  m_write_ped=1;
  m_write_lin=1;
  m_write_lut=1;
  m_write_wei=1;
  m_write_fgr=1;
  m_write_sli=1;
  m_write_spi=1; //modif-alex 21/01/11
  m_write_del=1; //modif-alex 21/01/11
  m_write_bxt=1;
  m_write_btt=1;
  m_write_bst=1;

  writeToDB_  = pSet.getParameter<bool>("writeToDB") ;
  DBEE_ = pSet.getParameter<bool>("allowDBEE") ;
  string DBsid    = pSet.getParameter<std::string>("DBsid") ;
  string DBuser   = pSet.getParameter<std::string>("DBuser") ;
  string DBpass   = pSet.getParameter<std::string>("DBpass") ;
  //uint32_t DBport = pSet.getParameter<unsigned int>("DBport") ;
  
  tag_   = pSet.getParameter<std::string>("TPGtag") ;
  version_ = pSet.getParameter<unsigned int>("TPGversion") ;

  m_write_ped = pSet.getParameter<unsigned int>("TPGWritePed") ;
  m_write_lin = pSet.getParameter<unsigned int>("TPGWriteLin") ;
  m_write_lut = pSet.getParameter<unsigned int>("TPGWriteLut") ;
  m_write_wei = pSet.getParameter<unsigned int>("TPGWriteWei") ;
  m_write_fgr = pSet.getParameter<unsigned int>("TPGWriteFgr") ;
  m_write_sli = pSet.getParameter<unsigned int>("TPGWriteSli") ;
  m_write_spi = pSet.getParameter<unsigned int>("TPGWriteSpi") ; //modif-alex 21/01/11
  m_write_del = pSet.getParameter<unsigned int>("TPGWriteDel") ; //modif-alex 21/01/11
  m_write_bxt = pSet.getParameter<unsigned int>("TPGWriteBxt") ;
  m_write_btt = pSet.getParameter<unsigned int>("TPGWriteBtt") ;
  m_write_bst = pSet.getParameter<unsigned int>("TPGWriteBst") ;

  btt_conf_id_=m_write_btt;
  bxt_conf_id_=m_write_bxt;
  bst_conf_id_=m_write_bst;

  if(m_write_ped != 0 && m_write_ped != 1 ) ped_conf_id_=m_write_ped;
  
  try {
    if (writeToDB_) std::cout << "data will be saved with tag and version="<< tag_<< ".version"<<version_<< endl;
    db_ = new EcalTPGDBApp(DBsid, DBuser, DBpass) ;
  } catch (exception &e) {
    cout << "ERROR:  " << e.what() << endl;
  } catch (...) {
    cout << "Unknown error caught" << endl;
  }

  writeToFiles_ =  pSet.getParameter<bool>("writeToFiles") ;
  if (writeToFiles_) {
    std::string outFile = pSet.getParameter<std::string>("outFile") ;
    out_file_ = new std::ofstream(outFile.c_str(), std::ios::out) ;  
  }
  geomFile_   = new std::ofstream("geomFile.txt", std::ios::out) ;  


  useTransverseEnergy_ = pSet.getParameter<bool>("useTransverseEnergy") ;
  
  Et_sat_EB_ = pSet.getParameter<double>("Et_sat_EB") ;
  Et_sat_EE_ = pSet.getParameter<double>("Et_sat_EE") ;
  sliding_ = pSet.getParameter<unsigned int>("sliding") ;
  weight_timeShift_ = pSet.getParameter<double>("weight_timeShift") ;
  sampleMax_ = pSet.getParameter<unsigned int>("weight_sampleMax") ;
  weight_unbias_recovery_ = pSet.getParameter<bool>("weight_unbias_recovery") ;

  forcedPedestalValue_ = pSet.getParameter<int>("forcedPedestalValue") ;
  forceEtaSlice_ = pSet.getParameter<bool>("forceEtaSlice") ;
    
  LUT_option_ = pSet.getParameter<std::string>("LUT_option") ;
  LUT_threshold_EB_ = pSet.getParameter<double>("LUT_threshold_EB") ;
  LUT_threshold_EE_ = pSet.getParameter<double>("LUT_threshold_EE") ;
  LUT_stochastic_EB_ = pSet.getParameter<double>("LUT_stochastic_EB") ;
  LUT_noise_EB_ =pSet.getParameter<double>("LUT_noise_EB") ;
  LUT_constant_EB_ =pSet.getParameter<double>("LUT_constant_EB") ;
  LUT_stochastic_EE_ = pSet.getParameter<double>("LUT_stochastic_EE") ;
  LUT_noise_EE_ =pSet.getParameter<double>("LUT_noise_EE") ;
  LUT_constant_EE_ =pSet.getParameter<double>("LUT_constant_EE") ;

  TTF_lowThreshold_EB_ = pSet.getParameter<double>("TTF_lowThreshold_EB") ;
  TTF_highThreshold_EB_ = pSet.getParameter<double>("TTF_highThreshold_EB") ;
  TTF_lowThreshold_EE_ = pSet.getParameter<double>("TTF_lowThreshold_EE") ;
  TTF_highThreshold_EE_ = pSet.getParameter<double>("TTF_highThreshold_EE") ;

  FG_lowThreshold_EB_ = pSet.getParameter<double>("FG_lowThreshold_EB") ;
  FG_highThreshold_EB_ = pSet.getParameter<double>("FG_highThreshold_EB") ;
  FG_lowRatio_EB_ = pSet.getParameter<double>("FG_lowRatio_EB") ;
  FG_highRatio_EB_ = pSet.getParameter<double>("FG_highRatio_EB") ;
  FG_lut_EB_ = pSet.getParameter<unsigned int>("FG_lut_EB") ;
  FG_Threshold_EE_ = pSet.getParameter<double>("FG_Threshold_EE") ;
  FG_lut_strip_EE_ = pSet.getParameter<unsigned int>("FG_lut_strip_EE") ;
  FG_lut_tower_EE_ = pSet.getParameter<unsigned int>("FG_lut_tower_EE") ;
  SFGVB_Threshold_ = pSet.getParameter<unsigned int>("SFGVB_Threshold") ;
  SFGVB_lut_ = pSet.getParameter<unsigned int>("SFGVB_lut") ;
  SFGVB_SpikeKillingThreshold_ = pSet.getParameter<int>("SFGVB_SpikeKillingThreshold") ; //modif-alex 21/01/11
  pedestal_offset_ = pSet.getParameter<unsigned int>("pedestal_offset") ;

  useInterCalibration_  = pSet.getParameter<bool>("useInterCalibration") ;
  H2_ = pSet.getUntrackedParameter<bool>("H2",false) ;

  //modif-alex-23/02/2011
  //convert the spike killing first from GeV to ADC (10 bits)                                                                                                                                                    
  //depending on the saturation scale: Et_sat_EB_                                                                                                                                                                
  if(SFGVB_SpikeKillingThreshold_ == -1 || (SFGVB_SpikeKillingThreshold_ > Et_sat_EB_)) 
    SFGVB_SpikeKillingThreshold_ = 1023; //nokilling
  else
    SFGVB_SpikeKillingThreshold_ = int(SFGVB_SpikeKillingThreshold_ * 1024/Et_sat_EB_);    
  std::cout << "INFO:SPIKE KILLING THRESHOLD (ADC)=" << SFGVB_SpikeKillingThreshold_ << std::endl;

  //modif-alex-02/02/11
  //TIMING information 
  TimingDelays_EB_ = pSet.getParameter<std::string>("timing_delays_EB") ;
  TimingDelays_EE_ = pSet.getParameter<std::string>("timing_delays_EE") ;
  TimingPhases_EB_ = pSet.getParameter<std::string>("timing_phases_EB") ;
  TimingPhases_EE_ = pSet.getParameter<std::string>("timing_phases_EE") ;

  std::cout << "INFO: READING timing files" << std::endl;
  std::ifstream delay_eb(TimingDelays_EB_.c_str()); 
  if(!delay_eb) std::cout << "ERROR: File " << TimingDelays_EB_.c_str() << " could not be opened" << std::endl;
  std::ifstream delay_ee(TimingDelays_EE_.c_str());
  if(!delay_ee) std::cout << "ERROR: File " << TimingDelays_EE_.c_str() << " could not be opened" << std::endl;
  std::ifstream phase_eb(TimingPhases_EB_.c_str());
  if(!phase_eb) std::cout << "ERROR: File " << TimingPhases_EB_.c_str() << " could not be opened" << std::endl;
  std::ifstream phase_ee(TimingPhases_EE_.c_str());
  if(!phase_ee) std::cout << "ERROR: File " << TimingPhases_EE_.c_str() << " could not be opened" << std::endl;
  
  char buf[1024];
  //READING DELAYS EB
  delay_eb.getline(buf,sizeof(buf),'\n');
  while( delay_eb ) {
    std::stringstream sin(buf);

    int tcc; 
    sin >> tcc;

    vector<int> vec_delays_eb;
    for(int ieb=0; ieb<68; ++ieb){
      int time_delay = -1;
      sin >> time_delay;
      vec_delays_eb.push_back(time_delay);
      if(time_delay==-1)  std::cout << "ERROR:Barrel timing delay -1, check file" << std::endl;
    }

    if(vec_delays_eb.size()!=68)
      std::cout << "ERROR:Barrel timing delay wrong, not enough towers, check file" << std::endl;

    if (delays_EB_.find(tcc) == delays_EB_.end())
      delays_EB_.insert(make_pair(tcc,vec_delays_eb));

    cout << tcc << " "; 
    for(unsigned int ieb=0; ieb<vec_delays_eb.size(); ++ieb)
      cout << vec_delays_eb[ieb] << " ";
    cout << endl;

    delay_eb.getline(buf,sizeof(buf),'\n');
  }//loop delay file EB
  delay_eb.close();

  //READING PHASES EB
  phase_eb.getline(buf,sizeof(buf),'\n');
  while( phase_eb ) {
    std::stringstream sin(buf);
    int tcc;
    sin >> tcc;

    vector<int> vec_phases_eb;
    for(unsigned int ieb=0; ieb<68; ++ieb){
      int time_phase = -1;
      sin >> time_phase;
      vec_phases_eb.push_back(time_phase);
      if(time_phase==-1)  std::cout << "ERROR:Barrel timing phase -1, check file" << std::endl;
    }

    if(vec_phases_eb.size()!=68)
      std::cout << "ERROR:Barrel timing phase wrong, not enough towers, check file" << std::endl;

    if (phases_EB_.find(tcc) == phases_EB_.end())
      phases_EB_.insert(make_pair(tcc,vec_phases_eb));

    cout << tcc << " ";
    for(unsigned int ieb=0; ieb<vec_phases_eb.size(); ++ieb)
      cout << vec_phases_eb[ieb] << " ";
    cout << endl;

    phase_eb.getline(buf,sizeof(buf),'\n');
  }//loop phase file EB
  phase_eb.close();

  //READING DELAYS EE//------------------------------------------------
  delay_ee.getline(buf,sizeof(buf),'\n');
  while( delay_ee ) {
    std::stringstream sin(buf);
    int tcc; 
    sin >> tcc;

    vector<int> vec_delays_ee;
    for(unsigned int iee=0; iee<48; ++iee){
      int time_delay = -1;
      sin >> time_delay;
      vec_delays_ee.push_back(time_delay);
      if(time_delay==-1)  std::cout << "ERROR:EE timing delay -1, check file" << std::endl;
    }

    if(vec_delays_ee.size()!=48)
      std::cout << "ERROR:EE timing delay wrong, not enough towers, check file" << std::endl;

    if (delays_EE_.find(tcc) == delays_EE_.end())
      delays_EE_.insert(make_pair(tcc,vec_delays_ee));

    cout << tcc << " "; 
    for(unsigned int iee=0; iee<vec_delays_ee.size(); ++iee)
      cout << vec_delays_ee[iee] << " ";
    cout << endl;

    delay_ee.getline(buf,sizeof(buf),'\n');
  }//loop delay file EE
  delay_ee.close();


  //READING PHASES EE
  phase_ee.getline(buf,sizeof(buf),'\n');
  while( phase_ee ) {
    std::stringstream sin(buf);
    int tcc;
    sin >> tcc;

    vector<int> vec_phases_ee;
    for(unsigned int iee=0; iee<48; ++iee){
      int time_phase = -1;
      sin >> time_phase;
      vec_phases_ee.push_back(time_phase);
      if(time_phase==-1)  std::cout << "ERROR:EE timing phase -1, check file" << std::endl;
    }

    if(vec_phases_ee.size()!=48)
      std::cout << "ERROR:EE timing phase wrong, not enough towers, check file" << std::endl;

    if (phases_EE_.find(tcc) == phases_EE_.end())
      phases_EE_.insert(make_pair(tcc,vec_phases_ee));

    cout << tcc << " ";
    for(unsigned int iee=0; iee<vec_phases_ee.size(); ++iee)
      cout << vec_phases_ee[iee] << " ";
    cout << endl;

    phase_ee.getline(buf,sizeof(buf),'\n');
  }//loop phase file EE
  phase_ee.close();

  std::cout << "INFO: DONE reading timing files for EB and EE" << std::endl;

}

EcalTPGParamBuilder::~EcalTPGParamBuilder

(

)

Definition at line 322 of file EcalTPGParamBuilder.cc.

References out_file_, and writeToFiles_.

{ 
  if (writeToFiles_) {
    (*out_file_ )<<"EOF"<<std::endl ;
    out_file_->close() ;
    delete out_file_ ;
  }
}

Member Function Documentation

void EcalTPGParamBuilder::analyze ( const edm::Event &  evt, const edm::EventSetup &  evtSetup  )

virtual

Implements edm::EDAnalyzer.

Definition at line 359 of file EcalTPGParamBuilder.cc.

References funct::abs(), newFWLiteAna::base, bst_conf_id_, btt_conf_id_, bxt_conf_id_, coeffStruc::calibCoeff_, EcalTriggerElectronicsId::channelId(), computeFineGrainEBParameters(), computeFineGrainEEParameters(), computeLinearizerParam(), computeLUT(), computeWeights(), gather_cfg::cout, dataset::dataset, db_, DBEE_, EcalElectronicsId::dccId(), EcalElectronicsMapping::DCCid(), TauDecayModes::dec, del_conf_id_, delays_EB_, delays_EE_, alignCSCRings::e, DetId::Ecal, EcalBarrel, EcalEndcap, Et_sat_EB_, Et_sat_EE_, PV3DBase< T, PVType, FrameType >::eta(), eTTmap_, cppFunctionSkipper::exception, EcalCondDBInterface::fetchDataSet(), FG_highRatio_EB_, FG_highThreshold_EB_, FG_lowRatio_EB_, FG_lowThreshold_EB_, FG_Threshold_EE_, fgr_conf_id_, plotBeamSpotDB::first, forcedPedestalValue_, forceEtaSlice_, g, callgraph::G, EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), coeffStruc::gainRatio_, edm::EventSetup::get(), FEConfigMainInfo::getBstId(), FEConfigMainInfo::getBttId(), FEConfigMainInfo::getBxtId(), getCoeff(), getCrate(), getDet(), EcalCondDBInterface::getEcalLogicID(), EcalCondDBInterface::getEcalLogicIDSetOrdered(), EcalElectronicsMapping::getElectronicsId(), FEConfigMainInfo::getFgrId(), getGCTRegionEta(), getGCTRegionPhi(), CaloSubdetectorGeometry::getGeometry(), EcalLogicID::getID1(), EcalLogicID::getID2(), EcalLogicID::getID3(), FEConfigMainInfo::getLinId(), EcalLogicID::getLogicID(), FEConfigMainInfo::getLUTId(), EcalCondObjectContainer< T >::getMap(), EcalLogicID::getName(), FEConfigMainInfo::getPedId(), FEConfigPedDat::getPedMeanG1(), FEConfigPedDat::getPedMeanG12(), FEConfigPedDat::getPedMeanG6(), CaloCellGeometry::getPosition(), FEConfigMainInfo::getSliId(), FEConfigMainInfo::getSpiId(), FEConfigMainInfo::getTimId(), EcalElectronicsMapping::getTriggerElectronicsId(), CaloSubdetectorGeometry::getValidDetIds(), FEConfigMainInfo::getWeiId(), H2_, ecalpyutils::hashedIndex(), i, EcalTrigTowerDetId::ieta(), EcalTrigTowerDetId::ietaAbs(), EcalCondObjectContainer< T >::insert(), EcalTrigTowerDetId::iphi(), EcalElectronicsMapping::iTT(), j, relval_steps::k, lin_conf_id_, lumiPlot::lut, lut_conf_id_, LUT_option_, lumiCalc2::lute, m_write_del, m_write_fgr, m_write_lin, m_write_lut, m_write_ped, m_write_sli, m_write_spi, m_write_wei, EcalPedestal::mean_x1, EcalPedestal::mean_x12, EcalPedestal::mean_x6, VarParsing::mult, linStruc::mult_, ntupleCrate_, ntupleDet_, ntupleInts_, EcalLogicID::NULLID, convertSQLiteXML::ok, dbtoconf::out, out_file_, AlCaHLTBitMon_ParallelJobs::p, ped_conf_id_, linStruc::pedestal_, pedestal_offset_, coeffStruc::pedestals_, phases_EB_, phases_EE_, PV3DBase< T, PVType, FrameType >::phi(), funct::pow(), edm::ESHandle< class >::product(), EcalTriggerElectronicsId::pseudoStripId(), EcalTriggerElectronicsId::rawId(), DetId::rawId(), realignBaseline(), EcalPedestal::rms_x1, EcalPedestal::rms_x12, EcalPedestal::rms_x6, compare_using_db::sample, IODConfig::setConfigTag(), FEConfigLinParamDat::setETSat(), FEConfigLUTParamDat::setETSat(), FEConfigFgrParamDat::setFGhighratio(), FEConfigFgrParamDat::setFGhighthresh(), FEConfigFgrParamDat::setFGlowratio(), FEConfigFgrParamDat::setFGlowthresh(), FEConfigFgrGroupDat::setFgrGroupId(), FEConfigFgrDat::setFgrGroupId(), FEConfigPedInfo::setId(), FEConfigFgrEEStripDat::setLutFgr(), FEConfigLUTGroupDat::setLUTGroupId(), FEConfigLUTDat::setLUTGroupId(), FEConfigFgrEETowerDat::setLutValue(), FEConfigLUTGroupDat::setLUTValue(), FEConfigFgrGroupDat::setLUTValue(), FEConfigLinDat::setMultX1(), FEConfigLinDat::setMultX12(), FEConfigLinDat::setMultX6(), FEConfigPedDat::setPedMeanG1(), FEConfigPedDat::setPedMeanG12(), FEConfigPedDat::setPedMeanG6(), FEConfigFgrGroupDat::setRatioHigh(), FEConfigFgrGroupDat::setRatioLow(), FEConfigLinDat::setShift1(), FEConfigLinDat::setShift12(), FEConfigLinDat::setShift6(), FEConfigSlidingDat::setSliding(), FEConfigSpikeDat::setSpikeThreshold(), FEConfigFgrGroupDat::setThreshHigh(), FEConfigFgrGroupDat::setThreshLow(), FEConfigFgrEEStripDat::setThreshold(), FEConfigTimingDat::setTimingPar1(), FEConfigTimingDat::setTimingPar2(), FEConfigLUTParamDat::setTTThreshhigh(), FEConfigLUTParamDat::setTTThreshlow(), FEConfigWeightGroupDat::setWeight0(), FEConfigWeightGroupDat::setWeight1(), FEConfigWeightGroupDat::setWeight2(), FEConfigWeightGroupDat::setWeight3(), FEConfigWeightGroupDat::setWeight4(), FEConfigWeightDat::setWeightGroupId(), FEConfigWeightGroupDat::setWeightGroupId(), SFGVB_lut_, SFGVB_SpikeKillingThreshold_, SFGVB_Threshold_, edm::shift, linStruc::shift_, funct::sin(), findQualityFiles::size, sli_conf_id_, sliding_, EBDetId::SMCRYSTALMODE, spi_conf_id_, AlCaHLTBitMon_QueryRunRegistry::string, EcalElectronicsId::stripId(), EcalTrigTowerDetId::subDet(), tag_, EcalElectronicsMapping::TCCid(), theBarrelGeometry_, theEndcapGeometry_, theMapping_, PV3DBase< T, PVType, FrameType >::theta(), theta(), dtDQMClient_cfg::threshold, EcalElectronicsId::towerId(), funct::true, TTF_highThreshold_EB_, TTF_highThreshold_EE_, TTF_lowThreshold_EB_, TTF_lowThreshold_EE_, useTransverseEnergy_, version_, wei_conf_id_, create_public_pileup_plots::weights, EcalTPGDBApp::writeToConfDB_Delay(), EcalTPGDBApp::writeToConfDB_Spike(), EcalTPGDBApp::writeToConfDB_TPGFgr(), EcalTPGDBApp::writeToConfDB_TPGLinearCoef(), EcalTPGDBApp::writeToConfDB_TPGLUT(), EcalTPGDBApp::writeToConfDB_TPGMain(), EcalTPGDBApp::writeToConfDB_TPGPedestals(), EcalTPGDBApp::writeToConfDB_TPGSliding(), EcalTPGDBApp::writeToConfDB_TPGWeight(), writeToDB_, writeToFiles_, x, xtal_LSB_EB_, xtal_LSB_EE_, EcalElectronicsId::xtalId(), EEDetId::XYMODE, detailsBasic3DVector::y, detailsBasic3DVector::z, EcalTrigTowerDetId::zside(), and ecaldqm::zside().

{
  using namespace edm;
  using namespace std;

  // geometry
  ESHandle<CaloGeometry> theGeometry;
  ESHandle<CaloSubdetectorGeometry> theEndcapGeometry_handle, theBarrelGeometry_handle;
  evtSetup.get<CaloGeometryRecord>().get( theGeometry );
  evtSetup.get<EcalEndcapGeometryRecord>().get("EcalEndcap",theEndcapGeometry_handle);
  evtSetup.get<EcalBarrelGeometryRecord>().get("EcalBarrel",theBarrelGeometry_handle);
  evtSetup.get<IdealGeometryRecord>().get(eTTmap_);
  theEndcapGeometry_ = &(*theEndcapGeometry_handle);
  theBarrelGeometry_ = &(*theBarrelGeometry_handle);

  // electronics mapping
  ESHandle< EcalElectronicsMapping > ecalmapping;
  evtSetup.get< EcalMappingRcd >().get(ecalmapping);
  theMapping_ = ecalmapping.product();

  
  // histo
  TFile saving ("EcalTPGParam.root","recreate") ;
  saving.cd () ;
  TH2F * ICEB = new TH2F("ICEB", "IC: Barrel", 360, 1, 361, 172, -86, 86) ;
  ICEB->GetYaxis()->SetTitle("eta index") ;
  ICEB->GetXaxis()->SetTitle("phi index") ;  
  TH2F * tpgFactorEB = new TH2F("tpgFactorEB", "tpgFactor: Barrel", 360, 1, 361, 172, -86, 86) ;
  tpgFactorEB->GetYaxis()->SetTitle("eta index") ;
  tpgFactorEB->GetXaxis()->SetTitle("phi index") ;  

  TH2F * ICEEPlus = new TH2F("ICEEPlus", "IC: Plus Endcap", 120, -9, 111, 120, -9, 111) ;
  ICEEPlus->GetYaxis()->SetTitle("y index") ;
  ICEEPlus->GetXaxis()->SetTitle("x index") ;
  TH2F * tpgFactorEEPlus = new TH2F("tpgFactorEEPlus", "tpgFactor: Plus Endcap", 120, -9, 111, 120, -9, 111) ;
  tpgFactorEEPlus->GetYaxis()->SetTitle("y index") ;
  tpgFactorEEPlus->GetXaxis()->SetTitle("x index") ;
  TH2F * ICEEMinus = new TH2F("ICEEMinus", "IC: Minus Endcap", 120, -9, 111, 120, -9, 111) ;
  ICEEMinus->GetYaxis()->SetTitle("y index") ;
  ICEEMinus->GetXaxis()->SetTitle("x index") ;
  TH2F * tpgFactorEEMinus = new TH2F("tpgFactorEEMinus", "tpgFactor: Minus Endcap", 120, -9, 111, 120, -9, 111) ;
  tpgFactorEEMinus->GetYaxis()->SetTitle("y index") ;
  tpgFactorEEMinus->GetXaxis()->SetTitle("x index") ;

  TH2F * IC = new TH2F("IC", "IC", 720, -acos(-1.), acos(-1.), 600, -3., 3.) ;
  IC->GetYaxis()->SetTitle("eta") ;
  IC->GetXaxis()->SetTitle("phi") ;  
  TH2F * tpgFactor = new TH2F("tpgFactor", "tpgFactor", 720, -acos(-1.), acos(-1.), 600, -3., 3.) ;
  tpgFactor->GetYaxis()->SetTitle("eta") ;
  tpgFactor->GetXaxis()->SetTitle("phi") ;  

  TH1F * hshapeEB = new TH1F("shapeEB", "shapeEB", 250, 0., 10.) ;
  TH1F * hshapeEE = new TH1F("shapeEE", "shapeEE", 250, 0., 10.) ;

  TTree * ntuple = new TTree("tpgmap","TPG geometry map") ;
  const std::string branchFloat[26] = {"fed","tcc","tower","stripInTower","xtalInStrip",
                       "CCU","VFE","xtalInVFE","xtalInCCU","ieta","iphi",
                       "ix","iy","iz","hashedId","ic","cmsswId","dbId","ietaTT","iphiTT",
                       "TCCch","TCCslot","SLBch","SLBslot","ietaGCT","iphiGCT"} ;
  ntupleInts_ = new Int_t[26] ;
  for (int i=0 ; i<26 ; i++) ntuple->Branch(branchFloat[i].c_str(),&ntupleInts_[i],(branchFloat[i]+string("/I")).c_str()) ;
  ntuple->Branch("det",ntupleDet_,"det/C") ;
  ntuple->Branch("crate",ntupleCrate_,"crate/C") ;


  TNtuple *ntupleSpike = new TNtuple("spikeParam","Spike parameters","gainId:theta:G:g:ped:pedLin") ;




  // Initialization section //
  list<uint32_t> towerListEB ;
  list<uint32_t> stripListEB ;
  list<uint32_t> towerListEE ;
  list<uint32_t> stripListEE ;
  list<uint32_t>::iterator itList ;
  
  map<int, uint32_t> stripMapEB ; // <EcalLogicId.hashed, strip elec id>
  map<uint32_t, uint32_t> stripMapEBsintheta ; // <strip elec id, sintheta>


  // Pedestals

  EcalPedestalsMap  pedMap ;

  if(m_write_ped == 1) {
    std::cout <<"Getting the pedestals from offline DB..."<<endl;

    ESHandle<EcalPedestals> pedHandle;
    evtSetup.get<EcalPedestalsRcd>().get( pedHandle );
    pedMap = pedHandle.product()->getMap() ;


    EcalPedestalsMapIterator pedIter ; 
    int nPed = 0 ;
    for (pedIter = pedMap.begin() ; pedIter != pedMap.end() && nPed<10 ; ++pedIter, nPed++) {
      EcalPedestals::Item aped = (*pedIter);
      std::cout<<aped.mean_x12<<", "<<aped.mean_x6<<", "<<aped.mean_x1<<std::endl ;
    }
  } else if(m_write_ped==0) {
    std::cout <<"Getting the pedestals from previous configuration"<<std::endl;
    
    EcalPedestals peds ;

    FEConfigMainInfo fe_main_info;
    fe_main_info.setConfigTag(tag_);
    try {
      std::cout << "trying to read previous tag if it exists tag="<< tag_<< ".version"<<version_<< endl;
      db_-> fetchConfigSet(&fe_main_info);
      if(fe_main_info.getPedId()>0 ) ped_conf_id_=fe_main_info.getPedId();
      
      FEConfigPedInfo fe_ped_info;
      fe_ped_info.setId(ped_conf_id_);
      db_-> fetchConfigSet(&fe_ped_info);
      std::map<EcalLogicID, FEConfigPedDat> dataset_TpgPed;
      db_->fetchDataSet(&dataset_TpgPed, &fe_ped_info);
      
      typedef std::map<EcalLogicID, FEConfigPedDat>::const_iterator CIfeped;
      EcalLogicID ecid_xt;
      FEConfigPedDat  rd_ped;
      int icells=0;
      for (CIfeped p = dataset_TpgPed.begin(); p != dataset_TpgPed.end(); p++) 
    {
      ecid_xt = p->first;
      rd_ped  = p->second;
      
      std::string ecid_name=ecid_xt.getName();
      
      // EB data
      if (ecid_name=="EB_crystal_number") {


        int sm_num=ecid_xt.getID1();
        int xt_num=ecid_xt.getID2();
        
        EBDetId ebdetid(sm_num,xt_num,EBDetId::SMCRYSTALMODE);
        EcalPedestals::Item item;
        item.mean_x1  =rd_ped.getPedMeanG1() ;
        item.mean_x6  =rd_ped.getPedMeanG6();
        item.mean_x12 =rd_ped.getPedMeanG12();
        item.rms_x1  =0.5 ;
        item.rms_x6  =1. ;
        item.rms_x12  =1.2 ;

        if(icells<10) std::cout << " copy the EB data " << " ped = "  << item.mean_x12<< std::endl;

        peds.insert(std::make_pair(ebdetid.rawId(),item));

        ++icells;
      }else if (ecid_name=="EE_crystal_number"){
              
        // EE data
        int z=ecid_xt.getID1();
        int x=ecid_xt.getID2();
        int y=ecid_xt.getID3();
        EEDetId eedetid(x,y,z,EEDetId::XYMODE);
        EcalPedestals::Item item;
        item.mean_x1  =rd_ped.getPedMeanG1() ;
        item.mean_x6  =rd_ped.getPedMeanG6();
        item.mean_x12 =rd_ped.getPedMeanG12();
        item.rms_x1  =0.5 ;
        item.rms_x6  =1. ;
        item.rms_x12  =1.2 ;
        
        peds.insert(std::make_pair(eedetid.rawId(),item));
        ++icells;
      }
    }
      
      pedMap = peds.getMap() ;
      
      EcalPedestalsMapIterator pedIter ; 
      int nPed = 0 ;
      for (pedIter = pedMap.begin() ; pedIter != pedMap.end() && nPed<10 ; ++pedIter, nPed++) {
    EcalPedestals::Item aped = (*pedIter);
    std::cout<<aped.mean_x12<<", "<<aped.mean_x6<<", "<<aped.mean_x1<<std::endl ;
      }
      
      
    } catch (exception &e) {
      cout << " error reading previous pedestals " << endl;
    } catch (...) {
      cout << "Unknown error reading previous pedestals " << endl;
    }
    

      

  } else if(m_write_ped >1) {
    std::cout <<"Getting the pedestals from configuration number"<< m_write_ped <<std::endl;



    EcalPedestals peds ;

      try {

      FEConfigPedInfo fe_ped_info;
      fe_ped_info.setId(m_write_ped);
      db_-> fetchConfigSet(&fe_ped_info);
      std::map<EcalLogicID, FEConfigPedDat> dataset_TpgPed;
      db_->fetchDataSet(&dataset_TpgPed, &fe_ped_info);
        
      typedef std::map<EcalLogicID, FEConfigPedDat>::const_iterator CIfeped;
      EcalLogicID ecid_xt;
      FEConfigPedDat  rd_ped;
      int icells=0;
      for (CIfeped p = dataset_TpgPed.begin(); p != dataset_TpgPed.end(); p++) 
        {
          ecid_xt = p->first;
          rd_ped  = p->second;
              
          std::string ecid_name=ecid_xt.getName();
              
          // EB data
          if (ecid_name=="EB_crystal_number") {
        if(icells<10) std::cout << " copy the EB data " << " icells = " << icells << std::endl;
        int sm_num=ecid_xt.getID1();
        int xt_num=ecid_xt.getID2();
              
        EBDetId ebdetid(sm_num,xt_num,EBDetId::SMCRYSTALMODE);
        EcalPedestals::Item item;
        item.mean_x1  =(unsigned int)rd_ped.getPedMeanG1() ;
        item.mean_x6  =(unsigned int)rd_ped.getPedMeanG6();
        item.mean_x12 =(unsigned int)rd_ped.getPedMeanG12();
        item.rms_x1  =0.5 ;
        item.rms_x6  =1. ;
        item.rms_x12  =1.2 ;
              
        peds.insert(std::make_pair(ebdetid.rawId(),item));
        ++icells;
          }else if (ecid_name=="EE_crystal_number"){
              
        // EE data
        int z=ecid_xt.getID1();
        int x=ecid_xt.getID2();
        int y=ecid_xt.getID3();
        EEDetId eedetid(x,y,z,EEDetId::XYMODE);
        EcalPedestals::Item item;
        item.mean_x1  =(unsigned int)rd_ped.getPedMeanG1();
        item.mean_x6  =(unsigned int)rd_ped.getPedMeanG6();
        item.mean_x12 =(unsigned int)rd_ped.getPedMeanG12();
        item.rms_x1  =0.5 ;
        item.rms_x6  =1. ;
        item.rms_x12  =1.2 ;
              
        peds.insert(std::make_pair(eedetid.rawId(),item));
        ++icells;
          }
        }

      pedMap = peds.getMap() ;

      EcalPedestalsMapIterator pedIter ; 
      int nPed = 0 ;
      for (pedIter = pedMap.begin() ; pedIter != pedMap.end() && nPed<10 ; ++pedIter, nPed++) {
        EcalPedestals::Item aped = (*pedIter);
        std::cout<<aped.mean_x12<<", "<<aped.mean_x6<<", "<<aped.mean_x1<<std::endl ;
      }
      

      } catch (exception &e) {
    cout << " error reading previous pedestals " << endl;
      } catch (...) {
    cout << "Unknown error reading previous pedestals " << endl;
      }


  }



  std::cout<<"...\n"<<std::endl ; 


  // Intercalib constants
  std::cout <<"Getting intercalib from offline DB..."<<endl;
  ESHandle<EcalIntercalibConstants> pIntercalib ;
  evtSetup.get<EcalIntercalibConstantsRcd>().get(pIntercalib) ;
  const EcalIntercalibConstants * intercalib = pIntercalib.product() ;
  const EcalIntercalibConstantMap & calibMap = intercalib->getMap() ;
  EcalIntercalibConstantMap::const_iterator calIter ;
  int nCal = 0 ;
  for (calIter = calibMap.begin() ; calIter != calibMap.end() && nCal<10 ; ++calIter, nCal++) {
    std::cout<<(*calIter)<<std::endl ;
  }  
  std::cout<<"...\n"<<std::endl ;
  float calibvec[1700] ;
  if (H2_) {
    std::cout<<"H2: overwriting IC coef with file"<<std::endl ;
    std::ifstream calibfile("calib_sm36.txt", std::ios::out) ;  
    int idata, icry ; 
    float fdata, fcali ;
    std::string strdata ;
    if (calibfile.is_open()) {
      calibfile >> strdata >>  strdata >>  strdata >>  strdata >>  strdata ;
      while(!calibfile.eof()) {
    calibfile >> idata >> icry >> fcali >> fdata >> fdata ;
    calibvec[icry-1] = fcali ;
    if(calibfile.eof()){
      break ; // avoid last line duplication
    }
      }
    }
  }    

  // Gain Ratios
  std::cout <<"Getting the gain ratios from offline DB..."<<endl;
  ESHandle<EcalGainRatios> pRatio;
  evtSetup.get<EcalGainRatiosRcd>().get(pRatio);
  const EcalGainRatioMap & gainMap = pRatio.product()->getMap();
  EcalGainRatioMap::const_iterator gainIter ;
  int nGain = 0 ;
  for (gainIter = gainMap.begin() ; gainIter != gainMap.end() && nGain<10 ; ++gainIter, nGain++) {
    const EcalMGPAGainRatio & aGain = (*gainIter) ;
    std::cout<<aGain.gain12Over6()<<", "<<aGain.gain6Over1() * aGain.gain12Over6()<<std::endl ;
  }  
  std::cout<<"...\n"<<std::endl ;    


  // ADCtoGeV
  std::cout <<"Getting the ADC to GEV from offline DB..."<<endl;
  ESHandle<EcalADCToGeVConstant> pADCToGeV ;
  evtSetup.get<EcalADCToGeVConstantRcd>().get(pADCToGeV) ;
  const EcalADCToGeVConstant * ADCToGeV = pADCToGeV.product() ;
  xtal_LSB_EB_ = ADCToGeV->getEBValue() ;
  xtal_LSB_EE_ = ADCToGeV->getEEValue() ;
  std::cout<<"xtal_LSB_EB_ = "<<xtal_LSB_EB_<<std::endl ;
  std::cout<<"xtal_LSB_EE_ = "<<xtal_LSB_EE_<<std::endl ;
  std::cout<<std::endl ;  

  
  vector<EcalLogicID> my_EcalLogicId;
  vector<EcalLogicID> my_TTEcalLogicId;
  vector<EcalLogicID> my_StripEcalLogicId;
  EcalLogicID my_EcalLogicId_EB;
    // Endcap identifiers
  EcalLogicID my_EcalLogicId_EE;
  vector<EcalLogicID> my_TTEcalLogicId_EE;
  vector<EcalLogicID> my_RTEcalLogicId_EE;
  vector<EcalLogicID> my_StripEcalLogicId1_EE;
  vector<EcalLogicID> my_StripEcalLogicId2_EE;
  vector<EcalLogicID> my_CrystalEcalLogicId_EE;
  vector<EcalLogicID> my_TTEcalLogicId_EB_by_TCC;
  vector<EcalLogicID> my_StripEcalLogicId_EE_strips_by_TCC;


  std::cout<<"going to get the ecal logic id set"<< endl;
  my_EcalLogicId_EB = db_->getEcalLogicID( "EB",EcalLogicID::NULLID,EcalLogicID::NULLID,EcalLogicID::NULLID,"EB");
  my_EcalLogicId_EE = db_->getEcalLogicID( "EE",EcalLogicID::NULLID,EcalLogicID::NULLID,EcalLogicID::NULLID,"EE");
  my_EcalLogicId = db_->getEcalLogicIDSetOrdered( "EB_crystal_number",
                          1, 36,
                          1, 1700,
                          EcalLogicID::NULLID,EcalLogicID::NULLID,
                          "EB_crystal_number",12 );
  my_TTEcalLogicId = db_->getEcalLogicIDSetOrdered( "EB_trigger_tower",
                            1, 36,
                            1, 68,
                            EcalLogicID::NULLID,EcalLogicID::NULLID,
                            "EB_trigger_tower",12 );
  my_StripEcalLogicId = db_->getEcalLogicIDSetOrdered( "EB_VFE",   1, 36,   1, 68,   1,5 ,  "EB_VFE",123 ); //last digi means ordered 1st by SM,then TT, then strip 
  std::cout<<"got the 3 ecal barrel logic id set"<< endl;

  // EE crystals identifiers
  my_CrystalEcalLogicId_EE = db_->getEcalLogicIDSetOrdered("EE_crystal_number",
                               -1, 1,
                               0, 200,
                               0, 200,
                               "EE_crystal_number",123 );
  
  // EE Strip identifiers
  // DCC=601-609 TT = ~40 EEstrip = 5
  my_StripEcalLogicId1_EE = db_->getEcalLogicIDSetOrdered( "ECAL_readout_strip",   
                               601, 609,   
                               1, 100,   
                               0,5 ,  
                               "ECAL_readout_strip",123 );
  // EE Strip identifiers
  // DCC=646-654 TT = ~40 EEstrip = 5
  my_StripEcalLogicId2_EE = db_->getEcalLogicIDSetOrdered( "ECAL_readout_strip",   
                               646, 654,   
                               1, 100,   
                               0,5 ,  
                               "ECAL_readout_strip",123 );
  // ----> modif here 31/1/2011   
  // EE Strip identifiers by TCC tower strip
  // TCC=1-108 TT = ~40 EEstrip = 1-5
  my_StripEcalLogicId_EE_strips_by_TCC = db_->getEcalLogicIDSetOrdered( "EE_trigger_strip",   
                               1, 108,   
                               1, 100,   
                               1,5 ,  
                               "EE_trigger_strip",123 );

  
  // ----> modif here 31/1/2011
  // ECAL Barrel trigger towers by TCC/tower identifiers
  // TTC=38-72 TT = 1-68
  my_TTEcalLogicId_EB_by_TCC = db_->getEcalLogicIDSetOrdered( "ECAL_trigger_tower",
                               1, 108,
                               1, 68,
                               EcalLogicID::NULLID,EcalLogicID::NULLID,
                               "ECAL_trigger_tower",12 );
  
  // EE TT identifiers
  // TTC=72 TT = 1440
  my_TTEcalLogicId_EE = db_->getEcalLogicIDSetOrdered( "EE_trigger_tower",
                               1, 108,
                               1, 40,
                               EcalLogicID::NULLID,EcalLogicID::NULLID,
                               "EE_trigger_tower",12 );
  
  // EE TT identifiers    
  // TTC=72 TT = 1440
  my_RTEcalLogicId_EE = db_->getEcalLogicIDSetOrdered( "EE_readout_tower",
                               1, 1000,
                               1, 100,
                               EcalLogicID::NULLID,EcalLogicID::NULLID,
                               "EE_readout_tower",12 );
  
  std::cout<<"got the end cap logic id set"<< endl;

  if (writeToDB_) {
    std::cout<<"Getting the latest ids for this tag (latest version) "<< endl;

    FEConfigMainInfo fe_main_info;
    fe_main_info.setConfigTag(tag_);
    try {
      std::cout << "trying to read previous tag if it exists tag="<< tag_<< ".version"<<version_<< endl;

      db_-> fetchConfigSet(&fe_main_info);
      if(fe_main_info.getPedId()>0 && ped_conf_id_ ==0 ) ped_conf_id_=fe_main_info.getPedId();
      lin_conf_id_=fe_main_info.getLinId();
      lut_conf_id_=fe_main_info.getLUTId();
      wei_conf_id_=fe_main_info.getWeiId();
      fgr_conf_id_=fe_main_info.getFgrId();
      sli_conf_id_=fe_main_info.getSliId();
      spi_conf_id_=fe_main_info.getSpiId(); //modif-alex 21/01/11
      del_conf_id_=fe_main_info.getTimId(); //modif-alex 21/01/11
      if(fe_main_info.getBxtId()>0 && bxt_conf_id_==0 ) bxt_conf_id_=fe_main_info.getBxtId();
      if(fe_main_info.getBttId()>0 && btt_conf_id_==0 ) btt_conf_id_=fe_main_info.getBttId();
      if(fe_main_info.getBstId()>0 && bst_conf_id_==0 ) bst_conf_id_=fe_main_info.getBstId();
      // those that are not written specifically in this program are propagated
      // from the previous record with the same tag and the highest version

      std::cout<<"got it "<< endl;

    } catch (exception &e) {
      cout << " tag did not exist a new tag will be created " << endl;
    } catch (...) {
      cout << "Unknown error caught" << endl;
    }

  }

  // Compute linearization coeff section //

  map<EcalLogicID, FEConfigPedDat> pedset ;
  map<EcalLogicID, FEConfigLinDat> linset ;
  map<EcalLogicID, FEConfigLinParamDat> linparamset ;
  map<EcalLogicID, FEConfigLUTParamDat> lutparamset ;
  map<EcalLogicID, FEConfigFgrParamDat> fgrparamset ;

  map<int, linStruc> linEtaSlice ;
  map< vector<int>, linStruc > linMap ;

  // count number of strip per tower
  int NbOfStripPerTCC[108][68] ; 
  for (int i=0 ; i<108 ; i++)
    for (int j=0 ; j<68 ; j++) 
      NbOfStripPerTCC[i][j] = 0 ;
  const std::vector<DetId> & ebCells = theBarrelGeometry_->getValidDetIds(DetId::Ecal, EcalBarrel);
  const std::vector<DetId> & eeCells = theEndcapGeometry_->getValidDetIds(DetId::Ecal, EcalEndcap);
  for (vector<DetId>::const_iterator it = ebCells.begin(); it != ebCells.end(); ++it) {
    EBDetId id(*it) ;
    const EcalTrigTowerDetId towid= id.tower();
    const EcalTriggerElectronicsId elId = theMapping_->getTriggerElectronicsId(id) ;
    int tccNb = theMapping_->TCCid(towid) ;
    int towerInTCC = theMapping_->iTT(towid) ;
    int stripInTower = elId.pseudoStripId() ;  
    if (stripInTower>NbOfStripPerTCC[tccNb-1][towerInTCC-1]) NbOfStripPerTCC[tccNb-1][towerInTCC-1] = stripInTower ;
  }
  for (vector<DetId>::const_iterator it = eeCells.begin(); it != eeCells.end(); ++it) {
    EEDetId id(*it) ;
    const EcalTrigTowerDetId towid= (*eTTmap_).towerOf(id) ;
    const EcalTriggerElectronicsId elId = theMapping_->getTriggerElectronicsId(id) ;
    int tccNb = theMapping_->TCCid(towid) ;
    int towerInTCC = theMapping_->iTT(towid) ; 
    int stripInTower = elId.pseudoStripId() ;
    if (stripInTower>NbOfStripPerTCC[tccNb-1][towerInTCC-1]) NbOfStripPerTCC[tccNb-1][towerInTCC-1] = stripInTower ;
  }



  // loop on EB xtals
  if (writeToFiles_) (*out_file_)<<"COMMENT ====== barrel crystals ====== "<<std::endl ;

  // special case of eta slices
  for (vector<DetId>::const_iterator it = ebCells.begin(); it != ebCells.end(); ++it) {
    EBDetId id(*it) ;
    double theta = theBarrelGeometry_->getGeometry(id)->getPosition().theta() ;
    if (!useTransverseEnergy_) theta = acos(0.) ;
    const EcalTrigTowerDetId towid= id.tower();
    const EcalTriggerElectronicsId elId = theMapping_->getTriggerElectronicsId(id) ;
    int dccNb = theMapping_->DCCid(towid) ;
    int tccNb = theMapping_->TCCid(towid) ;
    int towerInTCC = theMapping_->iTT(towid) ; // from 1 to 68 (EB)
    int stripInTower = elId.pseudoStripId() ;  // from 1 to 5
    int xtalInStrip = elId.channelId() ;       // from 1 to 5
    const EcalElectronicsId Id = theMapping_->getElectronicsId(id) ;
    int CCUid = Id.towerId() ;
    int VFEid = Id.stripId() ;
    int xtalInVFE = Id.xtalId() ;
    int xtalWithinCCUid = 5*(VFEid-1) + xtalInVFE -1 ; // Evgueni expects [0,24]
         
    (*geomFile_)<<"dccNb = "<<dccNb<<" tccNb = "<<tccNb<<" towerInTCC = "<<towerInTCC
        <<" stripInTower = "<<stripInTower<<" xtalInStrip = "<<xtalInStrip
        <<" CCUid = "<<CCUid<<" VFEid = "<<VFEid<<" xtalInVFE = "<<xtalInVFE
        <<" xtalWithinCCUid = "<<xtalWithinCCUid<<" ieta = "<<id.ieta()<<" iphi = "<<id.iphi()
        <<" xtalhashedId = "<<id.hashedIndex()<<" xtalNb = "<<id.ic()
        <<" ietaTT = "<<towid.ieta()<<" iphiTT = "<<towid.iphi()<<endl ;

    int TCCch = towerInTCC ;
    int SLBslot = int((towerInTCC-1)/8.) + 1 ;
    int SLBch = (towerInTCC-1)%8 + 1 ;
    int cmsswId = id.rawId() ;
    int ixtal=(id.ism()-1)*1700+(id.ic()-1);
    EcalLogicID logicId =my_EcalLogicId[ixtal];
    int dbId = logicId.getLogicID() ;
    int val[] = {dccNb+600,tccNb,towerInTCC,stripInTower,
         xtalInStrip,CCUid,VFEid,xtalInVFE,xtalWithinCCUid,id.ieta(),id.iphi(),
         -999,-999,towid.ieta()/abs(towid.ieta()),id.hashedIndex(),
         id.ic(), cmsswId, dbId, 
         towid.ieta(),towid.iphi(), TCCch, getCrate(tccNb).second, SLBch, SLBslot, 
         getGCTRegionEta(towid.ieta()),getGCTRegionPhi(towid.iphi())} ;
    for (int i=0 ; i<26 ; i++) ntupleInts_[i] = val[i] ;
    
    strcpy(ntupleDet_,getDet(tccNb).c_str()) ;
    strcpy(ntupleCrate_,getCrate(tccNb).first.c_str()) ;
    ntuple->Fill() ;
    
    
    if (tccNb == 37 && stripInTower == 3 && xtalInStrip == 3 && (towerInTCC-1)%4==0) {
      int etaSlice = towid.ietaAbs() ;
      coeffStruc coeff ;
      getCoeff(coeff, calibMap, id.rawId()) ;
      getCoeff(coeff, gainMap, id.rawId()) ;
      getCoeff(coeff, pedMap, id.rawId()) ;
      linStruc lin ;
      for (int i=0 ; i<3 ; i++) {
    int mult, shift ;
    bool ok = computeLinearizerParam(theta, coeff.gainRatio_[i], coeff.calibCoeff_, "EB", mult , shift) ;
    if (!ok) std::cout << "unable to compute the parameters for SM="<< id.ism()<<" xt="<< id.ic()<<" " <<dec<<id.rawId()<<std::endl ;  
    else {
      lin.pedestal_[i] = coeff.pedestals_[i] ;
      lin.mult_[i] = mult ;
      lin.shift_[i] = shift ;
    }
      }

      bool ok(true) ;
      if (forcedPedestalValue_ == -2) ok = realignBaseline(lin, 0) ;
      if (!ok) std::cout<<"SM="<< id.ism()<<" xt="<< id.ic()<<" " <<dec<<id.rawId()<<std::endl ; 
      linEtaSlice[etaSlice] = lin ; 
    }
  }    

  // general case
  for (vector<DetId>::const_iterator it = ebCells.begin(); it != ebCells.end(); ++it) {
    EBDetId id(*it) ;
    double theta = theBarrelGeometry_->getGeometry(id)->getPosition().theta() ;
    if (!useTransverseEnergy_) theta = acos(0.) ;
    const EcalTrigTowerDetId towid= id.tower();
    towerListEB.push_back(towid.rawId()) ;
    const EcalTriggerElectronicsId elId = theMapping_->getTriggerElectronicsId(id) ;
    stripListEB.push_back(elId.rawId() & 0xfffffff8) ;
    int dccNb = theMapping_->DCCid(towid) ;
    //int tccNb = theMapping_->TCCid(towid) ;
    int towerInTCC = theMapping_->iTT(towid) ; // from 1 to 68 (EB)
    //int stripInTower = elId.pseudoStripId() ;  // from 1 to 5
    //int xtalInStrip = elId.channelId() ;       // from 1 to 5
    const EcalElectronicsId Id = theMapping_->getElectronicsId(id) ;
    int CCUid = Id.towerId() ;
    int VFEid = Id.stripId() ;
    int xtalInVFE = Id.xtalId() ;
    int xtalWithinCCUid = 5*(VFEid-1) + xtalInVFE -1 ; // Evgueni expects [0,24]
    int etaSlice = towid.ietaAbs() ;

    // hashed index of strip EcalLogicID:
    int hashedStripLogicID = 68*5*(id.ism()-1) + 5*(towerInTCC-1) +  (VFEid-1) ;
    stripMapEB[hashedStripLogicID] = elId.rawId() & 0xfffffff8 ;
    stripMapEBsintheta[elId.rawId() & 0xfffffff8] = SFGVB_Threshold_ + abs(int(sin(theta)*pedestal_offset_)) ;

    //modif-debug
    /*
    FEConfigFgrEEStripDat stripdebug;
    EcalLogicID thestrip_debug = my_StripEcalLogicId[hashedStripLogicID] ;
    if(towid.ieta() == 6 && towid.iphi() == 7){
      std::cout << "xtal info=" << id << " VFE=" << VFEid << std::endl;
      std::cout << "TOWER DEBUG ieta=" << towid.ieta() << " iphi=" << towid.iphi() << " SFGVB=" << SFGVB_Threshold_ + abs(int(sin(theta)*pedestal_offset_)) 
        << " dbId=" << (elId.rawId() & 0xfffffff8) << " " << hashedStripLogicID << " " << thestrip_debug.getLogicID() << std::endl;  //modif-debug
    }//EB+3 TT24
    */
    //std::cout<<std::dec<<SFGVB_Threshold_ + abs(int(sin(theta)*pedestal_offset_))<<" "<<SFGVB_Threshold_<<" "<<abs(int(sin(theta)*pedestal_offset_))<<std::endl ;

    FEConfigPedDat pedDB ;
    FEConfigLinDat linDB ;
    if (writeToFiles_) (*out_file_)<<"CRYSTAL "<<dec<<id.rawId()<<std::endl ;
    //  if (writeToDB_) logicId = db_->getEcalLogicID ("EB_crystal_number", id.ism(), id.ic()) ;

    coeffStruc coeff ;
    getCoeff(coeff, calibMap, id.rawId()) ;
    if (H2_) coeff.calibCoeff_ = calibvec[id.ic()-1] ;
    getCoeff(coeff, gainMap, id.rawId()) ;
    getCoeff(coeff, pedMap, id.rawId()) ;
    ICEB->Fill(id.iphi(), id.ieta(), coeff.calibCoeff_) ;  
    IC->Fill(theBarrelGeometry_->getGeometry(id)->getPosition().phi(), theBarrelGeometry_->getGeometry(id)->getPosition().eta(), coeff.calibCoeff_) ;  

    vector<int> xtalCCU ;
    xtalCCU.push_back(dccNb+600) ; 
    xtalCCU.push_back(CCUid) ; 
    xtalCCU.push_back(xtalWithinCCUid) ;
    xtalCCU.push_back(id.rawId()) ;
    
    // compute and fill linearization parameters
    // case of eta slice
    if (forceEtaSlice_) {
      map<int, linStruc>::const_iterator itLin = linEtaSlice.find(etaSlice);
      if (itLin != linEtaSlice.end()) {
    linMap[xtalCCU] = itLin->second ;
    if (writeToFiles_) {
      for (int i=0 ; i<3 ; i++) 
        (*out_file_) << hex <<" 0x"<<itLin->second.pedestal_[i]<<" 0x"<<itLin->second.mult_[i]<<" 0x"<<itLin->second.shift_[i]<<std::endl;
    }
    if (writeToDB_) {
      for (int i=0 ; i<3 ; i++) {
        if (i==0)  {pedDB.setPedMeanG12(itLin->second.pedestal_[i]) ; linDB.setMultX12(itLin->second.mult_[i]) ; linDB.setShift12(itLin->second.shift_[i]) ; } 
        if (i==1)  {pedDB.setPedMeanG6(itLin->second.pedestal_[i]) ; linDB.setMultX6(itLin->second.mult_[i]) ; linDB.setShift6(itLin->second.shift_[i]) ; } 
        if (i==2)  {pedDB.setPedMeanG1(itLin->second.pedestal_[i]) ; linDB.setMultX1(itLin->second.mult_[i]) ; linDB.setShift1(itLin->second.shift_[i]) ; } 
      }
    }
    float factor = float(itLin->second.mult_[0])*pow(2.,-itLin->second.shift_[0])/xtal_LSB_EB_ ;
    tpgFactorEB->Fill(id.iphi(), id.ieta(), factor) ;
    tpgFactor->Fill(theBarrelGeometry_->getGeometry(id)->getPosition().phi(), 
            theBarrelGeometry_->getGeometry(id)->getPosition().eta(), factor) ;
      }
      else std::cout<<"current EtaSlice = "<<etaSlice<<" not found in the EtaSlice map"<<std::endl ;
    }
    else {
      // general case
      linStruc lin ;
      int forceBase12 = 0 ;
      for (int i=0 ; i<3 ; i++) {
    int mult, shift ;
    bool ok = computeLinearizerParam(theta, coeff.gainRatio_[i], coeff.calibCoeff_, "EB", mult , shift) ;
    if (!ok) std::cout << "unable to compute the parameters for SM="<< id.ism()<<" xt="<< id.ic()<<" " <<dec<<id.rawId()<<std::endl ;  
    else {
      //PP begin
      //  mult = 0 ; shift = 0 ;
      //  if (CCUid==1 && xtalWithinCCUid==21) {
      //    if (i==0) {mult = 0x80 ; shift = 0x3 ;}
      //    if (i==1) {mult = 0x80 ; shift = 0x2 ;}
      //    if (i==2) {mult = 0xc0 ; shift = 0x0 ;}
      //  } 
      //PP end
      double base = coeff.pedestals_[i] ;
      if (forcedPedestalValue_ == -3 && i==0) {
        double G = mult*pow(2.0,-(shift+2)) ;
        double g = G/sin(theta) ;
        // int pedestal = coeff.pedestals_[i] ;
        base = double(coeff.pedestals_[i]) - pedestal_offset_/g ;
        if (base<0.) base = 0 ;
        forceBase12 = int(base) ;
      }
      lin.pedestal_[i] = coeff.pedestals_[i] ;
      lin.mult_[i] = mult ;
      lin.shift_[i] = shift ;

//    if (xtalWithinCCUid != 14) {
//      forceBase12 = 0 ;
//      lin.pedestal_[i] = 0 ; 
//      lin.mult_[i] = 0 ;
//      lin.shift_[i] = 0 ;
//    }

    }
      }

      bool ok(true) ;
      if (forcedPedestalValue_ == -2) ok = realignBaseline(lin, 0) ;
      if (forcedPedestalValue_ == -3) ok = realignBaseline(lin, forceBase12) ;
      if (!ok) std::cout<<"SM="<< id.ism()<<" xt="<< id.ic()<<" " <<dec<<id.rawId()<<std::endl ; 
      
      for (int i=0 ; i<3 ; i++) {      
    if (writeToFiles_) (*out_file_) << hex <<" 0x"<<lin.pedestal_[i]<<" 0x"<<lin.mult_[i]<<" 0x"<<lin.shift_[i]<<std::endl; 
    if (writeToDB_) {
      if (i==0)  {pedDB.setPedMeanG12(lin.pedestal_[i]) ; linDB.setMultX12(lin.mult_[i]) ; linDB.setShift12(lin.shift_[i]) ; } 
      if (i==1)  {pedDB.setPedMeanG6(lin.pedestal_[i]) ; linDB.setMultX6(lin.mult_[i]) ; linDB.setShift6(lin.shift_[i]) ; } 
      if (i==2)  {pedDB.setPedMeanG1(lin.pedestal_[i]) ; linDB.setMultX1(lin.mult_[i]) ; linDB.setShift1(lin.shift_[i]) ; }
    }
    if (i==0) {
      float factor = float(lin.mult_[i])*pow(2.,-lin.shift_[i])/xtal_LSB_EB_ ;
      tpgFactorEB->Fill(id.iphi(), id.ieta(), factor) ;
      tpgFactor->Fill(theBarrelGeometry_->getGeometry(id)->getPosition().phi(), 
              theBarrelGeometry_->getGeometry(id)->getPosition().eta(), factor) ;               
    }
    double G = lin.mult_[i]*pow(2.0,-(lin.shift_[i]+2)) ;
    double g = G/sin(theta) ;
    float val[] = { float(i),
            float(theta),
            float(G),
            float(g),
            float(coeff.pedestals_[i]),
            float(lin.pedestal_[i])};// first arg = gainId (0 means gain12)
    ntupleSpike->Fill(val) ;
      }
      linMap[xtalCCU] = lin ;
    }
    if (writeToDB_) {
    // 1700 crystals/SM in the ECAL barrel
      int ixtal=(id.ism()-1)*1700+(id.ic()-1);
      EcalLogicID logicId =my_EcalLogicId[ixtal];
      pedset[logicId] = pedDB ;
      linset[logicId] = linDB ; 
    }
  } //ebCells

  if (writeToDB_) {
    // EcalLogicID  of the whole barrel is: my_EcalLogicId_EB
    FEConfigLinParamDat linparam ;
    linparam.setETSat(Et_sat_EB_); 
    linparamset[my_EcalLogicId_EB] = linparam ;

    FEConfigFgrParamDat fgrparam ;
    fgrparam.setFGlowthresh(FG_lowThreshold_EB_); 
    fgrparam.setFGhighthresh(FG_highThreshold_EB_); 
    fgrparam.setFGlowratio(FG_lowRatio_EB_); 
    fgrparam.setFGhighratio(FG_highRatio_EB_);
    fgrparamset[my_EcalLogicId_EB] = fgrparam ;


    FEConfigLUTParamDat lutparam ;
    lutparam.setETSat(Et_sat_EB_); 
    lutparam.setTTThreshlow(TTF_lowThreshold_EB_); 
    lutparam.setTTThreshhigh(TTF_highThreshold_EB_); 
    lutparamset[my_EcalLogicId_EB] = lutparam ;

  }


  // loop on EE xtals
  if (writeToFiles_) (*out_file_)<<"COMMENT ====== endcap crystals ====== "<<std::endl ;
  
  // special case of eta slices
  for (vector<DetId>::const_iterator it = eeCells.begin(); it != eeCells.end(); ++it) {
    EEDetId id(*it) ;
    double theta = theEndcapGeometry_->getGeometry(id)->getPosition().theta() ;
    if (!useTransverseEnergy_) theta = acos(0.) ;
    const EcalTrigTowerDetId towid= (*eTTmap_).towerOf(id) ;
    const EcalTriggerElectronicsId elId = theMapping_->getTriggerElectronicsId(id) ;
    const EcalElectronicsId Id = theMapping_->getElectronicsId(id) ;
    int dccNb = Id.dccId() ;
    int tccNb = theMapping_->TCCid(towid) ;
    int towerInTCC = theMapping_->iTT(towid) ; 
    int stripInTower = elId.pseudoStripId() ;
    int xtalInStrip = elId.channelId() ;
    int CCUid = Id.towerId() ;
    int VFEid = Id.stripId() ;
    int xtalInVFE = Id.xtalId() ;
    int xtalWithinCCUid = 5*(VFEid-1) + xtalInVFE -1 ; // Evgueni expects [0,24]

    (*geomFile_)<<"dccNb = "<<dccNb<<" tccNb = "<<tccNb<<" towerInTCC = "<<towerInTCC
        <<" stripInTower = "<<stripInTower<<" xtalInStrip = "<<xtalInStrip
        <<" CCUid = "<<CCUid<<" VFEid = "<<VFEid<<" xtalInVFE = "<<xtalInVFE
        <<" xtalWithinCCUid = "<<xtalWithinCCUid<<" ix = "<<id.ix()<<" iy = "<<id.iy()
        <<" xtalhashedId = "<<id.hashedIndex()<<" xtalNb = "<<id.isc()
        <<" ietaTT = "<<towid.ieta()<<" iphiTT = "<<towid.iphi()<<endl ;

    int TCCch = stripInTower ;
    int SLBslot, SLBch ;
    if (towerInTCC<5) {
      SLBslot = 1 ;
      SLBch = 4 + towerInTCC ;
    }
    else {
      SLBslot = int((towerInTCC-5)/8.) + 2 ;
      SLBch = (towerInTCC-5)%8 + 1 ;
    } 
    for (int j=0 ; j<towerInTCC-1 ; j++) TCCch += NbOfStripPerTCC[tccNb-1][j] ;

    int cmsswId = id.rawId() ;
    EcalLogicID logicId ;
    int iz = id.positiveZ() ;
    if (iz ==0) iz = -1 ;
    for (int k=0; k<(int)my_CrystalEcalLogicId_EE.size(); k++) {
      int z= my_CrystalEcalLogicId_EE[k].getID1() ;
      int x= my_CrystalEcalLogicId_EE[k].getID2() ;
      int y= my_CrystalEcalLogicId_EE[k].getID3() ;
      if (id.ix()==x && id.iy()==y && iz==z) logicId=my_CrystalEcalLogicId_EE[k];
    }
    int dbId = logicId.getLogicID() ;
    
    int val[] = {dccNb+600,tccNb,towerInTCC,stripInTower,
         xtalInStrip,CCUid,VFEid,xtalInVFE,xtalWithinCCUid,-999,-999,
         id.ix(),id.iy(),towid.ieta()/abs(towid.ieta()),id.hashedIndex(),
         id.ic(),cmsswId, dbId, 
         towid.ieta(),towid.iphi(),TCCch, getCrate(tccNb).second, SLBch, SLBslot, 
         getGCTRegionEta(towid.ieta()),getGCTRegionPhi(towid.iphi())} ;
    for (int i=0 ; i<26 ; i++) ntupleInts_[i] = val[i] ;
    strcpy(ntupleDet_,getDet(tccNb).c_str()) ;
    strcpy(ntupleCrate_,getCrate(tccNb).first.c_str()) ;
    ntuple->Fill() ;
     
    if ((tccNb == 76 || tccNb == 94) && stripInTower == 1 && xtalInStrip == 3 && (towerInTCC-1)%4==0) {
      int etaSlice = towid.ietaAbs() ;
      coeffStruc coeff ;
      getCoeff(coeff, calibMap, id.rawId()) ;
      getCoeff(coeff, gainMap, id.rawId()) ;
      getCoeff(coeff, pedMap, id.rawId()) ;
      linStruc lin ;
      for (int i=0 ; i<3 ; i++) {
    int mult, shift ;
    bool ok = computeLinearizerParam(theta, coeff.gainRatio_[i], coeff.calibCoeff_, "EE", mult , shift) ;
    if (!ok) std::cout << "unable to compute the parameters for Quadrant="<< id.iquadrant()<<" xt="<< id.ic()<<" " <<dec<<id.rawId()<<std::endl ;  
    else {
      lin.pedestal_[i] = coeff.pedestals_[i] ;
      lin.mult_[i] = mult ;
      lin.shift_[i] = shift ;
    }
      }

      bool ok(true) ;
      if (forcedPedestalValue_ == -2 || forcedPedestalValue_ == -3) ok = realignBaseline(lin, 0) ;
      if (!ok) std::cout<<"Quadrant="<< id.iquadrant()<<" xt="<< id.ic()<<" " <<dec<<id.rawId()<<std::endl ; 

      linEtaSlice[etaSlice] = lin ;
    }
  }

  // general case
  for (vector<DetId>::const_iterator it = eeCells.begin(); it != eeCells.end(); ++it) {
    EEDetId id(*it);
    double theta = theEndcapGeometry_->getGeometry(id)->getPosition().theta() ;
    if (!useTransverseEnergy_) theta = acos(0.) ;
    const EcalTrigTowerDetId towid= (*eTTmap_).towerOf(id) ;
    const EcalTriggerElectronicsId elId = theMapping_->getTriggerElectronicsId(id) ;
    const EcalElectronicsId Id = theMapping_->getElectronicsId(id) ;
    towerListEE.push_back(towid.rawId()) ;
    // special case of towers in inner rings of EE
    if (towid.ietaAbs() == 27 || towid.ietaAbs() == 28) {
      EcalTrigTowerDetId additionalTower(towid.zside(), towid.subDet(), towid.ietaAbs(), towid.iphi()+1) ;
      towerListEE.push_back(additionalTower.rawId()) ;
    }
    stripListEE.push_back(elId.rawId() & 0xfffffff8) ;
    int dccNb = Id.dccId() ;
    //int tccNb = theMapping_->TCCid(towid) ;
    //int towerInTCC = theMapping_->iTT(towid) ;
    //int stripInTower = elId.pseudoStripId() ;
    //int xtalInStrip = elId.channelId() ;
    int CCUid = Id.towerId() ;
    int VFEid = Id.stripId() ;
    int xtalInVFE = Id.xtalId() ;
    int xtalWithinCCUid = 5*(VFEid-1) + xtalInVFE -1 ; // Evgueni expects [0,24]
    int etaSlice = towid.ietaAbs() ;

    EcalLogicID logicId ;
    FEConfigPedDat pedDB ;
    FEConfigLinDat linDB ;
    if (writeToFiles_) (*out_file_)<<"CRYSTAL "<<dec<<id.rawId()<<std::endl ;
    if (writeToDB_ && DBEE_) {
      int iz = id.positiveZ() ;
      if (iz ==0) iz = -1 ;      
      for (int k=0; k<(int)my_CrystalEcalLogicId_EE.size(); k++) {
    int z= my_CrystalEcalLogicId_EE[k].getID1() ;
    int x= my_CrystalEcalLogicId_EE[k].getID2() ;
    int y= my_CrystalEcalLogicId_EE[k].getID3() ;   
    if (id.ix()==x && id.iy()==y && iz==z) logicId=my_CrystalEcalLogicId_EE[k];
      }
    }
    
    coeffStruc coeff ;
    getCoeff(coeff, calibMap, id.rawId()) ;
    getCoeff(coeff, gainMap, id.rawId()) ;
    getCoeff(coeff, pedMap, id.rawId()) ;
    if (id.zside()>0) ICEEPlus->Fill(id.ix(), id.iy(), coeff.calibCoeff_) ;  
    else ICEEMinus->Fill(id.ix(), id.iy(), coeff.calibCoeff_) ;  
    IC->Fill(theEndcapGeometry_->getGeometry(id)->getPosition().phi(), theEndcapGeometry_->getGeometry(id)->getPosition().eta(), coeff.calibCoeff_) ;  
  
    vector<int> xtalCCU ;
    xtalCCU.push_back(dccNb+600) ; 
    xtalCCU.push_back(CCUid) ; 
    xtalCCU.push_back(xtalWithinCCUid) ;
    xtalCCU.push_back(id.rawId()) ;

    // compute and fill linearization parameters
    // case of eta slice
    if (forceEtaSlice_) {
      map<int, linStruc>::const_iterator itLin = linEtaSlice.find(etaSlice);
      if (itLin != linEtaSlice.end()) {
    linMap[xtalCCU] = itLin->second ;
    if (writeToFiles_) {
      for (int i=0 ; i<3 ; i++) 
        (*out_file_) << hex <<" 0x"<<itLin->second.pedestal_[i]<<" 0x"<<itLin->second.mult_[i]<<" 0x"<<itLin->second.shift_[i]<<std::endl;
    }
    if (writeToDB_ && DBEE_) {
      for (int i=0 ; i<3 ; i++) {
        if (i==0)  {pedDB.setPedMeanG12(itLin->second.pedestal_[i]) ; linDB.setMultX12(itLin->second.mult_[i]) ; linDB.setShift12(itLin->second.shift_[i]) ; } 
        if (i==1)  {pedDB.setPedMeanG6(itLin->second.pedestal_[i]) ; linDB.setMultX6(itLin->second.mult_[i]) ; linDB.setShift6(itLin->second.shift_[i]) ; } 
        if (i==2)  {pedDB.setPedMeanG1(itLin->second.pedestal_[i]) ; linDB.setMultX1(itLin->second.mult_[i]) ; linDB.setShift1(itLin->second.shift_[i]) ; } 
      }
    }
    float factor = float(itLin->second.mult_[0])*pow(2.,-itLin->second.shift_[0])/xtal_LSB_EE_ ;
    if (id.zside()>0) tpgFactorEEPlus->Fill(id.ix(), id.iy(), factor) ;
    else tpgFactorEEMinus->Fill(id.ix(), id.iy(), factor) ;
    tpgFactor->Fill(theEndcapGeometry_->getGeometry(id)->getPosition().phi(), 
            theEndcapGeometry_->getGeometry(id)->getPosition().eta(), factor) ;
      }
      else std::cout<<"current EtaSlice = "<<etaSlice<<" not found in the EtaSlice map"<<std::endl ;      
    }
    else {
      // general case
      linStruc lin ;
      for (int i=0 ; i<3 ; i++) {
    int mult, shift ;
    bool ok = computeLinearizerParam(theta, coeff.gainRatio_[i], coeff.calibCoeff_, "EE", mult , shift) ;
    if (!ok) std::cout << "unable to compute the parameters for "<<dec<<id.rawId()<<std::endl ;  
    else {
      lin.pedestal_[i] = coeff.pedestals_[i] ;
      lin.mult_[i] = mult ;
      lin.shift_[i] = shift ;
    }
      }

      bool ok(true) ;
      if (forcedPedestalValue_ == -2 || forcedPedestalValue_ == -3) ok = realignBaseline(lin, 0) ;
      if (!ok) std::cout<<"Quadrant="<< id.iquadrant()<<" xt="<< id.ic()<<" " <<dec<<id.rawId()<<std::endl ; 

      for (int i=0 ; i<3 ; i++) {
    if (writeToFiles_) (*out_file_) << hex <<" 0x"<<lin.pedestal_[i]<<" 0x"<<lin.mult_[i]<<" 0x"<<lin.shift_[i]<<std::endl; 
    if (writeToDB_ && DBEE_) {
      if (i==0)  {pedDB.setPedMeanG12(lin.pedestal_[i]) ; linDB.setMultX12(lin.mult_[i]) ; linDB.setShift12(lin.shift_[i]) ; } 
      if (i==1)  {pedDB.setPedMeanG6(lin.pedestal_[i]) ; linDB.setMultX6(lin.mult_[i]) ; linDB.setShift6(lin.shift_[i]) ; } 
      if (i==2)  {pedDB.setPedMeanG1(lin.pedestal_[i]) ; linDB.setMultX1(lin.mult_[i]) ; linDB.setShift1(lin.shift_[i]) ; }
    }
    if (i==0) {
      float factor = float(lin.mult_[i])*pow(2.,-lin.shift_[i])/xtal_LSB_EE_ ;
      if (id.zside()>0) tpgFactorEEPlus->Fill(id.ix(), id.iy(), factor) ;
      else tpgFactorEEMinus->Fill(id.ix(), id.iy(), factor) ;       
      tpgFactor->Fill(theEndcapGeometry_->getGeometry(id)->getPosition().phi(), 
              theEndcapGeometry_->getGeometry(id)->getPosition().eta(), factor) ;                   
    }
      }
      linMap[xtalCCU] = lin ;
    }
    if (writeToDB_ && DBEE_) {
      pedset[logicId] = pedDB ;
      linset[logicId] = linDB ;
    }
  } //eeCells

  if (writeToDB_ ) {
    // EcalLogicID  of the whole barrel is: my_EcalLogicId_EB
    FEConfigLinParamDat linparam ;
    linparam.setETSat(Et_sat_EE_);
    linparamset[my_EcalLogicId_EE] = linparam ;

    FEConfigLUTParamDat lutparam ;
    lutparam.setETSat(Et_sat_EE_);
    lutparam.setTTThreshlow(TTF_lowThreshold_EE_); 
    lutparam.setTTThreshhigh(TTF_highThreshold_EE_);
    lutparamset[my_EcalLogicId_EE] = lutparam ;

 
    FEConfigFgrParamDat fgrparam ;
    fgrparam.setFGlowthresh(FG_Threshold_EE_); 
    fgrparam.setFGhighthresh(FG_Threshold_EE_); 
    fgrparamset[my_EcalLogicId_EE] = fgrparam ;

  }

  if (writeToDB_) {

    ostringstream ltag;
    ltag.str("EB_"); ltag<<Et_sat_EB_<<"_EE_"<<Et_sat_EE_;
    std::string lin_tag=ltag.str();
    std::cout<< " LIN tag "<<lin_tag<<endl;

    if(m_write_ped==1) {
      ped_conf_id_=db_->writeToConfDB_TPGPedestals(pedset, 1, "from_OfflineDB") ;
    } else {
      std::cout<< "the ped id ="<<ped_conf_id_<<" will be used for the pedestals "<<std::endl; 
    }

    if(m_write_lin==1) lin_conf_id_=db_->writeToConfDB_TPGLinearCoef(linset,linparamset, 1, lin_tag) ;
  }

  // Evgueni interface
  std::ofstream evgueni("TPG_hardcoded.hh", std::ios::out) ;  
  evgueni<<"void getLinParamTPG_hardcoded(int fed, int ccu, int xtal,"<<endl ;
  evgueni<<"                        int & mult12, int & shift12, int & base12,"<<endl ;
  evgueni<<"                        int & mult6, int & shift6, int & base6,"<<endl ;
  evgueni<<"                        int & mult1, int & shift1, int & base1)"<<endl ;
  evgueni<<"{"<<endl;
  evgueni<<"  mult12 = 0 ; shift12 = 0 ; base12 = 0 ; mult6 = 0 ; shift6 = 0 ; base6 = 0 ; mult1 = 0 ; shift1 = 0 ; base1 = 0 ;"<<endl ;
  map< vector<int>, linStruc>::const_iterator itLinMap ;
  for (itLinMap = linMap.begin() ; itLinMap != linMap.end() ; itLinMap++) {
    vector<int> xtalInCCU = itLinMap->first ;
    evgueni<<"  if (fed=="<<xtalInCCU[0]<<" && ccu=="<<xtalInCCU[1]<<" && xtal=="<<xtalInCCU[2]<<") {" ;
    evgueni<<"  mult12 = "<<itLinMap->second.mult_[0]<<" ; shift12 = "<<itLinMap->second.shift_[0]<<" ; base12 = "<<itLinMap->second.pedestal_[0]<<" ; " ;
    evgueni<<"  mult6 = "<<itLinMap->second.mult_[1]<<" ; shift6 = "<<itLinMap->second.shift_[1]<<" ; base6 = "<<itLinMap->second.pedestal_[1]<<" ; " ;
    evgueni<<"  mult1 = "<<itLinMap->second.mult_[2]<<" ; shift1 = "<<itLinMap->second.shift_[2]<<" ; base1 = "<<itLinMap->second.pedestal_[2]<<" ; " ;
    evgueni<<"  return ;}" <<endl ;
  }
  evgueni<<"}" <<endl ;
  evgueni.close() ;



  // Compute weights section //

  const int NWEIGROUPS = 2 ; 
  std::vector<unsigned int> weights[NWEIGROUPS] ;

  EBShape shapeEB ;
  EEShape shapeEE ;
  weights[0] = computeWeights(shapeEB, hshapeEB) ;
  weights[1] = computeWeights(shapeEE, hshapeEE) ;

  map<EcalLogicID, FEConfigWeightGroupDat> dataset;

  for (int igrp=0 ; igrp<NWEIGROUPS ; igrp++) {

    if (weights[igrp].size() == 5) {

      if (writeToFiles_) {
    (*out_file_) <<std::endl ;
    (*out_file_) <<"WEIGHT "<<igrp<<endl ;
    for (unsigned int sample=0 ; sample<5 ; sample++) (*out_file_) << "0x" <<hex<<weights[igrp][sample]<<" " ;
    (*out_file_)<<std::endl ; 
    (*out_file_) <<std::endl ;
      }
      if (writeToDB_) {
    std::cout<<"going to write the weights for groupe:"<<igrp<<endl;
    FEConfigWeightGroupDat gut;
    gut.setWeightGroupId(igrp);
    //PP WARNING: weights order is reverted when stored in the DB 
    gut.setWeight0(weights[igrp][4] );
    gut.setWeight1(weights[igrp][3]+ 0x80); //0x80 to identify the max of the pulse in the FENIX (doesn't exist in emulator)
    gut.setWeight2(weights[igrp][2] );
    gut.setWeight3(weights[igrp][1] );
    gut.setWeight4(weights[igrp][0] );
    EcalLogicID ecid = EcalLogicID( "DUMMY", igrp,igrp); //1 dummy ID per group
    // Fill the dataset
    dataset[ecid] = gut;
      }
    }
  }

  if (writeToDB_) {

    // now we store in the DB the correspondence btw channels and groups 
    map<EcalLogicID, FEConfigWeightDat> dataset2;

    // EB loop
    for (int ich=0; ich<(int)my_StripEcalLogicId.size() ; ich++){
      FEConfigWeightDat wut;
      int igroup = 0; // this group is for EB
      wut.setWeightGroupId(igroup);
      dataset2[my_StripEcalLogicId[ich]] = wut;
    }
    
    // EE loop
    for (int ich=0; ich<(int)my_StripEcalLogicId1_EE.size() ; ich++){
      FEConfigWeightDat wut;
      int igroup = 1; // this group is for EE
      wut.setWeightGroupId(igroup);
      // Fill the dataset
      dataset2[my_StripEcalLogicId1_EE[ich]] = wut;
    }
    // EE loop 2 (we had to split the ids of EE in 2 vectors to avoid crash!)
    for (int ich=0; ich<(int)my_StripEcalLogicId2_EE.size() ; ich++){
      FEConfigWeightDat wut;
      int igroup = 1; // this group is for EE
      wut.setWeightGroupId(igroup);
      // Fill the dataset
      dataset2[my_StripEcalLogicId2_EE[ich]] = wut;
    }

    // Insert the datasets
    ostringstream wtag;
    wtag.str(""); wtag<<"Shape_NGroups_"<<NWEIGROUPS;
    std::string weight_tag=wtag.str();
    std::cout<< " weight tag "<<weight_tag<<endl; 
    if (m_write_wei==1) wei_conf_id_=db_->writeToConfDB_TPGWeight(dataset, dataset2, NWEIGROUPS, weight_tag) ;
  }      


  // Compute FG section //

  // barrel
  unsigned int lowRatio, highRatio, lowThreshold, highThreshold, lutFG ;
  computeFineGrainEBParameters(lowRatio, highRatio, lowThreshold, highThreshold, lutFG) ;
  if (writeToFiles_) {
    (*out_file_) <<std::endl ;
    (*out_file_) <<"FG 0"<<std::endl ;
    (*out_file_)<<hex<<"0x"<<lowThreshold<<" 0x"<<highThreshold
          <<" 0x"<<lowRatio<<" 0x"<<highRatio<<" 0x"<<lutFG
          <<std::endl ;
  }

  // endcap
  unsigned int threshold, lut_tower ;
  unsigned int lut_strip;
  computeFineGrainEEParameters(threshold, lut_strip, lut_tower) ; 

  // and here we store the fgr part

  
  if (writeToDB_) {
    std::cout<<"going to write the fgr "<< endl;
      map<EcalLogicID, FEConfigFgrGroupDat> dataset;
      // we create 1 group
      int NFGRGROUPS =1; 
      for (int ich=0; ich<NFGRGROUPS; ich++){
    FEConfigFgrGroupDat gut;
    gut.setFgrGroupId(ich);
    gut.setThreshLow(lowRatio );
    gut.setThreshHigh(highRatio);
    gut.setRatioLow(lowThreshold);
    gut.setRatioHigh(highThreshold);
    gut.setLUTValue(lutFG);
    EcalLogicID ecid = EcalLogicID( "DUMMY", ich,ich);
    // Fill the dataset
    dataset[ecid] = gut; // we use any logic id but different, because it is in any case ignored... 
      }
      
      // now we store in the DB the correspondence btw channels and groups 
      map<EcalLogicID, FEConfigFgrDat> dataset2;
      // in this case I decide in a stupid way which channel belongs to which group 
      for (int ich=0; ich<(int)my_TTEcalLogicId.size() ; ich++){
    FEConfigFgrDat wut;
    int igroup=0;
    wut.setFgrGroupId(igroup);
    // Fill the dataset
    // the logic ids are ordered by SM (1,...36) and TT (1,...68)  
    // you have to calculate the right index here 
    dataset2[my_TTEcalLogicId[ich]] = wut;
      }

      // endcap loop
      for (int ich=0; ich<(int)my_RTEcalLogicId_EE.size() ; ich++){
    //  std::cout << " endcap FGR " << std::endl;
    FEConfigFgrDat wut;
    int igroup=0;
    wut.setFgrGroupId(igroup);
    // Fill the dataset
    // the logic ids are ordered by .... ?  
    // you have to calculate the right index here 
    dataset2[my_RTEcalLogicId_EE[ich]] = wut;
      }

      // endcap TT loop for the FEfgr EE Tower
      map<EcalLogicID, FEConfigFgrEETowerDat> dataset3;
      for (int ich=0; ich<(int)my_TTEcalLogicId_EE.size() ; ich++){
    FEConfigFgrEETowerDat fgreett;
    fgreett.setLutValue(lut_tower);
    dataset3[my_TTEcalLogicId_EE[ich]]=fgreett;
      }

      // endcap strip loop for the FEfgr EE strip
      // and barrel strip loop for the spike parameters (same structure than EE FGr)
      map<EcalLogicID, FEConfigFgrEEStripDat> dataset4;
      for (int ich=0; ich<(int)my_StripEcalLogicId1_EE.size() ; ich++){
    FEConfigFgrEEStripDat zut;
    zut.setThreshold(threshold);
    zut.setLutFgr(lut_strip);
    dataset4[my_StripEcalLogicId1_EE[ich]] = zut;
      }
      for (int ich=0; ich<(int)my_StripEcalLogicId2_EE.size() ; ich++){
    FEConfigFgrEEStripDat zut;
    zut.setThreshold(threshold);
    zut.setLutFgr(lut_strip);
    // Fill the dataset
    dataset4[my_StripEcalLogicId2_EE[ich]] = zut;
      }
      for (int ich=0; ich<(int)my_StripEcalLogicId.size() ; ich++){
    // EB
    FEConfigFgrEEStripDat zut;
    EcalLogicID thestrip = my_StripEcalLogicId[ich] ;
    uint32_t elStripId = stripMapEB[ich] ;
    map<uint32_t, uint32_t>::const_iterator it = stripMapEBsintheta.find(elStripId) ;
    if (it != stripMapEBsintheta.end()) zut.setThreshold(it->second);
    else {
      cout<<"ERROR: strip SFGVB threshold parameter not found for that strip:"<<thestrip.getID1()<<" "<<thestrip.getID3()<<" "<<thestrip.getID3()<<endl ; 
      cout<<" using value = "<<SFGVB_Threshold_+pedestal_offset_<<endl ;
      zut.setThreshold(SFGVB_Threshold_+pedestal_offset_) ;
    }
    zut.setLutFgr(SFGVB_lut_);
    // Fill the dataset
    dataset4[thestrip] = zut;
      }

      // Insert the dataset
      ostringstream wtag;
      wtag.str(""); wtag<<"FGR_"<<lutFG<<"_N_"<<NFGRGROUPS<<"_eb_"<<FG_lowThreshold_EB_<<"_EB_"<<FG_highThreshold_EB_;
      std::string weight_tag=wtag.str();
      std::cout<< " weight tag "<<weight_tag<<endl; 
      if(m_write_fgr==1) fgr_conf_id_=db_->writeToConfDB_TPGFgr(dataset, dataset2,  fgrparamset,dataset3, dataset4, NFGRGROUPS, weight_tag) ;

      //modif-alex 21/01/11
      map<EcalLogicID, FEConfigSpikeDat> datasetspike; //loob EB TT
      for (int ich=0; ich<(int)my_TTEcalLogicId.size() ; ich++){
        FEConfigSpikeDat spiketh;
        spiketh.setSpikeThreshold(SFGVB_SpikeKillingThreshold_);
        datasetspike[my_TTEcalLogicId[ich]] = spiketh;
      }//loop EB TT towers    

      //modif-alex 21/01/11
      ostringstream stag;
      stag.str(""); stag<<"SpikeTh"<<SFGVB_SpikeKillingThreshold_;
      std::string spike_tag=stag.str();
      std::cout<< " spike tag "<<spike_tag<<endl;
      if(m_write_spi==1) spi_conf_id_=db_->writeToConfDB_Spike(datasetspike, spike_tag) ; //modif-alex 21/01/11

      //modif-alex 31/01/11
      //DELAYS EB
      map<EcalLogicID, FEConfigTimingDat> datasetdelay; // the loop goes from TCC 38 to 72 and throught the towers from 1 to 68 
      for (int ich=0; ich<(int)my_TTEcalLogicId_EB_by_TCC.size() ; ich++){
        FEConfigTimingDat delay;

    EcalLogicID logiciddelay = my_TTEcalLogicId_EB_by_TCC[ich] ;
    int id1_tcc=my_TTEcalLogicId_EB_by_TCC[ich].getID1(); // the TCC
    int id2_tt =my_TTEcalLogicId_EB_by_TCC[ich].getID2(); // the tower
    std::map<int, vector<int> >::const_iterator ittEB  = delays_EB_.find(id1_tcc); 
    std::vector<int> TimingDelaysEB = ittEB->second;

    if (ittEB != delays_EB_.end()){
      if(TimingDelaysEB[id2_tt-1] == -1){
          std::cout << "ERROR: Barrel timing delay not specified, check file, putting default value 1" << std::endl;
          delay.setTimingPar1(1); 
        }
      else delay.setTimingPar1(TimingDelaysEB[id2_tt-1]);
    } else {
      std::cout << "ERROR:Barrel Could not find delay parameter for that trigger tower " << std::endl;
      std::cout << "Using default value = 1" << std::endl;
      delay.setTimingPar1(1);
    }

    std::map<int, vector<int> >::const_iterator ittpEB  = phases_EB_.find(id1_tcc);
    std::vector<int> TimingPhasesEB = ittpEB->second;

        if (ittpEB != phases_EB_.end()){
          if(TimingPhasesEB[id2_tt-1] == -1){
        std::cout << "ERROR: Barrel timing phase not specified, check file, putting default value 0" << std::endl;
        delay.setTimingPar2(0);
      }
          else delay.setTimingPar2(TimingPhasesEB[id2_tt-1]);
        } else {
      std::cout << "ERROR:Barrel Could not find phase parameter for that trigger tower " << std::endl;
      std::cout << "Using default value = 0" << std::endl;
          delay.setTimingPar2(0);
        }

    std::cout << ich << " tcc=" << id1_tcc << " TT=" << id2_tt << " logicId=" << logiciddelay.getLogicID() 
          << " delay=" << TimingDelaysEB[id2_tt-1] << " phase=" << TimingPhasesEB[id2_tt-1] <<  std::endl; 
    
        //delay.setTimingPar1(1);
        //delay.setTimingPar2(2);
        datasetdelay[my_TTEcalLogicId_EB_by_TCC[ich]] = delay;
      }//loop EB TT towers

      //DELAYS EE
      int stripindex = 0;
      int tccin = 1;
      for (int ich=0; ich<(int)my_StripEcalLogicId_EE_strips_by_TCC.size() ; ich++){
        FEConfigTimingDat delay;
    //int id1_strip=my_StripEcalLogicId_EE_strips_by_TCC[ich].getID1(); // the TCC
    //int id2_strip=my_StripEcalLogicId_EE_strips_by_TCC[ich].getID2(); // the Tower
    //int id3_strip=my_StripEcalLogicId_EE_strips_by_TCC[ich].getID3(); // the strip

        EcalLogicID logiciddelay = my_StripEcalLogicId_EE_strips_by_TCC[ich] ;
        int id1_tcc=my_StripEcalLogicId_EE_strips_by_TCC[ich].getID1(); // the TCC
        int id2_tt =my_StripEcalLogicId_EE_strips_by_TCC[ich].getID2(); // the tower
        int id3_st =my_StripEcalLogicId_EE_strips_by_TCC[ich].getID3(); // the strip

    //reset strip counter
    if(id1_tcc != tccin) {tccin = id1_tcc; stripindex=0;}

    std::map<int, vector<int> >::const_iterator ittEE  = delays_EE_.find(id1_tcc);
    std::vector<int> TimingDelaysEE = ittEE->second;

        if (ittEE != delays_EE_.end()){
          if(TimingDelaysEE[stripindex] == -1){
        std::cout << "ERROR: Endcap timing delay not specified, check file, putting default value 1" << std::endl;
        delay.setTimingPar1(1);
      }
          else delay.setTimingPar1(TimingDelaysEE[stripindex]);
        } else {
      std::cout << "ERROR:Endcap Could not find delay parameter for that trigger tower " << std::endl;
      std::cout << "Using default value = 1" << std::endl;
          delay.setTimingPar1(1);
        }

    std::map<int, vector<int> >::const_iterator ittpEE  = phases_EE_.find(id1_tcc);
    std::vector<int> TimingPhasesEE = ittpEE->second;

    if (ittpEE != phases_EE_.end()){
          if(TimingPhasesEE[stripindex] == -1){
        std::cout << "ERROR: Endcap timing phase not specified, check file, putting default value 0" << std::endl;
            delay.setTimingPar2(0);
          }
          else delay.setTimingPar2(TimingPhasesEE[stripindex]);
        } else {
      std::cout << "ERROR:Endcap Could not find phase parameter for that trigger tower " << std::endl;
      std::cout << "Using default value = 0" << std::endl;
          delay.setTimingPar2(0);
        }

    std::cout << ich << " stripindex=" << stripindex << " tcc=" << id1_tcc << " TT=" << id2_tt << " id3_st=" << id3_st 
          << " logicId=" << logiciddelay.getLogicID()
                  << " delay=" << TimingDelaysEE[stripindex] << " phase=" << TimingPhasesEE[stripindex] <<  std::endl;


        //delay.setTimingPar1(1);
        //delay.setTimingPar2(2);
        datasetdelay[my_StripEcalLogicId_EE_strips_by_TCC[ich]] = delay;
    stripindex++;
      }//loop EE strip towers

      ostringstream de_tag;
      de_tag.str(""); de_tag<<"DelaysFromFile";
      std::string delay_tag=de_tag.str();
      std::cout<< " delay tag "<<delay_tag<<endl;
      if(m_write_del==1) del_conf_id_=db_->writeToConfDB_Delay(datasetdelay, delay_tag) ; //modif-alex 31/01/11   


  } //write to DB

  if (writeToDB_) {
    std::cout<<"going to write the sliding "<< endl;
      map<EcalLogicID, FEConfigSlidingDat> dataset;
      // in this case I decide in a stupid way which channel belongs to which group 
      for (int ich=0; ich<(int)my_StripEcalLogicId.size() ; ich++){
    FEConfigSlidingDat wut;
    wut.setSliding(sliding_);
    // Fill the dataset
    // the logic ids are ordered by SM (1,...36) , TT (1,...68) and strip (1..5) 
    // you have to calculate the right index here 
    dataset[my_StripEcalLogicId[ich]] = wut;
      }

      // endcap loop
      for (int ich=0; ich<(int)my_StripEcalLogicId1_EE.size() ; ich++){
    FEConfigSlidingDat wut;
    wut.setSliding(sliding_);
    // Fill the dataset
    // the logic ids are ordered by fed tower strip
    // you have to calculate the right index here 
    dataset[my_StripEcalLogicId1_EE[ich]] = wut;
      }
      for (int ich=0; ich<(int)my_StripEcalLogicId2_EE.size() ; ich++){
    FEConfigSlidingDat wut;
    wut.setSliding(sliding_);
    // Fill the dataset
    // the logic ids are ordered by ... ? 
    // you have to calculate the right index here 
    dataset[my_StripEcalLogicId2_EE[ich]] = wut;
      }
      

      // Insert the dataset
      ostringstream wtag;
      wtag.str(""); wtag<<"Sliding_"<<sliding_;
      std::string justatag=wtag.str();
      std::cout<< " sliding tag "<<justatag<<endl;
      int iov_id=0; // just a parameter ... 
      if(m_write_sli==1) sli_conf_id_=db_->writeToConfDB_TPGSliding(dataset,iov_id, justatag) ;
  }

  



  // Compute LUT section //

  int lut_EB[1024], lut_EE[1024] ;

  // barrel
  computeLUT(lut_EB, "EB") ; 
  if (writeToFiles_) {
    (*out_file_) <<std::endl ;
    (*out_file_) <<"LUT 0"<<std::endl ;
    for (int i=0 ; i<1024 ; i++) (*out_file_)<<"0x"<<hex<<lut_EB[i]<<endl ;
    (*out_file_)<<endl ;
  }
  
  // endcap
  computeLUT(lut_EE, "EE") ;
  // check first if lut_EB and lut_EE are the same
  bool newLUT(false) ;
  for (int i=0 ; i<1024 ; i++) if (lut_EE[i] != lut_EB[i]) newLUT = true ;
  if (newLUT && writeToFiles_) { 
    (*out_file_) <<std::endl ;
    (*out_file_) <<"LUT 1"<<std::endl ;
    for (int i=0 ; i<1024 ; i++) (*out_file_)<<"0x"<<hex<<lut_EE[i]<<endl ;
    (*out_file_)<<endl ;
  }



  if (writeToDB_) {



    map<EcalLogicID, FEConfigLUTGroupDat> dataset;
    // we create 1 LUT group
    int NLUTGROUPS =0; 
    int ich=0;
      FEConfigLUTGroupDat lut;
      lut.setLUTGroupId(ich);
      for (int i=0; i<1024; i++){
    lut.setLUTValue(i, lut_EB[i] );
      }
      EcalLogicID ecid = EcalLogicID( "DUMMY", ich,ich);
      // Fill the dataset
      dataset[ecid] = lut; // we use any logic id but different, because it is in any case ignored... 
      
      ich++;
      
      FEConfigLUTGroupDat lute;
      lute.setLUTGroupId(ich);
      for (int i=0; i<1024; i++){
    lute.setLUTValue(i, lut_EE[i] );
      }
      EcalLogicID ecide = EcalLogicID( "DUMMY", ich,ich);
      // Fill the dataset
      dataset[ecide] = lute; // we use any logic id but different, because it is in any case ignored...

      ich++;

      NLUTGROUPS=ich;

    // now we store in the DB the correspondence btw channels and LUT groups 
    map<EcalLogicID, FEConfigLUTDat> dataset2;
    // in this case I decide in a stupid way which channel belongs to which group 
    for (int ich=0; ich<(int)my_TTEcalLogicId.size() ; ich++){
      FEConfigLUTDat lut;
      int igroup=0;
      lut.setLUTGroupId(igroup);
      // calculate the right TT - in the vector they are ordered by SM and by TT 
      // Fill the dataset
      dataset2[my_TTEcalLogicId[ich]] = lut;
    }

    // endcap loop 
    for (int ich=0; ich<(int)my_TTEcalLogicId_EE.size() ; ich++){
      FEConfigLUTDat lut;
      int igroup=1;
      lut.setLUTGroupId(igroup);
      // calculate the right TT  
      // Fill the dataset
      dataset2[my_TTEcalLogicId_EE[ich]] = lut;
    }    

    // Insert the dataset
    ostringstream ltag;
    ltag.str(""); ltag<<LUT_option_<<"_NGroups_"<<NLUTGROUPS;
    std::string lut_tag=ltag.str();
    std::cout<< " LUT tag "<<lut_tag<<endl; 
    if(m_write_lut==1) lut_conf_id_=db_->writeToConfDB_TPGLUT(dataset, dataset2,lutparamset, NLUTGROUPS, lut_tag) ;

  }

  // last we insert the FE_CONFIG_MAIN table 
 if (writeToDB_) {
   
   //int conf_id_=db_->writeToConfDB_TPGMain(ped_conf_id_,lin_conf_id_, lut_conf_id_, fgr_conf_id_, 
   //               sli_conf_id_, wei_conf_id_, bxt_conf_id_, btt_conf_id_, tag_, version_) ;
   int conf_id_=db_->writeToConfDB_TPGMain(ped_conf_id_,lin_conf_id_, lut_conf_id_, fgr_conf_id_, 
                       sli_conf_id_, wei_conf_id_, spi_conf_id_, del_conf_id_, bxt_conf_id_, btt_conf_id_, bst_conf_id_, tag_, version_) ; //modif-alex 21/01/11
   
   std::cout << "\n Conf ID = " << conf_id_ << std::endl;

 }


  // loop on strips and associate them with  values //

  // Barrel
  stripListEB.sort() ;
  stripListEB.unique() ;
  cout<<"Number of EB strips="<<dec<<stripListEB.size()<<endl ;
  if (writeToFiles_) {
    (*out_file_) <<std::endl ;
    for (itList = stripListEB.begin(); itList != stripListEB.end(); itList++ ) {
      (*out_file_) <<"STRIP_EB "<<dec<<(*itList)<<endl ;
      (*out_file_) << hex << "0x" <<sliding_<<std::endl ;
      (*out_file_) <<"0" <<std::endl ;
      (*out_file_) <<"0x"<<stripMapEBsintheta[(*itList)] << " 0x" << SFGVB_lut_ <<std::endl ;
    }
  }

  // Endcap
  stripListEE.sort() ;
  stripListEE.unique() ;
  cout<<"Number of EE strips="<<dec<<stripListEE.size()<<endl ;
  if (writeToFiles_) {
    (*out_file_) <<std::endl ;
    for (itList = stripListEE.begin(); itList != stripListEE.end(); itList++ ) {
      (*out_file_) <<"STRIP_EE "<<dec<<(*itList)<<endl ;
      (*out_file_) << hex << "0x" <<sliding_<<std::endl ;
      (*out_file_) <<" 0" << std::endl ;
      (*out_file_)<<hex<<"0x"<<threshold<<" 0x"<<lut_strip<<std::endl ;  
    }
  }


  // loop on towers and associate them with default values //

  // Barrel
  towerListEB.sort() ;
  towerListEB.unique() ;
  cout<<"Number of EB towers="<<dec<<towerListEB.size()<<endl ;
  if (writeToFiles_) {
    (*out_file_) <<std::endl ;
    for (itList = towerListEB.begin(); itList != towerListEB.end(); itList++ ) {
      (*out_file_) <<"TOWER_EB "<<dec<<(*itList)<<endl ;
      (*out_file_) <<" 0\n 0\n" ;
      (*out_file_) <<" " <<  SFGVB_SpikeKillingThreshold_ << std::endl; //modif-alex
    }
  }

  // Endcap
  towerListEE.sort() ;
  towerListEE.unique() ;
  cout<<"Number of EE towers="<<dec<<towerListEE.size()<<endl ;
  if (writeToFiles_) {
    (*out_file_) <<std::endl ;
    for (itList = towerListEE.begin(); itList != towerListEE.end(); itList++ ) {
      (*out_file_) <<"TOWER_EE "<<dec<<(*itList)<<endl ;
      if (newLUT) (*out_file_) <<" 1\n" ;
      else  (*out_file_) <<" 0\n" ;
      (*out_file_)<<hex<<"0x"<<lut_tower<<std::endl ;
    }
  }




  // store control histos //
  ICEB->Write() ;
  tpgFactorEB->Write() ;
  ICEEPlus->Write() ;
  tpgFactorEEPlus->Write() ;  
  ICEEMinus->Write() ;
  tpgFactorEEMinus->Write() ;
  IC->Write() ;
  tpgFactor->Write() ;
  hshapeEB->Write() ;
  hshapeEE->Write() ;
  ntuple->Write() ;
  ntupleSpike->Write() ;
  saving.Close () ;
}

void EcalTPGParamBuilder::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1941 of file EcalTPGParamBuilder.cc.

References gather_cfg::cout, create_header(), TauDecayModes::dec, DetId::Ecal, EcalBarrel, EcalEndcap, Et_sat_EB_, Et_sat_EE_, FG_highRatio_EB_, FG_highThreshold_EB_, FG_lowRatio_EB_, FG_lowThreshold_EB_, FG_Threshold_EE_, DetId::rawId(), TTF_highThreshold_EB_, TTF_highThreshold_EE_, TTF_lowThreshold_EB_, TTF_lowThreshold_EE_, and writeToFiles_.

{
  using namespace edm;
  using namespace std;

  std::cout<<"we are in beginJob"<<endl;

  create_header() ; 

  DetId eb(DetId::Ecal,EcalBarrel) ;
  DetId ee(DetId::Ecal,EcalEndcap) ;

  if (writeToFiles_) {
    (*out_file_)<<"PHYSICS_EB "<<dec<<eb.rawId()<<std::endl ;
    (*out_file_)<<Et_sat_EB_<<" "<<TTF_lowThreshold_EB_<<" "<<TTF_highThreshold_EB_<<std::endl ;
    (*out_file_)<<FG_lowThreshold_EB_<<" "<<FG_highThreshold_EB_<<" "
          <<FG_lowRatio_EB_<<" "<<FG_highRatio_EB_<<std::endl ;
    //(*out_file_) << SFGVB_SpikeKillingThreshold_ << std::endl; //modif-alex02/02/2011
    (*out_file_) <<std::endl ;

    (*out_file_)<<"PHYSICS_EE "<<dec<<ee.rawId()<<std::endl ;
    (*out_file_)<<Et_sat_EE_<<" "<<TTF_lowThreshold_EE_<<" "<<TTF_highThreshold_EE_<<std::endl ;
    (*out_file_)<<FG_Threshold_EE_<<" "<<-1<<" "
          <<-1<<" "<<-1<<std::endl ;
    (*out_file_) <<std::endl ;
  }

}

bool EcalTPGParamBuilder::checkIfOK

(

const EcalPedestals::Item &

item

)

Definition at line 332 of file EcalTPGParamBuilder.cc.

References EcalPedestal::mean_x1, EcalPedestal::mean_x12, EcalPedestal::mean_x6, query::result, EcalPedestal::rms_x1, EcalPedestal::rms_x12, and EcalPedestal::rms_x6.

{
  bool result=true;
  if( item.mean_x1 <150. || item.mean_x1 >250) result=false;
  if( item.mean_x6 <150. || item.mean_x6 >250) result=false;
  if( item.mean_x12 <150. || item.mean_x12 >250) result=false;
  if( item.rms_x1 <0 || item.rms_x1 > 2) result=false;
  if( item.rms_x6 <0 || item.rms_x6 > 3) result=false;
  if( item.rms_x12 <0 || item.rms_x12 > 5) result=false;
  return result; 
}

void EcalTPGParamBuilder::computeFineGrainEBParameters ( uint &  lowRatio, uint &  highRatio, uint &  lowThreshold, uint &  highThreshold, uint &  lut  )

private

Definition at line 2370 of file EcalTPGParamBuilder.cc.

References Et_sat_EB_, FG_highRatio_EB_, FG_highThreshold_EB_, FG_lowRatio_EB_, FG_lowThreshold_EB_, and FG_lut_EB_.

Referenced by analyze().

{
  lowRatio = int(0x80*FG_lowRatio_EB_ + 0.5) ;
  if (lowRatio>0x7f) lowRatio = 0x7f ;
  highRatio = int(0x80*FG_highRatio_EB_ + 0.5) ;
  if (highRatio>0x7f) highRatio = 0x7f ;
  
  // lsb at the stage of the FG calculation is:
  double lsb_FG = Et_sat_EB_/1024./4 ;
  lowThreshold = int(FG_lowThreshold_EB_/lsb_FG+0.5) ;
  if (lowThreshold>0xff) lowThreshold = 0xff ;
  highThreshold = int(FG_highThreshold_EB_/lsb_FG+0.5) ;
  if (highThreshold>0xff) highThreshold = 0xff ;

  // FG lut: FGVB response is LUT(adress) where adress is: 
  // bit3: maxof2/ET >= lowRatio, bit2: maxof2/ET >= highRatio, bit1: ET >= lowThreshold, bit0: ET >= highThreshold
  // FGVB =1 if jet-like (veto active), =0 if E.M.-like
  // the condition for jet-like is: ET>Threshold and  maxof2/ET < Ratio (only TT with enough energy are vetoed)

  // With the following lut, what matters is only max(TLow, Thigh) and max(Elow, Ehigh)
  // So, jet-like if maxof2/ettot<max(TLow, Thigh) && ettot >= max(Elow, Ehigh)
  if (FG_lut_EB_ == 0) lut = 0x0888 ; 
  else lut = FG_lut_EB_ ; // let's use the users value (hope he/she knows what he/she does!)
}

void EcalTPGParamBuilder::computeFineGrainEEParameters ( uint &  threshold, uint &  lut_strip, uint &  lut_tower  )

private

Definition at line 2396 of file EcalTPGParamBuilder.cc.

References Et_sat_EE_, FG_lut_strip_EE_, FG_lut_tower_EE_, and FG_Threshold_EE_.

Referenced by analyze().

{
  // lsb for EE:
  double lsb_FG = Et_sat_EE_/1024. ; // FIXME is it true????
  threshold = int(FG_Threshold_EE_/lsb_FG+0.5) ;
  lut_strip = FG_lut_strip_EE_  ;
  lut_tower = FG_lut_tower_EE_  ;
}

bool EcalTPGParamBuilder::computeLinearizerParam ( double  theta, double  gainRatio, double  calibCoeff, std::string  subdet, int &  mult, int &  shift  )

private

Definition at line 1972 of file EcalTPGParamBuilder.cc.

References gather_cfg::cout, Et_sat_EB_, Et_sat_EE_, edm::false, funct::sin(), sliding_, funct::true, xtal_LSB_EB_, and xtal_LSB_EE_.

Referenced by analyze().

{
  /*
    Linearization coefficient are determined in order to satisfy:
    tpg(ADC_sat) = 1024
    where: 
    tpg() is a model of the linearized tpg response on 10b 
    ADC_sat is the number of ADC count corresponding the Et_sat, the maximum scale of the transverse energy
    
    Since we have:
    Et_sat = xtal_LSB * ADC_sat * gainRatio * calibCoeff * sin(theta)
    and a simple model of tpg() being given by:
    tpg(X) = [ (X*mult) >> (shift+2) ] >> (sliding+shiftDet) 
    we must satisfy:
    [ (Et_sat/(xtal_LSB * gainRatio * calibCoeff * sin(theta)) * mult) >> (shift+2) ] >> (sliding+shiftDet) = 1024 
    that is:
    mult = 1024/Et_sat * xtal_LSB * gainRatio * calibCoeff * sin(theta) * 2^-(sliding+shiftDet+2) * 2^-shift
    mult = factor * 2^-shift
  */

  // case barrel:
  int shiftDet = 2 ; //fixed, due to FE FENIX TCP format
  double ratio = xtal_LSB_EB_/Et_sat_EB_ ;
  // case endcap:
  if (subdet=="EE") {
    shiftDet = 2 ; //applied in TCC-EE and not in FE FENIX TCP... This parameters is setable in the TCC-EE
    //shiftDet = 0 ; //was like this before with FE bug
    ratio = xtal_LSB_EE_/Et_sat_EE_ ;
  }



  double factor = 1024 * ratio * gainRatio * calibCoeff * sin(theta) * (1 << (sliding_ + shiftDet + 2)) ;
  // Let's try first with shift = 0 (trivial solution)
  mult = (int)(factor+0.5) ; 
  for (shift = 0 ; shift<15 ; shift++) {
    if (mult>=128  && mult<256) return true ;
    factor *= 2 ; 
    mult = (int)(factor+0.5) ;
  }
  std::cout << "too bad we did not manage to calculate the factor for calib="<<calibCoeff<<endl;
  return false ;
}

void EcalTPGParamBuilder::computeLUT ( int *  lut, std::string  det = "EB"  )

private

Definition at line 2244 of file EcalTPGParamBuilder.cc.

References Et_sat_EB_, Et_sat_EE_, i, j, LUT_constant_EB_, LUT_constant_EE_, LUT_noise_EB_, LUT_noise_EE_, LUT_option_, LUT_stochastic_EB_, LUT_stochastic_EE_, LUT_threshold_EB_, LUT_threshold_EE_, oneOverEtResolEt(), TTF_highThreshold_EB_, TTF_highThreshold_EE_, TTF_lowThreshold_EB_, and TTF_lowThreshold_EE_.

Referenced by analyze().

{
  double Et_sat = Et_sat_EB_ ;
  double LUT_threshold = LUT_threshold_EB_ ;
  double LUT_stochastic = LUT_stochastic_EB_ ;
  double LUT_noise = LUT_noise_EB_ ;
  double LUT_constant = LUT_constant_EB_ ;
  double TTF_lowThreshold = TTF_lowThreshold_EB_ ;
  double TTF_highThreshold = TTF_highThreshold_EB_ ;
  if (det == "EE") {
    Et_sat = Et_sat_EE_ ;
    LUT_threshold = LUT_threshold_EE_ ;
    LUT_stochastic = LUT_stochastic_EE_ ;
    LUT_noise = LUT_noise_EE_ ;
    LUT_constant = LUT_constant_EE_ ;
    TTF_lowThreshold = TTF_lowThreshold_EE_ ;
    TTF_highThreshold = TTF_highThreshold_EE_ ;
  }

  // initialisation with identity
  for (int i=0 ; i<1024 ; i++) {
    lut[i] = i ;
    if (lut[i]>0xff) lut[i] = 0xff ;
  }

  // case linear LUT
  if (LUT_option_ == "Linear") {
    int mylut = 0 ;
    for (int i=0 ; i<1024 ; i++) {
      lut[i] = mylut ;
      if ((i+1)%4 == 0 ) mylut++ ; 
      //if ((i+1)%8 == 0 ) mylut++ ;//modif-alex 16/12/2010 LSB==500MeV ONLY USED FOR BEAMV4 key 
    }
  }

  // case LUT following Ecal resolution
  if (LUT_option_ == "EcalResolution") {
    TF1 * func = new TF1("func",oneOverEtResolEt, 0., Et_sat,3) ;
    func->SetParameters(LUT_stochastic, LUT_noise, LUT_constant) ;
    double norm = func->Integral(0., Et_sat) ;
    for (int i=0 ; i<1024 ; i++) {   
      double Et = i*Et_sat/1024. ;
      lut[i] =  int(0xff*func->Integral(0., Et)/norm + 0.5) ;
    }
  }

  // Now, add TTF thresholds to LUT and apply LUT threshold if needed
  for (int j=0 ; j<1024 ; j++) {
    double Et_GeV = Et_sat/1024*(j+0.5) ;
    if (Et_GeV <= LUT_threshold) lut[j] = 0 ; // LUT threshold
    int ttf = 0x0 ;    
    if (Et_GeV >= TTF_highThreshold) ttf = 3 ;
    if (Et_GeV >= TTF_lowThreshold && Et_GeV < TTF_highThreshold) ttf = 1 ;
    ttf = ttf << 8 ;
    lut[j] += ttf ;
  }

}

std::vector< unsigned int > EcalTPGParamBuilder::computeWeights ( EcalShapeBase &  shape, TH1F *  histo  )

private

Definition at line 2128 of file EcalTPGParamBuilder.cc.

References complement2_, prof2calltree::count, gather_cfg::cout, cmsHarvester::index, bookConverter::max, min(), nSample_, funct::pow(), compare_using_db::sample, sampleMax_, cond::rpcobgas::time, HLT_25ns14e33_v1_cff::timeMax, EcalShapeBase::timeOfMax(), EcalShapeBase::timeOfThr(), uncodeWeight(), histoStyle::weight, weight_timeShift_, and weight_unbias_recovery_.

Referenced by analyze().

{
  std::cout<<"Computing Weights..."<<std::endl ;
  double timeMax = shape.timeOfMax() - shape.timeOfThr() ; // timeMax w.r.t begining of pulse
  double max = shape(timeMax) ;

  double sumf = 0. ;
  double sumf2 = 0. ;
  for (unsigned int sample = 0 ; sample<nSample_ ; sample++) {
    double time = timeMax - ((double)sampleMax_-(double)sample)*25. ;
    time -= weight_timeShift_ ;
    sumf += shape(time)/max ;
    sumf2 += shape(time)/max * shape(time)/max ;
    for (int subtime = 0 ; subtime<25 ; subtime++) histo->Fill(float(sample*25. + subtime)/25., shape(time+subtime)) ;
  }
  double lambda = 1./(sumf2-sumf*sumf/nSample_) ;
  double gamma = -lambda*sumf/nSample_ ;
  double * weight = new double[nSample_] ;
  for (unsigned int sample = 0 ; sample<nSample_ ; sample++) {
    double time = timeMax - ((double)sampleMax_-(double)sample)*25. ;
    time -= weight_timeShift_ ;
    weight[sample] = lambda*shape(time)/max + gamma ;
  }


  int * iweight = new int[nSample_] ;
  for (unsigned int sample = 0 ; sample<nSample_ ; sample++)   iweight[sample] = uncodeWeight(weight[sample], complement2_) ;

  // Let's check:  
  int isumw  = 0 ;  
  for (unsigned int sample = 0 ; sample<nSample_ ; sample++) isumw  += iweight[sample] ;
  unsigned int imax = (unsigned int)(pow(2.,int(complement2_))-1) ;
  isumw = (isumw & imax ) ;

  double ampl = 0. ;
  double sumw = 0. ;
  for (unsigned int sample = 0 ; sample<nSample_ ; sample++) {
     double time = timeMax - ((double)sampleMax_-(double)sample)*25. ;
     time -= weight_timeShift_ ;
     ampl += weight[sample]*shape(time) ;
     sumw += weight[sample] ;
     std::cout<<"weight="<<weight[sample]<<" shape="<<shape(time)<<std::endl ;
  }
  std::cout<<"Weights: sum="<<isumw<<" in float ="<<uncodeWeight(isumw, complement2_)<<" sum of floats ="<<sumw<<std::endl ;
  std::cout<<"Weights: sum (weight*shape) = "<<ampl<<std::endl ;



  // Let's correct for bias if any
  if (weight_unbias_recovery_) {
    int count = 0 ;
    while (isumw != 0 && count<10) {
      double min = 99. ;
      unsigned int index = 0 ;
      if ( (isumw & (1<<(complement2_-1))) != 0) {
    // add 1:
    std::cout<<"Correcting for bias: adding 1"<<std::endl ;
    for (unsigned int sample = 0 ; sample<nSample_ ; sample++) {
      int new_iweight = iweight[sample]+1 ; 
      double new_weight = uncodeWeight(new_iweight, complement2_) ;
      if (fabs(new_weight-weight[sample])<min) {
        min = fabs(new_weight-weight[sample]) ;
        index = sample ;
      }
    }
    iweight[index] ++ ; 
      } else {
    // Sub 1:
    std::cout<<"Correcting for bias: subtracting 1"<<std::endl ;
    for (unsigned int sample = 0 ; sample<nSample_ ; sample++) {
      int new_iweight = iweight[sample]-1 ;    
      double new_weight = uncodeWeight(new_iweight, complement2_) ;
      if (fabs(new_weight-weight[sample])<min) {
        min = fabs(new_weight-weight[sample]) ;
        index = sample ;
      }
    }
    iweight[index] -- ; 
      } 
      isumw  = 0 ;  
      for (unsigned int sample = 0 ; sample<nSample_ ; sample++) isumw  += iweight[sample] ; 
      imax = (unsigned int)(pow(2.,int(complement2_))-1) ;
      isumw = (isumw & imax ) ;
      std::cout<<"Correcting weight number: "<<index<<" sum weights = "<<isumw<<std::endl ;
      count ++ ;
    }
  }
  
  // let's check again
  isumw  = 0 ;  
  for (unsigned int sample = 0 ; sample<nSample_ ; sample++) isumw  += iweight[sample] ;
  imax = (unsigned int)(pow(2.,int(complement2_))-1) ;
  isumw = (isumw & imax ) ;
  ampl = 0. ;
  for (unsigned int sample = 0 ; sample<nSample_ ; sample++) {
     double time = timeMax - ((double)sampleMax_-(double)sample)*25. ;
     time -= weight_timeShift_ ;
     double new_weight = uncodeWeight(iweight[sample], complement2_) ;
     sumw += uncodeWeight(iweight[sample], complement2_) ;
     ampl += new_weight*shape(time) ;
     std::cout<<"weight unbiased after integer conversion="<<new_weight<<" shape="<<shape(time)<<std::endl ;
  }
  std::cout<<"Weights: sum="<<isumw<<" in float ="<<uncodeWeight(isumw, complement2_)<<" sum of floats ="<<sumw<<std::endl ;
  std::cout<<"Weights: sum (weight*shape) = "<<ampl<<std::endl ;



  std::vector<unsigned int> theWeights ;
  for (unsigned int sample = 0 ; sample<nSample_ ; sample++) theWeights.push_back(iweight[sample]) ;
  std::cout<<std::endl ;

  delete[] weight ;
  delete[] iweight ;
  return theWeights ;
}

void EcalTPGParamBuilder::create_header ( )

private

Definition at line 2016 of file EcalTPGParamBuilder.cc.

References mathSSE::return(), and writeToFiles_.

Referenced by beginJob().

{
  if (!writeToFiles_) return ;
  (*out_file_) <<"COMMENT put your comments here"<<std::endl ; 

  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT           physics EB structure"<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  (*out_file_) <<"COMMENT  EtSaturation (GeV), ttf_threshold_Low (GeV), ttf_threshold_High (GeV)"<<std::endl ;
  (*out_file_) <<"COMMENT  FG_lowThreshold (GeV), FG_highThreshold (GeV), FG_lowRatio, FG_highRatio"<<std::endl ;
  //(*out_file_) <<"COMMENT  SFGVB_SpikeKillingThreshold (GeV)"<<std::endl ; //modif-alex-02/02/2011
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;

  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT           physics EE structure"<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  (*out_file_) <<"COMMENT  EtSaturation (GeV), ttf_threshold_Low (GeV), ttf_threshold_High (GeV)"<<std::endl ;
  (*out_file_) <<"COMMENT  FG_Threshold (GeV), dummy, dummy, dummy"<<std::endl ;
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;

  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT           crystal structure (same for EB and EE)"<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  (*out_file_) <<"COMMENT  ped, mult, shift [gain12]"<<std::endl ;
  (*out_file_) <<"COMMENT  ped, mult, shift [gain6]"<<std::endl ;
  (*out_file_) <<"COMMENT  ped, mult, shift [gain1]"<<std::endl ;
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;

  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT           strip EB structure"<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  (*out_file_) <<"COMMENT  sliding_window"<<std::endl ;
  (*out_file_) <<"COMMENT  weightGroupId"<<std::endl ;
  (*out_file_) <<"COMMENT  threshold_sfg lut_sfg"<<std::endl ;
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;

  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT           strip EE structure"<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  (*out_file_) <<"COMMENT  sliding_window"<<std::endl ;
  (*out_file_) <<"COMMENT  weightGroupId"<<std::endl ;
  (*out_file_) <<"COMMENT  threshold_fg lut_fg"<<std::endl ;
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;

  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT           tower EB structure"<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  (*out_file_) <<"COMMENT  LUTGroupId"<<std::endl ;
  (*out_file_) <<"COMMENT  FgGroupId"<<std::endl ;
  (*out_file_) <<"COMMENT  spike_killing_threshold"<<std::endl ;//modif alex
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;

  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT           tower EE structure"<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  (*out_file_) <<"COMMENT  LUTGroupId"<<std::endl ;
  (*out_file_) <<"COMMENT  tower_lut_fg"<<std::endl ;
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT           Weight structure"<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  (*out_file_) <<"COMMENT  weightGroupId"<<std::endl ;
  (*out_file_) <<"COMMENT  w0, w1, w2, w3, w4"<<std::endl ;
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;

  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT           lut structure"<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  (*out_file_) <<"COMMENT  LUTGroupId"<<std::endl ;
  (*out_file_) <<"COMMENT  LUT[1-1024]"<<std::endl ;
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;

  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT           fg EB structure"<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;
  (*out_file_) <<"COMMENT  FgGroupId"<<std::endl ;
  (*out_file_) <<"COMMENT  el, eh, tl, th, lut_fg"<<std::endl ;
  (*out_file_) <<"COMMENT ================================="<<std::endl ;
  (*out_file_) <<"COMMENT"<<std::endl ;

  (*out_file_) <<std::endl ;
}

void EcalTPGParamBuilder::getCoeff ( coeffStruc &  coeff, const EcalIntercalibConstantMap &  calibMap, uint  rawId  )

private

Referenced by analyze().

void EcalTPGParamBuilder::getCoeff ( coeffStruc &  coeff, const EcalGainRatioMap &  gainMap, uint  rawId  )

private

void EcalTPGParamBuilder::getCoeff ( coeffStruc &  coeff, const EcalPedestalsMap &  pedMap, uint  rawId  )

private

void EcalTPGParamBuilder::getCoeff ( coeffStruc &  coeff, const std::map< EcalLogicID, MonPedestalsDat > &  pedMap, const EcalLogicID &  logicId  )

private

Definition at line 2352 of file EcalTPGParamBuilder.cc.

References gather_cfg::cout, EcalLogicID::getID1(), EcalLogicID::getID2(), EcalLogicID::getID3(), MonPedestalsDat::getPedMeanG1(), MonPedestalsDat::getPedMeanG12(), MonPedestalsDat::getPedMeanG6(), and coeffStruc::pedestals_.

{
  // get current pedestal
  coeff.pedestals_[0] = 0 ;
  coeff.pedestals_[1] = 0 ;
  coeff.pedestals_[2] = 0 ;

  map<EcalLogicID, MonPedestalsDat>::const_iterator it =  pedMap.find(logicId);
  if (it != pedMap.end()) {
    MonPedestalsDat ped = it->second ;
    coeff.pedestals_[0] = int(ped.getPedMeanG12() + 0.5) ; 
    coeff.pedestals_[1] = int(ped.getPedMeanG6() + 0.5) ; 
    coeff.pedestals_[2] = int(ped.getPedMeanG1() + 0.5) ; 
  } 
  else std::cout<<"getCoeff: "<<logicId.getID1()<<", "<<logicId.getID2()<<", "<<logicId.getID3()
        <<" not found in map<EcalLogicID, MonPedestalsDat>"<<std::endl ;
}

std::pair< std::string, int > EcalTPGParamBuilder::getCrate ( int  tcc )

private

Definition at line 2468 of file EcalTPGParamBuilder.cc.

References AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by analyze().

                                                                {
  std::stringstream crate ;
  std::string pos ;
  int slot = 0 ;

  crate<<"S2D" ;  
  if (tcc>=40 && tcc<=42) {crate<<"02d" ; slot = 5 + (tcc-40)*6 ;}
  if (tcc>=43 && tcc<=45) {crate<<"03d" ; slot = 5 + (tcc-43)*6 ;}
  if (tcc>=37 && tcc<=39) {crate<<"04d" ; slot = 5 + (tcc-37)*6 ;}
  if (tcc>=52 && tcc<=54) {crate<<"06d" ; slot = 5 + (tcc-52)*6 ;}
  if (tcc>=46 && tcc<=48) {crate<<"07d" ; slot = 5 + (tcc-46)*6 ;}
  if (tcc>=49 && tcc<=51) {crate<<"08d" ; slot = 5 + (tcc-49)*6 ;}
  if (tcc>=58 && tcc<=60) {crate<<"02h" ; slot = 5 + (tcc-58)*6 ;}
  if (tcc>=61 && tcc<=63) {crate<<"03h" ; slot = 5 + (tcc-61)*6 ;}
  if (tcc>=55 && tcc<=57) {crate<<"04h" ; slot = 5 + (tcc-55)*6 ;}
  if (tcc>=70 && tcc<=72) {crate<<"06h" ; slot = 5 + (tcc-70)*6 ;}
  if (tcc>=64 && tcc<=66) {crate<<"07h" ; slot = 5 + (tcc-64)*6 ;}
  if (tcc>=67 && tcc<=69) {crate<<"08h" ; slot = 5 + (tcc-67)*6 ;}

  if (tcc>=76 && tcc<=81) {
    crate<<"02l" ; 
    if (tcc%2==0) slot = 2 + (tcc-76)*3 ;
    else slot = 4 + (tcc-77)*3 ;
  }
  if (tcc>=94 && tcc<=99) {
    crate<<"02l" ; 
    if (tcc%2==0) slot = 3 + (tcc-94)*3 ;
    else slot = 5 + (tcc-95)*3 ;
  }

  if (tcc>=22 && tcc<=27) {
    crate<<"03l" ; 
    if (tcc%2==0) slot = 2 + (tcc-22)*3 ;
    else slot = 4 + (tcc-23)*3 ;
  }
  if (tcc>=4 && tcc<=9) {
    crate<<"03l" ; 
    if (tcc%2==0) slot = 3 + (tcc-4)*3 ;
    else slot = 5 + (tcc-5)*3 ;
  }

  if (tcc>=82 && tcc<=87) {
    crate<<"07l" ; 
    if (tcc%2==0) slot = 2 + (tcc-82)*3 ;
    else slot = 4 + (tcc-83)*3 ;
  }
  if (tcc>=100 && tcc<=105) {
    crate<<"07l" ; 
    if (tcc%2==0) slot = 3 + (tcc-100)*3 ;
    else slot = 5 + (tcc-101)*3 ;
  }

  if (tcc>=28 && tcc<=33) {
    crate<<"08l" ; 
    if (tcc%2==0) slot = 2 + (tcc-28)*3 ;
    else slot = 4 + (tcc-29)*3 ;
  }
  if (tcc>=10 && tcc<=15) {
    crate<<"08l" ; 
    if (tcc%2==0) slot = 3 + (tcc-10)*3 ;
    else slot = 5 + (tcc-11)*3 ;
  }

  if (tcc==34) {crate<<"04l" ; slot = 2 ;}
  if (tcc==16) {crate<<"04l" ; slot = 3 ;}
  if (tcc==35) {crate<<"04l" ; slot = 4 ;}
  if (tcc==17) {crate<<"04l" ; slot = 5 ;}
  if (tcc==36) {crate<<"04l" ; slot = 8 ;}
  if (tcc==18) {crate<<"04l" ; slot = 9 ;}
  if (tcc==19) {crate<<"04l" ; slot = 10 ;}
  if (tcc==1)  {crate<<"04l" ; slot = 11 ;}
  if (tcc==20) {crate<<"04l" ; slot = 14 ;}
  if (tcc==2)  {crate<<"04l" ; slot = 15 ;}
  if (tcc==21) {crate<<"04l" ; slot = 16 ;}
  if (tcc==3)  {crate<<"04l" ; slot = 17 ;}

  if (tcc==88)  {crate<<"06l" ; slot = 2 ;}
  if (tcc==106) {crate<<"06l" ; slot = 3 ;}
  if (tcc==89)  {crate<<"06l" ; slot = 4 ;}
  if (tcc==107) {crate<<"06l" ; slot = 5 ;}
  if (tcc==90)  {crate<<"06l" ; slot = 8 ;}
  if (tcc==108) {crate<<"06l" ; slot = 9 ;}
  if (tcc==73)  {crate<<"06l" ; slot = 10 ;}
  if (tcc==91)  {crate<<"06l" ; slot = 11 ;}
  if (tcc==74)  {crate<<"06l" ; slot = 14 ;}
  if (tcc==92)  {crate<<"06l" ; slot = 15 ;}
  if (tcc==75)  {crate<<"06l" ; slot = 16 ;}
  if (tcc==93)  {crate<<"06l" ; slot = 17 ;}

  return std::pair< std::string, int > (crate.str(),slot) ;
}

std::string EcalTPGParamBuilder::getDet ( int  tcc )

private

Definition at line 2444 of file EcalTPGParamBuilder.cc.

Referenced by analyze().

                                            {
  std::stringstream sdet ;

  if (tcc>36 && tcc<55) sdet<<"EB-"<<tcc-36 ;
  else if (tcc>=55 && tcc<73) sdet<<"EB+"<<tcc-54 ;
  else if (tcc<=36) sdet<<"EE-" ;
  else sdet<<"EE+" ;
  
  if (tcc<=36 || tcc>=73) {
    if (tcc>=73) tcc-=72 ;
    if (tcc==1 || tcc==18 || tcc==19 || tcc==36)  sdet<<7 ;
    else if (tcc==2 || tcc==3 || tcc==20 || tcc==21)  sdet<<8 ;
    else if (tcc==4 || tcc==5 || tcc==22 || tcc==23)  sdet<<9 ;
    else if (tcc==6 || tcc==7 || tcc==24 || tcc==25)  sdet<<1 ;
    else if (tcc==8 || tcc==9 || tcc==26 || tcc==27)  sdet<<2 ;
    else if (tcc==10 || tcc==11 || tcc==28 || tcc==29)  sdet<<3 ;
    else if (tcc==12 || tcc==13 || tcc==30 || tcc==31)  sdet<<4 ;
    else if (tcc==14 || tcc==15 || tcc==32 || tcc==33)  sdet<<5 ;
    else if (tcc==16 || tcc==17 || tcc==34 || tcc==35)  sdet<<6 ;
  }

  return sdet.str() ;
}

int EcalTPGParamBuilder::getEtaSlice ( int  tccId, int  towerInTCC  )

private

Definition at line 344 of file EcalTPGParamBuilder.cc.

{
  int etaSlice = (towerInTCC-1)/4+1 ;
  // barrel
  if (tccId>36 || tccId<73) return etaSlice ;
  //endcap
  else {
    if (tccId >=1 && tccId <= 18) etaSlice += 21 ; // inner -
    if (tccId >=19 && tccId <= 36) etaSlice += 17 ; // outer -
    if (tccId >=91 && tccId <= 108) etaSlice += 21 ; // inner +
    if (tccId >=73 && tccId <= 90) etaSlice += 17 ; // outer +
  }
  return etaSlice ;
}

int EcalTPGParamBuilder::getGCTRegionEta ( int  tteta )

private

Definition at line 2437 of file EcalTPGParamBuilder.cc.

Referenced by analyze().

                                                 {
  int gcteta = 0;
  if(tteta>0) gcteta = (tteta-1)/4 + 11;
  else if(tteta<0) gcteta = (tteta+1)/4 + 10;
  return gcteta;
}

int EcalTPGParamBuilder::getGCTRegionPhi ( int  ttphi )

private

Definition at line 2429 of file EcalTPGParamBuilder.cc.

Referenced by analyze().

                                                 {
  int gctphi=0;
  gctphi = (ttphi+1)/4;
  if(ttphi<=2) gctphi=0;
  if(ttphi>=71) gctphi=0;
  return gctphi;
}

bool EcalTPGParamBuilder::realignBaseline ( linStruc &  lin, float  forceBase12  )

private

Definition at line 2405 of file EcalTPGParamBuilder.cc.

References newFWLiteAna::base, gather_cfg::cout, edm::false, i, linStruc::mult_, convertSQLiteXML::ok, linStruc::pedestal_, funct::pow(), and linStruc::shift_.

Referenced by analyze().

{
  bool ok(true) ;
  float base[3] = {forceBase12, float(lin.pedestal_[1]), float(lin.pedestal_[2])} ;
  for (int i=1 ; i<3 ; i++) {
    if (lin.mult_[i]>0) {
      base[i] = float(lin.pedestal_[i]) - 
    float(lin.mult_[0])/float(lin.mult_[i])*pow(2., -(lin.shift_[0]-lin.shift_[i]))*(lin.pedestal_[0]-base[0]) ;
    } else base[i] = 0 ;
  }

  for (int i=0 ; i<3 ; i++) {
    lin.pedestal_[i] = base[i] ;
    //cout<<lin.pedestal_[i]<<" "<<base[i]<<endl ;
    if (base[i]<0 || lin.pedestal_[i]>1000) {
      cout<<"WARNING: base= "<<base[i]<<", "<<lin.pedestal_[i]<<" for gainId[0-2]="<<i<<" ==> forcing at 0"<<endl ;
      lin.pedestal_[i] = 0 ; 
      ok = false ;
    }
  }
  return ok ;

}

int EcalTPGParamBuilder::uncodeWeight ( double  weight, int  complement2 = 7  )

private

Definition at line 2110 of file EcalTPGParamBuilder.cc.

References bookConverter::max, and funct::pow().

Referenced by computeWeights().

{
  int iweight ;
  unsigned int max = (unsigned int)(pow(2.,complement2)-1) ;
  if (weight>0) iweight=int((1<<6)*weight+0.5) ; // +0.5 for rounding pb
  else iweight= max - int(-weight*(1<<6)+0.5) +1 ;
  iweight = iweight & max ;
  return iweight ;
}

double EcalTPGParamBuilder::uncodeWeight ( int  iweight, int  complement2 = 7  )

private

Definition at line 2120 of file EcalTPGParamBuilder.cc.

References funct::pow(), and histoStyle::weight.

{
  double weight = double(iweight)/pow(2., 6.) ;
  // test if negative weight:
  if ( (iweight & (1<<(complement2-1))) != 0) weight = (double(iweight)-pow(2., complement2))/pow(2., 6.) ;
  return weight ;
}

Member Data Documentation

int EcalTPGParamBuilder::bst_conf_id_

private

Definition at line 142 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::btt_conf_id_

private

Definition at line 141 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::bxt_conf_id_

private

Definition at line 140 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

unsigned int EcalTPGParamBuilder::complement2_

private

Definition at line 98 of file EcalTPGParamBuilder.h.

Referenced by computeWeights().

EcalTPGDBApp* EcalTPGParamBuilder::db_

private

Definition at line 127 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

bool EcalTPGParamBuilder::DBEE_

private

Definition at line 130 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::del_conf_id_

private

Definition at line 139 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

std::map<int, std::vector<int> > EcalTPGParamBuilder::delays_EB_

private

Definition at line 120 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

std::map<int, std::vector<int> > EcalTPGParamBuilder::delays_EE_

private

Definition at line 122 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::Et_sat_EB_

private

Definition at line 92 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeFineGrainEBParameters(), computeLinearizerParam(), computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::Et_sat_EE_

private

Definition at line 92 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeFineGrainEEParameters(), computeLinearizerParam(), computeLUT(), and EcalTPGParamBuilder().

edm::ESHandle<EcalTrigTowerConstituentsMap> EcalTPGParamBuilder::eTTmap_

private

Definition at line 87 of file EcalTPGParamBuilder.h.

Referenced by analyze().

double EcalTPGParamBuilder::FG_highRatio_EB_

private

Definition at line 105 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeFineGrainEBParameters(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::FG_highThreshold_EB_

private

Definition at line 105 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeFineGrainEBParameters(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::FG_lowRatio_EB_

private

Definition at line 105 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeFineGrainEBParameters(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::FG_lowThreshold_EB_

private

Definition at line 105 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeFineGrainEBParameters(), and EcalTPGParamBuilder().

unsigned int EcalTPGParamBuilder::FG_lut_EB_

private

Definition at line 106 of file EcalTPGParamBuilder.h.

Referenced by computeFineGrainEBParameters(), and EcalTPGParamBuilder().

unsigned int EcalTPGParamBuilder::FG_lut_strip_EE_

private

Definition at line 108 of file EcalTPGParamBuilder.h.

Referenced by computeFineGrainEEParameters(), and EcalTPGParamBuilder().

unsigned int EcalTPGParamBuilder::FG_lut_tower_EE_

private

Definition at line 108 of file EcalTPGParamBuilder.h.

Referenced by computeFineGrainEEParameters(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::FG_Threshold_EE_

private

Definition at line 107 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeFineGrainEEParameters(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::fgr_conf_id_

private

Definition at line 135 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::forcedPedestalValue_

private

Definition at line 109 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

bool EcalTPGParamBuilder::forceEtaSlice_

private

Definition at line 110 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

std::ofstream* EcalTPGParamBuilder::geomFile_

private

Definition at line 126 of file EcalTPGParamBuilder.h.

Referenced by EcalTPGParamBuilder().

bool EcalTPGParamBuilder::H2_

private

Definition at line 113 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::lin_conf_id_

private

Definition at line 133 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::lut_conf_id_

private

Definition at line 134 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::LUT_constant_EB_

private

Definition at line 101 of file EcalTPGParamBuilder.h.

Referenced by computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::LUT_constant_EE_

private

Definition at line 102 of file EcalTPGParamBuilder.h.

Referenced by computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::LUT_noise_EB_

private

Definition at line 101 of file EcalTPGParamBuilder.h.

Referenced by computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::LUT_noise_EE_

private

Definition at line 102 of file EcalTPGParamBuilder.h.

Referenced by computeLUT(), and EcalTPGParamBuilder().

std::string EcalTPGParamBuilder::LUT_option_

private

Definition at line 99 of file EcalTPGParamBuilder.h.

Referenced by analyze(), computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::LUT_stochastic_EB_

private

Definition at line 101 of file EcalTPGParamBuilder.h.

Referenced by computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::LUT_stochastic_EE_

private

Definition at line 102 of file EcalTPGParamBuilder.h.

Referenced by computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::LUT_threshold_EB_

private

Definition at line 100 of file EcalTPGParamBuilder.h.

Referenced by computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::LUT_threshold_EE_

private

Definition at line 100 of file EcalTPGParamBuilder.h.

Referenced by computeLUT(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_bst

private

Definition at line 155 of file EcalTPGParamBuilder.h.

Referenced by EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_btt

private

Definition at line 154 of file EcalTPGParamBuilder.h.

Referenced by EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_bxt

private

Definition at line 153 of file EcalTPGParamBuilder.h.

Referenced by EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_del

private

Definition at line 152 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_fgr

private

Definition at line 149 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_lin

private

Definition at line 146 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_lut

private

Definition at line 147 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_ped

private

Definition at line 145 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_sli

private

Definition at line 150 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_spi

private

Definition at line 151 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::m_write_wei

private

Definition at line 148 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

unsigned int EcalTPGParamBuilder::nSample_

private

Definition at line 97 of file EcalTPGParamBuilder.h.

Referenced by computeWeights().

Char_t EcalTPGParamBuilder::ntupleCrate_[10]

private

Definition at line 159 of file EcalTPGParamBuilder.h.

Referenced by analyze().

Char_t EcalTPGParamBuilder::ntupleDet_[10]

private

Definition at line 158 of file EcalTPGParamBuilder.h.

Referenced by analyze().

Int_t* EcalTPGParamBuilder::ntupleInts_

private

Definition at line 157 of file EcalTPGParamBuilder.h.

Referenced by analyze().

std::ofstream* EcalTPGParamBuilder::out_file_

private

Definition at line 125 of file EcalTPGParamBuilder.h.

Referenced by analyze(), EcalTPGParamBuilder(), and ~EcalTPGParamBuilder().

int EcalTPGParamBuilder::ped_conf_id_

private

Definition at line 132 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

unsigned int EcalTPGParamBuilder::pedestal_offset_

private

Definition at line 111 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

std::map<int, std::vector<int> > EcalTPGParamBuilder::phases_EB_

private

Definition at line 121 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

std::map<int, std::vector<int> > EcalTPGParamBuilder::phases_EE_

private

Definition at line 123 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

unsigned int EcalTPGParamBuilder::sampleMax_

private

Definition at line 94 of file EcalTPGParamBuilder.h.

Referenced by computeWeights(), and EcalTPGParamBuilder().

unsigned int EcalTPGParamBuilder::SFGVB_lut_

private

Definition at line 111 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::SFGVB_SpikeKillingThreshold_

private

Definition at line 112 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

unsigned int EcalTPGParamBuilder::SFGVB_Threshold_

private

Definition at line 111 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::sli_conf_id_

private

Definition at line 136 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

unsigned int EcalTPGParamBuilder::sliding_

private

Definition at line 93 of file EcalTPGParamBuilder.h.

Referenced by analyze(), computeLinearizerParam(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::spi_conf_id_

private

Definition at line 138 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

std::string EcalTPGParamBuilder::tag_

private

Definition at line 143 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

const CaloSubdetectorGeometry* EcalTPGParamBuilder::theBarrelGeometry_

private

Definition at line 86 of file EcalTPGParamBuilder.h.

Referenced by analyze().

const CaloSubdetectorGeometry* EcalTPGParamBuilder::theEndcapGeometry_

private

Definition at line 85 of file EcalTPGParamBuilder.h.

Referenced by analyze().

const EcalElectronicsMapping* EcalTPGParamBuilder::theMapping_

private

Definition at line 88 of file EcalTPGParamBuilder.h.

Referenced by analyze().

std::string EcalTPGParamBuilder::TimingDelays_EB_

private

Definition at line 116 of file EcalTPGParamBuilder.h.

Referenced by EcalTPGParamBuilder().

std::string EcalTPGParamBuilder::TimingDelays_EE_

private

Definition at line 117 of file EcalTPGParamBuilder.h.

Referenced by EcalTPGParamBuilder().

std::string EcalTPGParamBuilder::TimingPhases_EB_

private

Definition at line 118 of file EcalTPGParamBuilder.h.

Referenced by EcalTPGParamBuilder().

std::string EcalTPGParamBuilder::TimingPhases_EE_

private

Definition at line 119 of file EcalTPGParamBuilder.h.

Referenced by EcalTPGParamBuilder().

double EcalTPGParamBuilder::TTF_highThreshold_EB_

private

Definition at line 103 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::TTF_highThreshold_EE_

private

Definition at line 104 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::TTF_lowThreshold_EB_

private

Definition at line 103 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeLUT(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::TTF_lowThreshold_EE_

private

Definition at line 104 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), computeLUT(), and EcalTPGParamBuilder().

bool EcalTPGParamBuilder::useInterCalibration_

private

Definition at line 113 of file EcalTPGParamBuilder.h.

Referenced by EcalTPGParamBuilder().

bool EcalTPGParamBuilder::useTransverseEnergy_

private

Definition at line 90 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::version_

private

Definition at line 144 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

int EcalTPGParamBuilder::wei_conf_id_

private

Definition at line 137 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

double EcalTPGParamBuilder::weight_timeShift_

private

Definition at line 95 of file EcalTPGParamBuilder.h.

Referenced by computeWeights(), and EcalTPGParamBuilder().

bool EcalTPGParamBuilder::weight_unbias_recovery_

private

Definition at line 96 of file EcalTPGParamBuilder.h.

Referenced by computeWeights(), and EcalTPGParamBuilder().

bool EcalTPGParamBuilder::writeToDB_

private

Definition at line 128 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and EcalTPGParamBuilder().

bool EcalTPGParamBuilder::writeToFiles_

private

Definition at line 129 of file EcalTPGParamBuilder.h.

Referenced by analyze(), beginJob(), create_header(), EcalTPGParamBuilder(), and ~EcalTPGParamBuilder().

double EcalTPGParamBuilder::xtal_LSB_EB_

private

Definition at line 91 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and computeLinearizerParam().

double EcalTPGParamBuilder::xtal_LSB_EE_

private

Definition at line 91 of file EcalTPGParamBuilder.h.

Referenced by analyze(), and computeLinearizerParam().