AliDJetVReader.cxx

Go to the documentation of this file.

/**************************************************************************
 * Copyright(c) 1998-2016, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

#include <iostream>
#include <cstring>
#include <set>

#include <TMath.h>
#include <TFile.h>
#include <TCanvas.h>
#include <TLegend.h>
#include <TLatex.h>
#include <TRandom2.h>
#include <TH1D.h>
#include <TH2.h>
#include <TF1.h>
#include <TDatabasePDG.h>
#include <Riostream.h>
#include <TLine.h>

#include "AliDJetVReader.h"

ClassImp(AliDJetVReader);

AliDJetVReader::AliDJetVReader():
  TObject(),
  fpTmin(0),
  fpTmax(0),
  fnDbins(0),
  fDbinpTedges(nullptr),
  fnJetPtbins(0),
  fJetPtBinEdges(nullptr),
  fnJetzbins(0),
  fJetzBinEdges(nullptr),
  fDEffValues(nullptr),
  fMassRebin(1),
  fEfficiencyWeightSB(kFALSE),
  fMassPlot(nullptr),
  fMassVsJetPtPlot(nullptr),
  fMassVsJetzPlot(nullptr)
{
}

AliDJetVReader::AliDJetVReader(const AliDJetVReader &source):
  TObject(),
  fpTmin(source.fpTmin),
  fpTmax(source.fpTmax),
  fnDbins(0),
  fDbinpTedges(nullptr),
  fnJetPtbins(0),
  fJetPtBinEdges(nullptr),
  fnJetzbins(0),
  fJetzBinEdges(nullptr),
  fDEffValues(nullptr),
  fMassRebin(source.fMassRebin),
  fEfficiencyWeightSB(source.fEfficiencyWeightSB),
  fMassPlot(nullptr),
  fMassVsJetPtPlot(nullptr),
  fMassVsJetzPlot(nullptr)
{
  if (source.fnDbins > 0) {
    fnDbins = source.fnDbins;
    fDbinpTedges = new Double_t[fnDbins+1];
    memcpy(fDbinpTedges, source.fDbinpTedges, sizeof(Double_t)*(fnDbins+1));
    fDEffValues = new Double_t[fnDbins+1];
    memcpy(fDEffValues, source.fDEffValues, sizeof(Double_t)*(fnDbins+1));
  }
  if (source.fnJetPtbins > 0) {
    fnJetPtbins = source.fnJetPtbins;
    fJetPtBinEdges = new Double_t[fnJetPtbins+1];
    memcpy(fJetPtBinEdges, source.fJetPtBinEdges, sizeof(Double_t)*(fnJetPtbins+1));
  }
  if (source.fnJetzbins > 0) {
    fnJetzbins = source.fnJetzbins;
    fJetzBinEdges = new Double_t[fnJetzbins+1];
    memcpy(fJetzBinEdges, source.fJetzBinEdges, sizeof(Double_t)*(fnJetzbins+1));
  }
}

AliDJetVReader::~AliDJetVReader()
{
  if (fMassPlot) delete fMassPlot;
  if (fMassVsJetPtPlot) delete fMassVsJetPtPlot;
  if (fMassVsJetzPlot) delete fMassVsJetzPlot;
}

void AliDJetVReader::SetDmesonPtBins(Int_t nbins, Double_t* ptedges)
{
  fnDbins = nbins;
  if (fDbinpTedges) {
    delete[] fDbinpTedges;
    fDbinpTedges = nullptr;
  }
  if (nbins == 0) return;
  fDbinpTedges = new Double_t[fnDbins + 1];
  memcpy(fDbinpTedges, ptedges, sizeof(Double_t) * (fnDbins + 1));
}

void AliDJetVReader::SetJetPtBins(Int_t nbins, Double_t* ptedges)
{
  fnJetPtbins = nbins;
  if (fJetPtBinEdges) {
    delete[] fJetPtBinEdges;
    fJetPtBinEdges = nullptr;
  }
  if (nbins == 0) return;
  fJetPtBinEdges = new Double_t[fnJetPtbins + 1];
  memcpy(fJetPtBinEdges, ptedges, sizeof(Double_t) * (fnJetPtbins + 1));
}

void AliDJetVReader::SetJetzBins(Int_t nbins, Double_t* zedges)
{
  fnJetzbins = nbins;
  if (fJetzBinEdges) {
    delete[] fJetzBinEdges;
    fJetzBinEdges = nullptr;
  }
  if (nbins == 0) return;
  fJetzBinEdges = new Double_t[fnJetzbins + 1];
  memcpy(fJetzBinEdges, zedges, sizeof(Double_t) * (fnJetzbins + 1));
}

void AliDJetVReader::SetDmesonEfficiency(Double_t* effvalues)
{
  if (fDEffValues) {
    delete[] fDEffValues;
    fDEffValues = nullptr;
  }
  if (fnDbins == 0) return;
  fDEffValues = new Double_t[fnDbins];
  memcpy(fDEffValues, effvalues, sizeof(Double_t) * fnDbins);
}