ITS/ITSSDDQA.C

   1 void ITSSDDQA(char *iFile, Int_t MaxEvts=1000000, Int_t FirstEvt=0) {
   2
   3 cout << "SDD Quality Assurance Prototype Macro" << endl;
   4
   5 const Int_t nSDDmodules= 260;
   6 const Int_t imodoffset = 240;
   7 const Int_t modtotSDD  = nSDDmodules*2;
   8
   9 Float_t xi = 0.5;
  10 Float_t xf = xi + nSDDmodules;
  11 TH1F *modulePattern = new TH1F("patternModule","Modules pattern",nSDDmodules,xi,xf);
  12 xf = xi + modtotSDD;
  13 TH1F *moduleSidePattern = new TH1F("patternSide","Modules/Side pattern",modtotSDD,xi,xf);
  14
  15 TH2F *hismap[modtotSDD];  //260 dx e 260 sx  with A, T, Q
  16 TH2F *hispop[modtotSDD];  //260 dx e 260 sx  with A, T, Ncounts
  17
  18 Char_t *cindex = new Char_t[5];
  19 for(Int_t imod=0; imod<nSDDmodules;imod++){
  20   for(Int_t isid=0;isid<2;isid++){
  21     Int_t index=2*imod+isid;       //260*2 position
  22
  23     sprintf(cindex,"%d",index+1); // imod,isid);
  24     TString sindex((const char *) cindex);
  25
  26     TString histnam("chargeMap");
  27     TString histit("Total Charge, module number ");
  28     histnam.Append(sindex);
  29     histit.Append(sindex);
  30     hismap[index]=new TH2F(histnam.Data(),histit.Data(),256,-0.5,255.5,256,-0.5,255.5);
  31
  32     TString hisnam2("countsMap");
  33     TString histit2("Number of Counts, module number ");
  34     hisnam2.Append(sindex);
  35     histit2.Append(sindex);
  36     hispop[index]=new TH2F(hisnam2.Data(),histit2.Data(),256,-0.5,255.5,256,-0.5,255.5);
  37
  38     /*
  39       sprintf(hisnam,"hisprojX%03ds%d",imod,isid);
  40       sprintf(histitle,"layer , ladder, module position, channel, %d, %d", imod, isid);
  41       hisprojX[index]=new TH2F(hisnam,histitle,256,-0.5,255.5,256,-0.5,255.5);
  42
  43       sprintf(hisnam,"hisprojY%03ds%d",imod,isid);
  44       sprintf(histitle,"layer , ladder, module position, channel, %d, %d", imod, isid);
  45       hisprojY[index]=new TH2F(hisnam,histitle,256,-0.5,255.5,256,-0.5,255.5);
  46
  47       sprintf(hisnam,"hisprofX%03ds%d",imod,isid);
  48       sprintf(histitle,"layer , ladder, module position, channel, %d, %d", imod, isid);
  49       hisprofX[index]=new TH2F(hisnam,histitle,256,-0.5,255.5,256,-0.5,255.5);
  50
  51       sprintf(hisnam,"hisprofY%03ds%d",imod,isid);
  52       sprintf(histitle,"layer , ladder, module position, channel, %d, %d", imod, isid);
  53       hisprofY[index]=new TH2F(hisnam,histitle,256,-0.5,255.5,256,-0.5,255.5);
  54       */
  55
  56   }
  57 }
  58 delete [] cindex;
  59
  60   TString strFile = iFile;
  61   strFile += "?EventType=7";
  62   AliRawReader *rd = new AliRawReaderDate(strFile.Data(),FirstEvt);  // open run
  63   Int_t evCounter = 0;
  64
  65   //AliITS *itsRun = new AliITS();
  66   //  TGeoManager::Import("$ALICE_ROOT/EVE/alice-data/alice_fullgeo.root");
  67   /*
  68     AliITSInitGeometry *initgeom = new AliITSInitGeometry("AliITSvPPRasymmFMD",2);
  69     geom = initgeom->CreateAliITSgeom();
  70     delete initgeom;
  71   */
  72
  73   do{       // start loop on events
  74     if(++evCounter > MaxEvts) { cout << MaxEvts << " events read, stop" << endl; evCounter--; break; }
  75      cout << "Read Event: " << evCounter+FirstEvt-1 << endl;
  76
  77     rd->RequireHeader(kFALSE);
  78
  79
  80     rd->Reset();    // reset the current position to the beginning of the event
  81     AliITSRawStreamSDD s(rd);    //This class provides access to ITS SDD digits in raw data.
  82     Int_t iddl;
  83     Int_t isddmod;
  84     Int_t moduleSDD;
  85     while(s.Next()){                       //read the next raw digit; returns kFALSE if there is no digit left
  86
  87         iddl=rd->GetDDLID();
  88         isddmod=s.GetModuleNumber(iddl,s.GetCarlosId());  //this is the FEE Carlos
  89         //cout<<"DDLID= "<<iddl <<"; Module number= " <<isddmod  <<endl;
  90         modulePattern->Fill(isddmod-imodoffset+1);     // 1 to 260
  91         moduleSDD=2*(isddmod-imodoffset)+s.GetChannel();
  92         moduleSidePattern->Fill(moduleSDD+1);          // 1 to 520
  93         //cout << "anode " << s.GetCoord1() << ", time bin: " << s.GetCoord2() << ", charge: " << s.GetSignal() << endl;
  94         hismap[moduleSDD]->Fill(s.GetCoord1(), s.GetCoord2(),s.GetSignal() );
  95         hispop[moduleSDD]->Fill(s.GetCoord1(), s.GetCoord2() );
  96     }
  97
  98   } while(rd->NextEvent()); // end loop on events
  99   delete rd;
 100
 101   cout << "end after " << evCounter << " events" << endl;
 102
 103    TString oFileName(iFile);
 104    oFileName.Append(".root");
 105    TFile *oFile = TFile::Open(oFileName,"recreate");
 106    modulePattern->Write();
 107    moduleSidePattern->Write();
 108    for(Int_t i=0; i<modtotSDD; i++){
 109      hismap[i]->Write();
 110      hispop[i]->Write();
 111    }
 112    oFile->Close();
 113
 114    for(Int_t imod=0; imod<nSDDmodules;imod++){
 115      for(Int_t isid=0;isid<2;isid++){
 116        Int_t index=2*imod+isid;       //260*2 position
 117        delete hismap[index];
 118        delete hispop[index];
 119      }
 120    }
 121    delete modulePattern;
 122    delete moduleSidePattern;
 123
 124 }