source: trunk/MagicSoft/Mars/mtrigger/MTriggerPatternDecode.cc@ 9488

Last change on this file since 9488 was 9480, checked in by tbretz, 15 years ago
*** empty log message ***
File size: 5.2 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 to 00000000.
32// This can be changed using the information about 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->FindCreateObj("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 {
89 *fLog << err << "MRawEvtHeader not found... abort." << endl;
90 return kFALSE;
91 }
92
93 return kTRUE;
94}
95
96// --------------------------------------------------------------------------
97//
98// For files before file version 5 the trigger pattern is set to 00000000.
99// This can be changed using the information about the file-type.
100//
101Int_t MTriggerPatternDecode::Process()
102{
103 // No setting necessary because MTriggerPattern::reset() has
104 // been called just before
105 if (fRunHeader->GetFormatVersion()<5)
106 return kTRUE;
107
108 const UInt_t pattern = ~fEvtHeader->GetTriggerID();
109
110 // The the trigger pattern is currently written with inverted bits,
111 // but in the future this could be changed. In this case the file version
112 // number of the raw-data format must be checked. It tells you the
113 // file format (meaning of the bits)
114 // If for some reason the file-format has not been correctly changed
115 // use the run-number to decide whether the bits must be inverted or not.
116
117 // Trigger Pattern is a number of 16 bits, but if the machine
118 // is a 64 bits then the bit inversion the first 16 bits are set to 1,
119 // possibly giving problems with the filtering with the masks.
120 // Because UInt_t by definition is a 32 bit number on any
121 // machine (see root manual) no check is needed.
122 // The simplest workaround is:
123 // pattern &= 0xffffffff;
124
125 // For the moment the meaning of the bits in
126 // MRawEvtHeader::GetTriggerID and in MTriggerPattern::fTriggerPattern
127 // are identical - this may change in the future! In this case
128 // the decoding HERE must be changed - NOT MTriggerPattern.
129 // If an enhancement of the information stored in MTriggerPattern
130 // is necessary ADD new variables! Don't change the meaning of the old
131 // ones!
132 fPattern->fPrescaled = pattern & 0xff;
133 fPattern->fUnprescaled = (pattern>>8) & 0xff;
134
135 cout << (int)((pattern>>6)&1) << " " << fPattern->fPrescaled << endl;
136
137 fPattern->Print();
138
139 // This is a workaround for the new scheme in which L1TPU (the signal
140 // comming directly from the L1 is connected, but the L1 (routed
141 // over L2 is disconnected)
142 if (!fRunHeader->IsMonteCarloRun() && fRunHeader->GetTelescopeNumber()==1 &&
143 fRunHeader->GetRunNumber()>1006246)
144 {
145 fPattern->fPrescaled |= (pattern>> 6)&1;
146 fPattern->fUnprescaled |= (pattern>>14)&1;
147
148 fPattern->fPrescaled &= 0xbf;
149 fPattern->fUnprescaled &= 0xbf;
150 }
151
152 fPattern->Print();
153
154 return kTRUE;
155}
Note: See TracBrowser for help on using the repository browser.