more secure string operations
[u/mrichter/AliRoot.git] / TPC / AliTPCcalibAlignment.cxx
CommitLineData
9318a5b4 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
17///////////////////////////////////////////////////////////////////////////////
18// //
19// Class to make and monitor alignment of the TPC //
20//
21//
22
23
24#include "AliTracker.h"
25#include "AliTPCcalibAlignment.h"
26#include "AliTPCseed.h"
27#include "AliTPCclusterMI.h"
28#include "AliTPCClusterParam.h"
29#include "AliESDtrack.h"
30#include "TTreeStream.h"
31#include "AliTPCTracklet.h"
32
33#include <iostream>
34using namespace std;
35
36/* TEST
37
38 gSystem->Load("$ALICE_ROOT/TPC/TPCcalib/libTPCcalib.so");
39
40 .L AliXRDPROOFtoolkit.cxx+
41 AliXRDPROOFtoolkit tool;
42 TChain * chain = tool.MakeChain("listcosmic.txt","esdTree",0,100,0)
43//TChain * chain = tool.MakeChain("listpp.txt","esdTree",0,100,0)
44
45 AliTPCcalibTracks::AddInfo(chain,"TPCClusterParam.root");
46 AliTPCcalibTracks::AddCuts(chain,"LOWFLUX");
47
48// AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 0., 10., 2); // NO FIELD
49 AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 2); // FIELD
50// field->SetL3ConstField(0); //Using const. field in the barrel
51 AliTracker::SetFieldMap(field,kFALSE);
52
53 chain->SetBranchStatus("*",1);
54 chain->Process("$ALICE_ROOT/TPC/TPCcalib/AliTPCSelectorTracks.cxx+");
55
56 TFile f("AliTPCcalibAlignmentDebug.root")
57 TTree *t=f.Get("Alignment")
58
59 t->Draw("tracklet1.fInner.GetY()-trackletR1.fInner.GetY()")
60 t->Draw("tracklet1.fInner.GetSnp()-tracklet1.fOuter.GetSnp()")
61 t->Draw("common1.GetSnp()-common2.GetSnp()")
62 t->Draw("common1.GetY()-common2.GetY()")
63 t->Draw("common1.GetZ()-common2.GetZ()")
64
65 t->Draw("track1.GetY()-track2.GetY()","abs(track1.GetY()-track2.GetY())<0.2")
66 t->Draw("track1.GetZ()-track2.GetZ()","abs(track1.GetZ()-track2.GetZ())<0.2")
67 t->Draw("track1.GetY()-track2.GetY():track1.GetZ()-track2.GetZ()","abs(track1.GetY()-track2.GetY())<0.2&&abs(track1.GetZ()-track2.GetZ())<0.2")
68
69
70 */
71
72ClassImp(AliTPCcalibAlignment)
73
cec17745 74AliTPCcalibAlignment::AliTPCcalibAlignment():
75 fDebugStream(0)
76{
9318a5b4 77 //
78 // Constructor
79 //
80 fDebugStream=new TTreeSRedirector("AliTPCcalibAlignmentDebug.root");
81}
82
83AliTPCcalibAlignment::~AliTPCcalibAlignment() {
84 //
85 // Destructor
86 //
87 delete fDebugStream;
88}
89
90void AliTPCcalibAlignment::Process(AliTPCseed *track) {
91 //
92 // Track processing
93 //
94 TObjArray tracklets=
95 AliTPCTracklet::CreateTracklets(track,AliTPCTracklet::kKalman,kFALSE,20,2);
96 TObjArray trackletsL=
97 AliTPCTracklet::CreateTracklets(track,AliTPCTracklet::kLinear,kFALSE,20,2);
98 TObjArray trackletsQ=
99 AliTPCTracklet::CreateTracklets(track,AliTPCTracklet::kQuadratic,kFALSE,20,2);
100 TObjArray trackletsR=
101 AliTPCTracklet::CreateTracklets(track,AliTPCTracklet::kRiemann,kFALSE,20,2);
102 tracklets.SetOwner();
103 if (tracklets.GetEntries()==2) {
104 AliTPCTracklet* t1=static_cast<AliTPCTracklet*>(tracklets[0]);
105 AliTPCTracklet* t2=static_cast<AliTPCTracklet*>(tracklets[1]);
106 AliExternalTrackParam *common1=0,*common2=0;
107 if (AliTPCTracklet::PropagateToMeanX(*t1,*t2,common1,common2))
108 (*fDebugStream)<<"Alignment"
109 <<"common1.="<<common1
110 <<"common2.="<<common2
111 <<"tracklet1.="<<tracklets[0]
112 <<"tracklet2.="<<tracklets[1]
113 <<"trackletL1.="<<trackletsL[0]
114 <<"trackletL2.="<<trackletsL[1]
115 <<"trackletQ1.="<<trackletsQ[0]
116 <<"trackletQ2.="<<trackletsQ[1]
117 <<"trackletR1.="<<trackletsR[0]
118 <<"trackletR2.="<<trackletsR[1]
119 <<"\n";
120 delete common1;
121 delete common2;
122 }
123}
124
125