/data/refman/pasoursint/CMSSW_5_3_3/src/SimCalorimetry/EcalTrigPrimProducers/plugins/EcalTPCondAnalyzer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Class:      EcalTPDBAnalyzer
//
//
//
// Original Author:  Ursula Berthon
//         Created:  Wed Oct 15  11:38:38 CEST 2008
// $Id: EcalTPCondAnalyzer.cc,v 1.6 2011/10/24 15:27:54 ebecheva Exp $
//
//
//


// system include files
#include <memory>
#include <utility>
#include <iostream>
#include <fstream>

// user include files
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "CondFormats/DataRecord/interface/EcalTPGPhysicsConstRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGFineGrainEBGroupRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGLutGroupRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGLutIdMapRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGFineGrainEBIdMapRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGFineGrainTowerEERcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGLinearizationConstRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGPedestalsRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGSlidingWindowRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGWeightIdMapRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGWeightGroupRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGFineGrainStripEERcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGCrystalStatusRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGTowerStatusRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGStripStatusRcd.h"
#include "CondFormats/DataRecord/interface/EcalTPGSpikeRcd.h"

#include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"
#include "DataFormats/EcalDetId/interface/EcalTriggerElectronicsId.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"

#include "SimCalorimetry/EcalTrigPrimProducers/plugins/EcalTPCondAnalyzer.h"

EcalTPCondAnalyzer::EcalTPCondAnalyzer(const edm::ParameterSet&  iConfig)

{}

void EcalTPCondAnalyzer::beginRun(const edm::Run & run, edm::EventSetup const& evtSetup){
 // get geometry
  
  edm::ESHandle<CaloSubdetectorGeometry> theEndcapGeometry_handle, theBarrelGeometry_handle;
  evtSetup.get<EcalEndcapGeometryRecord>().get("EcalEndcap",theEndcapGeometry_handle);
  evtSetup.get<EcalBarrelGeometryRecord>().get("EcalBarrel",theBarrelGeometry_handle);
  theEndcapGeometry_ = &(*theEndcapGeometry_handle);
  theBarrelGeometry_ = &(*theBarrelGeometry_handle);

  cacheID_=this->getRecords(evtSetup);
}

void EcalTPCondAnalyzer::beginJob() 
{}


EcalTPCondAnalyzer::~EcalTPCondAnalyzer() {
}

unsigned long long  EcalTPCondAnalyzer::getRecords(edm::EventSetup const& setup) {
  //
  // get Eventsetup records and print them
  //
   printComment();

   edm::ESHandle<EcalTPGPhysicsConst> theEcalTPGPhysConst_handle;
  setup.get<EcalTPGPhysicsConstRcd>().get(theEcalTPGPhysConst_handle);
  const EcalTPGPhysicsConst * ecaltpPhysConst = theEcalTPGPhysConst_handle.product();
  printEcalTPGPhysicsConst(ecaltpPhysConst);
  // for EcalFenixStrip...

  // get parameter records for xtals
  edm::ESHandle<EcalTPGLinearizationConst> theEcalTPGLinearization_handle;
  setup.get<EcalTPGLinearizationConstRcd>().get(theEcalTPGLinearization_handle);
  const EcalTPGLinearizationConst * ecaltpLin = theEcalTPGLinearization_handle.product();

  edm::ESHandle<EcalTPGPedestals> theEcalTPGPedestals_handle;
  setup.get<EcalTPGPedestalsRcd>().get(theEcalTPGPedestals_handle);
  const EcalTPGPedestals * ecaltpPed = theEcalTPGPedestals_handle.product();
  printCRYSTAL(ecaltpPed,ecaltpLin );


  //weight
  edm::ESHandle<EcalTPGWeightIdMap> theEcalTPGWEightIdMap_handle;
  setup.get<EcalTPGWeightIdMapRcd>().get(theEcalTPGWEightIdMap_handle);
  const EcalTPGWeightIdMap * ecaltpgWeightMap = theEcalTPGWEightIdMap_handle.product();
  printWEIGHT(ecaltpgWeightMap);

   // .. and for EcalFenixTcp

  edm::ESHandle<EcalTPGFineGrainEBIdMap> theEcalTPGFineGrainEBIdMap_handle;
  setup.get<EcalTPGFineGrainEBIdMapRcd>().get(theEcalTPGFineGrainEBIdMap_handle);
  const EcalTPGFineGrainEBIdMap * ecaltpgFineGrainEB = theEcalTPGFineGrainEBIdMap_handle.product();
  printEcalTPGFineGrainEBIdMap(ecaltpgFineGrainEB);


  edm::ESHandle<EcalTPGLutIdMap> theEcalTPGLutIdMap_handle;
  setup.get<EcalTPGLutIdMapRcd>().get(theEcalTPGLutIdMap_handle);
  const EcalTPGLutIdMap * ecaltpgLut = theEcalTPGLutIdMap_handle.product();
  printEcalTPGLutIdMap(ecaltpgLut);

  //for strips
  edm::ESHandle<EcalTPGSlidingWindow> theEcalTPGSlidingWindow_handle;
  setup.get<EcalTPGSlidingWindowRcd>().get(theEcalTPGSlidingWindow_handle);
  const EcalTPGSlidingWindow * ecaltpgSlidW = theEcalTPGSlidingWindow_handle.product();
  edm::ESHandle<EcalTPGFineGrainStripEE> theEcalTPGFineGrainStripEE_handle;
  setup.get<EcalTPGFineGrainStripEERcd>().get(theEcalTPGFineGrainStripEE_handle);
  const EcalTPGFineGrainStripEE * ecaltpgFgStripEE = theEcalTPGFineGrainStripEE_handle.product();     
  edm::ESHandle<EcalTPGWeightGroup> theEcalTPGWEightGroup_handle;
  setup.get<EcalTPGWeightGroupRcd>().get(theEcalTPGWEightGroup_handle);
  const EcalTPGWeightGroup * ecaltpgWeightGroup = theEcalTPGWEightGroup_handle.product();
  printSTRIP(ecaltpgSlidW,ecaltpgWeightGroup,ecaltpgFgStripEE);

  // get parameter records for towers
  edm::ESHandle<EcalTPGLutGroup> theEcalTPGLutGroup_handle;
  setup.get<EcalTPGLutGroupRcd>().get(theEcalTPGLutGroup_handle);
  const EcalTPGLutGroup * ecaltpgLutGroup = theEcalTPGLutGroup_handle.product();

  edm::ESHandle<EcalTPGFineGrainEBGroup> theEcalTPGFineGrainEBGroup_handle;
  setup.get<EcalTPGFineGrainEBGroupRcd>().get(theEcalTPGFineGrainEBGroup_handle);
  const EcalTPGFineGrainEBGroup * ecaltpgFgEBGroup = theEcalTPGFineGrainEBGroup_handle.product();
  edm::ESHandle<EcalTPGSpike> theEcalTPGSpike_handle;
  setup.get<EcalTPGSpikeRcd>().get(theEcalTPGSpike_handle);
  const EcalTPGSpike * ecaltpgSpikeTh = theEcalTPGSpike_handle.product();

  printTOWEREB(ecaltpgSpikeTh, ecaltpgFgEBGroup,ecaltpgLutGroup);
  edm::ESHandle<EcalTPGFineGrainTowerEE> theEcalTPGFineGrainTowerEE_handle;
  setup.get<EcalTPGFineGrainTowerEERcd>().get(theEcalTPGFineGrainTowerEE_handle);
  const EcalTPGFineGrainTowerEE * ecaltpgFineGrainTowerEE = theEcalTPGFineGrainTowerEE_handle.product();

  printTOWEREE(ecaltpgFineGrainTowerEE, ecaltpgLutGroup);

  // get parameters for BadX
/*  edm::ESHandle<EcalTPGCrystalStatus> theEcalTPGCrystalStatus_handle;
  setup.get<EcalTPGCrystalStatusRcd>().get(theEcalTPGCrystalStatus_handle);
  const EcalTPGCrystalStatus * ecaltpgBadX = theEcalTPGCrystalStatus_handle.product();
  printBadX(ecaltpgBadX);

  // get parameters for BadTT
  edm::ESHandle<EcalTPGTowerStatus> theEcalTPGTowerStatus_handle;
  setup.get<EcalTPGTowerStatusRcd>().get(theEcalTPGTowerStatus_handle);
  const EcalTPGTowerStatus * ecaltpgBadTT = theEcalTPGTowerStatus_handle.product();
  printBadTT(ecaltpgBadTT);
*/  
  
  // get parameters for BadStrip
/*  edm::ESHandle<EcalTPGStripStatus> theEcalTPGStripStatus_handle;
  setup.get<EcalTPGStripStatusRcd>().get(theEcalTPGStripStatus_handle);
  const EcalTPGStripStatus * ecaltpgBadStrip = theEcalTPGStripStatus_handle.product();
  printBadStrip(ecaltpgBadStrip);
*/

  
  std::cout<<"EOF"<<std::endl;

  return setup.get<EcalTPGFineGrainTowerEERcd>().cacheIdentifier();
}

// ------------ method called to analyze the data  ------------
void
EcalTPCondAnalyzer::analyze(const edm::Event& iEvent, const  edm::EventSetup & iSetup){
}

void
EcalTPCondAnalyzer::endJob(){
}

void EcalTPCondAnalyzer::printEcalTPGPhysicsConst(const EcalTPGPhysicsConst *ecaltpgPhysConst) const {
  EcalTPGPhysicsConstMapIterator it;
  const EcalTPGPhysicsConstMap mymap=ecaltpgPhysConst->getMap();
  for (it=mymap.begin();it!=mymap.end();++it) {
    if (it==mymap.begin()) {
      std::cout<<"\nPHYSICS_EB "<<(*it).first<<std::endl;
    } else {
      std::cout<<"\nPHYSICS_EE "<<(*it).first<<std::endl;
    }
    std::cout<<(*it).second.EtSat<<" "<<(*it).second.ttf_threshold_Low<<" "<<(*it).second.ttf_threshold_High<<std::endl;
    std::cout<<(*it).second.FG_lowThreshold<<" "<<(*it).second.FG_highThreshold<<" "<<(*it).second.FG_lowRatio<<" "<<(*it).second.FG_highRatio<<std::endl;
  }
}

void  EcalTPCondAnalyzer::printSTRIP(const EcalTPGSlidingWindow *slWin,const EcalTPGWeightGroup *ecaltpgWeightGroup,const EcalTPGFineGrainStripEE * ecaltpgFgStripEE) const {
  // print STRIP information
  const EcalTPGSlidingWindowMap &slwinmap = slWin -> getMap();
  const EcalTPGFineGrainStripEEMap &fgstripEEmap= ecaltpgFgStripEE->getMap();
  EcalTPGSlidingWindowMapIterator it;
  const EcalTPGGroups::EcalTPGGroupsMap &gMap=ecaltpgWeightGroup->getMap();
  EcalTPGGroups::EcalTPGGroupsMapItr groupId;

  std::cout<<std::endl;
  for (int mysub=1;mysub<=2;++mysub) {
    std::cout<<std::endl;
    for (it=slwinmap.begin();it!=slwinmap.end();++it) {
      EcalTriggerElectronicsId elid((*it).first);
      groupId=gMap.find((*it).first);
      int subdet =elid.subdet();
      if (subdet==mysub) {
        if (subdet==1) {
          std::cout<<"STRIP_EB "<<std::dec<<(*it).first<<std::endl;
          std::cout << std::hex << "0x" <<(*it).second<<std::endl ;
          std::cout  <<"" <<(*groupId).second<< std::endl ; //weightgroupid
          EcalTPGFineGrainStripEEMapIterator it2=fgstripEEmap.find((*it).first);
          if (it2==fgstripEEmap.end()) {
            edm::LogWarning("EcalTPGCondAnalyzer") <<" could not find strip Id "<<(*it).first<<", given in sliding window, inside the EcalTPGFineGranStripEEMap!!!";
          } else {
            EcalTPGFineGrainStripEE::Item item=(*it2).second;
            std::cout<<std::hex<<"0x"<<item.threshold<<" 0x"<<item.lut<<std::endl ;  
          }
        }else if (subdet==2) {
          std::cout<<"STRIP_EE "<<std::dec<<(*it).first<<std::endl;
          std::cout << std::hex << "0x" <<(*it).second<<std::endl ;
          std::cout <<" " <<(*groupId).second<<std::endl ;//weightgroupid
          EcalTPGFineGrainStripEEMapIterator it2=fgstripEEmap.find((*it).first);
          if (it2==fgstripEEmap.end()) {
            edm::LogWarning("EcalTPGCondAnalyzer") <<" could not find strip Id "<<(*it).first<<", given in sliding window, inside the EcalTPGFineGranStripEEMap!!!";
          } else {
            EcalTPGFineGrainStripEE::Item item=(*it2).second;
            std::cout<<std::hex<<"0x"<<item.threshold<<" 0x"<<item.lut<<std::endl ;  
          }
        }
      }
    }
  }
}

void EcalTPCondAnalyzer::printWEIGHT(const EcalTPGWeightIdMap * ecaltpgWeightIdMap) const {
    
  std::cout<<std::endl;
    EcalTPGWeightIdMap::EcalTPGWeightMapItr it;
    uint32_t w0,w1,w2,w3,w4;
    const EcalTPGWeightIdMap::EcalTPGWeightMap map=ecaltpgWeightIdMap->getMap();
    for (it=map.begin();it!=map.end();++it) {
      std::cout <<"WEIGHT "<<(*it).first<<std::endl;
      (*it).second.getValues(w0,w1,w2,w3,w4);
      std::cout <<std::hex<<"0x"<<w0<<" 0x"<<w1<<" 0x"<<w2<<" 0x"<<w3<<" 0x"<<w4<<" "<<std::endl;
      std::cout<<std::endl;
      std::cout<<std::endl;
    }
}

void  EcalTPCondAnalyzer::printEcalTPGFineGrainEBIdMap(const EcalTPGFineGrainEBIdMap *ecaltpgFineGrainEB) const {
    EcalTPGFineGrainEBIdMap::EcalTPGFineGrainEBMapItr it;
    const EcalTPGFineGrainEBIdMap::EcalTPGFineGrainEBMap map=ecaltpgFineGrainEB->getMap();
    uint32_t ThresholdETLow, ThresholdETHigh, RatioLow, RatioHigh, LUT;

    //std::cout<<std::endl;
    for (it=map.begin();it!=map.end();++it) {
      std::cout <<"FG "<<(*it).first<<std::endl;
      (*it).second.getValues(ThresholdETLow, ThresholdETHigh, RatioLow, RatioHigh, LUT);
      std::cout <<std::hex<<"0x"<<ThresholdETLow<<" 0x"<<ThresholdETHigh<<" 0x"<<RatioLow<<" 0x"<<RatioHigh<<" 0x"<<LUT<<std::endl;
    }
}


void EcalTPCondAnalyzer::printEcalTPGLutIdMap(const EcalTPGLutIdMap *ecaltpgLut) const {
    EcalTPGLutIdMap::EcalTPGLutMapItr it;
    const EcalTPGLutIdMap::EcalTPGLutMap map=ecaltpgLut->getMap();

    std::cout<<std::endl;
    for (it=map.begin();it!=map.end();++it) {
      std::cout <<"LUT "<<(*it).first<<std::endl;
      const unsigned int * lut=(*it).second.getLut();
      for (unsigned int i=0;i<1024;++i)  std::cout <<std::hex<<"0x"<<*lut++<<std::endl;
    }
}

void EcalTPCondAnalyzer::printCRYSTAL(const EcalTPGPedestals * ecaltpPed, const EcalTPGLinearizationConst * ecaltpLin ) {

  std::cout<<std::endl;
  const EcalTPGPedestalsMap pedMap=ecaltpPed->getMap();
  const EcalTPGLinearizationConstMap linMap=ecaltpLin->getMap();

  const std::vector<DetId> & ebCells = theBarrelGeometry_->getValidDetIds(DetId::Ecal, EcalBarrel);

  std::cout<<"COMMENT ====== barrel crystals ====== "<<std::endl; 
  for (std::vector<DetId>::const_iterator it = ebCells.begin(); it != ebCells.end(); ++it) {
    EBDetId id(*it) ;
    std::cout <<"CRYSTAL "<<std::dec<<id.rawId()<<std::endl;
    const EcalTPGPedestal &ped=pedMap[id.rawId()];
    const EcalTPGLinearizationConstant &lin=linMap[id.rawId()];
    std::cout<<std::hex<<" 0x"<<ped.mean_x12<<" 0x"<<lin.mult_x12<<" 0x"<<lin.shift_x12<<std::endl;
    std::cout<<std::hex<<" 0x"<<ped.mean_x6 <<" 0x"<<lin.mult_x6 <<" 0x"<<lin.shift_x6<<std::endl;
    std::cout<<std::hex<<" 0x"<<ped.mean_x1 <<" 0x"<<lin.mult_x1 <<" 0x"<<lin.shift_x1<<std::endl;
  }

  const std::vector<DetId> & eeCells = theEndcapGeometry_->getValidDetIds(DetId::Ecal, EcalEndcap);
  std::cout<<"COMMENT ====== endcap crystals ====== "<<std::endl; 
  for (std::vector<DetId>::const_iterator it = eeCells.begin(); it != eeCells.end(); ++it) {
    EEDetId id(*it) ;
    std::cout <<"CRYSTAL "<<std::dec<<id.rawId()<<std::endl;
    const EcalTPGPedestal &ped=pedMap[id.rawId()];
    const EcalTPGLinearizationConstant &lin=linMap[id.rawId()];
    std::cout<<std::hex<<" 0x"<<ped.mean_x12<<" 0x"<<lin.mult_x12<<" 0x"<<lin.shift_x12<<std::endl;
    std::cout<<std::hex<<" 0x"<<ped.mean_x6 <<" 0x"<<lin.mult_x6 <<" 0x"<<lin.shift_x6<<std::endl;
    std::cout<<std::hex<<" 0x"<<ped.mean_x1 <<" 0x"<<lin.mult_x1 <<" 0x"<<lin.shift_x1<<std::endl;
  }
}
void EcalTPCondAnalyzer::printComment() const {
  std::cout<<"COMMENT put your comments here\n"<<
  "COMMENT =================================\n"<<
  "COMMENT           physics EB structure\n"<<
  "COMMENT\n"<<
  "COMMENT  EtSaturation (GeV), ttf_threshold_Low (GeV), ttf_threshold_High (GeV)\n"<<
  "COMMENT  FG_lowThreshold (GeV), FG_highThreshold (GeV), FG_lowRatio, FG_highRatio\n"<<
  "COMMENT =================================\n"<<
  "COMMENT\n"<<
  "COMMENT =================================\n"<<
  "COMMENT           physics EE structure\n"<<
  "COMMENT\n"<<
  "COMMENT  EtSaturation (GeV), ttf_threshold_Low (GeV), ttf_threshold_High (GeV)\n"<<
  "COMMENT  FG_Threshold (GeV), dummy, dummy, dummy\n"<<
  "COMMENT =================================\n"<<
  "COMMENT\n"<<
  "COMMENT =================================\n"<<
  "COMMENT           crystal structure (same for EB and EE)\n"<<
  "COMMENT\n"<<
  "COMMENT  ped, mult, shift [gain12]\n"<<
  "COMMENT  ped, mult, shift [gain6]\n"<<
  "COMMENT  ped, mult, shift [gain1]\n"<<
  "COMMENT =================================\n"<<
  "COMMENT\n"<<
  "COMMENT =================================\n"<<
  "COMMENT           strip EB structure\n"<<
  "COMMENT\n"<<
  "COMMENT  sliding_window\n"<<
  "COMMENT  weightGroupId\n"<<
  "COMMENT  threshold_sfg lut_sfg\n"<<
  "COMMENT =================================\n"<<
  "COMMENT\n"<<
  "COMMENT =================================\n"<<
  "COMMENT           strip EE structure\n"<<
  "COMMENT\n"<<
  "COMMENT  sliding_window\n"<<
  "COMMENT  weightGroupId\n"<<
  "COMMENT  threshold_fg lut_fg\n"<<
  "COMMENT =================================\n"<<
  "COMMENT\n"<<
  "COMMENT =================================\n"<<
  "COMMENT           tower EB structure\n"<<
  "COMMENT\n"<<
  "COMMENT  LUTGroupId\n"<<
  "COMMENT  FgGroupId\n"<<
  "COMMENT  spike_killing_threshold\n"<<
  "COMMENT =================================\n"<<
  "COMMENT\n"<<
  "COMMENT =================================\n"<<
  "COMMENT           tower EE structure\n"<<
  "COMMENT\n"<<
  "COMMENT  LUTGroupId\n"<<
  "COMMENT  tower_lut_fg\n"<<
  "COMMENT =================================\n"<<
  "COMMENT\n"<<
  "COMMENT =================================\n"<<
  "COMMENT           Weight structure\n"<<
  "COMMENT\n"<<
  "COMMENT  weightGroupId\n"<<
  "COMMENT  w0, w1, w2, w3, w4\n"<<
  "COMMENT =================================\n"<<
  "COMMENT\n"<<
  "COMMENT =================================\n"<<
  "COMMENT           lut structure\n"<<
  "COMMENT\n"<<
  "COMMENT  LUTGroupId\n"<<
  "COMMENT  LUT[1-1024]\n"<<
  "COMMENT =================================\n"<<
  "COMMENT\n"<<
  "COMMENT =================================\n"<<
  "COMMENT           fg EB structure\n"<<
  "COMMENT\n"<<
  "COMMENT  FgGroupId\n"<<
  "COMMENT  el, eh, tl, th, lut_fg\n"<<
  "COMMENT =================================\n"<<
  "COMMENT"<<std::endl;
}

/*void EcalTPCondAnalyzer::printTOWEREB(const EcalTPGFineGrainEBGroup *ecaltpgFgEBGroup,const EcalTPGLutGroup *ecaltpgLutGroup) const {

    const EcalTPGGroups::EcalTPGGroupsMap &lutMap=ecaltpgLutGroup->getMap();
    EcalTPGGroups::EcalTPGGroupsMapItr lutGroupId;
    const EcalTPGGroups::EcalTPGGroupsMap &fgMap=ecaltpgFgEBGroup->getMap();
    EcalTPGGroups::EcalTPGGroupsMapItr it;
    
    std::cout<<std::endl;
    for (it=fgMap.begin();it!=fgMap.end();++it) {
      std::cout <<"TOWER_EB "<<std::dec<<(*it).first<<std::endl;
      lutGroupId=lutMap.find((*it).first);
      std::cout <<" "<<(*it).second<<std::endl;
      std::cout <<" "<<(*lutGroupId).second<<std::endl;
    }
}
*/

void EcalTPCondAnalyzer::printTOWEREB(const EcalTPGSpike *ecaltpgSpikeTh, const EcalTPGFineGrainEBGroup *ecaltpgFgEBGroup,const EcalTPGLutGroup *ecaltpgLutGroup) const {
 
    const EcalTPGGroups::EcalTPGGroupsMap &lutMap=ecaltpgLutGroup->getMap();
    EcalTPGGroups::EcalTPGGroupsMapItr lutGroupId;
    const EcalTPGGroups::EcalTPGGroupsMap &fgMap=ecaltpgFgEBGroup->getMap();
    EcalTPGGroups::EcalTPGGroupsMapItr it;
    
    const EcalTPGSpike::EcalTPGSpikeMap spikeThMap = ecaltpgSpikeTh->getMap();
    EcalTPGSpike::EcalTPGSpikeMapIterator itSpikeTh;
    
    std::cout<<std::endl;
    for (it=fgMap.begin();it!=fgMap.end();++it) {
      std::cout <<"TOWER_EB "<<std::dec<<(*it).first<<std::endl;
      lutGroupId=lutMap.find((*it).first);
      itSpikeTh=spikeThMap.find((*it).first);
      std::cout <<" "<<(*it).second<<std::endl;
      std::cout <<" "<<(*lutGroupId).second<<std::endl;
      std::cout <<" "<<(*itSpikeTh).second<<std::endl; 
    }
}

void EcalTPCondAnalyzer::printTOWEREE(const EcalTPGFineGrainTowerEE *ecaltpgFineGrainTowerEE,const EcalTPGLutGroup *ecaltpgLutGroup) const {

    EcalTPGFineGrainTowerEEMapIterator it;
    const EcalTPGFineGrainTowerEEMap map=ecaltpgFineGrainTowerEE->getMap();
    const EcalTPGGroups::EcalTPGGroupsMap &lutMap=ecaltpgLutGroup->getMap();
    EcalTPGGroups::EcalTPGGroupsMapItr lutGroupId;

    std::cout<<std::endl;
    for (it=map.begin();it!=map.end();++it) {
      std::cout <<"TOWER_EE "<<std::dec<<(*it).first<<std::endl;
      lutGroupId=lutMap.find((*it).first);
      std::cout <<" "<<(*lutGroupId).second<<std::endl;
      std::cout <<std::hex<<"0x"<<(*it).second<<std::endl;
    }
}

/*void EcalTPCondAnalyzer::printBadX(const EcalTPGCrystalStatus *ecaltpgBadX) const {
  
  std::ofstream myfile;
  myfile.open("badXvalues.txt");
  
  const EcalTPGCrystalStatusMap & badXMap = ecaltpgBadX->getMap();

  const std::vector<DetId> & ebCells = theBarrelGeometry_->getValidDetIds(DetId::Ecal, EcalBarrel);

  myfile <<"COMMENT ====== barrel masked crystals ====== "<<std::endl; 
  myfile << "RawId         eta      phi   " << std::endl; 

  for (std::vector<DetId>::const_iterator it = ebCells.begin(); it != ebCells.end(); ++it) {

    EBDetId id(*it) ;
    const EcalTPGCrystalStatusCode &badXeb=badXMap[id.rawId()];
    // Print in the text file obly the masked crystals
    if (badXeb.getStatusCode() != 0){    
      myfile << "" << id.rawId() << "      " << id.ieta() << "      " << id.iphi() << std::endl;
    }
    
  }
  
  myfile << " " << std::endl;

  const std::vector<DetId> & eeCells = theEndcapGeometry_->getValidDetIds(DetId::Ecal, EcalEndcap);
  
  myfile <<"COMMENT ====== endcap masked crystals ====== "<<std::endl; 
  myfile << "RawId       x      y      z   " << std::endl; 
  
  for (std::vector<DetId>::const_iterator it = eeCells.begin(); it != eeCells.end(); ++it) {
    
    EEDetId id(*it) ;
    const EcalTPGCrystalStatusCode &badXee=badXMap[id.rawId()];
    // Print in the text file only the masked clystals
    if (badXee.getStatusCode() != 0){ 
      myfile << "" << id.rawId() << "   " << id.ix() << "      " << id.iy() << "     "<< id.zside() << std::endl;
    } 
  }
  
  myfile.close();
}
*/


/*void EcalTPCondAnalyzer::printBadTT(const EcalTPGTowerStatus *ecaltpgBadTT) const {
  std::ofstream myfilebadTT;
  myfilebadTT.open("badTTvalues.txt");
  int ieta = 0;
  int iphi = 0;
  
  const EcalTPGTowerStatusMap & badTTMap = ecaltpgBadTT -> getMap();
  EcalTPGTowerStatusMapIterator it;
  
  myfilebadTT <<"Barrel and endcap masked Trigger Towers"<<std::endl;
  myfilebadTT <<"RawId " << "     iphi " << "  ieta " << std::endl;
  myfilebadTT <<""<< std::endl;
  
  for (it=badTTMap.begin();it!=badTTMap.end();++it) {
    
    // Print in the text file only the masked barrel and endcap TTs
    if ((*it).second != 0){
    
      EcalTrigTowerDetId  ttId((*it).first);
      ieta = ttId.ieta();
      iphi = ttId.iphi();

      myfilebadTT <<""<< std::dec<<(*it).first << "  " << iphi << "     " << ieta << std::endl;
      
    }
  }
  
  
  myfilebadTT.close();

}
*/