AnalysisScripts/get_deadchannelMap_float.C at master · tklemenz/AnalysisScripts · 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
83
84
85
86
/// pedestalFile0 from first FEC setting
/// pedestalFile1 from second FEC setting

int get_deadchannelMap(TString pedestalFile0, TString pedestalFile1, TString filename="DeadChannelMap.root")
{
  using namespace o2::TPC;
  TFile f0(pedestalFile0);
  TFile f1(pedestalFile1);
  gROOT->cd();

  static const Mapper& mapper = Mapper::instance();
  ROC roc;

  CalDet<float> *pedestal0 = nullptr, *pedestal1 = nullptr;
  f0.GetObject("Pedestals", pedestal0);
  f1.GetObject("Pedestals", pedestal1);

  CalDet<float> *DCM = new CalDet<float>(PadSubset::ROC);
  //DCM->setName("DeadChannelMap");
  CalArray<float>& DeadChannelMap = DCM->getCalArray(roc);

  CalDet<float> *dummy = new CalDet<float>(PadSubset::ROC);
  CalArray<float>& dummyArr = dummy->getCalArray(roc);

  for (int irow=0; irow<63; ++irow) {
    int npads = mapper.getNumberOfPadsInRowROC(roc, irow);
    //cout<<irow<<"\t"<<npads<<endl;
    for (int ipad=0; ipad<npads; ++ipad) {
      //cout<<endl<<pad<<endl;
      const auto ped0 = pedestal0->getValue(roc, irow, ipad);
      const auto ped1 = pedestal1->getValue(roc, irow, ipad);
      if (ped0 != 0 || ped1 != 0) {
        dummyArr.setValue(irow, ipad, 1);
        DeadChannelMap.setValue(irow, ipad, 1);

        std::cout << DeadChannelMap.getValue(irow, ipad) << "true \n";
        //cout<<irow<<"\t"<<ipad<<endl;
      }
      else {
        dummyArr.setValue(irow, ipad, 0);
        DeadChannelMap.setValue(irow, ipad, 0);
        std::cout << DeadChannelMap.getValue(irow, ipad) << "false \n";
      }
    }
  }

  dummyArr.setValue(32, 51, 0);
  dummyArr.setValue(28, 60, 0);
  dummyArr.setValue(47, 57, 0);
  dummyArr.setValue(48, 71, 0);
  dummyArr.setValue(58, 63, 0);
  dummyArr.setValue(59, 61, 0);
  dummyArr.setValue(61, 69, 0);
  dummyArr.setValue(44, 57, 0);
  dummyArr.setValue(4, 42, 0);
  dummyArr.setValue(4, 33, 0);
  dummyArr.setValue(32, 33, 0);
  dummyArr.setValue(58, 35, 0);
  dummyArr.setValue(59, 34, 0);

  DeadChannelMap.setValue(32, 51, 0);
  DeadChannelMap.setValue(28, 60, 0);
  DeadChannelMap.setValue(47, 57, 0);
  DeadChannelMap.setValue(48, 71, 0);
  DeadChannelMap.setValue(58, 63, 0);
  DeadChannelMap.setValue(59, 61, 0);
  DeadChannelMap.setValue(61, 69, 0);
  DeadChannelMap.setValue(44, 57, 0);
  DeadChannelMap.setValue(4, 42, 0);
  DeadChannelMap.setValue(4, 33, 0);
  DeadChannelMap.setValue(32, 33, 0);
  DeadChannelMap.setValue(58, 35, 0);
  DeadChannelMap.setValue(59, 34, 0);

  auto dummyplot = Painter::getHistogram2D(dummyArr);

  auto cDummyDCM = new TCanvas("cDummyDCM", "DummyDCM");
  dummyplot->Draw("colz");

  TFile *f = TFile::Open(filename, "recreate");
  f->WriteObject(DCM, "DeadChannelMap");
  f->Close();
  delete f;

  return 0;
}