/[PAMELA software]/DarthVader/TrackerLevel2/src/TrkLevel0.cpp
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/TrkLevel0.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.3 by pam-fi, Thu Mar 15 12:17:10 2007 UTC revision 1.5 by pam-fi, Tue Nov 25 14:41:36 2008 UTC
# Line 4  Line 4 
4   */   */
5  #include <TrkLevel0.h>  #include <TrkLevel0.h>
6    
7    //......................................
8    // F77 routines
9    //......................................
10    extern "C" {    
11    
12        void filladc_(int*);
13        void evaluatecn_(int*);
14        void subtractped_(int*);
15    
16    
17    }
18  using namespace pamela::tracker;  using namespace pamela::tracker;
19    
20  /**  /**
21     *
22     *   Function to decode the tracker words.
23     *  
24     *   Tracker words are of three types:
25     *   - ADC data
26     *   - strip address
27     *   - end of ladder
28     *  
29     *   The function returns a struct variable (trkword)
30     *   that contains two quantities, type and decode, which are assigned
31     *   the following values:
32     *  
33     *   type                     decode              
34     *   ----------------------------------------------------
35     *   -1     error             0          
36     *    0     data              0-4095   ADC values
37     *    1     address           0-1023   strip address (relative to ladder)
38     *    2     end-of-ladder     1,2,3    ladder number     (compressed acq mode)      
39     *                            4,5,6    ladder number + 3 (full acq mode)
40     *      
41     **/
42    
43    trkword datadecode(int word){
44      int type =0;
45      int data =0;
46      int nodata = word&0xf000;
47      int zero = 0;
48    
49      trkword thisword;
50        
51      switch (nodata>>12){
52    
53      case 0:      
54            
55        thisword.decode = word;
56        thisword.type = 0;
57        //  cout << thisword.decode << "\n";
58        return (thisword);              //>>>>> ADC data (0)
59    
60      case 1:      
61    
62        type = word&0xC00;
63        data = word&0x3ff;
64    
65        switch(type>>10){
66    
67        case 0:  
68          thisword.decode = data;
69          thisword.type = 1; //>>> address (1)
70          return (thisword);          
71    
72        case 2:  
73          //            if(data>=4)data = data-3;
74          if(data>6){
75            printf("Error on data \n");
76            thisword.decode = zero;
77            thisword.type = -1;
78            return (thisword);     //>>>>> error (-1)
79          }
80          thisword.decode = data;      
81          thisword.type = 2;
82          return (thisword);       //>>> end-of-ladder
83    
84        default:
85          printf("Error on data \n");
86          thisword.decode = zero;
87          thisword.type = -1;
88          return (thisword);              //>>>>> error (-1)
89                
90        }
91    
92      default:
93    
94        printf("Error on data \n");
95        thisword.decode = zero;
96        thisword.type = -1;
97        return (thisword);              //>>>>> error (-1)
98      }
99    }
100    
101    
102    /**
103   * Pass values from the TrkLevel0 object to a struct cTrkLevel0 (to put data in F77 common).   * Pass values from the TrkLevel0 object to a struct cTrkLevel0 (to put data in F77 common).
104   */   */
105  //void TrkLevel0::GetCommonVar(cTrkLevel0 *l0) {  //void TrkLevel0::GetCommonVar(cTrkLevel0 *l0) {
# Line 89  int TrkLevel0::ProcessEvent(){ Line 182  int TrkLevel0::ProcessEvent(){
182      return 1;      return 1;
183    
184  }  }
185    /**
186     * Method to evaluate the raw signal (ADC) of each strip (it performs uncompression).
187     * (it calls F77 routine filladc )
188     */
189    bool TrkLevel0::FillADC(){
190    
191    //    cout << "int TrkLevel0::ProcessEvent()" << endl;
192        TrkParams::Load(6);
193        if( !TrkParams::IsLoaded(6) ){
194            cout << "int TrkLevel0::FillADC() -- ERROR -- VK-mask not loaded"<<endl;
195            return false;
196        };
197        TrkParams::LoadCalib( );
198        if( !TrkParams::CalibIsLoaded() ){
199            cout << "int TrkLevel0::FillADC() -- ERROR -- Calibration not loaded"<<endl;
200            return false;
201        };
202    
203        GetLevel0Struct();
204        int iflag=0;
205        filladc_(&iflag);
206        if(iflag!=0)return false;
207        return true;
208    
209    }
210    /**
211     * Method to evaluate calibrated signal (ADC-CN-PED) of each view
212     * (it calls F77 routines: filladc + evaluatecn + subtractped)
213     * @param iview view number (1-12)
214     * @param graph pointer to a TGraph to retrieve the output
215     * @param filladc if true calls F77 routine FillADC. otherwise must be called outside
216     */
217    bool TrkLevel0::GetCalibratedEvent(int iview, TGraph* graph, bool filladc){
218    
219        if ( iview<1 || iview>12 )return false;
220        if( filladc ){
221            if ( !FillADC() ) return false;
222        }
223        evaluatecn_(&iview);
224        subtractped_(&iview);
225    
226        if(graph){
227            graph->Set(NSTRIP);
228            for(int i=0; i<NSTRIP; i++)graph->SetPoint(i,(float)i,calibratedsignal_.value[i]);
229        }
230        
231    
232        return true;
233    
234    }
235    /**
236     * extract compressed/full data for a subset of channels
237     *
238     */
239    void TrkLevel0::Decode(int from,int to,bool& COMPRESSED,bool& FULL, int* datacomp, int* datafull){
240        //
241        int address = 0;
242        int ladder = 1;
243        Int_t whisto[3072]; for(int i=0; i<3072; i++)whisto[i] = -200;
244        //
245        COMPRESSED=false;
246        FULL=false;
247    
248        for(Int_t vi = from ; vi < to ; vi++){
249            int word = yodaobj->TrackerData.At(vi);    
250            trkword thisword = datadecode(word);
251            switch (thisword.type){
252                
253            case 0:  //ADC value
254                whisto[address] = thisword.decode;
255    //          cout << vi<<"    data " << thisword.decode << " @ "<<address << "\n";
256                address++;    
257                break;
258            
259            case 1:  //address
260                address = 1024*(ladder-1) + thisword.decode;
261    //          cout << vi<<"     adr " << address << "\n";
262                break;
263                
264            case 2:  //end-of-ladder
265                ladder = thisword.decode;
266    //          cout << vi<<" Ladder " << ladder << "\n";
267                if(ladder==3){
268                    //                  end of compressed data - FILL HISTO
269    //              cout << ">>> COMPRESSED data" << "\n";
270                    COMPRESSED=true;
271                    for(int ii = 0; ii < 3072; ii++){
272                        *(datacomp+ii) = whisto[ii];
273                        whisto[ii] = -200;
274                    }
275                    address = 0;
276                    ladder = 1;                
277                }else if(ladder==6){
278                    //                  end of full data - FILL HISTO
279    //              cout << ">>> FULL data" << "\n";
280                    FULL=true;
281                    for(int ii = 0; ii < 3072; ii++){
282                        *(datafull+ii) = whisto[ii];
283                        whisto[ii] = -200;
284                    }
285                    address = 0;
286                    ladder = 1;
287                }else{              
288                    if(ladder>3)    ladder=ladder-3;
289                    address= ladder*1024;          
290                    ladder = ladder + 1;    
291                }
292            }    
293        }    
294    
295    }
296    
297    /**
298     * Retrieve full event
299     *
300     */
301    bool TrkLevel0::GetFullEvent(int iview, TGraph* graph){
302    
303    
304        // retrieve the index of the dsp
305        int idsp=0;
306        do{
307            if ( yodaobj->DSPnumber[idsp] == iview )break;
308            idsp++;
309        }while(idsp<12);
310        if(idsp==12)return false;
311        //
312        // check DAQ mode:
313        //  9 full
314        // 10 compressed
315        // 11 compressed+full
316        //  8 special: compressed+full <--> compressed alternati
317        //
318        if( yodaobj->DAQmode[idsp] == 10 )return false;//compressed only
319        //
320        int offset =0;
321        for(int i=0; i<idsp; i++)offset+=yodaobj->DATAlength[i];
322        //
323        // decode data
324        //
325        Int_t whistocomp[3072];
326        Int_t whistofull[3072];
327        bool COMPRESSED=false;
328        bool FULL=false;
329        Decode(offset,offset+yodaobj->DATAlength[idsp],COMPRESSED,FULL,whistocomp,whistofull);
330        if(!FULL) return false;
331        
332        if(!graph)graph=new TGraph(3072);
333        for(int i=0; i<3072; i++)graph->SetPoint(i,(double)i,(double)whistofull[i]);
334        return true;
335    
336    }
337    
338    /**
339     * Retrieve compressed event
340     *
341     */
342    bool TrkLevel0::GetCompressedEvent(int iview, TGraph* graph){
343    
344    
345        // retrieve the index of the dsp
346        int idsp=0;
347        do{
348            if ( yodaobj->DSPnumber[idsp] == iview )break;
349            idsp++;
350        }while(idsp<12);
351        if(idsp==12)return false;
352        //
353        // check DAQ mode:
354        //  9 full
355        // 10 compressed
356        // 11 compressed+full
357        //  8 special: compressed+full <--> compressed alternati
358        //
359        if( yodaobj->DAQmode[idsp] == 9 )return false;
360        //
361        int offset =0;
362        for(int i=0; i<idsp; i++)offset+=yodaobj->DATAlength[i];
363        //
364        // decode data
365        //
366        Int_t whistocomp[3072];
367        Int_t whistofull[3072];
368        bool COMPRESSED=false;
369        bool FULL=false;
370        Decode(offset,offset+yodaobj->DATAlength[idsp],COMPRESSED,FULL,whistocomp,whistofull);
371        if(!COMPRESSED) return false;
372        
373        if(!graph)graph=new TGraph(3072);
374        for(int i=0; i<3072; i++)graph->SetPoint(i,(double)i,(double)whistocomp[i]);
375        return true;
376    
377    }
378    /**
379     * Retrieve uncompressed event event
380     */
381    bool TrkLevel0::GetUnCompressedEvent(int iview, TGraph* graph){
382    
383        if(!graph)graph=new TGraph(3072);
384        //
385        if( !GetCompressedEvent(iview,graph) )return false;
386        //  
387        double* adc;
388        double last;
389        adc = graph->GetY();
390        for(int i=0; i<3072; i++){
391            int ivk = (int)((float)i/128.);
392            int is = i%128;
393            if(*(adc+i)>0)last = adc[i] - (double)pedsigbad_.pedestal_t[is][ivk][iview-1];
394            else{
395                *(adc+i) = last + (double)pedsigbad_.pedestal_t[is][ivk][iview-1];
396            }
397            graph->SetPoint(i,(double)i,*(adc+i));
398        }
399        
400        return true;
401    
402    }
403    /**
404     * Retrieve sigmas
405     */
406    bool TrkLevel0::GetSigma(int iview, TGraph* graph){
407    
408        if(!graph)graph=new TGraph(3072);
409        //
410        for(int i=0; i<3072; i++){
411            int ivk = (int)((float)i/128.);
412            int is = i%128;
413            graph->SetPoint(i,(double)i,(double)pedsigbad_.sigma[is][ivk][iview-1]);
414        }
415        
416        return true;
417    
418    }
419    
420  ClassImp(TrkLevel0);  ClassImp(TrkLevel0);
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23