#include <TROOT.h>
#include <TSystem.h>
#include <TParticle.h>
-#include <TH1D.h>
#include <TF1.h>
+#include <TGraph.h>
+#include <TCanvas.h>
+#include <TStyle.h>
+//#include <TFrame.h>
// --- Standard library ---
#include "AliPHOSBeamTestEvent.h"
#include "AliPHOSGetter.h"
#include "AliPHOSLoader.h"
+#include "AliPHOSPulseGenerator.h"
#include "AliRunLoader.h"
#include "AliStack.h"
-#include "AliPHOSRawStream.h"
+#include "AliPHOSRawDecoder.h"
#include "AliRawReaderFile.h"
+#include "AliLog.h"
+#include "AliCDBLocal.h"
+#include "AliCDBStorage.h"
+#include "AliCDBManager.h"
+#include "AliPHOSRawDigiProducer.h"
ClassImp(AliPHOSGetter)
-AliPHOSGetter * AliPHOSGetter::fgObjGetter = 0 ;
-AliPHOSLoader * AliPHOSGetter::fgPhosLoader = 0;
+AliPHOSGetter * AliPHOSGetter::fgObjGetter = 0 ;
+AliPHOSLoader * AliPHOSGetter::fgPhosLoader = 0;
+AliPHOSCalibData* AliPHOSGetter::fgCalibData = 0;
Int_t AliPHOSGetter::fgDebug = 0;
// TFile * AliPHOSGetter::fgFile = 0 ;
+AliPHOSGetter::AliPHOSGetter() :
+ fBTE(0),
+ fLoadingStatus(),
+ fNPrimaries(0),
+ fPrimaries(0),
+ fESDFile(0),
+ fESDFileName(),
+ fESD(0),
+ fESDTree(0),
+ fRawDigits(kFALSE),
+ fcdb(0)
+{
+ // ctor: this is a singleton, the ctor should never be called but cint needs it as public
+ Fatal("ctor", "AliPHOSGetter is a singleton default ctor not callable") ;
+}
+
+
//____________________________________________________________________________
-AliPHOSGetter::AliPHOSGetter(const char* headerFile, const char* version, Option_t * openingOption)
+AliPHOSGetter::AliPHOSGetter(const char* headerFile, const char* version, Option_t * openingOption) :
+ fBTE(0),
+ fLoadingStatus(),
+ fNPrimaries(0),
+ fPrimaries(0),
+ fESDFile(0),
+ fESDFileName(),
+ fESD(0),
+ fESDTree(0),
+ fRawDigits(kFALSE),
+ fcdb(0)
{
// ctor only called by Instance()
if (rl->GetAliRun() == 0x0) {
rl->LoadgAlice();
gAlice = rl->GetAliRun(); // should be removed
+ rl->LoadHeader();
}
}
fgPhosLoader = dynamic_cast<AliPHOSLoader*>(rl->GetLoader("PHOSLoader"));
else
fgPhosLoader->SetTitle(version);
-
// initialize data members
SetDebug(0) ;
fBTE = 0 ;
fESDFile = 0 ;
fESD = 0 ;
fESDTree = 0 ;
+ fRawDigits = kFALSE ;
+
}
+AliPHOSGetter::AliPHOSGetter(const AliPHOSGetter & obj) :
+ TObject(obj),
+ fBTE(0),
+ fLoadingStatus(),
+ fNPrimaries(0),
+ fPrimaries(0),
+ fESDFile(0),
+ fESDFileName(),
+ fESD(0),
+ fESDTree(0),
+ fRawDigits(kFALSE),
+ fcdb(0)
+{
+ // cpy ctor requested by Coding Convention
+ Fatal("cpy ctor", "not implemented") ;
+}
+
+//____________________________________________________________________________
+AliPHOSGetter::AliPHOSGetter(Int_t /*i*/) :
+ fBTE(0),
+ fLoadingStatus(),
+ fNPrimaries(0),
+ fPrimaries(0),
+ fESDFile(0),
+ fESDFileName(),
+ fESD(0),
+ fESDTree(0),
+ fRawDigits(kFALSE),
+ fcdb(0)
+{
+ // special constructor for onflight
+}
+
+
//____________________________________________________________________________
AliPHOSGetter::~AliPHOSGetter()
{
// dtor
- delete fgPhosLoader ;
- fgPhosLoader = 0 ;
- delete fBTE ;
- fBTE = 0 ;
- fPrimaries->Delete() ;
- delete fPrimaries ;
- fgObjGetter = 0;
+ if(fgPhosLoader){
+ delete fgPhosLoader ;
+ fgPhosLoader = 0 ;
+ }
+ if(fBTE){
+ delete fBTE ;
+ fBTE = 0 ;
+ }
+ if(fPrimaries){
+ fPrimaries->Delete() ;
+ delete fPrimaries ;
+ }
if (fESD)
delete fESD ;
if (fESDTree)
- delete fESDTree ;
+ delete fESDTree ;
+
+ fgObjGetter = 0;
}
//____________________________________________________________________________
}
//____________________________________________________________________________
-TObjArray * AliPHOSGetter::CpvRecPoints()
+TObjArray * AliPHOSGetter::CpvRecPoints() const
{
// asks the Loader to return the CPV RecPoints container
}
//____________________________________________________________________________
-TClonesArray * AliPHOSGetter::Digits()
+TClonesArray * AliPHOSGetter::Digits() const
{
// asks the Loader to return the Digits container
//____________________________________________________________________________
-TObjArray * AliPHOSGetter::EmcRecPoints()
+TObjArray * AliPHOSGetter::EmcRecPoints() const
{
// asks the Loader to return the EMC RecPoints container
}
//____________________________________________________________________________
-TClonesArray * AliPHOSGetter::TrackSegments()
+TClonesArray * AliPHOSGetter::TrackSegments() const
{
// asks the Loader to return the TrackSegments container
}
//____________________________________________________________________________
-TClonesArray * AliPHOSGetter::RecParticles()
+TClonesArray * AliPHOSGetter::RecParticles() const
{
// asks the Loader to return the TrackSegments container
{
// Reads the content of all Tree's S, D and R
- if ( event >= MaxEvent() ) {
- Error("Event", "%d not found in TreeE !", event) ;
- return ;
- }
+// if ( event >= MaxEvent() ) {
+// Error("Event", "%d not found in TreeE !", event) ;
+// return ;
+// }
AliRunLoader * rl = AliRunLoader::GetRunLoader(PhosLoader()->GetTitle());
- // checks if we are dealing with test-beam data
- TBranch * btb = rl->TreeE()->GetBranch("AliPHOSBeamTestEvent") ;
- if(btb){
- if(!fBTE)
- fBTE = new AliPHOSBeamTestEvent() ;
- btb->SetAddress(&fBTE) ;
- btb->GetEntry(event) ;
- }
- else{
- if(fBTE){
- delete fBTE ;
- fBTE = 0 ;
- }
- }
+// // checks if we are dealing with test-beam data
+// TBranch * btb = rl->TreeE()->GetBranch("AliPHOSBeamTestEvent") ;
+// if(btb){
+// if(!fBTE)
+// fBTE = new AliPHOSBeamTestEvent() ;
+// btb->SetAddress(&fBTE) ;
+// btb->GetEntry(event) ;
+// }
+// else{
+// if(fBTE){
+// delete fBTE ;
+// fBTE = 0 ;
+// }
+// }
// Loads the type of object(s) requested
rl->GetEvent(event) ;
- if( strstr(opt,"X") || (strcmp(opt,"")==0) )
+ if(strstr(opt,"X") || (strcmp(opt,"")==0)){
ReadPrimaries() ;
-
- if(strstr(opt,"H") )
+ }
+
+ if(strstr(opt,"H") ){
ReadTreeH();
-
- if(strstr(opt,"S") )
+ }
+
+ if(strstr(opt,"S") ){
ReadTreeS() ;
-
- if( strstr(opt,"D") )
+ }
+
+ if(strstr(opt,"D") ){
ReadTreeD() ;
-
- if( strstr(opt,"R") )
+ }
+
+ if(strstr(opt,"R") ){
ReadTreeR() ;
+ }
- if( strstr(opt,"T") )
+ if( strstr(opt,"T") ){
ReadTreeT() ;
+ }
- if( strstr(opt,"P") )
+ if( strstr(opt,"P") ){
ReadTreeP() ;
+ }
- if( strstr(opt,"E") )
+ if( strstr(opt,"E") ){
ReadTreeE(event) ;
+ }
- if( strstr(opt,"W") )
- ReadRaw(event) ;
-
+}
-// if( strstr(opt,"Q") )
-// ReadTreeQA() ;
+
+//____________________________________________________________________________
+void AliPHOSGetter::Event(AliRawReader *rawReader, const char* opt, Bool_t isOldRCUFormat)
+{
+ // Reads the raw event from rawReader
+ // isOldRCUFormat defines whenever to assume
+ // the old RCU format or not
+
+ if( strstr(opt,"W") ){
+ rawReader->NextEvent();
+ ReadRaw(rawReader,isOldRCUFormat) ;
+ }
}
}
//____________________________________________________________________________
- TClonesArray * AliPHOSGetter::Hits()
+ TClonesArray * AliPHOSGetter::Hits() const
{
// asks the loader to return the Hits container
// Creates and returns the pointer of the unique instance
// Must be called only when the environment has changed
- //::Info("Instance","alirunFileName=%s version=%s openingOption=%s",alirunFileName,version,openingOption);
-
if(!fgObjGetter){ // first time the getter is called
fgObjGetter = new AliPHOSGetter(alirunFileName, version, openingOption) ;
}
}
}
else {
- AliRunLoader * rl = AliRunLoader::GetRunLoader(fgPhosLoader->GetTitle()) ;
+ rl = AliRunLoader::GetRunLoader(fgPhosLoader->GetTitle()) ;
if ( strstr(version, AliConfig::GetDefaultEventFolderName()) ) // false in case of merging
delete rl ;
fgObjGetter = new AliPHOSGetter(alirunFileName, version, openingOption) ;
AliPHOS * AliPHOSGetter:: PHOS() const
{
// returns the PHOS object
- AliPHOS * phos = dynamic_cast<AliPHOS*>(PhosLoader()->GetModulesFolder()->FindObject("PHOS")) ;
+ AliPHOSLoader * loader = 0;
+ static AliPHOSLoader * oldloader = 0;
+ static AliPHOS * phos = 0;
+
+ loader = PhosLoader();
+
+ if ( loader != oldloader) {
+ phos = dynamic_cast<AliPHOS*>(loader->GetModulesFolder()->FindObject("PHOS")) ;
+ oldloader = loader;
+ }
if (!phos)
if (fgDebug)
Warning("PHOS", "PHOS module not found in module folders: %s", PhosLoader()->GetModulesFolder()->GetName() ) ;
AliPHOSGeometry * rv = 0 ;
if (PHOS() )
rv = PHOS()->GetGeometry() ;
+ else {
+ rv = AliPHOSGeometry::GetInstance();
+ if (!rv) {
+ AliError("Could not find PHOS geometry! Loading the default one !");
+ rv = AliPHOSGeometry::GetInstance("IHEP","");
+ }
+ }
return rv ;
}
}
//____________________________________________________________________________
-Int_t AliPHOSGetter::ReadRaw(Int_t event)
+void AliPHOSGetter::FitRaw(Bool_t lowGainFlag, TGraph * gLowGain, TGraph * gHighGain, TF1* signalF, Double_t & energy, Double_t & time) const
{
- // reads the raw format data, converts it into digits format and store digits in Digits()
- // container.
+ // Fits the raw signal time distribution
+
+ const Int_t kNoiseThreshold = 0 ;
+ Double_t timezero1 = 0., timezero2 = 0., timemax = 0. ;
+ Double_t signal = 0., signalmax = 0. ;
+ time = 0. ;
+ energy = 0. ;
+
+ // Create a shaper pulse object which contains all the shaper parameters
+ AliPHOSPulseGenerator pulse;
+
+ if (lowGainFlag) {
+ timezero1 = timezero2 = signalmax = timemax = 0. ;
+ signalF->FixParameter(0, pulse.GetRawFormatLowCharge()) ;
+ signalF->FixParameter(1, pulse.GetRawFormatLowGain()) ;
+ Int_t index ;
+ for (index = 0; index < pulse.GetRawFormatTimeBins(); index++) {
+ gLowGain->GetPoint(index, time, signal) ;
+ if (signal > kNoiseThreshold && timezero1 == 0.)
+ timezero1 = time ;
+ if (signal <= kNoiseThreshold && timezero1 > 0. && timezero2 == 0.)
+ timezero2 = time ;
+ if (signal > signalmax) {
+ signalmax = signal ;
+ timemax = time ;
+ }
+ }
+ signalmax /=
+ pulse.RawResponseFunctionMax(pulse.GetRawFormatLowCharge(),
+ pulse.GetRawFormatLowGain()) ;
+ if ( timezero1 + pulse.GetRawFormatTimePeak() < pulse.GetRawFormatTimeMax() * 0.4 ) { // else its noise
+ signalF->SetParameter(2, signalmax) ;
+ signalF->SetParameter(3, timezero1) ;
+ gLowGain->Fit(signalF, "QRO", "", 0., timezero2); //, "QRON") ;
+ energy = signalF->GetParameter(2) ;
+ time = signalF->GetMaximumX() - pulse.GetRawFormatTimePeak() - pulse.GetRawFormatTimeTrigger() ;
+ }
+ } else {
+ timezero1 = timezero2 = signalmax = timemax = 0. ;
+ signalF->FixParameter(0, pulse.GetRawFormatHighCharge()) ;
+ signalF->FixParameter(1, pulse.GetRawFormatHighGain()) ;
+ Int_t index ;
+ for (index = 0; index < pulse.GetRawFormatTimeBins(); index++) {
+ gHighGain->GetPoint(index, time, signal) ;
+ if (signal > kNoiseThreshold && timezero1 == 0.)
+ timezero1 = time ;
+ if (signal <= kNoiseThreshold && timezero1 > 0. && timezero2 == 0.)
+ timezero2 = time ;
+ if (signal > signalmax) {
+ signalmax = signal ;
+ timemax = time ;
+ }
+ }
+ signalmax /= pulse.RawResponseFunctionMax(pulse.GetRawFormatHighCharge(),
+ pulse.GetRawFormatHighGain()) ;;
+ if ( timezero1 + pulse.GetRawFormatTimePeak() < pulse.GetRawFormatTimeMax() * 0.4 ) { // else its noise
+ signalF->SetParameter(2, signalmax) ;
+ signalF->SetParameter(3, timezero1) ;
+ gHighGain->Fit(signalF, "QRO", "", 0., timezero2) ;
+ energy = signalF->GetParameter(2) ;
+ time = signalF->GetMaximumX() - pulse.GetRawFormatTimePeak() - pulse.GetRawFormatTimeTrigger() ;
+ }
+ }
+ if (time == 0) energy = 0 ;
+}
- AliRawReaderFile rawReader(event) ;
- AliPHOSRawStream in(&rawReader);
+//____________________________________________________________________________
+Int_t AliPHOSGetter::CalibrateRaw(Double_t energy, Int_t *relId)
+{
+ // Convert energy into digitized amplitude for a cell relId
+ // It is a user responsilibity to open CDB and set
+ // AliPHOSCalibData object by the following operators:
+ //
+ // AliCDBLocal *loc = new AliCDBLocal("deCalibDB");
+ // AliPHOSCalibData* clb = (AliPHOSCalibData*)AliCDBStorage::Instance()
+ // ->Get(path_to_calibdata,run_number);
+ // AliPHOSGetter* gime = AliPHOSGetter::Instance("galice.root");
+ // gime->SetCalibData(clb);
+
+ if (CalibData() == 0)
+ Warning("CalibrateRaw","Calibration DB is not initiated!");
+
+ Int_t module = relId[0];
+ Int_t column = relId[3];
+ Int_t row = relId[2];
+
+ Float_t gainFactor = 0.0015; // width of one Emc ADC channel in GeV
+ Float_t pedestal = 0.005; // Emc pedestals
+
+ if(CalibData()) {
+ gainFactor = CalibData()->GetADCchannelEmc (module,column,row);
+ pedestal = CalibData()->GetADCpedestalEmc(module,column,row);
+ }
- const Int_t kHighGainIdOffset = PHOSGeometry()->GetNModules()
- * PHOSGeometry()->GetNPhi()
- * PHOSGeometry()->GetNZ()
- * 2 ;
+ Int_t amp = static_cast<Int_t>( (energy - pedestal) / gainFactor + 0.5 ) ;
+ return amp;
+}
+//____________________________________________________________________________
+Int_t AliPHOSGetter::ReadRaw(AliRawReader *rawReader,Bool_t isOldRCUFormat)
+{
+ // reads the raw format data, converts it into digits format and store digits in Digits()
+ // container.
+ // isOldRCUFormat = kTRUE in case of the old RCU
+ // format used in the raw data readout.
+ // Reimplemented by Boris Polichtchouk (Jul 2006)
+ // to make it working with the Jul-Aug 2006 beam test data.
+
+ //Create raw decoder.
- Bool_t first = kTRUE ;
-
- TH1D hLowGain("hLowGain", "Low Gain", 1000, 0., PHOS()->GetRawFormatTimeMax()) ;
- TH1D hHighGain("hHighGain", "High Gain", 1000, 0., PHOS()->GetRawFormatTimeMax()) ;
-
- // fit half the gaussian decay rather than AliPHOS::RawResponseFunction because thiswould give a floating point
- // exception during error calculation ... To solve...
- TF1 * gauss = new TF1("gauss", "gaus",
- PHOS()->GetRawFormatTimePeak(),
- PHOS()->GetRawFormatTimeMax() ) ;
+ AliPHOSRawDecoder dc(rawReader);
+ dc.SetOldRCUFormat(isOldRCUFormat);
+ dc.SubtractPedestals(kTRUE);
- Int_t relId[4], id ;
- Bool_t hgflag = kFALSE ;
-
TClonesArray * digits = Digits() ;
- digits->Clear() ;
- Int_t idigit = 0 ;
-
- while ( in.Next() ) { // PHOS entries loop
-
- if ( (in.IsNewRow() || in.IsNewColumn() || in.IsNewModule()) ) {
- if (!first) {
- hLowGain.Fit(gauss, "QRON") ;
- Int_t ampL = static_cast<Int_t>(gauss->Eval(gauss->GetParameter(2)) + 0.5) ;
- Double_t timeL = PHOS()->GetRawFormatTimePeak() - gauss->GetParameter(2) ;
- if (timeL < 0 ) // happens with noise
- timeL = PHOS()->GetRawFormatTimePeak() ;
- if (hgflag) {
- hHighGain.Fit(gauss, "QRON") ;
- Int_t ampH = static_cast<Int_t>(gauss->Eval(gauss->GetParameter(2)) + 0.5) ;
- Double_t timeH = PHOS()->GetRawFormatTimePeak() - gauss->GetParameter(2) ;
- if (timeH < 0 ) // happens with noise
- timeH = PHOS()->GetRawFormatTimePeak() ;
- new((*digits)[idigit]) AliPHOSDigit( -1, id+kHighGainIdOffset, ampH, timeH) ;
- idigit++ ;
- }
- else {
- new((*digits)[idigit]) AliPHOSDigit( -1, id, ampL, timeL) ;
- idigit++ ;
- }
- }
- first = kFALSE ;
- hLowGain.Reset() ;
- hHighGain.Reset() ;
- relId[0] = in.GetModule() ;
- if ( relId[0] >= PHOS()->GetRawFormatHighGainOffset() ) {
- relId[0] -= PHOS()->GetRawFormatHighGainOffset() ;
- hgflag = kTRUE ;
- } else
- hgflag = kFALSE ;
- relId[1] = 0 ;
- relId[2] = in.GetRow() ;
- relId[3] = in.GetColumn() ;
- PHOSGeometry()->RelToAbsNumbering(relId, id) ;
+ AliPHOSRawDigiProducer pr;
+ pr.MakeDigits(digits,&dc);
+
+ //!!!!for debug!!!
+ Int_t modMax=-111;
+ Int_t colMax=-111;
+ Int_t rowMax=-111;
+ Float_t eMax=-333;
+ //!!!for debug!!!
+
+ Int_t relId[4];
+ for(Int_t iDigit=0; iDigit<digits->GetEntries(); iDigit++) {
+ AliPHOSDigit* digit = (AliPHOSDigit*)digits->At(iDigit);
+ if(digit->GetEnergy()>eMax) {
+ PHOSGeometry()->AbsToRelNumbering(digit->GetId(),relId);
+ eMax=digit->GetEnergy();
+ modMax=relId[0];
+ rowMax=relId[2];
+ colMax=relId[3];
}
- if (hgflag)
- hHighGain.Fill(
- in.GetTime() * PHOS()->GetRawFormatTimeMax() / PHOS()->GetRawFormatTimeBins(),
- in. GetSignal() * PHOS()->GetRawFormatHighGainFactor() ) ;
- else
- hLowGain.Fill(
- in.GetTime() * PHOS()->GetRawFormatTimeMax() / PHOS()->GetRawFormatTimeBins(),
- in. GetSignal() ) ;
- } // PHOS entries loop
-
- delete gauss ;
-
- return Digits()->GetEntriesFast() ;
+ }
+
+ AliDebug(1,Form("Digit with max. energy: modMax %d colMax %d rowMax %d eMax %f\n\n",
+ modMax,colMax,rowMax,eMax));
+
+ return digits->GetEntriesFast() ;
}
//____________________________________________________________________________
-Int_t AliPHOSGetter::ReadTreeD()
+Int_t AliPHOSGetter::ReadTreeD() const
{
// Read the Digits
PhosLoader()->CleanDigits() ;
- // gets TreeD from the root file (PHOS.Digits.root)
- // if ( !IsLoaded("D") ) {
- PhosLoader()->LoadDigits("UPDATE") ;
- PhosLoader()->LoadDigitizer("UPDATE") ;
- // SetLoaded("D") ;
- //}
+ PhosLoader()->LoadDigits("UPDATE") ;
+ PhosLoader()->LoadDigitizer("UPDATE") ;
return Digits()->GetEntries() ;
}
//____________________________________________________________________________
-Int_t AliPHOSGetter::ReadTreeH()
+Int_t AliPHOSGetter::ReadTreeH() const
{
// Read the Hits
PhosLoader()->CleanHits() ;
}
//____________________________________________________________________________
-Int_t AliPHOSGetter::ReadTreeR()
+Int_t AliPHOSGetter::ReadTreeR() const
{
// Read the RecPoints
}
//____________________________________________________________________________
-Int_t AliPHOSGetter::ReadTreeT()
+Int_t AliPHOSGetter::ReadTreeT() const
{
// Read the TrackSegments
return TrackSegments()->GetEntries() ;
}
//____________________________________________________________________________
-Int_t AliPHOSGetter::ReadTreeP()
+Int_t AliPHOSGetter::ReadTreeP() const
{
// Read the RecParticles
return RecParticles()->GetEntries() ;
}
//____________________________________________________________________________
-Int_t AliPHOSGetter::ReadTreeS()
+Int_t AliPHOSGetter::ReadTreeS() const
{
// Read the SDigits
}
//____________________________________________________________________________
-TClonesArray * AliPHOSGetter::SDigits()
+TClonesArray * AliPHOSGetter::SDigits() const
{
// asks the Loader to return the Digits container
else
return 0 ;
}
+//____________________________________________________________________________
+
+AliPHOSCalibData* AliPHOSGetter::CalibData()
+{
+ // Check if the instance of AliPHOSCalibData exists, and return it
+
+ return fgCalibData;
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff