Major upgrade (P. Cerello)

[u/mrichter/AliRoot.git] / ITS / AliITSQASDDDataMakerSim.cxx

/**************************************************************************\r
 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *\r
 *                                                                        *\r
 * Author: The ALICE Off-line Project.                                    *\r
 * Contributors are mentioned in the code where appropriate.              *\r
 *                                                                        *\r
 * Permission to use, copy, modify and distribute this software and its   *\r
 * documentation strictly for non-commercial purposes is hereby granted   *\r
 * without fee, provided that the above copyright notice appears in all   *\r
 * copies and that both the copyright notice and this permission notice   *\r
 * appear in the supporting documentation. The authors make no claims     *\r
 * about the suitability of this software for any purpose. It is          *\r
 * provided "as is" without express or implied warranty.                  *\r
 **************************************************************************/\r
\r
/* $Id$ */\r
\r
//  *************************************************************\r
//  Checks the quality assurance \r
//  by comparing with reference data\r
//  contained in a DB\r
//  -------------------------------------------------------------\r
//  W. Ferrarese + P. Cerello Feb 2008\r
//  INFN Torino\r
\r
// --- ROOT system ---\r
#include <TTree.h>\r
// --- Standard library ---\r
\r
// --- AliRoot header files ---\r
#include "AliITSQASDDDataMakerSim.h"\r
#include "AliLog.h"\r
#include "AliQA.h"\r
#include "AliQAChecker.h"\r
#include "AliQADataMakerSim.h"\r
#include "AliITSQADataMakerSim.h"\r
#include "AliRawReader.h"\r
#include "AliITSdigit.h"\r
#include "AliITS.h"\r
#include "AliITSmodule.h"\r
#include "AliITShit.h"\r
#include "AliITSLoader.h"\r
#include "AliRunLoader.h"\r
#include "AliRun.h"\r
#include "AliITSsegmentationSDD.h"\r
#include "AliITSpList.h"\r
\r
ClassImp(AliITSQASDDDataMakerSim)\r
\r
//____________________________________________________________________________ \r
AliITSQASDDDataMakerSim::AliITSQASDDDataMakerSim(AliITSQADataMakerSim *aliITSQADataMakerSim) :\r
TObject(),\r
fAliITSQADataMakerSim(aliITSQADataMakerSim),\r
fSDDhDigits(0),\r
fSDDhSDigits(0),\r
fSDDhHits(0),\r
fDigitsOffset(0),\r
fSDigitsOffset(0),\r
fHitsOffset(0)\r
{\r
  //ctor used to discriminate OnLine-Offline analysis   \r
}\r
\r
//____________________________________________________________________________ \r
AliITSQASDDDataMakerSim::AliITSQASDDDataMakerSim(const AliITSQASDDDataMakerSim& qadm) :\r
TObject(),\r
fAliITSQADataMakerSim(qadm.fAliITSQADataMakerSim),\r
fSDDhDigits(qadm.fSDDhDigits),\r
fSDDhSDigits(qadm.fSDDhSDigits),\r
fSDDhHits(qadm.fSDDhHits),\r
fDigitsOffset(qadm.fDigitsOffset),\r
fSDigitsOffset(qadm.fSDigitsOffset),\r
fHitsOffset(qadm.fHitsOffset)\r
{\r
  //copy ctor \r
  fAliITSQADataMakerSim->SetName((const char*)qadm.fAliITSQADataMakerSim->GetName()) ; \r
  fAliITSQADataMakerSim->SetTitle((const char*)qadm.fAliITSQADataMakerSim->GetTitle());\r
  }\r
\r
//__________________________________________________________________\r
AliITSQASDDDataMakerSim& AliITSQASDDDataMakerSim::operator = (const AliITSQASDDDataMakerSim& qac )\r
{\r
  // Equal operator.\r
  this->~AliITSQASDDDataMakerSim();\r
  new(this) AliITSQASDDDataMakerSim(qac);\r
  return *this;\r
}\r
\r
//____________________________________________________________________________ \r
void AliITSQASDDDataMakerSim::StartOfDetectorCycle()\r
{\r
  //Detector specific actions at start of cycle\r
  AliDebug(1,"AliITSQADM::Start of SDD Cycle\n");\r
}\r
\r
//____________________________________________________________________________ \r
void AliITSQASDDDataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t /*task*/, TObjArray* /*list*/)\r
{\r
  // launch the QA checking\r
  AliDebug(1,"AliITSDM instantiates checker with Run(AliQA::kITS, task, list)\n"); \r
  //AliQAChecker::Instance()->Run( AliQA::kITS , task, list);\r
}\r
\r
//____________________________________________________________________________ \r
void AliITSQASDDDataMakerSim::InitDigits()\r
{ \r
  // Initialization for DIGIT data - SDD -\r
  fDigitsOffset = (fAliITSQADataMakerSim->fDigitsQAList)->GetEntries();\r
  //fSDDhDigits must be incremented by one unit every time a histogram is ADDED to the QA List\r
  //printf("AliITSQASDDDataMakerSim::InitDigits called \n");\r
  TH1F* h0=new TH1F("SDD DIGITS Module Pattern","SDD DIGITS Module Pattern",260,239.5,499.5);       //hmod\r
  h0->GetXaxis()->SetTitle("SDD Module Number");\r
  h0->GetYaxis()->SetTitle("# DIGITS");\r
  fAliITSQADataMakerSim->Add2DigitsList(h0,1+fDigitsOffset);\r
  fSDDhDigits ++;\r
  TH1F* h1=new TH1F("SDD Anode Distribution","DIGITS Anode Distribution",512,-0.5,511.5);      //hanocc\r
  h1->GetXaxis()->SetTitle("Anode Number");\r
  h1->GetYaxis()->SetTitle("# DIGITS");\r
  fAliITSQADataMakerSim->Add2DigitsList(h1,2+fDigitsOffset);\r
  fSDDhDigits ++;\r
  TH1F* h2=new TH1F("SDD Tbin Distribution","DIGITS Tbin Distribution",256,-0.5,255.5);      //htbocc\r
  h2->GetXaxis()->SetTitle("Tbin Number");\r
  h2->GetYaxis()->SetTitle("# DIGITS");\r
  fAliITSQADataMakerSim->Add2DigitsList(h2,3+fDigitsOffset);\r
  fSDDhDigits ++;\r
  TH1F* h3=new TH1F("SDD ADC Counts Distribution","DIGITS ADC Counts Distribution",200,0.,1024.);          //hsig\r
  h3->GetXaxis()->SetTitle("ADC Value");\r
  h3->GetYaxis()->SetTitle("# DIGITS");\r
  fAliITSQADataMakerSim->Add2DigitsList(h3,4+fDigitsOffset);\r
  fSDDhDigits ++;\r
  AliDebug(1,Form("%d SDD Digits histograms booked\n",fSDDhDigits));\r
}\r
\r
//____________________________________________________________________________\r
void AliITSQASDDDataMakerSim::MakeDigits(TTree * digits)\r
{ \r
  // Fill QA for DIGIT - SDD -\r
  //  printf("AliITSQASDDDataMakerSim::MakeDigits called \n");\r
  AliITS *fITS  = (AliITS*)gAlice->GetModule("ITS");\r
  fITS->SetTreeAddress();\r
  TClonesArray *iITSdigits  = fITS->DigitsAddress(1);\r
  for(Int_t i=0; i<260; i++){\r
    Int_t nmod=i+240;\r
    digits->GetEvent(nmod);\r
    Int_t ndigits = iITSdigits->GetEntries();\r
    fAliITSQADataMakerSim->GetDigitsData(1+fDigitsOffset)->Fill(nmod,ndigits);\r
    for (Int_t idig=0; idig<ndigits; idig++) {\r
      AliITSdigit *dig=(AliITSdigit*)iITSdigits->UncheckedAt(idig);\r
      Int_t iz=dig->GetCoord1();  // cell number z\r
      Int_t ix=dig->GetCoord2();  // cell number x\r
      Int_t sig=dig->GetSignal();\r
      fAliITSQADataMakerSim->GetDigitsData(2+fDigitsOffset)->Fill(iz);\r
      fAliITSQADataMakerSim->GetDigitsData(3+fDigitsOffset)->Fill(ix);\r
      fAliITSQADataMakerSim->GetDigitsData(4+fDigitsOffset)->Fill(sig);\r
    }\r
  }\r
}\r
\r
//____________________________________________________________________________ \r
void AliITSQASDDDataMakerSim::InitSDigits()\r
{ \r
  // Initialization for SDIGIT data - SDD -\r
  fSDigitsOffset = (fAliITSQADataMakerSim->fSDigitsQAList)->GetEntries();\r
  //printf("AliITSQASDDDataMakerSim::InitSDigits called \n");\r
  //fSDDhSDigits must be incremented by one unit every time a histogram is ADDED to the QA List\r
  TH1F* h0=new TH1F("SDD SDIGITS Module Pattern","SDIGITS SDD Module Pattern",260,239.5,499.5);       //hmod\r
  h0->GetXaxis()->SetTitle("SDD Module Number");\r
  h0->GetYaxis()->SetTitle("# SDIGITS");\r
  fAliITSQADataMakerSim->Add2SDigitsList(h0,1+fSDigitsOffset);\r
  fSDDhSDigits ++;\r
  TH1F* h1=new TH1F("SDD Anode Distribution","SDIGITS Anode Distribution",512,-0.5,511.5);      //hanocc\r
  h1->GetXaxis()->SetTitle("Anode Number");\r
  h1->GetYaxis()->SetTitle("# SDIGITS");\r
  fAliITSQADataMakerSim->Add2SDigitsList(h1,2+fSDigitsOffset);\r
  fSDDhSDigits ++;\r
  TH1F* h2=new TH1F("SDD Tbin Distribution","SDIGITS Tbin Distribution",256,-0.5,255.5);      //htbocc\r
  h2->GetXaxis()->SetTitle("Tbin Number");\r
  h2->GetYaxis()->SetTitle("# SDIGITS");\r
  fAliITSQADataMakerSim->Add2SDigitsList(h2,3+fSDigitsOffset);\r
  fSDDhSDigits ++;\r
  TH1F* h3=new TH1F("SDD ADC Counts Distribution","SDIGITS ADC Counts Distribution",200,0.,1024.);          //hsig\r
  h3->GetXaxis()->SetTitle("ADC Value");\r
  h3->GetYaxis()->SetTitle("# SDIGITS");\r
  fAliITSQADataMakerSim->Add2SDigitsList(h3,4+fSDigitsOffset);\r
  fSDDhSDigits ++;\r
\r
  AliDebug(1,Form("%d SDD SDigits histograms booked\n",fSDDhSDigits));\r
}\r
\r
//____________________________________________________________________________\r
void AliITSQASDDDataMakerSim::MakeSDigits(TTree * sdigits)\r
{ \r
  // Fill QA for SDIGIT - SDD -\r
  //  printf("AliITSQASDDDataMakerSim::MakeSDigits called \n");\r
  AliITSsegmentationSDD* seg = new AliITSsegmentationSDD();\r
  Int_t nan=seg->Npz();\r
  Int_t ntb=seg->Npx();\r
  Int_t scaleSize=4;\r
  AliITSpList* list=new AliITSpList(nan,ntb*scaleSize);\r
\r
  //AliITS *fITS  = (AliITS*)gAlice->GetModule("ITS");\r
  //fITS->SetTreeAddress();\r
  //TClonesArray *ITSdigits  = fITS->DigitsAddress(1);\r
  //TFile *sper = new TFile("sper.root","CREATE"); //agginto a mano x prova\r
  //digits->Write();\r
  //sper->Close();\r
\r
\r
  TBranch *brchSDigits = sdigits->GetBranch("ITS");\r
  for(Int_t id=0; id<260; id++){\r
    Int_t nmod=id+240;\r
    TClonesArray * sdig = new TClonesArray( "AliITSpListItem",1000 );\r
    brchSDigits->SetAddress( &sdig );\r
    brchSDigits->GetEvent(nmod);\r
    Int_t nsdig=sdig->GetEntries();\r
    fAliITSQADataMakerSim->GetSDigitsData(1+fSDigitsOffset)->Fill(nmod,nsdig);\r
    for(Int_t i=0;i<nsdig;i++){\r
      AliITSpListItem *cell=(AliITSpListItem*)sdig->At(i);\r
      Float_t sig=cell->GetSignal();\r
      Int_t idx=cell->GetIndex();\r
      Int_t ia,it;\r
      list->GetCell(idx,ia,it);\r
      fAliITSQADataMakerSim->GetSDigitsData(2+fSDigitsOffset)->Fill(ia);\r
      fAliITSQADataMakerSim->GetSDigitsData(3+fSDigitsOffset)->Fill(it);\r
      fAliITSQADataMakerSim->GetSDigitsData(4+fSDigitsOffset)->Fill(sig);\r
    }\r
    delete sdig;\r
  }\r
}\r
\r
//____________________________________________________________________________ \r
void AliITSQASDDDataMakerSim::InitHits()\r
{ \r
  // Initialization for HITS data - SDD -\r
  fHitsOffset = (fAliITSQADataMakerSim->fHitsQAList)->GetEntries();\r
  //fSDDhHits must be incremented by one unit every time a histogram is ADDED to the QA List\r
  //printf("AliITSQASDDDataMakerSim::InitHits called \n");\r
  TH1F *h0=new TH1F("SDD HITS Module Pattern","SDD HITS Module Pattern",260,239.5,499.5);  \r
  h0->GetXaxis()->SetTitle("SDD Module Number");\r
  h0->GetYaxis()->SetTitle("# HITS");\r
  fAliITSQADataMakerSim->Add2HitsList(h0,1+fHitsOffset);\r
  fSDDhHits ++;\r
  TH1F *h1=new TH1F("SDD HIT lenght along local Y Coord","HIT lenght along local Y Coord",200,0.,350.);\r
  h1->GetXaxis()->SetTitle("HIT lenght (um)");\r
  h1->GetYaxis()->SetTitle("# HITS");\r
  fAliITSQADataMakerSim->Add2HitsList(h1,2+fHitsOffset);\r
  fSDDhHits ++;\r
  TH1F *h2=new TH1F("SDD HIT lenght along local Y Coord - Zoom","SDD HIT lenght along local Y Coord - Zoom",200,250.,350.);\r
  h2->GetXaxis()->SetTitle("HIT lenght (um)");\r
  h2->GetYaxis()->SetTitle("# HITS");\r
  fAliITSQADataMakerSim->Add2HitsList(h2,3+fHitsOffset);\r
  fSDDhHits ++;\r
  TH1F *h3=new TH1F("SDD Deposited Energy Distribution (loc Y > 200um)","SDD HITS Deposited Energy Distribution (loc Y > 200um)",200,0.,350.);\r
  h3->GetXaxis()->SetTitle("ADC counts???");\r
  h3->GetYaxis()->SetTitle("# HITS");\r
  fAliITSQADataMakerSim->Add2HitsList(h3,4+fHitsOffset);\r
  fSDDhHits ++;\r
  //printf("%d SDD Hits histograms booked\n",fSDDhHits);\r
  AliDebug(1,Form("%d SDD Hits histograms booked\n",fSDDhHits));\r
}\r
\r
//____________________________________________________________________________\r
void AliITSQASDDDataMakerSim::MakeHits(TTree * hits)\r
{ \r
  // Fill QA for HITS - SDD -\r
  //printf("AliITSQASDDDataMakerSim::MakeHits called \n");\r
  AliITS *fITS  = (AliITS*)gAlice->GetModule("ITS");\r
  fITS->SetTreeAddress();\r
  Int_t nmodules;\r
  if(!(fITS->InitModules(-1,nmodules))){\r
    AliError("ITS geometry not available - nothing done");\r
    return;\r
  }\r
 \r
  fITS->FillModules(hits,0);\r
\r
  for(Int_t i=0; i<260; i++){\r
    Int_t nmod=i+240;\r
    AliITSmodule *modu = fITS->GetModule(nmod);\r
    TObjArray *arrHits = modu->GetHits();\r
    Int_t nhits = arrHits->GetEntriesFast();\r
    //printf("--w--AliITSQASDDDataMakerSim::MakeHits  nhits = %d\n",nhits);\r
    for (Int_t iHit=0;iHit<nhits;iHit++) {\r
      AliITShit *hit = (AliITShit*) arrHits->At(iHit);\r
      fAliITSQADataMakerSim->GetHitsData(1 + fHitsOffset)->Fill(nmod);\r
      Double_t xl,yl,zl,xl0,yl0,zl0;\r
      Double_t tof,tof0;\r
      hit->GetPositionL(xl,yl,zl,tof);\r
      hit->GetPositionL0(xl0,yl0,zl0,tof0);\r
      Float_t dyloc=TMath::Abs(yl-yl0)*10000.;\r
      fAliITSQADataMakerSim->GetHitsData(2 + fHitsOffset)->Fill(dyloc);\r
      Float_t edep=hit->GetIonization()*1000000;\r
      if(dyloc>200.){ \r
	fAliITSQADataMakerSim->GetHitsData(4 + fHitsOffset)->Fill(edep);\r
	fAliITSQADataMakerSim->GetHitsData(3 + fHitsOffset)->Fill(dyloc);\r
      }\r
    }\r
  }\r
}\r