* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-Revision 1.9 1999/11/14 14:33:25 fca
-Correct problems with distructors and pointers, thanks to I.Hrivnacova
-
-Revision 1.8 1999/09/29 09:24:19 fca
-Introduction of the Copyright and cvs Log
+/* $Id$ */
-*/
///////////////////////////////////////////////////////////////////////////////
-//
-// An overview of the basic philosophy of the ITS code development
-// and analysis is show in the figure below.
-//Begin_Html
-/*
+// //
+// An overview of the basic philosophy of the ITS code development //
+// and analysis is show in the figure below. //
+//Begin_Html //
+/*
<img src="picts/ITS/ITS_Analysis_schema.gif">
</pre>
<br clear=left>
// Modified and documented by Bjorn S. Nilsen
// July 11 1999
//
+// Version: 2
+// Modified and documented by A. Bologna
+// October 18 1999
+//
// AliITS is the general base class for the ITS. Also see AliDetector for
// futher information.
//
///////////////////////////////////////////////////////////////////////////////
-
-#include <TMath.h>
-#include <TRandom.h>
-#include <TVector.h>
-#include <TGeometry.h>
-#include <TNode.h>
-#include <TTUBE.h>
-#include "AliITSmodule.h"
+#include <stdlib.h>
+#include <TClonesArray.h>
+#include <TFile.h>
+#include <TParticle.h>
+#include <TString.h>
+#include <TTree.h>
+#include <TVirtualMC.h>
#include "AliDetector.h"
#include "AliITS.h"
-#include "TClonesArray.h"
-#include "TObjArray.h"
+#include "AliITSDetTypeSim.h"
+#include "AliITSDDLRawData.h"
+#include "AliITSLoader.h"
#include "AliITShit.h"
-#include "AliITSdigit.h"
+#include "AliITSmodule.h"
+#include "AliITSpListItem.h"
+#include "AliITSsimulation.h"
+#include "AliITSsimulationFastPoints.h"
+#include "AliMC.h"
+#include "AliITSDigitizer.h"
+#include "AliITSRecPoint.h"
+#include "AliITSsegmentationSPD.h"
+#include "AliITSsegmentationSDD.h"
+#include "AliITSsimulationSDD.h"
+#include "AliITSCalibrationSDD.h"
+#include "AliITSCalibrationSSD.h"
+#include "AliITSsegmentationSSD.h"
+#include "AliITSRawStreamSPD.h"
+#include "AliITSRawStreamSSD.h"
+#include "AliITSRawStreamSDD.h"
+#include "AliRawReader.h"
#include "AliRun.h"
+#include "AliLog.h"
+#include "AliITSInitGeometry.h"
+#include "AliITSFOSignalsSPD.h"
ClassImp(AliITS)
-
-//_____________________________________________________________________________
-AliITS::AliITS() {
- //
- // Default initialiser for ITS
- // The default constructor of the AliITS class. In addition to
+
+//______________________________________________________________________
+AliITS::AliITS() : AliDetector(),
+fDetTypeSim(0),
+fEuclidOut(0),
+fOpt("All"),
+fIdN(0),
+fIdSens(0),
+fIdName(0),
+fITSmodules(0),
+fTiming(kFALSE),
+fSimuParam(0),
+fModA(0),
+fpSDigits(0)
+{
+ // Default initializer for ITS
+ // The default constructor of the AliITS class. In addition to
// creating the AliITS class it zeros the variables fIshunt (a member
// of AliDetector class), fEuclidOut, and fIdN, and zeros the pointers
- // fITSpoints, fIdSens, and fIdName.
- //
- fITSpoints = 0;
- fIshunt = 0;
- fEuclidOut = 0;
- fIdN = 0;
- fIdName = 0;
- fIdSens = 0;
- fITSmodules = 0;
+ // fIdSens, and fIdName. The AliDetector default constructor
+ // is also called.
+
+// SetDetectors(); // default to fOpt="All". This variable not written out.
+//PH SetMarkerColor(kRed);
+}
+//______________________________________________________________________
+AliITS::AliITS(const Char_t *title):
+ AliDetector("ITS",title),
+ fDetTypeSim(0),
+ fEuclidOut(0),
+ fOpt("All"),
+ fIdN(0),
+ fIdSens(0),
+ fIdName(0),
+ fITSmodules(0),
+ fTiming(kFALSE),
+ fSimuParam(0),
+ fModA(0),
+ fpSDigits(0)
+{
+ // The standard Constructor for the ITS class.
+ // It also zeros the variables
+ // fIshunt (a member of AliDetector class), fEuclidOut, and zeros
+ // the pointers fIdSens and fIdName. To help in displaying hits via the
+ // ROOT macro display.C AliITS also sets the marker color to red. The
+ // variables passes with this constructor, const char *name and *title,
+ // are used by the constructor of AliDetector class. See AliDetector
+ // class for a description of these parameters and its constructor
+ // functions.
+ // Inputs:
+ // Char_t *title Simulation title for the ITS
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+
+ fHits = new TClonesArray("AliITShit",1560); // from AliDetector
+ if(gAlice->GetMCApp()) gAlice->GetMCApp()->AddHitList(fHits);
+ //fNhits=0; //done in AliDetector(name,title)
+ SetDetectors(); // default to fOpt="All". This variable not written out.
+ fDetTypeSim = new AliITSDetTypeSim();
+ //PH SetMarkerColor(kRed);
+ if(!fLoader) MakeLoader(AliConfig::GetDefaultEventFolderName());
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+}
+//______________________________________________________________________
+AliITS::AliITS(const char *name, const char *title):
+ AliDetector(name,title),
+ fDetTypeSim(0),
+ fEuclidOut(0),
+ fOpt("All"),
+ fIdN(0),
+ fIdSens(0),
+ fIdName(0),
+ fITSmodules(0),
+ fTiming(kFALSE),
+ fSimuParam(0),
+ fModA(0),
+ fpSDigits(0)
+{
+ // The standard Constructor for the ITS class.
+ // It also zeros the variables
+ // fIshunt (a member of AliDetector class), fEuclidOut, and zeros
+ // the pointers fIdSens and fIdName. To help in displaying hits via the
+ // ROOT macro display.C AliITS also sets the marker color to red. The
+ // variables passes with this constructor, const char *name and *title,
+ // are used by the constructor of AliDetector class. See AliDetector
+ // class for a description of these parameters and its constructor
+ // functions.
+
+ fHits = new TClonesArray("AliITShit",1560);
+ if(gAlice->GetMCApp()) gAlice->GetMCApp()->AddHitList(fHits);
+ //fNhits=0; //done in AliDetector(name,title)
+
+ SetDetectors(); // default to fOpt="All". This variable not written out.
+
+ fDetTypeSim = new AliITSDetTypeSim();
+ //PH SetMarkerColor(kRed);
+ if(!fLoader) MakeLoader(AliConfig::GetDefaultEventFolderName());
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
}
-//_____________________________________________________________________________
-AliITS::AliITS(const char *name, const char *title):AliDetector(name,title){
- //
- // Default initialiser for ITS
- // The constructor of the AliITS class. In addition to creating the
- // AliITS class, it allocates memory for the TClonesArrays fHits and
- // fDigits, and for the TObjArray fITSpoints. It also zeros the variables
- // fIshunt (a member of AliDetector class), fEuclidOut, and fIdN, and zeros
- // the pointers fIdSens and fIdName. To help in displaying hits via the ROOT
- // macro display.C AliITS also sets the marker color to red. The variables
- // passes with this constructor, const char *name and *title, are used by
- // the constructor of AliDetector class. See AliDetector class for a
- // description of these parameters and its constructor functions.
- //
+//______________________________________________________________________
+AliITS::~AliITS(){
+ // Default destructor for ITS.
+ // The default destructor of the AliITS class. In addition to deleting
+ // the AliITS class it deletes the memory pointed to by
+ // fIdSens, fIdName, fDetTypeSim and it's contents.
+ // Inputs:
+ // none.
+ // Outputs:
+ // none.
+ // Return:
+ // none.
- fHits = new TClonesArray("AliITShit", 1560);
- gAlice->AddHitList(fHits);
- fDigits = new TClonesArray("AliITSdigit",1000);
- fITSpoints = new TObjArray();
- fITSmodules = 0; //new AliITSmodules();
+ if (fHits) {
+ fHits->Delete();
+ delete fHits;
+ fHits=0;
+ }
+ if(fITSmodules) {
+ this->ClearModules();
+ delete fITSmodules;
+ fITSmodules = 0;
+ }// end if fITSmodules!=0
- fIshunt = 0;
- fEuclidOut = 0;
- fIdN = 0;
- fIdName = 0;
- fIdSens = 0;
+ delete[] fIdName; // Array of TStrings
+ delete[] fIdSens;
+ Int_t size = AliITSgeomTGeo::GetNModules();
+ if (fDetTypeSim){
+ delete fDetTypeSim;
+ fDetTypeSim = 0;
+ }
+ if(fSimuParam){
+ delete fSimuParam;
+ fSimuParam=0;
+ }
+ if(fModA){
+ if(size>0){
+ for(Int_t j=0; j<size; j++){
+ fModA[j]->Delete();
+ delete fModA[j];
+ }
+ }
+ delete []fModA;
+ }
+ if(fpSDigits){
+ fpSDigits->Delete();
+ delete fpSDigits;
+ }
+}
+//______________________________________________________________________
+AliDigitizer* AliITS::CreateDigitizer(AliRunDigitizer* manager)const{
+ // Creates the AliITSDigitizer in a standard way for use via AliModule.
+ // This function can not be included in the .h file because of problems
+ // with the order of inclusion (recursive).
+ // Inputs:
+ // AliRunDigitizer *manager The Manger class for Digitization
+ // Output:
+ // none.
+ // Return:
+ // A new AliITSRunDigitizer (cast as a AliDigitizer).
+
+ return new AliITSDigitizer(manager);
+}
+//______________________________________________________________________
+void AliITS::Init(){
+ // Initializer ITS after it has been built
+ // This routine initializes the AliITS class. It is intended to be
+ // called from the Init function in AliITSv?. Besides displaying a banner
+ // indicating that it has been called it initializes the array fIdSens
+ // and sets the default segmentation, response, digit and raw cluster
+ // classes therefore it should be called after a call to CreateGeometry.
+ // Inputs:
+ // none.
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+ Int_t i;
+ // Array of TStrings
+ if(gMC) for(i=0;i<fIdN;i++) fIdSens[i] = gMC->VolId(fIdName[i]);
+
+}
+//______________________________________________________________________
+void AliITS::SetDefaults(){
+ // sets the default segmentation, response, digit and raw cluster classes.
+ // Inputs:
+ // none.
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+ AliInfoClass("AliITS::Setting Defaults");
+ if(!fDetTypeSim) {
+ Error("SetDefaults()","fDetTypeSim is 0!");
+ return;
+ }
- SetMarkerColor(kRed);
+ fDetTypeSim->SetDefaults();
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
}
+//______________________________________________________________________
+void AliITS::SetDefaultSimulation(){
+ // sets the default simulation.
+ // Inputs:
+ // none.
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+ if(!fDetTypeSim) {
+ Error("SetDefaultSimulation()","fDetTypeSim is 0!");
+ return;
+ }
+
+ fDetTypeSim->SetDefaultSimulation();
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
-//_____________________________________________________________________________
-AliITS::~AliITS(){
- //
- // Default distructor for ITS
- // The default destructor of the AliITS class. In addition to deleting
- // the AliITS class it deletes the memory pointed to by the fHits, fDigits,
- // fIdSens, fIdName, and fITSpoints.
- //
- delete fHits;
- delete fDigits;
- if(fIdName!=0) delete[] fIdName;
- if(fIdSens!=0) delete[] fIdSens;
- delete fITSmodules;
- if(fITSpoints!=0) delete fITSpoints;
}
-//_____________________________________________________________________________
-void AliITS::AddDigit(Int_t *tracks, Int_t *digits){
- //
- // Add an ITS Digit
- // The function to add information to the AliITSdigits class. See the
- // AliITSdigits class for a full description. This function allocates the
- // necessary new space for the digits information and passes the pointers
- // *track and *digits to the AliITSdigits constructor function.
- //
- TClonesArray &ldigits = *fDigits;
- new(ldigits[fNdigits++]) AliITSdigit(tracks,digits);
+
+//______________________________________________________________________
+void AliITS::MakeBranch(Option_t* option){
+ // Creates Tree branches for the ITS.
+ // Inputs:
+ // Option_t *option String of Tree types S,D, and/or R.
+ // const char *file String of the file name where these branches
+ // are to be stored. If blank then these branches
+ // are written to the same tree as the Hits were
+ // read from.
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+ if(!fDetTypeSim) {
+ Error("MakeBranch","fDetTypeSim is 0!");
+ return;
+ }
+
+ Bool_t cH = (strstr(option,"H")!=0);
+ Bool_t cS = (strstr(option,"S")!=0);
+ Bool_t cD = (strstr(option,"D")!=0);
+
+ if(cH && (fHits == 0x0)) fHits = new TClonesArray("AliITShit", 1560);
+ AliDetector::MakeBranch(option);
+
+ if(cS) MakeBranchS(0);
+ if(cD) MakeBranchD(0);
+
+
}
+//___________________________________________________________________
+void AliITS::MakeBranchS(const char* fl){
+
+ // Creates Tree Branch for the ITS summable digits.
+ // Inputs:
+ // cont char *fl File name where SDigits branch is to be written
+ // to. If blank it write the SDigits to the same
+ // file in which the Hits were found.
-Int_t AliITS::AddDigit(AliITSdigit* d) {
+
+ if(!fDetTypeSim){
+ Error("MakeBranchS","fDetTypeSim is 0!");
+ }
+ Int_t buffersize = 4000;
+ char branchname[30];
+
+ // only one branch for SDigits.
+ sprintf(branchname,"%s",GetName());
+
+ if(fLoader->TreeS()){
+ TClonesArray* sdig = (TClonesArray*)fDetTypeSim->GetSDigits();
+ MakeBranchInTree(fLoader->TreeS(),branchname,&sdig,buffersize,fl);
+ }
+}
+//______________________________________________________________________
+void AliITS::MakeBranchD(const char* file){
+
+ //Make branch for digits
+ if(!fDetTypeSim) {
+ Warning("MakeBranchD","fDetTypeSim is 0!");
+ return;
+ }
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
+ MakeBranchInTreeD(fLoader->TreeD(),file);
+}
+
+//___________________________________________________________________
+void AliITS:: MakeBranchInTreeD(TTree* treeD, const char* file){
+ // Creates Tree branches for the ITS.
+
+ if(!fDetTypeSim){
+ Error("MakeBranchS","fDetTypeSim is 0!");
+ }
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
+
+ const Char_t *det[3] = {"SPD","SDD","SSD"};
+ const Char_t* digclass;
+ Int_t buffersize = 4000;
+ Char_t branchname[30];
+
+ if(!fDetTypeSim->GetDigits()){
+ fDetTypeSim->SetDigits(new TObjArray(fgkNTYPES));
+ }
+ for(Int_t i=0;i<fgkNTYPES;i++){
+ digclass = fDetTypeSim->GetDigitClassName(i);
+ TString classn = digclass;
+ if(!((fDetTypeSim->GetDigits())->At(i))){
+ (fDetTypeSim->GetDigits())->AddAt(new TClonesArray(classn.Data(),1000),i);
+ }
+ else ResetDigits(i);
+ if(fgkNTYPES==3) sprintf(branchname,"%sDigits%s",GetName(),det[i]);
+ else sprintf(branchname,"%sDigits%d",GetName(),i+1);
+ TObjArray* dig = DigitsAddress(i);
+ if(GetDigits() && treeD) AliDetector::MakeBranchInTree(treeD,branchname, &dig,buffersize,file);
+ }
- fDigits->Add(d);
- fNdigits = fDigits->GetEntriesFast();
- return fNdigits;
}
+//______________________________________________________________________
+void AliITS::SetTreeAddress(){
+ // Set branch address for the Trees.
+ // Inputs:
+ // none.
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+
+ if(!fDetTypeSim) {
+ Error("SetTreeAddress","fDetTypeSim is 0!");
+ return;
+ }
+
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
-//_____________________________________________________________________________
+ TTree *treeS = fLoader->TreeS();
+ TTree *treeD = fLoader->TreeD();
+ if (fLoader->TreeH() && (fHits == 0x0)) {
+ fHits = new TClonesArray("AliITShit", 1560);
+ }
+ AliDetector::SetTreeAddress();
+
+ fDetTypeSim->SetTreeAddressS(treeS, (Char_t*)GetName());
+ fDetTypeSim->SetTreeAddressD(treeD, (Char_t*)GetName());
+}
+//______________________________________________________________________
void AliITS::AddHit(Int_t track, Int_t *vol, Float_t *hits){
- //
- // Add an ITS hit
- // The function to add information to the AliITShit class. See the
- // AliITShit class for a full description. This function allocates the
- // necessary new space for the hit information and passes the variable
- // track, and the pointers *vol and *hits to the AliITShit constructor
- // function.
- //
+ // Add an ITS hit
+ // The function to add information to the AliITShit class. See the
+ // AliITShit class for a full description. This function allocates the
+ // necessary new space for the hit information and passes the variable
+ // track, and the pointers *vol and *hits to the AliITShit constructor
+ // function.
+ // Inputs:
+ // Int_t track Track number which produced this hit.
+ // Int_t *vol Array of Integer Hit information. See AliITShit.h
+ // Float_t *hits Array of Floating Hit information. see AliITShit.h
+ // Outputs:
+ // none.
+ // Return:
+ // none.
TClonesArray &lhits = *fHits;
new(lhits[fNhits++]) AliITShit(fIshunt,track,vol,hits);
}
-//_____________________________________________________________________________
-void AliITS::BuildGeometry(){
- //
- // Build ITS TNODE geometry for event display
- // This function builds a simple ITS geometry used by the ROOT macro
- // display.C. In general the geometry as coded is wrong.
- //
- TNode *Node, *Top;
- const int kColorITS=kYellow;
- //
- Top=gAlice->GetGeometry()->GetNode("alice");
-
- new TTUBE("S_layer1","Layer1 of ITS","void",3.9,3.9+0.05475,12.25);
- Top->cd();
- Node = new TNode("Layer1","Layer1","S_layer1",0,0,0,"");
- Node->SetLineColor(kColorITS);
- fNodes->Add(Node);
-
- new TTUBE("S_layer2","Layer2 of ITS","void",7.6,7.6+0.05475,16.3);
- Top->cd();
- Node = new TNode("Layer2","Layer2","S_layer2",0,0,0,"");
- Node->SetLineColor(kColorITS);
- fNodes->Add(Node);
-
- new TTUBE("S_layer3","Layer3 of ITS","void",14,14+0.05288,21.1);
- Top->cd();
- Node = new TNode("Layer3","Layer3","S_layer3",0,0,0,"");
- Node->SetLineColor(kColorITS);
- fNodes->Add(Node);
-
- new TTUBE("S_layer4","Layer4 of ITS","void",24,24+0.05288,29.6);
- Top->cd();
- Node = new TNode("Layer4","Layer4","S_layer4",0,0,0,"");
- Node->SetLineColor(kColorITS);
- fNodes->Add(Node);
-
- new TTUBE("S_layer5","Layer5 of ITS","void",40,40+0.05382,45.1);
- Top->cd();
- Node = new TNode("Layer5","Layer5","S_layer5",0,0,0,"");
- Node->SetLineColor(kColorITS);
- fNodes->Add(Node);
-
- new TTUBE("S_layer6","Layer6 of ITS","void",45,45+0.05382,50.4);
- Top->cd();
- Node = new TNode("Layer6","Layer6","S_layer6",0,0,0,"");
- Node->SetLineColor(kColorITS);
- fNodes->Add(Node);
-}
-//_____________________________________________________________________________
-void AliITS::CreateMaterials(){
- //
- // Create ITS materials
- // This function defines the default materials used in the Geant
- // Monte Carlo simulations. In general it is automatically replaced by
- // the CreatMaterials routine defined in AliITSv?. Should the function
- // CreateMaterials not exist for the geometry version you are using this
- // one is used. See the definition found in AliITSv5 or the other routine
- // for a complete definition.
- //
- // Water H2O
- Float_t awat[2] = { 1.00794,15.9994 };
- Float_t zwat[2] = { 1.,8. };
- Float_t wwat[2] = { 2.,1. };
- Float_t denswat = 1.;
- // Freon
- Float_t afre[2] = { 12.011,18.9984032 };
- Float_t zfre[2] = { 6.,9. };
- Float_t wfre[2] = { 5.,12. };
- Float_t densfre = 1.5;
- // Ceramics
- // 94.4% Al2O3 , 2.8% SiO2 , 2.3% MnO , 0.5% Cr2O3
- Float_t acer[5] = { 26.981539,15.9994,28.0855,54.93805,51.9961 };
- Float_t zcer[5] = { 13.,8.,14.,25., 24. };
- Float_t wcer[5] = { .49976,1.01233,.01307, .01782,.00342 };
- Float_t denscer = 3.6;
- //
- // 60% SiO2 , 40% G10FR4
- // PC board
- Float_t apcb[3] = { 28.0855,15.9994,17.749 };
- Float_t zpcb[3] = { 14.,8.,8.875 };
- Float_t wpcb[3] = { .28,.32,.4 };
- Float_t denspcb = 1.8;
- // POLYETHYL
- Float_t apoly[2] = { 12.01,1. };
- Float_t zpoly[2] = { 6.,1. };
- Float_t wpoly[2] = { .33,.67 };
- // SERVICES
- Float_t zserv[4] = { 1.,6.,26.,29. };
- Float_t aserv[4] = { 1.,12.,55.8,63.5 };
- Float_t wserv[4] = { .014,.086,.42,.48 };
+
+//______________________________________________________________________
+void AliITS::FillModules(Int_t evnt,Int_t bgrev,Int_t nmodules,
+ Option_t *option, const char *filename){
+ // fill the modules with the sorted by module hits; add hits from
+ // background if option=Add.
+
+ static TTree *trH1; //Tree with background hits
+ static Bool_t first=kTRUE;
+ static TFile *file;
+ const char *addBgr = strstr(option,"Add");
+
+ evnt = nmodules; // Dummy use of variables to remove warnings
+ if (addBgr ) {
+ if(first) {
+ file=new TFile(filename);
+ } // end if first
+ first=kFALSE;
+ file->cd();
+ file->ls();
+ // Get Hits Tree header from file
+ if(trH1) delete trH1;
+ trH1=0;
+
+ char treeName[20];
+ sprintf(treeName,"TreeH%d",bgrev);
+ trH1 = (TTree*)gDirectory->Get(treeName);
+ if (!trH1) {
+ Error("FillModules","cannot find Hits Tree for event:%d",bgrev);
+ } // end if !trH1
+ // Set branch addresses
+ } // end if addBgr
+
+ FillModules(fLoader->TreeH(),0); // fill from this file's tree.
+
+ if (addBgr ) {
+ FillModules(trH1,10000000); // Default mask 10M.
+ TTree *fAli=fLoader->GetRunLoader()->TreeK();
+ TFile *fileAli=0;
+ if (fAli) fileAli =fAli->GetCurrentFile();
+ fileAli->cd();
+ } // end if add
+
+
+}
+//______________________________________________________________________
+void AliITS::FillModules(TTree *treeH, Int_t mask) {
+ // fill the modules with the sorted by module hits;
+ // can be called many times to do a merging
+ // Inputs:
+ // TTree *treeH The tree containing the hits to be copied into
+ // the modules.
+ // Int_t mask The track number mask to indecate which file
+ // this hits came from.
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+
+ if (treeH == 0x0)
+ {
+ Error("FillModules","Tree is NULL");
+ }
+ Int_t lay,lad,det,index;
+ AliITShit *itsHit=0;
+ AliITSmodule *mod=0;
+ char branchname[20];
+ sprintf(branchname,"%s",GetName());
+ TBranch *branch = treeH->GetBranch(branchname);
+ if (!branch) {
+ Error("FillModules","%s branch in TreeH not found",branchname);
+ return;
+ } // end if !branch
+ branch->SetAddress(&fHits);
+ Int_t nTracks =(Int_t) treeH->GetEntries();
+ Int_t iPrimTrack,h;
+ for(iPrimTrack=0; iPrimTrack<nTracks; iPrimTrack++){
+ ResetHits();
+ Int_t nBytes = treeH->GetEvent(iPrimTrack);
+ if (nBytes <= 0) continue;
+ Int_t nHits = fHits->GetEntriesFast();
+ for(h=0; h<nHits; h++){
+ itsHit = (AliITShit *)fHits->UncheckedAt(h);
+ itsHit->GetDetectorID(lay,lad,det);
+ if (GetITSgeom()) {
+ index = GetITSgeom()->GetModuleIndex(lay,lad,det);
+ } else {
+ index=det-1; // This should not be used.
+ } // end if [You must have fITSgeom for this to work!]
+ mod = GetModule(index);
+ itsHit->SetTrack(itsHit->GetTrack()+mask); // Set track mask.
+ mod->AddHit(itsHit,iPrimTrack,h);
+ } // end loop over hits
+ } // end loop over tracks
+}
+
+//______________________________________________________________________
+Bool_t AliITS::InitModules(Int_t size,Int_t &nmodules){
+ // Initialize the modules array.
+ // Inputs:
+ // Int_t size Size of array of the number of modules to be
+ // created. If size <=0 then the number of modules
+ // is gotten from AliITSgeom class kept in fITSgeom.
+ // Outputs:
+ // Int_t &nmodules The number of modules existing.
+ // Return:
+ // none.
+
+ if(fITSmodules){
+ fITSmodules->Delete();
+ delete fITSmodules;
+ } // end fir fITSmoudles
+
+ if(!fDetTypeSim) {
+ Error("InitModules","fDetTypeSim is null!");
+ return kFALSE;
+ }
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
+
+ Int_t nl,indexMAX,index;
+
+ if(size<=0){ // default to using data stored in AliITSgeom
+ if(fDetTypeSim->GetITSgeom()==0) {
+ Error("InitModules","fITSgeom not defined");
+ return kFALSE;
+ } // end if fITSgeom==0
+ nl = fDetTypeSim->GetITSgeom()->GetNlayers();
+ indexMAX = fDetTypeSim->GetITSgeom()->GetIndexMax();
+ nmodules = indexMAX;
+ fITSmodules = new TObjArray(indexMAX);
+ for(index=0;index<indexMAX;index++){
+ fITSmodules->AddAt( new AliITSmodule(index),index);
+ } // end for index
+ }else{
+ fITSmodules = new TObjArray(size);
+ for(index=0;index<size;index++) {
+ fITSmodules->AddAt( new AliITSmodule(index),index);
+ } // end for index
+
+ nmodules = size;
+ } // end i size<=0
+ return kTRUE;
+}
+//______________________________________________________________________
+void AliITS::Hits2SDigits(){
+ // Standard Hits to summable Digits function.
+ // Inputs:
+ // none.
+ // Outputs:
+ // none.
+
+
+ if(!fDetTypeSim) {
+ Error("Hits2SDigits","fDetTypeSim is null!");
+ return;
+ }
+
+ SetDefaults();
+ fLoader->LoadHits("read");
+ fLoader->LoadSDigits("recreate");
+ AliRunLoader* rl = fLoader->GetRunLoader();
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
+
+ for (Int_t iEvent = 0; iEvent < rl->GetNumberOfEvents(); iEvent++) {
+ // Do the Hits to Digits operation. Use Standard input values.
+ // Event number from file, no background hit merging , use size from
+ // AliITSgeom class, option="All", input from this file only.
+ rl->GetEvent(iEvent);
+ if (!fLoader->TreeS()) fLoader->MakeTree("S");
+ MakeBranch("S");
+ SetTreeAddress();
+ HitsToPreDigits(iEvent,0,-1," ",fOpt," ");
+ } // end for iEvent
+
+ fLoader->UnloadHits();
+ fLoader->UnloadSDigits();
+
+}
+//______________________________________________________________________
+void AliITS::Hits2Digits(){
+
+ //Conversion from hits to digits
+ if(!fDetTypeSim) {
+ Error("Hits2SDigits","fDetTypeSim is 0!");
+ return;
+ }
+
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
+ SetDefaults();
+
+ fLoader->LoadHits("read");
+ fLoader->LoadDigits("recreate");
+ AliRunLoader* rl = fLoader->GetRunLoader();
+ for (Int_t iEvent = 0; iEvent < rl->GetNumberOfEvents(); iEvent++) {
+ rl->GetEvent(iEvent);
+ if (!fLoader->TreeD()) fLoader->MakeTree("D");
+ MakeBranch("D");
+ SetTreeAddress();
+ HitsToDigits(iEvent,0,-1," ",fOpt," ");
+ }
- Int_t ISXFLD = gAlice->Field()->Integ();
- Float_t SXMGMX = gAlice->Field()->Max();
+ fLoader->UnloadHits();
+ fLoader->UnloadDigits();
+}
+
+//______________________________________________________________________
+void AliITS::HitsToDigits(Int_t evNumber,Int_t bgrev,Int_t size,
+ Option_t *option,Option_t *opt,
+ const char *filename){
+ // Keep galice.root for signal and name differently the file for
+ // background when add! otherwise the track info for signal will be lost !
+ // the condition below will disappear when the geom class will be
+ // initialized for all versions - for the moment it is only for v5 !
+ // 7 is the SDD beam test version.
+ // Inputs:
+ // Int_t evnt Event to be processed.
+ // Int_t bgrev Background Hit tree number.
+ // Int_t nmodules Not used.
+ // Option_t *option String indicating if merging hits or not. To
+ // merge hits set equal to "Add". Otherwise no
+ // background hits are considered.
+ // Test_t *filename File name containing the background hits..
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+
+ if(!fDetTypeSim) {
+ Error("HitsToDigits","fDetTypeSim is null!");
+ return;
+ }
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
+ if(!GetITSgeom()) return; // need transformations to do digitization.
+ AliITSgeom *geom = GetITSgeom();
+
+ const char *all = strstr(opt,"All");
+ const char *det[3] = {strstr(opt,"SPD"),strstr(opt,"SDD"),
+ strstr(opt,"SSD")};
+ static Bool_t setDef=kTRUE;
+ if (setDef) SetDefaultSimulation();
+ setDef=kFALSE;
- // --- Define the various materials for GEANT ---
+ Int_t nmodules;
+ InitModules(size,nmodules);
+ FillModules(evNumber,bgrev,nmodules,option,filename);
+
+ // Reset Fast-OR signals for this event
+ fDetTypeSim->ResetFOSignals();
+
+ AliITSsimulation *sim = 0;
+ AliITSmodule *mod = 0;
+ Int_t id;
+ for(Int_t module=0;module<geom->GetIndexMax();module++){
+ id = geom->GetModuleType(module);
+ if (!all && !det[id]) continue;
+ sim = (AliITSsimulation*)fDetTypeSim->GetSimulationModel(id);
+ if (!sim) {
+ Error("HitsToDigits","The simulation class was not "
+ "instanciated for module %d type %s!",module,
+ geom->GetModuleTypeName(module));
+ exit(1);
+ } // end if !sim
+ mod = (AliITSmodule *)fITSmodules->At(module);
+ sim->DigitiseModule(mod,module,evNumber);
+ // fills all branches - wasted disk space
+ fLoader->TreeD()->Fill();
+ ResetDigits();
+ } // end for module
- // 200-224 --> Silicon Pixel Detectors (detectors, chips, buses, cooling,..)
+ ClearModules();
+
+ // Add random noise to FO signals
+ if (all || det[0]) { // SPD present
+ fDetTypeSim->ProcessNoiseForFastOr();
+ }
+
+ // Add Fast-OR signals to event (only one object per event)
+ if (all || det[0]) { // SPD present
+ fDetTypeSim->WriteFOSignals();
+ }
+
- AliMaterial(0, "SPD Si$", 28.0855, 14., 2.33, 9.36, 999);
- AliMaterial(1, "SPD Si chip$", 28.0855, 14., 2.33, 9.36, 999);
- AliMaterial(2, "SPD Si bus$", 28.0855, 14., 2.33, 9.36, 999);
- AliMaterial(3, "SPD C$", 12.011, 6., 2.265,18.8, 999);
- // v. dens
- AliMaterial(4, "SPD Air$", 14.61, 7.3, .001205, 30423., 999);
- AliMaterial(5, "SPD Vacuum$", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
- AliMaterial(6, "SPD Al$", 26.981539, 13., 2.6989, 8.9, 999);
- AliMixture( 7, "SPD Water $", awat, zwat, denswat, -2, wwat);
- AliMixture( 8, "SPD Freon$", afre, zfre, densfre, -2, wfre);
- // **
- AliMedium(0, "SPD Si$", 0, 1,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(1, "SPD Si chip$", 1, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(2, "SPD Si bus$", 2, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(3, "SPD C$", 3, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(4, "SPD Air$", 4, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(5, "SPD Vacuum$", 5, 0,ISXFLD,SXMGMX, 10.,1.00, .1, .100,10.00);
- AliMedium(6, "SPD Al$", 6, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(7, "SPD Water $", 7, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(8, "SPD Freon$", 8, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
+ fLoader->TreeD()->GetEntries();
+ fLoader->TreeD()->AutoSave();
+ // reset tree
+ fLoader->TreeD()->Reset();
+}
+//_____________________________________________________________________
+void AliITS::Hits2PreDigits(){
+ // Turn hits into SDigits
+
+ if(!fDetTypeSim) {
+ Error("Hits2SDigits","fDetTypeSim is 0!");
+ return;
+ }
+
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
+ SetDefaults();
- // 225-249 --> Silicon Drift Detectors (detectors, chips, buses, cooling,..)
+ HitsToPreDigits(fLoader->GetRunLoader()->GetEventNumber(),
+ 0,-1," ",fOpt," ");
+}
+
+//______________________________________________________________________
+void AliITS::HitsToPreDigits(Int_t evNumber,Int_t bgrev,Int_t size,
+ Option_t *option,Option_t *opt,
+ const char *filename){
+ // Keep galice.root for signal and name differently the file for
+ // background when add! otherwise the track info for signal will be lost !
+ // the condition below will disappear when the geom class will be
+ // initialized for all versions - for the moment it is only for v5 !
+ // 7 is the SDD beam test version.
+ // Inputs:
+ // Int_t evnt Event to be processed.
+ // Int_t bgrev Background Hit tree number.
+ // Int_t nmodules Not used.
+ // Option_t *option String indicating if merging hits or not. To
+ // merge hits set equal to "Add". Otherwise no
+ // background hits are considered.
+ // Test_t *filename File name containing the background hits..
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+
+
+ if(!fDetTypeSim) {
+ Error("HitsToPreDigits","fDetTypeSim is null!");
+ return;
+ }
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
+
+ if(!GetITSgeom()){
+ Error("HitsToPreDigits","fGeom is null!");
+ return; // need transformations to do digitization.
+ }
+ AliITSgeom *geom = GetITSgeom();
+
+ const char *all = strstr(opt,"All");
+ const char *det[3] = {strstr(opt,"SPD"),strstr(opt,"SDD"),
+ strstr(opt,"SSD")};
+ static Bool_t setDef=kTRUE;
+ if (setDef) SetDefaultSimulation();
+ setDef=kFALSE;
- AliMaterial(25, "SDD Si$", 28.0855, 14., 2.33, 9.36, 999);
- AliMaterial(26, "SDD Si chip$", 28.0855, 14., 2.33, 9.36, 999);
- AliMaterial(27, "SDD Si bus$", 28.0855, 14., 2.33, 9.36, 999);
- AliMaterial(28, "SDD C$", 12.011, 6., 2.265,18.8, 999);
- // v. dens
- AliMaterial(29, "SDD Air$", 14.61, 7.3, .001205, 30423., 999);
- AliMaterial(30, "SDD Vacuum$", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
- AliMaterial(31, "SDD Al$", 26.981539, 13., 2.6989, 8.9, 999);
- // After a call with ratios by number (negative number of elements),
- // the ratio array is changed to the ratio by weight, so all successive
- // calls with the same array must specify the number of elements as
- // positive
- AliMixture(32, "SDD Water $", awat, zwat, denswat, 2, wwat);
- // After a call with ratios by number (negative number of elements),
- // the ratio array is changed to the ratio by weight, so all successive
- // calls with the same array must specify the number of elements as
- // positive
- AliMixture( 33, "SDD Freon$", afre, zfre, densfre, 2, wfre);
- AliMixture( 34, "SDD PCB$", apcb, zpcb, denspcb, 3, wpcb);
- AliMaterial(35, "SDD Copper$", 63.546, 29., 8.96, 1.43, 999);
- AliMixture( 36, "SDD Ceramics$", acer, zcer, denscer, -5, wcer);
- AliMaterial(37, "SDD Kapton$", 12.011, 6., 1.3, 31.27, 999);
- // **
- // check A and Z
- AliMedium(25, "SDD Si$", 25, 1,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(26, "SDD Si chip$", 26, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(27, "SDD Si bus$", 27, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(28, "SDD C$", 28, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(29, "SDD Air$", 29, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(30, "SDD Vacuum$", 30, 0,ISXFLD,SXMGMX, 10.,1.00, .1, .100,10.00);
- AliMedium(31, "SDD Al$", 31, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(32, "SDD Water $", 32, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(33, "SDD Freon$", 33, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(34, "SDD PCB$", 34, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(35, "SDD Copper$", 35, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(36, "SDD Ceramics$",36, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(37, "SDD Kapton$", 37, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
+ Int_t nmodules;
+ InitModules(size,nmodules);
+ FillModules(evNumber,bgrev,nmodules,option,filename);
- // 250-274 --> Silicon Strip Detectors (detectors, chips, buses, cooling,..)
+
+ AliITSsimulation *sim = 0;
+ AliITSmodule *mod = 0;
+ Int_t id,module;
+ for(module=0;module<geom->GetIndexMax();module++){
+ id = geom->GetModuleType(module);
+ if (!all && !det[id]) continue;
+ sim = (AliITSsimulation*)GetSimulationModel(id);
+ if (!sim) {
+ Error("HitsToPreDigits","The simulation class was not "
+ "instanciated for module %d type %s!",module,
+ geom->GetModuleTypeName(module));
+ exit(1);
+ } // end if !sim
+ mod = (AliITSmodule *)fITSmodules->At(module);
+ sim->SDigitiseModule(mod,module,evNumber);
+ // fills all branches - wasted disk space
+ fLoader->TreeS()->Fill();
+ fDetTypeSim->ResetSDigits();
+ } // end for module
+
+ ClearModules();
+
- AliMaterial(50, "SSD Si$", 28.0855, 14., 2.33, 9.36, 999.);
- AliMaterial(51, "SSD Si chip$", 28.0855, 14., 2.33, 9.36, 999.);
- AliMaterial(52, "SSD Si bus$", 28.0855, 14., 2.33, 9.36, 999.);
- AliMaterial(53, "SSD C$", 12.011, 6., 2.265,18.8, 999.);
- // v. dens
- AliMaterial(54, "SSD Air$", 14.61, 7.3, .001205, 30423., 999);
- AliMaterial(55, "SSD Vacuum$", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
- AliMaterial(56, "SSD Al$", 26.981539, 13., 2.6989, 8.9, 999);
- // After a call with ratios by number (negative number of elements),
- // the ratio array is changed to the ratio by weight, so all successive
- // calls with the same array must specify the number of elements as
- // positive
- AliMixture(57, "SSD Water $", awat, zwat, denswat, 2, wwat);
- // After a call with ratios by number (negative number of elements),
- // the ratio array is changed to the ratio by weight, so all successive
- // calls with the same array must specify the number of elements as
- // positive
- AliMixture(58, "SSD Freon$", afre, zfre, densfre, 2, wfre);
- AliMixture(59, "SSD PCB$", apcb, zpcb, denspcb, 3, wpcb);
- AliMaterial(60, "SSD Copper$", 63.546, 29., 8.96, 1.43, 999.);
- // After a call with ratios by number (negative number of elements),
- // the ratio array is changed to the ratio by weight, so all successive
- // calls with the same array must specify the number of elements as
- // positive
- AliMixture( 61, "SSD Ceramics$", acer, zcer, denscer, 5, wcer);
- AliMaterial(62, "SSD Kapton$", 12.011, 6., 1.3, 31.27, 999.);
- // check A and Z
- AliMaterial(63, "SDD G10FR4$", 17.749, 8.875, 1.8, 21.822, 999.);
- // **
- AliMedium(50, "SSD Si$", 50, 1,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(51, "SSD Si chip$", 51, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(52, "SSD Si bus$", 52, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(53, "SSD C$", 53, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(54, "SSD Air$", 54, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(55, "SSD Vacuum$", 55, 0,ISXFLD,SXMGMX, 10.,1.00, .1, .100,10.00);
- AliMedium(56, "SSD Al$", 56, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(57, "SSD Water $", 57, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(58, "SSD Freon$", 58, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(59, "SSD PCB$", 59, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(60, "SSD Copper$", 60, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(61, "SSD Ceramics$",61, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(62, "SSD Kapton$", 62, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(63, "SSD G10FR4$", 63, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
+ fLoader->TreeS()->GetEntries();
+ fLoader->TreeS()->AutoSave();
+ fLoader->WriteSDigits("OVERWRITE");
+ // reset tree
+ fLoader->TreeS()->Reset();
+}
+
+//_____________________________________________________________________
+void AliITS::HitsToFastRecPoints(Int_t evNumber,Int_t bgrev,Int_t size,
+ Option_t *opt0,Option_t *opt1,
+ const char *flnm){
+ // keep galice.root for signal and name differently the file for
+ // background when add! otherwise the track info for signal will be lost !
+ // the condition below will disappear when the geom class will be
+ // initialized for all versions - for the moment it is only for v5 !
+ // Inputs:
+ // Int_t evnt Event to be processed.
+ // Int_t bgrev Background Hit tree number.
+ // Int_t size Size used by InitModules. See InitModules.
+ // Option_t *opt0 Option passed to FillModules. See FillModules.
+ // Option_t *opt1 String indicating if merging hits or not. To
+ // merge hits set equal to "Add". Otherwise no
+ // background hits are considered.
+ // Test_t *flnm File name containing the background hits..
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+
+
+
+ if(!GetITSgeom()){
+ Error("HitsToPreDigits","fGeom is null!");
+ return; // need transformations to do digitization.
+ }
+ AliITSgeom *geom = GetITSgeom();
+
+ AliITSLoader *pITSloader = (AliITSLoader*)fLoader;
+
+ const char *all = strstr(opt1,"All");
+ const char *det[3] ={strstr(opt1,"SPD"),strstr(opt1,"SDD"),
+ strstr(opt1,"SSD")};
+ Int_t nmodules;
+ InitModules(size,nmodules);
+ FillModules(evNumber,bgrev,nmodules,opt0,flnm);
+
+ AliITSsimulation *sim = 0;
+ AliITSmodule *mod = 0;
+ Int_t id,module;
+
+ TTree *lTR = pITSloader->TreeR();
+ if(!lTR) {
+ pITSloader->MakeTree("R");
+ lTR = pITSloader->TreeR();
+ }
- // 275-299 --> General (end-caps, frames, cooling, cables, etc.)
+ TClonesArray* ptarray = new TClonesArray("AliITSRecPoint",1000);
+ TBranch* branch = (TBranch*)lTR->Branch("ITSRecPointsF",&ptarray);
+ branch->SetAddress(&ptarray);
+ //m.b. : this change is nothing but a nice way to make sure
+ //the CPU goes up !
+ for(module=0;module<geom->GetIndexMax();module++){
+ id = geom->GetModuleType(module);
+ if (!all && !det[id]) continue;
+ sim = (AliITSsimulation*)GetSimulationModel(id);
+ if (!sim) {
+ Error("HitsToFastPoints","The simulation class was not "
+ "instantiated for module %d type %s!",module,
+ geom->GetModuleTypeName(module));
+ exit(1);
+ } // end if !sim
+ mod = (AliITSmodule *)fITSmodules->At(module);
+ sim->CreateFastRecPoints(mod,module,gRandom,ptarray);
+ lTR->Fill();
+ ptarray->Clear();
+ } // end for module
+
+ ClearModules();
+ fLoader->WriteRecPoints("OVERWRITE");
+ delete ptarray;
+}
+//_____________________________________________________________________
+Int_t AliITS::Hits2Clusters(TTree *hTree, TTree *cTree) {
+ //------------------------------------------------------------
+ // This function creates ITS clusters
+ //------------------------------------------------------------
+ if(!GetITSgeom()){
+ Error("HitsToPreDigits","fGeom is null!");
+ return 1; // need transformations to do digitization.
+ }
+ AliITSgeom *geom=GetITSgeom();
+ Int_t mmax=geom->GetIndexMax();
+
+ InitModules(-1,mmax);
+ FillModules(hTree,0);
+
+ TClonesArray *points = new TClonesArray("AliITSRecPoint",1000);
+ TBranch *branch=cTree->GetBranch("ITSRecPoints");
+ if (!branch) cTree->Branch("ITSRecPoints",&points);
+ else branch->SetAddress(&points);
+
+ AliITSsimulationFastPoints sim;
+ Int_t ncl=0;
+ for (Int_t m=0; m<mmax; m++) {
+ AliITSmodule *mod=GetModule(m);
+ sim.CreateFastRecPoints(mod,m,gRandom,points);
+ ncl+=points->GetEntriesFast();
+ cTree->Fill();
+ points->Clear();
+ }
+
+ AliDebug(1,Form("Number of found fast clusters : %d",ncl));
+
+ //cTree->Write();
+
+ delete points;
+ return 0;
+}
+
+//_____________________________________________________________________
+void AliITS::CheckLabels(Int_t lab[3]) const {
+ //------------------------------------------------------------
+ // Tries to find mother's labels
+ //------------------------------------------------------------
+
+ if(lab[0]<0 && lab[1]<0 && lab[2]<0) return; // In case of no labels just exit
+
+ Int_t ntracks = gAlice->GetMCApp()->GetNtrack();
+ for (Int_t i=0;i<3;i++){
+ Int_t label = lab[i];
+ if (label>=0 && label<ntracks) {
+ TParticle *part=(TParticle*)gAlice->GetMCApp()->Particle(label);
+ if (part->P() < 0.005) {
+ Int_t m=part->GetFirstMother();
+ if (m<0) {
+ continue;
+ }
+ if (part->GetStatusCode()>0) {
+ continue;
+ }
+ lab[i]=m;
+ }
+ }
+ }
- AliMaterial(75, "GEN C$", 12.011, 6., 2.265, 18.8, 999.);
- // verify density
- AliMaterial(76, "GEN Air$", 14.61, 7.3, .001205, 30423., 999);
- AliMaterial(77, "GEN Vacuum$", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
- AliMixture( 78, "GEN POLYETHYL$", apoly, zpoly, .95, -2, wpoly);
- AliMixture( 79, "GEN SERVICES$", aserv, zserv, 4.68, 4, wserv);
- AliMaterial(80, "GEN Copper$", 63.546, 29., 8.96, 1.43, 999.);
- // After a call with ratios by number (negative number of elements),
- // the ratio array is changed to the ratio by weight, so all successive
- // calls with the same array must specify the number of elements as
- // positive
- AliMixture(81, "GEN Water $", awat, zwat, denswat, 2, wwat);
- // **
- AliMedium(75,"GEN C$", 75, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(76,"GEN Air$", 76, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(77,"GEN Vacuum$", 77, 0,ISXFLD,SXMGMX, 10., .10, .1, .100,10.00);
- AliMedium(78,"GEN POLYETHYL$",78, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(79,"GEN SERVICES$", 79, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(80,"GEN Copper$", 80, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
- AliMedium(81,"GEN Water $", 81, 0,ISXFLD,SXMGMX, 10., .01, .1, .003, .003);
-}
-
-//_____________________________________________________________________________
-Int_t AliITS::DistancetoPrimitive(Int_t , Int_t ){
- //
- // Distance from mouse to ITS on the screen. Dummy routine
- // A dummy routine used by the ROOT macro display.C to allow for the
- // use of the mouse (pointing device) in the macro. In general this should
- // never be called. If it is it returns the number 9999 for any value of
- // x and y.
- //
- return 9999;
}
-//_____________________________________________________________________________
-void AliITS::Init(){
+//______________________________________________________________________
+void AliITS::SDigitsToDigits(Option_t *opt){
+ // Standard Summable digits to Digits function.
+ // Inputs:
+ // none.
+ // Outputs:
+ // none.
+ if (!fDetTypeSim) {
+ AliError("fDetTypeSim is 0!");
+ return;
+ }
+
+ const char *all = strstr(opt,"All");
+ const char *det[3] ={strstr(opt,"SPD"),strstr(opt,"SDD"),
+ strstr(opt,"SSD")};
+
+ // Reset Fast-OR signals for this event
+ fDetTypeSim->ResetFOSignals();
+
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+ SetDefaults();
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
+ fDetTypeSim->SDigitsToDigits(opt,(Char_t*)GetName());
+
+ // Add random noise to FO signals
+ if (all || det[0]) { // SPD present
+ fDetTypeSim->ProcessNoiseForFastOr();
+ }
+ // Add Fast-OR signals to event (only one object per event)
+ if (all || det[0]) { // SPD present
+ fDetTypeSim->WriteFOSignals();
+ }
+}
+
+//______________________________________________________________________
+void AliITS::ResetDigits(){
+ // Reset number of digits and the digits array for the ITS detector.
+ // Inputs:
+ // none.
+ // Outputs:
+ // none.
+ if(!fDetTypeSim) {
+ Error("ResetDigits","fDetTypeSim is 0!");
+ return;
+ }
+
+ fDetTypeSim->ResetDigits();
+
+
+}
+//______________________________________________________________________
+void AliITS::ResetDigits(Int_t branch){
+ // Reset number of digits and the digits array for this branch.
+ // Inputs:
+ // none.
+ // Outputs:
+ // none.
+
+ if(!fDetTypeSim) {
+ Error("ResetDigits","fDetTypeSim is 0!");
+ return;
+ }
+
+ fDetTypeSim->ResetDigits(branch);
+
+}
+//______________________________________________________________________
+void AliITS::AddSumDigit(AliITSpListItem &sdig){
+ // Adds the a module full of summable digits to the summable digits tree.
+ // Inputs:
+ // AliITSpListItem &sdig SDigit to be added to SDigits tree.
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+
+ if(!fDetTypeSim) {
+ Error("AddSumDigit","fDetTypeSim is 0!");
+ return;
+ }
+ fDetTypeSim->AddSumDigit(sdig);
+
+}
+//______________________________________________________________________
+void AliITS::AddSimDigit(Int_t branch, AliITSdigit *d){
+ // Add a simulated digit.
+ // Inputs:
+ // Int_t id Detector type number.
+ // AliITSdigit *d Digit to be added to the Digits Tree. See
+ // AliITSdigit.h
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+
+ if(!fDetTypeSim) {
+ Error("AddSimDigit","fDetTypeSim is 0!");
+ return;
+ }
+ fDetTypeSim->AddSimDigit(branch,d);
+
+}
+//______________________________________________________________________
+void AliITS::AddSimDigit(Int_t branch,Float_t phys,Int_t *digits,Int_t *tracks,
+ Int_t *hits,Float_t *charges, Int_t sigexpanded){
+ // Add a simulated digit to the list.
+ // Inputs:
+ // Int_t id Detector type number.
+ // Float_t phys Physics indicator. See AliITSdigits.h
+ // Int_t *digits Integer array containing the digits info. See
+ // AliITSdigit.h
+ // Int_t *tracks Integer array [AliITSdigitS?D::GetNTracks()]
+ // containing the track numbers that contributed to
+ // this digit.
+ // Int_t *hits Integer array [AliITSdigitS?D::GetNTracks()]
+ // containing the hit numbers, from AliITSmodule, that
+ // contributed to this digit.
+ // Float_t *charge Floating point array of the signals contributed
+ // to this digit by each track.
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+
+ if(!fDetTypeSim) {
+ Error("AddSimDigit","fDetTypeSim is 0!");
+ return;
+ }
+ fDetTypeSim->AddSimDigit(branch,phys,digits,tracks,hits,charges,sigexpanded);
+
+}
+//______________________________________________________________________
+void AliITS::Digits2Raw(){
+ // convert digits of the current event to raw data
+
+ if(!fDetTypeSim) {
+ Error("Digits2Raw","fDetTypeSim is 0!");
+ return;
+ }
+ fDetTypeSim->SetLoader((AliITSLoader*)fLoader);
+ SetDefaults();
+ if(fSimuParam) fDetTypeSim->SetSimuParam(fSimuParam);
+ fDetTypeSim->GetLoader()->LoadDigits();
+ TTree* digits = fDetTypeSim->GetLoader()->TreeD();
+ if (!digits) {
+ Error("Digits2Raw", "no digits tree");
+ return;
+ }
+ fDetTypeSim->SetTreeAddressD(digits,(Char_t*)GetName());
+
+ // Get the FO signals for this event
+ AliITSFOSignalsSPD* foSignals = NULL;
+ AliRunLoader* runLoader = AliRunLoader::Instance();
+ AliITSLoader* itsLoader = (AliITSLoader*) runLoader->GetLoader("ITSLoader");
+ if (!itsLoader) {
+ AliError("ITS loader is NULL.");
+ }
+ else {
+ if(!itsLoader->TreeD()) AliError(" !!! No TreeD available !!!");
+ foSignals = (AliITSFOSignalsSPD*)itsLoader->TreeD()->GetUserInfo()->FindObject("AliITSFOSignalsSPD");
+ if(!foSignals) AliError("FO signals not retrieved");
+ }
+
+ Bool_t deleteFOsignalsLater = kFALSE;
+ if (!foSignals) {
+ AliError("FO signals not available. No FO bits will be written.");
+ foSignals = new AliITSFOSignalsSPD(); // make a temporary dummy signals object
+ deleteFOsignalsLater = kTRUE;
+ }
+
+
+ AliITSDDLModuleMapSDD* ddlsdd=fDetTypeSim->GetDDLModuleMapSDD();
+ Char_t rawSDD=fDetTypeSim->GetSimuParam()->GetSDDRawDataFormat();
+ AliITSDDLRawData rawWriter;
+
+ rawWriter.SetSDDRawFormat(rawSDD);
+ //Verbose level
+ // 0: Silent
+ // 1: cout messages
+ // 2: txt files with digits
+ //BE CAREFUL, verbose level 2 MUST be used only for debugging and
+ //it is highly suggested to use this mode only for debugging digits files
+ //reasonably small, because otherwise the size of the txt files can reach
+ //quickly several MB wasting time and disk space.
+ rawWriter.SetVerbose(0);
+
+ //SILICON PIXEL DETECTOR
+ AliDebug(1,"Formatting raw data for SPD");
+ rawWriter.RawDataSPD(digits->GetBranch("ITSDigitsSPD"),foSignals);
+ if(deleteFOsignalsLater) delete foSignals;
+
+ //SILICON DRIFT DETECTOR
+ AliDebug(1,Form("Formatting raw data for SDD - Format code =%d",rawSDD));
+ rawWriter.RawDataSDD(digits->GetBranch("ITSDigitsSDD"),ddlsdd);
+
+ //SILICON STRIP DETECTOR
+ AliDebug(1,"Formatting raw data for SSD");
+ rawWriter.RawDataSSD(digits->GetBranch("ITSDigitsSSD"));
+
+ fLoader->UnloadDigits();
+}
+//______________________________________________________________________
+AliLoader* AliITS::MakeLoader(const char* topfoldername){
+ //builds ITSgetter (AliLoader type)
+ //if detector wants to use castomized getter, it must overload this method
+
+ AliDebug(1,Form("Creating AliITSLoader. Top folder is %s.",
+ topfoldername));
+ fLoader = new AliITSLoader(GetName(),topfoldername);
+ return fLoader;
+}
+//______________________________________________________________________
+Bool_t AliITS::Raw2SDigits(AliRawReader* rawReader)
+{
//
- // Initialise ITS after it has been built
- // This routine initializes the AliITS class. It is intended to be called
- // from the Init function in AliITSv?. Besides displaying a banner
- // indicating that it has been called it initializes the array fIdSens.
- // Therefore it should be called after a call to CreateGeometry.
+ // Converts RAW data to SDigits
//
- Int_t i;
+ // Get TreeS
//
- printf("\n");
- for(i=0;i<35;i++) printf("*");
- printf(" ITS_INIT ");
- for(i=0;i<35;i++) printf("*");
- printf("\n");
+ Int_t last = -1;
+ Int_t size = AliITSgeomTGeo::GetNModules();
+ if(!fModA) {
+ fModA = new TClonesArray*[size];
+ for (Int_t mod = 0; mod < size; mod++) fModA[mod] = new TClonesArray("AliITSpListItem", 10000);
+ }
+ AliLoader* loader = (AliRunLoader::Instance())->GetLoader("ITSLoader");
+ if (!loader){
+ Error("Open","Can not get ITS loader from Run Loader");
+ return kFALSE;
+ }
+
+ TTree* tree = 0;
+ tree = loader->TreeS();
+ if (!tree){
+ loader->MakeTree("S");
+ tree = loader->TreeS();
+ }
//
+ // Array for SDigits
+ //
+ if(!fpSDigits){
+ fpSDigits = new TClonesArray("AliITSpListItem",10000);
+ }
+ TClonesArray& aSDigits = *fpSDigits;
+ Int_t bufsize = 32000;
+ tree->Branch("ITS", &fpSDigits, bufsize);
+ Int_t npx = 0;
//
- for(i=0;i<fIdN;i++) fIdSens[i] = gMC->VolId(fIdName[i]);
+ // SPD
//
- for(i=0;i<80;i++) printf("*");
- printf("\n");
-}
+ AliITSsegmentationSPD* segSPD = (AliITSsegmentationSPD*) fDetTypeSim->GetSegmentationModel(0);
+ if(!segSPD){
+ AliWarning("Set AliITS defaults");
+ SetDefaults();
+ segSPD = (AliITSsegmentationSPD*) fDetTypeSim->GetSegmentationModel(0);
+ }
+ npx = segSPD->Npx();
+ Double_t thr, sigma;
+
+ AliITSRawStreamSPD inputSPD(rawReader);
+ while(1){
+ Bool_t next = inputSPD.Next();
+ if (!next) break;
+
+ Int_t module = inputSPD.GetModuleID();
+ Int_t column = inputSPD.GetColumn();
+ Int_t row = inputSPD.GetRow();
+ Int_t index = npx * column + row;
+
+ if (module >= size) continue;
+
+ last = (fModA[module])->GetEntries();
+ TClonesArray& dum = *fModA[module];
+ fDetTypeSim->GetSimuParam()->SPDThresholds(module,thr,sigma);
+ thr += 1.;
+ new (dum[last]) AliITSpListItem(-1, -1, module, index, thr);
+ }
+ rawReader->Reset();
-//_____________________________________________________________________________
-void AliITS::MakeBranch(Option_t* option){
//
- // Create Tree branches for the ITS.
- // Creates the TTree branch where the class AliITS is kept.
+ // SDD
+ //
+ AliITSsegmentationSDD* segSDD = (AliITSsegmentationSDD*) fDetTypeSim->GetSegmentationModel(1);
+ npx = segSDD->Npx();
+ Int_t scalef=AliITSsimulationSDD::ScaleFourier(segSDD);
+ Int_t firstSDD=AliITSgeomTGeo::GetModuleIndex(3,1,1);
+ Int_t firstSSD=AliITSgeomTGeo::GetModuleIndex(5,1,1);
+
+ AliITSRawStream* inputSDD=AliITSRawStreamSDD::CreateRawStreamSDD(rawReader);
+ for(Int_t iMod=firstSDD; iMod<firstSSD; iMod++){
+ AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)fDetTypeSim->GetCalibrationModel(iMod);
+ Bool_t isZeroSupp=cal->GetZeroSupp();
+ if(isZeroSupp){
+ for(Int_t iSid=0; iSid<2; iSid++) inputSDD->SetZeroSuppLowThreshold(iMod-firstSDD,iSid,cal->GetZSLowThreshold(iSid));
+ }else{
+ for(Int_t iSid=0; iSid<2; iSid++) inputSDD->SetZeroSuppLowThreshold(iMod-firstSDD,iSid,0);
+ }
+ }
+
+ AliITSDDLModuleMapSDD* ddlmap=fDetTypeSim->GetDDLModuleMapSDD();
+ inputSDD->SetDDLModuleMap(ddlmap);
+ while(inputSDD->Next()){
+ if(inputSDD->IsCompletedModule()==kFALSE &&
+ inputSDD->IsCompletedDDL()==kFALSE){
+
+ Int_t module = inputSDD->GetModuleID();
+ Int_t anode = inputSDD->GetCoord1()+segSDD->NpzHalf()*inputSDD->GetChannel();
+ Int_t time = inputSDD->GetCoord2();
+ Int_t signal10 = inputSDD->GetSignal();
+ Int_t index = AliITSpList::GetIndex(anode,time,scalef*npx);
+
+ if (module >= size) continue;
+ last = fModA[module]->GetEntries();
+ TClonesArray& dum = *fModA[module];
+ new (dum[last]) AliITSpListItem(-1, -1, module, index, Double_t(signal10));
+ ((AliITSpListItem*) dum.At(last))->AddSignalAfterElect(module, index, Double_t(signal10));
+ }
+ }
+ delete inputSDD;
+ rawReader->Reset();
+
//
- Int_t buffersize = 4000;
- char branchname[10];
- sprintf(branchname,"%s",GetName());
+ // SSD
+ //
+ AliITSsegmentationSSD* segSSD = (AliITSsegmentationSSD*) fDetTypeSim->GetSegmentationModel(2);
+ npx = segSSD->Npx();
+ AliITSRawStreamSSD inputSSD(rawReader);
+ while(1){
+ Bool_t next = inputSSD.Next();
+ if (!next) break;
- AliDetector::MakeBranch(option);
+ Int_t module = inputSSD.GetModuleID();
+ if(module<0)AliError(Form("Invalid SSD module %d \n",module));
+ if(module<0)continue;
+ Int_t side = inputSSD.GetSideFlag();
+ Int_t strip = inputSSD.GetStrip();
+ Int_t signal = inputSSD.GetSignal();
+ Int_t index = npx * side + strip;
+
+ if (module >= size) continue;
+
+ last = fModA[module]->GetEntries();
+ TClonesArray& dum = *fModA[module];
+ new (dum[last]) AliITSpListItem(-1, -1, module, index, Double_t(signal));
+ }
+ rawReader->Reset();
+ AliITSpListItem* sdig = 0;
+
+ Int_t firstssd = GetITSgeom()->GetStartDet(kSSD);
+ Double_t adcToEv = 1.;
+ for (Int_t mod = 0; mod < size; mod++)
+ {
+ if(mod>=firstssd) {
+ AliITSCalibrationSSD* calssd = (AliITSCalibrationSSD*)fDetTypeSim->GetCalibrationModel(mod);
+ adcToEv = 1./calssd->GetSSDDEvToADC(1.);
+ }
+ Int_t nsdig = fModA[mod]->GetEntries();
+ for (Int_t ie = 0; ie < nsdig; ie++) {
+ sdig = (AliITSpListItem*) (fModA[mod]->At(ie));
+ Double_t digsig = sdig->GetSignal();
+ if(mod>=firstssd) digsig*=adcToEv; // for SSD: convert back charge from ADC to electron
+ new (aSDigits[ie]) AliITSpListItem(-1, -1, mod, sdig->GetIndex(), digsig);
+ Float_t sig = sdig->GetSignalAfterElect();
+ if(mod>=firstssd) sig*=adcToEv;
+ if (sig > 0.) {
+ sdig = (AliITSpListItem*)aSDigits[ie];
+ sdig->AddSignalAfterElect(mod, sdig->GetIndex(), Double_t(sig));
+ }
+ }
+
+ tree->Fill();
+ aSDigits.Clear();
+ fModA[mod]->Clear();
+ }
+ loader->WriteSDigits("OVERWRITE");
+ return kTRUE;
+}
- char *D = strstr(option,"D");
-
- if (fDigits && gAlice->TreeD() && D) {
- gAlice->TreeD()->Branch(branchname,&fDigits, buffersize);
- printf("Making Branch %s for digits\n",branchname);
- } // end if
-}
-
-//____________________________________________________________________________
-void AliITS::Streamer(TBuffer &R__b){
- // Stream an object of class AliITS.
- Int_t i,j,l;
-
- if (R__b.IsReading()) {
- Version_t R__v = R__b.ReadVersion();
- if (R__v == 1) {
- AliDetector::Streamer(R__b);
- R__b >> fITSgeom;
-// R__b >> fITSmodules; //We do not write out modules so don't read them
- R__b >> fITSpoints;
- R__b >> fEuclidOut;
- R__b >> fIdN;
- if(fIdSens!=0) delete[] fIdSens;
- if(fIdName!=0) delete[] fIdName;
- fIdSens = new Int_t[fIdN];
- fIdName = new char*[fIdN];
- for(i=0;i<fIdN;i++) R__b >> fIdSens[i];
- for(i=0;i<fIdN;i++){
- R__b >> l;
- fIdName[i] = new char[l+1]; // add room for null character.
- for(j=0;j<l;j++) R__b >> fIdName[i][j];
- fIdName[i][l] = '\0'; // Null terminate this string.
- } // end for i
- R__b >> fMajorVersion;
- R__b >> fMinorVersion;
- } // end if (R__v)
- } else {
- R__b.WriteVersion(AliITS::IsA());
- AliDetector::Streamer(R__b);
- R__b << fITSgeom;
-// R__b << fITSmodules; //We don't want to write out the modules class.
- R__b << fITSpoints;
- R__b << fEuclidOut;
- R__b << fIdN;
- for(i=0;i<fIdN;i++) R__b <<fIdSens[i];
- for(i=0;i<fIdN;i++){
- l = strlen(fIdName[i]);
- R__b << l;
- for(j=0;j<l;j++) R__b << fIdName[i][j];
- } // end for i
- R__b << fMajorVersion;
- R__b << fMinorVersion;
- }
+//______________________________________________________________________
+void AliITS::UpdateInternalGeometry(){
+ //reads new geometry from TGeo
+// AliDebug(1,"Delete ITSgeom and create a new one reading TGeo");
+
+ AliITSVersion_t version = (AliITSVersion_t)IsVersion();
+ Int_t minor = 0;
+ if(version==kvPPRasymmFMD) minor=2; // default minor version for this geom.
+ AliITSInitGeometry initgeom;
+ AliITSgeom* geom = initgeom.CreateAliITSgeom(version,minor);
+ SetITSgeom(geom);
+}
+//______________________________________________________________________
+AliTriggerDetector* AliITS::CreateTriggerDetector() const {
+ // create an AliITSTrigger object (and set trigger conditions as input)
+ return new AliITSTrigger(fDetTypeSim->GetTriggerConditions());
}