myRawAna.cc

Go to the documentation of this file.

// myRawAna.cc
// Description:  Access raw data, specifically event fragment size per FED
// Author: Jim Hirschauer
// Date:  28 - July - 2009
// 
#include "RecoJets/JetAnalyzers/interface/myRawAna.h"

#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "FWCore/Framework/interface/Event.h"

#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include <TROOT.h>
#include <TSystem.h>
#include <TFile.h>
#include <TCanvas.h>
#include <cmath>

using namespace edm;
using namespace reco;
using namespace std;

#define DEBUG 1
#define MAXJETS 100

// Nothing passed from .cfg yet, but I leave the syntax here for
// future reference.

myRawAna::myRawAna( const ParameterSet & cfg ) {
}


// ************************
// ************************

void myRawAna::beginJob(  ) {

  edm::Service<TFileService> fs;

  fedSize = fs->make<TH2F>( "fedSize", "fedSize", 901, -0.5, 900.5, 20000, 0., 20000. );
  totFedSize = fs->make<TH1F>( "totFedSize", "totFedSize", 200, 0., 20000. );


}

// ************************
// ************************
void myRawAna::analyze( const edm::Event& evt, const edm::EventSetup& es ) {
  
  // **************************************************************
  // ** Access FED Information
  // **************************************************************


  edm::Handle<FEDRawDataCollection> theRaw;
  bool getFed = evt.getByLabel("source", theRaw);
  
  if ( ! getFed ) {
    std::cout << "fedRawData not available" << std::endl;
  } else { // got the fed raw data
    unsigned int totalFEDsize = 0 ; 
    
    // HCAL FEDs are 700-730
    unsigned int fedStart_ = 0;
    unsigned int fedStop_ = 900;
    
    for (unsigned int i=fedStart_; i<=fedStop_; ++i) {
      fedSize->Fill(i,theRaw->FEDData(i).size());
      totalFEDsize += theRaw->FEDData(i).size() ; 
    }
    totFedSize->Fill(totalFEDsize);
  }
  
  
}

// ***********************************
// ***********************************
void myRawAna::endJob() {

}
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(myRawAna);