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

/////////////////////////////////////////
// Class for SDD digits preprocessing  //
//                                     //
//                                     //
////////////////////////////////////////

#include "AliITSPreprocessorSDD.h"
#include "AliITSCalibrationSDD.h"
#include "AliShuttleInterface.h"
#include "AliCDBMetaData.h"
#include "TObjArray.h"
#include "AliLog.h"

const Int_t AliITSPreprocessorSDD::fgkNumberOfSDD = 260;
const Int_t AliITSPreprocessorSDD::fgkNumberOfChannels = 512;
const TString AliITSPreprocessorSDD::fgkNameHistoPedestals = "hpedestal";
const TString AliITSPreprocessorSDD::fgkNameHistoNoise = "hnoise";
ClassImp(AliITSPreprocessorSDD)


UInt_t AliITSPreprocessorSDD::Process(TMap*/* dcsAliasMap*/){

  //preprocessing. 


  AliCDBMetaData *md1= new AliCDBMetaData(); // metaData describing the object
  md1->SetObjectClassName("AliITSCalibration");
  md1->SetResponsible("Elisabetta Crescio, Francesco Prino");
  md1->SetBeamPeriod(0);
  md1->SetAliRootVersion("head 5 April 2007"); //root version
  md1->SetComment("This is a test");

  TObjArray respSDD(fgkNumberOfSDD);
  respSDD.SetOwner(kFALSE);
  
  Char_t filid[20];
  Float_t baseline,rawnoise,cmn,gain;
  Int_t isgoodan,i,im,is,isgoodmod;
  Int_t numOfBadChannels[fgkNumberOfSDD];
  for(Int_t imod=0;imod<fgkNumberOfSDD;imod++){
    AliITSCalibrationSDD *cal = new AliITSCalibrationSDD("simulated");
    numOfBadChannels[imod]=0;
    Int_t badch[fgkNumberOfChannels];
    for(Int_t isid=0;isid<=1;isid++){
      sprintf(filid,"DAQDAm%03ds%d",imod,isid);
      const char* filenamed= GetFile(kDAQ,filid,"GDC");
      FILE* filed = fopen(filenamed,"read");
      if (filed == 0) {
	AliWarning("File not found"); 
	continue;
      }
      
      fscanf(filed,"%d %d %d\n",&im,&is,&isgoodmod);
      if(!isgoodmod) cal->SetDead();
      for(Int_t ian=0;ian<(fgkNumberOfChannels/2);ian++){
	fscanf(filed,"%d %d %f %f %f %f\n",&i,&isgoodan,&baseline,&rawnoise,&cmn,&gain);
	Int_t ich=ian;
	if(isid==1) ich+=256;
	if(!isgoodan){ 
	  Int_t ibad=numOfBadChannels[imod];
	  numOfBadChannels[imod]++;
	  badch[ibad]=ich;
	}
	cal->SetBaseline(ich,baseline);
	cal->SetNoiseAfterElectronics(ich,rawnoise);	
	//	cal->SetGain(gain,isid,ian/4,ian%4);
      }
      cal->SetDeadChannels(numOfBadChannels[imod]);
      for(Int_t ibad=0;ibad<numOfBadChannels[imod];ibad++){
	cal->SetBadChannel(ibad,badch[ibad]);
      }
      fclose(filed);
    }
    respSDD.Add(cal);
  }

  Bool_t retcode = Store("Calib","Data",&respSDD,md1);


  if(retcode) return 0;
  else return 1;
}
Commit	Line	Data
e56160b8	1	/////////////////////////////////////////
e56160b8	2	// Class for SDD digits preprocessing //
7537d03c	3	// //
7537d03c	4	// //
e56160b8	5	////////////////////////////////////////
	6
	7	#include "AliITSPreprocessorSDD.h"
	8	#include "AliITSCalibrationSDD.h"
	9	#include "AliShuttleInterface.h"
	10	#include "AliCDBMetaData.h"
	11	#include "TObjArray.h"
e56160b8	12	#include "AliLog.h"
e56160b8	13
	14	const Int_t AliITSPreprocessorSDD::fgkNumberOfSDD = 260;
	15	const Int_t AliITSPreprocessorSDD::fgkNumberOfChannels = 512;
7537d03c	16	const TString AliITSPreprocessorSDD::fgkNameHistoPedestals = "hpedestal";
7537d03c	17	const TString AliITSPreprocessorSDD::fgkNameHistoNoise = "hnoise";
e56160b8	18	ClassImp(AliITSPreprocessorSDD)
	19
	20
	21	UInt_t AliITSPreprocessorSDD::Process(TMap/ dcsAliasMap*/){
	22
	23	//preprocessing.
	24
e56160b8	25
	26	AliCDBMetaData *md1= new AliCDBMetaData(); // metaData describing the object
	27	md1->SetObjectClassName("AliITSCalibration");
afa07040	28	md1->SetResponsible("Elisabetta Crescio, Francesco Prino");
e56160b8	29	md1->SetBeamPeriod(0);
afa07040	30	md1->SetAliRootVersion("head 5 April 2007"); //root version
e56160b8	31	md1->SetComment("This is a test");
e56160b8	32
afa07040	33	TObjArray respSDD(fgkNumberOfSDD);
e56160b8	34	respSDD.SetOwner(kFALSE);
e56160b8	35
afa07040	36	Char_t filid[20];
	37	Float_t baseline,rawnoise,cmn,gain;
	38	Int_t isgoodan,i,im,is,isgoodmod;
	39	Int_t numOfBadChannels[fgkNumberOfSDD];
e56160b8	40	for(Int_t imod=0;imod<fgkNumberOfSDD;imod++){
e56160b8	41	AliITSCalibrationSDD *cal = new AliITSCalibrationSDD("simulated");
afa07040	42	numOfBadChannels[imod]=0;
	43	Int_t badch[fgkNumberOfChannels];
	44	for(Int_t isid=0;isid<=1;isid++){
	45	sprintf(filid,"DAQDAm%03ds%d",imod,isid);
	46	const char* filenamed= GetFile(kDAQ,filid,"GDC");
	47	FILE* filed = fopen(filenamed,"read");
	48	if (filed == 0) {
	49	AliWarning("File not found");
	50	continue;
	51	}
	52
	53	fscanf(filed,"%d %d %d\n",&im,&is,&isgoodmod);
	54	if(!isgoodmod) cal->SetDead();
	55	for(Int_t ian=0;ian<(fgkNumberOfChannels/2);ian++){
	56	fscanf(filed,"%d %d %f %f %f %f\n",&i,&isgoodan,&baseline,&rawnoise,&cmn,&gain);
	57	Int_t ich=ian;
	58	if(isid==1) ich+=256;
	59	if(!isgoodan){
	60	Int_t ibad=numOfBadChannels[imod];
	61	numOfBadChannels[imod]++;
	62	badch[ibad]=ich;
	63	}
	64	cal->SetBaseline(ich,baseline);
	65	cal->SetNoiseAfterElectronics(ich,rawnoise);
	66	// cal->SetGain(gain,isid,ian/4,ian%4);
	67	}
	68	cal->SetDeadChannels(numOfBadChannels[imod]);
	69	for(Int_t ibad=0;ibad<numOfBadChannels[imod];ibad++){
	70	cal->SetBadChannel(ibad,badch[ibad]);
	71	}
	72	fclose(filed);
e56160b8	73	}
	74	respSDD.Add(cal);
	75	}
	76
afa07040	77	Bool_t retcode = Store("Calib","Data",&respSDD,md1);
e56160b8	78
e56160b8	79
afa07040	80	if(retcode) return 0;
afa07040	81	else return 1;
e56160b8	82	}