TPGCheck.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    TPGCheck
// Class:      TPGCheck
// 
//
// Original Author:  Muriel Cerutti
//         Created:  Thu Oct 26 10:47:17 CEST 2006
//


// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"

#include <vector>
#include <string>

#include "TH1I.h"
#include "TFile.h"

using namespace edm;
using namespace std;

//
// class declaration
//

class TPGCheck : public edm::EDAnalyzer {
   public:
      explicit TPGCheck(const edm::ParameterSet&);
      ~TPGCheck();


   private:
      virtual void beginJob() override ;
      virtual void analyze(const edm::Event&, const edm::EventSetup&) override;
      virtual void endJob() override ;

      // ----------member data ---------------------------
      TH1I *ecal_et_[2];
      TH1I *ecal_tt_[2];
      TH1I * ecal_fgvb_[2];
      
      TFile *histFile_;
      std::string label_;
      std::string producer_;
      std::vector<std::string> ecal_parts_;
  // fix for consumes
      edm::EDGetTokenT<EcalTrigPrimDigiCollection> ecal_tp_token_;
};

//
// constructors and destructor
//
TPGCheck::TPGCheck(const edm::ParameterSet& iConfig)
{
   //now do what ever initialization is needed
   
  ecal_parts_.push_back("Barrel");
  ecal_parts_.push_back("Endcap");

  histFile_=new TFile("histos.root","RECREATE");
  for (unsigned int i=0;i<2;++i) {
    // Energy
    char t[30];
    sprintf(t,"%s_energy",ecal_parts_[i].c_str());  
    ecal_et_[i]=new TH1I(t,"Et",255,0,255);
    
    // Trigger Tower flag
    char titleTTF[30];
    sprintf(titleTTF,"%s_ttf",ecal_parts_[i].c_str());
    ecal_tt_[i]=new TH1I(titleTTF,"TTF",10,0,10);
    
    // Fain Grain
    char titleFG[30];
    sprintf(titleFG,"%s_fgvb",ecal_parts_[i].c_str());
    ecal_fgvb_[i]=new TH1I(titleFG,"FGVB",10,0,10);
  }
  
  label_= iConfig.getParameter<std::string>("Label");
  producer_= iConfig.getParameter<std::string>("Producer");
  ecal_tp_token_ = consumes<EcalTrigPrimDigiCollection> (edm::InputTag(label_,producer_));
}


TPGCheck::~TPGCheck()
{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   
   histFile_->Write();
   histFile_->Close();

}


//
// member functions
//

// ------------ method called to for each event  ------------
void
TPGCheck::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  
  // Get input
  edm::Handle<EcalTrigPrimDigiCollection> tp;
  iEvent.getByToken(ecal_tp_token_,tp);
  for (unsigned int i=0;i<tp.product()->size();i++) {  
    EcalTriggerPrimitiveDigi d=(*(tp.product()))[i]; 
    int subdet=d.id().subDet()-1;
    // for endcap, regroup double TP-s that are generated for the 2 interior rings
    if (subdet==0) {
      ecal_et_[subdet]->Fill(d.compressedEt());
    }
    else {
      if (d.id().ietaAbs()==27 || d.id().ietaAbs()==28) {
    if (i%2) ecal_et_[subdet]->Fill(d.compressedEt()*2.);
      }
      else ecal_et_[subdet]->Fill(d.compressedEt());
    }
    ecal_tt_[subdet]->Fill(d.ttFlag());
    ecal_fgvb_[subdet]->Fill(d.fineGrain());
  }
}


// ------------ method called once each job just before starting event loop  ------------
void 
TPGCheck::beginJob()
{
}

// ------------ method called once each job just after ending the event loop  ------------
void 
TPGCheck::endJob() {
   for (unsigned int i=0;i<2;++i) {
    ecal_et_[i]->Write();
    ecal_tt_[i]->Write();
    ecal_fgvb_[i]->Write();
  }
}

//define this as a plug-in
DEFINE_FWK_MODULE(TPGCheck);