[u/mrichter/AliRoot.git] / TOF / CheckActiveChannels.C

CheckActiveChannels(const Char_t *);
CheckActiveChannels(Int_t nrun=-1);
CheckActiveChannelsFromCDBEntry(AliCDBEntry *);

//____________________________________________________________________

CheckActiveChannels(const Char_t *fileName)
{

  TFile *file = TFile::Open(fileName);
  AliCDBEntry *cdbe = (AliCDBEntry *)file->Get("AliCDBEntry");
  CheckActiveChannelsFromCDBEntry(cdbe);
  
}

//____________________________________________________________________

CheckActiveChannels(Int_t nrun)
{

  TGrid *alien = TGrid::Connect("alien://");
  if (!alien || !alien->IsConnected()) return;
  AliCDBManager *cdbm = AliCDBManager::Instance();
  cdbm->SetDefaultStorage("raw://");
  if (nrun==-1)
    cdbm->SetRun(80015);
  else
    cdbm->SetRun(nrun);
  AliCDBEntry *cdbe = cdbm->Get("TOF/Calib/Status");
  CheckActiveChannelsFromCDBEntry(cdbe);

}

//____________________________________________________________________

CheckActiveChannelsFromCDBEntry(AliCDBEntry *cdbe)
{

  AliTOFChannelOnlineStatusArray *array = (AliTOFChannelOnlineStatusArray *)cdbe->GetObject();
  TH1F *hStatus_hw_ok = new TH1F("hStatus_hw_ok", "HW status;index;status ok", array->GetSize(), 0., array->GetSize(););
  TH1F *hStatus_pulser_ok = new TH1F("hStatus_pulser_ok", "pulser status;index;status ok", array->GetSize(), 0., array->GetSize(););
  TH1F *hStatus_noise_ok = new TH1F("hStatus_noise_ok", "noise status;index;status ok", array->GetSize(), 0., array->GetSize(););
  TH1F *hStatus_active = new TH1F("hStatus_active", "active status;index;status ok", array->GetSize(), 0., array->GetSize(););
  Bool_t hw_ok, pulser_ok, noise_ok;
  for (Int_t i = 0; i <  array->GetSize(); i++) {
    hw_ok = array->GetHWStatus(i) & AliTOFChannelOnlineStatusArray::kTOFHWOk;
    pulser_ok = !(array->GetPulserStatus(i) & AliTOFChannelOnlineStatusArray::kTOFPulserBad);
    noise_ok = array->GetNoiseStatus(i) & AliTOFChannelOnlineStatusArray::kTOFNoiseOk;
    hStatus_hw_ok->SetBinContent(i + 1, hw_ok);
    hStatus_pulser_ok->SetBinContent(i + 1, pulser_ok);
    hStatus_noise_ok->SetBinContent(i + 1, noise_ok);
    hStatus_active->SetBinContent(i + 1, hw_ok && pulser_ok && noise_ok);
  }
  printf("%d active channels\n", hStatus_active->Integral());

  TFile *fileout = TFile::Open("ChannelStatus.root", "RECREATE");
  hStatus_hw_ok->Write();
  hStatus_pulser_ok->Write();
  hStatus_noise_ok->Write();
  hStatus_active->Write();
  fileout->Close();

}
Commit	Line	Data
	1	CheckActiveChannels(const Char_t *);
	2	CheckActiveChannels(Int_t nrun=-1);
	3	CheckActiveChannelsFromCDBEntry(AliCDBEntry *);
	4
	5	//____________________________________________________________________
	6
	7	CheckActiveChannels(const Char_t *fileName)
	8	{
	9
	10	TFile *file = TFile::Open(fileName);
	11	AliCDBEntry cdbe = (AliCDBEntry )file->Get("AliCDBEntry");
	12	CheckActiveChannelsFromCDBEntry(cdbe);
	13
	14	}
	15
	16	//____________________________________________________________________
	17
	18	CheckActiveChannels(Int_t nrun)
	19	{
	20
	21	TGrid *alien = TGrid::Connect("alien://");
	22	if (!alien \|\| !alien->IsConnected()) return;
	23	AliCDBManager *cdbm = AliCDBManager::Instance();
	24	cdbm->SetDefaultStorage("raw://");
	25	if (nrun==-1)
	26	cdbm->SetRun(80015);
	27	else
	28	cdbm->SetRun(nrun);
	29	AliCDBEntry *cdbe = cdbm->Get("TOF/Calib/Status");
	30	CheckActiveChannelsFromCDBEntry(cdbe);
	31
	32	}
	33
	34	//____________________________________________________________________
	35
	36	CheckActiveChannelsFromCDBEntry(AliCDBEntry *cdbe)
	37	{
	38
	39	AliTOFChannelOnlineStatusArray array = (AliTOFChannelOnlineStatusArray )cdbe->GetObject();
	40	TH1F *hStatus_hw_ok = new TH1F("hStatus_hw_ok", "HW status;index;status ok", array->GetSize(), 0., array->GetSize(););
	41	TH1F *hStatus_pulser_ok = new TH1F("hStatus_pulser_ok", "pulser status;index;status ok", array->GetSize(), 0., array->GetSize(););
	42	TH1F *hStatus_noise_ok = new TH1F("hStatus_noise_ok", "noise status;index;status ok", array->GetSize(), 0., array->GetSize(););
	43	TH1F *hStatus_active = new TH1F("hStatus_active", "active status;index;status ok", array->GetSize(), 0., array->GetSize(););
	44	Bool_t hw_ok, pulser_ok, noise_ok;
	45	for (Int_t i = 0; i < array->GetSize(); i++) {
	46	hw_ok = array->GetHWStatus(i) & AliTOFChannelOnlineStatusArray::kTOFHWOk;
	47	pulser_ok = !(array->GetPulserStatus(i) & AliTOFChannelOnlineStatusArray::kTOFPulserBad);
	48	noise_ok = array->GetNoiseStatus(i) & AliTOFChannelOnlineStatusArray::kTOFNoiseOk;
	49	hStatus_hw_ok->SetBinContent(i + 1, hw_ok);
	50	hStatus_pulser_ok->SetBinContent(i + 1, pulser_ok);
	51	hStatus_noise_ok->SetBinContent(i + 1, noise_ok);
	52	hStatus_active->SetBinContent(i + 1, hw_ok && pulser_ok && noise_ok);
	53	}
	54	printf("%d active channels\n", hStatus_active->Integral());
	55
	56	TFile *fileout = TFile::Open("ChannelStatus.root", "RECREATE");
	57	hStatus_hw_ok->Write();
	58	hStatus_pulser_ok->Write();
	59	hStatus_noise_ok->Write();
	60	hStatus_active->Write();
	61	fileout->Close();
	62
	63	}
	64
	65