1 #if !defined( __CINT__) || defined(__MAKECINT__)
6 #include <TProfile2D.h>
15 #include "AliCDBManager.h"
16 #include "AliCDBStorage.h"
17 #include "AliCDBEntry.h"
20 #include "../TRD/AliTRDarrayF.h"
21 #include "../TRD/AliTRDCalibPadStatus.h"
22 #include "../TRD/Cal/AliTRDCalPadStatus.h"
23 #include "../TRD/Cal/AliTRDCalDet.h"
24 #include "../TRD/Cal/AliTRDCalPad.h"
25 #include "../TRD/Cal/AliTRDCalROC.h"
26 #include "../TRD/AliTRDcalibDB.h"
32 //void PlotNoiseBaseline(Int_t run, Int_t sm, Int_t det, const char * pathdatabase="local:///d/alice12/bailhache/TestShuttle/database/", const char * pathreferencefile="local:///d/alice12/bailhache/TestShuttle/reference")
33 //void PlotNoiseBaseline(Int_t run=34529, Int_t sm=0, Int_t det=0, const char * pathdatabase="alien://Folder=/alice/data/2008/LHC08b/OCDB/", const char * pathreferencedatabase="alien://Folder=/alice/data/2008/LHC08b/Reference/")
34 //void PlotNoiseBaseline(Int_t run=1, Int_t sm=0, Int_t det=0, const char * pathdatabase="local:///d/alice12/bailhache/AliAnalysisTask/v4-13-Head/SHUTTLE/TestShuttle/TestCDB/", const char * pathreferencedatabase="local:///d/alice12/bailhache/AliAnalysisTask/v4-13-Head/SHUTTLE/TestShuttle/TestReference/")
35 void AliTRDplotNoiseBaseline(Int_t run=34529, Int_t sm=0, Int_t det=0, const char * pathdatabase="alien://Folder=/alice/data/2008/LHC08b/OCDB/", const char * pathreferencedatabase="alien://Folder=/alice/data/2008/LHC08b/Reference/")
38 //TGrid::Connect("alien://",0,0,"t");
40 AliCDBManager *CDB = AliCDBManager::Instance();
41 CDB->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
42 CDB->SetSpecificStorage("TRD/Calib/PadNoise",pathdatabase);
43 CDB->SetSpecificStorage("TRD/Calib/DetNoise",pathdatabase);
44 CDB->SetSpecificStorage("TRD/Calib/PadStatus",pathdatabase);
47 AliTRDcalibDB *cal = AliTRDcalibDB::Instance();
49 //const AliTRDCalDet *u = cal->GetNoiseDet();
51 AliTRDCalDet *u = new AliTRDCalDet("u","u");
52 for(Int_t k = 0; k < 540; k++){
57 //************************
58 gStyle->SetPalette(1);
59 gStyle->SetOptStat(1111);
60 gStyle->SetPadBorderMode(0);
61 gStyle->SetCanvasColor(10);
62 gStyle->SetPadLeftMargin(0.13);
63 gStyle->SetPadRightMargin(0.13);
66 //********************************
68 AliTRDCalPad *ki = new AliTRDCalPad("testnoise","testnoise");
69 for(Int_t k = 0; k < 540; k++){
70 ki->SetCalROC(k,(AliTRDCalROC *) cal->GetNoiseROC(k));
74 Int_t smn = (Int_t) det/30;
75 if((smn==0) || (smn==1) || (smn==2) || (smn==9) || (smn==10) || (smn==11)) smn = 1;
76 if((smn==3) || (smn==4) || (smn==5) || (smn==12) || (smn==13) || (smn==14)) smn = 2;
77 if((smn==6) || (smn==7) || (smn==8) || (smn==15) || (smn==16) || (smn==17)) smn = 3;
78 TString name("TRD/DAQData/PadStatus");
81 AliCDBEntry *entrypadstatus = AliCDBManager::Instance()->Get("TRD/Calib/PadStatus",run);
82 if(!entrypadstatus) return;
83 AliTRDCalPadStatus *lo = (AliTRDCalPadStatus *)entrypadstatus->GetObject();
84 AliCDBEntry *entryo = AliCDBManager::Instance()->GetStorage(pathreferencedatabase)->Get(name, run);
86 AliTRDCalibPadStatus *calpad = (AliTRDCalibPadStatus *) entryo->GetObject();
95 TCanvas *cnoise = new TCanvas((const char*)"noise1",(const char*)"noise1",50,50,600,800);
98 ((TH2F *)ki->MakeHisto2DSmPl(sm,0,u,0,0.0,3.5,-1))->Draw("colz");
100 ((TH2F *)ki->MakeHisto2DSmPl(sm,1,u,0,0.0,3.5,-1))->Draw("colz");
102 ((TH2F *)ki->MakeHisto2DSmPl(sm,2,u,0,0.0,3.5,-1))->Draw("colz");
104 ((TH2F *)ki->MakeHisto2DSmPl(sm,3,u,0,0.0,3.5,-1))->Draw("colz");
106 ((TH2F *)ki->MakeHisto2DSmPl(sm,4,u,0,0.0,3.5,-1))->Draw("colz");
108 ((TH2F *)ki->MakeHisto2DSmPl(sm,5,u,0,0.0,3.5,-1))->Draw("colz");
112 TCanvas *cpadstatus = new TCanvas((const char*)"padstatus",(const char*)"padstatus",50,50,600,800);
113 cpadstatus->Divide(3,2);
115 ((TH2F *)lo->MakeHisto2DSmPl(sm,0))->Draw("colz");
117 ((TH2F *)lo->MakeHisto2DSmPl(sm,1))->Draw("colz");
119 ((TH2F *)lo->MakeHisto2DSmPl(sm,2))->Draw("colz");
121 ((TH2F *)lo->MakeHisto2DSmPl(sm,3))->Draw("colz");
123 ((TH2F *)lo->MakeHisto2DSmPl(sm,4))->Draw("colz");
125 ((TH2F *)lo->MakeHisto2DSmPl(sm,5))->Draw("colz");
131 TCanvas *cpoui = new TCanvas((const char*)"cpoui",(const char*)"cpoui",50,50,600,800);
133 ((TH2F *)calpad->GetHisto(det))->Draw("lego");
136 AliTRDCalROC *ouip = calpad->GetCalRocMean(det);
137 TCanvas *cpouilo = new TCanvas((const char*)"cpouilo",(const char*)"cpouilo",50,50,600,800);
138 cpouilo->Divide(2,1);
140 ((TH1F *)ouip->MakeHisto1D(8.5,10.5,-1,10.0))->Draw();
141 //((TH1F *)ouip->MakeHisto1D(0.85,1.05,-1))->Draw();
143 ((TH2F *)ouip->MakeHisto2D(8.5,10.5,-1,10.0))->Draw("colz");
144 //((TH2F *)ouip->MakeHisto2D(0.85,1.05,-1))->Draw("colz");
146 AliTRDCalROC *ouiphy = calpad->GetCalRocRMS(det);
147 TCanvas *cpouiloh = new TCanvas((const char*)"cpouiloh",(const char*)"cpouiloh",50,50,600,800);
148 cpouiloh->Divide(2,1);
150 ((TH1F *)ouiphy->MakeHisto1D(0.1,4.5,-1,10.0))->Draw();
151 //((TH1F *)ouiphy->MakeHisto1D(0.01,0.45,-1))->Draw();
153 ((TH2F *)ouiphy->MakeHisto2D(0.1,4.5,-1,10.0))->Draw("colz");
154 //((TH2F *)ouiphy->MakeHisto2D(0.01,0.45,-1))->Draw("colz");