X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FAliPHOS.cxx;h=8c0e9e8fddc460ccdd7e0dbe4b28c7607744d97a;hb=34af53b5ec497551e2c61c5054e296aca0e44680;hp=e49219de24021d4f6f936e39939b47b9f41a43b7;hpb=65792ca0a6fa6abd7b334fa1e4afd0af7d9aa563;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/AliPHOS.cxx b/PHOS/AliPHOS.cxx index e49219de240..8c0e9e8fddc 100644 --- a/PHOS/AliPHOS.cxx +++ b/PHOS/AliPHOS.cxx @@ -13,6 +13,45 @@ * provided "as is" without express or implied warranty. * **************************************************************************/ /* $Id$ */ +/* History of cvs commits: + * + * $Log$ + * Revision 1.116 2007/10/10 09:05:10 schutz + * Changing name QualAss to QA + * + * Revision 1.115 2007/08/22 09:20:50 hristov + * Updated QA classes (Yves) + * + * Revision 1.114 2007/08/07 14:12:03 kharlov + * Quality assurance added (Yves Schutz) + * + * Revision 1.113 2007/07/18 16:29:54 policheh + * Raw Sdigits energy converted to GeV. + * + * Revision 1.112 2007/02/25 22:59:13 policheh + * Digits2Raw(): ALTRO buffer and mapping created per each DDL. + * + * Revision 1.111 2007/02/18 15:21:47 kharlov + * Corrections for memory leak in Digits2Raw due to AliAltroMapping + * + * Revision 1.110 2007/02/13 10:52:08 policheh + * Raw2SDigits() implemented + * + * Revision 1.109 2007/02/05 10:43:25 hristov + * Changes for correct initialization of Geant4 (Mihaela) + * + * Revision 1.108 2007/02/01 10:34:47 hristov + * Removing warnings on Solaris x86 + * + * Revision 1.107 2007/01/29 16:29:37 kharlov + * Digits2Raw(): special workaround for digits with time out of range + * + * Revision 1.106 2007/01/17 17:28:56 kharlov + * Extract ALTRO sample generation to a separate class AliPHOSPulseGenerator + * + * Revision 1.105 2007/01/12 21:44:29 kharlov + * Simulate and reconstruct two gains simulaneouslsy + */ //_________________________________________________________________________ // Base Class for PHOS description: @@ -27,7 +66,11 @@ // --- ROOT system --- class TFile; +#include #include +#include +#include +#include #include #include @@ -38,51 +81,54 @@ class TFile; #include "AliPHOS.h" #include "AliPHOSLoader.h" #include "AliRun.h" +#include "AliRawReader.h" #include "AliPHOSDigitizer.h" #include "AliPHOSSDigitizer.h" #include "AliPHOSDigit.h" #include "AliAltroBuffer.h" +#include "AliAltroMapping.h" +#include "AliCaloAltroMapping.h" +#include "AliLog.h" +#include "AliCDBManager.h" +#include "AliCDBEntry.h" +#include "AliCDBStorage.h" +#include "AliPHOSCalibData.h" +#include "AliPHOSPulseGenerator.h" +#include "AliDAQ.h" +#include "AliPHOSRawFitterv0.h" +#include "AliPHOSCalibData.h" +#include "AliPHOSRawDigiProducer.h" +#include "AliPHOSQAChecker.h" +#include "AliPHOSRecoParam.h" +#include "AliPHOSSimParam.h" ClassImp(AliPHOS) + //____________________________________________________________________________ -AliPHOS:: AliPHOS() : AliDetector() + AliPHOS:: AliPHOS() : AliDetector(),fgCalibData(0) { // Default ctor - fName="PHOS"; - fQATask = 0; - fTreeQA = 0; - fDebug = 0; + fName = "PHOS" ; + } //____________________________________________________________________________ -AliPHOS::AliPHOS(const char* name, const char* title): AliDetector(name, title) +AliPHOS::AliPHOS(const char* name, const char* title): AliDetector(name, title), +fgCalibData(0) { // ctor : title is used to identify the layout - - fQATask = 0; - fTreeQA = 0; - fDebug = 0; } //____________________________________________________________________________ AliPHOS::~AliPHOS() { - -} - -//____________________________________________________________________________ -void AliPHOS::Copy(AliPHOS & phos) -{ - // copy method to be used byy the cpy ctor - TObject::Copy(phos) ; - // fQATask = AliPHOSQAChecker::Copy(*(phos.fQATask)) ; - phos.fTreeQA = fTreeQA->CloneTree() ; + if(fgCalibData) delete fgCalibData ; } //____________________________________________________________________________ -AliDigitizer* AliPHOS::CreateDigitizer(AliRunDigitizer* manager) const +AliDigitizer* AliPHOS::CreateDigitizer(AliDigitizationInput* digInput) const { - return new AliPHOSDigitizer(manager); + return new AliPHOSDigitizer(digInput); } //____________________________________________________________________________ @@ -149,7 +195,7 @@ void AliPHOS::CreateMaterials() AliMixture(7, "Thermo Insul.$", aTI, zTI, dTI, -2, wTI) ; - // --- Textolitn --- + // --- Textolith --- Float_t aTX[4] = {16.0, 28.09, 12.011, 1.00794} ; Float_t zTX[4] = {8.0, 14.0, 6.0, 1.0} ; Float_t wTX[4] = {292.0, 68.0, 462.0, 736.0} ; @@ -158,12 +204,12 @@ void AliPHOS::CreateMaterials() AliMixture(8, "Textolit$", aTX, zTX, dTX, -4, wTX) ; //--- FR4 --- - Float_t aFR[3] = {28.0855, 15.9994, 17.749} ; - Float_t zFR[3] = {14., 8., 8.875} ; - Float_t wFR[3] = {.28, .32, .4} ; + Float_t aFR[4] = {16.0, 28.09, 12.011, 1.00794} ; + Float_t zFR[4] = {8.0, 14.0, 6.0, 1.0} ; + Float_t wFR[4] = {292.0, 68.0, 462.0, 736.0} ; Float_t dFR = 1.8 ; - AliMixture(9, "FR4$", aFR, zFR, dFR, -3, wFR) ; + AliMixture(9, "FR4$", aFR, zFR, dFR, -4, wFR) ; // --- The Composite Material for micromegas (so far polyetylene) --- Float_t aCM[2] = {12.01, 1.} ; @@ -180,7 +226,8 @@ void AliPHOS::CreateMaterials() Float_t aG10[4] = { 12., 1., 16., 28.} ; Float_t zG10[4] = { 6., 1., 8., 14.} ; Float_t wG10[4] = { .259, .288, .248, .205} ; - Float_t dG10 = 1.7 ; + Float_t dG10 = 1.7 ; + AliMixture(12, "G10$", aG10, zG10, dG10, -4, wG10); @@ -200,31 +247,17 @@ void AliPHOS::CreateMaterials() Float_t dAr = 0.001782 ; AliMaterial(15, "Ar$", 39.948, 18.0, dAr, 14.0, 0., 0, 0) ; - // ArCo2 - Char_t namate[21]=""; - Float_t aGM[2] ; - Float_t zGM[2] ; - Float_t wGM[2] ; - Float_t dGM ; - - Float_t absL, radL, density ; - Float_t buf[1] ; - Int_t nbuf ; - - gMC->Gfmate((*fIdmate)[15], namate, aGM[0], zGM[0], density, radL, absL, buf, nbuf) ; // Get properties of Ar - gMC->Gfmate((*fIdmate)[14], namate, aGM[1], zGM[1], density, radL, absL, buf, nbuf) ; // Get properties of CO2 - - - // Create gas mixture - - Float_t arContent = 0.80 ; // Ar-content of the Ar/CO2-mixture (80% / 20%) - - wGM[0] = arContent; - wGM[1] = 1. - arContent ; - dGM = wGM[0] * dAr + wGM[1] * dCO; + // Ar+CO2 Mixture (80% / 20%) + Float_t arContent = 0.80 ; // Ar-content of the ArCO2-mixture + Float_t aArCO[3] = {39.948, 12.0, 16.0} ; + Float_t zArCO[3] = {18.0 , 6.0, 8.0} ; + Float_t wArCO[3]; + wArCO[0] = arContent; + wArCO[1] = (1-arContent)*1; + wArCO[2] = (1-arContent)*2; + Float_t dArCO = arContent*dAr + (1-arContent)*dCO ; + AliMixture(16, "ArCO2$", aArCO, zArCO, dArCO, -3, wArCO) ; - - AliMixture(16, "ArCO2$", aGM, zGM, dGM, 2, wGM) ; // --- Stainless steel (let it be pure iron) --- AliMaterial(17, "Steel$", 55.845, 26, 7.87, 1.76, 0., 0, 0) ; @@ -260,9 +293,8 @@ void AliPHOS::CreateMaterials() // DEFINITION OF THE TRACKING MEDIA // for PHOS: idtmed[699->798] equivalent to fIdtmed[0->100] - Int_t * idtmed = fIdtmed->GetArray() - 699 ; - Int_t isxfld = gAlice->Field()->Integ() ; - Float_t sxmgmx = gAlice->Field()->Max() ; + Int_t isxfld = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ() ; + Float_t sxmgmx = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max() ; // The scintillator of the calorimeter made of PBW04 -> idtmed[699] AliMedium(0, "PHOS Xtal $", 0, 1, @@ -343,223 +375,271 @@ void AliPHOS::CreateMaterials() // Air -> idtmed[798] AliMedium(99, "Air $", 99, 0, isxfld, sxmgmx, 10.0, 1.0, 0.1, 0.1, 10.0, 0, 0) ; - - // --- Set decent energy thresholds for gamma and electron tracking - - // Tracking threshold for photons and electrons in the scintillator crystal - gMC->Gstpar(idtmed[699], "CUTGAM",0.5E-4) ; - gMC->Gstpar(idtmed[699], "CUTELE",1.0E-4) ; - - // --- Generate explicitly delta rays in the titan cover --- - gMC->Gstpar(idtmed[704], "LOSS",3.) ; - gMC->Gstpar(idtmed[704], "DRAY",1.) ; - // --- and in aluminium parts --- - gMC->Gstpar(idtmed[701], "LOSS",3.) ; - gMC->Gstpar(idtmed[701], "DRAY",1.) ; - // --- and in PIN diode - gMC->Gstpar(idtmed[705], "LOSS",3) ; - gMC->Gstpar(idtmed[705], "DRAY",1) ; - // --- and in the passive convertor - gMC->Gstpar(idtmed[712], "LOSS",3) ; - gMC->Gstpar(idtmed[712], "DRAY",1) ; - // Tracking threshold for photons and electrons in the gas ArC02 - gMC->Gstpar(idtmed[715], "CUTGAM",1.E-5) ; - gMC->Gstpar(idtmed[715], "CUTELE",1.E-5) ; - gMC->Gstpar(idtmed[715], "CUTNEU",1.E-5) ; - gMC->Gstpar(idtmed[715], "CUTHAD",1.E-5) ; - gMC->Gstpar(idtmed[715], "CUTMUO",1.E-5) ; - gMC->Gstpar(idtmed[715], "BCUTE",1.E-5) ; - gMC->Gstpar(idtmed[715], "BCUTM",1.E-5) ; - gMC->Gstpar(idtmed[715], "DCUTE",1.E-5) ; - gMC->Gstpar(idtmed[715], "DCUTM",1.E-5) ; - gMC->Gstpar(idtmed[715], "PPCUTM",1.E-5) ; - gMC->Gstpar(idtmed[715], "LOSS",2.) ; - gMC->Gstpar(idtmed[715], "DRAY",0.) ; - gMC->Gstpar(idtmed[715], "STRA",2.) ; - } -//____________________________________________________________________________ -void AliPHOS::Hits2SDigits() -{ -// create summable digits - - AliPHOSSDigitizer* phosDigitizer = - new AliPHOSSDigitizer(fLoader->GetRunLoader()->GetFileName().Data()) ; - phosDigitizer->SetEventRange(0, -1) ; // do all the events - phosDigitizer->ExecuteTask("all") ; -} - -//____________________________________________________________________________ -AliLoader* AliPHOS::MakeLoader(const char* topfoldername) +//_____________________________________________________________________________ +void AliPHOS::Init() { -//different behaviour than standard (singleton getter) -// --> to be discussed and made eventually coherent - fLoader = new AliPHOSLoader(GetName(),topfoldername); - return fLoader; -} - -//____________________________________________________________________________ -void AliPHOS::SetTreeAddress() -{ - // Links Hits in the Tree to Hits array - TBranch *branch; - char branchname[20]; - sprintf(branchname,"%s",GetName()); - - // Branch address for hit tree - TTree *treeH = TreeH(); - if (treeH) { - branch = treeH->GetBranch(branchname); - if (branch) - { - if (fHits == 0x0) fHits= new TClonesArray("AliPHOSHit",1000); - //Info("SetTreeAddress","<%s> Setting Hits Address",GetName()); - branch->SetAddress(&fHits); - } - } -} - -//____________________________________________________________________________ -void AliPHOS::WriteQA() -{ - - // Make TreeQA in the output file. - - if(fTreeQA == 0) - fTreeQA = new TTree("TreeQA", "QA Alarms") ; - // Create Alarms branches - Int_t bufferSize = 32000 ; - Int_t splitlevel = 0 ; - - TFolder* topfold = GetLoader()->GetTopFolder(); //get top aliroot folder; skowron - TString phosqafn(AliConfig::Instance()->GetQAFolderName()+"/"); //get name of QAaut folder relative to top event; skowron - phosqafn+=GetName(); //hard wired string!!! add the detector name to the pathname; skowron - TFolder * alarmsF = (TFolder*)topfold->FindObjectAny(phosqafn); //get the folder - - if (alarmsF == 0x0) - { - Error("WriteQA","Can not find folder with qa alarms"); - return; - } - TString branchName(alarmsF->GetName()); - TBranch * alarmsBranch = fTreeQA->Branch(branchName,"TFolder", &alarmsF, bufferSize, splitlevel); - TString branchTitle = branchName + " QA alarms" ; - alarmsBranch->SetTitle(branchTitle); - alarmsBranch->Fill() ; - - //fTreeQA->Fill() ; } //____________________________________________________________________________ void AliPHOS::Digits2Raw() { // convert digits of the current event to raw data + + if(AliPHOSSimParam::GetInstance()->IsEDigitizationOn()){ + AliError("Energy digitization should be OFF if use Digits2Raw") ; + } + + AliPHOSLoader * loader = static_cast(fLoader) ; // get the digits - ((AliPHOSLoader*) fLoader)->LoadDigits(); - TClonesArray* digits = ((AliPHOSLoader*) fLoader)->Digits(); + loader->LoadDigits(); + TClonesArray* digits = loader->Digits() ; + if (!digits) { - Error("Digits2Raw", "no digits"); + AliError(Form("No digits found !")); return; } // get the geometry AliPHOSGeometry* geom = GetGeometry(); if (!geom) { - Error("Digits2Raw", "no geometry"); + AliError(Form("No geometry found !")); return; } + // get mapping from OCDB + const TObjArray* maps = AliPHOSRecoParam::GetMappings(); + if(!maps) AliFatal("Cannot retrieve ALTRO mappings!!"); + // some digitization constants - const Int_t kDDLOffset = 0x600; - const Double_t kTimeMax = 1.28E-5; - const Int_t kTimeBins = 256; - const Double_t kTimePeak = 2.0E-6; - const Double_t kTimeRes = 1.5E-6; - const Int_t kThreshold = 3; - const Int_t kHighGainFactor = 40; - const Int_t kHighGainOffset = 0x200; - - AliAltroBuffer* buffer = NULL; + const Float_t kThreshold = 1.; // skip digits below 1 ADC channel + const Int_t kAdcThreshold = 1; // Lower ADC threshold to write to raw data + Int_t prevDDL = -1; - Int_t adcValuesLow[kTimeBins]; - Int_t adcValuesHigh[kTimeBins]; + + if(fgCalibData==0) + fgCalibData= new AliPHOSCalibData(-1) ; + + // Create a shaper pulse object + AliPHOSPulseGenerator *pulse = new AliPHOSPulseGenerator(); + + //Set Time step of sample + pulse->SetTimeStep(TMath::Abs(fgCalibData->GetSampleTimeStep())) ; + + Int_t *adcValuesLow = new Int_t[pulse->GetRawFormatTimeBins()]; + Int_t *adcValuesHigh= new Int_t[pulse->GetRawFormatTimeBins()]; + + const Int_t maxDDL = 20; + AliAltroBuffer *buffer[maxDDL]; + AliAltroMapping *mapping[maxDDL]; + + for(Int_t jDDL=0; jDDLGetEntries(); iDigit++) { - AliPHOSDigit* digit = (AliPHOSDigit*) digits->At(iDigit); - if (digit->GetAmp() < kThreshold) continue; + AliPHOSDigit* digit = static_cast(digits->At(iDigit)) ; + + // Skip small energy below treshold + if (digit->GetEnergy() < kThreshold) + continue; + Int_t relId[4]; geom->AbsToRelNumbering(digit->GetId(), relId); - Int_t module = relId[0]; - if (relId[1] != 0) continue; // ignore digits from CPV - Int_t iDDL = 4 * (module - 1) + (4 * (relId[2] - 1)) / geom->GetNPhi(); + Int_t module = 5-relId[0]; + + // Begin FIXME + if (relId[1] != 0) + continue; // ignore digits from CPV + // End FIXME + + Int_t row = relId[2]-1; + Int_t col = relId[3]-1; + + Int_t iRCU = -111; + + //RCU0 + if(0<=row&&row<16 && 0<=col&&col<56) iRCU=0; + + //RCU1 + if(16<=row&&row<32 && 0<=col&&col<56) iRCU=1; + + //RCU2 + if(32<=row&&row<48 && 0<=col&&col<56) iRCU=2; + + //RCU3 + if(48<=row&&row<64 && 0<=col&&col<56) iRCU=3; + + + // PHOS EMCA has 4 DDL per module. Splitting is based on the (row,column) numbers. + // here module already in PHOS online convention: 0Flush(); - buffer->WriteDataHeader(kFALSE, kFALSE); - delete buffer; - } - - // open new file and write dummy header - char fileName[256]; - sprintf(fileName, "PHOS_%d.ddl", iDDL + kDDLOffset); - buffer = new AliAltroBuffer(fileName, 1); - buffer->WriteDataHeader(kTRUE, kFALSE); //Dummy; + if (buffer[iDDL] == 0) { + // open new file and write dummy header + TString fileName = AliDAQ::DdlFileName("PHOS",iDDL); + mapping[iDDL] = (AliAltroMapping*)maps->At(iDDL); + buffer[iDDL] = new AliAltroBuffer(fileName.Data(),mapping[iDDL]); + buffer[iDDL]->WriteDataHeader(kTRUE, kFALSE); //Dummy; + } prevDDL = iDDL; } - // out of time range signal (?) - if (digit->GetTime() > kTimeMax) { - buffer->FillBuffer(digit->GetAmp()); - buffer->FillBuffer(kTimeBins); // time bin - buffer->FillBuffer(3); // bunch length - buffer->WriteTrailer(3, relId[3], relId[2], module); // trailer + AliDebug(2,Form("digit E=%.4f GeV, t=%g s, (mod,col,row)=(%d,%d,%d)\n", + digit->GetEnergy(),digit->GetTimeR(), + relId[0]-1,relId[3]-1,relId[2]-1)); + // if a signal is out of time range, write only trailer + if (digit->GetTimeR() > pulse->GetRawFormatTimeMax()*0.5 ) { + AliDebug(2,"Signal is out of time range.\n"); + buffer[iDDL]->FillBuffer(0); + buffer[iDDL]->FillBuffer(pulse->GetRawFormatTimeBins() ); // time bin + buffer[iDDL]->FillBuffer(3); // bunch length + buffer[iDDL]->WriteTrailer(3, relId[3]-1, relId[2]-1, 0); // trailer - // simulate linear rise and gaussian decay of signal + // calculate the time response function } else { - Bool_t highGain = kFALSE; - - // fill time bin values - for (Int_t iTime = 0; iTime < kTimeBins; iTime++) { - Double_t time = iTime * kTimeMax/kTimeBins; - Int_t signal = 0; - if (time < digit->GetTime() + kTimePeak) { - signal = Int_t(0.5 + digit->GetAmp() * - (time - digit->GetTime()) / kTimePeak); - } else { - signal = Int_t(0.5 + digit->GetAmp() * - TMath::Gaus(time, digit->GetTime() + kTimePeak, kTimeRes)); - } - if (signal < 0) signal = 0; - adcValuesLow[iTime] = signal; - if (signal > 0x3FF) adcValuesLow[iTime] = 0x3FF; - adcValuesHigh[iTime] = signal / kHighGainFactor; - if (adcValuesHigh[iTime] > 0) highGain = kTRUE; + Double_t energy = 0 ; + if (digit->GetId() <= geom->GetNModules() * geom->GetNCristalsInModule()) { + energy=digit->GetEnergy(); + if(energy>eMax) {eMax=energy; modMax=relId[0]; colMax=col; rowMax=row;} } - - // write low and eventually high gain channel - buffer->WriteChannel(relId[3], relId[2], module, - kTimeBins, adcValuesLow, kThreshold); - if (highGain) { - buffer->WriteChannel(relId[3], relId[2], module + kHighGainOffset, - kTimeBins, adcValuesHigh, 1); + else { + energy = 0; // CPV raw data format is now know yet } + pulse->SetAmplitude(energy); + pulse->SetTZero(digit->GetTimeR()); + Double_t r =fgCalibData->GetHighLowRatioEmc(relId[0],relId[3],relId[2]) ; + pulse->SetHG2LGRatio(r) ; + pulse->MakeSamples(); + pulse->GetSamples(adcValuesHigh, adcValuesLow) ; + + buffer[iDDL]->WriteChannel(relId[3]-1, relId[2]-1, 0, + pulse->GetRawFormatTimeBins(), adcValuesLow , kAdcThreshold); + buffer[iDDL]->WriteChannel(relId[3]-1, relId[2]-1, 1, + pulse->GetRawFormatTimeBins(), adcValuesHigh, kAdcThreshold); } } + delete [] adcValuesLow; + delete [] adcValuesHigh; // write real header and close last file - if (buffer) { - buffer->Flush(); - buffer->WriteDataHeader(kFALSE, kFALSE); - delete buffer; + for (Int_t iDDL=0; iDDLFlush(); + buffer[iDDL]->WriteDataHeader(kFALSE, kFALSE); + delete buffer[iDDL]; + //if (mapping[iDDL]) delete mapping[iDDL]; + } } + + AliDebug(1,Form("Digit with max. energy: modMax %d colMax %d rowMax %d eMax %f\n", + modMax,colMax,rowMax,eMax)); + + delete pulse; + loader->UnloadDigits(); +} + +//____________________________________________________________________________ +void AliPHOS::Hits2SDigits() +{ +// create summable digits + + AliPHOSSDigitizer phosDigitizer(fLoader->GetRunLoader()->GetFileName().Data()) ; + phosDigitizer.SetEventRange(0, -1) ; // do all the events + + phosDigitizer.Digitize("all") ; +} + - fLoader->UnloadDigits(); +//____________________________________________________________________________ +AliLoader* AliPHOS::MakeLoader(const char* topfoldername) +{ +//different behaviour than standard (singleton getter) +// --> to be discussed and made eventually coherent + fLoader = new AliPHOSLoader(GetName(),topfoldername); + return fLoader; +} + +//____________________________________________________________________________ +void AliPHOS::SetTreeAddress() +{ + // Links Hits in the Tree to Hits array + TBranch *branch; + char branchname[20]; + snprintf(branchname,20,"%s",GetName()); + // Branch address for hit tree + TTree *treeH = fLoader->TreeH(); + if (treeH) { + branch = treeH->GetBranch(branchname); + if (branch) + { + if (fHits == 0x0) fHits= new TClonesArray("AliPHOSHit",1000); + //AliInfo(Form("<%s> Setting Hits Address",GetName())); + branch->SetAddress(&fHits); + } + } } +//____________________________________________________________________________ +Bool_t AliPHOS::Raw2SDigits(AliRawReader* rawReader) +{ + + AliPHOSLoader * loader = static_cast(fLoader) ; + + TTree * tree = 0 ; + tree = loader->TreeS() ; + if ( !tree ) { + loader->MakeTree("S"); + tree = loader->TreeS() ; + } + + TClonesArray * sdigits = loader->SDigits() ; + if(!sdigits) { + loader->MakeSDigitsArray(); + sdigits = loader->SDigits(); + } + sdigits->Clear(); + + rawReader->Reset() ; + + const TObjArray* maps = AliPHOSRecoParam::GetMappings(); + if(!maps) AliFatal("Cannot retrieve ALTRO mappings!!"); + + AliAltroMapping *mapping[20]; + for(Int_t i = 0; i < 20; i++) { + mapping[i] = (AliAltroMapping*)maps->At(i); + } + + AliPHOSRawFitterv0 fitter; + + fitter.SubtractPedestals(AliPHOSSimParam::GetInstance()->EMCSubtractPedestals()); + fitter.SetAmpOffset(AliPHOSSimParam::GetInstance()->GetGlobalAltroOffset()); + fitter.SetAmpThreshold(AliPHOSSimParam::GetInstance()->GetGlobalAltroThreshold()); + + AliPHOSRawDigiProducer pr(rawReader,mapping); + pr.SetSampleQualityCut(AliPHOSSimParam::GetInstance()->GetEMCSampleQualityCut()); + pr.MakeDigits(sdigits,&fitter); + + Int_t bufferSize = 32000 ; + // TBranch * sdigitsBranch = tree->Branch("PHOS",&sdigits,bufferSize); + tree->Branch("PHOS",&sdigits,bufferSize); + tree->Fill(); + + fLoader->WriteSDigits("OVERWRITE"); + return kTRUE; + +}