BmnTriggerEfficiencyRun7.h

Go to the documentation of this file.

#include <BmnDataTriggerInfo.h>
#include <BmnLambdaAnalRun7.h>
#include <BmnTwoParticleDecay.h>
#include <TGraph.h>
 
using namespace std;
 
#ifndef BMNTRIGEFFRUN7_H
#define BMNTRIGEFFRUN7_H 1
 
class BmnTriggerEfficiencyRun7 : public BmnLambdaAnalRun7
{
  public:
    BmnTriggerEfficiencyRun7() {}
 
    BmnTriggerEfficiencyRun7(TString dir)
        : BmnLambdaAnalRun7(dir)
        , fEoSNode("root://ncm.jinr.ru/")
        , isAddTriggerCondition(kFALSE)
        , fSpectraFile("")
        , nFdCounts(1000)
        , nBdCounts(1000)
        , nMultBins(18)
    {
        // Preparing default multiplicity ranges (nTracks in prim. vertex) ...
        const Int_t minMult = 2;
        const Int_t maxMult = 20;
        Int_t multStep = (maxMult - minMult) / nMultBins;
 
        for (Int_t iBin = 0; iBin < nMultBins; iBin++) {
            Int_t step = iBin * multStep + minMult;
 
            fMultMap[iBin] = make_pair(step, step + multStep);
        }
    }
 
    virtual ~BmnTriggerEfficiencyRun7() { ; }
 
    // Setters ...
 
    void SetTrigger(TString trigg) { fTrigger = trigg; }
 
    void SetMultiplicityMap(vector<pair<Int_t, Int_t>> bins)
    {
        fMultMap.clear();
 
        nMultBins = bins.size();
 
        for (size_t iBin = 0; iBin < bins.size(); iBin++)
            fMultMap[iBin] = bins.at(iBin);
    }
 
    void SetSpectraFile(TString f) { fSpectraFile = f; }
 
    void SetEffFiles(TString f) { fEffFiles.push_back(f); }
 
    // Getters ...
 
    // Anals ...
    void triggerEfficiency();
 
  private:
    TString fEoSNode;
 
    TString fTrigger;
    Bool_t isAddTriggerCondition;
 
    TString fSpectraFile;
    vector<TString> fEffFiles;   // to be set if needed !!!
 
    enum TriggConditions
    {
        base = 0,
        severe
    };
 
    Int_t nFdCounts;
    Int_t nBdCounts;
    Int_t nMultBins;
 
    map<Int_t, pair<Int_t, Int_t>> fMultMap;
 
    // Spectrum histos ...
    TH1F*** hSpectra;   // [triggCond][multCondition]
 
    // Multiplicity histos ...
    TH1F** hMult;   // [triggCond]
 
    // Redefinition aimed at getting trigger condition when doing file list ...
    vector<TString> createFilelist(vector<TString>);
 
    // Redefinition aimed at getting trigger condition in anal ...
    void ReadFile(TString, BmnParticlePairCut* cut0);
 
    Int_t FinMultBin(Int_t multValue)
    {
        for (auto it : fMultMap) {
            Int_t bin = it.first;
            Int_t min = it.second.first;
            Int_t max = it.second.second;
 
            if (multValue >= min && multValue < max)
                return bin;
        }
        // If multBin not found ...
        return -1;
    }
 
  private:
    void DoNormalization(TH1F* h)
    {
        if (h->GetEntries() == 0)
            return;
 
        // Collecting all bin contents ...
        Double_t contentAll = 0.;
 
        for (Int_t iBin = 1; iBin < h->GetNbinsX() + 1; iBin++)
            contentAll += h->GetBinContent(iBin);
 
        // Normalizing histo ...
        for (Int_t iBin = 1; iBin < h->GetNbinsX() + 1; iBin++) {
            h->SetBinContent(iBin, h->GetBinContent(iBin) / contentAll);
            h->SetBinError(iBin, 0.);
        }
 
        h->GetYaxis()->SetRangeUser(0., 1.1 * h->GetMaximum());
    }
 
    ClassDef(BmnTriggerEfficiencyRun7, 0)
};
 
#endif