Trigger/TrigAnalysis/TrigInDetAnalysisUser/Resplot/src/main.cxx File Reference

#include <cmath>
#include <vector>
#include "TFile.h"
#include "TStyle.h"
#include "TMath.h"
#include "Resplot.h"
#include "Random.h"

Include dependency graph for Trigger/TrigAnalysis/TrigInDetAnalysisUser/Resplot/src/main.cxx:

Go to the source code of this file.

Functions
int	main ()

Detailed Description

Author

mark sutton

Date

Fri 11 Jan 2019 07:41:27 CET

Copyright (C) 2002-2019 CERN for the benefit of the ATLAS collaboration

Definition in file Trigger/TrigAnalysis/TrigInDetAnalysisUser/Resplot/src/main.cxx.

Function Documentation

main()

int main

(

)

getCoefficientMap(label, EigenIdxList)

input value:

1. label: falvor label in std::string format, could be one of B, C, T, Light
2. EigenIdxList is user defined vector containing all eigenvector index that user interested in. output: Map of format map<string, map<string, float>> containing decomposition coefficient of the list of eigenvectors defined by EigenIdxList.

getCoefficients(label, evIdx)

input value:

1. label: falvor label in std::string format, could be one of B, C, T, Light
2. evIdx: The index of eigenvector user interested in. output value: vector of coefficient values. The order is the same as output given by getListOfOriginalNuisanceParameters()

getListOfOriginalNuisanceParameters(label)

input value:

1. label: falvor label in std::string format, could be one of B, C, T, Light output value: List of original nuisance parameter names.

getNumEigenVectors(label)

input value:

1. label: falvor label in std::string format, could be one of B, C, T, Light return value: number of eigen vectors used for chosen label. Return 0 if error occured.

Definition at line 23 of file Trigger/TrigAnalysis/TrigInDetAnalysisUser/Resplot/src/main.cxx.

           { 
 
  std::cout << "main()" << std::endl;
 
  TFile f("out.root", "recreate");
 
 
  int     Nxbins   = 10;
  double  xlimitlo = -10;
  double  xlimithi =  10;
 
  int     Nybins    = 300;
  double  ylimitlo = -30;
  double  ylimithi =  30;
 
 
  Resplot  plot1("plot1", Nxbins, xlimitlo, xlimithi,   Nybins, ylimitlo, ylimithi );
  Resplot  plot2("plot2", Nxbins, xlimitlo, xlimithi,   Nybins, ylimitlo, ylimithi );
 
  std::cout << "running" << std::endl;
 
  Random& rnd = Random::RandomBuilder();
 
 
  // loop over some "events"
  for ( int i=0 ; i<10000 ; i++ ) { 
 
    if ( (i%1000)==0 ) std::cout << "process event " << i << std::endl;
 
    double x0 = 20*rnd.uniform()-10;
    double y0 = (2+x0*x0/20)*rnd.gauss();
    
    plot1.Fill(x0, y0);
    plot2.Fill(x0, y0);
  
  }
 
  std::cout << "finalising" << std::endl;
 
  plot1.Finalise(Resplot::FitNull95);
 
  std::cout << __LINE__ << std::endl;
 
  plot2.Finalise(Resplot::FitGaussian);
  
  plot1.Write();
  plot2.Write();
 
 
  f.Write();
 
 
  std::cout << "closing f" << std::endl; 
  f.Close();
 
 
  std::cout << "operations" << std::endl; 
 
  return 0;
}