// --- MUON header files ---
#include "AliMUONQADataMakerRec.h"
-#include "AliMUON2DMap.h"
-#include "AliMUONCluster.h"
-#include "AliMUONConstants.h"
-#include "AliMUONDDLTrigger.h"
-#include "AliMUONDarcHeader.h"
-#include "AliMUONDigitMaker.h"
-#include "AliMUONLocalStruct.h"
-#include "AliMUONLocalTrigger.h"
-#include "AliMUONRawStreamTracker.h"
-#include "AliMUONRawStreamTrigger.h"
-#include "AliMUONRegHeader.h"
-#include "AliMUONTrackerCalibratedDataMaker.h"
-#include "AliMUONTriggerDisplay.h"
-#include "AliMUONVCluster.h"
-#include "AliMUONVClusterStore.h"
-#include "AliMUONVDigit.h"
-#include "AliMUONVDigitStore.h"
-#include "AliMUONVTrackerData.h"
-#include "AliMUONVTriggerStore.h"
-#include "AliMUONTrack.h"
-#include "AliMUONTrackParam.h"
-#include "AliMUONESDInterface.h"
-#include "AliMUONCalibrationData.h"
-#include "AliMpBusPatch.h"
-#include "AliMpCDB.h"
-#include "AliMpConstants.h"
-#include "AliMpDDLStore.h"
-#include "AliMpDEIterator.h"
-#include "AliMpDEManager.h"
-#include "AliMpLocalBoard.h"
-#include "AliMpStationType.h"
-#include "AliMpTriggerCrate.h"
-#include "AliMpDCSNamer.h"
-#include "AliRawEventHeaderBase.h"
-
-// --- AliRoot header files ---
-#include "AliCDBManager.h"
-#include "AliCDBStorage.h"
-#include "AliESDEvent.h"
-#include "AliESDMuonTrack.h"
-#include "AliESDMuonCluster.h"
-#include "AliLog.h"
-#include "AliRawReader.h"
-#include "AliQAChecker.h"
-#include "AliCodeTimer.h"
-#include "AliDCSValue.h"
-
-// --- ROOT system ---
-#include <TClonesArray.h>
-#include <TFile.h>
-#include <TH1F.h>
-#include <TH1I.h>
-#include <TH2F.h>
-#include <TH3F.h>
-#include <Riostream.h>
-
//-----------------------------------------------------------------------------
/// \class AliMUONQADataMakerRec
///
/// MUON base class for quality assurance data (histo) maker
///
-/// \author C. Finck, D. Stocco, L. Aphecetche
+/// It is now only a steering class for the two subclasses AliMUONTrackerQADataMakerRec
+/// and AliMUONTriggerQADataMakerRec
+///
+/// \author C. Finck, D. Stocco, L. Aphecetche, A. Blanc
+
+#include "AliDAQ.h"
+#include "AliMUONTrackerQADataMakerRec.h"
+#include "AliMUONTriggerQADataMakerRec.h"
+#include "AliQAChecker.h"
+#include "AliRawReader.h"
+#include "AliRawEventHeaderBase.h"
/// \cond CLASSIMP
ClassImp(AliMUONQADataMakerRec)
/// \endcond
//____________________________________________________________________________
-AliMUONQADataMakerRec::AliMUONQADataMakerRec() :
-AliQADataMakerRec(AliQA::GetDetName(AliQA::kMUON), "MUON Quality Assurance Data Maker"),
-fIsInitRaws(kFALSE),
-fIsInitRecPointsTracker(kFALSE),
-fIsInitRecPointsTrigger(kFALSE),
-fIsInitESDs(kFALSE),
-fDigitStore(0x0),
-fTriggerStore(0x0),
-fDigitMaker(0x0),
-fClusterStore(0x0),
-fTrackerDataMaker(0x0),
-fhESDnTotClustersPerCh(0x0),
-fhESDnTotClustersPerDE(0x0),
-fhESDnTotFullClustersPerDE(0x0),
-fhESDSumClusterChargePerDE(0x0),
-fhESDSumClusterSizePerDE(0x0),
-fhESDSumResidualXPerDE(0x0),
-fhESDSumResidualYPerDE(0x0),
-fhESDSumResidualX2PerDE(0x0),
-fhESDSumResidualY2PerDE(0x0),
-fhESDSumLocalChi2XPerDE(0x0),
-fhESDSumLocalChi2YPerDE(0x0)
-{
- /// ctor
-
- AliDebug(1,"");
-
- Ctor();
-}
-
-//____________________________________________________________________________
-void
-AliMUONQADataMakerRec::Ctor()
+AliMUONQADataMakerRec::AliMUONQADataMakerRec(Bool_t tracker, Bool_t trigger) :
+AliQADataMakerRec(AliQAv1::GetDetName(AliQAv1::kMUON), "MUON Quality Assurance Data Maker"),
+fTracker(tracker ? new AliMUONTrackerQADataMakerRec(this) : 0x0),
+fTrigger(trigger ? new AliMUONTriggerQADataMakerRec(this) : 0x0)
{
- /// Init some members
- fDigitStore = AliMUONVDigitStore::Create("AliMUONDigitStoreV1");
- fDigitMaker = new AliMUONDigitMaker(kTRUE);
-}
-
-//____________________________________________________________________________
-AliMUONQADataMakerRec::AliMUONQADataMakerRec(const AliMUONQADataMakerRec& qadm) :
-AliQADataMakerRec(qadm),
-fIsInitRaws(kFALSE),
-fIsInitRecPointsTracker(kFALSE),
-fIsInitRecPointsTrigger(kFALSE),
-fIsInitESDs(kFALSE),
-fDigitStore(0x0),
-fTriggerStore(0x0),
-fDigitMaker(0x0),
-fClusterStore(0x0),
-fTrackerDataMaker(0x0),
-fhESDnTotClustersPerCh(0x0),
-fhESDnTotClustersPerDE(0x0),
-fhESDnTotFullClustersPerDE(0x0),
-fhESDSumClusterChargePerDE(0x0),
-fhESDSumClusterSizePerDE(0x0),
-fhESDSumResidualXPerDE(0x0),
-fhESDSumResidualYPerDE(0x0),
-fhESDSumResidualX2PerDE(0x0),
-fhESDSumResidualY2PerDE(0x0),
-fhESDSumLocalChi2XPerDE(0x0),
-fhESDSumLocalChi2YPerDE(0x0)
-{
- ///copy ctor
-
- AliDebug(1,"");
-
-
- SetName((const char*)qadm.GetName()) ;
- SetTitle((const char*)qadm.GetTitle());
-
- // Do not copy the digit store and digit maker, but create its own ones
-
- Ctor();
-
-}
-
-//__________________________________________________________________
-AliMUONQADataMakerRec& AliMUONQADataMakerRec::operator = (const AliMUONQADataMakerRec& qadm )
-{
- /// Assignment operator
-
- AliDebug(1,"");
-
- // check assignment to self
- if (this == &qadm) return *this;
-
- this->~AliMUONQADataMakerRec();
- new(this) AliMUONQADataMakerRec(qadm);
- return *this;
+ /// ctor
+ AliInfo(Form("tracker=%d trigger=%d",tracker,trigger));
}
//__________________________________________________________________
AliMUONQADataMakerRec::~AliMUONQADataMakerRec()
{
/// dtor
-
- AliDebug(1,"");
-
- AliCodeTimerAuto("");
-
- delete fDigitStore;
- delete fTriggerStore;
- delete fDigitMaker;
- delete fClusterStore;
- delete fTrackerDataMaker;
- delete fhESDnTotClustersPerCh;
- delete fhESDnTotClustersPerDE;
- delete fhESDnTotFullClustersPerDE;
- delete fhESDSumClusterChargePerDE;
- delete fhESDSumClusterSizePerDE;
- delete fhESDSumResidualXPerDE;
- delete fhESDSumResidualYPerDE;
- delete fhESDSumResidualX2PerDE;
- delete fhESDSumResidualY2PerDE;
- delete fhESDSumLocalChi2XPerDE;
- delete fhESDSumLocalChi2YPerDE;
+ delete fTracker;
+ delete fTrigger;
}
//____________________________________________________________________________
-void AliMUONQADataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray** list)
+Int_t AliMUONQADataMakerRec::Add2List(TH1 * hist, const Int_t index, AliQAv1::TASKINDEX_t task, const Bool_t expert, const Bool_t image, const Bool_t saveForCorr)
{
- ///Detector specific actions at end of cycle
-
- AliCodeTimerAuto("");
-
- // Display trigger histos in a more user friendly way
- DisplayTriggerInfo(task);
-
- for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
- SetEventSpecie(specie) ;
- if ( task == AliQA::kRAWS && fTrackerDataMaker )
- {
- TIter next(list[specie]);
- TObject* o;
- Bool_t alreadyThere(kFALSE);
- while ( ( o = next() ) && !alreadyThere )
- {
- TString classname(o->ClassName());
- if ( classname.Contains("TrackerData") ) alreadyThere = kTRUE;
- }
- if (!alreadyThere && fTrackerDataMaker)
- {
- AliInfo("Adding fTrackerDataMaker to the list of qa objects");
- list[specie]->AddAt(fTrackerDataMaker->Data(),(Int_t)kTrackerData);
- }
- if ( fTrackerDataMaker )
- {
- TH1* hbp = GetRawsData(kTrackerBusPatchOccupancy);
- hbp->Reset();
- TIter nextBP(AliMpDDLStore::Instance()->CreateBusPatchIterator());
- AliMpBusPatch* bp(0x0);
- AliMUONVTrackerData* data = fTrackerDataMaker->Data();
- Int_t occDim = 2;
-
- while ( ( bp = static_cast<AliMpBusPatch*>(nextBP())) )
- {
- Int_t busPatchId = bp->GetId();
- Int_t bin = hbp->FindBin(busPatchId);
- hbp->SetBinContent(bin,data->BusPatch(busPatchId,occDim));
- }
- }
- }
-
- // Normalize ESD histos
- if ( task == AliQA::kESDS ) {
-
- Double_t nTracks = GetESDsData(kESDnClustersPerTrack)->GetEntries();
- if (nTracks <= 0) continue;
-
- TH1* hESDnClustersPerCh = GetESDsData(kESDnClustersPerCh);
- TH1* hESDnClustersPerDE = GetESDsData(kESDnClustersPerDE);
- TH1* hESDClusterChargePerChMean = GetESDsData(kESDClusterChargePerChMean);
- TH1* hESDClusterChargePerChSigma = GetESDsData(kESDClusterChargePerChSigma);
- TH1* hESDClusterSizePerChMean = GetESDsData(kESDClusterSizePerChMean);
- TH1* hESDClusterSizePerChSigma = GetESDsData(kESDClusterSizePerChSigma);
- TH1* hESDResidualXPerChMean = GetESDsData(kESDResidualXPerChMean);
- TH1* hESDResidualXPerChSigma = GetESDsData(kESDResidualXPerChSigma);
- TH1* hESDResidualYPerChMean = GetESDsData(kESDResidualYPerChMean);
- TH1* hESDResidualYPerChSigma = GetESDsData(kESDResidualYPerChSigma);
- TH1* hESDLocalChi2XPerChMean = GetESDsData(kESDLocalChi2XPerChMean);
- TH1* hESDLocalChi2YPerChMean = GetESDsData(kESDLocalChi2YPerChMean);
- TH1* hESDClusterChargePerDE = GetESDsData(kESDClusterChargePerDE);
- TH1* hESDClusterSizePerDE = GetESDsData(kESDClusterSizePerDE);
- TH1* hESDResidualXPerDEMean = GetESDsData(kESDResidualXPerDEMean);
- TH1* hESDResidualXPerDESigma = GetESDsData(kESDResidualXPerDESigma);
- TH1* hESDResidualYPerDEMean = GetESDsData(kESDResidualYPerDEMean);
- TH1* hESDResidualYPerDESigma = GetESDsData(kESDResidualYPerDESigma);
- TH1* hESDLocalChi2XPerDEMean = GetESDsData(kESDLocalChi2XPerDEMean);
- TH1* hESDLocalChi2YPerDEMean = GetESDsData(kESDLocalChi2YPerDEMean);
-
- hESDnClustersPerCh->Reset();
- hESDnClustersPerDE->Reset();
- hESDnClustersPerCh->Add(fhESDnTotClustersPerCh, 1./nTracks);
- hESDnClustersPerDE->Add(fhESDnTotClustersPerDE, 1./nTracks);
-
- // loop over chambers
- for (Int_t iCh = 0; iCh < AliMUONConstants::NTrackingCh(); iCh++) {
-
- TH1* hESDClusterChargeInCh = GetESDsData(kESDClusterChargeInCh+iCh);
- Double_t sigmaCharge = hESDClusterChargeInCh->GetRMS();
- hESDClusterChargePerChMean->SetBinContent(iCh+1, hESDClusterChargeInCh->GetMean());
- hESDClusterChargePerChMean->SetBinError(iCh+1, hESDClusterChargeInCh->GetMeanError());
- hESDClusterChargePerChSigma->SetBinContent(iCh+1, sigmaCharge);
- hESDClusterChargePerChSigma->SetBinError(iCh+1, hESDClusterChargeInCh->GetRMSError());
-
- TH1* hESDClusterSizeInCh = GetESDsData(kESDClusterSizeInCh+iCh);
- Double_t sigmaSize = hESDClusterSizeInCh->GetRMS();
- hESDClusterSizePerChMean->SetBinContent(iCh+1, hESDClusterSizeInCh->GetMean());
- hESDClusterSizePerChMean->SetBinError(iCh+1, hESDClusterSizeInCh->GetMeanError());
- hESDClusterSizePerChSigma->SetBinContent(iCh+1, sigmaSize);
- hESDClusterSizePerChSigma->SetBinError(iCh+1, hESDClusterSizeInCh->GetRMSError());
-
- TH1* hESDResidualXInCh = GetESDsData(kESDResidualXInCh+iCh);
- Double_t sigmaResidualX = hESDResidualXInCh->GetRMS();
- hESDResidualXPerChMean->SetBinContent(iCh+1, hESDResidualXInCh->GetMean());
- hESDResidualXPerChMean->SetBinError(iCh+1, hESDResidualXInCh->GetMeanError());
- hESDResidualXPerChSigma->SetBinContent(iCh+1, sigmaResidualX);
- hESDResidualXPerChSigma->SetBinError(iCh+1, hESDResidualXInCh->GetRMSError());
-
- TH1* hESDResidualYInCh = GetESDsData(kESDResidualYInCh+iCh);
- Double_t sigmaResidualY = hESDResidualYInCh->GetRMS();
- hESDResidualYPerChMean->SetBinContent(iCh+1, hESDResidualYInCh->GetMean());
- hESDResidualYPerChMean->SetBinError(iCh+1, hESDResidualYInCh->GetMeanError());
- hESDResidualYPerChSigma->SetBinContent(iCh+1, sigmaResidualY);
- hESDResidualYPerChSigma->SetBinError(iCh+1, hESDResidualYInCh->GetRMSError());
-
- TH1* hESDLocalChi2XInCh = GetESDsData(kESDLocalChi2XInCh+iCh);
- Double_t sigmaLocalChi2X = hESDLocalChi2XInCh->GetRMS();
- hESDLocalChi2XPerChMean->SetBinContent(iCh+1, hESDLocalChi2XInCh->GetMean());
- hESDLocalChi2XPerChMean->SetBinError(iCh+1, hESDLocalChi2XInCh->GetMeanError());
-
- TH1* hESDLocalChi2YInCh = GetESDsData(kESDLocalChi2YInCh+iCh);
- Double_t sigmaLocalChi2Y = hESDLocalChi2YInCh->GetRMS();
- hESDLocalChi2YPerChMean->SetBinContent(iCh+1, hESDLocalChi2YInCh->GetMean());
- hESDLocalChi2YPerChMean->SetBinError(iCh+1, hESDLocalChi2YInCh->GetMeanError());
-
- // loop over DE into chamber iCh
- AliMpDEIterator it;
- it.First(iCh);
- while ( !it.IsDone()) {
-
- Int_t iDE = it.CurrentDEId();
-
- Double_t nClusters = fhESDnTotClustersPerDE->GetBinContent(iDE+1);
- if (nClusters > 1) {
-
- hESDClusterChargePerDE->SetBinContent(iDE+1, fhESDSumClusterChargePerDE->GetBinContent(iDE+1)/nClusters);
- hESDClusterChargePerDE->SetBinError(iDE+1, sigmaCharge/TMath::Sqrt(nClusters));
-
- Double_t meanResX = fhESDSumResidualXPerDE->GetBinContent(iDE+1)/nClusters;
- hESDResidualXPerDEMean->SetBinContent(iDE+1, meanResX);
- hESDResidualXPerDEMean->SetBinError(iDE+1, sigmaResidualX/TMath::Sqrt(nClusters));
- hESDResidualXPerDESigma->SetBinContent(iDE+1, TMath::Sqrt(fhESDSumResidualX2PerDE->GetBinContent(iDE+1)/nClusters - meanResX*meanResX));
- hESDResidualXPerDESigma->SetBinError(iDE+1, sigmaResidualX/TMath::Sqrt(2.*nClusters));
-
- Double_t meanResY = fhESDSumResidualYPerDE->GetBinContent(iDE+1)/nClusters;
- hESDResidualYPerDEMean->SetBinContent(iDE+1, meanResY);
- hESDResidualYPerDEMean->SetBinError(iDE+1, sigmaResidualY/TMath::Sqrt(nClusters));
- hESDResidualYPerDESigma->SetBinContent(iDE+1, TMath::Sqrt(fhESDSumResidualY2PerDE->GetBinContent(iDE+1)/nClusters - meanResY*meanResY));
- hESDResidualYPerDESigma->SetBinError(iDE+1, sigmaResidualY/TMath::Sqrt(2.*nClusters));
-
- hESDLocalChi2XPerDEMean->SetBinContent(iDE+1, fhESDSumLocalChi2XPerDE->GetBinContent(iDE+1)/nClusters);
- hESDLocalChi2XPerDEMean->SetBinError(iDE+1, sigmaLocalChi2X/TMath::Sqrt(nClusters));
-
- hESDLocalChi2YPerDEMean->SetBinContent(iDE+1, fhESDSumLocalChi2YPerDE->GetBinContent(iDE+1)/nClusters);
- hESDLocalChi2YPerDEMean->SetBinError(iDE+1, sigmaLocalChi2Y/TMath::Sqrt(nClusters));
-
- } else {
-
- hESDClusterChargePerDE->SetBinContent(iDE+1, fhESDSumClusterChargePerDE->GetBinContent(iDE+1));
- hESDClusterChargePerDE->SetBinError(iDE+1, hESDClusterChargeInCh->GetXaxis()->GetXmax());
-
- hESDResidualXPerDEMean->SetBinContent(iDE+1, fhESDSumResidualXPerDE->GetBinContent(iDE+1));
- hESDResidualXPerDEMean->SetBinError(iDE+1, hESDResidualXInCh->GetXaxis()->GetXmax());
- hESDResidualXPerDESigma->SetBinContent(iDE+1, 0.);
- hESDResidualXPerDESigma->SetBinError(iDE+1, hESDResidualXInCh->GetXaxis()->GetXmax());
-
- hESDResidualYPerDEMean->SetBinContent(iDE+1, fhESDSumResidualYPerDE->GetBinContent(iDE+1));
- hESDResidualYPerDEMean->SetBinError(iDE+1, hESDResidualYInCh->GetXaxis()->GetXmax());
- hESDResidualYPerDESigma->SetBinContent(iDE+1, 0.);
- hESDResidualYPerDESigma->SetBinError(iDE+1, hESDResidualYInCh->GetXaxis()->GetXmax());
-
- hESDLocalChi2XPerDEMean->SetBinContent(iDE+1, fhESDSumLocalChi2XPerDE->GetBinContent(iDE+1));
- hESDLocalChi2XPerDEMean->SetBinError(iDE+1, hESDLocalChi2XInCh->GetXaxis()->GetXmax());
-
- hESDLocalChi2YPerDEMean->SetBinContent(iDE+1, fhESDSumLocalChi2YPerDE->GetBinContent(iDE+1));
- hESDLocalChi2YPerDEMean->SetBinError(iDE+1, hESDLocalChi2YInCh->GetXaxis()->GetXmax());
-
- }
-
- Double_t nFullClusters = fhESDnTotFullClustersPerDE->GetBinContent(iDE+1);
- if (nFullClusters > 1) {
-
- hESDClusterSizePerDE->SetBinContent(iDE+1, fhESDSumClusterSizePerDE->GetBinContent(iDE+1)/nFullClusters);
- hESDClusterSizePerDE->SetBinError(iDE+1, sigmaSize/TMath::Sqrt(nFullClusters));
-
- } else {
-
- hESDClusterSizePerDE->SetBinContent(iDE+1, fhESDSumClusterSizePerDE->GetBinContent(iDE+1));
- hESDClusterSizePerDE->SetBinError(iDE+1, hESDClusterSizeInCh->GetXaxis()->GetXmax());
-
- }
-
- it.Next();
- }
-
- }
-
- }
-
+ TObjArray** list = GetList(task);
+ if (list)
+ {
+ return Add2List(hist,index,list,expert,image,saveForCorr);
}
-
- // do the QA checking
- AliQAChecker::Instance()->Run(AliQA::kMUON, task, list) ;
+ return -1;
}
//____________________________________________________________________________
-void AliMUONQADataMakerRec::InitRaws()
+void AliMUONQADataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
{
- /// create Raws histograms in Raws subdir
-
- AliCodeTimerAuto("");
-
- Bool_t forExpert(kTRUE);
-
- if ( ! AliCDBManager::Instance()->GetDefaultStorage() )
- {
- AliError("CDB default storage not set. Cannot work.");
- fIsInitRaws=kFALSE;
- }
-
- TH3F* h3 = new TH3F("hTriggerScalersBendPlane", "Trigger scalers in bending plane",
- 4, 10.5, 14.5,
- 234, 0.5, 234.5,
- 16, -0.5, 15.5);
- h3->GetXaxis()->SetTitle("Chamber");
- h3->GetYaxis()->SetTitle("Board");
- h3->GetZaxis()->SetTitle("Strip");
- Add2RawsList(h3, kTriggerScalersBP,forExpert);
-
- TH3F* h4 = new TH3F("hTriggerScalersNonBendPlane", "Trigger scalers in non-bending plane",
- 4, 10.5, 14.5,
- 234, 0.5, 234.5,
- 16, -0.5, 15.5);
- h4->GetXaxis()->SetTitle("Chamber");
- h4->GetYaxis()->SetTitle("Board");
- h4->GetZaxis()->SetTitle("Strip");
- Add2RawsList(h4, kTriggerScalersNBP,forExpert);
-
- AliMUONTriggerDisplay triggerDisplay;
- TString histoName, histoTitle;
- for(Int_t iCath=0; iCath<AliMpConstants::NofCathodes(); iCath++){
- TString cathName = ( iCath==0 ) ? "BendPlane" : "NonBendPlane";
- for(Int_t iChamber=0; iChamber<AliMpConstants::NofTriggerChambers(); iChamber++){
- histoName = Form("hScalers%sChamber%i", cathName.Data(), 11+iChamber);
- histoTitle = Form("Chamber %i: Scalers %s", 11+iChamber, cathName.Data());
- TH2F* h5 = (TH2F*)triggerDisplay.GetEmptyDisplayHisto(histoName, AliMUONTriggerDisplay::kDisplayStrips,
- iCath, iChamber, histoTitle);
- Add2RawsList(h5, kTriggerScalersDisplay + AliMpConstants::NofTriggerChambers()*iCath + iChamber,!forExpert);
- }
- }
-
- TH2F* h6 = new TH2F("hTriggerRPCI", "Trigger RPC currents",
- 4, 10.5, 14.5,
- 18, -0.5, 17.5);
- h6->GetXaxis()->SetTitle("Chamber");
- h6->GetYaxis()->SetTitle("RPC");
- Add2RawsList(h6, kTriggerRPCi, forExpert);
-
- TH2F* h7 = new TH2F("hTriggerRPCHV", "Trigger RPC HV",
- 4, 10.5, 14.5,
- 18, -0.5, 17.5);
- h7->GetXaxis()->SetTitle("Chamber");
- h7->GetYaxis()->SetTitle("RPC");
- Add2RawsList(h7, kTriggerRPChv, forExpert);
-
- for(Int_t iChamber=0; iChamber<AliMpConstants::NofTriggerChambers(); iChamber++){
- histoName = Form("hRPCIChamber%i", 11+iChamber);
- histoTitle = Form("Chamber %i: RPC Currents (#muA)", 11+iChamber);
- TH2F* h8 = (TH2F*)triggerDisplay.GetEmptyDisplayHisto(histoName, AliMUONTriggerDisplay::kDisplaySlats,
- 0, iChamber, histoTitle);
- Add2RawsList(h8, kTriggerIDisplay + iChamber, !forExpert);
- }
-
- for(Int_t iChamber=0; iChamber<AliMpConstants::NofTriggerChambers(); iChamber++){
- histoName = Form("hRPCHVChamber%i", 11+iChamber);
- histoTitle = Form("Chamber %i: RPC HV (V)", 11+iChamber);
- TH2F* h9 = (TH2F*)triggerDisplay.GetEmptyDisplayHisto(histoName, AliMUONTriggerDisplay::kDisplaySlats,
- 0, iChamber, histoTitle);
- Add2RawsList(h9, kTriggerHVDisplay + iChamber, !forExpert);
- }
-
-
- Int_t nbp(0);
- TIter next(AliMpDDLStore::Instance()->CreateBusPatchIterator());
- while (next())
+ /// Detector specific actions at end of cycle
+ //
+ ResetEventTrigClasses(); // RS
+ //
+ for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
{
- ++nbp;
- }
-
- TH1* hbp = new TH1F("hTrackerBusPatchOccupancy","Occupancy of bus patches",
- nbp,-0.5,nbp-0.5);
-
- Add2RawsList(hbp,kTrackerBusPatchOccupancy,!forExpert);
-
- const Bool_t histogram(kFALSE);
- const Bool_t fastDecoder(kTRUE);
-
- fTrackerDataMaker = new AliMUONTrackerCalibratedDataMaker(GetRecoParam(),
- AliCDBManager::Instance()->GetRun(),
- 0x0,
- AliCDBManager::Instance()->GetDefaultStorage()->GetURI(),
- "NOGAIN",
- histogram,
- 0.0,0.0,
- fastDecoder);
-
- fTrackerDataMaker->Data()->DisableChannelLevel(); // to save up disk space, we only store starting at the manu level
-
- fTrackerDataMaker->SetRunning(kTRUE);
-
- fIsInitRaws = kTRUE;
+ if (! IsValidEventSpecie(specie, list) ) continue;
+
+ SetEventSpecie(AliRecoParam::ConvertIndex(specie)); // needed by the GetXXXData methods
+
+ if ( task == AliQAv1::kRAWS )
+ {
+ if ( fTracker ) fTracker->EndOfDetectorCycleRaws(specie,list);
+ if ( fTrigger ) fTrigger->EndOfDetectorCycleRaws(specie,list);
+ }
+ else if ( task == AliQAv1::kRECPOINTS )
+ {
+ // normalize recpoints histograms
+ if ( fTracker ) fTracker->EndOfDetectorCycleRecPoints(specie,list);
+ if ( fTrigger ) fTrigger->EndOfDetectorCycleRecPoints(specie,list);
+ }
+ else if ( task == AliQAv1::kESDS )
+ {
+ // normalize esds histograms
+ if ( fTracker ) fTracker->EndOfDetectorCycleESDs(specie,list);
+ if ( fTrigger ) fTrigger->EndOfDetectorCycleESDs(specie,list);
+ }
+ else if ( task == AliQAv1::kDIGITSR )
+ {
+ if ( fTracker ) fTracker->EndOfDetectorCycleDigits(specie,list);
+ if ( fTrigger ) fTrigger->EndOfDetectorCycleDigits(specie,list);
+ }
+ else
+ {
+ AliFatal(Form("Not implemented for task %s",AliQAv1::GetTaskName(task).Data()));
+ }
+ } // loop on specie
+
+ // do the QA checking
+ AliQAChecker::Instance()->Run(AliQAv1::kMUON,task,list,const_cast<AliDetectorRecoParam*>(GetRecoParam()));
}
//____________________________________________________________________________
-void AliMUONQADataMakerRec::InitRecPoints()
+TObject* AliMUONQADataMakerRec::GetData(AliQAv1::TASKINDEX_t task, const Int_t index)
{
- /// create Reconstructed Points histograms in RecPoints subdir
-
- AliCodeTimerAuto("");
-
- InitRecPointsTrigger();
- InitRecPointsTracker();
+ TObjArray** list = GetList(task);
+ if (list) return GetData(list,index);
+ return 0x0;
}
//____________________________________________________________________________
-void AliMUONQADataMakerRec::InitRecPointsTracker()
+TObjArray** AliMUONQADataMakerRec::GetList(AliQAv1::TASKINDEX_t task)
{
- /// create Reconstructed Points histograms in RecPoints subdir for the
- /// MUON tracker subsystem.
-
- AliCodeTimerAuto("");
-
- Bool_t forExpert(kTRUE);
-
- AliMpDEIterator it;
-
- it.First();
-
- Int_t ndes(0);
-
- while ( !it.IsDone())
- {
- Int_t detElemId = it.CurrentDEId();
-
- it.Next();
-
- if ( AliMpDEManager::GetStationType(detElemId) != AliMp::kStationTrigger )
- {
- ndes = TMath::Max(ndes,detElemId);
-
- TH1* h = new TH1I(Form("hTrackerClusterMultiplicityForDE%04d",detElemId),
- Form("Multiplicity of the clusters in detection element %d",detElemId),
- 100,0,100);
-
- h->GetXaxis()->SetTitle("Detection Element Id");
-
- Add2RecPointsList(h,kTrackerClusterMultiplicityPerDE+detElemId,forExpert);
-
- h = new TH1I(Form("hTrackerClusterChargeForDE%04d",detElemId),
- Form("Charge of the clusters in detection element %d",detElemId),
- 100,0,1000);
-
- h->GetXaxis()->SetTitle("Detection Element Id");
-
- Add2RecPointsList(h,kTrackerClusterChargePerDE+detElemId,forExpert);
-
- }
-
- }
-
- TH1* h = new TH1I("hTrackerNumberOfClustersPerDE","Number of clusters per detection element",
- ndes, 0.5, ndes + 0.5);
-
- h->GetXaxis()->SetTitle("Detection Element Id");
-
- Add2RecPointsList(h, kTrackerNumberOfClustersPerDE,!forExpert);
-
- for ( Int_t i = 0; i < AliMpConstants::NofTrackingChambers(); ++i )
- {
- TH1* h1 = new TH1I("hTrackerNumberOfClustersPerChamber","Number of clusters per chamber",
- AliMpConstants::NofTrackingChambers(),-0.5,AliMpConstants::NofTrackingChambers()-0.5);
- Add2RecPointsList(h1,kTrackerNumberOfClustersPerChamber,forExpert);
- h1 = new TH1I(Form("hTrackerClusterMultiplicityForChamber%d",i),
- Form("Cluster multiplicity for chamber %d",i),
- 100,0,100);
- Add2RecPointsList(h1,kTrackerClusterMultiplicityPerChamber+i,forExpert);
- h1 = new TH1I(Form("hTrackerClusterChargeForChamber%d",i),
- Form("Cluster charge for chamber %d",i),
- 100,0,1000);
- Add2RecPointsList(h1,kTrackerClusterChargePerChamber+i,forExpert);
- }
-
- fIsInitRecPointsTracker=kTRUE;
+ // enum TASKINDEX_t {
+ // kNULLTASKINDEX=-1, kRAWS, kHITS, kSDIGITS, kDIGITS, kDIGITSR, kRECPOINTS, kTRACKSEGMENTS, kRECPARTICLES, kESDS, kNTASKINDEX };
+ if ( task == AliQAv1::kRAWS )
+ {
+ return fRawsQAList;
+ }
+ else if ( task == AliQAv1::kDIGITS || task == AliQAv1::kDIGITSR )
+ {
+ return fDigitsQAList;
+ }
+ else if ( task == AliQAv1::kRECPOINTS )
+ {
+ return fRecPointsQAList;
+ }
+ else
+ {
+ AliFatal(Form("task %s not supported here yet",AliQAv1::GetTaskName(task).Data()));
+ }
+ return 0x0;
}
//____________________________________________________________________________
-void AliMUONQADataMakerRec::InitRecPointsTrigger()
+void AliMUONQADataMakerRec::InitRaws()
{
- /// create Reconstructed Points histograms in RecPoints subdir for the
- /// MUON Trigger subsystem.
+ /// create Raws histograms in Raws subdir
- Bool_t forExpert(kTRUE);
-
- TH3F* h0 = new TH3F("hTriggerDigitsBendPlane", "Trigger digits in bending plane",
- 4, 10.5, 14.5,
- 234, 0.5, 234.5,
- 16, -0.5, 15.5);
- h0->GetXaxis()->SetTitle("Chamber");
- h0->GetYaxis()->SetTitle("Board");
- h0->GetZaxis()->SetTitle("Strip");
- Add2RecPointsList(h0, kTriggerDigitsBendPlane,forExpert);
-
- TH3F* h1 = new TH3F("hTriggerDigitsNonBendPlane", "Trigger digits in non-bending plane",
- 4, 10.5, 14.5,
- 234, 0.5, 234.5,
- 16, -0.5, 15.5);
- h1->GetXaxis()->SetTitle("Chamber");
- h1->GetYaxis()->SetTitle("Board");
- h1->GetZaxis()->SetTitle("Strip");
- Add2RecPointsList(h1, kTriggerDigitsNonBendPlane,forExpert);
-
- TH1F* h2 = new TH1F("hTriggeredBoards", "Triggered boards", 234, 0.5, 234.5);
- Add2RecPointsList(h2, kTriggeredBoards,forExpert);
+ if ( fTracker ) fTracker->InitRaws();
+ if ( fTrigger ) fTrigger->InitRaws();
+ //
+ ClonePerTrigClass(AliQAv1::kRAWS); // this should be the last line
+}
- AliMUONTriggerDisplay triggerDisplay;
- TString histoName, histoTitle;
- for(Int_t iCath=0; iCath<AliMpConstants::NofCathodes(); iCath++){
- TString cathName = ( iCath==0 ) ? "BendPlane" : "NonBendPlane";
- for(Int_t iChamber=0; iChamber<AliMpConstants::NofTriggerChambers(); iChamber++){
- histoName = Form("hTriggerDigits%sChamber%i", cathName.Data(), 11+iChamber);
- histoTitle = Form("Chamber %i: Fired pads %s", 11+iChamber, cathName.Data());
- TH2F* h3 = (TH2F*)triggerDisplay.GetEmptyDisplayHisto(histoName, AliMUONTriggerDisplay::kDisplayStrips,
- iCath, iChamber, histoTitle);
- Add2RecPointsList(h3, kTriggerDigitsDisplay + AliMpConstants::NofTriggerChambers()*iCath + iChamber,!forExpert);
- }
- }
+//__________________________________________________________________
+void AliMUONQADataMakerRec::InitDigits()
+{
+ /// Initialized Digits spectra
+ if ( fTracker ) fTracker->InitDigits();
+ if ( fTrigger ) fTrigger->InitDigits();
+ //
+ ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
+}
- TH2F* h4 = (TH2F*)triggerDisplay.GetEmptyDisplayHisto("hFiredBoardsDisplay", AliMUONTriggerDisplay::kDisplayBoards,
- 0, 0, "Fired boards");
- Add2RecPointsList(h4, kTriggerBoardsDisplay,!forExpert);
-
- fIsInitRecPointsTrigger = kTRUE;
+//____________________________________________________________________________
+void AliMUONQADataMakerRec::InitRecPoints()
+{
+ /// create Reconstructed Points histograms in RecPoints subdir
+ if ( fTracker ) fTracker->InitRecPoints();
+ if ( fTrigger ) fTrigger->InitRecPoints();
+ //
+ ClonePerTrigClass(AliQAv1::kRECPOINTS); // this should be the last line
}
void AliMUONQADataMakerRec::InitESDs()
{
///create ESDs histograms in ESDs subdir
-
- Bool_t forExpert(kTRUE);
-
- Int_t nCh = AliMUONConstants::NTrackingCh();
- Int_t nDE = 1100;
-
- // track info
- TH1F* hESDnTracks = new TH1F("hESDnTracks", "number of tracks;n_{tracks}", 20, 0., 20.);
- Add2ESDsList(hESDnTracks, kESDnTracks,!forExpert);
-
- TH1F* hESDMatchTrig = new TH1F("hESDMatchTrig", "number of tracks matched with trigger;n_{tracks}", 20, 0., 20.);
- Add2ESDsList(hESDMatchTrig, kESDMatchTrig,!forExpert);
-
- TH1F* hESDMomentum = new TH1F("hESDMomentum", "momentum distribution;p (GeV/c)", 300, 0., 300);
- Add2ESDsList(hESDMomentum, kESDMomentum,forExpert);
-
- TH1F* hESDPt = new TH1F("hESDPt", "transverse momentum distribution;p_{t} (GeV/c)", 200, 0., 50);
- Add2ESDsList(hESDPt, kESDPt,forExpert);
-
- TH1F* hESDRapidity = new TH1F("hESDRapidity", "rapidity distribution;rapidity", 200, -4.5, -2.);
- Add2ESDsList(hESDRapidity, kESDRapidity,forExpert);
-
- TH1F* hESDChi2 = new TH1F("hESDChi2", "normalized #chi^{2} distribution;#chi^{2} / ndf", 500, 0., 50.);
- Add2ESDsList(hESDChi2, kESDChi2,forExpert);
-
- TH1F* hESDProbChi2 = new TH1F("hESDProbChi2", "distribution of probability of #chi^{2};prob(#chi^{2})", 100, 0., 1.);
- Add2ESDsList(hESDProbChi2, kESDProbChi2,forExpert);
-
- // cluster info
- for (Int_t i = 0; i < nCh; i++) {
- Float_t rMax = AliMUONConstants::Rmax(i/2);
- TH2F* hESDClusterHitMap = new TH2F(Form("hESDClusterHitMap%d",i+1), Form("cluster position distribution in chamber %d;X (cm);Y (cm)",i+1),
- 100, -rMax, rMax, 100, -rMax, rMax);
- Add2ESDsList(hESDClusterHitMap, kESDClusterHitMap+i,forExpert);
- }
-
- TH1F* hESDnClustersPerTrack = new TH1F("hESDnClustersPerTrack", "number of associated clusters per track;n_{clusters}", 20, 0., 20.);
- Add2ESDsList(hESDnClustersPerTrack, kESDnClustersPerTrack,!forExpert);
-
- TH1F* hESDnClustersPerCh = new TH1F("hESDnClustersPerCh", "averaged number of clusters per chamber per track;chamber ID;<n_{clusters}>", nCh, 0, nCh);
- hESDnClustersPerCh->SetFillColor(kRed);
- Add2ESDsList(hESDnClustersPerCh, kESDnClustersPerCh,forExpert);
-
- TH1F* hESDnClustersPerDE = new TH1F("hESDnClustersPerDE", "averaged number of clusters per DE per track;DetElem ID;<n_{clusters}>", nDE, 0, nDE);
- hESDnClustersPerDE->SetFillColor(kRed);
- Add2ESDsList(hESDnClustersPerDE, kESDnClustersPerDE,forExpert);
-
- for (Int_t i = 0; i < nCh; i++) {
- TH1F* hESDClusterChargeInCh = new TH1F(Form("hESDClusterChargeInCh%d",i+1), Form("cluster charge distribution in chamber %d;charge (ADC counts)",i+1), 500, 0., 5000.);
- Add2ESDsList(hESDClusterChargeInCh, kESDClusterChargeInCh+i,forExpert);
- }
-
- TH1F* hESDClusterChargePerChMean = new TH1F("hESDClusterChargePerChMean", "cluster mean charge per chamber;chamber ID;<charge> (ADC counts)", nCh, 0, nCh);
- hESDClusterChargePerChMean->SetOption("P");
- hESDClusterChargePerChMean->SetMarkerStyle(kFullDotMedium);
- hESDClusterChargePerChMean->SetMarkerColor(kRed);
- Add2ESDsList(hESDClusterChargePerChMean, kESDClusterChargePerChMean,forExpert);
-
- TH1F* hESDClusterChargePerChSigma = new TH1F("hESDClusterChargePerChSigma", "cluster charge dispersion per chamber;chamber ID;#sigma_{charge} (ADC counts)", nCh, 0, nCh);
- hESDClusterChargePerChSigma->SetOption("P");
- hESDClusterChargePerChSigma->SetMarkerStyle(kFullDotMedium);
- hESDClusterChargePerChSigma->SetMarkerColor(kRed);
- Add2ESDsList(hESDClusterChargePerChSigma, kESDClusterChargePerChSigma,forExpert);
-
- TH1F* hESDClusterChargePerDE = new TH1F("hESDClusterChargePerDE", "cluster mean charge per DE;DetElem ID;<charge> (ADC counts)", nDE, 0, nDE);
- hESDClusterChargePerDE->SetOption("P");
- hESDClusterChargePerDE->SetMarkerStyle(kFullDotMedium);
- hESDClusterChargePerDE->SetMarkerColor(kRed);
- Add2ESDsList(hESDClusterChargePerDE, kESDClusterChargePerDE,forExpert);
-
- for (Int_t i = 0; i < nCh; i++) {
- TH1F* hESDClusterSizeInCh = new TH1F(Form("hESDClusterSizeInCh%d",i+1), Form("cluster size distribution in chamber %d;size (n_{pads})",i+1), 200, 0., 200.);
- Add2ESDsList(hESDClusterSizeInCh, kESDClusterSizeInCh+i,forExpert);
- }
-
- TH1F* hESDClusterSizePerChMean = new TH1F("hESDClusterSizePerChMean", "cluster mean size per chamber;chamber ID;<size> (n_{pads})", nCh, 0, nCh);
- hESDClusterSizePerChMean->SetOption("P");
- hESDClusterSizePerChMean->SetMarkerStyle(kFullDotMedium);
- hESDClusterSizePerChMean->SetMarkerColor(kRed);
- Add2ESDsList(hESDClusterSizePerChMean, kESDClusterSizePerChMean,forExpert);
-
- TH1F* hESDClusterSizePerChSigma = new TH1F("hESDClusterSizePerChSigma", "cluster size dispersion per chamber;chamber ID;#sigma_{size} (n_{pads})", nCh, 0, nCh);
- hESDClusterSizePerChSigma->SetOption("P");
- hESDClusterSizePerChSigma->SetMarkerStyle(kFullDotMedium);
- hESDClusterSizePerChSigma->SetMarkerColor(kRed);
- Add2ESDsList(hESDClusterSizePerChSigma, kESDClusterSizePerChSigma,forExpert);
-
- TH1F* hESDClusterSizePerDE = new TH1F("hESDClusterSizePerDE", "cluster mean size per DE;DetElem ID;<size> (n_{pads})", nDE, 0, nDE);
- hESDClusterSizePerDE->SetOption("P");
- hESDClusterSizePerDE->SetMarkerStyle(kFullDotMedium);
- hESDClusterSizePerDE->SetMarkerColor(kRed);
- Add2ESDsList(hESDClusterSizePerDE, kESDClusterSizePerDE,forExpert);
-
- // cluster - track info
- for (Int_t i = 0; i < nCh; i++) {
- TH1F* hESDResidualXInCh = new TH1F(Form("hESDResidualXInCh%d",i+1), Form("cluster-track residual-X distribution in chamber %d;#Delta_{X} (cm)",i+1), 1000, -5., 5.);
- Add2ESDsList(hESDResidualXInCh, kESDResidualXInCh+i,forExpert);
-
- TH1F* hESDResidualYInCh = new TH1F(Form("hESDResidualYInCh%d",i+1), Form("cluster-track residual-Y distribution in chamber %d;#Delta_{Y} (cm)",i+1), 1000, -1., 1.);
- Add2ESDsList(hESDResidualYInCh, kESDResidualYInCh+i,forExpert);
-
- TH1F* hESDLocalChi2XInCh = new TH1F(Form("hESDLocalChi2XInCh%d",i+1), Form("local chi2-X distribution in chamber %d;local #chi^{2}_{X}",i+1), 1000, 0., 25);
- Add2ESDsList(hESDLocalChi2XInCh, kESDLocalChi2XInCh+i,forExpert);
-
- TH1F* hESDLocalChi2YInCh = new TH1F(Form("hESDLocalChi2YInCh%d",i+1), Form("local chi2-Y distribution in chamber %d;local #chi^{2}_{Y}",i+1), 1000, 0., 25);
- Add2ESDsList(hESDLocalChi2YInCh, kESDLocalChi2YInCh+i,forExpert);
- }
-
- TH1F* hESDResidualXPerChMean = new TH1F("hESDResidualXPerChMean", "cluster-track residual-X per Ch: mean;chamber ID;<#Delta_{X}> (cm)", nCh, 0, nCh);
- hESDResidualXPerChMean->SetOption("P");
- hESDResidualXPerChMean->SetMarkerStyle(kFullDotMedium);
- hESDResidualXPerChMean->SetMarkerColor(kRed);
- Add2ESDsList(hESDResidualXPerChMean, kESDResidualXPerChMean,forExpert);
-
- TH1F* hESDResidualYPerChMean = new TH1F("hESDResidualYPerChMean", "cluster-track residual-Y per Ch: mean;chamber ID;<#Delta_{Y}> (cm)", nCh, 0, nCh);
- hESDResidualYPerChMean->SetOption("P");
- hESDResidualYPerChMean->SetMarkerStyle(kFullDotMedium);
- hESDResidualYPerChMean->SetMarkerColor(kRed);
- Add2ESDsList(hESDResidualYPerChMean, kESDResidualYPerChMean,forExpert);
-
- TH1F* hESDResidualXPerChSigma = new TH1F("hESDResidualXPerChSigma", "cluster-track residual-X per Ch: sigma;chamber ID;#sigma_{X} (cm)", nCh, 0, nCh);
- hESDResidualXPerChSigma->SetOption("P");
- hESDResidualXPerChSigma->SetMarkerStyle(kFullDotMedium);
- hESDResidualXPerChSigma->SetMarkerColor(kRed);
- Add2ESDsList(hESDResidualXPerChSigma, kESDResidualXPerChSigma,forExpert);
-
- TH1F* hESDResidualYPerChSigma = new TH1F("hESDResidualYPerChSigma", "cluster-track residual-Y per Ch: sigma;chamber ID;#sigma_{Y} (cm)", nCh, 0, nCh);
- hESDResidualYPerChSigma->SetOption("P");
- hESDResidualYPerChSigma->SetMarkerStyle(kFullDotMedium);
- hESDResidualYPerChSigma->SetMarkerColor(kRed);
- Add2ESDsList(hESDResidualYPerChSigma, kESDResidualYPerChSigma,forExpert);
-
- TH1F* hESDLocalChi2XPerChMean = new TH1F("hESDLocalChi2XPerChMean", "local chi2-X per Ch: mean;chamber ID;<local #chi^{2}_{X}>", nCh, 0, nCh);
- hESDLocalChi2XPerChMean->SetOption("P");
- hESDLocalChi2XPerChMean->SetMarkerStyle(kFullDotMedium);
- hESDLocalChi2XPerChMean->SetMarkerColor(kRed);
- Add2ESDsList(hESDLocalChi2XPerChMean, kESDLocalChi2XPerChMean,forExpert);
-
- TH1F* hESDLocalChi2YPerChMean = new TH1F("hESDLocalChi2YPerChMean", "local chi2-Y per Ch: mean;chamber ID;<local #chi^{2}_{Y}>", nCh, 0, nCh);
- hESDLocalChi2YPerChMean->SetOption("P");
- hESDLocalChi2YPerChMean->SetMarkerStyle(kFullDotMedium);
- hESDLocalChi2YPerChMean->SetMarkerColor(kRed);
- Add2ESDsList(hESDLocalChi2YPerChMean, kESDLocalChi2YPerChMean,forExpert);
-
- TH1F* hESDResidualXPerDEMean = new TH1F("hESDResidualXPerDEMean", "cluster-track residual-X per DE: mean;DetElem ID;<#Delta_{X}> (cm)", nDE, 0, nDE);
- hESDResidualXPerDEMean->SetOption("P");
- hESDResidualXPerDEMean->SetMarkerStyle(kFullDotMedium);
- hESDResidualXPerDEMean->SetMarkerColor(kRed);
- Add2ESDsList(hESDResidualXPerDEMean, kESDResidualXPerDEMean,forExpert);
-
- TH1F* hESDResidualYPerDEMean = new TH1F("hESDResidualYPerDEMean", "cluster-track residual-Y per DE: mean;DetElem ID;<#Delta_{Y}> (cm)", nDE, 0, nDE);
- hESDResidualYPerDEMean->SetOption("P");
- hESDResidualYPerDEMean->SetMarkerStyle(kFullDotMedium);
- hESDResidualYPerDEMean->SetMarkerColor(kRed);
- Add2ESDsList(hESDResidualYPerDEMean, kESDResidualYPerDEMean,forExpert);
-
- TH1F* hESDResidualXPerDESigma = new TH1F("hESDResidualXPerDESigma", "cluster-track residual-X per DE: sigma;DetElem ID;#sigma_{X} (cm)", nDE, 0, nDE);
- hESDResidualXPerDESigma->SetOption("P");
- hESDResidualXPerDESigma->SetMarkerStyle(kFullDotMedium);
- hESDResidualXPerDESigma->SetMarkerColor(kRed);
- Add2ESDsList(hESDResidualXPerDESigma, kESDResidualXPerDESigma,forExpert);
-
- TH1F* hESDResidualYPerDESigma = new TH1F("hESDResidualYPerDESigma", "cluster-track residual-Y per DE: sigma;DetElem ID;#sigma_{Y} (cm)", nDE, 0, nDE);
- hESDResidualYPerDESigma->SetOption("P");
- hESDResidualYPerDESigma->SetMarkerStyle(kFullDotMedium);
- hESDResidualYPerDESigma->SetMarkerColor(kRed);
- Add2ESDsList(hESDResidualYPerDESigma, kESDResidualYPerDESigma,forExpert);
-
- TH1F* hESDLocalChi2XPerDEMean = new TH1F("hESDLocalChi2XPerDEMean", "local chi2-X per DE: mean;DetElem ID;<local #chi^{2}_{X}>", nDE, 0, nDE);
- hESDLocalChi2XPerDEMean->SetOption("P");
- hESDLocalChi2XPerDEMean->SetMarkerStyle(kFullDotMedium);
- hESDLocalChi2XPerDEMean->SetMarkerColor(kRed);
- Add2ESDsList(hESDLocalChi2XPerDEMean, kESDLocalChi2XPerDEMean,forExpert);
-
- TH1F* hESDLocalChi2YPerDEMean = new TH1F("hESDLocalChi2YPerDEMean", "local chi2-Y per DE: mean;DetElem ID;<local #chi^{2}_{Y}>", nDE, 0, nDE);
- hESDLocalChi2YPerDEMean->SetOption("P");
- hESDLocalChi2YPerDEMean->SetMarkerStyle(kFullDotMedium);
- hESDLocalChi2YPerDEMean->SetMarkerColor(kRed);
- Add2ESDsList(hESDLocalChi2YPerDEMean, kESDLocalChi2YPerDEMean,forExpert);
-
- // temporary histograms
- fhESDnTotClustersPerCh = new TH1F("fhESDnTotClustersPerCh", "total number of associated clusters per chamber;chamber ID", nCh, 0, nCh);
- fhESDnTotClustersPerCh->SetDirectory(0);
- fhESDnTotClustersPerDE = new TH1F("fhESDnTotClustersPerDE", "total number of associated clusters per DE;DetElem ID", nDE, 0, nDE);
- fhESDnTotClustersPerDE->SetDirectory(0);
- fhESDnTotFullClustersPerDE = new TH1F("fhESDnTotFullClustersPerDE", "total number of associated clusters containing pad info per DE;DetElem ID", nDE, 0, nDE);
- fhESDnTotFullClustersPerDE->SetDirectory(0);
- fhESDSumClusterChargePerDE = new TH1F("fhESDSumClusterChargePerDE", "sum of cluster charge per DE;DetElem ID;ADC counts", nDE, 0, nDE);
- fhESDSumClusterChargePerDE->SetDirectory(0);
- fhESDSumClusterSizePerDE = new TH1F("fhESDSumClusterSizePerDE", "sum of cluster size per DE;DetElem ID;averaged number of associated pads", nDE, 0, nDE);
- fhESDSumClusterSizePerDE->SetDirectory(0);
- fhESDSumResidualXPerDE = new TH1F("fhESDSumResidualXPerDE", "sum of cluster-track residual-X per DE;DetElem ID", nDE, 0, nDE);
- fhESDSumResidualXPerDE->SetDirectory(0);
- fhESDSumResidualYPerDE = new TH1F("fhESDSumResidualYPerDE", "sum of cluster-track residual-Y per DE;DetElem ID", nDE, 0, nDE);
- fhESDSumResidualYPerDE->SetDirectory(0);
- fhESDSumResidualX2PerDE = new TH1F("fhESDSumResidualX2PerDE", "sum of cluster-track residual-X**2 per DE;DetElem ID", nDE, 0, nDE);
- fhESDSumResidualX2PerDE->SetDirectory(0);
- fhESDSumResidualY2PerDE = new TH1F("fhESDSumResidualY2PerDE", "sum of cluster-track residual-Y**2 per DE;DetElem ID", nDE, 0, nDE);
- fhESDSumResidualY2PerDE->SetDirectory(0);
- fhESDSumLocalChi2XPerDE = new TH1F("fhESDSumLocalChi2XPerDE", "sum of local chi2-X per DE;DetElem ID", nDE, 0, nDE);
- fhESDSumLocalChi2XPerDE->SetDirectory(0);
- fhESDSumLocalChi2YPerDE = new TH1F("fhESDSumLocalChi2YPerDE", "sum of local chi2-Y per DE;DetElem ID", nDE, 0, nDE);
- fhESDSumLocalChi2YPerDE->SetDirectory(0);
-
- fIsInitESDs = kTRUE;
+ if ( fTracker ) fTracker->InitESDs();
+ if ( fTrigger ) fTrigger->InitESDs();
+ //
+ ClonePerTrigClass(AliQAv1::kESDS); // this should be the last line
}
//____________________________________________________________________________
void AliMUONQADataMakerRec::MakeRaws(AliRawReader* rawReader)
{
- /// make QA for rawdata
-
- if ( ! fIsInitRaws ) {
- AliWarningStream()
- << "Skipping function due to a failure in Init" << endl;
- return;
- }
-
- if ( rawReader->GetType() == AliRawEventHeaderBase::kPhysicsEvent )
+ /// make QA for rawdata
+ /// Note that we do not call the sub-datamaker MakeRaws method
+ /// for events where the MCH or MTR is not part of the readout...
+
+ if ( !rawReader || !rawReader->GetDetectorPattern() ) return;
+
+ UInt_t clmask = rawReader->GetDetectorPattern()[0];
+
+ if ( fTracker && rawReader->GetType() == AliRawEventHeaderBase::kPhysicsEvent )
{
- rawReader->Reset();
- MakeRawsTracker(rawReader);
+ UInt_t mchMask = AliDAQ::DetectorPattern(" MUONTRK ");
+ if ( clmask & mchMask )
+ {
+ rawReader->Reset();
+ fTracker->MakeRaws(rawReader);
+ }
}
- rawReader->Reset();
- MakeRawsTrigger(rawReader);
+ if ( fTrigger && (rawReader->GetType() == AliRawEventHeaderBase::kPhysicsEvent ||
+ rawReader->GetType() == AliRawEventHeaderBase::kCalibrationEvent ) )
+ {
+ UInt_t mtrMask = AliDAQ::DetectorPattern(" MUONTRG ");
+ if ( clmask & mtrMask )
+ {
+ rawReader->Reset();
+ fTrigger->MakeRaws(rawReader);
+ }
+ }
+ //
+ IncEvCountCycleRaws();
+ IncEvCountTotalRaws();
+ //
}
-//____________________________________________________________________________
-void AliMUONQADataMakerRec::MakeRawsTracker(AliRawReader* rawReader)
+//__________________________________________________________________
+void AliMUONQADataMakerRec::MakeDigits()
{
- /// make QA for rawdata tracker
-
- ((AliMUONTrackerCalibratedDataMaker*)fTrackerDataMaker)->SetRawReader(rawReader);
-
- fTrackerDataMaker->ProcessEvent();
+ /// makes data from Digits
+
+ AliFatal("Not implemented");
}
-//____________________________________________________________________________
-void AliMUONQADataMakerRec::MakeRawsTrigger(AliRawReader* rawReader)
+//__________________________________________________________________
+void AliMUONQADataMakerRec::MakeDigits(TTree* digitsTree)
{
- /// make QA for rawdata trigger
-
- // Get trigger scalers
-
- Int_t loCircuit=0;
- AliMpCDB::LoadDDLStore();
-
- AliMUONRawStreamTrigger rawStreamTrig(rawReader);
- while (rawStreamTrig.NextDDL())
- {
- // If not a scaler event, do nothing
- Bool_t scalerEvent = rawReader->GetDataHeader()->GetL1TriggerMessage() & 0x1;
- if(!scalerEvent) break;
-
- AliMUONDDLTrigger* ddlTrigger = rawStreamTrig.GetDDLTrigger();
- AliMUONDarcHeader* darcHeader = ddlTrigger->GetDarcHeader();
-
- Int_t nReg = darcHeader->GetRegHeaderEntries();
-
- for(Int_t iReg = 0; iReg < nReg ;iReg++)
- { //reg loop
+ /// makes data from Digits
- // crate info
- AliMpTriggerCrate* crate = AliMpDDLStore::Instance()->
- GetTriggerCrate(rawStreamTrig.GetDDL(), iReg);
+ // Do nothing in case of calibration event
+ if ( GetEventSpecie() == AliRecoParam::kCalib ) return;
- AliMUONRegHeader* regHeader = darcHeader->GetRegHeaderEntry(iReg);
-
- // loop over local structures
- Int_t nLocal = regHeader->GetLocalEntries();
- for(Int_t iLocal = 0; iLocal < nLocal; iLocal++)
- {
- AliMUONLocalStruct* localStruct = regHeader->GetLocalEntry(iLocal);
-
- // if card exist
- if (!localStruct) continue;
-
- loCircuit = crate->GetLocalBoardId(localStruct->GetId());
-
- if ( !loCircuit ) continue; // empty slot
-
- AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(loCircuit, false);
-
- // skip copy cards
- if( !localBoard->IsNotified())
- continue;
-
- Int_t cathode = localStruct->GetComptXY()%2;
-
- ERaw hindex = (cathode==0) ? kTriggerScalersBP : kTriggerScalersNBP;
-
- // loop over strips
- for (Int_t ibitxy = 0; ibitxy < 16; ++ibitxy) {
- if(localStruct->GetXY1(ibitxy) > 0)
- ((TH3F*)GetRawsData(hindex))->Fill(11+0, loCircuit, ibitxy, 2*localStruct->GetXY1(ibitxy));
- if(localStruct->GetXY2(ibitxy) > 0)
- ((TH3F*)GetRawsData(hindex))->Fill(11+1, loCircuit, ibitxy, 2*localStruct->GetXY2(ibitxy));
- if(localStruct->GetXY3(ibitxy) > 0)
- ((TH3F*)GetRawsData(hindex))->Fill(11+2, loCircuit, ibitxy, 2*localStruct->GetXY3(ibitxy));
- if(localStruct->GetXY4(ibitxy) > 0)
- ((TH3F*)GetRawsData(hindex))->Fill(11+3, loCircuit, ibitxy, 2*localStruct->GetXY4(ibitxy));
- } // loop on strips
- } // iLocal
- } // iReg
- } // NextDDL
+ if ( fTracker ) fTracker->MakeDigits(digitsTree);
+ if ( fTrigger ) fTrigger->MakeDigits(digitsTree);
+ //
+ IncEvCountCycleDigits();
+ IncEvCountTotalDigits();
+ //
}
//____________________________________________________________________________
void AliMUONQADataMakerRec::MakeRecPoints(TTree* clustersTree)
{
/// Fill histograms from treeR
-
- if (fIsInitRecPointsTracker) MakeRecPointsTracker(clustersTree);
- if (fIsInitRecPointsTrigger) MakeRecPointsTrigger(clustersTree);
-}
-
-//____________________________________________________________________________
-void AliMUONQADataMakerRec::MakeRecPointsTracker(TTree* clustersTree)
-{
- /// Fill histograms related to tracker clusters
-
- // First things first : do we have clusters in the TreeR ?
- // In "normal" production mode, it should be perfectly normal
- // *not* to have them.
- // But if for some reason we de-activated the combined tracking,
- // then we have clusters in TreeR, so let's take that opportunity
- // to QA them...
-
- if (!fClusterStore)
- {
- AliCodeTimerAuto("ClusterStore creation");
- fClusterStore = AliMUONVClusterStore::Create(*clustersTree);
- if (!fClusterStore)
- {
- fIsInitRecPointsTracker = kFALSE;
- return;
- }
- }
-
- AliCodeTimerAuto("");
-
- fClusterStore->Connect(*clustersTree,kFALSE);
- clustersTree->GetEvent(0);
-
- TIter next(fClusterStore->CreateIterator());
- AliMUONVCluster* cluster;
-
- while ( ( cluster = static_cast<AliMUONVCluster*>(next()) ) )
- {
- Int_t detElemId = cluster->GetDetElemId();
- Int_t chamberId = AliMpDEManager::GetChamberId(detElemId);
-
- GetRecPointsData(kTrackerNumberOfClustersPerDE)->Fill(detElemId);
- GetRecPointsData(kTrackerClusterChargePerDE+detElemId)->Fill(cluster->GetCharge());
- GetRecPointsData(kTrackerClusterMultiplicityPerDE+detElemId)->Fill(cluster->GetNDigits());
- GetRecPointsData(kTrackerNumberOfClustersPerChamber)->Fill(chamberId);
- GetRecPointsData(kTrackerClusterChargePerChamber+chamberId)->Fill(cluster->GetCharge());
- GetRecPointsData(kTrackerClusterMultiplicityPerChamber+chamberId)->Fill(cluster->GetNDigits());
-
- }
+ // Do nothing in case of calibration event
+ if ( GetEventSpecie() == AliRecoParam::kCalib ) return;
- fClusterStore->Clear();
-}
-
-//____________________________________________________________________________
-void AliMUONQADataMakerRec::MakeRecPointsTrigger(TTree* clustersTree)
-{
- /// makes data from trigger response
-
- // Fired pads info
- fDigitStore->Clear();
-
- if (!fTriggerStore) fTriggerStore = AliMUONVTriggerStore::Create(*clustersTree);
- fTriggerStore->Clear();
- fTriggerStore->Connect(*clustersTree, false);
- clustersTree->GetEvent(0);
-
- AliMUONLocalTrigger* locTrg;
- TIter nextLoc(fTriggerStore->CreateLocalIterator());
-
- while ( ( locTrg = static_cast<AliMUONLocalTrigger*>(nextLoc()) ) )
- {
- if (locTrg->IsNull()) continue;
-
- TArrayS xyPattern[2];
- locTrg->GetXPattern(xyPattern[0]);
- locTrg->GetYPattern(xyPattern[1]);
-
- Int_t nBoard = locTrg->LoCircuit();
-
- Bool_t xTrig=locTrg->IsTrigX();
- Bool_t yTrig=locTrg->IsTrigY();
-
- if (xTrig && yTrig)
- ((TH1F*)GetRecPointsData(kTriggeredBoards))->Fill(nBoard);
-
- fDigitMaker->TriggerDigits(nBoard, xyPattern, *fDigitStore);
- }
-
- TIter nextDigit(fDigitStore->CreateIterator());
- AliMUONVDigit* mDigit;
- while ( ( mDigit = static_cast<AliMUONVDigit*>(nextDigit()) ) )
- {
- Int_t detElemId = mDigit->DetElemId();
- Int_t ch = detElemId/100;
- Int_t localBoard = mDigit->ManuId();
- Int_t channel = mDigit->ManuChannel();
- Int_t cathode = mDigit->Cathode();
- ERecPoints hindex
- = ( cathode == 0 ) ? kTriggerDigitsBendPlane : kTriggerDigitsNonBendPlane;
-
- ((TH3F*)GetRecPointsData(hindex))->Fill(ch, localBoard, channel);
- }
+ if ( fTracker ) fTracker->MakeRecPoints(clustersTree);
+ if ( fTrigger ) fTrigger->MakeRecPoints(clustersTree);
+ //
+ IncEvCountCycleRecPoints();
+ IncEvCountTotalRecPoints();
+ //
}
//____________________________________________________________________________
void AliMUONQADataMakerRec::MakeESDs(AliESDEvent* esd)
{
/// make QA data from ESDs
-
- if ( ! fIsInitESDs ) {
- AliWarningStream()
- << "Skipping function due to a failure in Init" << endl;
- return;
- }
-
- // use event specie from RecoParam
- AliRecoParam::EventSpecie_t savedEventSpecie = fEventSpecie;
- if (GetRecoParam()) SetEventSpecie(static_cast<AliRecoParam::EventSpecie_t>(GetRecoParam()->GetEventSpecie()));
-
- // load ESD event in the interface
- AliMUONESDInterface esdInterface;
- if (GetRecoParam()) AliMUONESDInterface::ResetTracker(GetRecoParam());
- else AliError("Unable to get recoParam: use default ones for residual calculation");
- esdInterface.LoadEvent(*esd);
-
- GetESDsData(kESDnTracks)->Fill(esdInterface.GetNTracks());
-
- Int_t nTrackMatchTrig = 0;
-
- // loop over tracks
- Int_t nTracks = (Int_t) esd->GetNumberOfMuonTracks();
- for (Int_t iTrack = 0; iTrack < nTracks; ++iTrack) {
-
- // get the ESD track and skip "ghosts"
- AliESDMuonTrack* esdTrack = esd->GetMuonTrack(iTrack);
- if (!esdTrack->ContainTrackerData()) continue;
-
- // get corresponding MUON track
- AliMUONTrack* track = esdInterface.FindTrack(esdTrack->GetUniqueID());
-
- if (esdTrack->ContainTriggerData()) nTrackMatchTrig++;
-
- GetESDsData(kESDMomentum)->Fill(esdTrack->P());
- GetESDsData(kESDPt)->Fill(esdTrack->Pt());
- GetESDsData(kESDRapidity)->Fill(esdTrack->Y());
- GetESDsData(kESDChi2)->Fill(track->GetNormalizedChi2());
- GetESDsData(kESDProbChi2)->Fill(TMath::Prob(track->GetGlobalChi2(),track->GetNDF()));
- GetESDsData(kESDnClustersPerTrack)->Fill(track->GetNClusters());
-
- // loop over clusters
- AliMUONTrackParam* trackParam = static_cast<AliMUONTrackParam*>(track->GetTrackParamAtCluster()->First());
- while (trackParam) {
-
- AliMUONVCluster* cluster = trackParam->GetClusterPtr();
- Int_t chId = cluster->GetChamberId();
- Int_t deID = cluster->GetDetElemId();
- Double_t residualX = cluster->GetX() - trackParam->GetNonBendingCoor();
- Double_t residualY = cluster->GetY() - trackParam->GetBendingCoor();
- Double_t localChi2X = (cluster->GetErrX2() > 0.) ? residualX*residualX/cluster->GetErrX2() : 0.;
- Double_t localChi2Y = (cluster->GetErrY2() > 0.) ? residualY*residualY/cluster->GetErrY2() : 0.;
-
- GetESDsData(kESDClusterHitMap+chId)->Fill(cluster->GetX(), cluster->GetY());
-
- fhESDnTotClustersPerCh->Fill(chId);
- fhESDnTotClustersPerDE->Fill(deID);
-
- GetESDsData(kESDClusterChargeInCh+chId)->Fill(cluster->GetCharge());
- fhESDSumClusterChargePerDE->Fill(deID, cluster->GetCharge());
-
- if (cluster->GetNDigits() > 0) { // discard clusters with pad not stored in ESD
- fhESDnTotFullClustersPerDE->Fill(deID);
- GetESDsData(kESDClusterSizeInCh+chId)->Fill(cluster->GetNDigits());
- fhESDSumClusterSizePerDE->Fill(deID, cluster->GetNDigits());
- }
-
- GetESDsData(kESDResidualXInCh+chId)->Fill(residualX);
- GetESDsData(kESDResidualYInCh+chId)->Fill(residualY);
- fhESDSumResidualXPerDE->Fill(deID, residualX);
- fhESDSumResidualYPerDE->Fill(deID, residualY);
- fhESDSumResidualX2PerDE->Fill(deID, residualX*residualX);
- fhESDSumResidualY2PerDE->Fill(deID, residualY*residualY);
-
- GetESDsData(kESDLocalChi2XInCh+chId)->Fill(localChi2X);
- GetESDsData(kESDLocalChi2YInCh+chId)->Fill(localChi2Y);
- fhESDSumLocalChi2XPerDE->Fill(deID, localChi2X);
- fhESDSumLocalChi2YPerDE->Fill(deID, localChi2Y);
-
- trackParam = static_cast<AliMUONTrackParam*>(track->GetTrackParamAtCluster()->After(trackParam));
- }
-
- }
-
- GetESDsData(kESDMatchTrig)->Fill(nTrackMatchTrig);
-
- // restore event specie
- SetEventSpecie(savedEventSpecie);
-
-}
-//____________________________________________________________________________
-void AliMUONQADataMakerRec::StartOfDetectorCycle()
-{
- /// Detector specific actions at start of cycle
+ // Do nothing in case of calibration event
+ if ( GetEventSpecie() == AliRecoParam::kCalib ) return;
-}
-
-//____________________________________________________________________________
-void AliMUONQADataMakerRec::DisplayTriggerInfo(AliQA::TASKINDEX_t task)
-{
+ if ( fTracker ) fTracker->MakeESDs(esd);
+ if ( fTrigger ) fTrigger->MakeESDs(esd);
//
- /// Display trigger information in a user-friendly way:
- /// from local board and strip numbers to their position on chambers
+ IncEvCountCycleESDs();
+ IncEvCountTotalESDs();
//
+ }
- if(task!=AliQA::kRECPOINTS && task!=AliQA::kRAWS) return;
-
- AliMUONTriggerDisplay triggerDisplay;
+//____________________________________________________________________________
+void AliMUONQADataMakerRec::ResetDetector(AliQAv1::TASKINDEX_t task)
+{
+ /// Reset internals
- TH3F* histoStrips=0x0;
- TH2F* histoDisplayStrips=0x0;
-
- for (Int_t iCath = 0; iCath < AliMpConstants::NofCathodes(); iCath++)
+ for (int spec = 0; spec < AliRecoParam::kNSpecies; ++spec)
{
- if(task==AliQA::kRECPOINTS){
- ERecPoints hindex
- = ( iCath == 0 ) ? kTriggerDigitsBendPlane : kTriggerDigitsNonBendPlane;
- histoStrips = (TH3F*)GetRecPointsData(hindex);
+ if (!AliQAv1::Instance()->IsEventSpecieSet(AliRecoParam::ConvertIndex(spec)))
+ continue;
+
+ if ( task == AliQAv1::kRAWS )
+ {
+ if (fTracker) fTracker->ResetDetectorRaws(fRawsQAList[spec]);
+ if (fTrigger) fTrigger->ResetDetectorRaws(fRawsQAList[spec]);
}
- else if(task==AliQA::kRAWS){
- ERaw hindex
- = ( iCath == 0 ) ? kTriggerScalersBP : kTriggerScalersNBP;
- histoStrips = (TH3F*)GetRawsData(hindex);
- if(histoStrips->GetEntries()==0) continue; // No scalers found
+ else if ( task == AliQAv1::kRECPOINTS )
+ {
+ if (fTracker) fTracker->ResetDetectorRecPoints(fRecPointsQAList[spec]);
+ if (fTrigger) fTrigger->ResetDetectorRecPoints(fRecPointsQAList[spec]);
}
-
- for (Int_t iChamber = 0; iChamber < AliMpConstants::NofTriggerChambers(); iChamber++)
+ else if ( task == AliQAv1::kESDS )
{
- if(task==AliQA::kRECPOINTS){
- histoDisplayStrips = (TH2F*)GetRecPointsData(kTriggerDigitsDisplay + AliMpConstants::NofTriggerChambers()*iCath + iChamber);
- }
- else if(task==AliQA::kRAWS){
- histoDisplayStrips = (TH2F*)GetRawsData(kTriggerScalersDisplay + AliMpConstants::NofTriggerChambers()*iCath + iChamber);
- }
- Int_t bin = histoStrips->GetXaxis()->FindBin(11+iChamber);
- histoStrips->GetXaxis()->SetRange(bin,bin);
- TH2F* inputHisto = (TH2F*)histoStrips->Project3D("zy");
- triggerDisplay.FillDisplayHistogram(inputHisto, histoDisplayStrips, AliMUONTriggerDisplay::kDisplayStrips, iCath, iChamber);
- } // iChamber
- } // iCath
-
- if(task==AliQA::kRAWS){
- TH2F* histoI = (TH2F*) GetRawsData(kTriggerRPCi);
- TH2F* histoHV = (TH2F*) GetRawsData(kTriggerRPChv);
- FillTriggerDCSHistos();
- for (Int_t iChamber = 0; iChamber < AliMpConstants::NofTriggerChambers(); iChamber++) {
- Int_t bin = histoI->GetXaxis()->FindBin(11+iChamber);
- TH2F* histoDisplayI = (TH2F*)GetRawsData(kTriggerIDisplay + iChamber);
- triggerDisplay.FillDisplayHistogram(histoI->ProjectionY("_px", bin, bin), histoDisplayI, AliMUONTriggerDisplay::kDisplaySlats, 0, iChamber);
- TH2F* histoDisplayHV = (TH2F*)GetRawsData(kTriggerHVDisplay + iChamber);
- bin = histoHV->GetXaxis()->FindBin(11+iChamber);
- triggerDisplay.FillDisplayHistogram(histoHV->ProjectionY("_px", bin, bin), histoDisplayHV, AliMUONTriggerDisplay::kDisplaySlats, 0, iChamber);
+ if (fTracker) fTracker->ResetDetectorESDs(fESDsQAList[spec]);
+ if (fTrigger) fTrigger->ResetDetectorESDs(fESDsQAList[spec]);
+ }
+ else if ( task == AliQAv1::kDIGITS )
+ {
+ if (fTracker) fTracker->ResetDetectorDigits(fDigitsQAList[spec]);
+ if (fTrigger) fTrigger->ResetDetectorDigits(fDigitsQAList[spec]);
+ }
+ else
+ {
+ AliFatal(Form("Not implemented for task %s",AliQAv1::GetTaskName(task).Data()));
}
}
-
- if(task==AliQA::kRECPOINTS){
- TH1F* histoBoards = (TH1F*)GetRecPointsData(kTriggeredBoards);
- TH2F* histoDisplayBoards = (TH2F*)GetRecPointsData(kTriggerBoardsDisplay);
- triggerDisplay.FillDisplayHistogram(histoBoards, histoDisplayBoards, AliMUONTriggerDisplay::kDisplayBoards, 0, 0);
- }
-}
-
-
-//_____________________________________________________________________________
-Bool_t
-AliMUONQADataMakerRec::FillTriggerDCSHistos()
-{
- /// Get HV and currents values for one trigger chamber
-
- AliCodeTimerAuto("");
-
- AliMpDEIterator deIt;
-
- deIt.First();
-
- AliMUONCalibrationData calibrationData(AliCDBManager::Instance()->GetRun());
-
- TMap* triggerDcsMap = calibrationData.TriggerDCS();
-
- AliMpDCSNamer triggerDcsNamer("TRIGGER");
-
- TH2* currHisto = 0x0;
-
- Bool_t error = kFALSE;
-
- while ( !deIt.IsDone() )
- {
- Int_t detElemId = deIt.CurrentDEId();
-
- if ( AliMpDEManager::GetStationType(detElemId) == AliMp::kStationTrigger) {
-
- Int_t iChamber = AliMpDEManager::GetChamberId(detElemId);
- Int_t slat = detElemId%100;
-
- for(Int_t iMeas=0; iMeas<AliMpDCSNamer::kNDCSMeas; iMeas++){
- TString currAlias = triggerDcsNamer.DCSChannelName(detElemId, 0, iMeas);
-
- AliDebug(2, Form("\nDetElemId %i dcsAlias %s", detElemId, currAlias.Data()));
-
- TPair* triggerDcsPair = static_cast<TPair*>(triggerDcsMap->FindObject(currAlias.Data()));
-
- if (!triggerDcsPair)
- {
- AliError(Form("Did not find expected alias (%s) for DE %d",
- currAlias.Data(),detElemId));
- error = kTRUE;
- }
- else
- {
- TObjArray* values = static_cast<TObjArray*>(triggerDcsPair->Value());
- if (!values)
- {
- AliError(Form("Could not get values for alias %s",currAlias.Data()));
- error = kTRUE;
- }
- else
- {
- TIter next(values);
- AliDCSValue* val;
-
- while ( ( val = static_cast<AliDCSValue*>(next()) ) )
- {
- Float_t hvi = val->GetFloat();
-
- AliDebug(2, Form("Value %f", hvi));
-
- switch(iMeas){
- case AliMpDCSNamer::kDCSI:
- currHisto = (TH2F*) GetRawsData(kTriggerRPCi);
- break;
- case AliMpDCSNamer::kDCSHV:
- currHisto = (TH2F*) GetRawsData(kTriggerRPChv);
- break;
- }
- Int_t binX = currHisto->GetXaxis()->FindBin(iChamber+1);
- Int_t binY = currHisto->GetYaxis()->FindBin(slat);
- currHisto->SetBinContent(binX, binY, hvi);
- } // loop on values
- } // if (!values)
- } // if (!triggerDcsPair)
- } // loop on measured types (HV and currents)
- } // if (stationType == kStationTrigger)
-
- deIt.Next();
- }
- return error;
}
//____________________________________________________________________________
-AliMUONVTrackerData* AliMUONQADataMakerRec::GetTrackerData() const
-{
-
- return fTrackerDataMaker->Data();
+void AliMUONQADataMakerRec::StartOfDetectorCycle()
+{
+ /// Detector specific actions at start of cycle
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff