/* $Id$ */
-/*
- example: fill pedestal with gausschen noise
- AliTRDCalibPadStatus ped;
- ped.TestEvent(numberofevent);
- // Method without histo
- //////////////////////////
- ped.Analyse();
- //Create the histo of the AliTRDCalROC
- TH2F * histo2dm = ped.GetCalRocMean(0,kFALSE)->MakeHisto2D();
- histo2dm->Scale(10.0);
- TH1F * histo1dm = ped.GetCalRocMean(0,kFALSE)->MakeHisto1D();
- histo1dm->Scale(10.0);
- TH2F * histo2ds = ped.GetCalRocSquares(0,kFALSE)->MakeHisto2D();
- histo2ds->Scale(10.0);
- TH1F * histo1ds = ped.GetCalRocSquares(0,kFALSE)->MakeHisto1D();
- histo1ds->Scale(10.0)
- //Draw output;
- TCanvas* c1 = new TCanvas;
- c1->Divide(2,2);
- c1->cd(1);
- histo2dm->Draw("colz");
- c1->cd(2);
- histo1dm->Draw();
- c1->cd(3);
- histo2ds->Draw("colz");
- c1->cd(4);
- histo1ds->Draw();
-// Method with histo
-/////////////////////////
-ped.AnalyseHisto();
-//Take the histo
-TH1F * histo = ped.GetHisto(31);
-histo->SetEntries(1);
-histo->Draw();
-
-*/
+////////////////////////////////////////////////////////////////////////////
+// //
+// Example: fill pedestal with Gaussian noise //
+// //
+// AliTRDCalibPadStatus ped; //
+// ped.TestEvent(numberofevent); //
+// //
+// // Method without histo //
+// ped.Analyse(); //
+// //
+// // Create the histo of the AliTRDCalROC //
+// TH2F * histo2dm = ped.GetCalRocMean(0,kFALSE)->MakeHisto2D(); //
+// histo2dm->Scale(10.0); //
+// TH1F * histo1dm = ped.GetCalRocMean(0,kFALSE)->MakeHisto1D(); //
+// histo1dm->Scale(10.0); //
+// TH2F * histo2ds = ped.GetCalRocSquares(0,kFALSE)->MakeHisto2D(); //
+// histo2ds->Scale(10.0); //
+// TH1F * histo1ds = ped.GetCalRocSquares(0,kFALSE)->MakeHisto1D(); //
+// histo1ds->Scale(10.0) //
+// //
+// // Draw output //
+// TCanvas* c1 = new TCanvas; //
+// c1->Divide(2,2); //
+// c1->cd(1); //
+// histo2dm->Draw("colz"); //
+// c1->cd(2); //
+// histo1dm->Draw(); //
+// c1->cd(3); //
+// histo2ds->Draw("colz"); //
+// c1->cd(4); //
+// histo1ds->Draw(); //
+// //
+// // Method with histo //
+// ped.AnalyseHisto(); //
+// //
+// // Take the histo //
+// TH1F *histo = ped.GetHisto(31); //
+// histo->SetEntries(1); //
+// histo->Draw(); //
+//
+// Authors:
+// R. Bailhache (R.Bailhache@gsi.de, rbailhache@ikf.uni-frankfurt.de)
+// J. Book (jbook@ikf.uni-frankfurt.de)
+// //
+////////////////////////////////////////////////////////////////////////////
+
//Root includes
#include <TObjArray.h>
-#include <TH1F.h>
#include <TH2F.h>
#include <TString.h>
#include <TMath.h>
-#include <TF1.h>
#include <TRandom.h>
+
+//#include <TRandom.h>
#include <TDirectory.h>
#include <TFile.h>
-#include <AliMathBase.h>
-#include "TTreeStream.h"
+
//AliRoot includes
+#include <AliMathBase.h>
#include "AliRawReader.h"
#include "AliRawReaderRoot.h"
#include "AliRawReaderDate.h"
-#include "AliTRDRawStream.h"
+
+//header file
+#include "AliLog.h"
+#include "AliTRDCalibPadStatus.h"
+#include "AliTRDgeometry.h"
+#include "AliTRDCommonParam.h"
#include "./Cal/AliTRDCalROC.h"
#include "./Cal/AliTRDCalPadStatus.h"
+#include "./Cal/AliTRDCalDet.h"
+#include "./Cal/AliTRDCalPad.h"
#include "./Cal/AliTRDCalSingleChamberStatus.h"
-#include "AliTRDarrayF.h"
-#include "AliTRDCommonParam.h"
+#include "AliTRDdigitsManager.h"
+#include "AliTRDdigitsParam.h"
+#include "AliTRDSignalIndex.h"
+#include "AliTRDarraySignal.h"
+#include "AliTRDarrayADC.h"
+#include "AliTRDfeeParam.h"
+
+#include "AliTRDrawStream.h"
#ifdef ALI_DATE
#include "event.h"
#endif
-//header file
-#include "AliTRDCalibPadStatus.h"
-
-
-
ClassImp(AliTRDCalibPadStatus) /*FOLD00*/
//_____________________________________________________________________
AliTRDCalibPadStatus::AliTRDCalibPadStatus() : /*FOLD00*/
TObject(),
+ fGeo(0),
fAdcMin(0),
- fAdcMax(20),
+ fAdcMax(21),
fDetector(-1),
fNumberOfTimeBins(0),
- fCalArrayEntries(540),
- fCalArrayMean(540),
- fCalArraySquares(540),
fCalRocArrayMean(540),
fCalRocArrayRMS(540),
- fHistoArray(540),
- fCalEntries(0x0),
- fCalMean(0x0),
- fCalSquares(0x0)
+ fCalRocArrayMeand(540),
+ fCalRocArrayRMSd(540),
+ fHistoArray(540)
{
//
// default constructor
//
+
+ fGeo = new AliTRDgeometry();
+
}
//_____________________________________________________________________
AliTRDCalibPadStatus::AliTRDCalibPadStatus(const AliTRDCalibPadStatus &ped) : /*FOLD00*/
TObject(ped),
+ fGeo(0),
fAdcMin(ped.GetAdcMin()),
fAdcMax(ped.GetAdcMax()),
fDetector(ped.fDetector),
fNumberOfTimeBins(ped.fNumberOfTimeBins),
- fCalArrayEntries(540),
- fCalArrayMean(540),
- fCalArraySquares(540),
fCalRocArrayMean(540),
fCalRocArrayRMS(540),
- fHistoArray(540),
- fCalEntries(0x0),
- fCalMean(0x0),
- fCalSquares(0x0)
+ fCalRocArrayMeand(540),
+ fCalRocArrayRMSd(540),
+ fHistoArray(540)
{
//
// copy constructor
//
for (Int_t idet = 0; idet < 540; idet++){
- const AliTRDarrayF *calEntries = (AliTRDarrayF*)ped.fCalArrayEntries.UncheckedAt(idet);
- const AliTRDarrayF *calMean = (AliTRDarrayF*)ped.fCalArrayMean.UncheckedAt(idet);
- const AliTRDarrayF *calSquares = (AliTRDarrayF*)ped.fCalArraySquares.UncheckedAt(idet);
const AliTRDCalROC *calRocMean = (AliTRDCalROC*)ped.fCalRocArrayMean.UncheckedAt(idet);
const AliTRDCalROC *calRocRMS = (AliTRDCalROC*)ped.fCalRocArrayRMS.UncheckedAt(idet);
+ const AliTRDCalROC *calRocMeand = (AliTRDCalROC*)ped.fCalRocArrayMeand.UncheckedAt(idet);
+ const AliTRDCalROC *calRocRMSd = (AliTRDCalROC*)ped.fCalRocArrayRMSd.UncheckedAt(idet);
const TH2F *hped = (TH2F*)ped.fHistoArray.UncheckedAt(idet);
- if ( calEntries != 0x0 ) fCalArrayEntries.AddAt(new AliTRDarrayF(*calEntries), idet);
- if ( calMean != 0x0 ) fCalArrayMean.AddAt(new AliTRDarrayF(*calMean), idet);
- if ( calSquares != 0x0 ) fCalArraySquares.AddAt(new AliTRDarrayF(*calSquares), idet);
if ( calRocMean != 0x0 ) fCalRocArrayMean.AddAt(new AliTRDCalROC(*calRocMean), idet);
if ( calRocRMS != 0x0 ) fCalRocArrayRMS.AddAt(new AliTRDCalROC(*calRocRMS), idet);
+ if ( calRocMeand != 0x0 ) fCalRocArrayMeand.AddAt(new AliTRDCalROC(*calRocMeand), idet);
+ if ( calRocRMSd != 0x0 ) fCalRocArrayRMSd.AddAt(new AliTRDCalROC(*calRocRMSd), idet);
+
if ( hped != 0x0 ){
TH2F *hNew = new TH2F(*hped);
hNew->SetDirectory(0);
}
}
+ if (fGeo) {
+ delete fGeo;
+ }
+ fGeo = new AliTRDgeometry();
}
+
//_____________________________________________________________________
AliTRDCalibPadStatus& AliTRDCalibPadStatus::operator = (const AliTRDCalibPadStatus &source)
{
//
// destructor
//
+ fCalRocArrayMean.Delete();
+ fCalRocArrayRMS.Delete();
+ fCalRocArrayMeand.Delete();
+ fCalRocArrayRMSd.Delete();
+ fHistoArray.Delete();
+ if (fGeo) {
+ delete fGeo;
+ }
}
//_____________________________________________________________________
-Int_t AliTRDCalibPadStatus::Update(const Int_t icdet, /*FOLD00*/
- const Int_t icRow,
- const Int_t icCol,
- const Int_t csignal,
- const Int_t crowMax)
+void AliTRDCalibPadStatus::Destroy()
{
- //
- // Signal filling methode
- //
- if ( (csignal>fAdcMax) || (csignal<fAdcMin) ) return 0;
-
- if(fDetector != icdet){
- fCalEntries = ((AliTRDarrayF *)GetCalEntries(icdet,kTRUE));
- fCalMean = ((AliTRDarrayF *)GetCalMean(icdet,kTRUE));
- fCalSquares = ((AliTRDarrayF *)GetCalSquares(icdet,kTRUE));
- }
-
- Float_t entries = fCalEntries->At(icRow+icCol*crowMax);
- Float_t mean = fCalMean->At(icRow+icCol*crowMax);
- Float_t squares = fCalSquares->At(icRow+icCol*crowMax);
-
- Float_t entriesn = entries+1.0;
- fCalEntries->AddAt(entriesn,(icRow+icCol*crowMax));
- Float_t meann = (mean*entries+((Float_t)(csignal+0.5)))/entriesn;
- fCalMean->AddAt(meann,icRow+icCol*crowMax);
- Float_t squaresn = ((squares*entries)+(((Float_t)(csignal+0.5))*((Float_t)(csignal+0.5))))/entriesn;
- fCalSquares->AddAt(squaresn,icRow+icCol*crowMax);
-
- //printf("icdet %d, icRow %d, icCol %d, csignal %d, crowMax %d\n",icdet,icRow,icCol,csignal,crowMax);
- //printf("entries %f, mean %f, squares %f\n",entriesn,meann,squaresn);
-
- fDetector = icdet;
-
- return 0;
+ //
+ // Destroy
+ //
+ fCalRocArrayMean.Delete();
+ fCalRocArrayRMS.Delete();
+ fCalRocArrayMeand.Delete();
+ fCalRocArrayRMSd.Delete();
+ fHistoArray.Delete();
}
//_____________________________________________________________________
Int_t AliTRDCalibPadStatus::UpdateHisto(const Int_t icdet, /*FOLD00*/
- const Int_t icRow,
- const Int_t icCol,
- const Int_t csignal,
- const Int_t crowMax)
+ const Int_t icRow,
+ const Int_t icCol,
+ const Int_t csignal,
+ const Int_t crowMax,
+ const Int_t ccold,
+ const Int_t icMcm,
+ const Int_t icRob
+ )
{
- //
- // Signal filling methode
- //
+ //
+ // Signal filling methode
+ //
Int_t nbchannel = icRow+icCol*crowMax;
+ Int_t mCm = icMcm%4;
+ Int_t rOb = icRob%2;
+
+ // now the case of double read channel
+ if(ccold > 0){
+ nbchannel = (((ccold-1)*8+ (mCm+rOb*4))*crowMax+icRow)+144*crowMax;
+ //printf("nbchannel %d, ccold %d, icMcm %d, crowMax %d, icRow %d\n",nbchannel,ccold,icMcm,crowMax,icRow);
+ }
// fast filling methode.
// Attention: the entry counter of the histogram is not increased
// this means that e.g. the colz draw option gives an empty plot
+
Int_t bin = 0;
- if ( !(((Int_t)csignal>fAdcMax ) || ((Int_t)csignal<fAdcMin)) )
- bin = (nbchannel+1)*(fAdcMax-fAdcMin+2)+((Int_t)csignal-fAdcMin+1);
+ if ( !(((Int_t)csignal>=fAdcMax ) || ((Int_t)csignal<fAdcMin)) )
+ bin = (nbchannel+1)*(fAdcMax-fAdcMin+2)+((Int_t)csignal-fAdcMin+1);
+
+ //GetHisto(icdet,kTRUE)->Fill(csignal,nbchannel);
+
GetHisto(icdet,kTRUE)->GetArray()[bin]++;
return 0;
}
//_____________________________________________________________________
-Bool_t AliTRDCalibPadStatus::ProcessEvent(AliTRDRawStream *rawStream, Bool_t nocheck)
+Int_t AliTRDCalibPadStatus::ProcessEvent(AliRawReader *rawReader)
{
//
- // Event Processing loop - AliTRDRawStream
+ // RawReader = AliTRDrawStream (Jochen Klein)
//
-
-
- Bool_t withInput = kFALSE;
-
- if(!nocheck) {
- while (rawStream->Next()) {
- Int_t rawversion = rawStream->GetRawVersion(); // current raw version
- if(rawversion != 2) return kFALSE;
- Int_t idetector = rawStream->GetDet(); // current detector
- Int_t iRow = rawStream->GetRow(); // current row
- Int_t iRowMax = rawStream->GetMaxRow(); // current rowmax
- Int_t iCol = rawStream->GetCol(); // current col
- Int_t iTimeBin = rawStream->GetTimeBin(); // current time bin
- Int_t *signal = rawStream->GetSignals(); // current ADC signal
- Int_t nbtimebin = rawStream->GetNumberOfTimeBins(); // number of time bins read from data
-
- if((fDetector != -1) && (nbtimebin != fNumberOfTimeBins)) return kFALSE;
- fNumberOfTimeBins = nbtimebin;
-
- Int_t fin = TMath::Min(nbtimebin,(iTimeBin+3));
- Int_t n = 0;
-
- for(Int_t k = iTimeBin; k < fin; k++){
- if(signal[n]>0) UpdateHisto(idetector,iRow,iCol,signal[n],iRowMax);
- n++;
- }
-
- withInput = kTRUE;
- }
- }
- else {
- while (rawStream->Next()) {
- Int_t idetector = rawStream->GetDet(); // current detector
- Int_t iRow = rawStream->GetRow(); // current row
- Int_t iRowMax = rawStream->GetMaxRow(); // current rowmax
- Int_t iCol = rawStream->GetCol(); // current col
- Int_t iTimeBin = rawStream->GetTimeBin(); // current time bin
- Int_t *signal = rawStream->GetSignals(); // current ADC signal
- Int_t nbtimebin = rawStream->GetNumberOfTimeBins(); // number of time bins read from data
-
- Int_t fin = TMath::Min(nbtimebin,(iTimeBin+3));
- Int_t n = 0;
-
- for(Int_t k = iTimeBin; k < fin; k++){
- if(signal[n]>0) UpdateHisto(idetector,iRow,iCol,signal[n],iRowMax);
- n++;
- }
-
- withInput = kTRUE;
- }
- }
-
- return withInput;
-}
-//_____________________________________________________________________
-Bool_t AliTRDCalibPadStatus::ProcessEvent(AliRawReader *rawReader, Bool_t nocheck)
-{
- //
- // Event processing loop - AliRawReader
+ // Event Processing loop - AliTRDRawStreamCosmic
+ // 0 time bin problem or zero suppression
+ // 1 no input
+ // 2 input
+ // Raw version number:
+ // [3,31] non zero suppressed
+ // 2,4 and [32,63] zero suppressed
//
+
+
+ Int_t withInput = 1;
- AliTRDRawStream rawStream(rawReader);
+ AliTRDdigitsManager *digitsManager = new AliTRDdigitsManager(kTRUE);
+ digitsManager->CreateArrays();
- rawReader->Select("TRD");
+ AliTRDrawStream *rawStream = new AliTRDrawStream(rawReader);
+ rawStream->SetDigitsManager(digitsManager);
+ //rawStream->SetNoErrorWarning();
+ //rawStream->SetSharedPadReadout(kTRUE);
+
+ AliTRDfeeParam *feeParam = AliTRDfeeParam::Instance();
- return ProcessEvent(&rawStream, nocheck);
-}
-//_________________________________________________________________________
-Bool_t AliTRDCalibPadStatus::ProcessEvent(
-#ifdef ALI_DATE
- eventHeaderStruct *event,
- Bool_t nocheck
-#else
- eventHeaderStruct* /*event*/,
- Bool_t /*nocheck*/
+ Int_t det = 0;
+ while ((det = rawStream->NextChamber(digitsManager, NULL, NULL)) >= 0) { //idetector
+ if (digitsManager->GetIndexes(det)->HasEntry()) {//QA
+ // printf("there is ADC data on this chamber!\n");
+
+ AliTRDarrayADC *digits = (AliTRDarrayADC *) digitsManager->GetDigits(det); //mod
+ if (digits->HasData()) { //array
+
+ AliTRDSignalIndex *indexes = digitsManager->GetIndexes(det);
+ if (indexes->IsAllocated() == kFALSE) {
+ AliError("Indexes do not exist!");
+ break;
+ }
+ Int_t iRow = 0;
+ Int_t iCol = 0;
+ indexes->ResetCounters();
+
+ while (indexes->NextRCIndex(iRow, iCol)) { //column,row
+
+ AliTRDdigitsParam *digitParam = (AliTRDdigitsParam *)digitsManager->GetDigitsParam();
+
+ Int_t mcm = 0; // MCM from AliTRDfeeParam
+ Int_t rob = 0; // ROB from AliTRDfeeParam
+ Int_t extCol = 0; // extended column from AliTRDfeeParam
+ mcm = feeParam->GetMCMfromPad(iRow,iCol);
+ rob = feeParam->GetROBfromPad(iRow,iCol);
+
+ Int_t idetector = det; // current detector
+ Int_t iRowMax = 16; // current rowmax
+ if(GetStack(det) == 2) iRowMax = 12;
+
+ Int_t adc = 20 - (iCol%18) -1; // current adc
+ Int_t col = 0; // col!=0 ->Shared Pad
+ extCol = feeParam->GetExtendedPadColFromADC(rob,mcm,adc);
+ //printf(" iCol %d iRow %d iRowMax %d rob %d mcm %d adc %d extCol %d\n",iCol,iRow,iRowMax,rob,mcm,adc,extCol);
+
+ // Signal for regular pads
+ Int_t nbtimebin = digitParam->GetNTimeBins(idetector); // number of time bins read from data
+ for(Int_t k = 0; k < nbtimebin; k++){
+ Short_t signal = 0;
+ signal = digits->GetData(iRow,iCol,k);
+
+ if(signal>0) {
+ UpdateHisto(idetector,iRow,iCol,signal,iRowMax,col,mcm,rob);
+ }
+ }
+
+
+
+ if((adc==3-1 || adc==20-1 || adc==19-1) && (iCol > 1 && iCol <142) ) { //SHARED PADS
-#endif
- )
-{
- //
- // process date event
- //
-#ifdef ALI_DATE
- AliRawReader *rawReader = new AliRawReaderDate((void*)event);
- Bool_t result=ProcessEvent(rawReader, nocheck);
- delete rawReader;
- return result;
-#else
- Fatal("AliTRDCalibPadStatus", "this class was compiled without DATE");
- return 0;
-#endif
-
-}
-//_____________________________________________________________________
-Bool_t AliTRDCalibPadStatus::TestEvent(Int_t nevent) /*FOLD00*/
-{
- //
- // Test event loop
- // fill one oroc and one iroc with random gaus
- //
-
- AliTRDCommonParam *comParam = AliTRDCommonParam::Instance();
- if (!comParam) {
- return 0;
- }
- gRandom->SetSeed(0);
-
- for (Int_t ism=0; ism<18; ism++){
- for (Int_t ich=0; ich < 5; ich++){
- for (Int_t ipl=0; ipl < 6; ipl++){
- for(Int_t irow = 0; irow < comParam->GetRowMax(ipl,ich,ism); irow++){
- for(Int_t icol = 0; icol < comParam->GetColMax(ipl); icol++){
- for (Int_t iTimeBin=0; iTimeBin<(30*nevent); iTimeBin++){
- Int_t signal=(Int_t)(ich+8+gRandom->Gaus(0,1.2));
- if ( signal>0 )Update((ipl+ich*6+ism*6*5),irow,icol,signal,comParam->GetRowMax(ipl,ich,ism));
- }
- }
- }
+ switch(adc) {
+ case 2:
+ adc = 20; //shared Pad adc
+ mcm = feeParam->GetMCMfromSharedPad(iRow,iCol); //shared Pad mcm
+ col = 1;
+ break;
+ case 19:
+ adc = 1; //shared Pad adc
+ mcm = feeParam->GetMCMfromSharedPad(iRow,iCol); //shared Pad mcm
+ col = 2;
+ break;
+ case 18:
+ adc = 0; //shared Pad adc
+ mcm = feeParam->GetMCMfromSharedPad(iRow,iCol); //shared Pad mcm
+ col = 3;
+ break;
}
- }
- }
- return kTRUE;
-}
-//_____________________________________________________________________
-Bool_t AliTRDCalibPadStatus::TestEventHisto(Int_t nevent) /*FOLD00*/
-{
- //
- // Test event loop
- // fill one oroc and one iroc with random gaus
- //
-
- AliTRDCommonParam *comParam = AliTRDCommonParam::Instance();
- if (!comParam) {
- return 0;
- }
- gRandom->SetSeed(0);
-
- for (Int_t ism=0; ism<18; ism++){
- for (Int_t ich=0; ich < 5; ich++){
- for (Int_t ipl=0; ipl < 6; ipl++){
- for(Int_t irow = 0; irow < comParam->GetRowMax(ipl,ich,ism); irow++){
- for(Int_t icol = 0; icol < comParam->GetColMax(ipl); icol++){
- for (Int_t iTimeBin=0; iTimeBin<(30*nevent); iTimeBin++){
- Int_t signal=(Int_t)(ich+8+gRandom->Gaus(0,1.2));
- if ( signal>0 )UpdateHisto((ipl+ich*6+ism*6*5),irow,icol,signal,comParam->GetRowMax(ipl,ich,ism));
- }
- }
+ rob = feeParam->GetROBfromSharedPad(iRow,iCol); //shared Pad rob
+
+
+ extCol = feeParam->GetExtendedPadColFromADC(rob,mcm,adc); //extended pad col via the shared pad rob,mcm and adc
+
+ //printf("SHARED PAD --- iCol %d iRow %d rob %d mcm %d adc %d extCol %d col %d\n",iCol,iRow,rob,mcm,adc,extCol,col);
+ for(Int_t k = 0; k < nbtimebin; k++){
+ Short_t signal = 0;
+ signal = digits->GetDataByAdcCol(iRow,extCol,k);
+
+ if(signal>0) {
+ UpdateHisto(idetector,iRow,iCol,signal,iRowMax,col,mcm,rob);
}
}
- }
- }
- return kTRUE;
+ } //shared pads end
+
+
+ withInput = 2;
+ }//column,row
+
+ }//array
+ }//QA
+ digitsManager->ClearArrays(det);
+ }//idetector
+ delete digitsManager;
+ delete rawStream;
+ return withInput;
+
}
+
//_____________________________________________________________________
TH2F* AliTRDCalibPadStatus::GetHisto(Int_t det, TObjArray *arr, /*FOLD00*/
Int_t nbinsY, Float_t ymin, Float_t ymax,
- Char_t *type, Bool_t force)
+ const Char_t *type, Bool_t force)
{
//
// return pointer to histogram
// if we are forced and histogram doesn't yes exist create it
Char_t name[255], title[255];
- sprintf(name,"hCalib%s%.2d",type,det);
- sprintf(title,"%s calibration histogram detector %.2d;ADC channel;Channel (pad)",type,det);
+ snprintf(name,255,"hCalib%s%.3d",type,det);
+ snprintf(title,255,"%s calibration histogram detector %.2d;ADC channel;Channel (pad)",type,det);
+
+
+ Int_t nbchannels = fGeo->GetRowMax(GetLayer(det),GetStack(det),GetSector(det))*fGeo->GetColMax(GetLayer(det));
+
+ // we will add 3*8*rowMax channels at the end for the double counted
+ nbchannels += 3*8*(fGeo->GetRowMax(GetLayer(det),GetStack(det),GetSector(det)));
- AliTRDCommonParam *comParam = AliTRDCommonParam::Instance();
- if (!comParam) {
- return 0x0;
- }
- Int_t nbchannels = comParam->GetRowMax(GetPlane(det),GetChamber(det),GetSector(det))*comParam->GetColMax(GetPlane(det));
// new histogram with calib information. One value for each pad!
TH2F* hist = new TH2F(name,title,
arr->AddAt(hist,det);
return hist;
}
+
//_____________________________________________________________________
TH2F* AliTRDCalibPadStatus::GetHisto(Int_t det, Bool_t force) /*FOLD00*/
{
// if force is true create a new histogram if it doesn't exist allready
//
TObjArray *arr = &fHistoArray;
- return GetHisto(det, arr, fAdcMax-fAdcMin, fAdcMin, fAdcMax, "Pedestal", force);
+ return GetHisto(det, arr, fAdcMax-fAdcMin, fAdcMin-0.5, fAdcMax-0.5, "Pedestal", force);
}
-//_____________________________________________________________________
-AliTRDarrayF* AliTRDCalibPadStatus::GetCal(Int_t det, TObjArray* arr, Bool_t force) /*FOLD00*/
-{
- //
- // return pointer to ROC Calibration
- // if force is true create a new AliTRDarrayF if it doesn't exist allready
- //
- if ( !force || arr->UncheckedAt(det) )
- return (AliTRDarrayF*)arr->UncheckedAt(det);
-
- // if we are forced and histogram doesn't yes exist create it
- AliTRDarrayF *croc = new AliTRDarrayF();
- AliTRDCommonParam *comParam = AliTRDCommonParam::Instance();
- if (!comParam) {
- return croc;
- }
- Int_t nbpad = comParam->GetRowMax(GetPlane(det),GetChamber(det),GetSector(det))*comParam->GetColMax(GetPlane(det));
- // new AliTRDCalROC. One value for each pad!
- croc->Expand(nbpad);
- for(Int_t k = 0; k < nbpad; k++){
- croc->AddAt(0.0,k);
- }
- arr->AddAt(croc,det);
- return croc;
-}
//_____________________________________________________________________
AliTRDCalROC* AliTRDCalibPadStatus::GetCalRoc(Int_t det, TObjArray* arr, Bool_t force) /*FOLD00*/
{
// if we are forced and histogram doesn't yes exist create it
// new AliTRDCalROC. One value for each pad!
- AliTRDCalROC *croc = new AliTRDCalROC(GetPlane(det),GetChamber(det));
+ AliTRDCalROC *croc = new AliTRDCalROC(GetLayer(det),GetStack(det));
arr->AddAt(croc,det);
return croc;
}
//_____________________________________________________________________
-AliTRDarrayF* AliTRDCalibPadStatus::GetCalEntries(Int_t det, Bool_t force) /*FOLD00*/
-{
- //
- // return pointer to Carge ROC Calibration
- // if force is true create a new histogram if it doesn't exist allready
- //
- TObjArray *arr = &fCalArrayEntries;
- return GetCal(det, arr, force);
-}
-//_____________________________________________________________________
-AliTRDarrayF* AliTRDCalibPadStatus::GetCalMean(Int_t det, Bool_t force) /*FOLD00*/
+AliTRDCalROC* AliTRDCalibPadStatus::GetCalRocMean(Int_t det, Bool_t force) /*FOLD00*/
{
//
// return pointer to Carge ROC Calibration
// if force is true create a new histogram if it doesn't exist allready
//
- TObjArray *arr = &fCalArrayMean;
- return GetCal(det, arr, force);
+ TObjArray *arr = &fCalRocArrayMean;
+ return GetCalRoc(det, arr, force);
}
+
//_____________________________________________________________________
-AliTRDarrayF* AliTRDCalibPadStatus::GetCalSquares(Int_t det, Bool_t force) /*FOLD00*/
+AliTRDCalROC* AliTRDCalibPadStatus::GetCalRocRMS(Int_t det, Bool_t force) /*FOLD00*/
{
//
// return pointer to Carge ROC Calibration
// if force is true create a new histogram if it doesn't exist allready
//
- TObjArray *arr = &fCalArraySquares;
- return GetCal(det, arr, force);
+ TObjArray *arr = &fCalRocArrayRMS;
+ return GetCalRoc(det, arr, force);
}
//_____________________________________________________________________
-AliTRDCalROC* AliTRDCalibPadStatus::GetCalRocMean(Int_t det, Bool_t force) /*FOLD00*/
+AliTRDCalROC* AliTRDCalibPadStatus::GetCalRocMeand(Int_t det, Bool_t force) /*FOLD00*/
{
//
// return pointer to Carge ROC Calibration
// if force is true create a new histogram if it doesn't exist allready
//
- TObjArray *arr = &fCalRocArrayMean;
+ TObjArray *arr = &fCalRocArrayMeand;
return GetCalRoc(det, arr, force);
}
+
//_____________________________________________________________________
-AliTRDCalROC* AliTRDCalibPadStatus::GetCalRocRMS(Int_t det, Bool_t force) /*FOLD00*/
+AliTRDCalROC* AliTRDCalibPadStatus::GetCalRocRMSd(Int_t det, Bool_t force) /*FOLD00*/
{
//
// return pointer to Carge ROC Calibration
// if force is true create a new histogram if it doesn't exist allready
//
- TObjArray *arr = &fCalRocArrayRMS;
+ TObjArray *arr = &fCalRocArrayRMSd;
return GetCalRoc(det, arr, force);
}
-//_________________________________________________________________________
-void AliTRDCalibPadStatus::Analyse() /*FOLD00*/
-{
- //
- // Calcul the rms properly
- //
-
- for(Int_t idet = 0; idet < 540; idet++){
- // Take the stuff
- fCalEntries = ((AliTRDarrayF *)GetCalEntries(idet));
- fCalMean = ((AliTRDarrayF *)GetCalMean(idet));
- fCalSquares = ((AliTRDarrayF *)GetCalSquares(idet));
-
- if(!fCalEntries) continue;
-
- AliTRDCalROC *calRocMean = ((AliTRDCalROC *)GetCalRocMean(idet,kTRUE));
- AliTRDCalROC *calRocRMS = ((AliTRDCalROC *)GetCalRocRMS(idet,kTRUE));
-
- // range channels
- Int_t channels = calRocMean->GetNchannels();
-
- for(Int_t ichannels = 0 ; ichannels < channels; ichannels++){
-
- Float_t entries = fCalEntries->At(ichannels);
- Float_t mean = fCalMean->At(ichannels);
- Float_t squares = fCalSquares->At(ichannels);
-
- Float_t rms = 0.0;
- if(entries > 0){
- Double_t rm = TMath::Abs(squares-(mean*mean));
- rms = TMath::Sqrt(rm);
- calRocRMS->SetValue(ichannels,rms/10.0);
- calRocMean->SetValue(ichannels,mean/10.0);
- }
- }
- }
-}
//_____________________________________________________________________
void AliTRDCalibPadStatus::AnalyseHisto() /*FOLD00*/
{
Int_t nbinsAdc = fAdcMax-fAdcMin;
- TVectorD param(3);
+ TVectorD param(4);
TMatrixD dummy(3,3);
- Float_t *array_hP=0;
+ Float_t *arrayHP=0;
for (Int_t idet=0; idet<540; idet++){
continue;
}
+ //printf("Entries for %d\n",idet);
+
AliTRDCalROC *rocMean = GetCalRocMean(idet,kTRUE);
AliTRDCalROC *rocRMS = GetCalRocRMS(idet,kTRUE);
- array_hP = hP->GetArray();
+ arrayHP = hP->GetArray();
Int_t nChannels = rocMean->GetNchannels();
for (Int_t iChannel=0; iChannel<nChannels; iChannel++){
Int_t offset = (nbinsAdc+2)*(iChannel+1)+1;
- Double_t ret = AliMathBase::FitGaus(array_hP+offset,nbinsAdc,fAdcMin,fAdcMax,¶m,&dummy);
- // if the fitting failed set noise and pedestal to 0
- if ((ret == -4) || (ret == -1)) {
+ Double_t ret = AliMathBase::FitGaus(arrayHP+offset,nbinsAdc,fAdcMin-0.5,fAdcMax-0.5,¶m,&dummy);
+ // if the fitting failed set noise and pedestal to 0
+ if ((ret==-4) || (ret==-1) || (ret==-2)) {
param[1]=0.0;
param[2]=0.0;
}
rocMean->SetValue(iChannel,param[1]/10.0);
rocRMS->SetValue(iChannel,param[2]/10.0);
}
+
+ // here we analyse doubled read channels
+
+ AliTRDCalROC *rocMeand = GetCalRocMeand(idet,kTRUE);
+ AliTRDCalROC *rocRMSd = GetCalRocRMSd(idet,kTRUE);
+
+ Int_t nrows = rocMeand->GetNrows();
+ Int_t shift = 144*nrows;
+ Int_t total = shift+3*8*nrows;
+
+ for (Int_t iChannel=shift; iChannel<total; iChannel++){
+ Int_t offset = (nbinsAdc+2)*(iChannel+1)+1;
+ Double_t ret = AliMathBase::FitGaus(arrayHP+offset,nbinsAdc,fAdcMin-0.5,fAdcMax-0.5,¶m,&dummy);
+ // if the fitting failed set noise and pedestal to 0
+ if ((ret==-4) || (ret==-1) || (ret==-2)) {
+ param[1]=0.0;
+ param[2]=0.0;
+ }
+ if((param[1]/10.0) > 65534.0) param[1] = 0.0;
+ if((param[2]/10.0) > 65534.0) param[2] = 0.0;
+
+ // here we have to recalculate backward
+ Int_t nb = iChannel-shift;
+ Int_t row = nb%nrows;
+ Int_t j = (Int_t)(nb/nrows);
+ Int_t imcm = j%8;
+ Int_t icol = (Int_t)(j/8);
+
+ Int_t finalcol = 18*imcm;
+ if(icol > 0) icol += 17;
+ else icol = -1;
+ finalcol += icol;
+
+ Int_t channel = row+finalcol*nrows;
+
+ //printf("iChannel %d, nrows %d, finalcol %d, row %d, channel %d\n",iChannel,nrows,finalcol,row,channel);
+ if((finalcol < 0) || (finalcol >= 144)) continue;
+
+ rocMeand->SetValue(channel,param[1]/10.0);
+ rocRMSd->SetValue(channel,param[2]/10.0);
+ }
+
}
-
+
}
+
//_______________________________________________________________________________________
AliTRDCalPadStatus* AliTRDCalibPadStatus::CreateCalPadStatus()
{
//
- //
+ // Create Pad Status out of Mean and RMS values
+ // The chamber without data are masked, this is the corrected in the preprocessor
//
AliTRDCalPadStatus* obj = new AliTRDCalPadStatus("padstatus", "padstatus");
{
AliTRDCalSingleChamberStatus *calROC = obj->GetCalROC(idet);
- //Take the stuff
- fCalEntries = ((AliTRDarrayF *)GetCalEntries(idet));
- AliTRDCalROC *calRocMean = ((AliTRDCalROC *)GetCalRocMean(idet));
- AliTRDCalROC *calRocRMS = ((AliTRDCalROC *)GetCalRocRMS(idet));
- if ( !calRocMean ) {
+ if ( !GetCalRocMean(idet)) {
for(Int_t k = 0; k < calROC->GetNchannels(); k++){
- calROC->SetStatus(k,AliTRDCalPadStatus::kMasked);
+ calROC->SetStatus(k,AliTRDCalPadStatus::kNotConnected);
}
- //printf("no fCalRocMean for %d\n",idet);
continue;
}
+
- //Range
- Int_t channels = calROC->GetNchannels();
-
-
- //Mean
- Float_t meanentries = 0.0;
- if(!fCalEntries){
- if(GetHisto(idet)){
- meanentries = GetHisto(idet)->GetEntries()/(channels);
- }
- }
- else meanentries = TMath::Mean(channels,((TArrayF *)fCalEntries)->GetArray());
- //Double_t meanmean = calRocMean->GetMean()*10.0;
- //Double_t meansquares = calRocRMS->GetMean()*10.0;
+ //Take the stuff
+ AliTRDCalROC *calRocMean = new AliTRDCalROC(*( (AliTRDCalROC *) GetCalRocMean(idet)));
+ AliTRDCalROC *calRocRMS = new AliTRDCalROC(*( (AliTRDCalROC *) GetCalRocRMS(idet)));
+
+ //Take the stuff second chance
+ AliTRDCalROC *calRocMeand = new AliTRDCalROC(*( (AliTRDCalROC *) GetCalRocMeand(idet)));
+ AliTRDCalROC *calRocRMSd = new AliTRDCalROC(*( (AliTRDCalROC *) GetCalRocRMSd(idet)));
+ calRocRMS->Unfold();
+ calRocRMSd->Unfold();
+
+
+ //Range
+ Int_t row = calROC->GetNrows();
+ Int_t col = calROC->GetNcols();
+
+ Double_t rmsmean = calRocMean->GetRMS()*10.0;
+ Double_t meanmean = calRocMean->GetMean()*10.0;
+ Double_t meansquares = calRocRMS->GetMean();
- for(Int_t ich = 0; ich < channels; ich++){
+
+ for(Int_t irow = 0; irow < row; irow++){
+
+ // for bridged pads
+ Float_t meanprevious = 0.0;
+ Float_t rmsprevious = 0.0;
+ Float_t mean = 0.0;
+ Float_t rms = 0.0;
- Float_t entries = 0.0;
- if(!fCalEntries){
- if(GetHisto(idet)){
- for(Int_t bin = 0; bin < (fAdcMax-fAdcMin); bin++){
- entries += GetHisto(idet)->GetArray()[(ich+1)*(fAdcMax-fAdcMin+2)+(bin+1)];
+ for(Int_t icol = 0; icol < col; icol++){
+
+ mean = calRocMean->GetValue(icol,irow)*10.0;
+ rms = calRocRMS->GetValue(icol,irow);
+
+ if(icol > 0) {
+ meanprevious = calRocMean->GetValue((icol -1),irow)*10.0;
+ rmsprevious = calRocRMS->GetValue((icol - 1),irow);
+ }
+
+ Bool_t pb = kFALSE;
+ // masked if two noisy
+ if((rms <= 0.0001) || (TMath::Abs(mean-meanmean)>(5*rmsmean)) || (TMath::Abs(rms)>(5.0*TMath::Abs(meansquares)))) {
+
+ pb = kTRUE;
+ // look at second chance
+ Float_t meand = calRocMeand->GetValue(icol,irow)*10.0;
+ Float_t rmsd = calRocRMSd->GetValue(icol,irow);
+
+ if((rmsd <= 0.0001) || (TMath::Abs(meand-meanmean)>(5*rmsmean)) || (TMath::Abs(rmsd)>(5.0*TMath::Abs(meansquares)))) {
+ if((rmsd <= 0.0001) && (rms <= 0.0001)) {
+ calROC->SetStatus(icol,irow,AliTRDCalPadStatus::kNotConnected);
+ }
+ else {
+ calROC->SetStatus(icol, irow, AliTRDCalPadStatus::kMasked);
+ }
+ }
+ else {
+ calROC->SetStatus(icol, irow, AliTRDCalPadStatus::kReadSecond);
}
}
+
+
+ // bridge if previous pad found something
+ if(!pb) {
+ if((TMath::Abs(meanprevious-mean)<0.000001) && (TMath::Abs(rmsprevious-rms)<0.000001) && (mean > 0.0001)) {
+ //printf("mean previous %f, mean %f, rms %f, rmsprevious %f, col %d\n",meanprevious,mean,rms,rmsprevious,icol);
+ calROC->SetStatus(icol -1 ,irow, AliTRDCalPadStatus::kPadBridgedRight);
+ calROC->SetStatus(icol ,irow, AliTRDCalPadStatus::kPadBridgedLeft);
+ }
+ }
+
}
- else entries = fCalEntries->At(ich);
- //Float_t mean = calRocMean->GetValue(ich)*10.0;
- Float_t rms = calRocRMS->GetValue(ich)*10.0;
- if(ich > 1720) printf("rms %f\n",rms);
-
- //if((entries < 0.3*meanentries) || (TMath::Abs(rms-meansquares) > ((Float_t)(fAdcMax-fAdcMin)/2.0)) || (rms == 0.0)) calROC->SetStatus(ich, AliTRDCalPadStatus::kMasked);
- if(rms <= 0.01) calROC->SetStatus(ich, AliTRDCalPadStatus::kMasked);
- }
+ }
+
+ delete calRocMean;
+ delete calRocRMS;
+ delete calRocMeand;
+ delete calRocRMSd;
+
+
+ }
+
+ return obj;
+
+}
+//_______________________________________________________________________________________
+AliTRDCalPad* AliTRDCalibPadStatus::CreateCalPad()
+{
+ //
+ // Create Pad Noise out of RMS values
+ //
+
+ AliTRDCalPad* obj = new AliTRDCalPad("PadNoise", "PadNoise");
+
+
+ for (Int_t det=0; det<AliTRDgeometry::kNdet; ++det) {
+
+ AliTRDCalROC *calROC22 = obj->GetCalROC(det);
+
+ AliTRDCalROC *calRocRMS = ((AliTRDCalROC *)GetCalRocRMS(det,kTRUE));
+
+ for(Int_t k = 0; k < calROC22->GetNchannels(); k++){
+ calROC22->SetValue(k,calRocRMS->GetValue(k));
}
+
+ }
return obj;
}
+
+//_______________________________________________________________________________________
+AliTRDCalDet* AliTRDCalibPadStatus::CreateCalDet() const
+{
+ //
+ // Create Det Noise correction factor
+ //
+
+ AliTRDCalDet* obj = new AliTRDCalDet("DetNoise", "DetNoise (correction factor)");
+
+ for(Int_t l = 0; l < 540; l++){
+ obj->SetValue(l,10.0);
+ }
+
+ return obj;
+
+}
+
//_____________________________________________________________________
void AliTRDCalibPadStatus::DumpToFile(const Char_t *filename, const Char_t *dir, Bool_t append) /*FOLD00*/
{
f.Close();
if ( backup ) backup->cd();
-}//_____________________________________________________________________________
-Int_t AliTRDCalibPadStatus::GetPlane(Int_t d) const
+}
+
+//_____________________________________________________________________
+void AliTRDCalibPadStatus::SetCalRocMean(AliTRDCalROC *mean, Int_t det) /*FOLD00*/
+{
+ //
+ // Put the AliTRDCalROC in the array fCalRocArrayMean
+ //
+
+
+ AliTRDCalROC *rocMean = GetCalRocMean(det,kTRUE);
+
+ Int_t nChannels = rocMean->GetNchannels();
+
+ for (Int_t iChannel=0; iChannel<nChannels; iChannel++){
+
+ rocMean->SetValue(iChannel,mean->GetValue(iChannel));
+
+ }
+
+}
+
+//_____________________________________________________________________
+void AliTRDCalibPadStatus::SetCalRocRMS(AliTRDCalROC *rms, Int_t det) /*FOLD00*/
+{
+ //
+ // Put the AliTRDCalROC in the array fCalRocArrayRMS
+ //
+
+
+ AliTRDCalROC *rocRms = GetCalRocRMS(det,kTRUE);
+
+ Int_t nChannels = rocRms->GetNchannels();
+
+ for (Int_t iChannel=0; iChannel<nChannels; iChannel++){
+
+ rocRms->SetValue(iChannel,rms->GetValue(iChannel));
+
+ }
+
+}
+//_____________________________________________________________________
+void AliTRDCalibPadStatus::SetCalRocMeand(AliTRDCalROC *mean, Int_t det) /*FOLD00*/
+{
+ //
+ // Put the AliTRDCalROC in the array fCalRocArrayMean
+ //
+
+
+ AliTRDCalROC *rocMean = GetCalRocMeand(det,kTRUE);
+
+ Int_t nChannels = rocMean->GetNchannels();
+
+ for (Int_t iChannel=0; iChannel<nChannels; iChannel++){
+
+ rocMean->SetValue(iChannel,mean->GetValue(iChannel));
+
+ }
+
+}
+
+//_____________________________________________________________________
+void AliTRDCalibPadStatus::SetCalRocRMSd(AliTRDCalROC *rms, Int_t det) /*FOLD00*/
+{
+ //
+ // Put the AliTRDCalROC in the array fCalRocArrayRMS
+ //
+
+
+ AliTRDCalROC *rocRms = GetCalRocRMSd(det,kTRUE);
+
+ Int_t nChannels = rocRms->GetNchannels();
+
+ for (Int_t iChannel=0; iChannel<nChannels; iChannel++){
+
+ rocRms->SetValue(iChannel,rms->GetValue(iChannel));
+
+ }
+
+}
+//_____________________________________________________________________________
+Int_t AliTRDCalibPadStatus::GetLayer(Int_t d) const
{
//
- // Reconstruct the plane number from the detector number
+ // Reconstruct the layer number from the detector number
//
return ((Int_t) (d % 6));
}
//_____________________________________________________________________________
-Int_t AliTRDCalibPadStatus::GetChamber(Int_t d) const
+Int_t AliTRDCalibPadStatus::GetStack(Int_t d) const
{
//
// Reconstruct the chamber number from the detector number
return ((Int_t) (d % 30) / 6);
}
+
//_____________________________________________________________________________
Int_t AliTRDCalibPadStatus::GetSector(Int_t d) const
{
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff