+/**************************************************************************
+ * 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. *
+ **************************************************************************/
+
+/* $Id$ */
+
////////////////////////////////////////////////////////////
// Author: Henrik Tydesjo //
// This class is used in the detector algorithm framework //
// to process the data stored in special container files //
// (see AliITSOnlineSPDscan). For instance, minimum //
-// threshold values can be calculated. //
+// threshold values can be extracted. //
////////////////////////////////////////////////////////////
#include "AliITSOnlineSPDscanAnalyzer.h"
#include "AliITSRawStreamSPD.h"
#include <TStyle.h>
#include <TMath.h>
+#include <TLine.h>
#include <TF1.h>
#include <TGraph.h>
#include <TH2F.h>
+#include <TError.h>
+#include <iostream>
+#include <fstream>
+
+using std::ifstream;
Double_t itsSpdErrorf(Double_t *x, Double_t *par){
if (par[2]<0) par[2]=0;
//Double_t itsSpdErrorfOrig(Double_t *x, Double_t *par){
// return 0.5+0.5*TMath::Erf((x[0]-par[0])/par[1]/sqrt(2.));
//}
+//_________________________________________________________________________
+Double_t itsSpdScurveForMeanTh(Double_t *x, Double_t *par){
+ if (par[2]<0) par[2]=0;
+ Double_t val = 1.- par[2]*(1.-TMath::Erf((x[0]-par[0])/par[1]/sqrt(2.)));
+// Double_t val = par[2]+(0.12*256*32-par[2])*(0.5+0.5*TMath::Erf((x[0]-par[0])/par[1]/sqrt(2.)));
+ return val;
+}
-
-AliITSOnlineSPDscanAnalyzer::AliITSOnlineSPDscanAnalyzer(Char_t *fileName) :
- fType(99),fDacId(99),fScanObj(NULL),fTriggers(NULL),
+//_________________________________________________________________________
+AliITSOnlineSPDscanAnalyzer::AliITSOnlineSPDscanAnalyzer(const Char_t *fileName, AliITSOnlineCalibrationSPDhandler *handler, Bool_t readFromGridFile) :
+ fType(99),fDacId(99),fFileName(fileName),fScanObj(NULL),fHandler(handler),fTriggers(NULL),fTPeff(0),fTPeffHS(NULL),fDeadPixel(0),fDeadPixelHS(NULL),fNoisyPixel(0),fNoisyPixelHS(NULL),
fOverWrite(kFALSE),fNoiseThreshold(0.01),fNoiseMinimumEvents(100),
- fMinNrStepsBeforeIncrease(5),fMinIncreaseFromBaseLine(2),fStepDownDacSafe(2),fMaxBaseLineLevel(10)
+ fMinNrStepsBeforeIncrease(5),fMinIncreaseFromBaseLine(2),fStepDownDacSafe(5),fMaxBaseLineLevel(10)
{
// constructor
- sprintf(fFileName,"%s",fileName);
for (UInt_t chipNr=0; chipNr<11; chipNr++) {
for (UInt_t hs=0; hs<6; hs++) {
fMeanMultiplicity[hs][chipNr]=NULL;
fHitEventEfficiency[hs][chipNr]=NULL;
}
}
- for (Int_t module=0; module<240; module++) {
- fHandler[module]=NULL;
+ for (UInt_t hs=0; hs<6; hs++) {
+ fTPeffChip[hs]=NULL;
+ fDeadPixelChip[hs]=NULL;
+ fNoisyPixelChip[hs]=NULL;
}
- Init();
-}
+ for (UInt_t mod=0; mod<240; mod++) {
+ fbModuleScanned[mod]=kFALSE;
+ }
+ Init(readFromGridFile);
+}
+//_________________________________________________________________________
AliITSOnlineSPDscanAnalyzer::AliITSOnlineSPDscanAnalyzer(const AliITSOnlineSPDscanAnalyzer& handle) :
- fType(99),fDacId(99),fScanObj(NULL),fTriggers(NULL),
+ fType(99),fDacId(99),fFileName("."),fScanObj(NULL),fHandler(NULL),fTriggers(NULL),fTPeff(0),fTPeffHS(NULL),fDeadPixel(0),fDeadPixelHS(NULL),fNoisyPixel(0),fNoisyPixelHS(NULL),
fOverWrite(kFALSE),fNoiseThreshold(0.01),fNoiseMinimumEvents(100),
- fMinNrStepsBeforeIncrease(5),fMinIncreaseFromBaseLine(2),fStepDownDacSafe(2),fMaxBaseLineLevel(10)
+ fMinNrStepsBeforeIncrease(5),fMinIncreaseFromBaseLine(2),fStepDownDacSafe(5),fMaxBaseLineLevel(10)
{
- // copy constructor, only copies the filename (not the processed data)
- sprintf(fFileName,"%s",handle.fFileName);
+ // copy constructor, only copies the filename and params (not the processed data)
+ fFileName=handle.fFileName;
+ fOverWrite=handle.fOverWrite;
+ fNoiseThreshold=handle.fNoiseThreshold;
+ fNoiseMinimumEvents=handle.fNoiseMinimumEvents;
+ fMinNrStepsBeforeIncrease=handle.fMinNrStepsBeforeIncrease;
+ fMinIncreaseFromBaseLine=handle.fMinIncreaseFromBaseLine;
+ fStepDownDacSafe=handle.fStepDownDacSafe;
+ fMaxBaseLineLevel=handle.fMaxBaseLineLevel;
- fScanObj=NULL;
- fType=99;
- fDacId=99;
for (UInt_t chipNr=0; chipNr<11; chipNr++) {
for (UInt_t hs=0; hs<6; hs++) {
fMeanMultiplicity[hs][chipNr]=NULL;
fHitEventEfficiency[hs][chipNr]=NULL;
}
}
- fTriggers=NULL;
- for (Int_t module=0; module<240; module++) {
- fHandler[module]=NULL;
+ for (UInt_t hs=0; hs<6; hs++) {
+ fTPeffChip[hs]=NULL;
+ fDeadPixelChip[hs]=NULL;
+ fNoisyPixelChip[hs]=NULL;
+ }
+
+ for (UInt_t mod=0; mod<240; mod++) {
+ fbModuleScanned[mod]=kFALSE;
}
Init();
}
-
+//_________________________________________________________________________
AliITSOnlineSPDscanAnalyzer::~AliITSOnlineSPDscanAnalyzer() {
// destructor
for (UInt_t hs=0; hs<6; hs++) {
for (UInt_t chipNr=0; chipNr<11; chipNr++) {
if (fMeanMultiplicity[hs][chipNr]!=NULL) {
delete fMeanMultiplicity[hs][chipNr];
+ fMeanMultiplicity[hs][chipNr]=NULL;
}
if (fHitEventEfficiency[hs][chipNr]!=NULL) {
delete fHitEventEfficiency[hs][chipNr];
+ fHitEventEfficiency[hs][chipNr]=NULL;
}
}
}
- if (fTriggers!=NULL) delete fTriggers;
- if (fScanObj!=NULL) delete fScanObj;
- for (Int_t module=0; module<240; module++) {
- if (fHandler[module]!=NULL) {
- delete fHandler[module];
- }
+
+ if (fTriggers!=NULL) {
+ delete fTriggers;
+ fTriggers=NULL;
}
-}
+ DeleteUniformityHistograms();
+
+ if (fScanObj!=NULL) {
+ delete fScanObj;
+ fScanObj=NULL;
+ }
+}
+//_________________________________________________________________________
AliITSOnlineSPDscanAnalyzer& AliITSOnlineSPDscanAnalyzer::operator=(const AliITSOnlineSPDscanAnalyzer& handle) {
- // assignment operator, only copies the filename (not the processed data)
+ // assignment operator, only copies the filename and params (not the processed data)
if (this!=&handle) {
for (UInt_t hs=0; hs<6; hs++) {
for (UInt_t chipNr=0; chipNr<11; chipNr++) {
}
}
}
- if (fTriggers!=NULL) delete fTriggers;
- if (fScanObj!=NULL) delete fScanObj;
- for (Int_t module=0; module<240; module++) {
- if (fHandler[module]!=NULL) {
- delete fHandler[module];
- }
+ if (fTriggers!=NULL) {
+ delete fTriggers;
+ fTriggers=NULL;
+ }
+
+ DeleteUniformityHistograms();
+
+ if (fScanObj!=NULL) {
+ delete fScanObj;
+ fScanObj=NULL;
}
- sprintf(fFileName,"%s",handle.fFileName);
+ fFileName=handle.fFileName;
+ fOverWrite=handle.fOverWrite;
+ fNoiseThreshold=handle.fNoiseThreshold;
+ fNoiseMinimumEvents=handle.fNoiseMinimumEvents;
+ fMinNrStepsBeforeIncrease=handle.fMinNrStepsBeforeIncrease;
+ fMinIncreaseFromBaseLine=handle.fMinIncreaseFromBaseLine;
+ fStepDownDacSafe=handle.fStepDownDacSafe;
+ fMaxBaseLineLevel=handle.fMaxBaseLineLevel;
- fScanObj=NULL;
- fType=99;
- fDacId=99;
for (UInt_t chipNr=0; chipNr<11; chipNr++) {
for (UInt_t hs=0; hs<6; hs++) {
fMeanMultiplicity[hs][chipNr]=NULL;
fHitEventEfficiency[hs][chipNr]=NULL;
}
}
- fTriggers=NULL;
- for (Int_t module=0; module<240; module++) {
- fHandler[module]=NULL;
+ for (UInt_t mod=0; mod<240; mod++) {
+ fbModuleScanned[mod]=kFALSE;
}
+
+ fHandler=NULL;
+ fType=99;
+ fDacId=99;
+
Init();
}
return *this;
}
-
-void AliITSOnlineSPDscanAnalyzer::Init() {
+//_________________________________________________________________________
+void AliITSOnlineSPDscanAnalyzer::Init(Bool_t readFromGridFile) {
// first checks type of container and then initializes container obj
- FILE* fp0 = fopen(fFileName, "r");
- if (fp0 == NULL) {
- return;
- }
- else {
- fclose(fp0);
+ if (!readFromGridFile) {
+ FILE* fp0 = fopen(fFileName.Data(), "r");
+ if (fp0 == NULL) {
+ return;
+ }
+ else {
+ fclose(fp0);
+ }
}
- fScanObj = new AliITSOnlineSPDscan(fFileName);
+
+ fScanObj = new AliITSOnlineSPDscan(fFileName.Data(),readFromGridFile);
fType = fScanObj->GetType();
delete fScanObj;
switch(fType) {
case kUNIMA:
case kNOISE:
- fScanObj = new AliITSOnlineSPDscanSingle(fFileName);
+ fScanObj = new AliITSOnlineSPDscanSingle(fFileName.Data(),readFromGridFile);
break;
case kMINTH:
case kDAC:
case kDELAY:
- fScanObj = new AliITSOnlineSPDscanMultiple(fFileName);
+ fScanObj = new AliITSOnlineSPDscanMultiple(fFileName.Data(),readFromGridFile);
fDacId = ((AliITSOnlineSPDscanMultiple*)fScanObj)->GetDacId();
break;
case kMEANTH:
- fScanObj = new AliITSOnlineSPDscanMeanTh(fFileName);
+ fScanObj = new AliITSOnlineSPDscanMeanTh(fFileName.Data(),readFromGridFile);
fDacId = ((AliITSOnlineSPDscanMeanTh*)fScanObj)->GetDacId();
break;
default:
- printf("Type %d not defined!\n",fType);
+ Error("AliITSOnlineSPDscanAnalyzer::Init","Type %d not defined!",fType);
fScanObj=NULL;
return;
break;
}
- // set some default values (these should later be read from text file)
- fOverWrite=kFALSE;
- fNoiseThreshold=0.01;
- fNoiseMinimumEvents=100;
- fMinNrStepsBeforeIncrease=6;
- fMinIncreaseFromBaseLine=2;
- fStepDownDacSafe=2;
- fMaxBaseLineLevel=10;
-
}
-
-
-Bool_t AliITSOnlineSPDscanAnalyzer::ProcessDeadPixels(Char_t *oldcalibDir) {
- // process dead pixel data
+//_________________________________________________________________________
+void AliITSOnlineSPDscanAnalyzer::SetParam(const Char_t *pname, const Char_t *pval) {
+ // set a parameter
+ TString name = pname;
+ TString val = pval;
+ if (name.CompareTo("fOverWrite")==0) {
+ if (val.CompareTo("YES")==0 || val.CompareTo("1")==0) {
+ fOverWrite = kTRUE;
+ }
+ else fOverWrite = kFALSE;
+ }
+ else if (name.CompareTo("fNoiseThreshold")==0) {
+ fNoiseThreshold = val.Atof();
+ }
+ else if (name.CompareTo("fNoiseMinimumEvents")==0) {
+ fNoiseMinimumEvents = val.Atoi();
+ }
+ else if (name.CompareTo("fMinNrStepsBeforeIncrease")==0) {
+ fMinNrStepsBeforeIncrease = val.Atoi();
+ }
+ else if (name.CompareTo("fMinIncreaseFromBaseLine")==0) {
+ fMinIncreaseFromBaseLine = val.Atof();
+ }
+ else if (name.CompareTo("fStepDownDacSafe")==0) {
+ fStepDownDacSafe = val.Atoi();
+ }
+ else if (name.CompareTo("fMaxBaseLineLevel")==0) {
+ fMaxBaseLineLevel = val.Atof();
+ }
+ else {
+ Error("AliITSOnlineSPDscanAnalyzer::SetParam","Parameter %s in configuration file unknown.",name.Data());
+ }
+}
+//_________________________________________________________________________
+void AliITSOnlineSPDscanAnalyzer::ReadParamsFromLocation(const Char_t *dirName) {
+ // opens file (default name) in dir dirName and reads parameters from it
+ TString paramsFileName = Form("%s/standal_params.txt",dirName);
+ ifstream paramsFile;
+ paramsFile.open(paramsFileName, ifstream::in);
+ if (paramsFile.fail()) {
+ printf("No config file (%s) present. Using default tuning parameters.\n",paramsFileName.Data());
+ }
+ else {
+ while(1) {
+ Char_t paramN[50];
+ Char_t paramV[50];
+ paramsFile >> paramN;
+ if (paramsFile.eof()) break;
+ paramsFile >> paramV;
+ SetParam(paramN,paramV);
+ if (paramsFile.eof()) break;
+ }
+ paramsFile.close();
+ }
+}
+//_________________________________________________________________________
+Bool_t AliITSOnlineSPDscanAnalyzer::IsChipPresent(UInt_t hs, UInt_t chipNr) {
+ // is the chip present?
+ if (fScanObj==NULL) {
+ Warning("AliITSOnlineSPDscanAnalyzer::IsChipPresent","No data!");
+ return kFALSE;
+ }
+ return fScanObj->GetChipPresent(hs,chipNr);
+}
+//_________________________________________________________________________
+Bool_t AliITSOnlineSPDscanAnalyzer::ProcessDeadPixels(/*Char_t *oldcalibDir*/) {
+ // process dead pixel data, for uniformity scan,
+ // NB: This will not be the general way of finding dead pixels.
if (fScanObj==NULL) {
- printf("No data!\n");
+ Warning("AliITSOnlineSPDscanAnalyzer::ProcessDeadPixels","No data!");
return kFALSE;
}
// should be type kUNIMA
if (fType!=kUNIMA) {
- printf("Dead pixels only for scan type %d\n",kUNIMA);
+ Warning("AliITSOnlineSPDscanAnalyzer::ProcessDeadPixels","Dead pixels only for scan type %d.",kUNIMA);
return kFALSE;
}
-
- Int_t nrDead[240];
- for (Int_t i=0; i<240; i++) {
- nrDead[i]=0;
+ // handler should be initialized
+ if (fHandler==NULL) {
+ Error("AliITSOnlineSPDscanAnalyzer::ProcessDeadPixels","Calibration handler is not initialized!");
+ return kFALSE;
}
+
UInt_t routerNr = fScanObj->GetRouterNr();
UInt_t rowStart = fScanObj->GetRowStart();
UInt_t rowEnd = fScanObj->GetRowEnd();
for (UInt_t hs=0; hs<6; hs++) {
for (UInt_t chipNr=0; chipNr<10; chipNr++) {
if (fScanObj->GetChipPresent(hs,chipNr) && fScanObj->GetAverageMultiplicity(0,hs,chipNr)>0) { // check the status of the chippresent parameter in the mood header!!!!!!!!!!!!!!!!!!!!!!!!!!!!
- UInt_t module = AliITSRawStreamSPD::GetModuleNumber(routerNr,hs*2+chipNr/5);
- if (!fHandler[module]) {
- fHandler[module] = new AliITSOnlineCalibrationSPDhandler(module);
+ if (fOverWrite) {fHandler->ResetDeadForChip(routerNr,hs,chipNr);}
+ for (UInt_t col=0; col<32; col++) {
+ for (UInt_t row=rowStart; row<=rowEnd; row++) {
+ if (col!=1 && col!=9 && col!=17 && col!=25) { //exclude test columns!!!
+ if (fScanObj->GetHits(0,hs,chipNr,col,row)==0) {
+ fHandler->SetDeadPixel(routerNr,hs,chipNr,col,row);
+ }
+ }
+ }
}
- fHandler[module]->SetFileLocation(oldcalibDir);
- fHandler[module]->ReadFromFile();
- if (fOverWrite) {fHandler[module]->ResetDead();}
+ }
+ }
+ }
+ return kTRUE;
+}
+//_________________________________________________________________________
+Bool_t AliITSOnlineSPDscanAnalyzer::ProcessUniformity() {
+ // process uniformity scan data (thanks to Roberta Ferretti for providing this method)
+ if (fScanObj==NULL) {
+ Warning("AliITSOnlineSPDscanAnalyzer::ProcessUniformity","No data!");
+ return kFALSE;
+ }
+ // should be type kUNIMA
+ if (fType!=kUNIMA) {
+ Warning("AliITSOnlineSPDscanAnalyzer::ProcessUniformity","Only for scan type %d.",kUNIMA);
+ return kFALSE;
+ }
+
+ CreateUniformityHistograms(); // create all histograms that will be filled here
+
+ // UInt_t routerNr = fScanObj->GetRouterNr();
+ UInt_t rowStart = fScanObj->GetRowStart();
+ UInt_t rowEnd = fScanObj->GetRowEnd();
+ UInt_t nrTriggers = fScanObj->GetTriggers(0)/(rowEnd-rowStart+1);
+
+ Float_t pixel100=0;
+ Float_t zeri=0;
+ Float_t pixelN=0;
+ UInt_t numChipsActive=0;
+
+ for (UInt_t hs=0; hs<6; hs++) {
+ Float_t pixel100hs=0;
+ Float_t zerihs=0;
+ Float_t pixelNhs=0;
+ UInt_t numChipsActiveHS=0;
+
+ for (UInt_t chipNr=0; chipNr<10; chipNr++) {
+ Float_t pixel100chip=0;
+ Float_t zerichip=0;
+ Float_t pixelNchip=0;
+
+ if (fScanObj->GetChipPresent(hs,chipNr)) { // check the status of the chippresent parameter in the mood header!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+ numChipsActive++;
+ numChipsActiveHS++;
+
for (UInt_t col=0; col<32; col++) {
for (UInt_t row=rowStart; row<=rowEnd; row++) {
if (col!=1 && col!=9 && col!=17 && col!=25) { //exclude test columns!!!
+
+ if (fScanObj->GetHits(0,hs,chipNr,col,row)==nrTriggers) {
+ pixel100++;
+ pixel100hs++;
+ pixel100chip++;
+ }
if (fScanObj->GetHits(0,hs,chipNr,col,row)==0) {
- Int_t newCol = 32*(chipNr%5) + col;
- if (fHandler[module]->SetDeadPixel(newCol,row)) {
- nrDead[module]++;
- }
+ zeri++;
+ zerihs++;
+ zerichip++;
+ }
+ if (fScanObj->GetHits(0,hs,chipNr,col,row)>nrTriggers) {
+ pixelN++;
+ pixelNhs++;
+ pixelNchip++;
}
}
}
}
+
+ Float_t tPeffChip=(pixel100chip/(28*(rowEnd-rowStart+1)))*100;
+ fTPeffChip[hs]->Fill(chipNr,tPeffChip);
+
+ Float_t deadPixelChip=(zerichip/(28*(rowEnd-rowStart+1)))*100;
+ fDeadPixelChip[hs]->Fill(chipNr,deadPixelChip);
+
+ Float_t noisyPixelChip=(pixelNchip/(28*(rowEnd-rowStart+1)))*100;
+ fNoisyPixelChip[hs]->Fill(chipNr,noisyPixelChip);
}
}
+
+ Float_t tPeffHS=(pixel100hs/(28*numChipsActiveHS*(rowEnd-rowStart+1)))*100;
+ fTPeffHS->Fill(hs,tPeffHS);
+
+ Float_t deadPixelHS=(zerihs/(28*numChipsActiveHS*(rowEnd-rowStart+1)))*100;
+ fDeadPixelHS->Fill(hs,deadPixelHS);
+
+ Float_t noisyPixelHS=(pixelNhs/(28*numChipsActiveHS*(rowEnd-rowStart+1)))*100;
+ fNoisyPixelHS->Fill(hs,noisyPixelHS);
}
+
+ fTPeff=(pixel100/(28*numChipsActive*(rowEnd-rowStart+1)))*100;
+ fDeadPixel=(zeri/(28*numChipsActive*(rowEnd-rowStart+1)))*100;
+ fNoisyPixel=(pixelN/(28*numChipsActive*(rowEnd-rowStart+1)))*100;
+
return kTRUE;
}
+//_________________________________________________________________________
+void AliITSOnlineSPDscanAnalyzer::CreateUniformityHistograms() {
+ // create uniformity histograms to be filled by "ProcessUniformity" method
+ DeleteUniformityHistograms(); // make sure no old histograms are lying around...
+ UInt_t eq = GetRouterNr();
+ TString label;
+
+ label = Form("Ratio of 'Good' Pixels Per HS (eq %d)",eq);
+ fTPeffHS = new TH1F(label.Data(),label.Data(),6,-0.5,5.5);
+ fTPeffHS->SetXTitle("hs");
+ fTPeffHS->SetYTitle("ratio [%]");
+ fTPeffHS->SetFillColor(kBlue);
+ fTPeffHS->SetStats(0);
+ label = Form("Ratio of 'Dead' Pixels Per HS (eq %d)",eq);
+ fDeadPixelHS = new TH1F(label.Data(),label.Data(),6,-0.5,5.5);
+ fDeadPixelHS->SetXTitle("hs");
+ fDeadPixelHS->SetYTitle("ratio [%]");
+ fDeadPixelHS->SetFillColor(kBlue);
+ fDeadPixelHS->SetStats(0);
+
+ label = Form("Ratio of 'Noisy' Pixels Per HS (eq %d)",eq);
+ fNoisyPixelHS = new TH1F(label.Data(),label.Data(),6,-0.5,5.5);
+ fNoisyPixelHS->SetXTitle("hs");
+ fNoisyPixelHS->SetYTitle("ratio [%]");
+ fNoisyPixelHS->SetFillColor(kBlue);
+ fNoisyPixelHS->SetStats(0);
+
+ for (UInt_t hs=0; hs<6; hs++) {
+ label = Form("Ratio of 'Good' Pixels Per Chip (eq %d, hs %d)",eq,hs);
+ fTPeffChip[hs] = new TH1F(label.Data(),label.Data(),10,-0.5,9.5);
+ fTPeffChip[hs]->SetXTitle("chip");
+ fTPeffChip[hs]->SetYTitle("ratio [%]");
+ fTPeffChip[hs]->SetFillColor(kBlue);
+ fTPeffChip[hs]->SetStats(0);
+ label = Form("Ratio of 'Dead' Pixels Per Chip (eq %d, hs %d)",eq,hs);
+ fDeadPixelChip[hs] = new TH1F(label.Data(),label.Data(),10,-0.5,9.5);
+ fDeadPixelChip[hs]->SetXTitle("chip");
+ fDeadPixelChip[hs]->SetYTitle("ratio [%]");
+ fDeadPixelChip[hs]->SetFillColor(kBlue);
+ fDeadPixelChip[hs]->SetStats(0);
-Bool_t AliITSOnlineSPDscanAnalyzer::ProcessNoisyPixels(Char_t *oldcalibDir) {
+ label = Form("Ratio of 'Noisy' Pixels Per Chip (eq %d, hs %d)",eq,hs);
+ fNoisyPixelChip[hs] = new TH1F(label.Data(),label.Data(),10,-0.5,9.5);
+ fNoisyPixelChip[hs]->SetXTitle("chip");
+ fNoisyPixelChip[hs]->SetYTitle("ratio [%]");
+ fNoisyPixelChip[hs]->SetFillColor(kBlue);
+ fNoisyPixelChip[hs]->SetStats(0);
+ }
+
+}
+//_________________________________________________________________________
+void AliITSOnlineSPDscanAnalyzer::DeleteUniformityHistograms() {
+ // remove uniformity histograms if they are created
+ if (fTPeffHS!=NULL) {
+ delete fTPeffHS;
+ fTPeffHS=NULL;
+ }
+ if (fDeadPixelHS!=NULL) {
+ delete fDeadPixelHS;
+ fDeadPixelHS=NULL;
+ }
+ if (fNoisyPixelHS!=NULL) {
+ delete fNoisyPixelHS;
+ fNoisyPixelHS=NULL;
+ }
+ for (UInt_t hs=0; hs<6; hs++) {
+ if (fTPeffChip[hs]!=NULL) {
+ delete fTPeffChip[hs];
+ fTPeffChip[hs]=NULL;
+ }
+ if (fDeadPixelChip[hs]!=NULL) {
+ delete fDeadPixelChip[hs];
+ fDeadPixelChip[hs]=NULL;
+ }
+ if (fNoisyPixelChip[hs]!=NULL) {
+ delete fNoisyPixelChip[hs];
+ fNoisyPixelChip[hs]=NULL;
+ }
+ }
+}
+//_________________________________________________________________________
+Bool_t AliITSOnlineSPDscanAnalyzer::ProcessNoisyPixels(/*Char_t *oldcalibDir*/) {
// process noisy pixel data
if (fScanObj==NULL) {
- printf("No data!\n");
+ Warning("AliITSOnlineSPDscanAnalyzer::ProcessNoisyPixels","No data!");
return kFALSE;
}
// should be type kNOISE
if (fType != kNOISE) {
- printf("Noisy pixels only for scan type %d\n",kNOISE);
+ Warning("AliITSOnlineSPDscanAnalyzer::ProcessNoisyPixels","Noisy pixels only for scan type %d.",kNOISE);
return kFALSE;
}
-
- Int_t nrNoisy[240];
- for (Int_t i=0; i<240; i++) {
- nrNoisy[i]=0;
+ // handler should be initialized
+ if (fHandler==NULL) {
+ Error("AliITSOnlineSPDscanAnalyzer::ProcessNoisyPixels","Calibration handler is not initialized!");
+ return kFALSE;
}
+ // check if enough statistics
if (fScanObj->GetTriggers(0)<fNoiseMinimumEvents) {
- printf("Error: Process noisy: Too few events.\n");
+ Warning("AliITSOnlineSPDscanAnalyzer::ProcessNoisyPixels","Process noisy: Too few events.");
return kFALSE;
}
+
UInt_t routerNr = fScanObj->GetRouterNr();
for (UInt_t hs=0; hs<6; hs++) {
for (UInt_t chipNr=0; chipNr<10; chipNr++) {
if (fScanObj->GetChipPresent(hs,chipNr)) { // check the status of the chippresent parameter in the mood header!!!!!!!!!!!!!!!!!!!!!!!!!!!!
- UInt_t module = AliITSRawStreamSPD::GetModuleNumber(routerNr,hs*2+chipNr/5);
+ if (fOverWrite) {fHandler->ResetNoisyForChip(routerNr,hs,chipNr);}
for (UInt_t col=0; col<32; col++) {
for (UInt_t row=0; row<256; row++) {
if (fScanObj->GetHitsEfficiency(0,hs,chipNr,col,row)>fNoiseThreshold) {
- if (!fHandler[module]) {
- fHandler[module] = new AliITSOnlineCalibrationSPDhandler(module);
- }
- fHandler[module]->SetFileLocation(oldcalibDir);
- fHandler[module]->ReadFromFile();
- if (fOverWrite) {fHandler[module]->ResetNoisy();}
- Int_t newCol = 32*(chipNr%5) + col;
- if (fHandler[module]->SetNoisyPixel(newCol,row)) {
- nrNoisy[module]++;
- }
+ fHandler->SetNoisyPixel(routerNr,hs,chipNr,col,row);
}
}
}
}
return kTRUE;
}
-
-Bool_t AliITSOnlineSPDscanAnalyzer::SaveDeadNoisyPixels(UInt_t module, Char_t *calibDir) {
- // save dead and noisy pixels to file in dir calibDir
- if (fHandler[module]!=NULL) {
- fHandler[module]->SetFileLocation(calibDir);
- fHandler[module]->WriteToFile();
- return kTRUE;
- }
- return kFALSE;
-}
-
+//_________________________________________________________________________
Int_t AliITSOnlineSPDscanAnalyzer::GetDelay(UInt_t hs, UInt_t chipNr) {
// get delay
if (hs>=6 || chipNr>10) return -1;
if (fScanObj==NULL) {
- printf("No data!\n");
+ Warning("AliITSOnlineSPDscanAnalyzer::GetDelay","No data!");
return -1;
}
// should be type kDELAY or kDAC with id 42 (delay_ctrl)
if (fType!=kDELAY && (fType!=kDAC || fDacId!=42)) {
- printf("Delay only for scan type %d or %d and dac_id 42\n",kDELAY,kDAC);
+ Warning("AliITSOnlineSPDscanAnalyzer::GetDelay","Delay only for scan type %d or %d and dac_id 42.",kDELAY,kDAC);
return -1;
}
if (fMeanMultiplicity[hs][chipNr]==NULL) {
}
}
- Char_t funcName[30];
- sprintf(funcName,"Fit delay func HS%d CHIP%d",hs,chipNr);
- Int_t minDac = ((AliITSOnlineSPDscanMultiple*)fScanObj)->GetDacValue(0);
- Int_t maxDac = ((AliITSOnlineSPDscanMultiple*)fScanObj)->GetDacValue(fScanObj->GetNSteps()-1);
- TF1* delayFunc = new TF1(funcName,"gaus",minDac,maxDac);
- fMeanMultiplicity[hs][chipNr]->Fit(funcName,"Q0");
- Double_t mean = delayFunc->GetParameter(1);
- // Double_t sigma = fDelayFunc[hs][chipNr]->GetParameter(2);
- delete delayFunc;
- if (mean>minDac && mean<maxDac) {
- return (Int_t) (mean+0.5);
+ UInt_t maxStep=0;
+ Float_t maxVal=0;
+ for (UInt_t step=0; step<fScanObj->GetNSteps(); step++) {
+ Double_t thisDac;
+ Double_t thisMult;
+ fMeanMultiplicity[hs][chipNr]->GetPoint(step,thisDac,thisMult);
+ if (thisMult > maxVal) {
+ maxVal = thisMult;
+ maxStep = step;
+ }
+ }
+
+ if (maxVal>0) {
+ return ((AliITSOnlineSPDscanMultiple*)fScanObj)->GetDacValue(maxStep);
}
else {
return -1;
}
-}
+}
+//_________________________________________________________________________
Int_t AliITSOnlineSPDscanAnalyzer::GetNrNoisyUnima(UInt_t hs, UInt_t chipNr) {
// in case of a uniformity scan, returns the nr of noisy pixels, (here > 200 hits)
if (hs>=6 || chipNr>10) return -1;
if (fScanObj==NULL) {
- printf("No data!\n");
+ Error("AliITSOnlineSPDscanAnalyzer::GetNrNoisyUnima","No data!");
return kFALSE;
}
// should be type kUNIMA
if (fType != kUNIMA) {
- printf("Noisy pixels Unima only for scan type %d\n",kUNIMA);
+ Error("AliITSOnlineSPDscanAnalyzer::GetNrNoisyUnima","Noisy pixels Unima only for scan type %d.",kUNIMA);
return kFALSE;
}
if (fScanObj->GetTriggers(0)!=25600) {
- printf("Error: Process noisy unima: Incorrect number of events (!=25600.\n");
+ Error("AliITSOnlineSPDscanAnalyzer::GetNrNoisyUnima","Process noisy unima: Incorrect number of events (!=25600.");
return kFALSE;
}
}
return nrNoisy;
}
-
+//_________________________________________________________________________
Int_t AliITSOnlineSPDscanAnalyzer::FindLastMinThDac(UInt_t hs, UInt_t chipNr) {
// returns dac value where fMinIncreaseFromBaseLine reached
if (hs>=6 || chipNr>10) return -1;
}
return step;
}
-
+//_________________________________________________________________________
Float_t AliITSOnlineSPDscanAnalyzer::GetCompareLine(UInt_t step, UInt_t hs, UInt_t chipNr, Float_t basePar2) {
// returns value to compare mean mult with (when finding min th)
if (hs>=6 || chipNr>10) return -1;
// calculates and returns the minimum threshold
if (hs>=6 || chipNr>10) return -1;
if (fScanObj==NULL) {
- printf("No data!\n");
+ Error("AliITSOnlineSPDscanAnalyzer::GetMinTh","No data!");
return -1;
}
// should be type kMINTH or kDAC with id 39 (pre_vth)
if (fType!=kMINTH && (fType!=kDAC || fDacId!=39)) {
- printf("MinTh only for scan type %d OR %d with dac_id 39\n",kMINTH,kDAC);
+ Error("AliITSOnlineSPDscanAnalyzer::GetMinTh","MinTh only for scan type %d OR %d with dac_id 39.",kMINTH,kDAC);
return -1;
}
if (fMeanMultiplicity[hs][chipNr]==NULL) {
Int_t lastDac = FindLastMinThDac(hs,chipNr);
if (lastDac==-1) {
- printf("GetMinTh: HS%d, Chip%d: Increase of Mean Multiplicity by %1.2f never reached.\n",hs,chipNr,fMinIncreaseFromBaseLine);
+ Warning("AliITSOnlineSPDscanAnalyzer::GetMinTh","HS%d, Chip%d: Increase of Mean Multiplicity by %1.2f never reached.",hs,chipNr,fMinIncreaseFromBaseLine);
return -1;
}
Int_t minDac = ((AliITSOnlineSPDscanMultiple*)fScanObj)->GetDacValue(0);
- Char_t funcName[30];
- sprintf(funcName,"Fit minth func HS%d CHIP%d",hs,chipNr);
- TF1 *minThFunc = new TF1(funcName,itsSpdErrorf,100,500,3);
+ TString funcName = Form("Fit minth func HS%d CHIP%d",hs,chipNr);
+ TF1 *minThFunc = new TF1(funcName.Data(),itsSpdErrorf,100,500,3);
minThFunc->SetParameter(0,lastDac+10);
minThFunc->SetParameter(1,2);
minThFunc->SetParameter(2,0);
delete minThFunc;
if (baseLine>fMaxBaseLineLevel) {
- printf("GetMinTh: HS%d, Chip%d: BaseLine too large (%1.2f>%1.2f)\n",hs,chipNr,baseLine,fMaxBaseLineLevel);
+ Warning("AliITSOnlineSPDscanAnalyzer::GetMinTh","HS%d, Chip%d: BaseLine too large (%1.2f>%1.2f).",hs,chipNr,baseLine,fMaxBaseLineLevel);
return -1;
}
UInt_t step=FindClosestLowerStep(lastDac);
Float_t compareLine = GetCompareLine(step,hs,chipNr,baseLine);
if (compareLine==-1) {
- printf("GetMinTh: HS%d, Chip%d: Not enough steps (%d<%d) before increase to get a compare line.\n",hs,chipNr,step,fMinNrStepsBeforeIncrease);
+ Warning("AliITSOnlineSPDscanAnalyzer::GetMinTh","HS%d, Chip%d: Not enough steps (%d<%d) before increase to get a compare line.",hs,chipNr,step,fMinNrStepsBeforeIncrease);
return -1;
}
return minth;
}
else {
- printf("GetMinTh: HS%d, Chip%d: Did not find a point below the compare line (%f).\n",hs,chipNr,compareLine);
+ Warning("AliITSOnlineSPDscanAnalyzer::GetMinTh","HS%d, Chip%d: Did not find a point below the compare line (%f).",hs,chipNr,compareLine);
return -1;
}
}
+//_________________________________________________________________________
+TArrayI AliITSOnlineSPDscanAnalyzer::GetMeanTh(UInt_t hs, UInt_t chipNr) {
+ // calculates and returns the mean threshold
+ TArrayI fitresults(4);
+ if (hs>=6 || chipNr>10) return fitresults;
+ if (fScanObj==NULL) {
+ Error("AliITSOnlineSPDscanAnalyzer::GetMeanTh","No data!");
+ return fitresults;
+ }
+ // should be type kMEANTH or kDAC with id 39
+ if (fType!=kMEANTH && (fType!=kDAC || fDacId!=105)) {
+ Error("AliITSOnlineSPDscanAnalyzer::GetMeanTh","MeanTh only for scan type %d OR %d with dac_id 105.",kMEANTH,kDAC);
+ return fitresults;
+ }
+ if (fHitEventEfficiency[hs][chipNr]==NULL) {
+ printf("processing hit efficiency \n");
+ if (!ProcessHitEventEfficiency()) {
+ printf("...not processed!!\n");
+ return fitresults;
+ }
+ }
+ Double_t x,y;
+ fHitEventEfficiency[hs][chipNr]->GetPoint(fHitEventEfficiency[hs][chipNr]->GetN()-1,x,y);
+ Double_t min = x;
+ fHitEventEfficiency[hs][chipNr]->GetPoint(0,x,y);
+ Double_t max = x;
+ Double_t mean = 0.5*(min+max);
+ TString funcName = Form("Fit meanth func HS%d CHIP%d",hs,chipNr);
+ TF1 *minThFunc = new TF1(funcName.Data(),itsSpdScurveForMeanTh,min,max,3);
+ minThFunc->SetParameter(0,mean);
+ minThFunc->SetParameter(1,264); // 4 (mV) * 66 (el)
+ minThFunc->SetParLimits(1,100,1000); // minimum value is 1 mV (-> 100 in TPAmplitude)
+ minThFunc->SetParameter(2,0.4);
+ minThFunc->SetParName(0,"Mean");
+ minThFunc->SetParName(1,"Sigma");
+ minThFunc->SetParName(2,"Half");
+ minThFunc->SetLineWidth(1);
+ fHitEventEfficiency[hs][chipNr]->Fit(funcName,"Q","",min,max);
+
+ fitresults.AddAt((Int_t)minThFunc->GetParameter(0),0);
+ fitresults.AddAt((Int_t)minThFunc->GetParError(0),1);
+ fitresults.AddAt((Int_t)minThFunc->GetParameter(1),2);
+ fitresults.AddAt((Int_t)minThFunc->GetParError(1),3);
+ TLine *ly = new TLine((Double_t)min,0.5,(Double_t)fitresults.At(0),0.5); ly->SetLineStyle(6);
+ ly->Draw("same");
+ TLine *lx = new TLine((Double_t)fitresults.At(0),0.,(Double_t)fitresults.At(0),0.5);
+ lx->SetLineStyle(6);
+ lx->Draw("same");
+ delete minThFunc;
+
+ return fitresults;
+}
+//_________________________________________________________________________
Bool_t AliITSOnlineSPDscanAnalyzer::ProcessMeanMultiplicity() {
// process mean multiplicity data
if (fScanObj==NULL) {
- printf("No data!\n");
+ Error("AliITSOnlineSPDscanAnalyzer::ProcessMeanMultiplicity","No data!");
return kFALSE;
}
for (UInt_t step=0; step<fScanObj->GetNSteps(); step++) {
}
return kTRUE;
}
-
+//_________________________________________________________________________
TGraph* AliITSOnlineSPDscanAnalyzer::GetMeanMultiplicityG(UInt_t hs, UInt_t chipNr) {
// returns mean multiplicity graph
if (hs>=6 || chipNr>10) return NULL;
}
return fMeanMultiplicity[hs][chipNr];
}
-
+//_________________________________________________________________________
Bool_t AliITSOnlineSPDscanAnalyzer::ProcessHitEventEfficiency() {
// process hit event efficiency
if (fScanObj==NULL) {
- printf("No data!\n");
+ Error("AliITSOnlineSPDscanAnalyzer::ProcessHitEventEfficiency","No data!");
return kFALSE;
}
for (UInt_t step=0; step<fScanObj->GetNSteps(); step++) {
fHitEventEfficiency[hs][chipNr] = new TGraph();
}
Float_t efficiency=fScanObj->GetHitEventsEfficiency(step,hs,chipNr);
- if (fType==kMINTH || fType==kMEANTH || fType==kDAC || fType==kDELAY) {
+ if (fType==kMINTH || fType==kDAC || fType==kDELAY) {
fHitEventEfficiency[hs][chipNr]->SetPoint(step,((AliITSOnlineSPDscanMultiple*)fScanObj)->GetDacValue(step),efficiency);
- }
- else {
+ } else if(fType==kMEANTH){
+ fHitEventEfficiency[hs][chipNr]->SetPoint(step,((AliITSOnlineSPDscanMeanTh*)fScanObj)->GetTPAmp(step,hs),efficiency);
+ } else {
fHitEventEfficiency[hs][chipNr]->SetPoint(step,0,efficiency);
}
}
}
return kTRUE;
}
-
+//_________________________________________________________________________
TGraph* AliITSOnlineSPDscanAnalyzer::GetHitEventEfficiencyG(UInt_t hs, UInt_t chipNr) {
// returns hit event efficiency graph
if (hs>=6 || chipNr>10) return NULL;
}
return fHitEventEfficiency[hs][chipNr];
}
-
-
+//_________________________________________________________________________
Bool_t AliITSOnlineSPDscanAnalyzer::ProcessNrTriggers() {
// process nr of triggers data
if (fScanObj==NULL) {
- printf("No data!\n");
+ Error("AliITSOnlineSPDscanAnalyzer::ProcessNrTriggers","No data!");
return kFALSE;
}
for (UInt_t step=0; step<fScanObj->GetNSteps(); step++) {
}
return kTRUE;
}
-
+//_________________________________________________________________________
TGraph* AliITSOnlineSPDscanAnalyzer::GetNrTriggersG() {
// returns nr of triggers graph
if (fTriggers==NULL) {
}
return fTriggers;
}
-
+//_________________________________________________________________________
Bool_t AliITSOnlineSPDscanAnalyzer::GetHalfStavePresent(UInt_t hs) {
// returns half stave present info
if (hs<6 && fScanObj!=NULL) {
}
return kFALSE;
}
-
-AliITSOnlineCalibrationSPDhandler* AliITSOnlineSPDscanAnalyzer::GetOnlineCalibrationHandler(UInt_t module) {
- // returns a pointer to the AliITSOnlineCalibrationSPDhandler
- if (module<240) return fHandler[module];
- else return NULL;
-}
-
+//_________________________________________________________________________
UInt_t AliITSOnlineSPDscanAnalyzer::GetRouterNr() {
// returns the router nr of scan obj
if (fScanObj!=NULL) return fScanObj->GetRouterNr();
else return 99;
}
-
+//_________________________________________________________________________
TH2F* AliITSOnlineSPDscanAnalyzer::GetHitMapTot(UInt_t step) {
// creates and returns a pointer to a hitmap histo (half sector style a la spdmood)
if (fScanObj==NULL) {
- printf("No data!\n");
+ Error("AliITSOnlineSPDscanAnalyzer::GetHitMapTot","No data!");
return NULL;
}
- Char_t histoname[50];
+ TString histoname;
if (fType==kMINTH || fType==kMEANTH || fType==kDAC) {
- sprintf(histoname,"Router %d , DAC %d",GetRouterNr(),((AliITSOnlineSPDscanMultiple*)fScanObj)->GetDacValue(step));
+ histoname = Form("Router %d , DAC %d",GetRouterNr(),((AliITSOnlineSPDscanMultiple*)fScanObj)->GetDacValue(step));
}
else {
- sprintf(histoname,"Router %d ",GetRouterNr());
+ histoname = Form("Router %d ",GetRouterNr());
}
- TH2F* fHitMapTot = new TH2F(histoname,histoname,32*10,-0.5,32*10-0.5,256*6,-0.5,256*6-0.5);
+ TH2F* fHitMapTot = new TH2F(histoname.Data(),histoname.Data(),32*10,-0.5,32*10-0.5,256*6,-0.5,256*6-0.5);
fHitMapTot->SetNdivisions(-10,"X");
fHitMapTot->SetNdivisions(-006,"Y");
fHitMapTot->SetTickLength(0,"X");
}
return fHitMapTot;
}
+//_________________________________________________________________________
+TH2F* AliITSOnlineSPDscanAnalyzer::GetPhysicalHitMapTot(UInt_t step) {
+ // creates and returns a pointer to a hitmap histo (-> 1 equpment opened up)
+
+ if (fScanObj==NULL) {
+ Error("AliITSOnlineSPDscanAnalyzer::GetHitMapTot","No data!");
+ return NULL;
+ }
+ TString histoname;
+ if (fType==kMINTH || fType==kMEANTH || fType==kDAC) {
+ histoname = Form("Router %d , DAC %d",GetRouterNr(),((AliITSOnlineSPDscanMultiple*)fScanObj)->GetDacValue(step));
+ }
+ else {
+ histoname = Form("Router %d ",GetRouterNr());
+ }
+ TH2F* hPhysicalHitMapTot = new TH2F(histoname.Data(),histoname.Data(),32*10,-0.5,32*10-0.5,256*6,-0.5,256*6-0.5);
+ hPhysicalHitMapTot->SetNdivisions(-10,"X");
+ hPhysicalHitMapTot->SetNdivisions(-006,"Y");
+ hPhysicalHitMapTot->SetTickLength(0,"X");
+ hPhysicalHitMapTot->SetTickLength(0,"Y");
+ hPhysicalHitMapTot->GetXaxis()->SetLabelColor(gStyle->GetCanvasColor());
+ hPhysicalHitMapTot->GetYaxis()->SetLabelColor(gStyle->GetCanvasColor());
+ Int_t correctChip = -1;
+ for (UInt_t hs=0; hs<6; hs++) {
+ for (UInt_t chipNr=0; chipNr<10; chipNr++) {
+ if(GetRouterNr()<10) correctChip = 9-chipNr;
+ else correctChip=chipNr;
+ for (UInt_t col=0; col<32; col++) {
+ for (UInt_t row=0; row<256; row++) {
+ if(hs<2) hPhysicalHitMapTot->Fill(correctChip*32+col,(5-hs)*256+row,fScanObj->GetHits(step,hs,chipNr,col,row));
+ else hPhysicalHitMapTot->Fill(correctChip*32+(31-col),(5-hs)*256+row,fScanObj->GetHits(step,hs,chipNr,col,row));
+ }
+ }
+ }
+ }
+ return hPhysicalHitMapTot;
+}
+
+//_________________________________________________________________________
+TH2F* AliITSOnlineSPDscanAnalyzer::GetHitMapChip(UInt_t step, UInt_t hs, UInt_t chip) {
+ // creates and returns a pointer to a hitmap histo (chip style a la spdmood)
+ if (fScanObj==NULL) {
+ Error("AliITSOnlineSPDscanAnalyzer::GetHitMapChip","No data!");
+ return NULL;
+ }
+
+ TString histoName;
+ TString histoTitle;
+ histoName = Form("fChipHisto_%d_%d_%d", GetRouterNr(), hs, chip);
+ histoTitle = Form("Eq ID %d, Half Stave %d, Chip %d ", GetRouterNr(), hs, chip);
+
+ TH2F *returnHisto = new TH2F(histoName.Data(), histoTitle.Data(), 32, -0.5, 31.5, 256, -0.5, 255.5);
+ returnHisto->SetMinimum(0);
+ for (UInt_t col=0; col<32; col++) {
+ for (UInt_t row=0; row<256; row++) {
+ if(hs<2) returnHisto->Fill(31-col,row,fScanObj->GetHits(step,hs,chip,col,row));
+ else returnHisto->Fill(col,row,fScanObj->GetHits(step,hs,chip,col,row));
+ }
+ }
+
+ return returnHisto;
+}
+
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff