#include "AliPMDdigit.h"
#include "AliPMDClusterFinder.h"
#include "AliPMDClustering.h"
+#include "AliPMDClusteringV1.h"
#include "AliPMDcluster.h"
#include "AliPMDrecpoint1.h"
+#include "AliPMDrechit.h"
#include "AliPMDRawStream.h"
+#include "AliPMDCalibData.h"
+#include "AliPMDPedestal.h"
+#include "AliPMDddldata.h"
+
+#include "AliDAQ.h"
+#include "AliCDBManager.h"
+#include "AliCDBEntry.h"
AliPMDClusterFinder::AliPMDClusterFinder():
fRunLoader(0),
fPMDLoader(0),
+ fCalibGain(GetCalibGain()),
+ fCalibPed(GetCalibPed()),
fTreeD(0),
fTreeR(0),
fDigits(new TClonesArray("AliPMDdigit", 1000)),
fRecpoints(new TClonesArray("AliPMDrecpoint1", 1000)),
+ fRechits(new TClonesArray("AliPMDrechit", 1000)),
fNpoint(0),
+ fNhit(0),
+ fDetNo(0),
fEcut(0.)
{
//
AliPMDClusterFinder::AliPMDClusterFinder(AliRunLoader* runLoader):
fRunLoader(runLoader),
fPMDLoader(runLoader->GetLoader("PMDLoader")),
+ fCalibGain(GetCalibGain()),
+ fCalibPed(GetCalibPed()),
fTreeD(0),
fTreeR(0),
fDigits(new TClonesArray("AliPMDdigit", 1000)),
fRecpoints(new TClonesArray("AliPMDrecpoint1", 1000)),
+ fRechits(new TClonesArray("AliPMDrechit", 1000)),
fNpoint(0),
+ fNhit(0),
+ fDetNo(0),
fEcut(0.)
{
//
//
}
// ------------------------------------------------------------------------- //
+AliPMDClusterFinder::AliPMDClusterFinder(const AliPMDClusterFinder & finder):
+ TObject(finder),
+ fRunLoader(0),
+ fPMDLoader(0),
+ fCalibGain(GetCalibGain()),
+ fCalibPed(GetCalibPed()),
+ fTreeD(0),
+ fTreeR(0),
+ fDigits(NULL),
+ fRecpoints(NULL),
+ fRechits(NULL),
+ fNpoint(0),
+ fNhit(0),
+ fDetNo(0),
+ fEcut(0.)
+{
+ // copy constructor
+ AliError("Copy constructor not allowed");
+}
+// ------------------------------------------------------------------------- //
+AliPMDClusterFinder &AliPMDClusterFinder::operator=(const AliPMDClusterFinder & /*finder*/)
+{
+ // assignment op
+ AliError("Assignment Operator not allowed");
+ return *this;
+}
+// ------------------------------------------------------------------------- //
AliPMDClusterFinder::~AliPMDClusterFinder()
{
// Destructor
delete fRecpoints;
fRecpoints=0;
}
+ if (fRechits)
+ {
+ fRechits->Delete();
+ delete fRechits;
+ fRechits=0;
+ }
}
// ------------------------------------------------------------------------- //
// Converts digits to recpoints after running clustering
// algorithm on CPV plane and PREshower plane
//
+
Int_t det = 0,smn = 0;
Int_t xpos,ypos;
Float_t adc;
Int_t ismn;
Int_t idet;
- Float_t clusdata[5];
+ Float_t clusdata[6];
TObjArray *pmdcont = new TObjArray();
- AliPMDClustering *pmdclust = new AliPMDClustering();
+ AliPMDClustering *pmdclust = new AliPMDClusteringV1();
pmdclust->SetEdepCut(fEcut);
}
Int_t bufsize = 16000;
- fTreeR->Branch("PMDRecpoint", &fRecpoints, bufsize);
+ TBranch * branch1 = fTreeR->Branch("PMDRecpoint", &fRecpoints, bufsize);
+ TBranch * branch2 = fTreeR->Branch("PMDRechit", &fRechits, bufsize);
Int_t nmodules = (Int_t) fTreeD->GetEntries();
xpos = pmddigit->GetRow();
ypos = pmddigit->GetColumn();
adc = pmddigit->GetADC();
+
+ // CALIBRATION
+ Float_t gain = fCalibGain->GetGainFact(det,smn,xpos,ypos);
+ // printf("adc = %d gain = %f\n",adc,gain);
+
+ adc = adc*gain;
+
//Int_t trno = pmddigit->GetTrackNumber();
fCellADC[xpos][ypos] = (Double_t) adc;
}
clusdata[1] = pmdcl->GetClusY();
clusdata[2] = pmdcl->GetClusADC();
clusdata[3] = pmdcl->GetClusCells();
- clusdata[4] = pmdcl->GetClusRadius();
+ clusdata[4] = pmdcl->GetClusSigmaX();
+ clusdata[5] = pmdcl->GetClusSigmaY();
AddRecPoint(idet,ismn,clusdata);
+
+ Int_t ncell = (Int_t) clusdata[3];
+ for(Int_t ihit = 0; ihit < ncell; ihit++)
+ {
+ Int_t celldataX = pmdcl->GetClusCellX(ihit);
+ Int_t celldataY = pmdcl->GetClusCellY(ihit);
+ AddRecHit(celldataX, celldataY);
+ }
+ branch2->Fill();
+ ResetRechit();
}
pmdcont->Clear();
- fTreeR->Fill();
+ branch1->Fill();
ResetRecpoint();
} // modules
}
// ------------------------------------------------------------------------- //
+void AliPMDClusterFinder::Digits2RecPoints(TTree *digitsTree,
+ TTree *clustersTree)
+{
+ // Converts digits to recpoints after running clustering
+ // algorithm on CPV plane and PREshower plane
+ //
+
+ Int_t det = 0,smn = 0;
+ Int_t xpos,ypos;
+ Float_t adc;
+ Int_t ismn;
+ Int_t idet;
+ Float_t clusdata[6];
+
+ TObjArray *pmdcont = new TObjArray();
+ AliPMDClustering *pmdclust = new AliPMDClusteringV1();
+
+ pmdclust->SetEdepCut(fEcut);
+
+ AliPMDdigit *pmddigit;
+ TBranch *branch = digitsTree->GetBranch("PMDDigit");
+ branch->SetAddress(&fDigits);
+
+ ResetRecpoint();
+
+ Int_t bufsize = 16000;
+ TBranch * branch1 = clustersTree->Branch("PMDRecpoint", &fRecpoints, bufsize);
+ TBranch * branch2 = clustersTree->Branch("PMDRechit", &fRechits, bufsize);
+
+ Int_t nmodules = (Int_t) digitsTree->GetEntries();
+
+ for (Int_t imodule = 0; imodule < nmodules; imodule++)
+ {
+ ResetCellADC();
+ digitsTree->GetEntry(imodule);
+ Int_t nentries = fDigits->GetLast();
+ for (Int_t ient = 0; ient < nentries+1; ient++)
+ {
+ pmddigit = (AliPMDdigit*)fDigits->UncheckedAt(ient);
+
+ det = pmddigit->GetDetector();
+ smn = pmddigit->GetSMNumber();
+ xpos = pmddigit->GetRow();
+ ypos = pmddigit->GetColumn();
+ adc = pmddigit->GetADC();
+
+ // Pedestal Subtraction
+ Int_t pedmeanrms = fCalibPed->GetPedMeanRms(det,smn,xpos,ypos);
+ Int_t pedrms1 = (Int_t) pedmeanrms%1000;
+ Float_t pedrms = (Float_t)pedrms1/10.;
+ Float_t pedmean = (Float_t) (pedmeanrms - pedrms1)/1000.0;
+ //printf("%f %f\n",pedmean, pedrms);
+
+ Float_t adc1 = adc - (pedmean + 3.0*pedrms);
+
+ // CALIBRATION
+ Float_t gain = fCalibGain->GetGainFact(det,smn,xpos,ypos);
+ // printf("adc = %d gain = %f\n",adc,gain);
+
+ adc = adc1*gain;
+
+ //Int_t trno = pmddigit->GetTrackNumber();
+ fCellADC[xpos][ypos] = (Double_t) adc;
+ }
+
+ idet = det;
+ ismn = smn;
+ pmdclust->DoClust(idet,ismn,fCellADC,pmdcont);
+
+ Int_t nentries1 = pmdcont->GetEntries();
+
+ AliDebug(1,Form("Total number of clusters/module = %d",nentries1));
+
+ for (Int_t ient1 = 0; ient1 < nentries1; ient1++)
+ {
+ AliPMDcluster *pmdcl = (AliPMDcluster*)pmdcont->UncheckedAt(ient1);
+ idet = pmdcl->GetDetector();
+ ismn = pmdcl->GetSMN();
+ clusdata[0] = pmdcl->GetClusX();
+ clusdata[1] = pmdcl->GetClusY();
+ clusdata[2] = pmdcl->GetClusADC();
+ clusdata[3] = pmdcl->GetClusCells();
+ clusdata[4] = pmdcl->GetClusSigmaX();
+ clusdata[5] = pmdcl->GetClusSigmaY();
+
+ AddRecPoint(idet,ismn,clusdata);
+
+ Int_t ncell = (Int_t) clusdata[3];
+ for(Int_t ihit = 0; ihit < ncell; ihit++)
+ {
+ Int_t celldataX = pmdcl->GetClusCellX(ihit);
+ Int_t celldataY = pmdcl->GetClusCellY(ihit);
+ AddRecHit(celldataX, celldataY);
+ }
+ branch2->Fill();
+ ResetRechit();
+ }
+ pmdcont->Clear();
+
+ branch1->Fill();
+ ResetRecpoint();
+
+ } // modules
+
+ ResetCellADC();
+
+ // delete the pointers
+ delete pmdclust;
+ delete pmdcont;
+
+}
+// ------------------------------------------------------------------------- //
+
void AliPMDClusterFinder::Digits2RecPoints(AliRawReader *rawReader,
TTree *clustersTree)
{
// Converts RAW data to recpoints after running clustering
// algorithm on CPV and PREshower plane
//
+ // This method is called at the time of reconstruction
+
- Float_t clusdata[5];
+ Float_t clusdata[6];
+ TObjArray pmdddlcont;
TObjArray *pmdcont = new TObjArray();
- AliPMDClustering *pmdclust = new AliPMDClustering();
+ AliPMDClustering *pmdclust = new AliPMDClusteringV1();
pmdclust->SetEdepCut(fEcut);
ResetRecpoint();
Int_t bufsize = 16000;
- clustersTree->Branch("PMDRecpoint", &fRecpoints, bufsize);
+ TBranch *branch1 = clustersTree->Branch("PMDRecpoint", &fRecpoints, bufsize);
- const Int_t kDDL = 6;
+ TBranch * branch2 = clustersTree->Branch("PMDRechit", &fRechits, bufsize);
+
+ const Int_t kDDL = AliDAQ::NumberOfDdls("PMD");
const Int_t kRow = 48;
const Int_t kCol = 96;
Int_t idet = 0;
Int_t iSMN = 0;
-
+
+
for (Int_t indexDDL = 0; indexDDL < kDDL; indexDDL++)
{
if (indexDDL < 4)
ResetCellADC();
rawReader->Reset();
AliPMDRawStream pmdinput(rawReader);
- rawReader->Select(12, indexDDL, indexDDL);
- while(pmdinput.Next())
+
+ rawReader->Select("PMD", indexDDL, indexDDL);
+
+ pmdinput.DdlData(indexDDL,&pmdddlcont);
+
+ Int_t indexsmn = 0;
+ Int_t ientries = pmdddlcont.GetEntries();
+ for (Int_t ient = 0; ient < ientries; ient++)
{
- Int_t det = pmdinput.GetDetector();
- Int_t smn = pmdinput.GetSMN();
- //Int_t mcm = pmdinput.GetMCM();
- //Int_t chno = pmdinput.GetChannel();
- Int_t row = pmdinput.GetRow();
- Int_t col = pmdinput.GetColumn();
- Int_t sig = pmdinput.GetSignal();
+ AliPMDddldata *pmdddl = (AliPMDddldata*)pmdddlcont.UncheckedAt(ient);
- Int_t indexsmn = 0;
+ Int_t det = pmdddl->GetDetector();
+ Int_t smn = pmdddl->GetSMN();
+ //Int_t mcm = pmdddl->GetMCM();
+ //Int_t chno = pmdddl->GetChannel();
+ Int_t row = pmdddl->GetRow();
+ Int_t col = pmdddl->GetColumn();
+ Int_t sig = pmdddl->GetSignal();
+
+ // Pedestal Subtraction
+ Int_t pedmeanrms = fCalibPed->GetPedMeanRms(det,smn,row,col);
+ Int_t pedrms1 = (Int_t) pedmeanrms%1000;
+ Float_t pedrms = (Float_t)pedrms1/10.;
+ Float_t pedmean = (Float_t) (pedmeanrms - pedrms1)/1000.0;
+
+ //printf("%f %f\n",pedmean, pedrms);
+
+ // Float_t sig1 = (Float_t) sig;
+ Float_t sig1 = (Float_t) sig - (pedmean + 3.0*pedrms);
+
+ // CALIBRATION
+ Float_t gain = fCalibGain->GetGainFact(det,smn,row,col);
+ //printf("sig = %d gain = %f\n",sig,gain);
+ sig = (Int_t) (sig1*gain);
if (indexDDL < 4)
{
{
indexsmn = smn;
}
- else if (smn >= 12 && smn < 18)
+ else if (smn >= 18 && smn < 24)
{
- indexsmn = smn - 6;
+ indexsmn = smn - 12;
}
}
else if (indexDDL == 5)
{
indexsmn = smn - 6;
}
- else if (smn >= 18 && smn < 24)
+ else if (smn >= 12 && smn < 18)
{
- indexsmn = smn - 12;
+ indexsmn = smn - 6;
}
}
precpvADC[indexsmn][row][col] = sig;
- } // while loop
+ }
+
+ pmdddlcont.Clear();
Int_t ismn = 0;
for (Int_t indexsmn = 0; indexsmn < iSMN; indexsmn++)
(Double_t) precpvADC[indexsmn][irow][icol];
} // row
} // col
+
if (indexDDL < 4)
{
ismn = indexsmn + indexDDL * 6;
}
else if (indexsmn >= 6 && indexsmn < 12)
{
- ismn = indexsmn + 6;
+ ismn = indexsmn + 12;
}
idet = 1;
}
}
else if (indexsmn >= 6 && indexsmn < 12)
{
- ismn = indexsmn + 12;
+ ismn = indexsmn + 6;
}
idet = 1;
}
-
pmdclust->DoClust(idet,ismn,fCellADC,pmdcont);
Int_t nentries1 = pmdcont->GetEntries();
clusdata[1] = pmdcl->GetClusY();
clusdata[2] = pmdcl->GetClusADC();
clusdata[3] = pmdcl->GetClusCells();
- clusdata[4] = pmdcl->GetClusRadius();
+ clusdata[4] = pmdcl->GetClusSigmaX();
+ clusdata[5] = pmdcl->GetClusSigmaY();
AddRecPoint(idet,ismn,clusdata);
+
+ Int_t ncell = (Int_t) clusdata[3];
+ for(Int_t ihit = 0; ihit < ncell; ihit++)
+ {
+ Int_t celldataX = pmdcl->GetClusCellX(ihit);
+ Int_t celldataY = pmdcl->GetClusCellY(ihit);
+ AddRecHit(celldataX, celldataY);
+ }
+ branch2->Fill();
+ ResetRechit();
+
}
pmdcont->Clear();
- //fTreeR->Fill();
- clustersTree->Fill();
+ branch1->Fill();
ResetRecpoint();
// algorithm on CPV and PREshower plane
//
- Float_t clusdata[5];
-
+ Float_t clusdata[6];
+ TObjArray pmdddlcont;
TObjArray *pmdcont = new TObjArray();
- AliPMDClustering *pmdclust = new AliPMDClustering();
+
+ AliPMDClustering *pmdclust = new AliPMDClusteringV1();
pmdclust->SetEdepCut(fEcut);
fTreeR = fPMDLoader->TreeR();
}
Int_t bufsize = 16000;
- fTreeR->Branch("PMDRecpoint", &fRecpoints, bufsize);
+ TBranch *branch1 = fTreeR->Branch("PMDRecpoint", &fRecpoints, bufsize);
+ TBranch *branch2 = fTreeR->Branch("PMDRechit", &fRechits, bufsize);
- const Int_t kDDL = 6;
+ const Int_t kDDL = AliDAQ::NumberOfDdls("PMD");
const Int_t kRow = 48;
const Int_t kCol = 96;
}
ResetCellADC();
rawReader->Reset();
+ rawReader->Select("PMD", indexDDL, indexDDL);
+
AliPMDRawStream pmdinput(rawReader);
- rawReader->Select(12, indexDDL, indexDDL);
- while(pmdinput.Next())
+ pmdinput.DdlData(indexDDL,&pmdddlcont);
+
+ Int_t indexsmn = 0;
+ Int_t ientries = pmdddlcont.GetEntries();
+ for (Int_t ient = 0; ient < ientries; ient++)
{
- Int_t det = pmdinput.GetDetector();
- Int_t smn = pmdinput.GetSMN();
- //Int_t mcm = pmdinput.GetMCM();
- //Int_t chno = pmdinput.GetChannel();
- Int_t row = pmdinput.GetRow();
- Int_t col = pmdinput.GetColumn();
- Int_t sig = pmdinput.GetSignal();
+ AliPMDddldata *pmdddl = (AliPMDddldata*)pmdddlcont.UncheckedAt(ient);
- Int_t indexsmn = 0;
+ Int_t det = pmdddl->GetDetector();
+ Int_t smn = pmdddl->GetSMN();
+ //Int_t mcm = pmdddl->GetMCM();
+ //Int_t chno = pmdddl->GetChannel();
+ Int_t row = pmdddl->GetRow();
+ Int_t col = pmdddl->GetColumn();
+ Int_t sig = pmdddl->GetSignal();
+
+ // Pedestal Subtraction
+ Int_t pedmeanrms = fCalibPed->GetPedMeanRms(det,smn,row,col);
+ Int_t pedrms1 = (Int_t) pedmeanrms%1000;
+ Float_t pedrms = (Float_t)pedrms1/10.;
+ Float_t pedmean = (Float_t) (pedmeanrms - pedrms1)/1000.0;
+
+ //printf("%f %f\n",pedmean, pedrms);
+
+ //Float_t sig1 = (Float_t) sig;
+ Float_t sig1 = (Float_t) sig - (pedmean + 3.0*pedrms);
+ // CALIBRATION
+ Float_t gain = fCalibGain->GetGainFact(det,smn,row,col);
+
+ //printf("sig = %d gain = %f\n",sig,gain);
+ sig = (Int_t) (sig1*gain);
+
if (indexDDL < 4)
{
{
indexsmn = smn;
}
- else if (smn >= 12 && smn < 18)
+ else if (smn >= 18 && smn < 24)
{
- indexsmn = smn - 6;
+ indexsmn = smn - 12;
}
}
else if (indexDDL == 5)
{
indexsmn = smn - 6;
}
- else if (smn >= 18 && smn < 24)
+ else if (smn >= 12 && smn < 18)
{
- indexsmn = smn - 12;
+ indexsmn = smn - 6;
}
}
precpvADC[indexsmn][row][col] = sig;
- } // while loop
+
+ }
+
+ pmdddlcont.Clear();
Int_t ismn = 0;
for (Int_t indexsmn = 0; indexsmn < iSMN; indexsmn++)
(Double_t) precpvADC[indexsmn][irow][icol];
} // row
} // col
+
+
if (indexDDL < 4)
{
ismn = indexsmn + indexDDL * 6;
}
else if (indexsmn >= 6 && indexsmn < 12)
{
- ismn = indexsmn + 6;
+ ismn = indexsmn + 12;
}
idet = 1;
}
}
else if (indexsmn >= 6 && indexsmn < 12)
{
- ismn = indexsmn + 12;
+ ismn = indexsmn + 6;
}
idet = 1;
}
-
pmdclust->DoClust(idet,ismn,fCellADC,pmdcont);
Int_t nentries1 = pmdcont->GetEntries();
clusdata[1] = pmdcl->GetClusY();
clusdata[2] = pmdcl->GetClusADC();
clusdata[3] = pmdcl->GetClusCells();
- clusdata[4] = pmdcl->GetClusRadius();
+ clusdata[4] = pmdcl->GetClusSigmaX();
+ clusdata[5] = pmdcl->GetClusSigmaY();
AddRecPoint(idet,ismn,clusdata);
+
+ Int_t ncell = (Int_t) clusdata[3];
+ for(Int_t ihit = 0; ihit < ncell; ihit++)
+ {
+ Int_t celldataX = pmdcl->GetClusCellX(ihit);
+ Int_t celldataY = pmdcl->GetClusCellY(ihit);
+ AddRecHit(celldataX, celldataY);
+ }
+ branch2->Fill();
+ ResetRechit();
+
}
pmdcont->Clear();
- fTreeR->Fill();
+ branch1->Fill();
ResetRecpoint();
delete newrecpoint;
}
// ------------------------------------------------------------------------- //
+void AliPMDClusterFinder::AddRecHit(Int_t celldataX,Int_t celldataY)
+{
+ // Add associated cell hits to the Reconstructed points
+ //
+ TClonesArray &lrechits = *fRechits;
+ AliPMDrechit *newrechit;
+ newrechit = new AliPMDrechit(celldataX, celldataY);
+ new(lrechits[fNhit++]) AliPMDrechit(newrechit);
+ delete newrechit;
+}
+// ------------------------------------------------------------------------- //
void AliPMDClusterFinder::ResetCellADC()
{
// Reset the individual cell ADC value to zero
if (fRecpoints) fRecpoints->Clear();
}
// ------------------------------------------------------------------------- //
+void AliPMDClusterFinder::ResetRechit()
+{
+ // Clear the list of reconstructed points
+ fNhit = 0;
+ if (fRechits) fRechits->Clear();
+}
+// ------------------------------------------------------------------------- //
void AliPMDClusterFinder::Load()
{
// Load all the *.root files
fPMDLoader->UnloadRecPoints();
}
// ------------------------------------------------------------------------- //
+
+AliPMDCalibData* AliPMDClusterFinder::GetCalibGain() const
+{
+ // The run number will be centralized in AliCDBManager,
+ // you don't need to set it here!
+ // Added by ZA
+ AliCDBEntry *entry = AliCDBManager::Instance()->Get("PMD/Calib/Gain");
+
+ if(!entry) AliFatal("Calibration object retrieval failed! ");
+
+ AliPMDCalibData *calibdata=0;
+ if (entry) calibdata = (AliPMDCalibData*) entry->GetObject();
+
+ if (!calibdata) AliFatal("No calibration data from calibration database !");
+
+ return calibdata;
+}
+
+// ------------------------------------------------------------------------- //
+
+AliPMDPedestal* AliPMDClusterFinder::GetCalibPed() const
+{
+ // The run number will be centralized in AliCDBManager,
+ // you don't need to set it here!
+ AliCDBEntry *entry = AliCDBManager::Instance()->Get("PMD/Calib/Ped");
+
+ if(!entry) AliFatal("Pedestal object retrieval failed!");
+
+ AliPMDPedestal *pedestal = 0;
+ if (entry) pedestal = (AliPMDPedestal*) entry->GetObject();
+
+ if (!pedestal) AliFatal("No pedestal data from pedestal database !");
+
+ return pedestal;
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff