#include <fstream>
#include <sstream>
#include <iostream>
+#include <stdexcept>
#include <cmath>
#include "AliRawReader.h"
TObject(),
fPedestalLowGain(),
fPedestalHighGain(),
- fPedestalRMSLowGain(),
- fPedestalRMSHighGain(),
+ fPedestalLEDRefLowGain(),
+ fPedestalLEDRefHighGain(),
fPeakMinusPedLowGain(),
fPeakMinusPedHighGain(),
fPeakMinusPedHighGainHisto(),
fPedestalLowGainDiff(),
fPedestalHighGainDiff(),
+ fPedestalLEDRefLowGainDiff(),
+ fPedestalLEDRefHighGainDiff(),
fPeakMinusPedLowGainDiff(),
fPeakMinusPedHighGainDiff(),
fPedestalLowGainRatio(),
fPedestalHighGainRatio(),
+ fPedestalLEDRefLowGainRatio(),
+ fPedestalLEDRefHighGainRatio(),
fPeakMinusPedLowGainRatio(),
fPeakMinusPedHighGainRatio(),
fDeadMap(),
fDetType(kNone),
fColumns(0),
fRows(0),
+ fLEDRefs(0),
fModules(0),
fRowMin(0),
fRowMax(0),
if (detectorType == kPhos) {
fColumns = fgkPhosCols;
fRows = fgkPhosRows;
+ fLEDRefs = fgkPhosLEDRefs;
fModules = fgkPhosModules;
fCaloString = "PHOS";
fRowMin = -1*fRows;
//case, if someone intentionally gives another number
fColumns = AliEMCALGeoParams::fgkEMCALCols;
fRows = AliEMCALGeoParams::fgkEMCALRows;
+ fLEDRefs = AliEMCALGeoParams::fgkEMCALLEDRefs;
fModules = AliEMCALGeoParams::fgkEMCALModules;
fCaloString = "EMCAL";
fRowMin = 0;
fRowMultiplier = 1;
}
fDetType = detectorType;
-
+
+ // ValidateProfiles(); // not to be done in ctor; info from Axel N.
+}
+
+//_____________________________________________________________________
+void AliCaloCalibPedestal::ValidateProfiles()
+{
+ //Make sure the basic histos exist
+ if (!fPedestalLowGain.IsEmpty()) return; //The profiles already exist. We just check one, because they're all created at
+ //the same time
+
//Then, loop for the requested number of modules
TString title, name;
for (int i = 0; i < fModules; i++) {
fPedestalHighGain.Add(new TProfile2D(name, title,
fColumns, 0.0, fColumns,
fRows, fRowMin, fRowMax,"s"));
- //All Samples, low gain
- name = "hPedestalRMSlowgain";
+
+ //LED Ref/Mon pedestals, low gain
+ name = "hPedestalLEDReflowgain";
name += i;
- title = "Pedestal RMS, low gain, module ";
+ title = "Pedestal LEDRef, low gain, module ";
title += i;
- fPedestalRMSLowGain.Add(new TProfile2D(name, title,
- fColumns, 0.0, fColumns,
- fRows, fRowMin, fRowMax,"s"));
-
- //All Samples, high gain
- name = "hPedestalRMShighgain";
+ fPedestalLEDRefLowGain.Add(new TProfile(name, title,
+ fLEDRefs, 0.0, fLEDRefs, "s"));
+
+ //LED Ref/Mon pedestals, high gain
+ name = "hPedestalLEDRefhighgain";
name += i;
- title = "Pedestal RMS, high gain, module ";
+ title = "Pedestal LEDRef, high gain, module ";
title += i;
- fPedestalRMSHighGain.Add(new TProfile2D(name, title,
- fColumns, 0.0, fColumns,
- fRows, fRowMin, fRowMax,"s"));
+ fPedestalLEDRefHighGain.Add(new TProfile(name, title,
+ fLEDRefs, 0.0, fLEDRefs, "s"));
//Peak-Pedestals, low gain
name = "hPeakMinusPedlowgain";
fRows, fRowMin, fRowMax));
}//end for nModules create the histograms
-
+
+ CompressAndSetOwner();
+}
+
+//_____________________________________________________________________
+void AliCaloCalibPedestal::CompressAndSetOwner()
+{
//Compress the arrays, in order to remove the empty objects (a 16 slot array is created by default)
fPedestalLowGain.Compress();
fPedestalHighGain.Compress();
- fPedestalRMSLowGain.Compress();
- fPedestalRMSHighGain.Compress();
+ fPedestalLEDRefLowGain.Compress();
+ fPedestalLEDRefHighGain.Compress();
fPeakMinusPedLowGain.Compress();
fPeakMinusPedHighGain.Compress();
fPeakMinusPedHighGainHisto.Compress();
fDeadMap.Compress();
- //Make them the owners of the profiles, so we don't need to care about deleting them
- //fPedestalLowGain.SetOwner();
- //fPedestalHighGain.SetOwner();
- //fPeakMinusPedLowGain.SetOwner();
- //fPeakMinusPedHighGain.SetOwner();
-
+
+ // set owner ship for everyone
+ fPedestalLowGain.SetOwner(kTRUE);
+ fPedestalHighGain.SetOwner(kTRUE);
+ fPedestalLEDRefLowGain.SetOwner(kTRUE);
+ fPedestalLEDRefHighGain.SetOwner(kTRUE);
+ fPeakMinusPedLowGain.SetOwner(kTRUE);
+ fPeakMinusPedHighGain.SetOwner(kTRUE);
+ fPeakMinusPedHighGainHisto.SetOwner(kTRUE);
+ fPedestalLowGainDiff.SetOwner(kTRUE);
+ fPedestalHighGainDiff.SetOwner(kTRUE);
+ fPedestalLEDRefLowGainDiff.SetOwner(kTRUE);
+ fPedestalLEDRefHighGainDiff.SetOwner(kTRUE);
+ fPeakMinusPedLowGainDiff.SetOwner(kTRUE);
+ fPeakMinusPedHighGainDiff.SetOwner(kTRUE);
+ fPedestalLowGainRatio.SetOwner(kTRUE);
+ fPedestalHighGainRatio.SetOwner(kTRUE);
+ fPedestalLEDRefLowGainRatio.SetOwner(kTRUE);
+ fPedestalLEDRefHighGainRatio.SetOwner(kTRUE);
+ fPeakMinusPedLowGainRatio.SetOwner(kTRUE);
+ fPeakMinusPedHighGainRatio.SetOwner(kTRUE);
+ fDeadMap.SetOwner(kTRUE);
}
// dtor
//_____________________________________________________________________
AliCaloCalibPedestal::~AliCaloCalibPedestal()
{
+ //dtor
+
if (fReference) delete fReference;//Delete the reference object, if it has been loaded
- //TObjArray will delete the histos/profiles when it is deleted.
+
+ // delete also TObjArray's
+ fPedestalLowGain.Delete();
+ fPedestalHighGain.Delete();
+ fPedestalLEDRefLowGain.Delete();
+ fPedestalLEDRefHighGain.Delete();
+ fPeakMinusPedLowGain.Delete();
+ fPeakMinusPedHighGain.Delete();
+ fPeakMinusPedHighGainHisto.Delete();
+ fPedestalLowGainDiff.Delete();
+ fPedestalHighGainDiff.Delete();
+ fPedestalLEDRefLowGainDiff.Delete();
+ fPedestalLEDRefHighGainDiff.Delete();
+ fPeakMinusPedLowGainDiff.Delete();
+ fPeakMinusPedHighGainDiff.Delete();
+ fPedestalLowGainRatio.Delete();
+ fPedestalHighGainRatio.Delete();
+ fPedestalLEDRefLowGainRatio.Delete();
+ fPedestalLEDRefHighGainRatio.Delete();
+ fPeakMinusPedLowGainRatio.Delete();
+ fPeakMinusPedHighGainRatio.Delete();
+ fDeadMap.Delete();
+
}
// copy ctor
//_____________________________________________________________________
-AliCaloCalibPedestal::AliCaloCalibPedestal(const AliCaloCalibPedestal &ped) :
+AliCaloCalibPedestal::AliCaloCalibPedestal(AliCaloCalibPedestal &ped) :
TObject(ped),
fPedestalLowGain(),
fPedestalHighGain(),
- fPedestalRMSLowGain(),
- fPedestalRMSHighGain(),
+ fPedestalLEDRefLowGain(),
+ fPedestalLEDRefHighGain(),
fPeakMinusPedLowGain(),
fPeakMinusPedHighGain(),
fPeakMinusPedHighGainHisto(),
fPedestalLowGainDiff(),
fPedestalHighGainDiff(),
+ fPedestalLEDRefLowGainDiff(),
+ fPedestalLEDRefHighGainDiff(),
fPeakMinusPedLowGainDiff(),
fPeakMinusPedHighGainDiff(),
fPedestalLowGainRatio(),
fPedestalHighGainRatio(),
+ fPedestalLEDRefLowGainRatio(),
+ fPedestalLEDRefHighGainRatio(),
fPeakMinusPedLowGainRatio(),
fPeakMinusPedHighGainRatio(),
fDeadMap(),
fDetType(ped.GetDetectorType()),
fColumns(ped.GetColumns()),
fRows(ped.GetRows()),
+ fLEDRefs(ped.GetLEDRefs()),
fModules(ped.GetModules()),
fRowMin(ped.GetRowMin()),
fRowMax(ped.GetRowMax()),
for (int i = 0; i < fModules; i++) {
fPedestalLowGain.Add( ped.GetPedProfileLowGain(i) );
fPedestalHighGain.Add( ped.GetPedProfileHighGain(i) );
- fPedestalRMSLowGain.Add( ped.GetPedRMSProfileLowGain(i) );
- fPedestalRMSHighGain.Add( ped.GetPedRMSProfileHighGain(i) );
+ fPedestalLEDRefLowGain.Add( ped.GetPedLEDRefProfileLowGain(i) );
+ fPedestalLEDRefHighGain.Add( ped.GetPedLEDRefProfileHighGain(i) );
fPeakMinusPedLowGain.Add( ped.GetPeakProfileLowGain(i) );
fPeakMinusPedHighGain.Add( ped.GetPeakProfileHighGain(i) );
fPeakMinusPedHighGainHisto.Add( ped.GetPeakHighGainHisto(i) );
fDeadMap.Add( ped.GetDeadMap(i) );
}//end for nModules
- //Compress the arrays, in order to remove the empty objects (a 16 slot array is created by default)
- fPedestalLowGain.Compress();
- fPedestalHighGain.Compress();
- fPedestalRMSLowGain.Compress();
- fPedestalRMSHighGain.Compress();
- fPeakMinusPedLowGain.Compress();
- fPeakMinusPedHighGain.Compress();
- fPeakMinusPedHighGainHisto.Compress();
-
- fDeadMap.Compress();
+ CompressAndSetOwner();
}
// assignment operator; use copy ctor to make life easy..
//_____________________________________________________________________
-AliCaloCalibPedestal& AliCaloCalibPedestal::operator = (const AliCaloCalibPedestal &source)
+AliCaloCalibPedestal& AliCaloCalibPedestal::operator = (AliCaloCalibPedestal &source)
{
// assignment operator; use copy ctor
if (&source == this) return *this;
//_____________________________________________________________________
void AliCaloCalibPedestal::Reset()
-{
- // Reset all arrays/histograms
+{ // Reset all arrays/histograms
+ ValidateProfiles(); // make sure histos/profiles exist
for (int i = 0; i < fModules; i++) {
GetPedProfileLowGain(i)->Reset();
GetPedProfileHighGain(i)->Reset();
+ GetPedLEDRefProfileLowGain(i)->Reset();
+ GetPedLEDRefProfileHighGain(i)->Reset();
GetPeakProfileLowGain(i)->Reset();
GetPeakProfileHighGain(i)->Reset();
GetPeakHighGainHisto(i)->Reset();
GetPedProfileLowGainDiff(i)->Reset();
GetPedProfileHighGainDiff(i)->Reset();
+ GetPedLEDRefProfileLowGainDiff(i)->Reset();
+ GetPedLEDRefProfileHighGainDiff(i)->Reset();
GetPeakProfileLowGainDiff(i)->Reset();
GetPeakProfileHighGainDiff(i)->Reset();
GetPedProfileLowGainRatio(i)->Reset();
GetPedProfileHighGainRatio(i)->Reset();
+ GetPedLEDRefProfileLowGainRatio(i)->Reset();
+ GetPedLEDRefProfileHighGainRatio(i)->Reset();
GetPeakProfileLowGainRatio(i)->Reset();
GetPeakProfileHighGainRatio(i)->Reset();
}
s >> keyValue;
// check stream status
- if( s.rdstate() & ios::failbit ) break;
+ if( ( s.rdstate() & ios::failbit ) == ios::failbit) break;
// skip rest of line if comments found
if( keyValue.substr( 0, 2 ) == "//" ) break;
}
//_____________________________________________________________________
-Bool_t AliCaloCalibPedestal::AddInfo(const AliCaloCalibPedestal *ped)
+Bool_t AliCaloCalibPedestal::AddInfo(AliCaloCalibPedestal *ped)
{
// just do this for the basic histograms/profiles that get filled in ProcessEvent
// may not have data for all modules, but let's just Add everything..
+ ValidateProfiles(); // make sure histos/profiles exist
+
for (int i = 0; i < fModules; i++) {
GetPedProfileLowGain(i)->Add( ped->GetPedProfileLowGain(i) );
GetPedProfileHighGain(i)->Add( ped->GetPedProfileHighGain(i) );
+ GetPedLEDRefProfileLowGain(i)->Add( ped->GetPedLEDRefProfileLowGain(i) );
+ GetPedLEDRefProfileHighGain(i)->Add( ped->GetPedLEDRefProfileHighGain(i) );
GetPeakProfileLowGain(i)->Add( ped->GetPeakProfileLowGain(i) );
GetPeakProfileHighGain(i)->Add( ped->GetPeakProfileHighGain(i) );
GetPeakHighGainHisto(i)->Add( ped->GetPeakHighGainHisto(i) );
-
}//end for nModules
+ // We should also copy other pieces of info: counters and parameters
+ // (not number of columns and rows etc which should be the same)
+ // note that I just assign them here rather than Add them, but we
+ // normally just Add (e.g. in Preprocessor) one object so this should be fine.
+ fNEvents = ped->GetNEvents();
+ fNChanFills = ped->GetNChanFills();
+ fDeadTowers = ped->GetDeadTowerCount();
+ fNewDeadTowers = ped->GetDeadTowerNew();
+ fResurrectedTowers = ped->GetDeadTowerResurrected();
+ fRunNumber = ped->GetRunNumber();
+ fSelectPedestalSamples = ped->GetSelectPedestalSamples();
+ fFirstPedestalSample = ped->GetFirstPedestalSample();
+ fLastPedestalSample = ped->GetLastPedestalSample();
+ fDeadThreshold = ped->GetDeadThreshold();
+ fWarningThreshold = ped->GetWarningThreshold();
+ fWarningFraction = ped->GetWarningFraction();
+ fHotSigma = ped->GetHotSigma();
+
// DeadMap; Diff profiles etc would need to be redone after this operation
return kTRUE;//We succesfully added info from the supplied object
{
// Method to process=analyze one event in the data stream
if (!in) return kFALSE; //Return right away if there's a null pointer
+ in->Reset(); // just in case the next customer forgets to check if the stream was reset..
+
fNEvents++; // one more event
+
+ if (fNEvents==1) ValidateProfiles(); // 1st event, make sure histos/profiles exist
// indices for the reading
int sample = 0;
- int gain = 0;
int time = 0;
int i = 0; // sample counter
int startBin = 0;
int max = AliEMCALGeoParams::fgkSampleMin, min = AliEMCALGeoParams::fgkSampleMax; // min and max sample values
int nsamples = 0;
- // for the pedestal calculation
- int sampleSum = 0; // sum of samples
- int squaredSampleSum = 0; // sum of samples squared
- int nSum = 0; // number of samples in sum
- // calc. quantities
- double mean = 0, squaredMean = 0, rms = 0;
-
+ // pedestal samples
+ int nPed = 0;
+ vector<int> pedSamples;
+
while (in->NextBunch()) {
const UShort_t *sig = in->GetSignals();
startBin = in->GetStartTimeBin();
// should we add it for the pedestal calculation?
if ( (fFirstPedestalSample<=time && time<=fLastPedestalSample) || // sample time in range
!fSelectPedestalSamples ) { // or we don't restrict the sample range.. - then we'll take all
- sampleSum += sample;
- squaredSampleSum += sample*sample;
- nSum++;
+ pedSamples.push_back( sig[i] );
+ nPed++;
}
} // loop over samples in bunch
if (nsamples > 0) { // this check is needed for when we have zero-supp. on, but not sparse readout
- // calculate pedesstal estimate: mean of possibly selected samples
- if (nSum > 0) {
- mean = sampleSum / (1.0 * nSum);
- squaredMean = squaredSampleSum / (1.0 * nSum);
- // The variance (rms squared) is equal to the mean of the squares minus the square of the mean..
- rms = sqrt(squaredMean - mean*mean);
- }
- else {
- mean = 0;
- squaredMean = 0;
- rms = 0;
- }
-
- // we're done with the calc. for this channel; let's prepare to fill histo
- gain = -1; // init to not valid value
- if ( in->IsLowGain() ) {
- gain = 0;
- }
- else if ( in->IsHighGain() ) {
- gain = 1;
- }
-
// it should be enough to check the SuperModule info for each DDL really, but let's keep it here for now
int arrayPos = in->GetModule(); //The modules are numbered starting from 0
if (arrayPos >= fModules) {
fNChanFills++; // one more channel found, and profile to be filled
//NOTE: coordinates are (column, row) for the profiles
- if (gain == 0) {
+ if ( in->IsLowGain() ) {
//fill the low gain histograms
((TProfile2D*)fPeakMinusPedLowGain[arrayPos])->Fill(in->GetColumn(), fRowMultiplier*in->GetRow(), max - min);
- if (nSum>0) { // only fill pedestal info in case it could be calculated
- ((TProfile2D*)fPedestalLowGain[arrayPos])->Fill(in->GetColumn(), fRowMultiplier*in->GetRow(), mean);
- ((TProfile2D*)fPedestalRMSLowGain[arrayPos])->Fill(in->GetColumn(), fRowMultiplier*in->GetRow(), rms);
+ if (nPed>0) { // only fill pedestal info in case it could be calculated
+ for ( i=0; i<nPed; i++) {
+ ((TProfile2D*)fPedestalLowGain[arrayPos])->Fill(in->GetColumn(), fRowMultiplier*in->GetRow(), pedSamples[i]);
+ }
}
}
- else if (gain == 1) {
+ else if ( in->IsHighGain() ) {
//fill the high gain ones
((TProfile2D*)fPeakMinusPedHighGain[arrayPos])->Fill(in->GetColumn(), fRowMultiplier*in->GetRow(), max - min);
- if (nSum>0) { // only fill pedestal info in case it could be calculated
- ((TProfile2D*)fPedestalHighGain[arrayPos])->Fill(in->GetColumn(), fRowMultiplier*in->GetRow(), mean);
- ((TProfile2D*)fPedestalRMSHighGain[arrayPos])->Fill(in->GetColumn(), fRowMultiplier*in->GetRow(), rms);
+ if (nPed>0) { // only fill pedestal info in case it could be calculated
+ for ( i=0; i<nPed; i++) {
+ ((TProfile2D*)fPedestalHighGain[arrayPos])->Fill(in->GetColumn(), fRowMultiplier*in->GetRow(), pedSamples[i]);
+ }
}
// for warning checks
int idx = in->GetRow() + fRows * in->GetColumn();
((TH2F*)fPeakMinusPedHighGainHisto[arrayPos])->Fill(idx, max - min);
- }//end if valid gain
+ }
+ else if ( in->IsLEDMonData() ) {
+ // for LED Mon data, the mapping class holds the gain info in the Row variable
+ // and the Strip number in the Column..
+ int gain = in->GetRow();
+ int stripId = in->GetColumn();
+ if (nPed>0 && stripId<fLEDRefs) {
+ if (gain == 0) {
+ for ( i=0; i<nPed; i++) {
+ ((TProfile*)fPedestalLEDRefLowGain[arrayPos])->Fill(stripId, pedSamples[i]);
+ }
+ }
+ else {
+ for ( i=0; i<nPed; i++) {
+ ((TProfile*)fPedestalLEDRefHighGain[arrayPos])->Fill(stripId, pedSamples[i]);
+ }
+ }
+ }
+ }
} // nsamples>0 check, some data found for this channel; not only trailer/header
}// end while over channel
}//end while over DDL's, of input stream
- in->Reset(); // just in case the next customer forgets to check if the stream was reset..
return kTRUE;
}
Bool_t AliCaloCalibPedestal::SaveHistograms(TString fileName, Bool_t saveEmptyHistos)
{
//Saves all the histograms (or profiles, to be accurate) to the designated file
-
+ ValidateProfiles(); // make sure histos/profiles exist
TFile destFile(fileName, "recreate");
if (destFile.IsZombie()) {
if( ((TProfile2D *)fPedestalHighGain[i])->GetEntries() || saveEmptyHistos) {
fPedestalHighGain[i]->Write();
}
- if( ((TProfile2D *)fPedestalRMSLowGain[i])->GetEntries() || saveEmptyHistos) {
- fPedestalRMSLowGain[i]->Write();
+ if( ((TProfile *)fPedestalLEDRefLowGain[i])->GetEntries() || saveEmptyHistos) {
+ fPedestalLEDRefLowGain[i]->Write();
}
- if( ((TProfile2D *)fPedestalRMSHighGain[i])->GetEntries() || saveEmptyHistos) {
- fPedestalRMSHighGain[i]->Write();
+ if( ((TProfile *)fPedestalLEDRefHighGain[i])->GetEntries() || saveEmptyHistos) {
+ fPedestalLEDRefHighGain[i]->Write();
}
if( ((TH2F *)fPeakMinusPedHighGainHisto[i])->GetEntries() || saveEmptyHistos) {
fPeakMinusPedHighGainHisto[i]->Write();
}
delete sourceFile;
-
+
//Reset the histogram ownership behaviour. NOTE: a better workaround would be good, since this may accidentally set AddDirectory to true, even
//if we are called by someone who has set it to false...
TH1::AddDirectory(kTRUE);
return kTRUE;//We succesfully loaded the object
}
+
+//_____________________________________________________________________
+Bool_t AliCaloCalibPedestal::SetReference(AliCaloCalibPedestal *ref)
+{ // set reference object
+ if (fReference) delete fReference;//Delete the reference object, if it already exists
+ fReference = 0;
+
+ fReference = ref;
+
+ if (!fReference || (fReference->GetDetectorType() != fDetType)) {
+ if (fReference) delete fReference;//Delete the object, in case we had an object of the wrong type
+ fReference = 0;
+ return kFALSE;
+ }
+
+ return kTRUE;//We succesfully loaded the object
+}
+
//_____________________________________________________________________
void AliCaloCalibPedestal::ValidateComparisonProfiles()
{
fColumns, 0.0, fColumns,
fRows, fRowMin, fRowMax,"s"));
+ //LED Ref/Mon pedestals, low gain
+ name = "hPedestalLEDReflowgainDiff";
+ name += i;
+ title = "Pedestal difference LEDRef, low gain, module ";
+ title += i;
+ fPedestalLEDRefLowGainDiff.Add(new TProfile(name, title,
+ fLEDRefs, 0.0, fLEDRefs, "s"));
+
+ //LED Ref/Mon pedestals, high gain
+ name = "hPedestalLEDRefhighgainDiff";
+ name += i;
+ title = "Pedestal difference LEDRef, high gain, module ";
+ title += i;
+ fPedestalLEDRefHighGainDiff.Add(new TProfile(name, title,
+ fLEDRefs, 0.0, fLEDRefs, "s"));
+
//Peak-Pedestals, high gain
name = "hPeakMinusPedhighgainDiff";
name += i;
fPeakMinusPedHighGainDiff.Add(new TProfile2D(name, title,
fColumns, 0.0, fColumns,
fRows, fRowMin, fRowMax,"s"));
+
+ //Peak-Pedestals, low gain
+ name = "hPeakMinusPedlowgainDiff";
+ name += i;
+ title = "Peak-Pedestal difference, low gain, module ";
+ title += i;
+ fPeakMinusPedLowGainDiff.Add(new TProfile2D(name, title,
+ fColumns, 0.0, fColumns,
+ fRows, fRowMin, fRowMax,"s"));
//Pedestals, low gain
name = "hPedlowgainRatio";
fPedestalHighGainRatio.Add(new TProfile2D(name, title,
fColumns, 0.0, fColumns,
fRows, fRowMin, fRowMax,"s"));
+
+ //LED Ref/Mon pedestals, low gain
+ name = "hPedestalLEDReflowgainRatio";
+ name += i;
+ title = "Pedestal ratio LEDRef, low gain, module ";
+ title += i;
+ fPedestalLEDRefLowGainRatio.Add(new TProfile(name, title,
+ fLEDRefs, 0.0, fLEDRefs, "s"));
+
+ //LED Ref/Mon pedestals, high gain
+ name = "hPedestalLEDRefhighgainRatio";
+ name += i;
+ title = "Pedestal ratio LEDRef, high gain, module ";
+ title += i;
+ fPedestalLEDRefHighGainRatio.Add(new TProfile(name, title,
+ fLEDRefs, 0.0, fLEDRefs, "s"));
//Peak-Pedestals, low gain
name = "hPeakMinusPedlowgainRatio";
//_____________________________________________________________________
void AliCaloCalibPedestal::ComputeDiffAndRatio()
-{
- // calculate differences and ratios relative to a reference
+{ // calculate differences and ratios relative to a reference
+ ValidateProfiles(); // make sure histos/profiles exist
ValidateComparisonProfiles();//Make sure the comparison histos exist
if (!fReference) {
return;//Return if the reference object isn't loaded
}
+ int bin = 0;
+ double diff = 0;
+ double ratio = 1;
for (int i = 0; i < fModules; i++) {
- //Compute the ratio of the histograms
-
- ((TProfile2D*)fPedestalLowGainRatio[i])->Divide(GetPedProfileLowGain(i), fReference->GetPedProfileLowGain(i), 1.0, 1.0);
- ((TProfile2D*)fPedestalHighGainRatio[i])->Divide(GetPedProfileHighGain(i), fReference->GetPedProfileHighGain(i), 1.0, 1.0);
- ((TProfile2D*)fPeakMinusPedLowGainRatio[i])->Divide(GetPeakProfileLowGain(i), fReference->GetPeakProfileLowGain(i), 1.0, 1.0);
- ((TProfile2D*)fPeakMinusPedHighGainRatio[i])->Divide(GetPeakProfileHighGain(i), fReference->GetPeakProfileHighGain(i), 1.0, 1.0);
-
//For computing the difference, we cannot simply do TProfile2D->Add(), because that subtracts the sum of all entries,
//which means that the mean of the new profile will not be the difference of the means. So do it by hand:
- for (int j = 0; j <= fColumns; j++) {
- for (int k = 0; k <= fRows; k++) {
- int bin = ((TProfile2D*)fPeakMinusPedHighGainDiff[i])->GetBin(j+1, k+1);//Note that we assume here that all histos have the same structure...
- double diff = fReference->GetPeakProfileHighGain(i)->GetBinContent(bin) - GetPeakProfileHighGain(i)->GetBinContent(bin);
- ((TProfile2D*)fPeakMinusPedHighGainDiff[i])->SetBinContent(j+1, k+1, diff);
- ((TProfile2D*)fPeakMinusPedHighGainDiff[i])->SetBinEntries(bin, 1);
-
- diff = fReference->GetPeakProfileLowGain(i)->GetBinContent(bin) - GetPeakProfileLowGain(i)->GetBinContent(bin);
- ((TProfile2D*)fPeakMinusPedLowGainDiff[i])->SetBinContent(j+1, k+1, diff);
- ((TProfile2D*)fPeakMinusPedLowGainDiff[i])->SetBinEntries(bin, 1);
-
- diff = fReference->GetPedProfileHighGain(i)->GetBinContent(bin) - GetPedProfileHighGain(i)->GetBinContent(bin);
- ((TProfile2D*)fPedestalHighGainDiff[i])->SetBinContent(j+1, k+1, diff);
- ((TProfile2D*)fPedestalHighGainDiff[i])->SetBinEntries(bin, 1);
-
- diff = fReference->GetPedProfileLowGain(i)->GetBinContent(bin) - GetPedProfileLowGain(i)->GetBinContent(bin);
- ((TProfile2D*)fPedestalLowGainDiff[i])->SetBinContent(j+1, k+1, diff);
- ((TProfile2D*)fPedestalLowGainDiff[i])->SetBinEntries(bin, 1);
-
+ for (int j = 0; j < fColumns; j++) {
+ for (int k = 0; k < fRows; k++) {
+ bin = ((TProfile2D*)fPeakMinusPedHighGainDiff[i])->GetBin(j+1, k+1);//Note that we assume here that all histos have the same structure...
+
+ if (fReference->GetPeakProfileHighGain(i)->GetBinContent(bin) > 0) {
+ diff = GetPeakProfileHighGain(i)->GetBinContent(bin) - fReference->GetPeakProfileHighGain(i)->GetBinContent(bin);
+ ((TProfile2D*)fPeakMinusPedHighGainDiff[i])->SetBinContent(bin, diff);
+ ((TProfile2D*)fPeakMinusPedHighGainDiff[i])->SetBinEntries(bin, 1);
+ ratio = GetPeakProfileHighGain(i)->GetBinContent(bin) / fReference->GetPeakProfileHighGain(i)->GetBinContent(bin);
+ ((TProfile2D*)fPeakMinusPedHighGainRatio[i])->SetBinContent(bin, ratio);
+ ((TProfile2D*)fPeakMinusPedHighGainRatio[i])->SetBinEntries(bin, 1);
+ }
+
+ if (fReference->GetPeakProfileLowGain(i)->GetBinContent(bin) > 0) {
+ diff = GetPeakProfileLowGain(i)->GetBinContent(bin) - fReference->GetPeakProfileLowGain(i)->GetBinContent(bin);
+ ((TProfile2D*)fPeakMinusPedLowGainDiff[i])->SetBinContent(bin, diff);
+ ((TProfile2D*)fPeakMinusPedLowGainDiff[i])->SetBinEntries(bin, 1);
+ ratio = GetPeakProfileLowGain(i)->GetBinContent(bin) / fReference->GetPeakProfileLowGain(i)->GetBinContent(bin);
+ ((TProfile2D*)fPeakMinusPedLowGainRatio[i])->SetBinContent(bin, ratio);
+ ((TProfile2D*)fPeakMinusPedLowGainRatio[i])->SetBinEntries(bin, 1);
+ }
+
+ if (fReference->GetPedProfileHighGain(i)->GetBinContent(bin) > 0) {
+ diff = GetPedProfileHighGain(i)->GetBinContent(bin) - fReference->GetPedProfileHighGain(i)->GetBinContent(bin);
+ ((TProfile2D*)fPedestalHighGainDiff[i])->SetBinContent(bin, diff);
+ ((TProfile2D*)fPedestalHighGainDiff[i])->SetBinEntries(bin, 1);
+ ratio = GetPedProfileHighGain(i)->GetBinContent(bin) / fReference->GetPedProfileHighGain(i)->GetBinContent(bin);
+ ((TProfile2D*)fPedestalHighGainRatio[i])->SetBinContent(bin, ratio);
+ ((TProfile2D*)fPedestalHighGainRatio[i])->SetBinEntries(bin, 1);
+ }
+
+ if (fReference->GetPedProfileLowGain(i)->GetBinContent(bin) > 0) {
+ diff = GetPedProfileLowGain(i)->GetBinContent(bin) - fReference->GetPedProfileLowGain(i)->GetBinContent(bin);
+ ((TProfile2D*)fPedestalLowGainDiff[i])->SetBinContent(bin, diff);
+ ((TProfile2D*)fPedestalLowGainDiff[i])->SetBinEntries(bin, 1);
+ ratio = GetPedProfileLowGain(i)->GetBinContent(bin) / fReference->GetPedProfileLowGain(i)->GetBinContent(bin);
+ ((TProfile2D*)fPedestalLowGainRatio[i])->SetBinContent(bin, ratio);
+ ((TProfile2D*)fPedestalLowGainRatio[i])->SetBinEntries(bin, 1);
+ }
+
} // rows
} // columns
-
+
+ // same for LED Ref/Mon channels
+ for (int j = 0; j <= fLEDRefs; j++) {
+ bin = j+1;//Note that we assume here that all histos have the same structure...
+
+ if (fReference->GetPedLEDRefProfileHighGain(i)->GetBinContent(bin) > 0) {
+ diff = GetPedLEDRefProfileHighGain(i)->GetBinContent(bin) - fReference->GetPedLEDRefProfileHighGain(i)->GetBinContent(bin);
+ ((TProfile*)fPedestalLEDRefHighGainDiff[i])->SetBinContent(bin, diff);
+ ((TProfile*)fPedestalLEDRefHighGainDiff[i])->SetBinEntries(bin, 1);
+ ratio = GetPedLEDRefProfileHighGain(i)->GetBinContent(bin) / fReference->GetPedLEDRefProfileHighGain(i)->GetBinContent(bin);
+ ((TProfile*)fPedestalLEDRefHighGainRatio[i])->SetBinContent(bin, ratio);
+ ((TProfile*)fPedestalLEDRefHighGainRatio[i])->SetBinEntries(bin, 1);
+ }
+
+ if (fReference->GetPedLEDRefProfileLowGain(i)->GetBinContent(bin) > 0) {
+ diff = GetPedLEDRefProfileLowGain(i)->GetBinContent(bin) - fReference->GetPedLEDRefProfileLowGain(i)->GetBinContent(bin);
+ ((TProfile*)fPedestalLEDRefLowGainDiff[i])->SetBinContent(bin, diff);
+ ((TProfile*)fPedestalLEDRefLowGainDiff[i])->SetBinEntries(bin, 1);
+ ratio = GetPedLEDRefProfileLowGain(i)->GetBinContent(bin) / fReference->GetPedLEDRefProfileLowGain(i)->GetBinContent(bin);
+ ((TProfile*)fPedestalLEDRefLowGainRatio[i])->SetBinContent(bin, ratio);
+ ((TProfile*)fPedestalLEDRefLowGainRatio[i])->SetBinEntries(bin, 1);
+ }
+
+ }
+
} // modules
}
//_____________________________________________________________________
void AliCaloCalibPedestal::ComputeHotAndWarningTowers(const char * hotMapFile)
{ // look for hot/noisy towers
+ ValidateProfiles(); // make sure histos/profiles exist
ofstream * fout = 0;
char name[512];//Quite a long temp buffer, just in case the filename includes a path
}
}
- //...and then we fit it
- hPeakFit->Fit("gaus", "OQ", "", hPeakFit->GetMean() - 3*hPeakFit->GetRMS(),
- hPeakFit->GetMean() + 3*hPeakFit->GetRMS());
+ //...and then we try to fit it
+ double mean = hPeakFit->GetMean();
+ double sigma = hPeakFit->GetRMS();
+ try {
+ hPeakFit->Fit("gaus", "OQ", "", mean - 3*sigma, mean + 3*sigma);
+ mean = hPeakFit->GetFunction("gaus")->GetParameter(1);
+ sigma = hPeakFit->GetFunction("gaus")->GetParameter(2);
+ }
+ catch (const std::exception & e) {
+ printf("AliCaloCalibPedestal: TH1D PeakFit exception %s", e.what());
+ }
//hPeakFit->Draw();
- double mean = hPeakFit->GetFunction("gaus")->GetParameter(1);
- double sigma = hPeakFit->GetFunction("gaus")->GetParameter(2);
delete hPeakFit;
//_____________________________________________________________________
void AliCaloCalibPedestal::ComputeDeadTowers(const char * deadMapFile)
{
+ ValidateProfiles(); // make sure histos/profiles exist
//Computes the number of dead towers etc etc into memory, after this you can call the GetDead... -functions
int countTot = 0;
int countNew = 0;
}
if (!diff->is_open()) {
delete diff;
- fout = 0;//Set the pointer to empty if the file was not opened
+ diff = 0;//Set the pointer to empty if the file was not opened
}
}
//_____________________________________________________________________
Bool_t AliCaloCalibPedestal::IsBadChannel(int imod, int icol, int irow) const
{
+ // Check if status info histo exists
+ if (!fDeadMap[imod]) {
+ return kFALSE;
+ }
+
//Check if channel is dead or hot.
Int_t status = (Int_t) ( ((TH2D*)fDeadMap[imod])->GetBinContent(icol,irow) );
if(status == kAlive)
//_____________________________________________________________________
void AliCaloCalibPedestal::SetChannelStatus(int imod, int icol, int irow, int status)
{
+ ValidateProfiles(); // make sure histos/profiles exist
//Set status of channel dead, hot, alive ...
((TH2D*)fDeadMap[imod])->SetBinContent(icol, irow, status);
}