#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, gen::k, nisopixeltrackL3, npixelclusters, npixeltracksL3, 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");
}