ROOT logo
void CreateNoiseRunHisto(Int_t iLDC=0){

  // macro to create histogram as if from pulser data

  AliTOFGeometry * geom = new AliTOFGeometry();
  AliTOFDecoder * decoder = new AliTOFDecoder();
  Int_t array0[2160];
  Int_t array1[2208];
  Int_t array2[2160];
  Int_t array3[2208];
  Int_t array4[2160];
  Int_t array5[2208];
  // simulating data from LDC reading DDL # iLDC*6+{0,1,2,3,4,5}
  Int_t iDDL[6];
  for (Int_t j=0;j<6;j++){
    iDDL[j]=j+iLDC*6;
  }
  decoder->GetArrayDDL(array0,iDDL[0]);
  cout << " DDL " << iDDL[0] << " done " << endl;
  decoder->GetArrayDDL(array1,iDDL[1]);
  cout << " DDL " << iDDL[1] << " done " << endl;
  decoder->GetArrayDDL(array2,iDDL[2]);
  cout << " DDL " << iDDL[2] << " done " << endl;
  decoder->GetArrayDDL(array3,iDDL[3]);
  cout << " DDL " << iDDL[3] << " done " << endl;
  decoder->GetArrayDDL(array4,iDDL[4]);
  cout << " DDL " << iDDL[4] << " done " << endl;
  decoder->GetArrayDDL(array5,iDDL[5]);
  cout << " DDL " << iDDL[5] << " done " << endl;

  // debugging printings
  for (Int_t i=0;i<2160;i++){
    cout << " array0[" << i << "] = " << array0[i] << endl;
    cout << " array2[" << i << "] = " << array2[i] << endl;
    cout << " array4[" << i << "] = " << array4[i] << endl;
  }
  for (Int_t i=0;i<2208;i++){
    cout << " array1[" << i << "] = " << array1[i] << endl;
    cout << " array3[" << i << "] = " << array3[i] << endl;
    cout << " array5[" << i << "] = " << array5[i] << endl;
  }

  static const Int_t size = AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors();
  cout << " size = " << size << endl;
  TH1F::AddDirectory(0);
  TH1F * htofNoise = new TH1F("hTOFnoise","histo with signals on TOF during noise", size,-0.5,size-0.5);
  for (Int_t ibin=1;ibin<=size;ibin++){
    htofNoise->SetBinContent(ibin,-1);
  }

  for (Int_t i=0;i<2160;i++){
    if (htofNoise->GetBinContent(array0[i]+1)<0) htofNoise->SetBinContent(array0[i]+1,0);
    if (htofNoise->GetBinContent(array2[i]+1)<0) htofNoise->SetBinContent(array2[i]+1,0);
    if (htofNoise->GetBinContent(array4[i]+1)<0) htofNoise->SetBinContent(array4[i]+1,0);
  }

  for (Int_t i=0;i<2208;i++){
    if (htofNoise->GetBinContent(array1[i]+1)<0) htofNoise->SetBinContent(array1[i]+1,0);
    if (htofNoise->GetBinContent(array3[i]+1)<0) htofNoise->SetBinContent(array3[i]+1,0);
    if (htofNoise->GetBinContent(array5[i]+1)<0) htofNoise->SetBinContent(array5[i]+1,0);
  }

  const Int_t ndead =3;
  Int_t idead[18]; // indexes of dead channels
  Float_t ifloatdead[ndead]; // indexes of dead channels
  TRandom *rnd = new TRandom(4357);
  TDatime *dt=new TDatime();
  TDatime ciccio;
  rnd->SetSeed(ciccio.Get());
  rnd->RndmArray(ndead,ifloatdead);

  for (Int_t i=0;i<ndead;i++){
    Int_t index = Int_t(ifloatdead[i]*2160);
    idead[i*3] = array0[index];
    idead[i*3+1] = array2[index];
    idead[i*3+2] = array4[index];
    cout << " channels " << idead[i*3] << ", " << idead[i*3+1] << ", " << idead[i*3+2] << ", " << " will be noisy " << endl;
  }

  for (Int_t i=0;i<ndead;i++){
    Int_t index = Int_t(ifloatdead[i]*2208);
    idead[i*3+9] = array1[index];
    idead[i*3+9+1] = array3[index];
    idead[i*3+9+2] = array5[index];
    cout << " channels " << idead[i*3+9] << ", " << idead[i*3+9+1] << ", " << idead[i*3+9+2] << ", " << " will be noisy " << endl;
  }

  for (Int_t i=0;i<ndead*6;i++){
    cout << " noisy[" << i << "] = " << idead[i] << endl;
  }

  Bool_t tobeskipped=kFALSE;
  for (Int_t i=0;i<size;i++){
    tobeskipped = kFALSE;
    for (Int_t j=0;j<ndead*6;j++){
      if (i==idead[j]) {
	tobeskipped=kTRUE;
	cout << " noisy channel " << i << " with idead = " << idead[j] << endl;
	break;
      }
    }
    if (tobeskipped && (htofNoise->GetBinContent(i+1)!=-1)) htofNoise->SetBinContent(i+1,2.2);
    else {
      //      cout << " filling channel " << i << endl;
      //      htofNoise->Fill(i,10000);
    }
  }
  TCanvas *c = new TCanvas("c","c",-2,30,500,500);
  htofNoise->Draw();
  char filename[100];
  sprintf(filename,"$ALICE_ROOT/TOF/ShuttleInput/TOFoutNoiseLDC_%02i.root",iLDC);
  TFile *fileout = new TFile(filename,"RECREATE");
  htofNoise->Write("");
  fileout->Close();
}
 CreateNoiseRunHisto.C:1
 CreateNoiseRunHisto.C:2
 CreateNoiseRunHisto.C:3
 CreateNoiseRunHisto.C:4
 CreateNoiseRunHisto.C:5
 CreateNoiseRunHisto.C:6
 CreateNoiseRunHisto.C:7
 CreateNoiseRunHisto.C:8
 CreateNoiseRunHisto.C:9
 CreateNoiseRunHisto.C:10
 CreateNoiseRunHisto.C:11
 CreateNoiseRunHisto.C:12
 CreateNoiseRunHisto.C:13
 CreateNoiseRunHisto.C:14
 CreateNoiseRunHisto.C:15
 CreateNoiseRunHisto.C:16
 CreateNoiseRunHisto.C:17
 CreateNoiseRunHisto.C:18
 CreateNoiseRunHisto.C:19
 CreateNoiseRunHisto.C:20
 CreateNoiseRunHisto.C:21
 CreateNoiseRunHisto.C:22
 CreateNoiseRunHisto.C:23
 CreateNoiseRunHisto.C:24
 CreateNoiseRunHisto.C:25
 CreateNoiseRunHisto.C:26
 CreateNoiseRunHisto.C:27
 CreateNoiseRunHisto.C:28
 CreateNoiseRunHisto.C:29
 CreateNoiseRunHisto.C:30
 CreateNoiseRunHisto.C:31
 CreateNoiseRunHisto.C:32
 CreateNoiseRunHisto.C:33
 CreateNoiseRunHisto.C:34
 CreateNoiseRunHisto.C:35
 CreateNoiseRunHisto.C:36
 CreateNoiseRunHisto.C:37
 CreateNoiseRunHisto.C:38
 CreateNoiseRunHisto.C:39
 CreateNoiseRunHisto.C:40
 CreateNoiseRunHisto.C:41
 CreateNoiseRunHisto.C:42
 CreateNoiseRunHisto.C:43
 CreateNoiseRunHisto.C:44
 CreateNoiseRunHisto.C:45
 CreateNoiseRunHisto.C:46
 CreateNoiseRunHisto.C:47
 CreateNoiseRunHisto.C:48
 CreateNoiseRunHisto.C:49
 CreateNoiseRunHisto.C:50
 CreateNoiseRunHisto.C:51
 CreateNoiseRunHisto.C:52
 CreateNoiseRunHisto.C:53
 CreateNoiseRunHisto.C:54
 CreateNoiseRunHisto.C:55
 CreateNoiseRunHisto.C:56
 CreateNoiseRunHisto.C:57
 CreateNoiseRunHisto.C:58
 CreateNoiseRunHisto.C:59
 CreateNoiseRunHisto.C:60
 CreateNoiseRunHisto.C:61
 CreateNoiseRunHisto.C:62
 CreateNoiseRunHisto.C:63
 CreateNoiseRunHisto.C:64
 CreateNoiseRunHisto.C:65
 CreateNoiseRunHisto.C:66
 CreateNoiseRunHisto.C:67
 CreateNoiseRunHisto.C:68
 CreateNoiseRunHisto.C:69
 CreateNoiseRunHisto.C:70
 CreateNoiseRunHisto.C:71
 CreateNoiseRunHisto.C:72
 CreateNoiseRunHisto.C:73
 CreateNoiseRunHisto.C:74
 CreateNoiseRunHisto.C:75
 CreateNoiseRunHisto.C:76
 CreateNoiseRunHisto.C:77
 CreateNoiseRunHisto.C:78
 CreateNoiseRunHisto.C:79
 CreateNoiseRunHisto.C:80
 CreateNoiseRunHisto.C:81
 CreateNoiseRunHisto.C:82
 CreateNoiseRunHisto.C:83
 CreateNoiseRunHisto.C:84
 CreateNoiseRunHisto.C:85
 CreateNoiseRunHisto.C:86
 CreateNoiseRunHisto.C:87
 CreateNoiseRunHisto.C:88
 CreateNoiseRunHisto.C:89
 CreateNoiseRunHisto.C:90
 CreateNoiseRunHisto.C:91
 CreateNoiseRunHisto.C:92
 CreateNoiseRunHisto.C:93
 CreateNoiseRunHisto.C:94
 CreateNoiseRunHisto.C:95
 CreateNoiseRunHisto.C:96
 CreateNoiseRunHisto.C:97
 CreateNoiseRunHisto.C:98
 CreateNoiseRunHisto.C:99
 CreateNoiseRunHisto.C:100
 CreateNoiseRunHisto.C:101
 CreateNoiseRunHisto.C:102
 CreateNoiseRunHisto.C:103
 CreateNoiseRunHisto.C:104
 CreateNoiseRunHisto.C:105
 CreateNoiseRunHisto.C:106
 CreateNoiseRunHisto.C:107
 CreateNoiseRunHisto.C:108
 CreateNoiseRunHisto.C:109
 CreateNoiseRunHisto.C:110
 CreateNoiseRunHisto.C:111
 CreateNoiseRunHisto.C:112
 CreateNoiseRunHisto.C:113
 CreateNoiseRunHisto.C:114
 CreateNoiseRunHisto.C:115
 CreateNoiseRunHisto.C:116