]>
Commit | Line | Data |
---|---|---|
c9ae910e | 1 | Float_t gpTMin_T1 = 6.0; |
2 | Float_t gpTMax_T1 = 14.0; | |
3 | Float_t gpTMin_T2 = 4.0; | |
4 | Float_t gpTMax_T2 = 10.0; | |
5 | Float_t gpTMin_assoc = 1.0; | |
6 | Float_t gpTMax_assoc = 6.0; | |
7 | Float_t gZVtxRange = -1; | |
8 | Float_t gAxis = 4; | |
9 | ||
10 | void loadlibs() | |
11 | { | |
12 | gSystem->Load("libANALYSIS"); | |
13 | gSystem->Load("libANALYSISalice"); | |
14 | gSystem->Load("libCORRFW"); | |
15 | gSystem->Load("libPWGTools"); | |
16 | gSystem->Load("libPWGCFCorrelationsBase"); | |
17 | } | |
18 | ||
19 | void* events = 0; | |
20 | ||
21 | void* GetTwoPlusOne(const char* fileName, TList** listRef = 0, Bool_t mixed = kFALSE, const char* tag = "") | |
22 | { | |
23 | file = TFile::Open(fileName); | |
24 | if (!file) | |
25 | return 0; | |
26 | ||
27 | list = (TList*) gFile->Get("PWGCF_TwoPlusOne/histosTwoPlusOne"); | |
28 | ||
29 | if (!list) | |
30 | return 0; | |
31 | ||
32 | if (listRef) | |
33 | *listRef = list; | |
34 | ||
35 | return events = list->FindObject("AliTwoPlusOneContainer"); | |
36 | } | |
37 | ||
38 | void PlotQA(const char* fileName, const char* tag = "") | |
39 | { | |
40 | loadlibs(); | |
41 | ||
42 | TFile::Open(fileName); | |
43 | ||
44 | AliTwoPlusOneContainer* h = (AliTwoPlusOneContainer*) GetTwoPlusOne(fileName, 0, kFALSE, tag); | |
45 | ||
46 | if (h->GetData()->GetTrackHist(0)->GetGrid(6)->GetGrid()->GetNbins() == 0) | |
47 | { | |
48 | Printf("We have %d axes", ((AliTHn*) h->GetData()->GetTrackHist(0)->GetNVar())); | |
49 | ||
50 | ((AliTHn*) h->GetData()->GetTrackHist(0))->FillParent(); | |
51 | ((AliTHn*) h->GetData()->GetTrackHist(0))->DeleteContainers(); | |
52 | } | |
53 | ||
54 | TCanvas* c1 = new TCanvas("can1", "can1", 1200, 800); | |
55 | c1->Divide(2, 1); | |
56 | ||
57 | c1->cd(1); | |
58 | AliCFGridSparse* near_plot = h->GetData()->GetTrackHist(0)->GetGrid((AliUEHist::CFStep) AliTwoPlusOneContainer::kSameNS); | |
59 | //near_plot->SetRangeUser(2, gpTMin_T1, gpTMax_T1); | |
60 | //near_plot->SetRangeUser(6, gpTMin_T2, gpTMax_T2); | |
61 | //near_plot->SetRangeUser(1, gpTMin_assoc, gpTMax_assoc); | |
62 | TH1D* tracks_near = near_plot->Project(gAxis); | |
63 | ||
64 | tracks_near->DrawCopy(); | |
65 | c1->cd(2); | |
66 | ||
67 | Printf("\n output_TwoPlusOne: entries %f \n", h->GetData()->GetTrackHist(AliUEHist::kToward)->GetGrid((AliUEHist::CFStep) AliTwoPlusOneContainer::kSameNS)->Project(1)->GetEntries()); | |
68 | ||
69 | AliCFGridSparse* away_plot = h->GetData()->GetTrackHist(0)->GetGrid((AliUEHist::CFStep) AliTwoPlusOneContainer::kMixedNS); | |
70 | //away_plot->SetRangeUser(2, gpTMin_T1, gpTMax_T1); | |
71 | //away_plot->SetRangeUser(6, gpTMin_T2, gpTMax_T2); | |
72 | //away_plot->SetRangeUser(1, gpTMin_assoc, gpTMax_assoc); | |
73 | TH1D* tracks_away = away_plot->Project(gAxis); | |
74 | tracks_away->DrawCopy(); | |
75 | ||
76 | AliCFGridSparse* trigger = (AliCFGridSparse*) h->GetData()->GetEventHist()->GetGrid((AliUEHist::CFStep) AliTwoPlusOneContainer::kSameNS); | |
77 | //trigger->SetRangeUser(0, gpTMin_T1, gpTMax_T1); | |
78 | ||
79 | TCanvas* c2 = new TCanvas("can2", "can2", 900, 600); | |
80 | TH1D* trigger_hist = trigger->Project(0); | |
81 | trigger_hist->DrawCopy(); | |
82 | ||
83 | ||
84 | TCanvas* c3 = new TCanvas("can3", "can3", 900, 600); | |
85 | TH2* h2_near = h->GetData()->GetSumOfRatios2(h->GetData(), (AliUEHist::CFStep) AliTwoPlusOneContainer::kSameNS, 0, 4.0, 14.0, 0, 100, kFALSE, (AliUEHist::CFStep) AliTwoPlusOneContainer::kMixedNS); | |
86 | ||
87 | h2_near->DrawCopy("surf1"); | |
88 | } | |
89 | ||
90 | ||
91 | ||
92 |