Script to get # of dead channels from OCDB
[u/mrichter/AliRoot.git] / FMD / scripts / dqm / SpectraRaw.C
CommitLineData
e064ab4a 1//____________________________________________________________________
2//
3// $Id: PatternDigits.C 28055 2008-08-18 00:33:20Z cholm $
4//
5// Draw hits in the specialised FMD event display
6//
7/** Display hits
8 @ingroup FMD_script
9 */
ddbc3f82 10Bool_t
11CheckFile(const char* prefix, int number, TString& f)
12{
13 f = (Form("%s%d.csv", prefix, number));
14 std::cout << "Checking if " << f << " exists ... " << std::flush;
15 f = gSystem->Which("$(HOME)/calib/", f.Data());
16 std::cout << '"' << f << '"' << std::endl;
17 return !f.IsNull();
18}
e064ab4a 19void
20SpectraRaw(const char* file="raw.root", Int_t runno=0)
21{
e064ab4a 22 gSystem->Load("libFMDutil.so");
23
24 AliCDBManager* cdb = AliCDBManager::Instance();
25 const char* cdbUri = gSystem->Getenv("AMORE_CDB_URI");
26 cdb->SetDefaultStorage(cdbUri);
27 cdb->SetRun(runno);
28
ddbc3f82 29 AliFMDCalibStripRange* range = new AliFMDCalibStripRange;
30 AliFMDCalibSampleRate* rate = new AliFMDCalibSampleRate;
31 AliFMDCalibPedestal* peds = new AliFMDCalibPedestal;
32 AliFMDCalibGain* gains = new AliFMDCalibGain;
33 Bool_t gotConds = kFALSE;
34 Bool_t gotPeds = kFALSE;
35 Bool_t gotGains = kFALSE;
36 for (Int_t i = 1; i <= 3; i++) {
37 TString f;
38 if (CheckFile("conditions", i, f)) {
39 gotConds = kTRUE;
40 std::cout << "Reading conditions for FMD" <<i<< " from " <<f<< std::endl;
41 std::ifstream in(f.Data());
42 range->ReadFromFile(in);
43 rate->ReadFromFile(in);
44 }
45 if (CheckFile("peds", i, f)) {
46 gotPeds = kTRUE;
47 std::cout << "Reading pedestals for FMD" <<i<< " from " <<f<< std::endl;
48 std::ifstream in(f.Data());
49 peds->ReadFromFile(in);
50 }
51 if (CheckFile("gains", i, f)) {
52 gotGains = kTRUE;
53 std::cout << "Reading gains for FMD" <<i<< " from " <<f<< std::endl;
54 std::ifstream in(f.Data());
55 gains->ReadFromFile(in);
56 }
57 }
58
59 Int_t mask = (AliFMDParameters::kDeadMap|
60 AliFMDParameters::kZeroSuppression|
61 AliFMDParameters::kAltroMap);
62
63
64 if (!gotConds) mask |= AliFMDParameters::kStripRange;
65 if (!gotConds) mask |= AliFMDParameters::kSampleRate;
66 if (!gotPeds) mask |= AliFMDParameters::kPedestal;
67 if (!gotGains) mask |= AliFMDParameters::kPulseGain;
68
69 AliFMDParameters* pars = AliFMDParameters::Instance();
70 pars->Init(kFALSE, mask);
71
72 if (gotConds) pars->SetStripRange(range);
73 if (gotConds) pars->SetSampleRate(rate);
74 if (gotPeds) pars->SetPedestal(peds);
75 if (gotGains) pars->SetGain(gains);
e064ab4a 76
77 AliFMDSpectraDisplay* d = new AliFMDSpectraDisplay;
78 d->AddLoad(AliFMDInput::kRaw);
79 d->SetRawFile(file);
80 d->SetName("raw");
81 d->SetTitle("Raw");
82 d->Run();
83}
84
85//____________________________________________________________________
86//
87// EOF
88//