Utilities.cc File Reference

#include "L1Trigger/L1TMuonEndCap/interface/Utilities.h"
#include "TRandom3.h"
#include "TStopwatch.h"
#include "TTree.h"
#include "TNtuple.h"
#include "TFile.h"
#include "TChain.h"
#include "TMath.h"

Go to the source code of this file.

Functions
void	mergeNtuples (const char *ntuplename, const char *filestomerge, const char *outputfile)
float	processPrediction (float BDTPt, int Quality, float PrelimFit)
void	sortNtupleByEvent (const char *ntuplename, const char *filenametosort, const char *outputfile)

Variables
const Double_t	ptscale [31]
const std::vector< Double_t >	ptScale = std::vector<Double_t>(ptscale, ptscale + sizeof ptscale / sizeof ptscale[0])
const Double_t	twoJets_scale [16]
const std::vector< Double_t >	twoJetsScale = std::vector<Double_t>(twoJets_scale, twoJets_scale + sizeof twoJets_scale / sizeof twoJets_scale[0])

Function Documentation

void mergeNtuples	(	const char *	ntuplename,
		const char *	filestomerge,
		const char *	outputfile
	)

Definition at line 85 of file Utilities.cc.

{
  TChain chain(ntuplename);
  chain.Add(filestomerge);
  chain.Merge(outputfile);
}

float processPrediction	(	float	BDTPt,
		int	Quality,
		float	PrelimFit
	)

Definition at line 47 of file Utilities.cc.

References RPCpg::pts, and ptscale.

{
// Discretize and scale the BDTPt prediction 

  // Fix terrible predictions
  if(BDTPt < 0) BDTPt = PrelimFit;
  if(BDTPt > 250) BDTPt = PrelimFit;
  
  float BDTPt1 = BDTPt;
  float scaleF = 1.0;
  
  // Scale based upon quality
  if (Quality == 3) scaleF = 1.15;
  if (Quality == 2) scaleF = 1.3;
  if (Quality == 1) scaleF = 1.7;
  
  BDTPt1 = scaleF*BDTPt1;
  
  // Discretize based upon ptscale
  for (int pts=0; pts<31; pts++)
    {
      if (ptscale[pts]<=BDTPt1 && ptscale[pts+1]>BDTPt1)
    {
      BDTPt1 = ptscale[pts];
      break;
    }
    }
  
  if (BDTPt1 > 140) BDTPt1 = 140;
  if (BDTPt1 < 0) BDTPt1 = 0;
  
  return BDTPt1;
}

void sortNtupleByEvent	(	const char *	ntuplename,
		const char *	filenametosort,
		const char *	outputfile
	)

Definition at line 96 of file Utilities.cc.

References f, validate-o2o-wbm::f2, i, cmsHarvester::index, and MainPageGenerator::tree.

{
  //TFile f("../../all_test_redux_post.root");
  TFile f(filenametosort);
  TNtuple *tree = (TNtuple*)f.Get(ntuplename);
  Int_t nentries = (Int_t)tree->GetEntries();
  //Drawing variable pz with no graphics option.
  //variable pz stored in array fV1 (see TTree::Draw)
  tree->Draw("Event","","goff");
  Int_t *index = new Int_t[nentries];
  //sort array containing pz in decreasing order
  //The array index contains the entry numbers in decreasing order
  TMath::Sort(nentries,tree->GetV1(),index);
  
  //open new file to store the sorted Tree
  //TFile f2("../../test_events_sorted.root","recreate");
  TFile f2(outputfile,"recreate");
  
  //Create an empty clone of the original tree
  TTree *tsorted = (TTree*)tree->CloneTree(0);
  for (Int_t i=0;i<nentries;i++) {
    tree->GetEntry(index[i]);
    tsorted->Fill();
  }
  tsorted->Write();
  delete [] index;
}

Variable Documentation

const Double_t ptscale[31]

Initial value:

=  
  { 0,     1.5,   2.0,   2.5,   3.0,   3.5,   4.0,
    4.5,   5.0,   6.0,   7.0,   8.0,  10.0,  12.0,  14.0,
    16.0,  18.0,  20.0,  25.0,  30.0,  35.0,  40.0,  45.0,
    50.0,  60.0,  70.0,  80.0,  90.0, 100.0, 120.0, 140.0 }

Definition at line 26 of file Utilities.cc.

Referenced by L1TGMT::bookHistograms(), and processPrediction().

const std::vector<Double_t> ptScale = std::vector<Double_t>(ptscale, ptscale + sizeof ptscale / sizeof ptscale[0])

Definition at line 32 of file Utilities.cc.

Referenced by CSCTFPtLUT::CSCTFPtLUT(), and CSCTFTrackProducer::produce().

const Double_t twoJets_scale[16]

Initial value:

=  
  { 0,     0.5,   1.0,    2.0,    3.0,   
    4.0,   5.0,   10.0,   20.0,   50.0,
    100,   500,   1000,   5000,   7500,  50000}

Definition at line 35 of file Utilities.cc.

const std::vector<Double_t> twoJetsScale = std::vector<Double_t>(twoJets_scale, twoJets_scale + sizeof twoJets_scale / sizeof twoJets_scale[0])

Definition at line 40 of file Utilities.cc.