EcalPreshowerNoiseDistrib.cc

Go to the documentation of this file.

/*
 * \file EcalPreshowerNoiseDistrib.cc
 *
*/

#include <Validation/EcalDigis/interface/EcalPreshowerNoiseDistrib.h>

EcalPreshowerNoiseDistrib::EcalPreshowerNoiseDistrib(const edm::ParameterSet& ps):
  ESdigiCollectionToken_( consumes<ESDigiCollection>( ps.getParameter<edm::InputTag>( "ESdigiCollection" ) ) )
{
  
  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
    
  // histos
  meESDigiMultiplicity_=0;
  for (int ii=0; ii<3; ii++ ) { meESDigiADC_[ii] = 0; }
  
}

void EcalPreshowerNoiseDistrib::bookHistograms(DQMStore::IBooker &ibooker, edm::Run const&, edm::EventSetup const&) {
  Char_t histo[200];

  sprintf (histo, "multiplicity" ) ;
  meESDigiMultiplicity_ = ibooker.book1D(histo, histo, 1000, 0., 137728);

  for ( int ii = 0; ii < 3 ; ii++ ) {
    sprintf (histo, "esRefHistos%02d", ii) ;
    meESDigiADC_[ii] = ibooker.book1D(histo, histo, 35, 983.5, 1018.5) ;
  }

  for ( int ii = 0; ii < 3 ; ii++ ) {
    sprintf (histo, "esRefHistosCorr%02d", ii) ;
    meESDigiCorr_[ii] = ibooker.book2D(histo, histo, 35, 983.5, 1018.5, 35, 983.5, 1018.5) ;
  }

  meESDigi3D_ = ibooker.book3D("meESDigi3D_", "meESDigi3D_", 35, 983.5, 1018.5, 35, 983.5, 1018.5, 35, 983.5, 1018.5) ;
}


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

  edm::Handle<ESDigiCollection> EcalDigiES;
  
  e.getByToken( ESdigiCollectionToken_ , EcalDigiES );

  // retrun if no data
  if( !EcalDigiES.isValid() ) return;
  
  // loop over Digis
  const ESDigiCollection * preshowerDigi = EcalDigiES.product () ;
  
  std::vector<double> esADCCounts ;
  esADCCounts.reserve(ESDataFrame::MAXSAMPLES);
  
  int nDigis = 0;

  for (unsigned int digis=0; digis<EcalDigiES->size(); ++digis) {
    nDigis++;
    ESDataFrame esdf=(*preshowerDigi)[digis];
    int nrSamples=esdf.size();
    for (int sample = 0 ; sample < nrSamples; ++sample) {
      ESSample mySample = esdf[sample];
      if (meESDigiADC_[sample]) { meESDigiADC_[sample] ->Fill(mySample.adc()); }
    }

    // to study correlations
    if(meESDigiCorr_[0]){ meESDigiCorr_[0]->Fill(esdf[0].adc(),esdf[1].adc()); } 
    if(meESDigiCorr_[1]){ meESDigiCorr_[1]->Fill(esdf[0].adc(),esdf[2].adc()); } 
    if(meESDigiCorr_[2]){ meESDigiCorr_[2]->Fill(esdf[1].adc(),esdf[2].adc()); } 

    // reference histo: sample0, sample1, sample2
    if ( meESDigi3D_ ) meESDigi3D_ -> Fill(esdf[0].adc(),esdf[1].adc(),esdf[2].adc());
  }
  
  if ( meESDigiMultiplicity_ ) meESDigiMultiplicity_->Fill(nDigis);
  
}