X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FAliPHOS.cxx;h=8c0e9e8fddc460ccdd7e0dbe4b28c7607744d97a;hb=5b494540960a8144022b4c2431cc60a842b25183;hp=ff91cd13ad42972613237c4ac325f06d8731a2af;hpb=0bc3b8eda6a9ea194e4af879b7bf9d7e2ed553a7;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/AliPHOS.cxx b/PHOS/AliPHOS.cxx index ff91cd13ad4..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,42 +66,69 @@ // --- ROOT system --- class TFile; -#include "TROOT.h" -#include "TTree.h" -#include "TFolder.h" +#include +#include +#include +#include +#include +#include +#include // --- Standard library --- // --- AliRoot header files --- +#include "AliMagF.h" #include "AliPHOS.h" +#include "AliPHOSLoader.h" #include "AliRun.h" -#include "AliMagF.h" -#include "AliPHOSGeometry.h" -#include "AliPHOSQAChecker.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; + 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; } //____________________________________________________________________________ AliPHOS::~AliPHOS() { - + if(fgCalibData) delete fgCalibData ; +} + +//____________________________________________________________________________ +AliDigitizer* AliPHOS::CreateDigitizer(AliDigitizationInput* digInput) const +{ + return new AliPHOSDigitizer(digInput); } //____________________________________________________________________________ @@ -129,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} ; @@ -138,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.} ; @@ -160,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); @@ -180,30 +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) ; @@ -228,18 +282,19 @@ void AliPHOS::CreateMaterials() AliMixture(19, "Cables $", aCA, zCA, dCA, -4, wCA) ; - - // --- Air --- - AliMaterial(99, "Air$", 14.61, 7.3, 0.001205, 30420., 67500., 0, 0) ; - + Float_t aAir[4]={12.0107,14.0067,15.9994,39.948}; + Float_t zAir[4]={6.,7.,8.,18.}; + Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827}; + Float_t dAir = 1.20479E-3; + AliMixture(99, "Air$", aAir, zAir, dAir, 4, wAir) ; + // 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, @@ -320,84 +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) ; +} + +//_____________________________________________________________________________ +void AliPHOS::Init() +{ +} + +//____________________________________________________________________________ +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 + loader->LoadDigits(); + TClonesArray* digits = loader->Digits() ; + + if (!digits) { + AliError(Form("No digits found !")); + return; + } + + // get the geometry + AliPHOSGeometry* geom = GetGeometry(); + if (!geom) { + 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 Float_t kThreshold = 1.; // skip digits below 1 ADC channel + const Int_t kAdcThreshold = 1; // Lower ADC threshold to write to raw data - // --- Set decent energy thresholds for gamma and electron tracking + Int_t prevDDL = -1; - // 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) ; + 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 = 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 = 5-relId[0]; - // --- 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.) ; + // 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: 0At(iDDL); + buffer[iDDL] = new AliAltroBuffer(fileName.Data(),mapping[iDDL]); + buffer[iDDL]->WriteDataHeader(kTRUE, kFALSE); //Dummy; + } + prevDDL = iDDL; + } + + 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 + + // calculate the time response function + } else { + 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;} + } + 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 + 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(); } //____________________________________________________________________________ -AliPHOSGeometry * AliPHOS::GetGeometry() const -{ - // gets the pointer to the AliPHOSGeometry unique instance - - return AliPHOSGeometry::GetInstance(GetTitle(),"") ; +void AliPHOS::Hits2SDigits() +{ +// create summable digits + + AliPHOSSDigitizer phosDigitizer(fLoader->GetRunLoader()->GetFileName().Data()) ; + phosDigitizer.SetEventRange(0, -1) ; // do all the events + + phosDigitizer.Digitize("all") ; +} + +//____________________________________________________________________________ +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() { - // Set tree address to hit array + // Links Hits in the Tree to Hits array TBranch *branch; char branchname[20]; - sprintf(branchname,"%s",GetName()); - + snprintf(branchname,20,"%s",GetName()); // Branch address for hit tree - TTree *treeH = gAlice->TreeH(); - if (treeH && fHits) { + TTree *treeH = fLoader->TreeH(); + if (treeH) { branch = treeH->GetBranch(branchname); - if (branch) branch->SetAddress(&fHits); + if (branch) + { + if (fHits == 0x0) fHits= new TClonesArray("AliPHOSHit",1000); + //AliInfo(Form("<%s> Setting Hits Address",GetName())); + branch->SetAddress(&fHits); + } } } -//____________________________________________________________________________ -void AliPHOS::WriteQA() -{ +//____________________________________________________________________________ +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); + } - // 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 * alarmsF = (TFolder*)gROOT->FindObjectAny("Folders/Run/Conditions/QA/PHOS") ; - 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() ; + 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; + } -