]>
Commit | Line | Data |
---|---|---|
325500ed | 1 | // |
2 | // Procedures to create a cluster parametererization | |
3 | // Currently only pseudo code - once will be stable - to be updated as a "normal" macro | |
4 | // | |
5 | ||
6 | void PseudoCode(){ | |
7 | // | |
8 | // That rough sequence to update a Cluster param calibration using debug streamers | |
9 | // | |
10 | ||
11 | // | |
12 | //0. Load libraries | |
13 | gSystem->Load("libANALYSIS"); | |
14 | gSystem->Load("libSTAT"); | |
15 | gSystem->Load("libTPCcalib"); | |
16 | ||
17 | ||
18 | //1. Load Parameters to be modified: | |
19 | //e.g: | |
20 | AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT"); | |
21 | AliCDBManager::Instance()->SetRun(0) ; | |
22 | AliTPCClusterParam * param = AliTPCcalibDB::Instance()->GetClusterParam(); | |
23 | // | |
24 | //2. Load chain from debug streamers | |
25 | // | |
26 | //e.g | |
27 | gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/macros"); | |
28 | gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+"); | |
29 | AliXRDPROOFtoolkit tool; | |
30 | TChain * chainres = tool.MakeChain("tracks.txt","ResolCl",0,10200); | |
31 | chainres->Lookup(); | |
32 | TChain * chainGain = tool.MakeChain("gain.txt","dEdx",0,1000000); | |
33 | chainGain->Lookup(); | |
34 | TChain * chainCosmic = tool.MakeChain("cosmic.txt","Track0",0,1000000); | |
35 | chainCosmic->Lookup(); | |
36 | // | |
37 | // 3. Do fits - pos correction | |
38 | // | |
39 | AliTPCcalibTracks::MakeQPosNormAll(chainres,param,200000,0) ; | |
40 | ||
41 | // | |
42 | //4. Do fits gain | |
43 | // | |
44 | param->fPosQTnorm[0] = AliTPCcalibTracksGain::MakeQPosNorm(chainGain,0,kFALSE,200000,kTRUE); | |
45 | param->fPosQTnorm[1] = AliTPCcalibTracksGain::MakeQPosNorm(chainGain,1,kFALSE,200000,kTRUE); | |
46 | param->fPosQTnorm[2] = AliTPCcalibTracksGain::MakeQPosNorm(chainGain,2,kFALSE,200000,kTRUE); | |
47 | // | |
48 | param->fPosQMnorm[0] = AliTPCcalibTracksGain::MakeQPosNorm(chainGain,0,kTRUE,200000,kTRUE); | |
49 | param->fPosQMnorm[1] = AliTPCcalibTracksGain::MakeQPosNorm(chainGain,1,kTRUE,200000,kTRUE); | |
50 | param->fPosQMnorm[2] = AliTPCcalibTracksGain::MakeQPosNorm(chainGain,2,kTRUE,200000,kTRUE); | |
51 | // | |
52 | //5. Store results | |
53 | // | |
54 | TFile f("paramout.root","recreate"); | |
55 | param->Write("clusterParam"); | |
56 | f.Close(); | |
57 | ||
58 | ||
59 | //6. Verification | |
60 | TFile f2("paramout.root"); | |
61 | AliTPCClusterParam *param2 = (AliTPCClusterParam*)f2.Get("clusterParam"); | |
62 | param2->SetInstance(param2); | |
63 | // | |
64 | // position correction | |
65 | chainres->Draw("AliTPCClusterParam::SPosCorrection(1,0,Cl.fPad,Cl.fTimeBin,Cl.fZ,Cl.fSigmaY2,Cl.fSigmaZ2,Cl.fMax)","Cl.fDetector<36","",10000); // should be line | |
66 | // gain correction | |
67 | TCut cutA("dedge>3&&fraction2<0.7"); | |
68 | chainGain->Draw("(Cl.fMax/gain/dedxM.fElements[2]):AliTPCClusterParam::SQnormPos(0,1,Cl.fPad,Cl.fTimeBin,Cl.fZ,Cl.fSigmaY2,Cl.fSigmaZ2,Cl.fMax,Cl.fQ)","IPad==0"+cutA,"prof",100000); | |
69 | ||
70 | ||
71 | } | |
72 | ||
73 | ||
74 | void UpdateParam(){ | |
75 | // | |
76 | // Pseudo code -to update cluster params | |
77 | // | |
78 | .L $ALICE_ROOT/TPC/Cal/AliTPCCreateDummyCDB.C | |
79 | TFile f2("paramout.root"); | |
80 | AliTPCClusterParam *param2 = (AliTPCClusterParam*)f2.Get("clusterParam"); | |
81 | AliCDBManager *man = AliCDBManager::Instance(); | |
82 | gStorLoc = man->GetStorage(gCDBpath); | |
83 | metaData = CreateMetaObject("AliTPCClusterParam"); | |
84 | StoreObject("TPC/Calib/ClusterParam/", param2, metaData); | |
85 | ||
86 | } | |
87 | ||
88 | ||
89 | void TestChainCosmicDedx(){ | |
90 | // | |
91 | // pseudo cose | |
92 | // | |
93 | gSystem->Load("libANALYSIS"); | |
94 | gSystem->Load("libSTAT"); | |
95 | gSystem->Load("libTPCcalib"); | |
96 | AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT"); | |
97 | AliCDBManager::Instance()->SetRun(0) ; | |
98 | AliTPCClusterParam * param = AliTPCcalibDB::Instance()->GetClusterParam(); | |
99 | param->SetInstance(param); | |
100 | //e.g | |
101 | gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/macros"); | |
102 | gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+"); | |
103 | AliXRDPROOFtoolkit tool; | |
104 | TChain * chainCosmic = tool.MakeChain("cosmic.txt","Track0",0,1000000); | |
105 | chainCosmic->Lookup(); | |
106 | ||
107 | ||
108 | ||
109 | TCut cutT("cutT","abs(Tr1.fP[3]+Tr0.fP[3])<0.01"); // OK | |
110 | TCut cutD("cutD","abs(Tr0.fP[0]+Tr1.fP[0])<2"); // OK | |
111 | TCut cutP1("cutP1","abs(Tr0.fP[1]-Tr1.fP[1])<3"); // OK | |
112 | TCut cutPt("cutPt","abs(Tr1.fP[4]+Tr0.fP[4])<0.1&&abs(Tr0.fP[4])+abs(Tr1.fP[4])<10"); | |
113 | TCut cutN("cutN","min(Orig0.fTPCncls,Orig1.fTPCncls)>100"); | |
114 | TCut cutA=cutT+cutD+cutPt+cutN+cutP1; | |
115 | ||
116 | TCut cutS("cutS","Orig0.fIp.fP[1]*Orig1.fIp.fP[1]>0"); | |
117 | ||
118 | if (gProof) chainCosmic->SetProof(kTRUE); | |
119 | // | |
120 | // | |
121 | // | |
122 | // | |
123 | chainCosmic->Draw("Seed0.CookdEdxNorm(0.01,0.65,1,0,159,0)/Seed1.CookdEdxNorm(0.01,0.65,1,0,159,0)","min(Orig0.fTPCncls,Orig1.fTPCncls)>120","",10000); | |
124 | ||
125 | // | |
126 | // Q normalization - short - medium - long | |
127 | // | |
128 | ||
129 | TVectorD paramT(3); | |
130 | TVectorD paramM(3); | |
131 | TH1F * his =0; | |
132 | chainCosmic->Draw("Seed0.CookdEdxNorm(0.01,0.65,0,0,64,0,kTRUE)/Seed0.CookdEdxNorm(0.01,0.65,0,0,159,0,kTRUE)>>his(100,0.5,1.5)","min(Orig0.fTPCncls,Orig1.fTPCncls)>140","",100000); | |
133 | TH1 * his = (TH1*)gROOT->FindObject("his"); | |
134 | paramT[0]=his->GetMean(); | |
135 | // | |
136 | chainCosmic->Draw("Seed0.CookdEdxNorm(0.01,0.65,0,64,129,0,kTRUE)/Seed0.CookdEdxNorm(0.01,0.65,0,0,159,0,kTRUE)>>his(100,0.5,1.5)","min(Orig0.fTPCncls,Orig1.fTPCncls)>140","",100000); | |
137 | TH1 * his = (TH1*)gROOT->FindObject("his"); | |
138 | paramT[1]=his->GetMean(); | |
139 | // | |
140 | chainCosmic->Draw("Seed0.CookdEdxNorm(0.01,0.65,0,129,159,0,kTRUE)/Seed0.CookdEdxNorm(0.01,0.65,0,0,159,0,kTRUE)>>his(100,0.5,2.0)","min(Orig0.fTPCncls,Orig1.fTPCncls)>140","",100000); | |
141 | TH1 * his = (TH1*)gROOT->FindObject("his"); | |
142 | paramT[2]=his->GetMean(); | |
143 | // | |
144 | // | |
145 | chainCosmic->Draw("Seed0.CookdEdxNorm(0.01,0.65,1,0,64,0,kTRUE)/Seed0.CookdEdxNorm(0.01,0.65,1,0,159,0,kTRUE)>>his(100,0.5,1.5)","min(Orig0.fTPCncls,Orig1.fTPCncls)>140","",100000); | |
146 | TH1 * his = (TH1*)gROOT->FindObject("his"); | |
147 | paramM[0]=his->GetMean(); | |
148 | // | |
149 | chainCosmic->Draw("Seed0.CookdEdxNorm(0.01,0.65,1,64,129,0,kTRUE)/Seed0.CookdEdxNorm(0.01,0.65,1,0,159,0,kTRUE)>>his(100,0.5,1.5)","min(Orig0.fTPCncls,Orig1.fTPCncls)>140","",100000); | |
150 | TH1 * his = (TH1*)gROOT->FindObject("his"); | |
151 | paramM[1]=his->GetMean(); | |
152 | // | |
153 | chainCosmic->Draw("Seed0.CookdEdxNorm(0.01,0.65,1,129,159,0,kTRUE)/Seed0.CookdEdxNorm(0.01,0.65,1,0,159,0,kTRUE)>>his(100,0.5,2.0)","min(Orig0.fTPCncls,Orig1.fTPCncls)>140","",100000); | |
154 | TH1 * his = (TH1*)gROOT->FindObject("his"); | |
155 | paramM[2]=his->GetMean(); | |
156 | ||
157 | param->fQpadTnorm=(TVectorD*)paramT->Clone(); | |
158 | param->fQpadMnorm=(TVectorD*)paramM->Clone(); | |
159 | ||
160 | TFile f("paramout.root","recreate"); | |
161 | param->Write("clusterParam"); | |
162 | f.Close(); | |
163 | ||
164 | } | |
ee59d24d | 165 | |
166 | ||
167 | void dedxDemo(){ | |
168 | chainCosmic->Draw("(Seed0.CookdEdxNorm(0.01,0.65,1,0,159,0,kTRUE,kTRUE)+Seed1.CookdEdxNorm(0.01,0.65,1,0,159,0,kTRUE,kTRUE))*0.5:log(Tr0.P())>>his2(100,-2,5,500,0,200)","min(Orig0.fTPCncls,Orig1.fTPCncls)>110"); | |
169 | } |