X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=TPC%2FTPCPEDESTALda.cxx;h=988a93ce4e7e688cbc47fab701e2995f8b190e43;hp=b3235cf18e072d82c44e76ba81548295a0b0288c;hb=3f9bda933a3fbd28852007a5c0b4d1df070d2c90;hpb=ac940b58176ad80ff4113b908922658388386f0d diff --git a/TPC/TPCPEDESTALda.cxx b/TPC/TPCPEDESTALda.cxx index b3235cf18e0..988a93ce4e7 100644 --- a/TPC/TPCPEDESTALda.cxx +++ b/TPC/TPCPEDESTALda.cxx @@ -26,6 +26,7 @@ TPCda_pedestal.cxx - calibration algorithm for TPC pedestal runs 18/09/2008 christian.lippmann@cern.ch : Noisy channels are output to ASCII file. Use max noise in ALTRO. 19/09/2008 J.Wiechula@gsi.de: Added export of the calibration data to the AMORE data base. Added support for configuration files. +31/01/2011 Christian.Lippmann@cern.ch : Updates for changed setup at P2 with 2 LDCs per sector contact: marian.ivanov@cern.ch @@ -38,6 +39,12 @@ and save results in a file (named from RESULT_FILE define - see below). #define FILE_ID "pedestals" #define MAPPING_FILE "tpcMapping.root" #define CONFIG_FILE "TPCPEDESTALda.conf" +#define PED_FILE "tpcPedestals.data" +#define NOISE_FILE "tpcNoise.data" +#define PEDMEM_FILE "tpcPedestalMem.data" +#define NOISY_FILE "tpcNoisyChannels.data" +#define VERYNOISY_FILE "tpcVeryNoisyChannels.data" +#define DEAD_FILE "tpcDeadChannels.data" #define AliDebugLevel() -1 extern "C" { @@ -102,6 +109,20 @@ int main(int argc, char **argv) { printf("Wrong number of arguments\n"); return -1; } + + TString daterolename(gSystem->Getenv("DATE_ROLE_NAME")); + if ( daterolename == "" ) { + printf("Error: Variable DATE_ROLE_NAME not defined! Exiting ...\n"); + return -1; + } + bool inner; + if ( daterolename.EndsWith("-0") ) inner = true; + else if ( daterolename.EndsWith("-1") ) inner = false; + else { + printf("Error: Variable DATE_ROLE_NAME neither ends with -0 nor -1 (E.g. ldc-TPC-C12-1)! Exiting ...\n"); + return -1; + } + AliLog::SetClassDebugLevel("AliTPCRawStream",-5); AliLog::SetClassDebugLevel("AliRawReaderDate",-5); AliLog::SetClassDebugLevel("AliTPCAltroMapping",-5); @@ -109,10 +130,10 @@ int main(int argc, char **argv) { /* magic line */ gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo", - "*", - "TStreamerInfo", - "RIO", - "TStreamerInfo()"); + "*", + "TStreamerInfo", + "RIO", + "TStreamerInfo()"); /* declare monitoring program */ int i, status; @@ -164,18 +185,28 @@ int main(int argc, char **argv) { } AliTPCConfigDA config(CONFIG_FILE); // check configuration + Bool_t skipAmore=kFALSE; + if ( (Int_t)config.GetValue("NoTimeAnalysis") == 1 ) { printf("WARNING: Time analysis was switched off in the configuration file!\n"); timeAnalysis=kFALSE; } + if ( (Int_t)config.GetValue("UseFastDecoder") == 1 ){ printf("Info: The fast decoder will be used for the processing.\n"); fastDecoding=kTRUE; } - + + if ( config.GetConfigurationMap()->GetValue("SkipAmore") ) { + skipAmore=((TObjString*)config.GetConfigurationMap()->GetValue("SkipAmore"))->GetString().Atoi(); + printf("TPCPEDESTALda: Skip Amore set in config\n"); + } + + // create calibration object AliTPCCalibPedestal calibPedestal(config.GetConfigurationMap()); // pedestal and noise calibration calibPedestal.SetAltroMapping(mapping->GetAltroMapping()); // Use altro mapping we got from daqDetDb + calibPedestal.SetTimeAnalysis(timeAnalysis); // pedestal(t) calibration //===========================// // loop over RAW data files // @@ -210,22 +241,19 @@ int main(int argc, char **argv) { break; } - /* skip start/end of run events */ - if ( (event->eventType != physicsEvent) && (event->eventType != calibrationEvent) ) - continue; - /* retry if got no event */ if (event==NULL) continue; + /* skip start/end of run events */ + if ( (event->eventType != physicsEvent) && (event->eventType != calibrationEvent) ) + continue; + nevents++; // get the run number runNb = event->eventRunNb; // Pedestal calibration - AliRawReader *rawReader = new AliRawReaderDate((void*)event); - if ( fastDecoding ) calibPedestal.ProcessEventFast(rawReader); - else calibPedestal.ProcessEvent(rawReader); - delete rawReader; + calibPedestal.ProcessEvent(event); /* free resources */ free(event); @@ -252,42 +280,48 @@ int main(int argc, char **argv) { // //Send objects to the AMORE DB // - printf ("AMORE part\n"); - const char *amoreDANameorig=gSystem->Getenv("AMORE_DA_NAME"); - //cheet a little -- temporary solution (hopefully) - // - //currently amoreDA uses the environment variable AMORE_DA_NAME to create the mysql - //table in which the calib objects are stored. This table is dropped each time AmoreDA - //is initialised. This of course makes a problem if we would like to store different - //calibration entries in the AMORE DB. Therefore in each DA which writes to the AMORE DB - //the AMORE_DA_NAME env variable is overwritten. - const char *role=gSystem->Getenv("DATE_ROLE_NAME"); - gSystem->Setenv("AMORE_DA_NAME",Form("%s-%s",role,FILE_ID)); - // - // end cheet - if ( role ){ - TDatime time; - TObjString info(Form("Run: %u; Date: %s",runNb,time.AsString())); - - amore::da::AmoreDA amoreDA(amore::da::AmoreDA::kSender); - Int_t status=0; - status+=amoreDA.Send("Pedestals",calibPedestal.GetCalPadPedestal()); - status+=amoreDA.Send("Noise",calibPedestal.GetCalPadRMS()); - status+=amoreDA.Send("Info",&info); - if ( status ) - printf("Waring: Failed to write one of the calib objects to the AMORE database\n"); + if (!skipAmore){ + printf ("AMORE part\n"); + const char *amoreDANameorig=gSystem->Getenv("AMORE_DA_NAME"); + //cheet a little -- temporary solution (hopefully) + // + //currently amoreDA uses the environment variable AMORE_DA_NAME to create the mysql + //table in which the calib objects are stored. This table is dropped each time AmoreDA + //is initialised. This of course makes a problem if we would like to store different + //calibration entries in the AMORE DB. Therefore in each DA which writes to the AMORE DB + //the AMORE_DA_NAME env variable is overwritten. + + //find processed sector + Char_t sideName='A'; + Int_t sector = -1; + for ( Int_t roc = 0; roc < 72; roc++ ) { + if ( !calibPedestal.GetCalRocPedestal(roc) ) continue; + if (mapping->GetSideFromRoc(roc)==1) sideName='C'; + sector = mapping->GetSectorFromRoc(roc); + } + // gSystem->Setenv("AMORE_DA_NAME",Form("TPC-%c%02d-%s",sideName,sector,FILE_ID)); + gSystem->Setenv("AMORE_DA_NAME",Form("%s-%s",gSystem->Getenv("DATE_ROLE_NAME"),FILE_ID)); + + // + // end cheet + if (sector>-1){ + TDatime time; + TObjString info(Form("Run: %u; Date: %s",runNb,time.AsSQLString())); + + amore::da::AmoreDA amoreDA(amore::da::AmoreDA::kSender); + Int_t statusDA=0; + statusDA+=amoreDA.Send("Pedestals",calibPedestal.GetCalPadPedestal()); + statusDA+=amoreDA.Send("Noise",calibPedestal.GetCalPadRMS()); + statusDA+=amoreDA.Send("Info",&info); + if ( statusDA ) + printf("Warning: Failed to write one of the calib objects to the AMORE database\n"); + } else { + printf("Warning: No data found!\n"); + } // reset env var if (amoreDANameorig) gSystem->Setenv("AMORE_DA_NAME",amoreDANameorig); - } else { - printf ("Warning: environment variable 'AMORE_DA_NAME' not set. Cannot write to the AMORE database\n"); - } - - - if ( !timeAnalysis ){ - printf("Skipping ASCII file preparation for local ALTRO configuration, as requested\n"); - return status; } - + // // Now prepare ASCII files for local ALTRO configuration through DDL. // @@ -295,16 +329,15 @@ int main(int argc, char **argv) { ofstream noisefile; ofstream pedmemfile; ofstream noisychannelfile; + ofstream verynoisychannelfile; + ofstream deadchannelfile; - char filename[255]; - sprintf(filename,"tpcPedestals.data"); - pedfile.open(filename); - sprintf(filename,"tpcNoise.data"); - noisefile.open(filename); - sprintf(filename,"tpcPedestalMem.data"); - pedmemfile.open(filename); - sprintf(filename,"tpcNoisyChannels.data"); - noisychannelfile.open(filename); + pedfile.open(PED_FILE); + noisefile.open(NOISE_FILE); + pedmemfile.open(PEDMEM_FILE); + noisychannelfile.open(NOISY_FILE); + verynoisychannelfile.open(VERYNOISY_FILE); + deadchannelfile.open(DEAD_FILE); TArrayF **timePed = calibPedestal.GetTimePedestals(); // pedestal values for each time bin @@ -312,38 +345,47 @@ int main(int argc, char **argv) { Int_t ctr_altro = 0; Int_t ctr_pattern = 0; Int_t ctr_noisy = 0; - - pedfile << 10 << std::endl; // Mark file to contain PEDESTALS per channel - noisefile << 11 << std::endl; // Mark file to contain NOISE per altro - pedmemfile << 12 << std::endl; // Mark file to contain PEDESTALs per time bin - noisychannelfile << 14 << std::endl; // Mark file to contain NOISY CHANNELS - + Int_t ctr_vnoisy = 0; + Int_t ctr_dead = 0; + + pedfile << 10 << std::endl; // Mark file to contain PEDESTALS per channel + noisefile << 11 << std::endl; // Mark file to contain NOISE per altro + pedmemfile << 12 << std::endl; // Mark file to contain PEDESTALs per time bin + noisychannelfile << 14 << std::endl; // Mark file to contain NOISY or DEAD CHANNELS + verynoisychannelfile << 14 << std::endl; // Mark file to contain NOISY or DEAD CHANNELS + deadchannelfile << 14 << std::endl; // Mark file to contain NOISY or DEAD CHANNELS + + // inner==true : calROC from ldc-0 contains: rcus 0,1 for IROC and rcu 2 for OROC + // inner==false: calROC from ldc-1 contains: nothing for IROC and rcus 3,4,5 for OROC for ( Int_t roc = 0; roc < 72; roc++ ) { if ( !calibPedestal.GetCalRocPedestal(roc) ) continue; + bool isIROC = mapping->IsIROC(roc); Int_t side = mapping->GetSideFromRoc(roc); Int_t sector = mapping->GetSectorFromRoc(roc); - //printf("Analysing ROC %d (side %d, sector %d) ...\n", roc, side, sector); - Int_t nru = mapping->IsIROC(roc) ? 2 : 4; - for ( int rcu = 0; rcu < nru; rcu++ ) { - Int_t patch = mapping->IsIROC(roc) ? rcu : rcu+2; + Int_t minrcu, maxrcu; + if ( isIROC ) { minrcu=0; maxrcu=1; } + else if ( inner ) { minrcu=2; maxrcu=2; } + else { minrcu=3; maxrcu=5; } + for ( int rcu = minrcu; rcu <= maxrcu; rcu++ ) { + //Int_t patch = mapping->IsIROC(roc) ? rcu : rcu+2; for ( int branch = 0; branch < 2; branch++ ) { - for ( int fec = 0; fec < mapping->GetNfec(patch, branch); fec++ ) { + for ( int fec = 0; fec < mapping->GetNfec(rcu, branch); fec++ ) { for ( int altro = 0; altro < 8; altro++ ) { Float_t rms = 0.; Float_t maxrms = 0.; Float_t ctr_altrochannel = 0.; for ( int channel = 0; channel < 16; channel++ ) { Int_t hwadd = mapping->CodeHWAddress(branch, fec, altro, channel); - Int_t row = mapping->GetPadRow(patch, hwadd); // row in a ROC - Int_t globalrow = mapping->GetGlobalPadRow(patch, hwadd); // row in full sector - Int_t pad = mapping->GetPad(patch, hwadd); + Int_t row = mapping->GetPadRow(rcu, hwadd); // row in a ROC + Int_t globalrow = mapping->GetGlobalPadRow(rcu, hwadd); // row in full sector + Int_t pad = mapping->GetPad(rcu, hwadd); Float_t ped = calibPedestal.GetCalRocPedestal(roc)->GetValue(row,pad); // fixed pedestal - pedfile << ctr_channel++ << "\t" << side << "\t" << sector << "\t" << patch << "\t" - << hwadd << "\t" << ped << std::endl; - // pedestal(t) + pedfile << ctr_channel++ << "\t" << side << "\t" << sector << "\t" << rcu << "\t" + << hwadd << "\t" << ped << std::endl; + // pedestal(t)=pedestal memories if ( timePed && fabs(timePed[globalrow][pad].GetSum()) > 1e-10 ) { - pedmemfile << ctr_pattern++ << "\t" << side << "\t" << sector << "\t" << patch + pedmemfile << ctr_pattern++ << "\t" << side << "\t" << sector << "\t" << rcu << "\t" << hwadd; for ( Int_t timebin = 0; timebin < 1024; timebin++ ) pedmemfile << "\t" << timePed[globalrow][pad].At(timebin); @@ -352,34 +394,38 @@ int main(int argc, char **argv) { // rms=noise Float_t rms2 = calibPedestal.GetCalRocRMS(roc)->GetValue(row,pad); if ( fabs(ped) < 1.e-10 ) { // dead channel - noisychannelfile << ctr_noisy++ << "\t" << side << "\t" << sector << "\t" - << patch << "\t" << hwadd << "\t" << rms2 << std::endl; + deadchannelfile << ctr_dead++ << "\t" << side << "\t" << sector << "\t" + << rcu << "\t" << hwadd << "\t" << rms2 << std::endl; } else if ( (ped > 1.e-10) && (rms2 > 1.e-10) ) { // not dead - // Find noisy channels - if ( ((roc<36) && (rms2 > 2.0)) || // IROC - ((roc>35) && (row<65) && (rms2 > 2.0)) || // OROC, small pads - ((roc>35) && (row>64) && (rms2 > 3.0)) ) { // OROC, large pads (50% more signal) + // Find noisy channels + if ( rms2 > 6.0 ) { // VERY noisy + verynoisychannelfile << ctr_vnoisy++ << "\t" << side << "\t" << sector << "\t" + << rcu << "\t" << hwadd << "\t" << rms2 << std::endl; + + } else if ( ((roc<36) && (rms2 > 2.0)) || // IROC + ((roc>35) && (row<65) && (rms2 > 2.0)) || // OROC, small pads + ((roc>35) && (row>64) && (rms2 > 3.0)) ) { // OROC, large pads (50% more signal) noisychannelfile << ctr_noisy++ << "\t" << side << "\t" << sector << "\t" - << patch << "\t" << hwadd << "\t" << rms2 << std::endl; - } else { - // Not noisy. Get average and maximum noise in this ALTRO - rms += rms2; - ctr_altrochannel += 1.; - if (rms2 > maxrms) maxrms = rms2; - } // end if noisy + << rcu << "\t" << hwadd << "\t" << rms2 << std::endl; + } else { + // Not noisy. Get average and maximum noise in this ALTRO + rms += rms2; + ctr_altrochannel += 1.; + if (rms2 > maxrms) maxrms = rms2; + } // end if noisy } // end if some signal } // end channel for loop Int_t hwadd = mapping->CodeHWAddress(branch, fec, altro, 0); // ALTRO address - // Noise data (rms) averaged over all channels in this ALTRO. + // Noise data (rms) averaged over all channels in this ALTRO. if ( ctr_altrochannel > 1.e-10 ) { - /* - // average noise of this ALTRO (excluding high-noise channels) - noisefile << ctr_altro << "\t" << side << "\t" << sector << "\t" << patch << "\t" + /* + // average noise of this ALTRO (excluding high-noise channels) + noisefile << ctr_altro << "\t" << side << "\t" << sector << "\t" << rcu << "\t" << hwadd << "\t" << rms/ctr_altrochannel << std::endl; - */ - // maximum noise of this ALTRO (excluding high-noise channels) - noisefile << ctr_altro << "\t" << side << "\t" << sector << "\t" << patch << "\t" - << hwadd << "\t" << maxrms << std::endl; + */ + // maximum noise of this ALTRO (excluding high-noise channels) + noisefile << ctr_altro << "\t" << side << "\t" << sector << "\t" << rcu << "\t" + << hwadd << "\t" << maxrms << std::endl; ctr_altro++; } } // end altro for loop @@ -392,7 +438,9 @@ int main(int argc, char **argv) { noisefile.close(); pedmemfile.close(); noisychannelfile.close(); - printf("Wrote ASCII files. Found %d noisy channels.\n", ctr_noisy); + verynoisychannelfile.close(); + deadchannelfile.close(); + printf("Wrote ASCII files. Found %d noisy, %d very noisy and %d dead channels.\n", ctr_noisy, ctr_vnoisy, ctr_dead); return status;