BTaggingCommission/runCode4ttbar.C at master · alphatsai/BTaggingCommission · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
{
    const int i=runN;
    //const int i=1;

    // Options for CommPlotProducer4ttbar::Loop
    int trig_data=3;
    float PtMin_Cut=30;  // useless now
    float PtMax_Cut=500;
    TString syst="";

    // Samples defination
    //Data13TeV_MuonEG_2016B
    //MC13TeV_DY10to50
    //MC13TeV_DY50toInf
    //MC13TeV_SingleT_tW
    //MC13TeV_SingleTbar_tW
    //MC13TeV_TTJets_powheg_pythia8
    //MC13TeV_WW
    //MC13TeV_WZ
    //MC13TeV_ZZ
    //
    const int Nsamples=10;
    float luminosity = 6260; // /pb see usage.txt
    //float luminosity = 2600; // /pb see usage.txt
    //float luminosity = 218.042; // /pb see usage.txt
    //float luminosity = 589.3;
    //float luminosity = 804.2;
    TString iname[Nsamples] = {
                                //"Data13TeV_MuonEG_2016B",       // 0
                                //"Data13TeV_MuonEG_2016B_589pb", // 0
                                //"Data13TeV_MuonEG_2016B_804pb",  // 0
                                //"Data13TeV_MuonEG_2016B_2fb",  // 0
                                "Data13TeV_MuonEG_2016B_6fb",  // 0
                                "Data13TeV_MuonEG_2016C_6fb",  // 1
                                "MC13TeV_TTJets_powheg_pythia8", // 2
                                "MC13TeV_DY10to50",              // 3
                                "MC13TeV_DY50toInf",             // 4
                                "MC13TeV_SingleT_tW",            // 5
                                "MC13TeV_SingleTbar_tW",         // 6
                                "MC13TeV_WW",                    // 7
                                "MC13TeV_WZ",                    // 8
                                "MC13TeV_ZZ"                     // 9
                                };
    TString oname[Nsamples] = { "output_dataB_mueg",
                                "output_dataC_mueg",
                                "output_ttbar",
                                "output_dy1",
                                "output_dy2",
                                "output_st1",
                                "output_st2",
                                "output_ww",
                                "output_wz",
                                "output_zz"
                                };

    //float sf_dy=1.3; // in Run2015
    float sf_dy=1;
    int datatype[Nsamples] = {0, 0, 1,     2,    2,     3,     3,      4,     5,    6 };
    float xsec[Nsamples]   = {0, 0, 831.77, 18610*sf_dy, 6025*sf_dy, 35.85, 35.85, 118.7, 47.13, 16.5 };

    //TString samplepath="root://eoscms.cern.ch//eos/cms/store/group/phys_btag/Commissioning/TTbar/Run2016/May2016_v0";
    //TString samplepath="root://eoscms.cern.ch//eos/cms/store/group/phys_btag/Commissioning/TTbar/Run2016/May2016_v1";
    TString samplepath="root://eoscms.cern.ch//eos/cms/store/group/phys_btag/Commissioning/TTbar/Run2016/July12_v0";
    TString sampleroot="MergedJetTree.root";

    // Run...
    TFile* f = TFile::Open(samplepath+"/"+iname[i]+"/"+sampleroot);
    TChain* tree = (TChain*)f->Get("btagana/ttree");
    CommPlotProducer4ttbar *run = new CommPlotProducer4ttbar(tree,1,1);

    if( xsec[i] != 0 )
    {
        TH1F* inputWeight = (TH1F*)f->Get("ttbarselectionproducer/wgtcounter");
        double wgtcounter = inputWeight->GetBinContent(1);
        run->SetNorm(xsec[i]*luminosity/wgtcounter);
        run->Loop( datatype[i], trig_data, PtMin_Cut, PtMax_Cut, oname[i], inputWeight, syst);
    }
    else
    {
        run->Loop( datatype[i], trig_data, PtMin_Cut, PtMax_Cut, oname[i], 0, syst);
    }
}