Adding the number of fired chips in the ITS fast OR (Domenico, Massimo)
[u/mrichter/AliRoot.git] / ITS / AliITSReconstructor.cxx
CommitLineData
121a60bd 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///////////////////////////////////////////////////////////////////////////////
19// //
20// class for ITS reconstruction //
21// //
22///////////////////////////////////////////////////////////////////////////////
23
e62c1aea 24#include "Riostream.h"
121a60bd 25#include "AliITSReconstructor.h"
00a7cc50 26#include "AliRun.h"
113c12f1 27#include "AliRawReader.h"
00a7cc50 28#include "AliITSDetTypeRec.h"
6cae184e 29#include "AliITSgeom.h"
7d62fb64 30#include "AliITSLoader.h"
e43c066c 31#include "AliITStrackerMI.h"
c0342e67 32#include "AliITStrackerV2.h"
121a60bd 33#include "AliITStrackerSA.h"
34#include "AliITSVertexerIons.h"
35#include "AliITSVertexerFast.h"
c8735dd3 36#include "AliITSVertexerFixed.h"
3acc14d5 37#include "AliITSVertexer3D.h"
121a60bd 38#include "AliITSVertexerZ.h"
3acc14d5 39#include "AliITSVertexerCosmics.h"
af885e0f 40#include "AliESDEvent.h"
121a60bd 41#include "AliITSpidESD.h"
e62c1aea 42#include "AliITSpidESD1.h"
43#include "AliITSpidESD2.h"
023ae34b 44#include "AliITSInitGeometry.h"
121a60bd 45
44347160 46
121a60bd 47ClassImp(AliITSReconstructor)
44347160 48
e62c1aea 49//___________________________________________________________________________
94631b2f 50AliITSReconstructor::AliITSReconstructor() : AliReconstructor(),
6cae184e 51fItsPID(0),
52fDetTypeRec(0)
94631b2f 53{
e62c1aea 54 // Default constructor
e62c1aea 55}
56 //___________________________________________________________________________
57AliITSReconstructor::~AliITSReconstructor(){
58// destructor
59 delete fItsPID;
6cae184e 60 if(fDetTypeRec) delete fDetTypeRec;
e62c1aea 61}
62//______________________________________________________________________
94631b2f 63AliITSReconstructor::AliITSReconstructor(const AliITSReconstructor &ob) :AliReconstructor(ob),
5c47990a 64fItsPID(ob.fItsPID),
65fDetTypeRec(ob.fDetTypeRec)
6cae184e 66
94631b2f 67{
e62c1aea 68 // Copy constructor
e62c1aea 69}
70
71//______________________________________________________________________
94631b2f 72AliITSReconstructor& AliITSReconstructor::operator=(const AliITSReconstructor& ob ){
e62c1aea 73 // Assignment operator
94631b2f 74 this->~AliITSReconstructor();
75 new(this) AliITSReconstructor(ob);
e62c1aea 76 return *this;
77}
6cae184e 78
023ae34b 79//______________________________________________________________________
d76c31f4 80void AliITSReconstructor::Init() {
023ae34b 81 // Initalize this constructor bet getting/creating the objects
82 // nesseary for a proper ITS reconstruction.
83 // Inputs:
d76c31f4 84 // none.
023ae34b 85 // Output:
86 // none.
87 // Return:
88 // none.
00a7cc50 89
108bd0fe 90 AliITSInitGeometry initgeom;
91 AliITSgeom *geom = initgeom.CreateAliITSgeom();
92 AliInfo(Form("Geometry name: %s",(initgeom.GetGeometryName()).Data()));
6cae184e 93
94 fDetTypeRec = new AliITSDetTypeRec();
95 fDetTypeRec->SetITSgeom(geom);
96 fDetTypeRec->SetDefaults();
ed446fa3 97
023ae34b 98 return;
99}
6cae184e 100
121a60bd 101//_____________________________________________________________________________
6cae184e 102void AliITSReconstructor::Reconstruct(TTree *digitsTree, TTree *clustersTree) const
121a60bd 103{
104// reconstruct clusters
105
9f9cae94 106 Int_t cluFindOpt = GetRecoParam()->GetClusterFinder();
107 Bool_t useV2=kTRUE; // Default: V2 cluster finder
108 if(cluFindOpt==1) useV2=kFALSE;
121a60bd 109
6cae184e 110 fDetTypeRec->SetTreeAddressD(digitsTree);
111 fDetTypeRec->MakeBranch(clustersTree,"R");
112 fDetTypeRec->SetTreeAddressR(clustersTree);
9f9cae94 113 fDetTypeRec->DigitsToRecPoints(digitsTree,clustersTree,0,"All",useV2);
121a60bd 114}
115
00a7cc50 116//_________________________________________________________________
6cae184e 117void AliITSReconstructor::Reconstruct(AliRawReader* rawReader, TTree *clustersTree) const
113c12f1 118{
6cae184e 119 // reconstruct clusters from raw data
00a7cc50 120
6cae184e 121 fDetTypeRec->SetDefaultClusterFindersV2(kTRUE);
122 fDetTypeRec->DigitsToRecPoints(rawReader,clustersTree);
113c12f1 123}
124
125//_____________________________________________________________________________
d76c31f4 126AliTracker* AliITSReconstructor::CreateTracker() const
121a60bd 127{
128// create a ITS tracker
129
9f9cae94 130 Int_t trackerOpt = GetRecoParam()->GetTracker();
e62c1aea 131 AliTracker* tracker;
9f9cae94 132 if (trackerOpt==1) {
e341247d 133 tracker = new AliITStrackerMI(0);
23197852 134 AliITStrackerMI *mit=(AliITStrackerMI*)tracker;
135 mit->SetDetTypeRec(fDetTypeRec);
c0342e67 136 }
9f9cae94 137 else if (trackerOpt==2) {
c0342e67 138 tracker = new AliITStrackerV2(0);
e62c1aea 139 }
140 else {
cc088660 141 tracker = new AliITStrackerSA(0); // inherits from AliITStrackerMI
1966f03c 142 AliITStrackerSA *sat=(AliITStrackerSA*)tracker;
23197852 143 sat->SetDetTypeRec(fDetTypeRec);
9f9cae94 144 if(GetRecoParam()->GetTrackerSAOnly()) sat->SetSAFlag(kTRUE);
1966f03c 145 if(sat->GetSAFlag())AliDebug(1,"Tracking Performed in ITS only\n");
9f9cae94 146 sat->SetOuterStartLayer(GetRecoParam()->GetOuterStartLayerSA());
e62c1aea 147 }
7d62fb64 148
9f9cae94 149 Int_t pidOpt = GetRecoParam()->GetPID();
150
d76c31f4 151 AliITSReconstructor* nc = const_cast<AliITSReconstructor*>(this);
9f9cae94 152 if(pidOpt==1){
d76c31f4 153 Info("FillESD","ITS LandauFitPID option has been selected\n");
154 nc->fItsPID = new AliITSpidESD2((AliITStrackerMI*)tracker);
e62c1aea 155 }
156 else{
d76c31f4 157 Info("FillESD","ITS default PID\n");
7a8614f3 158 Double_t parITS[] = {79.,0.13, 5.}; //IB: this is "pp tuning"
d76c31f4 159 nc->fItsPID = new AliITSpidESD1(parITS);
e62c1aea 160 }
38013a1a 161
e62c1aea 162 return tracker;
7d62fb64 163
121a60bd 164}
165
166//_____________________________________________________________________________
d76c31f4 167AliVertexer* AliITSReconstructor::CreateVertexer() const
121a60bd 168{
169// create a ITS vertexer
170
ff44c37c 171 AliITSVertexer *vptr;
9f9cae94 172 Int_t vtxOpt = GetRecoParam()->GetVertexer();
173 if(vtxOpt==3){
121a60bd 174 Info("CreateVertexer","a AliITSVertexerIons object has been selected\n");
ff44c37c 175 vptr = new AliITSVertexerIons();
121a60bd 176 }
ff44c37c 177 else if(vtxOpt==4){
121a60bd 178 Double_t smear[3]={0.005,0.005,0.01};
179 Info("CreateVertexer","a AliITSVertexerFast object has been selected\n");
ff44c37c 180 vptr = new AliITSVertexerFast(smear);
121a60bd 181 }
ff44c37c 182 else if(vtxOpt==1){
99f38bf6 183 Info("CreateVertexer","a AliITSVertexerZ object has been selected\n");
ff44c37c 184 vptr = new AliITSVertexerZ();
3acc14d5 185 }
ff44c37c 186 else if(vtxOpt==2){
3acc14d5 187 Info("CreateVertexer","a AliITSVertexerCosmics object has been selected\n");
ff44c37c 188 vptr = new AliITSVertexerCosmics();
121a60bd 189 }
ff44c37c 190 else if(vtxOpt==5){
c8735dd3 191 Info("CreateVertexer","vertex is fixed in the position of the TDI\n");
ff44c37c 192 vptr = new AliITSVertexerFixed("TDI");
c8735dd3 193 }
ff44c37c 194 else if(vtxOpt==6){
c8735dd3 195 Info("CreateVertexer","vertex is fixed in the position of the TED\n");
ff44c37c 196 vptr = new AliITSVertexerFixed("TED");
c8735dd3 197 }
ff44c37c 198 else {
99f38bf6 199 // by default an AliITSVertexer3D object is instatiated
ff44c37c 200 Info("CreateVertexer","a AliITSVertexer3D object has been selected\n");
201 AliITSVertexer3D* vtxr = new AliITSVertexer3D();
202 Float_t dzw=GetRecoParam()->GetVertexer3DWideFiducialRegionZ();
203 Float_t drw=GetRecoParam()->GetVertexer3DWideFiducialRegionR();
204 vtxr->SetWideFiducialRegion(dzw,drw);
205 Float_t dzn=GetRecoParam()->GetVertexer3DNarrowFiducialRegionZ();
206 Float_t drn=GetRecoParam()->GetVertexer3DNarrowFiducialRegionR();
207 vtxr->SetNarrowFiducialRegion(dzn,drn);
208 Float_t dphil=GetRecoParam()->GetVertexer3DLooseDeltaPhiCut();
209 Float_t dphit=GetRecoParam()->GetVertexer3DTightDeltaPhiCut();
210 vtxr->SetDeltaPhiCuts(dphil,dphit);
211 Float_t dcacut=GetRecoParam()->GetVertexer3DDCACut();
212 vtxr->SetDCACut(dcacut);
213 vptr = vtxr;
214 }
215 vptr->SetDetTypeRec(fDetTypeRec);
216 return vptr;
121a60bd 217}
218
219//_____________________________________________________________________________
d76c31f4 220void AliITSReconstructor::FillESD(TTree * /*digitsTree*/, TTree *clustersTree,
af885e0f 221 AliESDEvent* esd) const
121a60bd 222{
e62c1aea 223// make PID, find V0s and cascade
d76c31f4 224 if(fItsPID!=0) {
9f9cae94 225 Int_t pidOpt = GetRecoParam()->GetPID();
226 if(pidOpt==1){
d76c31f4 227 fItsPID->MakePID(clustersTree,esd);
9f9cae94 228 }else{
d76c31f4 229 fItsPID->MakePID(esd);
230 }
e62c1aea 231 }
232 else {
233 Error("FillESD","!! cannot do the PID !!\n");
234 }
121a60bd 235}