#include <HLTTrack.h>

Public Member Functions
void	analyze (const edm::Handle< reco::IsolatedPixelTrackCandidateCollection > &IsoPixelTrackL3, const edm::Handle< reco::IsolatedPixelTrackCandidateCollection > &IsoPixelTrackL2, const edm::Handle< reco::VertexCollection > &pixelVertices, const edm::Handle< reco::RecoChargedCandidateCollection > &PixelTracksL3, const edm::Handle< FEDRawDataCollection > hfedraw, const edm::Handle< edmNew::DetSetVector< SiPixelCluster > > &pixelClusters, TTree *tree)

	HLTTrack ()

void	setup (const edm::ParameterSet &pSet, TTree *tree)

Private Attributes
bool	_Debug

bool	_Monte

int	evtCounter

float *	isopixeltrackL2dXY

float *	isopixeltrackL2eta

float *	isopixeltrackL2pt

float *	isopixeltrackL3energy

float *	isopixeltrackL3eta

float *	isopixeltrackL3maxptpxl

float *	isopixeltrackL3phi

float *	isopixeltrackL3pt

int	nisopixeltrackL3

int	npixelclusters

int	npixeltracksL3

int	pixelfedsize

float *	pixeltracksL3eta

float *	pixeltracksL3phi

float *	pixeltracksL3pt

float *	pixeltracksL3vz

Detailed Description

Definition at line 35 of file HLTTrack.h.

Constructor & Destructor Documentation

HLTTrack::HLTTrack ( )

Definition at line 15 of file HLTTrack.cc.

References _Debug, _Monte, and evtCounter.

                    {
   evtCounter=0;
 
   //set parameter defaults 
   _Monte=false;
   _Debug=false;
 }

Member Function Documentation

void HLTTrack::analyze	(	const edm::Handle< reco::IsolatedPixelTrackCandidateCollection > &	IsoPixelTrackL3,
		const edm::Handle< reco::IsolatedPixelTrackCandidateCollection > &	IsoPixelTrackL2,
		const edm::Handle< reco::VertexCollection > &	pixelVertices,
		const edm::Handle< reco::RecoChargedCandidateCollection > &	PixelTracksL3,
		const edm::Handle< FEDRawDataCollection >	hfedraw,
		const edm::Handle< edmNew::DetSetVector< SiPixelCluster > > &	pixelClusters,
		TTree *	tree
	)

Analyze the Data

Definition at line 77 of file HLTTrack.cc.

References deltaR(), FEDRawDataCollection::FEDData(), i, edm::Ref< C, T, F >::isNonnull(), isopixeltrackL2dXY, isopixeltrackL2eta, isopixeltrackL2pt, isopixeltrackL3energy, isopixeltrackL3eta, isopixeltrackL3maxptpxl, isopixeltrackL3phi, isopixeltrackL3pt, edm::HandleBase::isValid(), j, relval_2017::k, HLT_25ns10e33_v2_cff::minDR, nisopixeltrackL3, npixelclusters, npixeltracksL3, HLT_25ns10e33_v2_cff::pixelClusters, pixelfedsize, pixeltracksL3eta, pixeltracksL3phi, pixeltracksL3pt, pixeltracksL3vz, PixelVertexes_cfi::pixelVertices, and FEDRawData::size().

Referenced by HLTAnalyzer::analyze().

                                {
 
   //isoPixel
   if (IsoPixelTrackL3.isValid()) { 
     // Ref to Candidate object to be recorded in filter object 
     edm::Ref<reco::IsolatedPixelTrackCandidateCollection> candref; 
     
     nisopixeltrackL3 = IsoPixelTrackL3->size();
     
     for (unsigned int i=0; i<IsoPixelTrackL3->size(); i++) 
       { 
         candref = edm::Ref<reco::IsolatedPixelTrackCandidateCollection>(IsoPixelTrackL3, i); 
         
         isopixeltrackL3maxptpxl[i] = candref->maxPtPxl();
         isopixeltrackL3pt[i] = candref->pt();
         isopixeltrackL3eta[i] = candref->track()->eta();
         isopixeltrackL3phi[i] = candref->phi();   
         isopixeltrackL3energy[i] = (candref->pt())*cosh(candref->track()->eta());
       } 
   }
   else {nisopixeltrackL3 = 0;} 
 
   //minBiasPixel  
   if (PixelTracksL3.isValid()) { 
     // Ref to Candidate object to be recorded in filter object 
     edm::Ref<reco::RecoChargedCandidateCollection> candref; 
     
     npixeltracksL3 = PixelTracksL3->size();
     
     for (unsigned int i=0; i<PixelTracksL3->size(); i++) 
       { 
     candref = edm::Ref<reco::RecoChargedCandidateCollection>(PixelTracksL3, i); 
     
     pixeltracksL3pt[i] = candref->pt();
     pixeltracksL3eta[i] = candref->eta();
         pixeltracksL3phi[i] = candref->phi();   
     pixeltracksL3vz[i] = candref->vz();
     
     if (IsoPixelTrackL2.isValid()) {
       double minDR=100;
       edm::Ref<reco::IsolatedPixelTrackCandidateCollection> candrefl2;
       edm::Ref<reco::IsolatedPixelTrackCandidateCollection> candrefl2matched;
       for (unsigned int j=0; j<IsoPixelTrackL2->size(); j++) 
         { 
           candrefl2 = edm::Ref<reco::IsolatedPixelTrackCandidateCollection>(IsoPixelTrackL2, j);
           double drL3L2 = deltaR(candrefl2->eta(), candrefl2->phi(),candref->eta(), candref->phi());
           if (drL3L2<minDR)
         {
           candrefl2matched=candrefl2;
           minDR=drL3L2;
         }
         }
       if (candrefl2matched.isNonnull())
         {
           isopixeltrackL2pt[i]=candrefl2matched->pt();
           isopixeltrackL2eta[i]=candrefl2matched->eta();
           if (pixelVertices.isValid()) 
         {
           double minDZ=100;
           edm::Ref<reco::VertexCollection> vertref; 
           edm::Ref<reco::VertexCollection> vertrefMatched;
           for (unsigned int k=0; k<pixelVertices->size(); k++)
             {
               vertref=edm::Ref<reco::VertexCollection>(pixelVertices, k);
               double dz=fabs(candrefl2matched->track()->dz(vertref->position()));
               if (dz<minDZ)
             {
               minDZ=dz;
               vertrefMatched=vertref;
             }
             }
           if (vertrefMatched.isNonnull()) isopixeltrackL2dXY[i] = candrefl2matched->track()->dxy(vertref->position());
           else isopixeltrackL2dXY[i]=0;
         }
         }
     }
       } 
   }
   else {npixeltracksL3 = 0;} 
 
   //Pixel FED activity
   const FEDRawDataCollection theRaw = * hfedraw;
   int sumoffeds = 0; 
 
   for (unsigned int i = 0; i <= 39; i++)
     {
       sumoffeds = sumoffeds + (theRaw.FEDData(i).size());
     }
 
   pixelfedsize = sumoffeds;
 
   npixelclusters = pixelClusters->dataSize();
 
   
 }

void HLTTrack::setup	(	const edm::ParameterSet &	pSet,
		TTree *	tree
	)

Definition at line 24 of file HLTTrack.cc.

References _Debug, _Monte, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNames(), isopixeltrackL2dXY, isopixeltrackL2eta, isopixeltrackL2pt, isopixeltrackL3energy, isopixeltrackL3eta, isopixeltrackL3maxptpxl, isopixeltrackL3phi, isopixeltrackL3pt, nisopixeltrackL3, npixelclusters, npixeltracksL3, pixelfedsize, pixeltracksL3eta, pixeltracksL3phi, pixeltracksL3pt, and pixeltracksL3vz.

Referenced by HLTAnalyzer::HLTAnalyzer().

                                                                 {
 
   edm::ParameterSet myEmParams = pSet.getParameter<edm::ParameterSet>("RunParameters") ;
   std::vector<std::string> parameterNames = myEmParams.getParameterNames() ;
   
   for ( std::vector<std::string>::iterator iParam = parameterNames.begin();
     iParam != parameterNames.end(); iParam++ ){
     if  ( (*iParam) == "Monte" ) _Monte =  myEmParams.getParameter<bool>( *iParam );
     else if ( (*iParam) == "Debug" ) _Debug =  myEmParams.getParameter<bool>( *iParam );
   }
 
    //common
   const int kMaxTrackL3 = 10000;
    //isoPixel
   isopixeltrackL3pt = new float[kMaxTrackL3];
   isopixeltrackL3eta = new float[kMaxTrackL3];
   isopixeltrackL3phi = new float[kMaxTrackL3]; 
   isopixeltrackL3maxptpxl = new float[kMaxTrackL3];
   isopixeltrackL3energy = new float[kMaxTrackL3]; 
   isopixeltrackL2pt = new float[kMaxTrackL3];
   isopixeltrackL2eta = new float[kMaxTrackL3];
   isopixeltrackL2dXY = new float[kMaxTrackL3];
 
    //minBiasPixel
   pixeltracksL3pt = new float[kMaxTrackL3];
   pixeltracksL3eta = new float[kMaxTrackL3];
   pixeltracksL3phi = new float[kMaxTrackL3]; 
   pixeltracksL3vz = new float[kMaxTrackL3];
 
   // Track-specific branches of the tree 
    //isoPixel
   HltTree->Branch("NohIsoPixelTrackL3",&nisopixeltrackL3,"NohIsoPixelTrackL3/I");
   HltTree->Branch("ohIsoPixelTrackL3Pt",isopixeltrackL3pt,"ohIsoPixelTrackL3Pt[NohIsoPixelTrackL3]/F");
   HltTree->Branch("ohIsoPixelTrackL3Eta",isopixeltrackL3eta,"ohIsoPixelTrackL3Eta[NohIsoPixelTrackL3]/F");
   HltTree->Branch("ohIsoPixelTrackL3Phi",isopixeltrackL3phi,"ohIsoPixelTrackL3Phi[NohIsoPixelTrackL3]/F"); 
   HltTree->Branch("ohIsoPixelTrackL3MaxPtPxl",isopixeltrackL3maxptpxl,"ohIsoPixelTrackL3MaxPtPxl[NohIsoPixelTrackL3]/F");
   HltTree->Branch("ohIsoPixelTrackL3Energy",isopixeltrackL3energy,"ohIsoPixelTrackL3Energy[NohIsoPixelTrackL3]/F");
   HltTree->Branch("ohIsoPixelTrackL2pt",isopixeltrackL2pt,"ohIsoPixelTrackL2pt[NohIsoPixelTrackL3]/F");
   HltTree->Branch("ohIsoPixelTrackL2eta",isopixeltrackL2eta,"ohIsoPixelTrackL2eta[NohIsoPixelTrackL3]/F");
   HltTree->Branch("ohIsoPixelTrackL2dXY",isopixeltrackL2dXY,"ohIsoPixelTrackL2dXY[NohIsoPixelTrackL3]/F");
 
    //minBiasPixel
   HltTree->Branch("NohPixelTracksL3",&npixeltracksL3,"NohPixelTracksL3/I");
   HltTree->Branch("ohPixelTracksL3Pt",pixeltracksL3pt,"ohPixelTracksL3Pt[NohPixelTracksL3]/F");
   HltTree->Branch("ohPixelTracksL3Eta",pixeltracksL3eta,"ohPixelTracksL3Eta[NohPixelTracksL3]/F");
   HltTree->Branch("ohPixelTracksL3Phi",pixeltracksL3phi,"ohPixelTracksL3Phi[NohPixelTracksL3]/F"); 
   HltTree->Branch("ohPixelTracksL3Vz",pixeltracksL3vz,"ohPixelTracksL3Vz[NohPixelTracksL3]/F");
 
   HltTree->Branch("ohPixelFEDSize",&pixelfedsize,"ohPixelFEDSize/I");
   HltTree->Branch("NohPixelClusters",&npixelclusters,"NohPixelClusters/I");
 }

Member Data Documentation

bool HLTTrack::_Debug

private

Definition at line 64 of file HLTTrack.h.

Referenced by HLTTrack(), and setup().

bool HLTTrack::_Monte

private

Definition at line 64 of file HLTTrack.h.

Referenced by HLTTrack(), and setup().

int HLTTrack::evtCounter

private

Definition at line 66 of file HLTTrack.h.

Referenced by HLTTrack().

float * HLTTrack::isopixeltrackL2dXY

private