///////////////////////////////////////////////////////////////////////////////
// //
-// TPC calibration class for parameters which saved per pad //
+// TPC calibration class for parameters which are saved per pad //
+// Each AliTPCCalPad consists of 72 AliTPCCalROC-objects //
// //
///////////////////////////////////////////////////////////////////////////////
#include "TTreeStream.h"
#include "TFile.h"
#include "TKey.h"
+#include <TFormula.h>
+#include <TString.h>
+#include <TObjString.h>
+#include <iostream>
+#include <AliLog.h>
+
+//graphic includes
+#include <TTree.h>
+#include <TH1.h>
+#include <TCanvas.h>
+#include <TLegend.h>
+#include <TCut.h>
+#include <TVirtualPad.h>
+
ClassImp(AliTPCCalPad)
TObject::Copy(c);
}
+
+void AliTPCCalPad::SetCalROC(AliTPCCalROC* roc, Int_t sector){
+ //
+ // Set AliTPCCalROC copies values from 'roc'
+ // if sector == -1 the sector specified in 'roc' is used
+ // else sector specified in 'roc' is ignored and specified sector is filled
+ //
+ if (sector == -1) sector = roc->GetSector();
+ if (!fROC[sector]) fROC[sector] = new AliTPCCalROC(sector);
+ for (UInt_t ichannel = 0; ichannel < roc->GetNchannels(); ichannel++)
+ fROC[sector]->SetValue(ichannel, roc->GetValue(ichannel));
+}
+
+
+
//_____________________________________________________________________________
void AliTPCCalPad::Add(Float_t c1)
{
//
- // add constant for all channels of all ROCs
+ // add constant c1 to all channels of all ROCs
//
for (Int_t isec = 0; isec < kNsec; isec++) {
//_____________________________________________________________________________
void AliTPCCalPad::Multiply(Float_t c1)
{
- //
- // multiply constant for all channels of all ROCs
- //
+ //
+ // multiply each channel of all ROCs with c1
+ //
for (Int_t isec = 0; isec < kNsec; isec++) {
if (fROC[isec]){
fROC[isec]->Multiply(c1);
//_____________________________________________________________________________
void AliTPCCalPad::Add(const AliTPCCalPad * pad, Double_t c1)
{
- //
- // add calpad channel by channel multiplied by c1 - all ROCs
- //
+ //
+ // multiply AliTPCCalPad 'pad' by c1 and add each channel to the coresponing channel in all ROCs
+ // - pad by pad -
+ //
for (Int_t isec = 0; isec < kNsec; isec++) {
- if (fROC[isec]){
+ if (fROC[isec] && pad->GetCalROC(isec)){
fROC[isec]->Add(pad->GetCalROC(isec),c1);
}
}
//_____________________________________________________________________________
void AliTPCCalPad::Multiply(const AliTPCCalPad * pad)
{
- //
- // multiply calpad channel by channel - all ROCs
- //
- for (Int_t isec = 0; isec < kNsec; isec++) {
+ //
+ // multiply each channel of all ROCs with the coresponding channel of 'pad'
+ // - pad by pad -
+ //
+ for (Int_t isec = 0; isec < kNsec; isec++) {
if (fROC[isec]){
fROC[isec]->Multiply(pad->GetCalROC(isec));
}
//_____________________________________________________________________________
void AliTPCCalPad::Divide(const AliTPCCalPad * pad)
{
- //
- // divide calpad channel by channel - all ROCs
- //
+ //
+ // divide each channel of all ROCs by the coresponding channel of 'pad'
+ // - pad by pad -
+ //
for (Int_t isec = 0; isec < kNsec; isec++) {
if (fROC[isec]){
fROC[isec]->Divide(pad->GetCalROC(isec));
// type=1 - mean
// 2 - median
// 3 - LTM
+ //
Int_t npoints = 0;
for (Int_t i=0;i<72;i++) if (fROC[i]) npoints++;
TGraph * graph = new TGraph(npoints);
Double_t AliTPCCalPad::GetMeanRMS(Double_t &rms)
{
//
- // Calculate mean an RMS of all rocs
+ // Calculates mean and RMS of all ROCs
//
Double_t sum = 0, sum2 = 0, n=0, val=0;
for (Int_t isec = 0; isec < kNsec; isec++) {
}
//_____________________________________________________________________________
-TH1F * AliTPCCalPad::MakeHisto1D(Float_t min, Float_t max,Int_t type){
+TH1F * AliTPCCalPad::MakeHisto1D(Float_t min, Float_t max,Int_t type, Int_t side){
//
// make 1D histo
// type -1 = user defined range
// 0 = nsigma cut nsigma=min
+ //
if (type>=0){
if (type==0){
// nsigma range
max = mean+sigma;
}
}
- char name[1000];
- sprintf(name,"%s Pad 1D",GetTitle());
- TH1F * his = new TH1F(name,name,100, min,max);
+ TString name=Form("%s Pad 1D",GetTitle());
+ TH1F * his = new TH1F(name.Data(),name.Data(),100, min,max);
for (Int_t isec = 0; isec < kNsec; isec++) {
+ if (side==1 && isec%36>18) continue;
+ if (side==-1 && isec%36<18) continue;
if (fROC[isec]){
for (UInt_t irow=0; irow<fROC[isec]->GetNrows(); irow++){
UInt_t npads = (Int_t)fROC[isec]->GetNPads(irow);
// Make 2D graph
// side - specify the side A = 0 C = 1
// type - used types of determination of boundaries in z
+ //
Float_t kEpsilon = 0.000000000001;
TH2F * his = new TH2F(GetName(), GetName(), 250,-250,250,250,-250,250);
AliTPCROC * roc = AliTPCROC::Instance();
}
+AliTPCCalPad* AliTPCCalPad::LocalFit(const char* padName, Int_t rowRadius, Int_t padRadius, AliTPCCalPad* PadOutliers, Bool_t robust, Double_t chi2Threshold, Double_t robustFraction, Bool_t printCurrentSector) const {
+ //
+ // Loops over all AliTPCCalROCs and performs a localFit in each ROC
+ // AliTPCCalPad with fit-data is returned
+ // rowRadius and padRadius specifies a window around a given pad in one sector.
+ // The data of this window are fitted with a parabolic function.
+ // This function is evaluated at the pad's position.
+ // At the edges the window is shifted, so that the specified pad is not anymore in the center of the window.
+ // rowRadius - radius - rows to be used for smoothing
+ // padradius - radius - pads to be used for smoothing
+ // ROCoutlier - map of outliers - pads not to be used for local smoothing
+ // robust - robust method of fitting - (much slower)
+ // chi2Threshold: Threshold for chi2 when EvalRobust is called
+ // robustFraction: Fraction of data that will be used in EvalRobust
+ //
+ //
+ AliTPCCalPad* pad = new AliTPCCalPad(padName, padName);
+ for (Int_t isec = 0; isec < 72; isec++){
+ if (printCurrentSector) std::cout << "LocalFit in sector " << isec << "\r" << std::flush;
+ if (PadOutliers)
+ pad->SetCalROC(GetCalROC(isec)->LocalFit(rowRadius, padRadius, PadOutliers->GetCalROC(isec), robust, chi2Threshold, robustFraction));
+ else
+ pad->SetCalROC(GetCalROC(isec)->LocalFit(rowRadius, padRadius, 0, robust, chi2Threshold, robustFraction));
+ }
+ return pad;
+}
+
+AliTPCCalPad* AliTPCCalPad::GlobalFit(const char* padName, AliTPCCalPad* PadOutliers, Bool_t robust, Int_t fitType, Double_t chi2Threshold, Double_t robustFraction, Double_t err, TObjArray *fitParArr, TObjArray *fitCovArr){
+ //
+ // Loops over all AliTPCCalROCs and performs a globalFit in each ROC
+ // AliTPCCalPad with fit-data is returned
+ // chi2Threshold: Threshold for chi2 when EvalRobust is called
+ // robustFraction: Fraction of data that will be used in EvalRobust
+ // chi2Threshold: Threshold for chi2 when EvalRobust is called
+ // robustFraction: Fraction of data that will be used in EvalRobust
+ // err: error of the data points
+ // if fitParArr and/or fitCovArr is given, write fitParameters and/or covariance Matrices into the array
+ //
+ AliTPCCalPad* pad = new AliTPCCalPad(padName, padName);
+ TVectorD fitParam(0);
+ TMatrixD covMatrix(0,0);
+ Float_t chi2 = 0;
+ for (Int_t isec = 0; isec < 72; isec++){
+ if (PadOutliers)
+ GetCalROC(isec)->GlobalFit(PadOutliers->GetCalROC(isec), robust, fitParam, covMatrix, chi2, fitType, chi2Threshold, robustFraction, err);
+ else
+ GetCalROC(isec)->GlobalFit(0, robust, fitParam, covMatrix, chi2, fitType, chi2Threshold, robustFraction, err);
+
+ AliTPCCalROC *roc=AliTPCCalROC::CreateGlobalFitCalROC(fitParam, isec);
+ pad->SetCalROC(roc);
+ delete roc;
+ if ( fitParArr ) fitParArr->AddAtAndExpand(new TVectorD(fitParam), isec);
+ if ( fitCovArr ) fitCovArr->AddAtAndExpand(new TMatrixD(covMatrix), isec);
+ }
+ return pad;
+}
+//_____________________________________________________________________________
+TObjArray* AliTPCCalPad::CreateFormulaArray(const char *fitFormula)
+{
+ //
+ // create an array of TFormulas for the each parameter of the fit function
+ //
-void AliTPCCalPad::MakeTree(const char * fileName, TObjArray * array, const char * mapFileName, AliTPCCalPad* outlierPad, Float_t ltmFraction) {
- //
- // Write tree with all available information
- //
- AliTPCROC* tpcROCinstance = AliTPCROC::Instance();
-
- TObjArray* mapIROCs = 0;
- TObjArray* mapOROCs = 0;
- TVectorF *mapIROCArray = 0;
- TVectorF *mapOROCArray = 0;
- Int_t mapEntries = 0;
- TString* mapNames = 0;
-
- if (mapFileName) {
- TFile mapFile(mapFileName, "read");
-
- TList* listOfROCs = mapFile.GetListOfKeys();
- mapEntries = listOfROCs->GetEntries()/2;
- mapIROCs = new TObjArray(mapEntries*2);
- mapOROCs = new TObjArray(mapEntries*2);
- mapIROCArray = new TVectorF[mapEntries];
- mapOROCArray = new TVectorF[mapEntries];
-
- mapNames = new TString[mapEntries];
- for (Int_t ivalue = 0; ivalue < mapEntries; ivalue++) {
- TString ROCname(((TKey*)(listOfROCs->At(ivalue*2)))->GetName());
- ROCname.Remove(ROCname.Length()-4, 4);
- mapIROCs->AddAt((AliTPCCalROC*)mapFile.Get((ROCname + "IROC").Data()), ivalue);
- mapOROCs->AddAt((AliTPCCalROC*)mapFile.Get((ROCname + "OROC").Data()), ivalue);
- mapNames[ivalue].Append(ROCname);
- }
-
- for (Int_t ivalue = 0; ivalue < mapEntries; ivalue++) {
- mapIROCArray[ivalue].ResizeTo(tpcROCinstance->GetNChannels(0));
- mapOROCArray[ivalue].ResizeTo(tpcROCinstance->GetNChannels(36));
-
- for (UInt_t ichannel = 0; ichannel < tpcROCinstance->GetNChannels(0); ichannel++)
- (mapIROCArray[ivalue])[ichannel] = ((AliTPCCalROC*)(mapIROCs->At(ivalue)))->GetValue(ichannel);
- for (UInt_t ichannel = 0; ichannel < tpcROCinstance->GetNChannels(36); ichannel++)
- (mapOROCArray[ivalue])[ichannel] = ((AliTPCCalROC*)(mapOROCs->At(ivalue)))->GetValue(ichannel);
- }
+ // split fit string in single parameters
+ // find dimension of the fit:
+ TString fitString(fitFormula);
+ fitString.ReplaceAll("++","#");
+ fitString.ReplaceAll(" ","");
+ TObjArray *arrFitParams = fitString.Tokenize("#");
+ Int_t ndim = arrFitParams->GetEntries();
+ //create array of TFormulas to evaluate the parameters
+ TObjArray *arrFitFormulas = new TObjArray(ndim);
+ arrFitFormulas->SetOwner(kTRUE);
+ for (Int_t idim=0;idim<ndim;++idim){
+ TString s=((TObjString*)arrFitParams->At(idim))->GetString();
+ s.ReplaceAll("gx","[0]");
+ s.ReplaceAll("gy","[1]");
+ s.ReplaceAll("lx","[2]");
+ s.ReplaceAll("ly","[3]");
+ s.ReplaceAll("sector","[4]");
+ arrFitFormulas->AddAt(new TFormula(Form("param%02d",idim),s.Data()),idim);
+ }
+ delete arrFitParams;
+
+ return arrFitFormulas;
+}
+//_____________________________________________________________________________
+void AliTPCCalPad::EvalFormulaArray(const TObjArray &arrFitFormulas, TVectorD &results,
+ const Int_t sec, const Int_t row, const Int_t pad)
+{
+ //
+ // evaluate the fit formulas
+ //
+ Int_t ndim=arrFitFormulas.GetEntries();
+ results.ResizeTo(ndim);
+
+ AliTPCROC* tpcROCinstance = AliTPCROC::Instance(); // to calculate the pad's position
+ Float_t localXYZ[3];
+ Float_t globalXYZ[3];
+ tpcROCinstance->GetPositionLocal(sec, row, pad, localXYZ);
+ tpcROCinstance->GetPositionGlobal(sec, row, pad, globalXYZ);
+ //calculate parameter values
+ for (Int_t idim=0;idim<ndim;++idim){
+ TFormula *f=(TFormula*)arrFitFormulas.At(idim);
+ f->SetParameters(globalXYZ[0],globalXYZ[1],localXYZ[0],localXYZ[1],sec);
+ results[idim]=f->Eval(0);
+ }
+}
+//_____________________________________________________________________________
+void AliTPCCalPad::GlobalSidesFit(const AliTPCCalPad* PadOutliers, const char* fitFormula, TVectorD &fitParamSideA, TVectorD &fitParamSideC,TMatrixD &covMatrixSideA, TMatrixD &covMatrixSideC, Float_t & chi2SideA, Float_t & chi2SideC, AliTPCCalPad *pointError, Bool_t robust, Double_t robustFraction){
+ //
+ // Performs a fit on both sides.
+ // Valid information for the fitFormula are the variables
+ // - gx, gy, lx ,ly: meaning global x, global y, local x, local y value of the padName
+ // - sector: the sector number.
+ // eg. a formula might look 'gy' or '(sector<36) ++ gy' or 'gx ++ gy' or 'gx ++ gy ++ lx ++ lx^2' and so on
+ //
+ // PadOutliers - pads with value !=0 are not used in fitting procedure
+ // chi2Threshold: Threshold for chi2 when EvalRobust is called
+ // robustFraction: Fraction of data that will be used in EvalRobust
+ //
- } // if (mapFileName)
+ TObjArray* arrFitFormulas=CreateFormulaArray(fitFormula);
+ Int_t ndim = arrFitFormulas->GetEntries();
+ //resize output data arrays
+ fitParamSideA.ResizeTo(ndim+1);
+ fitParamSideC.ResizeTo(ndim+1);
+ covMatrixSideA.ResizeTo(ndim+1,ndim+1);
+ covMatrixSideC.ResizeTo(ndim+1,ndim+1);
+ // create linear fitter for A- and C- Side
+ TLinearFitter* fitterGA = new TLinearFitter(ndim+1,Form("hyp%d",ndim));
+ TLinearFitter* fitterGC = new TLinearFitter(ndim+1,Form("hyp%d",ndim));
+ fitterGA->StoreData(kTRUE);
+ fitterGC->StoreData(kTRUE);
+ //parameter values
+ TVectorD parValues(ndim);
+
+ AliTPCCalROC *rocErr=0x0;
- TTreeSRedirector cstream(fileName);
- Int_t arrayEntries = array->GetEntries();
-
- TString* names = new TString[arrayEntries];
- for (Int_t ivalue = 0; ivalue < arrayEntries; ivalue++)
- names[ivalue].Append(((AliTPCCalPad*)array->At(ivalue))->GetName());
-
- for (UInt_t isector = 0; isector < tpcROCinstance->GetNSectors(); isector++) {
- //
- // get statistic for given sector
- //
- TVectorF median(arrayEntries);
- TVectorF mean(arrayEntries);
- TVectorF rms(arrayEntries);
- TVectorF ltm(arrayEntries);
- TVectorF ltmrms(arrayEntries);
- TVectorF medianWithOut(arrayEntries);
- TVectorF meanWithOut(arrayEntries);
- TVectorF rmsWithOut(arrayEntries);
- TVectorF ltmWithOut(arrayEntries);
- TVectorF ltmrmsWithOut(arrayEntries);
-
- TVectorF *vectorArray = new TVectorF[arrayEntries];
- for (Int_t ivalue = 0; ivalue < arrayEntries; ivalue++)
- vectorArray[ivalue].ResizeTo(tpcROCinstance->GetNChannels(isector));
-
- for (Int_t ivalue = 0; ivalue < arrayEntries; ivalue++) {
- AliTPCCalPad* calPad = (AliTPCCalPad*) array->At(ivalue);
- AliTPCCalROC* calROC = calPad->GetCalROC(isector);
- AliTPCCalROC* outlierROC = 0;
- if (outlierPad) outlierROC = outlierPad->GetCalROC(isector);
- if (calROC) {
- median[ivalue] = calROC->GetMedian();
- mean[ivalue] = calROC->GetMean();
- rms[ivalue] = calROC->GetRMS();
- Double_t ltmrmsValue = 0;
- ltm[ivalue] = calROC->GetLTM(<mrmsValue, ltmFraction);
- ltmrms[ivalue] = ltmrmsValue;
- if (outlierROC) {
- medianWithOut[ivalue] = calROC->GetMedian(outlierROC);
- meanWithOut[ivalue] = calROC->GetMean(outlierROC);
- rmsWithOut[ivalue] = calROC->GetRMS(outlierROC);
- ltmrmsValue = 0;
- ltmWithOut[ivalue] = calROC->GetLTM(<mrmsValue, ltmFraction, outlierROC);
- ltmrmsWithOut[ivalue] = ltmrmsValue;
- }
- }
- else {
- median[ivalue] = 0.;
- mean[ivalue] = 0.;
- rms[ivalue] = 0.;
- ltm[ivalue] = 0.;
- ltmrms[ivalue] = 0.;
- medianWithOut[ivalue] = 0.;
- meanWithOut[ivalue] = 0.;
- rmsWithOut[ivalue] = 0.;
- ltmWithOut[ivalue] = 0.;
- ltmrmsWithOut[ivalue] = 0.;
- }
- }
-
- //
- // fill vectors of variable per pad
- //
- TVectorF *posArray = new TVectorF[8];
- for (Int_t ivalue = 0; ivalue < 8; ivalue++)
- posArray[ivalue].ResizeTo(tpcROCinstance->GetNChannels(isector));
-
- Float_t posG[3] = {0};
- Float_t posL[3] = {0};
- Int_t ichannel = 0;
- for (UInt_t irow = 0; irow < tpcROCinstance->GetNRows(isector); irow++) {
- for (UInt_t ipad = 0; ipad < tpcROCinstance->GetNPads(isector, irow); ipad++) {
- tpcROCinstance->GetPositionLocal(isector, irow, ipad, posL);
- tpcROCinstance->GetPositionGlobal(isector, irow, ipad, posG);
- posArray[0][ichannel] = irow;
- posArray[1][ichannel] = ipad;
- posArray[2][ichannel] = posL[0];
- posArray[3][ichannel] = posL[1];
- posArray[4][ichannel] = posG[0];
- posArray[5][ichannel] = posG[1];
- posArray[6][ichannel] = (Int_t)(ipad - (Double_t)(tpcROCinstance->GetNPads(isector, irow))/2);
- posArray[7][ichannel] = ichannel;
-
- // loop over array containing AliTPCCalPads
- for (Int_t ivalue = 0; ivalue < arrayEntries; ivalue++) {
- AliTPCCalPad* calPad = (AliTPCCalPad*) array->At(ivalue);
- AliTPCCalROC* calROC = calPad->GetCalROC(isector);
- if (calROC)
- (vectorArray[ivalue])[ichannel] = calROC->GetValue(irow, ipad);
- else
- (vectorArray[ivalue])[ichannel] = 0;
- }
- ichannel++;
- }
+ for (UInt_t isec = 0; isec<kNsec; ++isec){
+ AliTPCCalROC *rocOut=PadOutliers->GetCalROC(isec);
+ AliTPCCalROC *rocData=GetCalROC(isec);
+ if (pointError) rocErr=pointError->GetCalROC(isec);
+ if (!rocData) continue;
+ for (UInt_t irow = 0; irow < GetCalROC(isec)->GetNrows(); irow++) {
+ for (UInt_t ipad = 0; ipad < GetCalROC(isec)->GetNPads(irow); ipad++) {
+ //check for outliers
+ if (rocOut && rocOut->GetValue(irow,ipad)) continue;
+ //calculate parameter values
+ EvalFormulaArray(*arrFitFormulas,parValues,isec,irow,ipad);
+ //get value
+ Float_t value=rocData->GetValue(irow,ipad);
+ //point error
+ Int_t err=1;
+ if (rocErr) {
+ err=TMath::Nint(rocErr->GetValue(irow,ipad));
+ if (err==0) err=1;
+ }
+ //add points to the fitters
+ if (isec/18%2==0){
+ fitterGA->AddPoint(parValues.GetMatrixArray(),value,err);
+ }else{
+ fitterGC->AddPoint(parValues.GetMatrixArray(),value,err);
+ }
}
-
- cstream << "calPads" <<
- "sector=" << isector;
-
- for (Int_t ivalue = 0; ivalue < arrayEntries; ivalue++) {
- cstream << "calPads" <<
- (Char_t*)((names[ivalue] + "_Median=").Data()) << median[ivalue] <<
- (Char_t*)((names[ivalue] + "_Mean=").Data()) << mean[ivalue] <<
- (Char_t*)((names[ivalue] + "_RMS=").Data()) << rms[ivalue] <<
- (Char_t*)((names[ivalue] + "_LTM=").Data()) << ltm[ivalue] <<
- (Char_t*)((names[ivalue] + "_RMS_LTM=").Data()) << ltmrms[ivalue];
- if (outlierPad) {
- cstream << "calPads" <<
- (Char_t*)((names[ivalue] + "_Median_OutlierCutted=").Data()) << medianWithOut[ivalue] <<
- (Char_t*)((names[ivalue] + "_Mean_OutlierCutted=").Data()) << meanWithOut[ivalue] <<
- (Char_t*)((names[ivalue] + "_RMS_OutlierCutted=").Data()) << rmsWithOut[ivalue] <<
- (Char_t*)((names[ivalue] + "_LTM_OutlierCutted=").Data()) << ltmWithOut[ivalue] <<
- (Char_t*)((names[ivalue] + "_RMS_LTM_OutlierCutted=").Data()) << ltmrmsWithOut[ivalue];
- }
+ }
+ }
+ if (robust){
+ fitterGA->EvalRobust(robustFraction);
+ fitterGC->EvalRobust(robustFraction);
+ } else {
+ fitterGA->Eval();
+ fitterGC->Eval();
+ }
+ chi2SideA=fitterGA->GetChisquare()/(fitterGA->GetNpoints()-(ndim+1));
+ chi2SideC=fitterGC->GetChisquare()/(fitterGC->GetNpoints()-(ndim+1));
+ fitterGA->GetParameters(fitParamSideA);
+ fitterGC->GetParameters(fitParamSideC);
+ fitterGA->GetCovarianceMatrix(covMatrixSideA);
+ fitterGC->GetCovarianceMatrix(covMatrixSideC);
+
+ delete arrFitFormulas;
+ delete fitterGA;
+ delete fitterGC;
+
+}
+//
+AliTPCCalPad *AliTPCCalPad::CreateCalPadFit(const char* fitFormula, const TVectorD &fitParamSideA, const TVectorD &fitParamSideC)
+{
+ //
+ //
+ //
+ TObjArray *arrFitFormulas=CreateFormulaArray(fitFormula);
+ Int_t ndim = arrFitFormulas->GetEntries();
+ //check if dimension of fit formula and fit parameters agree
+ if (ndim!=fitParamSideA.GetNrows()||ndim!=fitParamSideC.GetNrows()){
+ printf("AliTPCCalPad::CreateCalPadFit: Dimensions of fit formula and fit Parameters does not match!");
+ return 0;
+ }
+ //create cal pad
+ AliTPCCalPad *pad=new AliTPCCalPad("fitResultPad",Form("Fit result: %s",fitFormula));
+ //fill cal pad with fit results if requested
+ for (UInt_t isec = 0; isec<kNsec; ++isec){
+ AliTPCCalROC *roc=pad->GetCalROC(isec);
+ for (UInt_t irow = 0; irow < roc->GetNrows(); irow++) {
+ for (UInt_t ipad = 0; ipad < roc->GetNPads(irow); ipad++) {
+ const TVectorD *fitPar=0;
+ TVectorD fitResArray;
+ if (isec/18%2==0){
+ fitPar=&fitParamSideA;
+ }else{
+ fitPar=&fitParamSideC;
+ }
+ EvalFormulaArray(*arrFitFormulas,fitResArray, isec, irow, ipad);
+ for (Int_t idim=0;idim<ndim;++idim)
+ fitResArray(idim)*=(*fitPar)(idim);
+ roc->SetValue(irow,ipad,fitResArray.Sum());
}
+ }
+ }
+ delete arrFitFormulas;
+ return pad;
+}
- for (Int_t ivalue = 0; ivalue < arrayEntries; ivalue++) {
- cstream << "calPads" <<
- (Char_t*)((names[ivalue] + ".=").Data()) << &vectorArray[ivalue];
- }
- if (mapFileName) {
- for (Int_t ivalue = 0; ivalue < mapEntries; ivalue++) {
- if (isector < 36)
- cstream << "calPads" <<
- (Char_t*)((mapNames[ivalue] + ".=").Data()) << &mapIROCArray[ivalue];
- else
- cstream << "calPads" <<
- (Char_t*)((mapNames[ivalue] + ".=").Data()) << &mapOROCArray[ivalue];
- }
- }
- cstream << "calPads" <<
- "row.=" << &posArray[0] <<
- "pad.=" << &posArray[1] <<
- "lx.=" << &posArray[2] <<
- "ly.=" << &posArray[3] <<
- "gx.=" << &posArray[4] <<
- "gy.=" << &posArray[5] <<
- "rpad.=" << &posArray[6] <<
- "channel.=" << &posArray[7];
-
- cstream << "calPads" <<
- "\n";
-
- delete[] posArray;
- delete[] vectorArray;
- }
-
- delete[] names;
- if (mapFileName) {
- delete mapIROCs;
- delete mapOROCs;
- delete[] mapIROCArray;
- delete[] mapOROCArray;
- delete[] mapNames;
- }
+TCanvas * AliTPCCalPad::MakeReportPadSector(TTree *chain, const char* varName, const char*varTitle, const char *axisTitle, Float_t min, Float_t max, const char *cutUser){
+ //
+ // Make a report - cal pads per sector
+ // mean valeus per sector and local X
+ //
+ TH1* his=0;
+ TLegend *legend = 0;
+ TCanvas *canvas = new TCanvas(Form("Sector: %s",varTitle),Form("Sector: %s",varTitle),1500,1100);
+
+ canvas->Divide(2);
+ chain->SetAlias("lX","lx.fElements");
+ //
+ canvas->cd(1);
+ TString strDraw=varName;
+ strDraw+=":lX";
+ legend = new TLegend(0.5,0.50,0.9,0.9, Form("%s TPC A side", varTitle));
+ for (Int_t isec=-1; isec<18; isec+=1){
+ TCut cutSec=Form("sector%%36==%d",isec);
+ cutSec+=cutUser;
+ if (isec==-1) cutSec="sector%36<18";
+ chain->SetMarkerColor(1+(isec+2)%5);
+ chain->SetLineColor(1+(isec+2)%5);
+ chain->SetMarkerStyle(25+(isec+2)%4);
+ //
+ chain->Draw(strDraw.Data(),cutSec,"profgoff");
+ his=(TH1*)chain->GetHistogram()->Clone();
+ delete chain->GetHistogram();
+ his->SetMaximum(max);
+ his->SetMinimum(min);
+ his->GetXaxis()->SetTitle("R (cm)");
+ his->GetYaxis()->SetTitle(axisTitle);
+ his->SetTitle(Form("%s- sector %d",varTitle, isec));
+ his->SetName(Form("%s- sector %d",varTitle, isec));
+ if (isec==-1) his->SetTitle(Form("%s A side",varTitle));
+ if (isec==-1) his->Draw();
+ his->Draw("same");
+ legend->AddEntry(his);
+ }
+ legend->Draw();
+ canvas->cd(2);
+ //
+ legend = new TLegend(0.5,0.50,0.9,0.9, Form("%s TPC C side", varTitle));
+ for (Int_t isec=-1; isec<18; isec+=1){
+ TCut cutSec=Form("(sector+18)%%36==%d",isec);
+ cutSec+=cutUser;
+ if (isec==-1) cutSec="sector%36>18";
+ chain->SetMarkerColor(1+(isec+2)%5);
+ chain->SetLineColor(1+(isec+2)%5);
+ chain->SetMarkerStyle(25+isec%4);
+ //
+ chain->Draw(strDraw.Data(),cutSec,"profgoff");
+ his=(TH1*)chain->GetHistogram()->Clone();
+ delete chain->GetHistogram();
+ his->SetMaximum(max);
+ his->SetMinimum(min);
+ his->GetXaxis()->SetTitle("R (cm)");
+ his->GetYaxis()->SetTitle(axisTitle);
+ his->SetTitle(Form("%s- sector %d",varTitle,isec));
+ his->SetName(Form("%s- sector %d",varTitle,isec));
+ if (isec==-1) his->SetTitle(Form("%s C side",varTitle));
+ if (isec==-1) his->Draw();
+ his->Draw("same");
+ legend->AddEntry(his);
+ }
+ legend->Draw();
+ //
+ //
+ return canvas;
}
+TCanvas * AliTPCCalPad::MakeReportPadSector2D(TTree *chain, const char* varName, const char*varTitle, const char *axisTitle, Float_t min, Float_t max, const char *cutUser){
+ //
+ // Make a report - cal pads per sector
+ // 2D view
+ // Input tree should be created using AliPreprocesorOnline before
+ //
+ TH1* his=0;
+ TCanvas *canvas = new TCanvas(Form("%s2D",varTitle),Form("%s2D",varTitle),1500,1100);
+ canvas->Divide(2);
+ //
+ TString strDraw=varName;
+ strDraw+=":gy.fElements:gx.fElements>>his(250,-250,250,250,-250,250)";
+ //
+ TVirtualPad * pad=0;
+ pad=canvas->cd(1);
+ pad->SetMargin(0.15,0.15,0.15,0.15);
+ TCut cut=cutUser;
+ chain->Draw(strDraw.Data(),"sector%36<18"+cut,"profgoffcolz2");
+ his=(TH1*)chain->GetHistogram()->Clone();
+ delete chain->GetHistogram();
+ his->SetMaximum(max);
+ his->SetMinimum(min);
+ his->GetXaxis()->SetTitle("x (cm)");
+ his->GetYaxis()->SetTitle("y (cm)");
+ his->GetZaxis()->SetTitle(axisTitle);
+ his->SetTitle(Form("%s A side",varTitle));
+ his->SetName(Form("%s A side",varTitle));
+ his->Draw("colz2");
+ //
+ pad=canvas->cd(2);
+ pad->SetMargin(0.15,0.15,0.15,0.15);
+
+ chain->Draw(strDraw.Data(),"sector%36>=18"+cut,"profgoffcolz2");
+ his=(TH1*)chain->GetHistogram()->Clone();
+ delete chain->GetHistogram();
+ his->SetMaximum(max);
+ his->SetMinimum(min);
+ his->GetXaxis()->SetTitle("x (cm)");
+ his->GetYaxis()->SetTitle("y (cm)");
+ his->GetZaxis()->SetTitle(axisTitle);
+ his->SetTitle(Form("%s C side",varTitle));
+ his->SetName(Form("%s C side",varTitle));
+ his->Draw("colz2");
+ //
+ //
+ return canvas;
+}
+
+void AliTPCCalPad::Draw(Option_t* option){
+ //
+ // Draw function - standard 2D view
+ //
+ TH1* his=0;
+ TCanvas *canvas = new TCanvas(Form("%s2D",GetTitle()),Form("%s2D",GetTitle()),900,900);
+ canvas->Divide(2,2);
+ //
+ //
+ TVirtualPad * pad=0;
+ pad=canvas->cd(1);
+ pad->SetMargin(0.15,0.15,0.15,0.15);
+ his=MakeHisto2D(0);
+ his->GetXaxis()->SetTitle("x (cm)");
+ his->GetYaxis()->SetTitle("y (cm)");
+ his->GetZaxis()->SetTitle(GetTitle());
+ his->SetTitle(Form("%s A side",GetTitle()));
+ his->SetName(Form("%s A side",GetTitle()));
+ his->Draw(option);
+ //
+ pad=canvas->cd(2);
+ pad->SetMargin(0.15,0.15,0.15,0.15);
+ his=MakeHisto2D(1);
+ his->GetXaxis()->SetTitle("x (cm)");
+ his->GetYaxis()->SetTitle("y (cm)");
+ his->GetZaxis()->SetTitle(GetTitle());
+ his->SetTitle(Form("%s C side",GetTitle()));
+ his->SetName(Form("%s C side",GetTitle()));
+ his->Draw(option);
+ //
+ pad=canvas->cd(3);
+ pad->SetMargin(0.15,0.15,0.15,0.15);
+ his=MakeHisto1D(-8,8,0,1);
+ his->GetXaxis()->SetTitle(GetTitle());
+ his->SetTitle(Form("%s A side",GetTitle()));
+ his->SetName(Form("%s A side",GetTitle()));
+ his->Draw("err");
+ //
+ pad=canvas->cd(4);
+ pad->SetMargin(0.15,0.15,0.15,0.15);
+ his=MakeHisto1D(-8,8,0,-1);
+ his->GetXaxis()->SetTitle(GetTitle());
+ his->SetTitle(Form("%s C side",GetTitle()));
+ his->SetName(Form("%s C side",GetTitle()));
+ his->Draw("err");
+
+
+}
+
+
+AliTPCCalPad * AliTPCCalPad::MakeCalPadFromHistoRPHI(TH2 * hisA, TH2* hisC){
+ //
+ // Make cal pad from r-phi histograms
+ //
+ AliTPCROC *proc= AliTPCROC::Instance();
+ AliTPCCalPad *calPad = new AliTPCCalPad("his","his");
+ Float_t globalPos[3];
+ for (Int_t isec=0; isec<72; isec++){
+ AliTPCCalROC* calRoc = calPad->GetCalROC(isec);
+ TH2 * his = ((isec%36<18) ? hisA:hisC);
+ for (UInt_t irow=0; irow<calRoc->GetNrows(); irow+=1){
+ Int_t jrow=irow;
+ if (isec>=36) jrow+=63;
+ for (UInt_t ipad=0;ipad<proc->GetNPads(isec,irow);ipad+=1){
+ proc->GetPositionGlobal(isec,irow,ipad, globalPos);
+ Double_t phi=TMath::ATan2(globalPos[1],globalPos[0]);
+ //if (phi<0) phi+=TMath::Pi()*2;
+ Int_t bin=his->FindBin(phi,jrow);
+ Float_t value= his->GetBinContent(bin);
+ calRoc->SetValue(irow,ipad,value);
+ }
+ }
+ }
+ return calPad;
+}