MUON/AliMUONTrackerIO.cxx

   1 /**************************************************************************
   2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3 *                                                                        *
   4 * Author: The ALICE Off-line Project.                                    *
   5 * Contributors are mentioned in the code where appropriate.              *
   6 *                                                                        *
   7 * Permission to use, copy, modify and distribute this software and its   *
   8 * documentation strictly for non-commercial purposes is hereby granted   *
   9 * without fee, provided that the above copyright notice appears in all   *
  10 * copies and that both the copyright notice and this permission notice   *
  11 * appear in the supporting documentation. The authors make no claims     *
  12 * about the suitability of this software for any purpose. It is          *
  13 * provided "as is" without express or implied warranty.                  *
  14 **************************************************************************/
  15
  16 // $Id$
  17
  18 #include <cstdlib>
  19 #include "AliMUONTrackerIO.h"
  20
  21 /// \class AliMUONTrackerIO
  22 ///
  23 /// Reader class for ASCII calibration files for MUON tracker :
  24 /// converts those ASCII files into AliMUONVStore (suitable to e.g. feed
  25 /// the OCDB).
  26 ///
  27 /// \author Laurent Aphecetche, Subatech
  28
  29 /// \cond CLASSIMP
  30 ClassImp(AliMUONTrackerIO)
  31 /// \endcond
  32
  33 #include "AliDCSValue.h"
  34 #include "AliLog.h"
  35 #include "AliMUONCalibParamND.h"
  36 #include "AliMUONCalibParamNF.h"
  37 #include "AliMUONVStore.h"
  38 #include "AliMpConstants.h"
  39 #include "AliMpDDLStore.h"
  40 #include "AliMpDEManager.h"
  41 #include "AliMpDetElement.h"
  42 #include <Riostream.h>
  43 #include <TClass.h>
  44 #include <TObjString.h>
  45 #include <TSystem.h>
  46 #include <sstream>
  47
  48 //_____________________________________________________________________________
  49 AliMUONTrackerIO::AliMUONTrackerIO()
  50 {
  51   /// ctor
  52 }
  53
  54 //_____________________________________________________________________________
  55 AliMUONTrackerIO::~AliMUONTrackerIO()
  56 {
  57   /// dtor
  58 }
  59
  60 //_____________________________________________________________________________
  61 Int_t
  62 AliMUONTrackerIO::ReadOccupancy(const char* filename,AliMUONVStore& occupancyMap)
  63 {
  64   /// Read occupancy file created by online DA
  65   /// and append values to the occupancyMap store.
  66   /// Expected format of the file is :
  67   /// busPatchId manuId sumofn nevt
  68
  69   TString sFilename(gSystem->ExpandPathName(filename));
  70
  71   std::ifstream in(sFilename.Data());
  72   if (!in.good())
  73   {
  74     return kCannotOpenFile;
  75   }
  76
  77   TString datastring;
  78   ostringstream stream;
  79   char line[1024];
  80   while ( in.getline(line,1024) )
  81         stream << line << "\n";
  82
  83   in.close();
  84
  85   return DecodeOccupancy(stream.str().c_str(),occupancyMap);
  86
  87 }
  88
  89 //_____________________________________________________________________________
  90 Int_t
  91 AliMUONTrackerIO::DecodeOccupancy(const char* data, AliMUONVStore& occupancyMap)
  92 {
  93   /// Decode occupancy string created append values to the occupancyMap store.
  94   /// Expected format of the file is :
  95   /// busPatchId manuId sumofn nevt
  96
  97   if ( ! AliMpDDLStore::Instance(kFALSE) )
  98   {
  99     AliErrorClass("Mapping not loaded. Cannot work");
 100     return kNoMapping;
 101   }
 102
 103   char line[1024];
 104   istringstream in(data);
 105
 106   Int_t n(0);
 107
 108   while ( in.getline(line,1024) )
 109   {
 110     AliDebugClass(3,Form("line=%s",line));
 111     if ( line[0] == '/' && line[1] == '/' ) continue;
 112     std::istringstream sin(line);
 113
 114     Int_t busPatchId, manuId;
 115     Int_t numberOfEvents;
 116     Double_t sumn;
 117
 118     sin >> busPatchId >> manuId >> sumn >> numberOfEvents;
 119
 120     if ( busPatchId == -1 && manuId == -1 && sumn == 0 && numberOfEvents == 0 )
 121     {
 122       /// DA could not produce information (because run failed somehow).
 123       /// That's OK, but let the world know about it
 124       return kNoInfoFile;
 125     }
 126
 127     Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromBus(busPatchId);
 128
 129     AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
 130
 131     Int_t numberOfChannelsInManu = -1;
 132
 133     if (de) numberOfChannelsInManu = de->NofChannelsInManu(manuId);
 134
 135     if ( numberOfChannelsInManu <= 0 )
 136     {
 137       AliErrorClass(Form("BP %5d DE %5d MANU %5d nchannels=%d",busPatchId,detElemId,manuId,numberOfChannelsInManu));
 138       continue;
 139     }
 140
 141     AliMUONVCalibParam* occupancy =
 142     static_cast<AliMUONVCalibParam*>(occupancyMap.FindObject(detElemId,manuId));
 143     if (occupancy)
 144     {
 145       AliErrorClass(Form("DE %5d MANU %5d is already there ?!",detElemId,manuId));
 146       continue;
 147     }
 148
 149     occupancy = new AliMUONCalibParamND(5,1,detElemId,manuId,0);
 150
 151     occupancyMap.Add(occupancy);
 152
 153     occupancy->SetValueAsDouble(0,0,sumn);
 154     occupancy->SetValueAsDouble(0,1,sumn); // with only 0 and 1s, sumw = sumw2 = sumn
 155     occupancy->SetValueAsDouble(0,2,sumn);
 156     occupancy->SetValueAsInt(0,3,numberOfChannelsInManu);
 157     occupancy->SetValueAsInt(0,4,numberOfEvents);
 158     ++n;
 159   }
 160
 161   return n;
 162 }
 163
 164 //_____________________________________________________________________________
 165 Int_t
 166 AliMUONTrackerIO::ReadPedestals(const char* filename, AliMUONVStore& pedStore)
 167 {
 168   /// Read pedestal file (produced by the MUONTRKda.exe program for instance)
 169   /// and append the read values into the given VStore
 170   /// To be used when the input is a file (for instance when reading data
 171   /// from the OCDB).
 172
 173   TString sFilename(gSystem->ExpandPathName(filename));
 174
 175   std::ifstream in(sFilename.Data());
 176   if (!in.good())
 177   {
 178     return kCannotOpenFile;
 179   }
 180
 181   ostringstream stream;
 182   char line[1024];
 183   while ( in.getline(line,1024) )
 184         stream << line << "\n";
 185
 186   in.close();
 187
 188   return DecodePedestals(stream.str().c_str(),pedStore);
 189
 190 }
 191
 192 //_____________________________________________________________________________
 193 Int_t
 194 AliMUONTrackerIO::DecodePedestals(const char* data, AliMUONVStore& pedStore)
 195 {
 196   /// Read pedestal Data (produced by the MUONTRKda.exe program for instance)
 197   /// and append the read values into the given VStore
 198   /// To be used when the input is a TString (for instance when getting data
 199   /// from AMORE DB).
 200
 201   char line[1024];
 202   Int_t busPatchID, manuID, manuChannel;
 203   Float_t pedMean, pedSigma;
 204   Int_t n(0);
 205   istringstream in(data);
 206
 207   while ( in.getline(line,1024) )
 208   {
 209     AliDebugClass(3,Form("line=%s",line));
 210     if ( line[0] == '/' && line[1] == '/' ) continue;
 211     std::istringstream sin(line);
 212     sin >> busPatchID >> manuID >> manuChannel >> pedMean >> pedSigma;
 213     Int_t detElemID = AliMpDDLStore::Instance()->GetDEfromBus(busPatchID);
 214
 215     if ( !AliMpDEManager::IsValidDetElemId(detElemID) )
 216     {
 217       AliErrorClass(Form("Got an invalid DE = %d from busPatchId=%d manuId=%d",
 218                          detElemID,busPatchID,manuID));
 219       continue;
 220     }
 221
 222     AliDebugClass(3,Form("BUSPATCH %3d DETELEMID %4d MANU %3d CH %3d MEAN %7.2f SIGMA %7.2f",
 223                     busPatchID,detElemID,manuID,manuChannel,pedMean,pedSigma));
 224
 225     AliMUONVCalibParam* ped =
 226       static_cast<AliMUONVCalibParam*>(pedStore.FindObject(detElemID,manuID));
 227     if (!ped)
 228     {
 229       ped = new AliMUONCalibParamNF(2,AliMpConstants::ManuNofChannels(),
 230                                     detElemID,manuID,
 231                                     AliMUONVCalibParam::InvalidFloatValue());
 232       pedStore.Add(ped);
 233     }
 234     ped->SetValueAsFloat(manuChannel,0,pedMean);
 235     ped->SetValueAsFloat(manuChannel,1,pedSigma);
 236     ++n;
 237   }
 238
 239   return n;
 240 }
 241
 242 //_____________________________________________________________________________
 243 Int_t
 244 AliMUONTrackerIO::ReadGains(const char* filename, AliMUONVStore& gainStore,
 245                             TString& comment)
 246 {
 247   /// Read gain file (produced by the MUONTRKda.exe program for instance)
 248   /// and append the read values into the given VStore
 249   /// To be used when the input is a file (for instance when reading data
 250   /// from the OCDB).
 251
 252   comment = "";
 253
 254   TString sFilename(gSystem->ExpandPathName(filename));
 255
 256   std::ifstream in(sFilename.Data());
 257   if (!in.good())
 258   {
 259     return kCannotOpenFile;
 260   }
 261
 262   ostringstream stream;
 263   char line[1024];
 264   while ( in.getline(line,1024) )
 265         stream << line << "\n";
 266
 267   in.close();
 268
 269   return DecodeGains(stream.str().c_str(),gainStore,comment);
 270
 271 }
 272
 273 //_____________________________________________________________________________
 274 Int_t
 275 AliMUONTrackerIO::DecodeGains(const char* data, AliMUONVStore& gainStore,
 276                               TString& comment)
 277 {
 278   /// Read gain file (produced by the MUONTRKda.exe program for instance)
 279   /// and append the read values into the given VStore
 280   /// To be used when the input is a string (for instance when getting data
 281   /// from AMORE DB).
 282
 283   char line[1024];
 284   istringstream in(data);
 285   Int_t busPatchID, manuID, manuChannel;
 286   Float_t a0, a1;
 287   Int_t thres;
 288   UInt_t qual;
 289   const Int_t kSaturation(3000); // FIXME: how to get this number ?
 290   Int_t n(0);
 291   Int_t runNumber(-1);
 292   Int_t* runs(0x0);
 293   Int_t* dac(0x0);
 294   Int_t nDAC(0);
 295   Int_t nInit(0),f1nbp(0),f2nbp(0);
 296
 297   while ( in.getline(line,1024) )
 298   {
 299     if ( strlen(line) < 10 ) continue;
 300     if ( line[0] == '/' && line[1] == '/' )
 301     {
 302       TString sline(line);
 303       if ( sline.Contains("DUMMY") )
 304       {
 305         AliDebugClass(1,"Got a dummy file here");
 306         delete [] runs;
 307         delete [] dac;
 308         return kDummyFile;
 309       }
 310       if ( sline.Contains("* Run") )
 311       {
 312         TObjArray* a = sline.Tokenize(":");
 313         if ( a->GetLast() >= 1 )
 314         {
 315           TString s = static_cast<TObjString*>(a->At(1))->String();
 316           runNumber = s.Atoi();
 317           AliDebugClass(1,Form("runNumber is %d",runNumber));
 318         }
 319       }
 320       if ( sline.Contains("DAC values") )
 321       {
 322         nDAC = TString(sline(2,sline.Length()-2)).Atoi();
 323         AliDebugClass(1,Form("# of DAC values = %d",nDAC));
 324         if ( nDAC > 0 )
 325         {
 326           if ( nDAC < 100 )
 327           {
 328             delete[] runs;
 329             delete[] dac;
 330             runs = new Int_t[nDAC];
 331             dac = new Int_t[nDAC];
 332             // skip two lines
 333             in.getline(line,1024);
 334             sline = line;
 335             if (!sline.Contains("*  nInit ="))
 336             {
 337               AliErrorClass("Improper format : was expecting nInit= line...");
 338             }
 339             else
 340             {
 341               sscanf(line,"//   *  nInit = %d  *  f1nbp = %d  *  f2nbp = %d",&nInit,&f1nbp,&f2nbp);
 342               AliDebugClass(1,Form("nInit = %d",nInit));
 343               AliDebugClass(1,Form("f1nbp = %d",f1nbp));
 344               AliDebugClass(1,Form("f2nbp = %d",f2nbp));
 345             }
 346             in.getline(line,1024);
 347             in.getline(line,1024);
 348             // then get run and dac values
 349             Int_t iDAC(0);
 350             for ( Int_t i = 0; i < nDAC; ++i )
 351             {
 352               in.getline(line,1024);
 353               Int_t a,b;
 354               sscanf(line,"// %d %d",&a,&b);
 355               runs[iDAC] = a;
 356               dac[iDAC] = b;
 357               AliDebugClass(1,Form("RUN %d is DAC %d",runs[iDAC],dac[iDAC]));
 358               ++iDAC;
 359             }
 360           }
 361           else
 362           {
 363             AliErrorClass(Form("Something went wrong, as I get too big nDAC = %d",nDAC));
 364             nDAC = 0;
 365             delete [] runs;
 366             delete [] dac;
 367             return kFormatError;
 368           }
 369         }
 370       }
 371       continue;
 372     }
 373
 374     sscanf(line,"%d %d %d %f %f %d %x",&busPatchID,&manuID,&manuChannel,
 375            &a0,&a1,&thres,&qual);
 376     AliDebugClass(3,Form("line=%s",line));
 377     Int_t detElemID = AliMpDDLStore::Instance()->GetDEfromBus(busPatchID);
 378
 379     if ( !AliMpDEManager::IsValidDetElemId(detElemID) )
 380     {
 381       AliErrorClass(Form("Got an invalid DE = %d from busPatchId=%d manuId=%d",
 382                          detElemID,busPatchID,manuID));
 383       continue;
 384     }
 385
 386     AliDebugClass(3,Form("BUSPATCH %3d DETELEMID %4d MANU %3d CH %3d A0 %7.2f "
 387                     "A1 %e THRES %5d QUAL %x",
 388                     busPatchID,detElemID,manuID,manuChannel,a0,a1,thres,qual));
 389     if ( qual == 0 ) continue;
 390
 391     AliMUONVCalibParam* gain =
 392       static_cast<AliMUONVCalibParam*>(gainStore.FindObject(detElemID,manuID));
 393
 394    if (!gain)
 395     {
 396       gain = new AliMUONCalibParamNF(5,AliMpConstants::ManuNofChannels(),detElemID,manuID,0);
 397       gainStore.Add(gain);
 398     }
 399     gain->SetValueAsFloat(manuChannel,0,a0);
 400     gain->SetValueAsFloat(manuChannel,1,a1);
 401     gain->SetValueAsInt(manuChannel,2,thres);
 402     gain->SetValueAsInt(manuChannel,3,qual);
 403     gain->SetValueAsInt(manuChannel,4,kSaturation);
 404     ++n;
 405   }
 406
 407   comment = "";
 408
 409   if ( runNumber > 0 )
 410   {
 411     comment = Form("RUN %d",runNumber);
 412   }
 413
 414   for ( Int_t i = 0; i < nDAC; ++i )
 415   {
 416     comment += Form(";(RUN %d = DAC %d)",runs[i],dac[i]);
 417   }
 418   comment += Form(";(nDAC = %d)",nDAC);
 419   comment += Form(";(nInit = %d)",nInit);
 420   comment += Form(";(f1nbp = %d)",f1nbp);
 421   comment += Form(";(f2nbp = %d)",f2nbp);
 422
 423   delete[] runs;
 424   delete[] dac;
 425
 426   return n;
 427 }
 428
 429 //_____________________________________________________________________________
 430 Int_t
 431 AliMUONTrackerIO::ReadCapacitances(const char* filename, AliMUONVStore& capaStore)
 432 {
 433   /// Read capacitance file
 434   /// and append the read values into the given VStore
 435
 436   TString sFilename(gSystem->ExpandPathName(filename));
 437
 438   std::ifstream in(sFilename.Data());
 439   if (!in.good())
 440   {
 441     return kCannotOpenFile;
 442   }
 443
 444   ostringstream stream;
 445   char line[1024];
 446   while ( in.getline(line,1024) )
 447         stream << line << "\n";
 448
 449   in.close();
 450
 451   return DecodeCapacitances(stream.str().c_str(),capaStore);
 452 }
 453
 454 //_____________________________________________________________________________
 455 Int_t
 456 AliMUONTrackerIO::DecodeCapacitances(const char* data, AliMUONVStore& capaStore)
 457 {
 458   /// Read capacitances and append the read values into the given VStore
 459   /// To be used when the input is a string (for instance when getting data
 460   /// from AMORE DB).
 461
 462   Int_t ngenerated(0);
 463
 464   char line[1024];
 465   Int_t serialNumber(-1);
 466   AliMUONVCalibParam* param(0x0);
 467   istringstream in(data);
 468
 469   while ( in.getline(line,1024,'\n') )
 470   {
 471     if ( isdigit(line[0]) )
 472     {
 473       serialNumber = atoi(line);
 474       param = static_cast<AliMUONVCalibParam*>(capaStore.FindObject(serialNumber));
 475       if (param)
 476       {
 477         AliErrorClass(Form("serialNumber %d appears several times !",serialNumber));
 478         continue;
 479       }
 480       param = new AliMUONCalibParamNF(2,AliMpConstants::ManuNofChannels(),serialNumber,0,1.0);
 481       Bool_t ok = capaStore.Add(param);
 482       if (!ok)
 483       {
 484         AliErrorClass(Form("Could not set serialNumber=%d",serialNumber));
 485         continue;
 486       }
 487       continue;
 488     }
 489
 490     if (!param) continue;
 491
 492     Int_t channel;
 493     Float_t capaValue;
 494     Float_t injectionGain;
 495     sscanf(line,"%d %f %f",&channel,&capaValue,&injectionGain);
 496     AliDebugClass(1,Form("SerialNumber %10d Channel %3d Capa %f injectionGain %f",
 497                     serialNumber,channel,capaValue,injectionGain));
 498     param->SetValueAsFloat(channel,0,capaValue);
 499     param->SetValueAsFloat(channel,1,injectionGain);
 500     ++ngenerated;
 501   }
 502
 503   return ngenerated;
 504 }
 505
 506 //_____________________________________________________________________________
 507 Int_t
 508 AliMUONTrackerIO::ReadConfig(const char* filename, AliMUONVStore& confStore)
 509 {
 510   /// Read config file (produced by the MUONTRKda.exe program for instance)
 511   /// and append the read values into the given VStore
 512   /// To be used when the input is a file (for instance when reading data
 513   /// from the OCDB).
 514
 515   TString sFilename(gSystem->ExpandPathName(filename));
 516
 517   std::ifstream in(sFilename.Data());
 518   if (!in.good())
 519   {
 520     return kCannotOpenFile;
 521   }
 522
 523   ostringstream stream;
 524   char line[1024];
 525   while ( in.getline(line,1024) )
 526         stream << line << "\n";
 527
 528   in.close();
 529
 530   return DecodeConfig(stream.str().c_str(),confStore);
 531 }
 532
 533 //_____________________________________________________________________________
 534 Int_t
 535 AliMUONTrackerIO::DecodeConfig(const char* data, AliMUONVStore& confStore)
 536 {
 537   /// Read config data (produced by the MUONTRKda.exe program for instance)
 538   /// and append the read values into the given VStore
 539   /// To be used when the input is a TString (for instance when getting data
 540   /// from AMORE DB).
 541
 542   char line[1024];
 543   Int_t busPatchID, manuID;
 544   Int_t n(0);
 545   istringstream in(data);
 546
 547   while ( in.getline(line,1024) )
 548   {
 549     AliDebugClass(3,Form("line=%s",line));
 550     if ( line[0] == '#' )
 551     {
 552       continue;
 553     }
 554     std::istringstream sin(line);
 555     sin >> busPatchID >> manuID;
 556
 557     Int_t detElemID = AliMpDDLStore::Instance()->GetDEfromBus(busPatchID);
 558
 559     if ( !AliMpDEManager::IsValidDetElemId(detElemID) )
 560     {
 561       AliErrorClass(Form("Got an invalid DE = %d from busPatchId=%d manuId=%d",
 562                          detElemID,busPatchID,manuID));
 563       continue;
 564     }
 565
 566     AliMUONVCalibParam* conf =
 567     static_cast<AliMUONVCalibParam*>(confStore.FindObject(detElemID,manuID));
 568     if (!conf)
 569     {
 570       conf = new AliMUONCalibParamNF(1,1,detElemID,manuID,1);
 571       confStore.Add(conf);
 572     }
 573     ++n;
 574   }
 575
 576   return n;
 577 }
 578
 579 //_____________________________________________________________________________
 580 Int_t
 581 AliMUONTrackerIO::WriteConfig(ofstream& out, const AliMUONVStore& confStore)
 582 {
 583   /// Write the conf store as an ASCII file
 584   /// Note that we are converting (back) the detElemId into a busPatchId
 585   /// Return the number of lines written
 586
 587   if ( !AliMpDDLStore::Instance() )
 588   {
 589     cout << "ERROR: mapping not loaded. Cannot work" << endl;
 590     return 0;
 591   }
 592
 593   TIter next(confStore.CreateIterator());
 594   AliMUONVCalibParam* param;
 595   Int_t n(0);
 596
 597   while ( (param=static_cast<AliMUONVCalibParam*>(next())) )
 598   {
 599     Int_t detElemId = param->ID0();
 600     Int_t manuId = param->ID1();
 601
 602     Int_t busPatchId = AliMpDDLStore::Instance()->GetBusPatchId(detElemId,manuId);
 603     ++n;
 604
 605     out << busPatchId << " " << manuId << endl;
 606   }
 607   return n;
 608 }
 609