EcalBarrelDigisValidation Class Reference

#include <EcalBarrelDigisValidation.h>

Inheritance diagram for EcalBarrelDigisValidation:

Public Member Functions
void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override
	EcalBarrelDigisValidation (const edm::ParameterSet &ps)
	Constructor. More...
	~EcalBarrelDigisValidation () override
	Destructor. More...
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &c) override
	Analyze. More...
void	checkCalibrations (edm::EventSetup const &c)
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override

Private Types
typedef std::map< uint32_t, float, std::less< uint32_t > >	MapType

Private Attributes
double	barrelADCtoGeV_
edm::EDGetTokenT< EBDigiCollection >	EBdigiCollection_
double	endcapADCtoGeV_
std::map< int, double, std::less< int > >	gainConv_
MonitorElement *	meEBDigiADCAnalog_ [10]
MonitorElement *	meEBDigiADCg12_ [10]
MonitorElement *	meEBDigiADCg1_ [10]
MonitorElement *	meEBDigiADCg6_ [10]
MonitorElement *	meEBDigiADCGlobal_
MonitorElement *	meEBDigiADCgS_ [10]
MonitorElement *	meEBDigiGain_ [10]
MonitorElement *	meEBDigiMultiplicity_
MonitorElement *	meEBDigiOccupancy_
MonitorElement *	meEBMaximumgt100ADC_
MonitorElement *	meEBMaximumgt10ADC_
MonitorElement *	meEBnADCafterSwitch_
MonitorElement *	meEBPedestal_
std::string	outputFile_
bool	verbose_

Detailed Description

Definition at line 36 of file EcalBarrelDigisValidation.h.

Member Typedef Documentation

typedef std::map<uint32_t,float,std::less<uint32_t> > EcalBarrelDigisValidation::MapType

private

Definition at line 38 of file EcalBarrelDigisValidation.h.

Constructor & Destructor Documentation

EcalBarrelDigisValidation::EcalBarrelDigisValidation

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 16 of file EcalBarrelDigisValidation.cc.

References barrelADCtoGeV_, endcapADCtoGeV_, gainConv_, edm::ParameterSet::getUntrackedParameter(), mps_fire::i, meEBDigiADCAnalog_, meEBDigiADCg12_, meEBDigiADCg1_, meEBDigiADCg6_, meEBDigiADCGlobal_, meEBDigiADCgS_, meEBDigiGain_, meEBDigiMultiplicity_, meEBDigiOccupancy_, meEBMaximumgt100ADC_, meEBMaximumgt10ADC_, meEBnADCafterSwitch_, meEBPedestal_, and verbose_.

                                                                          :
  EBdigiCollection_(consumes<EBDigiCollection>(ps.getParameter<edm::InputTag>("EBdigiCollection")))
{
  
  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
    
                                                                                                                                          
  // get hold of back-end interface
                                                                                                                                          

  gainConv_[1] = 1.;
  gainConv_[2] = 2.;
  gainConv_[3] = 12.;
  gainConv_[0] = 12.;   // saturated channels
  barrelADCtoGeV_ = 0.035;
  endcapADCtoGeV_ = 0.06;
 
  meEBDigiOccupancy_ = nullptr;

  meEBDigiMultiplicity_ = nullptr;

  meEBDigiADCGlobal_ = nullptr;

  for (int i = 0; i < 10 ; i++ ) {
    meEBDigiADCAnalog_[i] = nullptr;
    meEBDigiADCgS_[i]  = nullptr;
    meEBDigiADCg1_[i]  = nullptr;
    meEBDigiADCg6_[i]  = nullptr;
    meEBDigiADCg12_[i] = nullptr;
    meEBDigiGain_[i] = nullptr;
  }

  meEBPedestal_ = nullptr;
                                 
  meEBMaximumgt100ADC_ = nullptr; 
                                 
  meEBMaximumgt10ADC_ = nullptr; 

  meEBnADCafterSwitch_ = nullptr;
 
}

EcalBarrelDigisValidation::~EcalBarrelDigisValidation ( )

override

Destructor.

Definition at line 59 of file EcalBarrelDigisValidation.cc.

                                                     {
 
}

Member Function Documentation

void EcalBarrelDigisValidation::analyze ( edm::Event const &  e, edm::EventSetup const &  c  )

overrideprotected

Analyze.

Definition at line 121 of file EcalBarrelDigisValidation.cc.

References EcalMGPASample::adc(), barrelADCtoGeV_, EBdigiCollection_, MonitorElement::Fill(), gainConv_, EcalMGPASample::gainId(), edm::Event::getByToken(), mps_fire::i, EBDataFrame::id(), EBDetId::ieta(), createfilelist::int, EBDetId::iphi(), edm::HandleBase::isValid(), LogDebug, EcalDataFrame::MAXSAMPLES, meEBDigiADCAnalog_, meEBDigiADCg12_, meEBDigiADCg1_, meEBDigiADCg6_, meEBDigiADCGlobal_, meEBDigiADCgS_, meEBDigiGain_, meEBDigiMultiplicity_, meEBDigiOccupancy_, meEBMaximumgt100ADC_, meEBMaximumgt10ADC_, meEBnADCafterSwitch_, meEBPedestal_, dataAnalyzerFineBiningParameters_cff::Pmax, edm::Handle< T >::product(), DetId::rawId(), edm::DataFrameContainer::reserve(), simplePhotonAnalyzer_cfi::sample, EcalDataFrame::size(), and edm::DataFrameContainer::size().

                                                                            {

  //LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();

  Handle<EBDigiCollection> EcalDigiEB;

  e.getByToken( EBdigiCollection_ , EcalDigiEB );

  //Return if no Barrel data 
  if( !EcalDigiEB.isValid() ) return;

  // BARREL

  // loop over Digis

  const EBDigiCollection * barrelDigi = EcalDigiEB.product () ;

  std::vector<double> ebAnalogSignal ;
  std::vector<double> ebADCCounts ;
  std::vector<double> ebADCGains ;
  ebAnalogSignal.reserve(EBDataFrame::MAXSAMPLES);
  ebADCCounts.reserve(EBDataFrame::MAXSAMPLES);
  ebADCGains.reserve(EBDataFrame::MAXSAMPLES);

  int nDigis = 0;
  
  for (unsigned int digis=0; digis<EcalDigiEB->size(); ++digis) 
    {
      
      EBDataFrame ebdf = (*barrelDigi)[digis];
      int nrSamples = ebdf.size();
      
      EBDetId ebid = ebdf.id () ;

      nDigis++;
      if (meEBDigiOccupancy_) meEBDigiOccupancy_->Fill( ebid.iphi(), ebid.ieta() );
      
      double Emax = 0. ;
      int Pmax = 0 ;
      double pedestalPreSample = 0.;
      double pedestalPreSampleAnalog = 0.;
      int countsAfterGainSwitch = -1;
      double higherGain = 1.;
      int higherGainSample = 0;
      
      for (int sample = 0 ; sample < nrSamples; ++sample) {
        ebAnalogSignal[sample] = 0.;
        ebADCCounts[sample] = 0.;
        ebADCGains[sample] = 0.;
      }
      
      for (int sample = 0 ; sample < nrSamples; ++sample) 
        {
      EcalMGPASample thisSample = ebdf[sample];
      
          ebADCCounts[sample] = (thisSample.adc());
          ebADCGains[sample]  = (thisSample.gainId());
          ebAnalogSignal[sample] = (ebADCCounts[sample]*gainConv_[(int)ebADCGains[sample]]*barrelADCtoGeV_);
      
          if (Emax < ebAnalogSignal[sample] ) {
            Emax = ebAnalogSignal[sample] ;
            Pmax = sample ;
          }
      
          if ( sample < 3 ) {
            pedestalPreSample += ebADCCounts[sample] ;
            pedestalPreSampleAnalog += ebADCCounts[sample]*gainConv_[(int)ebADCGains[sample]]*barrelADCtoGeV_ ;
          }
      
          if ( sample > 0 && ( ((ebADCGains[sample] > ebADCGains[sample-1]) && (ebADCGains[sample-1]!=0)) || (countsAfterGainSwitch<0 && ebADCGains[sample]==0)) ) {
            higherGain = ebADCGains[sample];
            higherGainSample = sample;
            countsAfterGainSwitch = 1;
          }
      
          if ( (higherGain > 1 && (higherGainSample != sample) && (ebADCGains[sample] == higherGain)) || (higherGain==3 && (higherGainSample != sample) && (ebADCGains[sample]==0)) || (higherGain==0 && (higherGainSample != sample) && ((ebADCGains[sample] == 3) || (ebADCGains[sample]==0))) ) { countsAfterGainSwitch++ ; }
        }
      
      pedestalPreSample /= 3. ; 
      pedestalPreSampleAnalog /= 3. ; 
      
      LogDebug("DigiInfo") << "Barrel Digi for EBDetId = " << ebid.rawId() << " eta,phi " << ebid.ieta() << " " << ebid.iphi() ;
      for ( int i = 0; i < 10 ; i++ ) {
    LogDebug("DigiInfo") << "sample " << i << " ADC = " << ebADCCounts[i] << " gain = " << ebADCGains[i] << " Analog = " << ebAnalogSignal[i];  
      }
      LogDebug("DigiInfo") << "Maximum energy = " << Emax << " in sample " << Pmax << " Pedestal from pre-sample = " << pedestalPreSampleAnalog;
      if ( countsAfterGainSwitch > 0 ) LogDebug("DigiInfo") << "Counts after switch " << countsAfterGainSwitch;
      
      if ( countsAfterGainSwitch > 0 && countsAfterGainSwitch < 5 ) {
    edm::LogWarning("DigiWarning") << "Wrong number of counts after gain switch before next switch! " << countsAfterGainSwitch ;
    for ( int i = 0; i < 10 ; i++ ) {
      edm::LogWarning("DigiWarning") << "sample " << i << " ADC = " << ebADCCounts[i] << " gain = " << ebADCGains[i] << " Analog = " << ebAnalogSignal[i];
    } 
      }
      
      for ( int i = 0 ; i < 10 ; i++ ) {
    if (meEBDigiADCGlobal_ && (Emax-pedestalPreSampleAnalog*gainConv_[(int)ebADCGains[Pmax]]) > 100.*barrelADCtoGeV_) meEBDigiADCGlobal_->Fill( i , ebAnalogSignal[i] ) ;
    if (meEBDigiADCAnalog_[i]) meEBDigiADCAnalog_[i]->Fill( ebAnalogSignal[i] ) ;
    
    if ( ebADCGains[i] == 0) {
      if (meEBDigiADCgS_[i]) meEBDigiADCgS_[i]->Fill( ebADCCounts[i] ) ;
    }
    else if ( ebADCGains[i] == 3 ) {
      if (meEBDigiADCg1_[i]) meEBDigiADCg1_[i]->Fill( ebADCCounts[i] ) ;
    }
    else if ( ebADCGains[i] == 2 ) {
      if (meEBDigiADCg6_[i]) meEBDigiADCg6_[i]->Fill( ebADCCounts[i] ) ;
    }
    else if ( ebADCGains[i] == 1 ) {
      if (meEBDigiADCg12_[i]) meEBDigiADCg12_[i]->Fill( ebADCCounts[i] ) ;
    }
    if (meEBDigiGain_[i]) meEBDigiGain_[i]->Fill( ebADCGains[i] ) ;
      }
      
      if (meEBPedestal_) meEBPedestal_->Fill ( pedestalPreSample ) ;
      if (meEBMaximumgt10ADC_ && (Emax-pedestalPreSampleAnalog*gainConv_[(int)ebADCGains[Pmax]]) > 10.*barrelADCtoGeV_) meEBMaximumgt10ADC_->Fill( Pmax ) ;
      if (meEBMaximumgt100ADC_ && (Emax-pedestalPreSampleAnalog*gainConv_[(int)ebADCGains[Pmax]]) > 100.*barrelADCtoGeV_) meEBMaximumgt100ADC_->Fill( Pmax ) ;
      if (meEBnADCafterSwitch_) meEBnADCafterSwitch_->Fill( countsAfterGainSwitch ) ;
      
    }
  
  if ( meEBDigiMultiplicity_ ) meEBDigiMultiplicity_->Fill(nDigis);
    
}

void EcalBarrelDigisValidation::bookHistograms ( DQMStore::IBooker &  i, edm::Run const &  , edm::EventSetup const &    )

override

Definition at line 63 of file EcalBarrelDigisValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), trackerHits::histo, mps_fire::i, meEBDigiADCAnalog_, meEBDigiADCg12_, meEBDigiADCg1_, meEBDigiADCg6_, meEBDigiADCGlobal_, meEBDigiADCgS_, meEBDigiGain_, meEBDigiMultiplicity_, meEBDigiOccupancy_, meEBMaximumgt100ADC_, meEBMaximumgt10ADC_, meEBnADCafterSwitch_, meEBPedestal_, and DQMStore::IBooker::setCurrentFolder().

                                                                                                           {

    Char_t histo[200];

    ibooker.setCurrentFolder("EcalDigisV/EcalDigiTask");
  
    sprintf (histo, "EcalDigiTask Barrel occupancy" ) ;
    meEBDigiOccupancy_ = ibooker.book2D(histo, histo, 360, 0., 360., 170, -85., 85.);

    sprintf (histo, "EcalDigiTask Barrel digis multiplicity" ) ;
    meEBDigiMultiplicity_ = ibooker.book1D(histo, histo, 612, 0., 61200);
  
    sprintf (histo, "EcalDigiTask Barrel global pulse shape" ) ;
    meEBDigiADCGlobal_ = ibooker.bookProfile(histo, histo, 10, 0, 10, 10000, 0., 1000.) ;
    
    for (int i = 0; i < 10 ; i++ ) {

      sprintf (histo, "EcalDigiTask Barrel analog pulse %02d", i+1) ;
      meEBDigiADCAnalog_[i] = ibooker.book1D(histo, histo, 4000, 0., 400.);

      sprintf (histo, "EcalDigiTask Barrel ADC pulse %02d Gain 0 - Saturated", i+1) ;
      meEBDigiADCgS_[i] = ibooker.book1D(histo, histo, 4096, -0.5, 4095.5);

      sprintf (histo, "EcalDigiTask Barrel ADC pulse %02d Gain 1", i+1) ;
      meEBDigiADCg1_[i] = ibooker.book1D(histo, histo, 4096, -0.5, 4095.5);

      sprintf (histo, "EcalDigiTask Barrel ADC pulse %02d Gain 6", i+1) ;
      meEBDigiADCg6_[i] = ibooker.book1D(histo, histo, 4096, -0.5, 4095.5);

      sprintf (histo, "EcalDigiTask Barrel ADC pulse %02d Gain 12", i+1) ;
      meEBDigiADCg12_[i] = ibooker.book1D(histo, histo, 4096, -0.5, 4095.5);

      sprintf (histo, "EcalDigiTask Barrel gain pulse %02d", i+1) ;
      meEBDigiGain_[i] = ibooker.book1D(histo, histo, 4, 0, 4);

    }
    
    sprintf (histo, "EcalDigiTask Barrel pedestal for pre-sample" ) ;
    meEBPedestal_ = ibooker.book1D(histo, histo, 4096, -0.5, 4095.5) ;

    sprintf (histo, "EcalDigiTask Barrel maximum position gt 100 ADC" ) ;
    meEBMaximumgt100ADC_ = ibooker.book1D(histo, histo, 10, 0., 10.) ;

    sprintf (histo, "EcalDigiTask Barrel maximum position gt 10 ADC" ) ;
    meEBMaximumgt10ADC_ = ibooker.book1D(histo, histo, 10, 0., 10.) ;

    sprintf (histo, "EcalDigiTask Barrel ADC counts after gain switch" ) ;
    meEBnADCafterSwitch_ = ibooker.book1D(histo, histo, 10, 0., 10.) ;

}

void EcalBarrelDigisValidation::checkCalibrations ( edm::EventSetup const &  c )

protected

Definition at line 246 of file EcalBarrelDigisValidation.cc.

References barrelADCtoGeV_, endcapADCtoGeV_, EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), gainConv_, edm::EventSetup::get(), EcalADCToGeVConstant::getEBValue(), EcalADCToGeVConstant::getEEValue(), LogDebug, and edm::ESHandle< T >::product().

Referenced by dqmBeginRun().

  {
    
  // ADC -> GeV Scale
  edm::ESHandle<EcalADCToGeVConstant> pAgc;
  eventSetup.get<EcalADCToGeVConstantRcd>().get(pAgc);
  const EcalADCToGeVConstant* agc = pAgc.product();
  
  EcalMGPAGainRatio * defaultRatios = new EcalMGPAGainRatio();

  gainConv_[1] = 1.;
  gainConv_[2] = defaultRatios->gain12Over6() ;
  gainConv_[3] = gainConv_[2]*(defaultRatios->gain6Over1()) ;
  gainConv_[0] = gainConv_[2]*(defaultRatios->gain6Over1()) ;

  LogDebug("EcalDigi") << " Gains conversions: " << "\n" << " g0 = " << gainConv_[0] << "\n" << " g1 = " << gainConv_[1] << "\n" << " g2 = " << gainConv_[2]  << "\n" << " g3 = " << gainConv_[3]; 

  delete defaultRatios;

  const double barrelADCtoGeV_  = agc->getEBValue();
  LogDebug("EcalDigi") << " Barrel GeV/ADC = " << barrelADCtoGeV_;
  const double endcapADCtoGeV_ = agc->getEEValue();
  LogDebug("EcalDigi") << " Endcap GeV/ADC = " << endcapADCtoGeV_;

}

void EcalBarrelDigisValidation::dqmBeginRun ( edm::Run const &  , edm::EventSetup const &  c  )

overrideprotected

Definition at line 114 of file EcalBarrelDigisValidation.cc.

References checkCalibrations().

                                                                                {

  checkCalibrations(c);

}

Member Data Documentation

double EcalBarrelDigisValidation::barrelADCtoGeV_

private

Definition at line 68 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), checkCalibrations(), and EcalBarrelDigisValidation().

edm::EDGetTokenT<EBDigiCollection> EcalBarrelDigisValidation::EBdigiCollection_

private

Definition at line 64 of file EcalBarrelDigisValidation.h.

Referenced by analyze().

double EcalBarrelDigisValidation::endcapADCtoGeV_

private

Definition at line 69 of file EcalBarrelDigisValidation.h.

Referenced by checkCalibrations(), and EcalBarrelDigisValidation().

std::map<int, double, std::less<int> > EcalBarrelDigisValidation::gainConv_

private

Definition at line 66 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), checkCalibrations(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBDigiADCAnalog_[10]

private

Definition at line 77 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBDigiADCg12_[10]

private

Definition at line 82 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBDigiADCg1_[10]

private

Definition at line 80 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBDigiADCg6_[10]

private

Definition at line 81 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBDigiADCGlobal_

private

Definition at line 75 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBDigiADCgS_[10]

private

Definition at line 79 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBDigiGain_[10]

private

Definition at line 84 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBDigiMultiplicity_

private

Definition at line 73 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBDigiOccupancy_

private

Definition at line 71 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBMaximumgt100ADC_

private

Definition at line 88 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBMaximumgt10ADC_

private

Definition at line 90 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBnADCafterSwitch_

private

Definition at line 92 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

MonitorElement* EcalBarrelDigisValidation::meEBPedestal_

private

Definition at line 86 of file EcalBarrelDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalBarrelDigisValidation().

std::string EcalBarrelDigisValidation::outputFile_

private

Definition at line 62 of file EcalBarrelDigisValidation.h.

bool EcalBarrelDigisValidation::verbose_

private

Definition at line 60 of file EcalBarrelDigisValidation.h.

Referenced by EcalBarrelDigisValidation().