atlasresearch/temp.c at master · donjpierce/atlasresearch · 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
{

	Int_t k = 0;
	for (int q = 0; q < 4; q++)
	{
		for (int l = 0; l < 4; l++)
		{
			if (q == l)	continue;
			correlationgraph[k]->SetTitle("2016 Muons (L1KFmuontriggers...48Runs-002) for L1 > 50GeV");
			k++;
		}
	}


}

	/*

	TFile *muonFile = TFile::Open("../PhysicsMain.L1KFmuontriggers.2016.f731f758_m1659m1710.Run309759.48Runs-002.root");
	TTree* muonTree = (TTree*)muonFile->Get("tree");
	Int_t passmuonmed, passmuonvarmed, muonrecal, muonclean;
	Float_t ml1, muonint;
	muonTree->SetBranchAddress("passmu26med", &passmuonmed);
	muonTree->SetBranchAddress("passmu26varmed", &passmuonvarmed);
	muonTree->SetBranchAddress("metl1", &ml1);
	muonTree->SetBranchAddress("actint", &muonint);
	muonTree->SetBranchAddress("recalbroke", &muonrecal);
	muonTree->SetBranchAddress("passcleancuts", &muonclean);

	TFile *zbFile = TFile::Open("../PhysicsMain.L1KFnoalgXEtriggers.2016.f731f758_m1659m1710.Run309759.48Runs-001.root");
	TTree *zbTree = (TTree*)zbFile->Get("tree");
	Int_t zbl1gt10, zbl1gt30, zbl1gt40, zbl1gt45;
	Float_t zbl1, zbint;
	zbTree->SetBranchAddress("passnoalgL1XE10", &zbl1gt10);
	zbTree->SetBranchAddress("passnoalgL1XE30", &zbl1gt30);
	zbTree->SetBranchAddress("passnoalgL1XE40", &zbl1gt40);
	zbTree->SetBranchAddress("passnoalgL1XE45", &zbl1gt45);
	zbTree->SetBranchAddress("actint", &zbint);

	//___Calculate Tail Events Based on Resolutions___
	TString metalgName[4] = {"metl1", "metl1kf", "metl1mht", "metl1mhtkf"};
	TString setalgName[4] = {"setl1", "setl1kf", "setl1mht", "setl1mhtkf"};

	// create arrays for MET and SET branches
	Float_t met[4]; Float_t set[4];
	for (Int_t i = 0; i < 4; i++)
	{
		zbTree->SetBranchAddress(metalgName[i], &met[i]);
		zbTree->SetBranchAddress(setalgName[i], &set[i]);
	}

		TH2F *correlationgraph = new TH2F("correlationgraph", "",100, 0., 100., 100, 0., 100.);

		TTree* runtree = zbTree;

		Int_t n = 0;
		Long64_t nentries = runtree->GetEntries();
		for (Int_t i = 0; i < nentries/10; i++)
		{
			runtree->GetEntry(i);
			//std::cout << zbl1gt10 << zbl1gt30 << zbl1gt40 << zbl1gt45  << std::endl;

			Double_t x[4]; // x = bulkmet and y = tailmet will be calculated for each algorithm
			Double_t y[4];

			if ( (zbl1gt10 > 0.1 || zbl1gt30 > 0.1 || zbl1gt40 > 0.1 || zbl1gt45 > 0.1) && met[0] > 50.)
			{
				x[0] = met[0];
				y[0] = met[1];
				correlationgraph->Fill(x[0], y[0]); // and populate the appropraite correlationgraph
			}
		}

}