EcalSignalGenerator.cc

Go to the documentation of this file.

#include "SimCalorimetry/EcalSimAlgos/interface/EcalSignalGenerator.h"


  template <>
  CaloSamples EcalSignalGenerator<EBDigitizerTraits>::samplesInPE(const DIGI & digi)
  {
    // calibration, for future reference:  (same block for all Ecal types)
    //EcalDetId cell = digi.id();
    //         const EcalCalibrations& calibrations=conditions->getEcalCalibrations(cell);
    //const EcalQIECoder* channelCoder = theConditions->getEcalCoder (cell);
    //const EcalQIEShape* channelShape = theConditions->getEcalShape (cell);
    //EcalCoderDb coder (*channelCoder, *channelShape);
    //CaloSamples result;
    //coder.adc2fC(digi, result);
    //fC2pe(result);

    DetId detId = digi.id();

    double Emax = fullScaleEnergy(detId); 
    double LSB[NGAINS+1]; 

    //double icalconst = findIntercalibConstant( detId );

    double icalconst = 1.;  // find the correct value.

    const EcalIntercalibConstantMCMap &icalMap = ical->getMap();
    EcalIntercalibConstantMCMap::const_iterator icalit = icalMap.find(detId);
    if( icalit!=icalMap.end() )
    {
        icalconst = (*icalit);
    }

    double peToA = peToAConversion ( detId ) ;

    const std::vector<float> gainRatios = GetGainRatios(detId);

    for( unsigned int igain ( 0 ); igain <= NGAINS ; ++igain ) 
      {
    LSB[igain] = 0.;
    if ( igain > 0 ) LSB[igain]= Emax/(MAXADC*gainRatios[igain]);
      }

    //    std::cout << " intercal, LSBs, egains " << icalconst << " " << LSB[0] << " " << LSB[1] << " " << gainRatios[0] << " " << gainRatios[1] << " " << Emax << std::endl;

    CaloSamples result(detId, digi.size());

    for(int isample = 0; isample<digi.size(); ++isample){

      int gainId = digi[isample].gainId();
      //int gainId = 1;

      if(gainId == 1) {
        result[isample] = float(digi[isample].adc())/1000./peToA; // special coding                   
      }
      else if(gainId > 1) {
        result[isample] = float(digi[isample].adc())*LSB[gainId-1]*icalconst/peToA;
      }  // gain = 0                                                                                  
      else { result[isample] = float(digi[isample].adc())*LSB[gainId]*icalconst/peToA;}

      // old version:
      //result[isample] = float(digi[isample].adc())*LSB[gainId]*icalconst/peToA;
    }

    //std::cout << " EcalSignalGenerator:EB noise input " << digi << std::endl;

    //std::cout << " converted noise sample " << std::endl;
    //for(int isample = 0; isample<digi.size(); ++isample){
    //  std::cout << " " << result[isample] ;
    //}
    //std::cout << std::endl;

    return result;
  }


  template <>
  CaloSamples EcalSignalGenerator<EEDigitizerTraits>::samplesInPE(const DIGI & digi)
  {
    // calibration, for future reference:  (same block for all Ecal types)
    //EcalDetId cell = digi.id();
    //         const EcalCalibrations& calibrations=conditions->getEcalCalibrations(cell);
    //const EcalQIECoder* channelCoder = theConditions->getEcalCoder (cell);
    //const EcalQIEShape* channelShape = theConditions->getEcalShape (cell);
    //EcalCoderDb coder (*channelCoder, *channelShape);
    //CaloSamples result;
    //coder.adc2fC(digi, result);
    //fC2pe(result);

    DetId detId = digi.id();

    double Emax = fullScaleEnergy(detId); 
    double LSB[NGAINS+1]; 

    double icalconst = 1.; //findIntercalibConstant( detId );

    const EcalIntercalibConstantMCMap &icalMap = ical->getMap();
    EcalIntercalibConstantMCMap::const_iterator icalit = icalMap.find(detId);
    if( icalit!=icalMap.end() )
    {
        icalconst = (*icalit);
    }


    double peToA = peToAConversion ( detId ) ;

    const std::vector<float> gainRatios = GetGainRatios(detId);

    for( unsigned int igain ( 0 ); igain <= NGAINS ; ++igain ) 
      {
    LSB[igain] = 0.;
    if ( igain > 0 ) LSB[igain]= Emax/(MAXADC*gainRatios[igain]);
      }

    //    std::cout << " intercal, LSBs, egains " << icalconst << " " << LSB[0] << " " << LSB[1] << " " << gainRatios[0] << " " << gainRatios[1] << " " << Emax << std::endl;

    CaloSamples result(detId, digi.size());

    for(int isample = 0; isample<digi.size(); ++isample){

      int gainId = digi[isample].gainId();
      //int gainId = 1;

      if(gainId == 1) {
        result[isample] = float(digi[isample].adc())/1000./peToA; // special coding                   
      }
      else if(gainId > 1) {
        result[isample] = float(digi[isample].adc())*LSB[gainId-1]*icalconst/peToA;
      }  // gain = 0                                                                                  
      else { result[isample] = float(digi[isample].adc())*LSB[gainId]*icalconst/peToA;}

      // old version
      //result[isample] = float(digi[isample].adc())*LSB[gainId]*icalconst/peToA;
    }

    //std::cout << " EcalSignalGenerator:EE noise input " << digi << std::endl;

    //std::cout << " converted noise sample " << std::endl;
    //for(int isample = 0; isample<digi.size(); ++isample){
    //  std::cout << " " << result[isample] ;
    // }
    //std::cout << std::endl;

    return result;
  }

  template <>
  CaloSamples EcalSignalGenerator<ESDigitizerTraits>::samplesInPE(const DIGI & digi)
  {
    // calibration, for future reference:  (same block for all Ecal types)
    //EcalDetId cell = digi.id();
    //         const EcalCalibrations& calibrations=conditions->getEcalCalibrations(cell);
    //const EcalQIECoder* channelCoder = theConditions->getEcalCoder (cell);
    //const EcalQIEShape* channelShape = theConditions->getEcalShape (cell);
    //EcalCoderDb coder (*channelCoder, *channelShape);
    //CaloSamples result;
    //coder.adc2fC(digi, result);
    //fC2pe(result);

    DetId detId = digi.id();

    double icalconst = 1.; //findIntercalibConstant( detId );

    const ESIntercalibConstantMap &icalMap = esmips->getMap();
    ESIntercalibConstantMap::const_iterator icalit = icalMap.find(detId);
    if( icalit!=icalMap.end() )
      {
    icalconst = double (*icalit);
      }

    CaloSamples result(detId, digi.size());

    for(int isample = 0; isample<digi.size(); ++isample){
      result[isample] = float(digi[isample].adc())/icalconst*ESMIPToGeV;
    }

    //std::cout << " EcalSignalGenerator:ES noise input " << digi << std::endl;

    //std::cout << " converted noise sample " << std::endl;
    //for(int isample = 0; isample<digi.size(); ++isample){
    //  std::cout << " " << result[isample] ;
    //}
    //std::cout << std::endl;

    return result;
  }