]>
Commit | Line | Data |
---|---|---|
c12208b8 | 1 | /* |
2 | ||
3 | TPCda_pedestal.cxx - calibration algorithm for TPC pedestal runs | |
4 | ||
5 | 10/06/2007 sylvain.chapeland@cern.ch : first version - clean skeleton based on DAQ DA case1 | |
fee88714 | 6 | 19/10/2007 christian.lippmann@cern.ch : Possibility to write output to ASCII file |
7 | 24/10/2007 christian.lippmann@cern.ch : Including pedestal calibration for time bins | |
8 | 23/11/2007 christian.lippmann@cern.ch : Fix in order to avoid streamer problems in case of | |
9 | invalid ROOTSTYS. The famous magic line provided by Rene. | |
10 | 28/11/2007 christian.lippmann@cern.ch : TPC mapping file is read from DaqDetDB | |
c12208b8 | 11 | |
12 | contact: marian.ivanov@cern.ch | |
13 | ||
c12208b8 | 14 | This process reads RAW data from the files provided as command line arguments |
15 | and save results in a file (named from RESULT_FILE define - see below). | |
16 | ||
17 | */ | |
18 | ||
fee88714 | 19 | #define RESULT_FILE "tpcPedestal.root" |
20 | #define MAPPING_FILE "tpcMapping.root" | |
c12208b8 | 21 | |
22 | ||
23 | extern "C" { | |
24 | #include <daqDA.h> | |
25 | } | |
26 | #include "event.h" | |
27 | #include "monitor.h" | |
fee88714 | 28 | |
29 | #include "stdio.h" | |
30 | #include "stdlib.h" | |
31 | #include <fstream> | |
c12208b8 | 32 | |
33 | // | |
34 | //Root includes | |
35 | // | |
fee88714 | 36 | #include "TFile.h" |
37 | #include "TArrayF.h" | |
8b1405a7 | 38 | #include "TROOT.h" |
39 | #include "TPluginManager.h" | |
c12208b8 | 40 | |
41 | // | |
42 | //AliRoot includes | |
43 | // | |
44 | #include "AliRawReader.h" | |
45 | #include "AliRawReaderDate.h" | |
87d57858 | 46 | #include "AliTPCmapper.h" |
c12208b8 | 47 | #include "AliTPCRawStream.h" |
48 | #include "AliTPCROC.h" | |
49 | #include "AliTPCCalROC.h" | |
50 | #include "AliTPCCalPad.h" | |
51 | #include "AliMathBase.h" | |
52 | #include "TTreeStream.h" | |
53 | ||
54 | // | |
55 | // TPC calibration algorithm includes | |
56 | // | |
57 | #include "AliTPCCalibPedestal.h" | |
58 | ||
bdf99a93 | 59 | /* |
60 | Main routine, TPC pedestal detector algorithm to be run on TPC LDC | |
61 | Arguments: list of DATE raw data files | |
c12208b8 | 62 | */ |
bdf99a93 | 63 | |
c12208b8 | 64 | int main(int argc, char **argv) { |
bdf99a93 | 65 | // |
66 | // Main for TPC pedestal detector algorithm | |
67 | // | |
8b1405a7 | 68 | |
bdf99a93 | 69 | Bool_t timeAnalysis = kTRUE; |
c12208b8 | 70 | |
c12208b8 | 71 | if (argc<2) { |
72 | printf("Wrong number of arguments\n"); | |
73 | return -1; | |
74 | } | |
75 | ||
fee88714 | 76 | /* magic line */ |
77 | gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo", | |
78 | "*", | |
79 | "TStreamerInfo", | |
80 | "RIO", | |
81 | "TStreamerInfo()"); | |
82 | int i, status; | |
83 | ||
c12208b8 | 84 | /* log start of process */ |
85 | printf("TPC DA started - %s\n",__FILE__); | |
86 | ||
c12208b8 | 87 | /* declare monitoring program */ |
88 | status=monitorDeclareMp( __FILE__ ); | |
89 | if (status!=0) { | |
90 | printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status)); | |
91 | return -1; | |
92 | } | |
93 | ||
fee88714 | 94 | /* copy locally the mapping file from daq detector config db */ |
95 | status = daqDA_DB_getFile(MAPPING_FILE,"./tpcMapping.root"); | |
96 | if (status) { | |
97 | printf("Failed to get mapping file (%s) from DAQdetDB, status=%d\n", MAPPING_FILE, status); | |
4958b652 | 98 | printf("Continue anyway ... maybe it works?\n"); // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
99 | //return -1; // temporarily uncommented for testing on pcald47 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! | |
fee88714 | 100 | } |
101 | ||
102 | /* open the mapping file and retrieve mapping object */ | |
fee88714 | 103 | AliTPCmapper *mapping = 0; // The TPC mapping |
4958b652 | 104 | TFile *fileMapping = new TFile(MAPPING_FILE, "read"); |
fee88714 | 105 | mapping = (AliTPCmapper*) fileMapping->Get("tpcMapping"); |
106 | if (mapping == 0) { | |
107 | printf("Failed to get mapping object from %s. Exiting ...\n", MAPPING_FILE); | |
108 | delete fileMapping; | |
109 | return -1; | |
110 | } else { | |
111 | printf("Got mapping object from %s\n", MAPPING_FILE); | |
112 | } | |
fee88714 | 113 | |
4958b652 | 114 | AliTPCCalibPedestal calibPedestal; // pedestal and noise calibration |
115 | calibPedestal.SetTimeAnalysis(timeAnalysis); // pedestal(t) calibration | |
116 | calibPedestal.SetAltroMapping(mapping->GetAltroMapping()); // Use altro mapping we got from daqDetDb | |
fee88714 | 117 | |
c12208b8 | 118 | /* loop over RAW data files */ |
119 | int nevents=0; | |
bdf99a93 | 120 | for ( i=1; i<argc; i++ ) { |
c12208b8 | 121 | |
122 | /* define data source : this is argument i */ | |
123 | printf("Processing file %s\n", argv[i]); | |
124 | status=monitorSetDataSource( argv[i] ); | |
125 | if (status!=0) { | |
fee88714 | 126 | printf("monitorSetDataSource() failed. Error=%s. Exiting ...\n", monitorDecodeError(status)); |
127 | delete fileMapping; | |
c12208b8 | 128 | return -1; |
129 | } | |
130 | ||
131 | /* read until EOF */ | |
bdf99a93 | 132 | for ( ; ; ) { |
c12208b8 | 133 | struct eventHeaderStruct *event; |
134 | ||
135 | /* check shutdown condition */ | |
136 | if (daqDA_checkShutdown()) {break;} | |
137 | ||
138 | /* get next event (blocking call until timeout) */ | |
139 | status=monitorGetEventDynamic((void **)&event); | |
140 | if (status==MON_ERR_EOF) { | |
bdf99a93 | 141 | printf ("End of File %d (%s) detected\n", i, argv[i]); |
142 | break; /* end of monitoring file has been reached */ | |
c12208b8 | 143 | } |
c12208b8 | 144 | if (status!=0) { |
bdf99a93 | 145 | printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status)); |
146 | break; | |
c12208b8 | 147 | } |
148 | ||
bdf99a93 | 149 | /* skip start/end of run events */ |
150 | if ( (event->eventType != physicsEvent) && (event->eventType != calibrationEvent) ) | |
151 | continue; | |
152 | ||
c12208b8 | 153 | /* retry if got no event */ |
4958b652 | 154 | if (event==NULL) |
bdf99a93 | 155 | continue; |
bdf99a93 | 156 | |
c12208b8 | 157 | nevents++; |
158 | ||
159 | // Pedestal calibration | |
160 | AliRawReader *rawReader = new AliRawReaderDate((void*)event); | |
4958b652 | 161 | //calibPedestal.ProcessEvent(rawReader); |
162 | calibPedestal.ProcessEventFast(rawReader); // fast data reader | |
c12208b8 | 163 | delete rawReader; |
164 | ||
165 | /* free resources */ | |
166 | free(event); | |
167 | } | |
168 | } | |
169 | ||
fee88714 | 170 | // |
171 | // Analyse pedestals and write them to rootfile | |
172 | // | |
173 | ||
174 | calibPedestal.Analyse(); | |
175 | calibPedestal.AnalyseTime(nevents); | |
176 | printf ("%d physics/calibration events processed.\n",nevents); | |
c12208b8 | 177 | |
bdf99a93 | 178 | TFile *fileTPC = new TFile(RESULT_FILE, "recreate"); |
c12208b8 | 179 | calibPedestal.Write("calibPedestal"); |
180 | delete fileTPC; | |
fee88714 | 181 | printf("Wrote %s.\n",RESULT_FILE); |
c12208b8 | 182 | |
bdf99a93 | 183 | // |
fee88714 | 184 | // Now prepare ASCII files for local ALTRO configuration through DDL. |
bdf99a93 | 185 | // |
87d57858 | 186 | |
bdf99a93 | 187 | ofstream pedfile; |
188 | ofstream noisefile; | |
189 | ofstream pedmemfile; | |
87d57858 | 190 | char filename[255]; |
bdf99a93 | 191 | sprintf(filename,"tpcPedestals.data"); |
192 | pedfile.open(filename); | |
193 | sprintf(filename,"tpcNoise.data"); | |
194 | noisefile.open(filename); | |
195 | sprintf(filename,"tpcPedestalMem.data"); | |
196 | pedmemfile.open(filename); | |
197 | ||
fee88714 | 198 | TArrayF **timePed = calibPedestal.GetTimePedestals(); // pedestal values for each time bin |
87d57858 | 199 | |
bdf99a93 | 200 | Int_t ctr_channel = 0; |
201 | Int_t ctr_altro = 0; | |
202 | Int_t ctr_pattern = 0; | |
87d57858 | 203 | |
fee88714 | 204 | pedfile << 10 << std::endl; // mark file to contain PEDESTALS per channel |
205 | noisefile << 11 << std::endl; // mark file to contain NOISE per altro | |
206 | pedmemfile << 12 << std::endl; // mark file to contain PEDESTALs per time bin | |
87d57858 | 207 | |
bdf99a93 | 208 | for ( Int_t roc = 0; roc < 72; roc++ ) { |
209 | if ( !calibPedestal.GetCalRocPedestal(roc) ) continue; | |
fee88714 | 210 | Int_t side = mapping->GetSideFromRoc(roc); |
211 | Int_t sector = mapping->GetSectorFromRoc(roc); | |
bdf99a93 | 212 | //printf("Analysing ROC %d (side %d, sector %d) ...\n", roc, side, sector); |
fee88714 | 213 | Int_t nru = mapping->IsIROC(roc) ? 2 : 4; |
bdf99a93 | 214 | for ( int rcu = 0; rcu < nru; rcu++ ) { |
fee88714 | 215 | Int_t patch = mapping->IsIROC(roc) ? rcu : rcu+2; |
bdf99a93 | 216 | for ( int branch = 0; branch < 2; branch++ ) { |
fee88714 | 217 | for ( int fec = 0; fec < mapping->GetNfec(patch, branch); fec++ ) { |
bdf99a93 | 218 | for ( int altro = 0; altro < 8; altro++ ) { |
219 | Float_t rms = 0.; | |
220 | Float_t ctr = 0.; | |
221 | for ( int channel = 0; channel < 16; channel++ ) { | |
4958b652 | 222 | Int_t hwadd = mapping->CodeHWAddress(branch, fec, altro, channel); |
223 | Int_t row = mapping->GetPadRow(patch, hwadd); // row in a ROC | |
224 | Int_t globalrow = mapping->GetGlobalPadRow(patch, hwadd); // row in full sector | |
225 | Int_t pad = mapping->GetPad(patch, hwadd); | |
226 | Float_t ped = calibPedestal.GetCalRocPedestal(roc)->GetValue(row,pad); | |
bdf99a93 | 227 | // fixed pedestal |
228 | if ( ped > 1.e-10 ) { | |
fee88714 | 229 | pedfile << ctr_channel << "\t" << side << "\t" << sector << "\t" << patch << "\t" |
230 | << hwadd << "\t" << ped << std::endl; | |
bdf99a93 | 231 | ctr_channel++; |
232 | } | |
233 | // pedestal(t) | |
4958b652 | 234 | if ( timePed && fabs(timePed[globalrow][pad].GetSum()) > 1e-10 ) { |
fee88714 | 235 | pedmemfile << ctr_pattern << "\t" << side << "\t" << sector << "\t" << patch |
236 | << "\t" << hwadd; | |
bdf99a93 | 237 | for ( Int_t timebin = 0; timebin < 1024; timebin++ ) |
4958b652 | 238 | pedmemfile << "\t" << timePed[globalrow][pad].At(timebin); |
bdf99a93 | 239 | pedmemfile << std::endl; |
240 | ctr_pattern++; | |
241 | } | |
242 | // rms=noise | |
243 | Float_t rms2 = calibPedestal.GetCalRocRMS(roc)->GetValue(row,pad); | |
244 | if ( rms2 > 1.e-10 ) { rms += rms2; ctr += 1.; } | |
245 | } // end channel for loop | |
246 | // noise data (rms) averaged over all channels in this ALTRO. | |
fee88714 | 247 | Int_t hwadd = mapping->CodeHWAddress(branch, fec, altro, 0); |
bdf99a93 | 248 | if ( ctr > 1.e-10 ) { |
fee88714 | 249 | noisefile << ctr_altro << "\t" << side << "\t" << sector << "\t" << patch << "\t" |
250 | << hwadd << "\t" << rms/ctr << std::endl; | |
bdf99a93 | 251 | ctr_altro++; |
252 | } | |
253 | } // end altro for loop | |
254 | } // end fec for loop | |
255 | } // end branch for loop | |
256 | } // end rcu for loop | |
257 | } // end roc loop | |
258 | ||
259 | pedfile.close(); | |
260 | noisefile.close(); | |
261 | pedmemfile.close(); | |
fee88714 | 262 | printf("Wrote ASCII files.\n"); |
bdf99a93 | 263 | |
fee88714 | 264 | delete fileMapping; |
87d57858 | 265 | |
c12208b8 | 266 | return status; |
267 | } |