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 87 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 98 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 39 of file Utilities.cc.