EcalBxOrbitNumberGrapher.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:   EcalBxOrbitNumberGrapher 
// Class:     EcalBxOrbitNumberGrapher 
// 
//
// Original Author:  Seth COOPER
//         Created:  Th Nov 22 5:46:22 CEST 2007
//
//

#include "CaloOnlineTools/EcalTools/plugins/EcalBxOrbitNumberGrapher.h"

using namespace cms;
using namespace edm;
using namespace std;

#include <iostream>

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
EcalBxOrbitNumberGrapher::EcalBxOrbitNumberGrapher(const edm::ParameterSet& iConfig) :
  digiProducer_(iConfig.getParameter<std::string>("RawDigis")),
  runNum_(-1),
  fileName_ (iConfig.getUntrackedParameter<std::string>("fileName", std::string("ecalURechHitHists")))
{
  
  
}


EcalBxOrbitNumberGrapher::~EcalBxOrbitNumberGrapher()
{
}


//
// member functions
//

// ------------ method called to for each event  ------------
void
EcalBxOrbitNumberGrapher::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  using namespace edm;
   using namespace cms;
  //int ievt = iEvent.id().event();
  int orbit = -100;
  int bx = -100;
  int numorbiterrors = 0;
  bool orbiterror = false;

  
  edm::Handle<EcalRawDataCollection> DCCHeaders;
  iEvent.getByLabel(digiProducer_, DCCHeaders);
  if (!DCCHeaders.isValid()) {
    edm::LogError("BxOrbitNumber") << "can't get the product for EcalRawDataCollection";
  }

  //-----------------BX STuff here
   for ( EcalRawDataCollection::const_iterator headerItr= DCCHeaders->begin();headerItr != DCCHeaders->end(); 
      ++headerItr ) {
    headerItr->getEventSettings();
    int myorbit = headerItr->getOrbit();
    int mybx = headerItr->getBX();
    
    if ( orbit == -100 )
    {
      orbit = myorbit;
    }
    else if (orbit != myorbit)
    {
       std::cout << " NOOOO This header has a conflicting orbit OTHER " << orbit << " new " << myorbit  << std::endl;
       orbiterror = true; 
       numorbiterrors++;
       orbitErrorBxDiffPlot_->Fill(myorbit - orbit);
    }
    
    if ( bx == -100 )
    {
      bx = mybx;
    }
    else if (bx != mybx)
    {
       std::cout << " NOOOO This header has a conflicting bx OTHER " << bx << " new " << mybx  << std::endl; 
    }
    //LogDebug("EcalTimingCosmic") << " Lambda " << lambda; //hmm... this isn't good, I should keep a record of the wavelength in the headers as an inactive SM might have a different wavelength for this field and make this not go through.
  }
  
  if ( (bx != -100) & ( orbit != -100) )
  {
     std::cout << " Interesting event Orbit " << orbit << " BX " << bx << std::endl;
     bxnumberPlot_->Fill(bx);
     if (orbiterror) 
     {
        orbitErrorPlot_->Fill(bx);  
     }
  }
  numberofOrbitDiffPlot_->Fill(numorbiterrors);
  
  
  if(runNum_==-1)
  {
    runNum_ = iEvent.id().run();
  }
}


// insert the hist map into the map keyed by FED number
void EcalBxOrbitNumberGrapher::initHists(int FED)
{
 
  
}

// ------------ method called once each job just before starting event loop  ------------
void 
EcalBxOrbitNumberGrapher::beginJob()
{
  bxnumberPlot_ = new TH1F("bxnumber", "BX number of interexting events",3600, 0., 3600.);
  orbitErrorPlot_ = new TH1F("bxOfOrbitDiffs", "BX number of interexting events with orbit changes",3600, 0., 3600.);
  orbitErrorBxDiffPlot_ = new TH1F("orbitErrorDiffPlot","Orbit Difference of those HEADERS that have a difference",20,-10.,10.);
  numberofOrbitDiffPlot_ = new TH1F("numberOfOrbitDiffsPlot","Number of Orbit Differences",54,0., 54.);
}

// ------------ method called once each job just after ending the event loop  ------------
void 
EcalBxOrbitNumberGrapher::endJob()
{
  using namespace std;
  fileName_ += ".bx.root";

  TFile root_file_(fileName_.c_str() , "RECREATE");

  bxnumberPlot_->Write();
  orbitErrorPlot_->Write();
  numberofOrbitDiffPlot_->Write();
  orbitErrorBxDiffPlot_->Write();
  root_file_.Close();

}