MUON/MUONGenerateBusPatch.C

   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 #if !defined(__CINT__) || defined(__MAKECINT__)
  19 // ROOT includes
  20 #include <Riostream.h>
  21 #include <TROOT.h>
  22 #include <TSystem.h>
  23 #include <TString.h>
  24
  25 // MUON includes
  26 #include "AliMpBusPatch.h"
  27
  28 #endif
  29
  30 /// \ingroup macros
  31 /// \file MUONGenerateBusPatch.C
  32 /// \brief Generates buspatch id and DDL id for given detection element id
  33 ///
  34 /// - station 1 & 2 assuming 24 buspatches per quadrant
  35 /// - station345, reading from DetElemIdToSlatType.dat file and calculates
  36 ///   the number of bus patches per slat (and also number of Translator and Bridge Boards).
  37 /// Generates an output file DetElemIdToBusPatch.dat.out, preserve from overwriting
  38 ///
  39 /// List of changes:
  40 /// - July 05, first version
  41 /// - Nov. 05, added DDL
  42 /// - June 06, correction for St123
  43 /// - Feb. 07, add 1st manu list for St12 (starting on NB !) and new ddl sharing for station 3)
  44 /// - June. 07, new numbering of station 345, and buspatch starts at 1)
  45 ///
  46 /// \author Ch. Finck, July 05
  47
  48
  49 void MUONGenerateBusPatch()
  50 {
  51
  52   TString dirPath2 = gSystem->Getenv("ALICE_ROOT");
  53   dirPath2 += "/MUON/mapping/data/";
  54
  55   TString dirPath1 = dirPath2 + "station345/";
  56
  57   TString infile =  dirPath1 + "DetElemIdToSlatType.dat";
  58   TString outfile = dirPath2 + "DetElemIdToBusPatch.dat.out";
  59
  60   ifstream in(infile, ios::in);
  61   ofstream out(outfile, ios::out);
  62   char line[80];
  63   char nameSlat[82];
  64   Int_t idDE;
  65   Int_t i;
  66   Int_t cursor = 0;
  67   Int_t begin[900];
  68   Int_t end[900];
  69   Int_t nbBusPatch;
  70   Int_t nbTB = 0;
  71   Int_t nbBB = 0;
  72   Int_t idSt12[] = {100, 101, 102, 103,
  73                     200, 201, 202, 203,
  74                     300, 301, 302, 303,
  75                     400, 401, 402, 403};
  76
  77   Char_t manuListSt1[]
  78       = " 1,27,53,79,105,131,157,183,201,214,224,232,1025,1051,1077,1103,1129,1155,1181,1207,1225,1238,1249,1257";
  79
  80   Char_t manuListSt2[]
  81       = " 1,27,53,79,105,131,157,183,201,214,226,246,1025,1051,1077,1103,1129,1155,1181,1207,1225,1238,1251,1269";
  82
  83
  84
  85   Int_t iDDL = 0;
  86   // station 1 & 2
  87   nbBusPatch = 24;
  88   Int_t nbHalfBusPatch =  nbBusPatch/2;
  89   cout << "#DE BusPatch DDL SlatName" << endl;
  90   out << "#DE BusPatch DDL  1st manu in buspatch" << endl;
  91
  92   for (Int_t j = 0; j < 16; j++) {
  93
  94     idDE = idSt12[j];
  95     if (idDE % 100 == 0) {
  96       iDDL++;
  97       begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
  98       cout << "# Chamber " << idDE/100 << endl;
  99       out  << "# Chamber " << idDE/100 << endl;
 100     }
 101     if (idDE % 100 == 1) {
 102       iDDL++;
 103       begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
 104     }
 105
 106     if (idDE % 100 == 3) {
 107       iDDL--;
 108       begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1) + nbBusPatch;
 109     }
 110
 111     end[cursor]     = begin[cursor] + nbBusPatch - 1;
 112     ++cursor; // it seems that vec[++i] does not work in Cint
 113     begin[cursor] = end[cursor-1] + 1;
 114
 115     cout << idDE << " " << begin[cursor-1] + nbHalfBusPatch << "-" <<end[cursor-1] << ";" <<
 116         begin[cursor-1] << "-" << end[cursor-1] - nbHalfBusPatch << " " << iDDL-1 << endl;
 117
 118     if (idDE < 300 )
 119         out << idDE  << " " << begin[cursor-1] + nbHalfBusPatch << "-" <<end[cursor-1] << ";" <<
 120             begin[cursor-1]     << "-" << end[cursor-1] - nbHalfBusPatch << " " << iDDL-1 <<
 121             manuListSt1 << endl;
 122     else
 123         out << idDE  << " " << begin[cursor-1] + nbHalfBusPatch << "-" <<end[cursor-1] << ";" <<
 124             begin[cursor-1]     << "-" << end[cursor-1] - nbHalfBusPatch << " " <<  iDDL-1 <<
 125             manuListSt2  << endl;
 126     if (idDE % 100 == 3) iDDL++;
 127
 128   }
 129
 130   // station 345
 131   nbBusPatch = 0;
 132   Int_t idCh5swp1 = 501; // 1/4 chamber for DDL on horizontal
 133   Int_t idCh5swp2 = 505; // 1/4 chamber for DDL on horizontal
 134   Int_t idCh5swp3 = 510; // 1/4 chamber for DDL on horizontal
 135   Int_t idCh5swp4 = 514; // 1/4 chamber for DDL on horizontal
 136
 137   Int_t idCh6swp1 = 600; // 1/4 chamber for DDL on horizontal
 138   Int_t idCh6swp2 = 605; // 1/4 chamber for DDL on horizontal
 139   Int_t idCh6swp3 = 609; // 1/4 chamber for DDL on horizontal
 140   Int_t idCh6swp4 = 614; // 1/4 chamber for DDL on horizontal
 141
 142   Int_t idSt45swp1 = 7;  // half chamber for DDL in vertical cutting twice the official numbering
 143   Int_t idSt45swp2 = 20;
 144
 145   Int_t offsetBusCh5swp1 = 13; // number of buspatches between DE 501-504
 146   Int_t offsetBusCh5swp2 = 17; // number of buspatches before DE 510
 147
 148   Int_t offsetBusCh6swp1 = 17; // number of buspatches before DE 600-604
 149   Int_t offsetBusCh6swp2 = 13; // number of buspatches between DE 614-617
 150
 151   Int_t offsetBusSt4 = 25;     // number of buspatches before DE 700-800
 152   Int_t offsetBusSt4swp = 46;  // number of buspatches between 801-806 (701-706)
 153
 154   Int_t offsetBusSt5 = 27;     // number of buspatches before DE 900-1000
 155   Int_t offsetBusSt5swp = 50;  // number of buspatches between 901-906 (1001-1006)
 156
 157   // reads from file
 158   while ( in.getline(line,80) ) {
 159
 160     if ( line[0] == '#' || line[0] == '\0' ) continue;
 161     i = 0;
 162     sscanf(line, "%d %s", &idDE, nameSlat);
 163
 164     // number of PCB's
 165     while(nameSlat[i] < '4' && nameSlat[i] !='0') i++;
 166     Int_t len = strlen(nameSlat);
 167
 168     // number of buspatch
 169     if (i == 2 || i == 3)
 170       nbBusPatch = 2;
 171     else
 172       nbBusPatch = 4;
 173
 174     // calculate the number of TB & BB
 175     // R1: 2 bridges more per PCB
 176     // R2 & R3: 1 translator more per PCB
 177     if (i == 2) nbBB += 2;
 178     if (i == 3) nbBB += 4;
 179     if (i == 4) nbBB += 4;
 180     if (i == 5) nbBB += 6;
 181     if (i == 6) nbBB += 8;
 182
 183     nbTB += nbBusPatch;
 184     if (nameSlat[len-1] == '1') {
 185       nbTB += 2;
 186       nbBusPatch+=2;
 187     }
 188     if (nameSlat[len-1] == '2' || nameSlat[len-1] == '3') {
 189       nbTB += 1;
 190       nbBusPatch++;
 191     }
 192
 193     // station 3
 194     // for buspatch length reasons, one ddl connects one 1/4 of two chambers
 195     // really messy isn't it ?
 196     // much more with the new DDL sharing for station 3
 197
 198     if (idDE < 700 ) {
 199
 200       if (idDE % 100 == 0) {
 201         cout << "# Chamber " << idDE/100 << endl;
 202         out  << "# Chamber " << idDE/100 << endl;
 203       }
 204
 205       // chamber 5
 206       if (idDE  == 500 ) {
 207         iDDL = (idDE/100)*2+2;
 208         begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
 209
 210         end[cursor]   = begin[cursor] + nbBusPatch - 1;
 211         ++cursor;
 212         begin[cursor] = end[cursor-1] + 1;
 213
 214         cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 215         out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
 216       }
 217
 218       if (idDE  >= idCh5swp1 && idDE < idCh5swp2) {
 219
 220         if ( idDE ==  idCh5swp1) {
 221           iDDL = (idDE/100)*2 ;
 222           begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp1, iDDL-1);
 223         }
 224
 225         end[cursor]   = begin[cursor] - nbBusPatch + 1;
 226         ++cursor;
 227         begin[cursor] = end[cursor-1] - 1;
 228
 229         cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 230         out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
 231       }
 232
 233       if (idDE >=  idCh5swp2  && idDE  < idCh5swp3) {
 234
 235         if (idDE ==  idCh5swp2) {
 236           iDDL = (idDE/100)*2-1;
 237               begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
 238         }
 239
 240         end[cursor]   = begin[cursor] + nbBusPatch - 1;
 241         ++cursor;
 242         begin[cursor] = end[cursor-1] + 1;
 243
 244         cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 245         out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
 246       }
 247
 248       if (idDE >=  idCh5swp3  && idDE  < idCh5swp4) {
 249
 250         if (idDE ==  idCh5swp3) {
 251           iDDL = (idDE/100)*2+1;
 252               begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
 253         }
 254
 255         end[cursor]   = begin[cursor] + nbBusPatch - 1;
 256         ++cursor;
 257         begin[cursor] = end[cursor-1] + 1;
 258
 259         cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 260         out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
 261       }
 262
 263      if (idDE  >= idCh5swp4 && idDE < 600) {
 264
 265         if ( idDE ==  idCh5swp4) {
 266           iDDL = (idDE/100)*2+2 ;
 267           begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp2, iDDL-1);
 268         }
 269
 270         end[cursor]   = begin[cursor] - nbBusPatch + 1;
 271         ++cursor;
 272         begin[cursor] = end[cursor-1] - 1;
 273
 274         cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 275         out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
 276       }
 277
 278      // chamber 6
 279       if (idDE  >= idCh6swp1 && idDE < idCh6swp2) {
 280
 281         if ( idDE ==  idCh6swp1) {
 282           iDDL = (idDE/100)*2-2 ;
 283           begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp1+offsetBusCh6swp1, iDDL-1);
 284         }
 285
 286         end[cursor]   = begin[cursor] - nbBusPatch + 1;
 287         ++cursor;
 288         begin[cursor] = end[cursor-1] - 1;
 289
 290         cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 291         out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
 292       }
 293
 294       if (idDE >=  idCh6swp2  && idDE  < idCh6swp3) {
 295
 296         if (idDE ==  idCh6swp2) {
 297           iDDL = (idDE/100)*2-3;
 298               begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp2+1, iDDL-1);
 299         }
 300
 301         end[cursor]   = begin[cursor] + nbBusPatch - 1;
 302         ++cursor;
 303         begin[cursor] = end[cursor-1] + 1;
 304
 305         cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 306         out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
 307       }
 308
 309       if (idDE >=  idCh6swp3  && idDE  < idCh6swp4) {
 310
 311         if (idDE ==  idCh6swp3) {
 312           iDDL = (idDE/100)*2-1;
 313               begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp1+1, iDDL-1);
 314         }
 315
 316         end[cursor]   = begin[cursor] + nbBusPatch - 1;
 317         ++cursor;
 318         begin[cursor] = end[cursor-1] + 1;
 319
 320         cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 321         out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
 322       }
 323
 324      if (idDE  >= idCh6swp4 && idDE < 700) {
 325
 326         if ( idDE ==  idCh6swp4) {
 327           iDDL = (idDE/100)*2 ;
 328           begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp2+offsetBusCh6swp2, iDDL-1);
 329         }
 330
 331         end[cursor]   = begin[cursor] - nbBusPatch + 1;
 332         ++cursor;
 333         begin[cursor] = end[cursor-1] - 1;
 334
 335         cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 336         out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
 337       }
 338
 339     }
 340
 341     // station 4 & 5
 342     // back to normal one ddl connects one 1/2 chamber
 343     // finally !
 344
 345     if (idDE >= 700) {
 346
 347       if (idDE % 100 == 0) {
 348         cout << "# Chamber " << idDE/100 << endl;
 349         out  << "# Chamber " << idDE/100 << endl;
 350       }
 351
 352       if ((idDE % 100) >= idSt45swp1 && (idDE % 100) < idSt45swp2) {
 353         if ((idDE % 100) == idSt45swp1) {
 354           iDDL = (idDE/100)*2-1;
 355           begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
 356         }
 357         end[cursor]   = begin[cursor] + nbBusPatch - 1;
 358         ++cursor;
 359         begin[cursor] = end[cursor-1] + 1;
 360
 361         cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 362         out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
 363       }
 364
 365       if (idDE % 100 >= idSt45swp2) {
 366
 367         if ((idDE % 100) == idSt45swp2 ) {
 368           iDDL = (idDE/100)*2 ;
 369           if (idDE/100 == 7 || idDE/100 == 8)
 370               begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt4swp, iDDL-1);
 371           else
 372               begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt5swp, iDDL-1);
 373
 374         }
 375         end[cursor]   = begin[cursor] - nbBusPatch + 1;
 376         ++cursor;
 377         begin[cursor] = end[cursor-1] - 1;
 378
 379         cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 380         out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
 381       }
 382
 383       if ((idDE % 100) >= 0 && (idDE % 100) < idSt45swp1) {
 384
 385         if ((idDE % 100) == 0 ) {
 386           iDDL = (idDE/100)*2;
 387           if (idDE/100 == 7 || idDE/100 == 8)
 388               begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt4, iDDL-1);
 389           else
 390               begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt5, iDDL-1);
 391
 392         }
 393         end[cursor]   = begin[cursor] - nbBusPatch + 1;
 394         ++cursor;
 395         begin[cursor] = end[cursor-1] - 1;
 396
 397         cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
 398         out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
 399       }
 400
 401     }
 402
 403   }
 404   printf(" Slat: number of TB %d and BB %d\n", nbTB, nbBB);
 405   in.close();
 406   out.close();
 407 }