[u/mrichter/AliRoot.git] / ITS / AliITSOnlineSDDTP.cxx

/**************************************************************************
 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
#include <TFile.h>
#include "AliITSOnlineSDDTP.h"
#include "AliLog.h"
#include <TH2F.h>
#include <TMath.h>


///////////////////////////////////////////////////////////////////
//                                                               //
// Implemetation of the class SDD Test Pulse analysis            //
// Origin: F.Prino, Torino, prino@to.infn.it                     //
//                                                               //
///////////////////////////////////////////////////////////////////


ClassImp(AliITSOnlineSDDTP)

//______________________________________________________________________
AliITSOnlineSDDTP::AliITSOnlineSDDTP():AliITSOnlineSDD(),fDAC(0.),fLowThreshold(0),fHighThreshold(0),fNSigmaGain(0.),fNSigmaNoise(0.)
{
  // default constructor
  Reset();
  SetNSigmaGain();
  SetNSigmaNoise();
}
//______________________________________________________________________
AliITSOnlineSDDTP::AliITSOnlineSDDTP(Int_t nddl, Int_t ncarlos, Int_t sid, Float_t xDAC):AliITSOnlineSDD(nddl,ncarlos,sid),fDAC(xDAC),fLowThreshold(0),fHighThreshold(0),fNSigmaGain(0.),fNSigmaNoise(0.)
{
  // standard constructor
  Reset();
  SetNSigmaGain();
  SetNSigmaNoise();
}
//______________________________________________________________________
AliITSOnlineSDDTP::~AliITSOnlineSDDTP(){
  // Destructor
}
//______________________________________________________________________
void AliITSOnlineSDDTP::Reset(){
  //
  for(Int_t i=0;i<fgkNAnodes;i++){
    fNEvents[i]=0;
    fGoodAnode[i]=1;
    fBaseline[i]=0.;
    fCMN[i]=0.;
    fRawNoise[i]=0.;
    fCorrNoise[i]=0.;
    fSumTPPeak[i]=0.;
    fTPPos[i]=0.;
  }
  ReadBaselines();
}

//______________________________________________________________________
void AliITSOnlineSDDTP::AddEvent(TH2F* hrawd){
  // 
  for(Int_t ian=0;ian<fgkNAnodes;ian++){
    Float_t auxmax=0.;
    Int_t auxtb=0;
    if(!fGoodAnode[ian]) continue;
    for(Int_t itb=fFirstGoodTB;itb<=fLastGoodTB;itb++){
      Float_t cnt=hrawd->GetBinContent(itb+1,ian+1);
      if(cnt>auxmax){ 
	auxmax=cnt;
	auxtb=itb;
      }
    }
    if(auxmax>fBaseline[ian]+fNSigmaNoise*fRawNoise[ian]){
      fSumTPPeak[ian]+=auxmax-fBaseline[ian];
      fTPPos[ian]+=auxtb;
      fNEvents[ian]++;
    }
  }
}
//______________________________________________________________________
void AliITSOnlineSDDTP::ReadBaselines(){
  // assume baselines and good anodes are taken from previous run
  Char_t basfilnam[100];
  sprintf(basfilnam,"SDDbase_step2_ddl%02dc%02d_sid%d.data",fDDL,fCarlos,fSide);
  FILE* basf=fopen(basfilnam,"r");
  if(basf==0){
    AliWarning(Form("Baseline file not present (ddl %d  carlos %d side %d, Set all baselines to 20",fDDL,fCarlos,fSide));
    for(Int_t ian=0;ian<fgkNAnodes;ian++){ 
      fBaseline[ian]=20.;
      fEqBaseline[ian]=20;
      fOffsetBaseline[ian]=0;
      fGoodAnode[ian]=1;
    }
    return;
  }
  fscanf(basf,"%d\n",&fHighThreshold);
  fscanf(basf,"%d\n",&fLowThreshold);
  Int_t n,ok,eqbase,offbase;
  Float_t base,rms,cmn,corrnoi;
  for(Int_t ian=0;ian<fgkNAnodes;ian++){
    fscanf(basf,"%d %d %f %d %d %f %f %f\n",&n,&ok,&base,&eqbase,&offbase,&rms,&cmn,&corrnoi);
    fBaseline[ian]=base;
    fEqBaseline[ian]=eqbase;
    fOffsetBaseline[ian]=offbase;
    fGoodAnode[ian]=ok;
    fRawNoise[ian]=rms;
    fCMN[ian]=cmn;
    fCorrNoise[ian]=corrnoi;
  }
  fclose(basf);
}

//______________________________________________________________________
Bool_t AliITSOnlineSDDTP::IsModuleGood() const{
  //
  // Check if there is at least 1 good anode
  //
  for(Int_t ian=0;ian<fgkNAnodes;ian++){
    if(fGoodAnode[ian]) return kTRUE;
  }
  return kFALSE;
}
//______________________________________________________________________
void AliITSOnlineSDDTP::ValidateAnodes(){
  //
  Float_t meang,rmsg;
  StatGain(meang,rmsg);
  Float_t lowlim=meang-fNSigmaGain*rmsg;
  Float_t hilim=meang+fNSigmaGain*rmsg;

  for(Int_t ian=0;ian<fgkNAnodes;ian++){
    if(!fGoodAnode[ian]) continue;
    if(GetChannelGain(ian)<lowlim||GetChannelGain(ian)>hilim) fGoodAnode[ian]=0;
  }
}


//______________________________________________________________________
void AliITSOnlineSDDTP::StatGain(Float_t &mean, Float_t  &rms){
  //
  Float_t sum=0.,sumq=0.;
  Int_t cnt=0;
  for(Int_t ian=0;ian<fgkNAnodes;ian++){
    if(!fGoodAnode[ian]) continue;
    if(fNEvents[ian]==0) continue;
    sum+=GetChannelGain(ian);
    sumq+=TMath::Power(GetChannelGain(ian),2);
    cnt++;
  }
  if(cnt>0){ 
    mean=sum/(Float_t)cnt;
    Float_t variance=sumq/(Float_t)cnt-mean*mean;
    if(variance>0.) rms=TMath::Sqrt(variance);
    else rms = 0;
  }else{ 
    mean=0.;
    rms=0.;
  }
  return;
}

//______________________________________________________________________
void AliITSOnlineSDDTP::WriteToASCII(){
  //
  Char_t outfilnam[100];
  sprintf(outfilnam,"SDDbase_ddl%02dc%02d_sid%d.data",fDDL,fCarlos,fSide);
  FILE* outf=fopen(outfilnam,"w");
  fprintf(outf,"%d %d %d\n",fCarlos,fSide,IsModuleGood());
  fprintf(outf,"%d\n",fHighThreshold);
  fprintf(outf,"%d\n",fLowThreshold);
  for(Int_t ian=0;ian<fgkNAnodes;ian++){
    fprintf(outf,"%d %d %8.3f %d %d %8.3f %8.3f %8.3f %8.3f\n",ian,IsAnodeGood(ian),GetAnodeBaseline(ian),GetAnodeEqualizedBaseline(ian),GetAnodeBaselineOffset(ian),GetAnodeRawNoise(ian),GetAnodeCommonMode(ian),GetAnodeCorrNoise(ian),GetChannelGain(ian));
  }
  fclose(outf);  
}
//______________________________________________________________________
Bool_t AliITSOnlineSDDTP::WriteToROOT(TFile *fil){
  //
  if(fil==0){ 
    AliWarning("Invalid pointer to ROOT file");
    return kFALSE;    
  }
  Char_t hisnam[20];
  fil->cd();
  sprintf(hisnam,"hgood%02dc%02ds%d",fDDL,fCarlos,fSide);
  TH1F hgood(hisnam,"",256,-0.5,255.5);
  sprintf(hisnam,"hbase%02dc%02ds%d",fDDL,fCarlos,fSide);
  TH1F hbase(hisnam,"",256,-0.5,255.5);
  sprintf(hisnam,"hnois%02dc%02ds%d",fDDL,fCarlos,fSide);
  TH1F hnois(hisnam,"",256,-0.5,255.5);
  sprintf(hisnam,"hcmn%02dc%02ds%d",fDDL,fCarlos,fSide);
  TH1F hcmn(hisnam,"",256,-0.5,255.5);
  sprintf(hisnam,"hcorn%02dc%02ds%d",fDDL,fCarlos,fSide);
  TH1F hcorn(hisnam,"",256,-0.5,255.5);
  sprintf(hisnam,"hgain%02dc%02ds%d",fDDL,fCarlos,fSide);
  TH1F hgain(hisnam,"",256,-0.5,255.5);
  sprintf(hisnam,"htptb%02dc%02ds%d",fDDL,fCarlos,fSide);
  TH1F htptb(hisnam,"",256,-0.5,255.5);
  for(Int_t ian=0;ian<fgkNAnodes;ian++){
    hgood.SetBinContent(ian+1,float(IsAnodeGood(ian)));
    hbase.SetBinContent(ian+1,GetAnodeBaseline(ian));
    hnois.SetBinContent(ian+1,GetAnodeRawNoise(ian));
    hcmn.SetBinContent(ian+1,GetAnodeCommonMode(ian));
    hcorn.SetBinContent(ian+1,GetAnodeCorrNoise(ian));
    hgain.SetBinContent(ian+1,GetChannelGain(ian));
    htptb.SetBinContent(ian+1,GetTimeBinTPPeak(ian));
  }
  hgood.Write();
  hbase.Write();
  hnois.Write();
  hcmn.Write();
  hcorn.Write();
  hgain.Write();
  htptb.Write();
  return kTRUE;
}
Commit	Line	Data
348f80b7	1	/**************************************************************************
	2	* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
3083967f	15	#include <TFile.h>
348f80b7	16	#include "AliITSOnlineSDDTP.h"
7f347565	17	#include "AliLog.h"
348f80b7	18	#include <TH2F.h>
	19	#include <TMath.h>
	20
	21
	22	///////////////////////////////////////////////////////////////////
	23	// //
	24	// Implemetation of the class SDD Test Pulse analysis //
	25	// Origin: F.Prino, Torino, prino@to.infn.it //
	26	// //
	27	///////////////////////////////////////////////////////////////////
	28
	29
	30	ClassImp(AliITSOnlineSDDTP)
	31
	32	//______________________________________________________________________
e7610d19	33	AliITSOnlineSDDTP::AliITSOnlineSDDTP():AliITSOnlineSDD(),fDAC(0.),fLowThreshold(0),fHighThreshold(0),fNSigmaGain(0.),fNSigmaNoise(0.)
348f80b7	34	{
	35	// default constructor
	36	Reset();
	37	SetNSigmaGain();
750296dd	38	SetNSigmaNoise();
348f80b7	39	}
348f80b7	40	//______________________________________________________________________
e7610d19	41	AliITSOnlineSDDTP::AliITSOnlineSDDTP(Int_t nddl, Int_t ncarlos, Int_t sid, Float_t xDAC):AliITSOnlineSDD(nddl,ncarlos,sid),fDAC(xDAC),fLowThreshold(0),fHighThreshold(0),fNSigmaGain(0.),fNSigmaNoise(0.)
348f80b7	42	{
	43	// standard constructor
	44	Reset();
	45	SetNSigmaGain();
750296dd	46	SetNSigmaNoise();
348f80b7	47	}
	48	//______________________________________________________________________
	49	AliITSOnlineSDDTP::~AliITSOnlineSDDTP(){
	50	// Destructor
	51	}
	52	//______________________________________________________________________
	53	void AliITSOnlineSDDTP::Reset(){
beb262b4	54	//
348f80b7	55	for(Int_t i=0;i<fgkNAnodes;i++){
750296dd	56	fNEvents[i]=0;
348f80b7	57	fGoodAnode[i]=1;
348f80b7	58	fBaseline[i]=0.;
3083967f	59	fCMN[i]=0.;
	60	fRawNoise[i]=0.;
	61	fCorrNoise[i]=0.;
348f80b7	62	fSumTPPeak[i]=0.;
	63	fTPPos[i]=0.;
	64	}
	65	ReadBaselines();
	66	}
	67
	68	//______________________________________________________________________
	69	void AliITSOnlineSDDTP::AddEvent(TH2F* hrawd){
	70	//
348f80b7	71	for(Int_t ian=0;ian<fgkNAnodes;ian++){
	72	Float_t auxmax=0.;
	73	Int_t auxtb=0;
	74	if(!fGoodAnode[ian]) continue;
750296dd	75	for(Int_t itb=fFirstGoodTB;itb<=fLastGoodTB;itb++){
348f80b7	76	Float_t cnt=hrawd->GetBinContent(itb+1,ian+1);
	77	if(cnt>auxmax){
	78	auxmax=cnt;
	79	auxtb=itb;
	80	}
	81	}
750296dd	82	if(auxmax>fBaseline[ian]+fNSigmaNoise*fRawNoise[ian]){
	83	fSumTPPeak[ian]+=auxmax-fBaseline[ian];
	84	fTPPos[ian]+=auxtb;
	85	fNEvents[ian]++;
	86	}
348f80b7	87	}
	88	}
	89	//______________________________________________________________________
	90	void AliITSOnlineSDDTP::ReadBaselines(){
	91	// assume baselines and good anodes are taken from previous run
	92	Char_t basfilnam[100];
979b5a5f	93	sprintf(basfilnam,"SDDbase_step2_ddl%02dc%02d_sid%d.data",fDDL,fCarlos,fSide);
348f80b7	94	FILE* basf=fopen(basfilnam,"r");
348f80b7	95	if(basf==0){
2f0c3b74	96	AliWarning(Form("Baseline file not present (ddl %d carlos %d side %d, Set all baselines to 20",fDDL,fCarlos,fSide));
348f80b7	97	for(Int_t ian=0;ian<fgkNAnodes;ian++){
2f0c3b74	98	fBaseline[ian]=20.;
2f0c3b74	99	fEqBaseline[ian]=20;
d8303ac8	100	fOffsetBaseline[ian]=0;
348f80b7	101	fGoodAnode[ian]=1;
	102	}
	103	return;
	104	}
e7610d19	105	fscanf(basf,"%d\n",&fHighThreshold);
e7610d19	106	fscanf(basf,"%d\n",&fLowThreshold);
d8303ac8	107	Int_t n,ok,eqbase,offbase;
348f80b7	108	Float_t base,rms,cmn,corrnoi;
348f80b7	109	for(Int_t ian=0;ian<fgkNAnodes;ian++){
d8303ac8	110	fscanf(basf,"%d %d %f %d %d %f %f %f\n",&n,&ok,&base,&eqbase,&offbase,&rms,&cmn,&corrnoi);
348f80b7	111	fBaseline[ian]=base;
d8303ac8	112	fEqBaseline[ian]=eqbase;
d8303ac8	113	fOffsetBaseline[ian]=offbase;
348f80b7	114	fGoodAnode[ian]=ok;
3083967f	115	fRawNoise[ian]=rms;
	116	fCMN[ian]=cmn;
	117	fCorrNoise[ian]=corrnoi;
348f80b7	118	}
	119	fclose(basf);
	120	}
	121
3083967f	122	//______________________________________________________________________
	123	Bool_t AliITSOnlineSDDTP::IsModuleGood() const{
	124	//
	125	// Check if there is at least 1 good anode
	126	//
	127	for(Int_t ian=0;ian<fgkNAnodes;ian++){
	128	if(fGoodAnode[ian]) return kTRUE;
	129	}
	130	return kFALSE;
	131	}
348f80b7	132	//______________________________________________________________________
348f80b7	133	void AliITSOnlineSDDTP::ValidateAnodes(){
beb262b4	134	//
348f80b7	135	Float_t meang,rmsg;
348f80b7	136	StatGain(meang,rmsg);
348f80b7	137	Float_t lowlim=meang-fNSigmaGain*rmsg;
	138	Float_t hilim=meang+fNSigmaGain*rmsg;
	139
	140	for(Int_t ian=0;ian<fgkNAnodes;ian++){
	141	if(!fGoodAnode[ian]) continue;
	142	if(GetChannelGain(ian)<lowlim\|\|GetChannelGain(ian)>hilim) fGoodAnode[ian]=0;
	143	}
	144	}
	145
	146
	147	//______________________________________________________________________
	148	void AliITSOnlineSDDTP::StatGain(Float_t &mean, Float_t &rms){
beb262b4	149	//
348f80b7	150	Float_t sum=0.,sumq=0.;
	151	Int_t cnt=0;
	152	for(Int_t ian=0;ian<fgkNAnodes;ian++){
	153	if(!fGoodAnode[ian]) continue;
750296dd	154	if(fNEvents[ian]==0) continue;
348f80b7	155	sum+=GetChannelGain(ian);
	156	sumq+=TMath::Power(GetChannelGain(ian),2);
	157	cnt++;
	158	}
	159	if(cnt>0){
	160	mean=sum/(Float_t)cnt;
9b8481e1	161	Float_t variance=sumq/(Float_t)cnt-mean*mean;
	162	if(variance>0.) rms=TMath::Sqrt(variance);
	163	else rms = 0;
348f80b7	164	}else{
	165	mean=0.;
	166	rms=0.;
	167	}
	168	return;
	169	}
	170
	171	//______________________________________________________________________
e44f571c	172	void AliITSOnlineSDDTP::WriteToASCII(){
348f80b7	173	//
348f80b7	174	Char_t outfilnam[100];
979b5a5f	175	sprintf(outfilnam,"SDDbase_ddl%02dc%02d_sid%d.data",fDDL,fCarlos,fSide);
348f80b7	176	FILE* outf=fopen(outfilnam,"w");
979b5a5f	177	fprintf(outf,"%d %d %d\n",fCarlos,fSide,IsModuleGood());
e7610d19	178	fprintf(outf,"%d\n",fHighThreshold);
e7610d19	179	fprintf(outf,"%d\n",fLowThreshold);
348f80b7	180	for(Int_t ian=0;ian<fgkNAnodes;ian++){
d8303ac8	181	fprintf(outf,"%d %d %8.3f %d %d %8.3f %8.3f %8.3f %8.3f\n",ian,IsAnodeGood(ian),GetAnodeBaseline(ian),GetAnodeEqualizedBaseline(ian),GetAnodeBaselineOffset(ian),GetAnodeRawNoise(ian),GetAnodeCommonMode(ian),GetAnodeCorrNoise(ian),GetChannelGain(ian));
348f80b7	182	}
	183	fclose(outf);
	184	}
3083967f	185	//______________________________________________________________________
	186	Bool_t AliITSOnlineSDDTP::WriteToROOT(TFile *fil){
	187	//
	188	if(fil==0){
	189	AliWarning("Invalid pointer to ROOT file");
	190	return kFALSE;
	191	}
	192	Char_t hisnam[20];
	193	fil->cd();
979b5a5f	194	sprintf(hisnam,"hgood%02dc%02ds%d",fDDL,fCarlos,fSide);
3083967f	195	TH1F hgood(hisnam,"",256,-0.5,255.5);
979b5a5f	196	sprintf(hisnam,"hbase%02dc%02ds%d",fDDL,fCarlos,fSide);
3083967f	197	TH1F hbase(hisnam,"",256,-0.5,255.5);
979b5a5f	198	sprintf(hisnam,"hnois%02dc%02ds%d",fDDL,fCarlos,fSide);
3083967f	199	TH1F hnois(hisnam,"",256,-0.5,255.5);
979b5a5f	200	sprintf(hisnam,"hcmn%02dc%02ds%d",fDDL,fCarlos,fSide);
3083967f	201	TH1F hcmn(hisnam,"",256,-0.5,255.5);
979b5a5f	202	sprintf(hisnam,"hcorn%02dc%02ds%d",fDDL,fCarlos,fSide);
3083967f	203	TH1F hcorn(hisnam,"",256,-0.5,255.5);
979b5a5f	204	sprintf(hisnam,"hgain%02dc%02ds%d",fDDL,fCarlos,fSide);
3083967f	205	TH1F hgain(hisnam,"",256,-0.5,255.5);
a7996467	206	sprintf(hisnam,"htptb%02dc%02ds%d",fDDL,fCarlos,fSide);
a7996467	207	TH1F htptb(hisnam,"",256,-0.5,255.5);
3083967f	208	for(Int_t ian=0;ian<fgkNAnodes;ian++){
	209	hgood.SetBinContent(ian+1,float(IsAnodeGood(ian)));
	210	hbase.SetBinContent(ian+1,GetAnodeBaseline(ian));
	211	hnois.SetBinContent(ian+1,GetAnodeRawNoise(ian));
	212	hcmn.SetBinContent(ian+1,GetAnodeCommonMode(ian));
	213	hcorn.SetBinContent(ian+1,GetAnodeCorrNoise(ian));
	214	hgain.SetBinContent(ian+1,GetChannelGain(ian));
a7996467	215	htptb.SetBinContent(ian+1,GetTimeBinTPPeak(ian));
3083967f	216	}
	217	hgood.Write();
	218	hbase.Write();
	219	hnois.Write();
	220	hcmn.Write();
	221	hcorn.Write();
	222	hgain.Write();
a7996467	223	htptb.Write();
3083967f	224	return kTRUE;
3083967f	225	}
348f80b7	226