[u/mrichter/AliRoot.git] / HMPID / AliHMPIDCalib.cxx

//  **************************************************************************
//  * Copyright(c) 1998-1999, 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 "AliHMPIDCalib.h" //class header
#include "TTreePlayer.h"
#include <fstream>
#include <TTree.h>
//#include "AliHMPIDDigit.h"

ClassImp(AliHMPIDCalib) 


//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
AliHMPIDCalib::AliHMPIDCalib():
   fPedTree(0x0),
   fa(-1),
   fd(-1),
   fr(-1),
   fq(-1),
   fPedHisto(0) 
{
  Init();
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
AliHMPIDCalib::~AliHMPIDCalib()
{
  //destructor
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliHMPIDCalib::Init()
{
  for(Int_t iddl=0;iddl<11;iddl++) faddl[iddl]=kFALSE;
  //
  // Called to initialize the TTree in which the raw data words will be stored 
  //
  fPedHisto=new TH1F("fPedHisto","Temporary pedestal",4096,-0.5,4095.5);                                //init pedestal histo
  fPedTree = new TTree("PedTree","HMPID Pedestal Tree");                                                //init pedestal tree
  fPedTree->Branch("diladd",&fa,"diladd/I");
  fPedTree->Branch("dilnum",&fd,"dilnum/I");
  fPedTree->Branch("dilrow",&fr,"dilrow/I");
  fPedTree->Branch("qdc"   ,&fq,"qdc/F   ");
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliHMPIDCalib::FillPedestal(Int_t ddl,Int_t row, Int_t dil,Int_t adr,Int_t q)
{
  //
  //Called from the HMPIDda and fills the pedestal tree 
  //
  fq=q;
  fr=row;
  fd=dil;
  fa=adr;
  faddl[ddl]=kTRUE; 
  if(fq>-1) fPedTree->Fill();    
}//FillPedestal()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Bool_t AliHMPIDCalib::CalcPedestal(Int_t ddl, Char_t* name)    
{
  //
  //Calculate pedestal for each pad  
  //
  
  ofstream out;                                           //to write the pedestal text files
  Double_t mean,sigma;
  Int_t inhard;
  if(faddl[ddl]==kFALSE) return kFALSE;                   //if ddl is missing no ped file is created (and also for LDC selection). Check with Paolo what he checks for?!
    out.open(name);
    for(Int_t row=1; row < 25; row++){
     for(Int_t dil=1; dil < 11; dil++){
      for(Int_t adr=0; adr < 48; adr++){
          fPedHisto->Reset();
          fPedTree->Draw("fq>>fPedHisto",Form("fr==%d&&fd==%d&&fa==%d",row,dil,adr));
          mean=fPedHisto->GetMean();
          sigma=fPedHisto->GetRMS();
          inhard=((Int_t(mean))<<9)+Int_t(mean+3*sigma);
	  out << Form("%2i %2i %2i %5.2f %5.2f %x\n",row,dil,adr,mean,sigma,inhard);
        }//adr
      }//dil
    }//row
    out.close();                                          //write pedestal file
  return kTRUE;
}//CaclPedestal()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Commit	Line	Data
21f61e25	1	// **************************************************************************
	2	// * Copyright(c) 1998-1999, 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	// **************************************************************************
	15
	16	#include "AliHMPIDCalib.h" //class header
	17	#include "TTreePlayer.h"
	18	#include <fstream>
	19	#include <TTree.h>
	20	//#include "AliHMPIDDigit.h"
	21
	22	ClassImp(AliHMPIDCalib)
	23
	24
	25	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	26	AliHMPIDCalib::AliHMPIDCalib():
	27	fPedTree(0x0),
	28	fa(-1),
	29	fd(-1),
	30	fr(-1),
	31	fq(-1),
	32	fPedHisto(0)
	33	{
	34	Init();
	35	}
	36	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	37	AliHMPIDCalib::~AliHMPIDCalib()
	38	{
	39	//destructor
	40	}
	41	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	42	void AliHMPIDCalib::Init()
	43	{
	44	for(Int_t iddl=0;iddl<11;iddl++) faddl[iddl]=kFALSE;
	45	//
	46	// Called to initialize the TTree in which the raw data words will be stored
	47	//
	48	fPedHisto=new TH1F("fPedHisto","Temporary pedestal",4096,-0.5,4095.5); //init pedestal histo
	49	fPedTree = new TTree("PedTree","HMPID Pedestal Tree"); //init pedestal tree
	50	fPedTree->Branch("diladd",&fa,"diladd/I");
	51	fPedTree->Branch("dilnum",&fd,"dilnum/I");
	52	fPedTree->Branch("dilrow",&fr,"dilrow/I");
	53	fPedTree->Branch("qdc" ,&fq,"qdc/F ");
	54	}
	55	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	56	void AliHMPIDCalib::FillPedestal(Int_t ddl,Int_t row, Int_t dil,Int_t adr,Int_t q)
	57	{
	58	//
	59	//Called from the HMPIDda and fills the pedestal tree
	60	//
	61	fq=q;
	62	fr=row;
	63	fd=dil;
	64	fa=adr;
65	faddl[ddl]=kTRUE;
66	if(fq>-1) fPedTree->Fill();
67	}//FillPedestal()
68	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
69	Bool_t AliHMPIDCalib::CalcPedestal(Int_t ddl, Char_t* name)
70	{
71	//
72	//Calculate pedestal for each pad
73	//
74
75	ofstream out; //to write the pedestal text files
76	Double_t mean,sigma;
77	Int_t inhard;
78	if(faddl[ddl]==kFALSE) return kFALSE; //if ddl is missing no ped file is created (and also for LDC selection). Check with Paolo what he checks for?!
79	out.open(name);
80	for(Int_t row=1; row < 25; row++){
81	for(Int_t dil=1; dil < 11; dil++){
82	for(Int_t adr=0; adr < 48; adr++){
83	fPedHisto->Reset();
84	fPedTree->Draw("fq>>fPedHisto",Form("fr==%d&&fd==%d&&fa==%d",row,dil,adr));
85	mean=fPedHisto->GetMean();
86	sigma=fPedHisto->GetRMS();
87	inhard=((Int_t(mean))<<9)+Int_t(mean+3*sigma);
88	out << Form("%2i %2i %2i %5.2f %5.2f %x\n",row,dil,adr,mean,sigma,inhard);
89	}//adr
90	}//dil
91	}//row
92	out.close(); //write pedestal file
93	return kTRUE;
94	}//CaclPedestal()
95	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++