source: trunk/Mars/mtrigger/MTriggerPatternDecode.cc@ 10100

Last change on this file since 10100 was 9531, checked in by tbretz, 15 years ago
*** empty log message ***
File size: 5.4 KB
Line 
1/* ======================================================================== *\
2!
3! *
4! * This file is part of MARS, the MAGIC Analysis and Reconstruction
5! * Software. It is distributed to you in the hope that it can be a useful
6! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
7! * It is distributed WITHOUT ANY WARRANTY.
8! *
9! * Permission to use, copy, modify and distribute this software and its
10! * documentation for any purpose is hereby granted without fee,
11! * provided that the above copyright notice appear in all copies and
12! * that both that copyright notice and this permission notice appear
13! * in supporting documentation. It is provided "as is" without express
14! * or implied warranty.
15! *
16!
17!
18! Author(s): Nicola Galante 12/2004 <mailto:nicola.galante@pi.infn.it>
19! Author(s): Thomas Bretz 12/2004 <mailto:tbretz@astro.uni-wuerzburg.de>
20!
21! Copyright: MAGIC Software Development, 2004
22!
23\* ======================================================================== */
24
25/////////////////////////////////////////////////////////////////////////////
26//
27// MTriggerPatternDecode
28//
29// Decodes the trigger pattern from MRawEvtData into MTriggerPattern.
30//
31// For files before file version 5 the trigger pattern is set according
32// to the file type
33//
34// For more details see: MTriggerPattern
35//
36// Input:
37// MRawEvtData
38//
39// Output:
40// MTriggerPattern
41//
42/////////////////////////////////////////////////////////////////////////////
43#include "MTriggerPatternDecode.h"
44
45#include "MLog.h"
46#include "MLogManip.h"
47
48#include "MParList.h"
49#include "MRawEvtHeader.h"
50#include "MRawRunHeader.h"
51#include "MTriggerPattern.h"
52
53ClassImp(MTriggerPatternDecode);
54
55using namespace std;
56
57// --------------------------------------------------------------------------
58//
59// Default constructor
60//
61MTriggerPatternDecode::MTriggerPatternDecode(const char *name, const char *title)
62 : fRunHeader(0), fEvtHeader(0), fPattern(0)
63{
64 fName = name ? name : "MTriggerPatternDecode";
65 fTitle = title ? title : "Task to decode Trigger Pattern";
66}
67
68// --------------------------------------------------------------------------
69//
70Int_t MTriggerPatternDecode::PreProcess(MParList *pList)
71{
72 fRunHeader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
73 if (!fRunHeader)
74 {
75 *fLog << err << "MRawRunHeader not found... abort." << endl;
76 return kFALSE;
77 }
78
79 fEvtHeader = (MRawEvtHeader*)pList->FindObject("MRawEvtHeader");
80 if (!fEvtHeader)
81 {
82 *fLog << err << "MRawEvtHeader not found... abort." << endl;
83 return kFALSE;
84 }
85
86 fPattern = (MTriggerPattern*)pList->FindCreateObj("MTriggerPattern");
87 if (!fPattern)
88 return kFALSE;
89
90 return kTRUE;
91}
92
93// --------------------------------------------------------------------------
94//
95//
96Int_t MTriggerPatternDecode::Process()
97{
98 if (fRunHeader->GetFormatVersion()<5)
99 {
100 switch (fRunHeader->GetRunType()&0xff)
101 {
102 case MRawRunHeader::kRTData:
103 fPattern->fPrescaled = MTriggerPattern::kTriggerLvl1;
104 fPattern->fUnprescaled = MTriggerPattern::kTriggerLvl1;
105 return kTRUE;
106
107 case MRawRunHeader::kRTPedestal:
108 fPattern->fPrescaled = MTriggerPattern::kPedestal;
109 fPattern->fUnprescaled = MTriggerPattern::kPedestal;
110 return kTRUE;
111
112 case MRawRunHeader::kRTCalibration:
113 fPattern->fPrescaled = MTriggerPattern::kCalibration;
114 fPattern->fUnprescaled = MTriggerPattern::kCalibration;
115 return kTRUE;
116 }
117 return kTRUE;
118 }
119
120 const UInt_t pattern = ~fEvtHeader->GetTriggerID();
121
122 // The the trigger pattern is currently written with inverted bits,
123 // but in the future this could be changed. In this case the file version
124 // number of the raw-data format must be checked. It tells you the
125 // file format (meaning of the bits)
126 // If for some reason the file-format has not been correctly changed
127 // use the run-number to decide whether the bits must be inverted or not.
128
129 // Trigger Pattern is a number of 16 bits, but if the machine
130 // is a 64 bits then the bit inversion the first 16 bits are set to 1,
131 // possibly giving problems with the filtering with the masks.
132 // Because UInt_t by definition is a 32 bit number on any
133 // machine (see root manual) no check is needed.
134 // The simplest workaround is:
135 // pattern &= 0xffffffff;
136
137 // For the moment the meaning of the bits in
138 // MRawEvtHeader::GetTriggerID and in MTriggerPattern::fTriggerPattern
139 // are identical - this may change in the future! In this case
140 // the decoding HERE must be changed - NOT MTriggerPattern.
141 // If an enhancement of the information stored in MTriggerPattern
142 // is necessary ADD new variables! Don't change the meaning of the old
143 // ones!
144 fPattern->fPrescaled = pattern & 0xff;
145 fPattern->fUnprescaled = (pattern>>8) & 0xff;
146
147 // This is a workaround for the new scheme in which L1TPU (the signal
148 // comming directly from the L1 is connected, but the L1 (routed
149 // over L2 is disconnected)
150 if (!fRunHeader->IsMonteCarloRun() && fRunHeader->GetTelescopeNumber()==1 &&
151 fRunHeader->GetRunNumber()>1006246)
152 {
153 fPattern->fPrescaled |= (pattern>> 6)&1;
154 fPattern->fUnprescaled |= (pattern>>14)&1;
155
156 fPattern->fPrescaled &= 0xbf;
157 fPattern->fUnprescaled &= 0xbf;
158 }
159
160 return kTRUE;
161}
Note: See TracBrowser for help on using the repository browser.