EcalPreshowerDigisValidation.cc

Go to the documentation of this file.

/*
 * \file EcalPreshowerDigisValidation.cc
 *
 * \author F. Cossutti
 *
*/

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

EcalPreshowerDigisValidation::EcalPreshowerDigisValidation(const edm::ParameterSet& ps):
  ESdigiCollectionToken_( consumes<ESDigiCollection>( ps.getParameter<edm::InputTag>( "ESdigiCollection" ) ) )
{
  
  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
                                                                                                                                           
  meESDigiMultiplicity_=nullptr;

  for (int i = 0; i < 3 ; i++ ) {
    meESDigiADC_[i] = nullptr;
  }
 
}

void EcalPreshowerDigisValidation::bookHistograms(DQMStore::IBooker &ibooker, edm::Run const&, edm::EventSetup const&){

  Char_t histo[200];
 
  ibooker.setCurrentFolder("EcalDigisV/EcalDigiTask");

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

  for ( int i = 0; i < 3 ; i++ ) {

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

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

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

  edm::Handle<ESDigiCollection> EcalDigiES;

  e.getByToken( ESdigiCollectionToken_ , EcalDigiES );

  // Return if no preshower data
  if( !EcalDigiES.isValid() ) return;

  // PRESHOWER
  
  // 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) {

    ESDataFrame esdf=(*preshowerDigi)[digis];
    int nrSamples=esdf.size();
    
    ESDetId esid = esdf.id () ;
    
    nDigis++;
    
    for (int sample = 0 ; sample < nrSamples; ++sample) {
      esADCCounts[sample] = 0.;
    }
    
    for (int sample = 0 ; sample < nrSamples; ++sample) {
      ESSample mySample = esdf[sample];
      esADCCounts[sample] = (mySample.adc()) ;
    }
    if (verbose_) {
      LogDebug("DigiInfo") << "Preshower Digi for ESDetId: z side " << esid.zside() << "  plane " << esid.plane() << esid.six() << ',' << esid.siy() << ':' << esid.strip();
      for ( int i = 0; i < 3 ; i++ ) {
    LogDebug("DigiInfo") << "sample " << i << " ADC = " << esADCCounts[i];
      }
    }
    
    for ( int i = 0 ; i < 3 ; i++ ) {
      if (meESDigiADC_[i]) meESDigiADC_[i]->Fill( esADCCounts[i] ) ;
    }
    
  } 
  
  if ( meESDigiMultiplicity_ ) meESDigiMultiplicity_->Fill(nDigis);
  
}