InputAnalyzer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    InputAnalyzer
// Class:      InputAnalyzer
// 
//
// Original Author:  Emilia Lubenova Becheva
//         Created:  Mon Apr 20 13:43:06 CEST 2009
// $Id: InputAnalyzer.cc,v 1.3 2011/11/15 21:57:47 gowdy Exp $
//
//


// system include files
#include <memory>

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

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"


#include "DataFormats/Provenance/interface/ProductID.h"
#include "DataFormats/Common/interface/Handle.h"

#include "FWCore/Utilities/interface/InputTag.h"

#include "InputAnalyzer.h"

#include "SimDataFormats/Track/interface/SimTrackContainer.h"
#include "SimDataFormats/CrossingFrame/interface/PCrossingFrame.h"

//
// constructors and destructor
//
namespace edm
{

InputAnalyzer::InputAnalyzer(const edm::ParameterSet& iConfig)
{

 dataStep2_ = iConfig.getParameter<bool>("dataStep2");
 
 if (dataStep2_)
   // The data file contain the PCrossingFrame<SimTrack> 
   label_   = iConfig.getParameter<edm::InputTag>("collPCF");
 else 
   // The data file contain the SimTrack
   label_   = iConfig.getParameter<edm::InputTag>("collSimTrack");
}


InputAnalyzer::~InputAnalyzer()
{
}


//
// member functions
//

// ------------ method called to for each event  ------------
void
InputAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  std::cout << " dataStep2_ = " << dataStep2_ << std::endl;
  
  if (!dataStep2_){
  // Get the SimTrack collection
  
    //double simPt=0;
   int i=0;
   
   // Get the SimTrack collection from the event
   edm::Handle<SimTrackContainer> simTracks;
   bool gotTracks = iEvent.getByLabel(label_,simTracks);
   
   if (!gotTracks)
   {  
      std::cout<<"-> Could not read SimTracks !!!!"<<std::endl;
   }
   else{
      std::cout<<"-> Could read SimTracks !!!!"<<std::endl;
      
   }

  
   // Loop over the tracks
   SimTrackContainer::const_iterator simTrack;
   for (simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack){
    i++;
    
    //simPt=(*simTrack).momentum().Pt();
    //std::cout << " # i = " << i << " simPt = " << simPt << std::endl;
   
   }
   
   
 }
 else{
 // Get the PCrossingFrame collection given as signal
   
    edm::Handle<PCrossingFrame<SimTrack> > cf_simtrack;
    bool gotTracks = iEvent.getByLabel("CFWriter","g4SimHits",cf_simtrack);
    
    if (!gotTracks)
    {  
      std::cout<<"-> Could not read PCrossingFrame<SimTracks> !!!!"<<std::endl;
    }
    else
      std::cout<<"-> Could read PCrossingFrame<SimTracks> !!!!"<<std::endl;
      
   }
}


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

// ------------ method called once each job just after ending the event loop  ------------
void 
InputAnalyzer::endJob() {
}

}//edm