/* History of cvs commits:
*
* $Log$
+ * Revision 1.99 2006/06/28 11:36:09 cvetan
+ * New detector numbering scheme (common for DAQ/HLT/Offline). All the subdetectors shall use the AliDAQ class for the sim and rec of the raw data. The AliDAQ and raw reader classes now provide all the necessary interfaces to write and select the detector specific raw-data payload. Look into the AliDAQ.h and AliRawReader.h for more details.
+ *
+ * Revision 1.98 2006/05/11 11:30:48 cvetan
+ * Major changes in AliAltroBuffer. Now it can be used only for writing of raw data. All the corresponding read method are removed. It is based now on AliFstream in order to avoid endianess problems. The altro raw data is written always with little endian
+ *
+ * Revision 1.97 2006/04/22 10:30:17 hristov
+ * Add fEnergy to AliPHOSDigit and operate with EMC amplitude in energy units (Yu.Kharlov)
+ *
+ * Revision 1.96 2006/04/07 08:41:59 hristov
+ * Follow AliAlignObj framework and remove AliPHOSAlignData (Yu.Kharlov)
+ *
+ * Revision 1.95 2006/03/14 19:40:41 kharlov
+ * Remove De-digitizing of raw data and digitizing the raw data fit
+ *
+ * Revision 1.94 2006/03/07 18:56:25 kharlov
+ * CDB is passed via environment variable
+ *
+ * Revision 1.93 2005/11/22 08:45:11 kharlov
+ * Calibration is read from CDB if any (Boris Polichtchouk)
+ *
+ * Revision 1.92 2005/11/03 13:09:19 hristov
+ * Removing meaningless const declarations (linuxicc)
+ *
+ * Revision 1.91 2005/07/27 15:08:53 kharlov
+ * Mixture ArCO2 is corrected
+ *
* Revision 1.90 2005/06/17 07:39:07 hristov
* Removing GetDebug and SetDebug from AliRun and AliModule. Using AliLog for the messages
*
#include "AliPHOSDigit.h"
#include "AliAltroBuffer.h"
#include "AliLog.h"
+#include "AliCDBManager.h"
+#include "AliCDBEntry.h"
+#include "AliCDBStorage.h"
+#include "AliPHOSCalibData.h"
+#include "AliDAQ.h"
ClassImp(AliPHOS)
Double_t AliPHOS::fgTimePeak = 4.1E-6 ; // 4 micro seconds
Double_t AliPHOS::fgTimeTrigger = 100E-9 ; // 100ns, just for a reference
+Double_t AliPHOS::fgHighCharge = 8.2; // adjusted for a high gain range of 5.12 GeV (10 bits)
+Double_t AliPHOS::fgHighGain = 6.64;
+Double_t AliPHOS::fgHighLowGainFactor = 16.; // adjusted for a low gain range of 82 GeV (10 bits)
+
//____________________________________________________________________________
AliPHOS:: AliPHOS() : AliDetector()
{
AliPHOS::AliPHOS(const char* name, const char* title): AliDetector(name, title)
{
// ctor : title is used to identify the layout
-
- fHighCharge = 8.2 ; // adjusted for a high gain range of 5.12 GeV (10 bits)
- fHighGain = 6.64 ;
- fHighLowGainFactor = 16. ; // adjusted for a low gain range of 82 GeV (10 bits)
- fLowGainOffset = GetGeometry()->GetNModules() + 1 ;
- // offset added to the module id to distinguish high and low gain data
}
//____________________________________________________________________________
//____________________________________________________________________________
void AliPHOS::Copy(TObject &obj)const
{
- // copy method to be used byy the cpy ctor
+ // copy method to be used by the cpy ctor
TObject::Copy(obj);
-
- AliPHOS &phos = static_cast<AliPHOS &>(obj);
-
- phos.fHighCharge = fHighCharge ;
- phos.fHighGain = fHighGain ;
- phos.fHighLowGainFactor = fHighLowGainFactor ;
- phos.fLowGainOffset = fLowGainOffset;
}
//____________________________________________________________________________
return;
}
- // get the digitizer
- loader->LoadDigitizer();
- AliPHOSDigitizer * digitizer = dynamic_cast<AliPHOSDigitizer *>(loader->Digitizer()) ;
-
// get the geometry
AliPHOSGeometry* geom = GetGeometry();
if (!geom) {
}
// some digitization constants
- const Int_t kDDLOffset = 0x600; // assigned to PHOS
- const Int_t kThreshold = 1; // skip digits below this threshold
+// const Int_t kThreshold = 1; // skip digits below this threshold // YVK
+ const Float_t kThreshold = 0.001; // skip digits below 1 MeV
+ const Int_t kAdcThreshold = 1; // Lower ADC threshold to write to raw data
AliAltroBuffer* buffer = NULL;
Int_t prevDDL = -1;
// loop over digits (assume ordered digits)
for (Int_t iDigit = 0; iDigit < digits->GetEntries(); iDigit++) {
AliPHOSDigit* digit = dynamic_cast<AliPHOSDigit *>(digits->At(iDigit)) ;
- if (digit->GetAmp() < kThreshold)
+ if (digit->GetEnergy() < kThreshold)
continue;
Int_t relId[4];
geom->AbsToRelNumbering(digit->GetId(), relId);
}
// open new file and write dummy header
- TString fileName("PHOS_") ;
- fileName += (iDDL + kDDLOffset) ;
- fileName += ".ddl" ;
- buffer = new AliAltroBuffer(fileName.Data(), 1);
+ TString fileName = AliDAQ::DdlFileName("PHOS",iDDL);
+ buffer = new AliAltroBuffer(fileName.Data());
buffer->WriteDataHeader(kTRUE, kFALSE); //Dummy;
prevDDL = iDDL;
// out of time range signal (?)
if (digit->GetTimeR() > GetRawFormatTimeMax() ) {
- buffer->FillBuffer(digit->GetAmp());
+ buffer->FillBuffer((Int_t)digit->GetEnergy());
buffer->FillBuffer(GetRawFormatTimeBins() ); // time bin
buffer->FillBuffer(3); // bunch length
buffer->WriteTrailer(3, relId[3], relId[2], module); // trailer
// calculate the time response function
} else {
- Double_t energy = 0 ;
- if ( digit->GetId() <= geom->GetNModules() * geom->GetNCristalsInModule())
- energy = digit->GetAmp() * digitizer->GetEMCchannel() + digitizer->GetEMCpedestal() ;
- else
- energy = digit->GetAmp() * digitizer->GetCPVchannel() + digitizer->GetCPVpedestal() ;
-
+ Double_t energy = 0 ;
+ Int_t module = relId[0];
+ if ( digit->GetId() <= geom->GetNModules() * geom->GetNCristalsInModule()) {
+ energy=digit->GetEnergy();
+ }
+ else {
+// energy = digit->GetAmp()*digitizer->GetCPVchannel()+digitizer->GetCPVpedestal();
+ energy = 0; // CPV raw data format is now know yet
+ }
Bool_t lowgain = RawSampledResponse(digit->GetTimeR(), energy, adcValuesHigh, adcValuesLow) ;
- if (lowgain)
- buffer->WriteChannel(relId[3], relId[2], module + fLowGainOffset,
- GetRawFormatTimeBins(), adcValuesLow, kThreshold);
+ if (lowgain)
+ buffer->WriteChannel(relId[3], relId[2], module + GetGeometry()->GetNModules() + 1,
+ GetRawFormatTimeBins(), adcValuesLow , kAdcThreshold);
else
buffer->WriteChannel(relId[3], relId[2], module,
- GetRawFormatTimeBins(), adcValuesHigh, kThreshold);
+ GetRawFormatTimeBins(), adcValuesHigh, kAdcThreshold);
}
}
}
//__________________________________________________________________
-Bool_t AliPHOS::RawSampledResponse(const Double_t dtime, const Double_t damp, Int_t * adcH, Int_t * adcL) const
+Bool_t AliPHOS::RawSampledResponse(Double_t dtime, Double_t damp, Int_t * adcH, Int_t * adcL) const
{
// for a start time dtime and an amplitude damp given by digit,
// calculates the raw sampled response AliPHOS::RawResponseFunction
+ // Input: dtime - signal start time
+ // damp - signal amplitude (energy)
+ // Output: adcH - array[fkTimeBins] of 10-bit samples for high-gain channel
+ // adcL - array[fkTimeBins] of 10-bit samples for low-gain channel
const Int_t kRawSignalOverflow = 0x3FF ;
Bool_t lowGain = kFALSE ;