[u/mrichter/AliRoot.git] / PHOS / AliPHOSRcuDA1.cxx

#include "AliPHOSRcuDA1.h"
#include "TString.h"

ClassImp(AliPHOSRcuDA1)

//----------------------------------------------------------------
AliPHOSRcuDA1::AliPHOSRcuDA1(Int_t module, Int_t rcu) : TNamed(),
	   fHistoFile(0),fMod(module),fRCU(rcu),fWriteToFile(kTRUE),fHistoArray()

{
  // Create DA1 ("Calibration DA") object.
  // module is the PHOS module number (0..4).
  // Checks existence of histograms which might have been left
  // from the previous runs to continue their filling.
  // Histogram names: module_iX_iZ_gain for TH2 and  module_iX_iZ for TH1.
  // Root file name: PHOS_ModuleX_RCUY_Calib.root, where X - module number,
  // Y - RCU number. If no RCU specified (rcu<0), file name is simply
  // PHOS_ModuleX_Calib.root.
  
  char name[128]; TString sname;

  if(rcu<0) { 
    sname="PHOS_Module%d_Calib";
    snprintf(name,sname.Length(),sname.Data(),fMod);
  }
  else {
    sname="PHOS_Module%d_RCU%d_Calib";
    snprintf(name,sname.Length(),sname.Data(),fMod,fRCU);
  }

  SetName(name);
  SetTitle("Calibration Detector Algorithm");

  char rootname[128]; 
  TString srootname="%s.root";
  snprintf(rootname,srootname.Length(),srootname.Data(),GetName());

  fHistoFile =  new TFile(rootname,"update");
  fHistoArray.SetName("histo_container");

  char hname[128];
  TH1F* hist1=0;
  TH2F* hist2=0;
  TString shname;

  for(Int_t iX=0; iX<64; iX++) {
    for(Int_t iZ=0; iZ<56; iZ++) {
      
      shname = "%d_%d_%d";
      snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ);

      hist1 = (TH1F*)fHistoFile->Get(hname);
      if(hist1) { 
	fHgLgRatio[iX][iZ] = hist1;
	fHistoArray.Add(hist1);
      }
      else
	fHgLgRatio[iX][iZ] = 0;

      for(Int_t iGain=0; iGain<2; iGain++) {
	
	shname = "%d_%d_%d_%d";
	snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ,iGain);
	
	hist2 = (TH2F*)fHistoFile->Get(hname);
	if(hist2) {
	  fTimeEnergy[iX][iZ][iGain] = hist2;
	  fHistoArray.Add(hist2);
	}
	else
	  fTimeEnergy[iX][iZ][iGain] = 0;
      }

    }
  }
 
}

//-------------------------------------------------------------------
AliPHOSRcuDA1::AliPHOSRcuDA1(Int_t module, Int_t rcu, TObjArray* oldTimeEnergy) :
  TNamed(),fHistoFile(0),fMod(module),fRCU(rcu),fWriteToFile(kFALSE),
  fHistoArray()
{
  // Constructor. 
  // oldTimeEnergy is an array of histograms kept from the previous run.
  // By default the final histograms will not be saved to the root file.

  char name[128]; TString sname;

  if(rcu<0) {
    sname="PHOS_Module%d_Calib";
    snprintf(name,sname.Length(),sname.Data(),fMod);
  }
  else {
    sname="PHOS_Module%d_RCU%d_Calib";
    snprintf(name,sname.Length(),sname.Data(),fMod,fRCU);
  }
  
  SetName(name);
  SetTitle("Calibration Detector Algorithm");
  
  fHistoArray.SetName("histo_container");
  
  char hname[128];
  TH1F* hist1=0;
  TH2F* hist2=0;
  TString shname;

  for(Int_t iX=0; iX<64; iX++) {
    for(Int_t iZ=0; iZ<56; iZ++) {
      
      shname = "%d_%d_%d";
      snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ);
      
      if(oldTimeEnergy)
	hist1 = (TH1F*)oldTimeEnergy->FindObject(hname);
      if(hist1) {
        fHgLgRatio[iX][iZ] = hist1;
	fHistoArray.Add(hist1);
      }
      else
        fHgLgRatio[iX][iZ] = 0;
      
      for(Int_t iGain=0; iGain<2; iGain++) {
	
        shname = "%d_%d_%d_%d";
        snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ,iGain);
	
	if(oldTimeEnergy)
	  hist2 = (TH2F*)oldTimeEnergy->FindObject(hname);
        if(hist2) {
          fTimeEnergy[iX][iZ][iGain] = hist2;
	  fHistoArray.Add(hist2);
        }
        else
          fTimeEnergy[iX][iZ][iGain] = 0;
      }

    }
  }

}

//-------------------------------------------------------------------
AliPHOSRcuDA1::~AliPHOSRcuDA1()
{
  // Destructor
  
  UpdateHistoFile();
  if(fHistoFile) delete fHistoFile;
  
}

//-------------------------------------------------------------------
void AliPHOSRcuDA1::FillHistograms(Float_t e[64][56][2], Float_t t[64][56][2]) 
{
  // Fill energy vs time-of-flight and HG/LG ratio histograms in one event.
  // Energy and TOF data are encoded as e[X][Z][gain] and t[X][Z][gain], 
  // where X(0..63) and Z(0..55) - crystal position in the module, 
  // gain=0 - low gain, gain=1 - high gain.
  // Energy in ADC counts, time in samples.
  // If no energy or time read for particular channel, 
  // the correspondent array entry should be filled by zero.
  // WARNING: this function should be called once per event!

  TString hname;
  TString htitle;

  for(Int_t iX=0; iX<64; iX++) {
    for (Int_t iZ=0; iZ<56; iZ++) {
      
      // HG/LG
      if(e[iX][iZ][0]>10. && e[iX][iZ][1]>10. && e[iX][iZ][1]<900.) {
	
	if(fHgLgRatio[iX][iZ]) {
	  //printf("iX=%d iZ=%d,e[iX][iZ][1]=%.3f,e[iX][iZ][0]=%.3f, t1=%.3f,t0=%.3f\n",
	  // 		 iX,iZ,e[iX][iZ][1],e[iX][iZ][0], t[iX][iZ][1],t[iX][iZ][0]);
	  fHgLgRatio[iX][iZ]->Fill(e[iX][iZ][1]/e[iX][iZ][0]); 
	}
	else
	  {
	    hname.Clear(); htitle.Clear();
	    hname += fMod; hname += "_"; hname += iX; hname += "_"; hname += iZ;
	    htitle += "HG/LG ratio for crystal "; htitle += hname;
	    
	    fHgLgRatio[iX][iZ] = new TH1F(hname,htitle,400,14.,18.);
// 	    printf("iX=%d iZ=%d,e[iX][iZ][1]=%.3f,e[iX][iZ][0]=%.3f\n",
// 		   iX,iZ,e[iX][iZ][1],e[iX][iZ][0]);
	    fHgLgRatio[iX][iZ]->Fill(e[iX][iZ][1]/e[iX][iZ][0]);
	    fHistoArray.Add(fHgLgRatio[iX][iZ]);
	  }
      }
	  
      // Energy vs TOF
      for(Int_t iGain=0; iGain<2; iGain++) {

	if(e[iX][iZ][iGain]<10) continue;
	if(!(t[iX][iZ][iGain]>0)) continue;

	if(fTimeEnergy[iX][iZ][iGain]) 
	  fTimeEnergy[iX][iZ][iGain]->Fill(e[iX][iZ][iGain],t[iX][iZ][iGain]);
	else 
	  {
	    hname.Clear(); htitle.Clear();
	    hname += fMod; hname += "_"; hname += iX; hname += "_"; hname += iZ;
	    htitle = "Energy vs TOF for crystal "; htitle += hname; htitle += " and gain "; htitle += iGain;
	    hname += "_"; hname += iGain;
	    
	    fTimeEnergy[iX][iZ][iGain] = new TH2F(hname,htitle,100,0.,1024.,50,0.,50.);
	    fTimeEnergy[iX][iZ][iGain]->Fill(e[iX][iZ][iGain],t[iX][iZ][iGain]);
	    fHistoArray.Add(fTimeEnergy[iX][iZ][iGain]);
	  }
      }
      
    }
  }
  
}

//-------------------------------------------------------------------
void AliPHOSRcuDA1::UpdateHistoFile()
{
  // Write histograms to file

  if(!fWriteToFile) return;
  if(!fHistoFile) return;
  if(!fHistoFile->IsOpen()) return;

  fHistoArray.Write();
  fHistoFile->Purge();

}

//-----------------------------------------------------------------------
void AliPHOSRcuDA1::SetWriteToFile(Bool_t write)
{
  if(!write) { 
    fWriteToFile = write;
    return;
  }

  if(!fHistoFile) {
    TString rootname(GetName());
    rootname += ".root";
    fHistoFile =  new TFile(rootname.Data(),"update");
  } 
  
  fWriteToFile = write;
}
Commit	Line	Data
354c9df9	1	#include "AliPHOSRcuDA1.h"
	2	#include "TString.h"
	3
	4	ClassImp(AliPHOSRcuDA1)
	5
	6	//----------------------------------------------------------------
	7	AliPHOSRcuDA1::AliPHOSRcuDA1(Int_t module, Int_t rcu) : TNamed(),
	8	fHistoFile(0),fMod(module),fRCU(rcu),fWriteToFile(kTRUE),fHistoArray()
	9
	10	{
	11	// Create DA1 ("Calibration DA") object.
	12	// module is the PHOS module number (0..4).
	13	// Checks existence of histograms which might have been left
	14	// from the previous runs to continue their filling.
	15	// Histogram names: module_iX_iZ_gain for TH2 and module_iX_iZ for TH1.
4056a2e7	16	// Root file name: PHOS_ModuleX_RCUY_Calib.root, where X - module number,
	17	// Y - RCU number. If no RCU specified (rcu<0), file name is simply
	18	// PHOS_ModuleX_Calib.root.
354c9df9	19
6ba1dd81	20	char name[128]; TString sname;
	21
	22	if(rcu<0) {
	23	sname="PHOS_Module%d_Calib";
	24	snprintf(name,sname.Length(),sname.Data(),fMod);
	25	}
	26	else {
	27	sname="PHOS_Module%d_RCU%d_Calib";
	28	snprintf(name,sname.Length(),sname.Data(),fMod,fRCU);
	29	}
354c9df9	30
6ba1dd81	31	SetName(name);
354c9df9	32	SetTitle("Calibration Detector Algorithm");
354c9df9	33
6ba1dd81	34	char rootname[128];
	35	TString srootname="%s.root";
	36	snprintf(rootname,srootname.Length(),srootname.Data(),GetName());
354c9df9	37
	38	fHistoFile = new TFile(rootname,"update");
	39	fHistoArray.SetName("histo_container");
	40
354c9df9	41	char hname[128];
	42	TH1F* hist1=0;
	43	TH2F* hist2=0;
6ba1dd81	44	TString shname;
354c9df9	45
	46	for(Int_t iX=0; iX<64; iX++) {
	47	for(Int_t iZ=0; iZ<56; iZ++) {
6ba1dd81	48
	49	shname = "%d_%d_%d";
	50	snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ);
354c9df9	51
52fccb4e	52	hist1 = (TH1F*)fHistoFile->Get(hname);
354c9df9	53	if(hist1) {
	54	fHgLgRatio[iX][iZ] = hist1;
	55	fHistoArray.Add(hist1);
	56	}
	57	else
	58	fHgLgRatio[iX][iZ] = 0;
	59
	60	for(Int_t iGain=0; iGain<2; iGain++) {
6ba1dd81	61
	62	shname = "%d_%d_%d_%d";
	63	snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ,iGain);
	64
52fccb4e	65	hist2 = (TH2F*)fHistoFile->Get(hname);
354c9df9	66	if(hist2) {
	67	fTimeEnergy[iX][iZ][iGain] = hist2;
	68	fHistoArray.Add(hist2);
	69	}
	70	else
	71	fTimeEnergy[iX][iZ][iGain] = 0;
	72	}
	73
	74	}
	75	}
	76
	77	}
	78
	79	//-------------------------------------------------------------------
	80	AliPHOSRcuDA1::AliPHOSRcuDA1(Int_t module, Int_t rcu, TObjArray* oldTimeEnergy) :
	81	TNamed(),fHistoFile(0),fMod(module),fRCU(rcu),fWriteToFile(kFALSE),
	82	fHistoArray()
	83	{
	84	// Constructor.
	85	// oldTimeEnergy is an array of histograms kept from the previous run.
	86	// By default the final histograms will not be saved to the root file.
	87
6ba1dd81	88	char name[128]; TString sname;
	89
	90	if(rcu<0) {
	91	sname="PHOS_Module%d_Calib";
	92	snprintf(name,sname.Length(),sname.Data(),fMod);
	93	}
	94	else {
	95	sname="PHOS_Module%d_RCU%d_Calib";
	96	snprintf(name,sname.Length(),sname.Data(),fMod,fRCU);
	97	}
354c9df9	98
6ba1dd81	99	SetName(name);
354c9df9	100	SetTitle("Calibration Detector Algorithm");
6ba1dd81	101
354c9df9	102	fHistoArray.SetName("histo_container");
	103
	104	char hname[128];
	105	TH1F* hist1=0;
	106	TH2F* hist2=0;
6ba1dd81	107	TString shname;
354c9df9	108
	109	for(Int_t iX=0; iX<64; iX++) {
	110	for(Int_t iZ=0; iZ<56; iZ++) {
	111
6ba1dd81	112	shname = "%d_%d_%d";
	113	snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ);
	114
354c9df9	115	if(oldTimeEnergy)
	116	hist1 = (TH1F*)oldTimeEnergy->FindObject(hname);
	117	if(hist1) {
	118	fHgLgRatio[iX][iZ] = hist1;
	119	fHistoArray.Add(hist1);
	120	}
	121	else
	122	fHgLgRatio[iX][iZ] = 0;
	123
	124	for(Int_t iGain=0; iGain<2; iGain++) {
6ba1dd81	125
	126	shname = "%d_%d_%d_%d";
	127	snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ,iGain);
	128
354c9df9	129	if(oldTimeEnergy)
	130	hist2 = (TH2F*)oldTimeEnergy->FindObject(hname);
	131	if(hist2) {
	132	fTimeEnergy[iX][iZ][iGain] = hist2;
	133	fHistoArray.Add(hist2);
	134	}
	135	else
	136	fTimeEnergy[iX][iZ][iGain] = 0;
	137	}
	138
	139	}
	140	}
	141
	142	}
	143
	144	//-------------------------------------------------------------------
	145	AliPHOSRcuDA1::~AliPHOSRcuDA1()
	146	{
	147	// Destructor
	148
	149	UpdateHistoFile();
	150	if(fHistoFile) delete fHistoFile;
	151
	152	}
	153
	154	//-------------------------------------------------------------------
	155	void AliPHOSRcuDA1::FillHistograms(Float_t e[64][56][2], Float_t t[64][56][2])
	156	{
	157	// Fill energy vs time-of-flight and HG/LG ratio histograms in one event.
	158	// Energy and TOF data are encoded as e[X][Z][gain] and t[X][Z][gain],
	159	// where X(0..63) and Z(0..55) - crystal position in the module,
	160	// gain=0 - low gain, gain=1 - high gain.
	161	// Energy in ADC counts, time in samples.
	162	// If no energy or time read for particular channel,
	163	// the correspondent array entry should be filled by zero.
	164	// WARNING: this function should be called once per event!
	165
7e6cd815	166	TString hname;
7e6cd815	167	TString htitle;
354c9df9	168
	169	for(Int_t iX=0; iX<64; iX++) {
	170	for (Int_t iZ=0; iZ<56; iZ++) {
	171
	172	// HG/LG
2d03aabf	173	if(e[iX][iZ][0]>10. && e[iX][iZ][1]>10. && e[iX][iZ][1]<900.) {
354c9df9	174
	175	if(fHgLgRatio[iX][iZ]) {
	176	//printf("iX=%d iZ=%d,e[iX][iZ][1]=%.3f,e[iX][iZ][0]=%.3f, t1=%.3f,t0=%.3f\n",
	177	// iX,iZ,e[iX][iZ][1],e[iX][iZ][0], t[iX][iZ][1],t[iX][iZ][0]);
	178	fHgLgRatio[iX][iZ]->Fill(e[iX][iZ][1]/e[iX][iZ][0]);
	179	}
	180	else
	181	{
7e6cd815	182	hname.Clear(); htitle.Clear();
	183	hname += fMod; hname += "_"; hname += iX; hname += "_"; hname += iZ;
	184	htitle += "HG/LG ratio for crystal "; htitle += hname;
	185
	186	fHgLgRatio[iX][iZ] = new TH1F(hname,htitle,400,14.,18.);
354c9df9	187	// printf("iX=%d iZ=%d,e[iX][iZ][1]=%.3f,e[iX][iZ][0]=%.3f\n",
	188	// iX,iZ,e[iX][iZ][1],e[iX][iZ][0]);
	189	fHgLgRatio[iX][iZ]->Fill(e[iX][iZ][1]/e[iX][iZ][0]);
	190	fHistoArray.Add(fHgLgRatio[iX][iZ]);
	191	}
	192	}
	193
	194	// Energy vs TOF
	195	for(Int_t iGain=0; iGain<2; iGain++) {
	196
	197	if(e[iX][iZ][iGain]<10) continue;
	198	if(!(t[iX][iZ][iGain]>0)) continue;
	199
	200	if(fTimeEnergy[iX][iZ][iGain])
	201	fTimeEnergy[iX][iZ][iGain]->Fill(e[iX][iZ][iGain],t[iX][iZ][iGain]);
7e6cd815	202	else
	203	{
	204	hname.Clear(); htitle.Clear();
	205	hname += fMod; hname += "_"; hname += iX; hname += "_"; hname += iZ;
	206	htitle = "Energy vs TOF for crystal "; htitle += hname; htitle += " and gain "; htitle += iGain;
	207	hname += "_"; hname += iGain;
	208
	209	fTimeEnergy[iX][iZ][iGain] = new TH2F(hname,htitle,100,0.,1024.,50,0.,50.);
	210	fTimeEnergy[iX][iZ][iGain]->Fill(e[iX][iZ][iGain],t[iX][iZ][iGain]);
	211	fHistoArray.Add(fTimeEnergy[iX][iZ][iGain]);
	212	}
354c9df9	213	}
7e6cd815	214
354c9df9	215	}
354c9df9	216	}
7e6cd815	217
354c9df9	218	}
	219
	220	//-------------------------------------------------------------------
	221	void AliPHOSRcuDA1::UpdateHistoFile()
	222	{
	223	// Write histograms to file
	224
	225	if(!fWriteToFile) return;
	226	if(!fHistoFile) return;
	227	if(!fHistoFile->IsOpen()) return;
	228
f39fcb60	229	fHistoArray.Write();
354c9df9	230	fHistoFile->Purge();
	231
	232	}
	233
	234	//-----------------------------------------------------------------------
	235	void AliPHOSRcuDA1::SetWriteToFile(Bool_t write)
	236	{
	237	if(!write) {
	238	fWriteToFile = write;
	239	return;
	240	}
	241
	242	if(!fHistoFile) {
7e6cd815	243	TString rootname(GetName());
	244	rootname += ".root";
	245	fHistoFile = new TFile(rootname.Data(),"update");
354c9df9	246	}
	247
	248	fWriteToFile = write;
	249	}