PamOffLineSW/physics/TriggerReader.cpp

/** @file
 * $Source: /repository/PamOffLineSW/physics/TriggerReader.cpp,v $
 * $Id: TriggerReader.cpp,v 1.4 2008-03-05 11:37:37 messineo Exp $
 * $Author: messineo $
 * 
 * Implementation of the CalorimeterReader class.
 */

#include <iostream>
#include <string>
#include "TriggerReader.h"

extern "C" {
  

  extern struct {
    int         evcount;
    int         pmtpl[3];
    int         trigrate[6];
    int         dltime[2];
    int         s4calcount[2]; 
    int         pmtcount1[24];
    int         pmtcount2[24];
    int         patternbusy[3];
    int         patterntrig[6];
    int         trigconf;
  } trig_; 

  void triggerunpack_(unsigned char[], long int*, int*);

    //Struct per il passaggio di dati da e verso la chiamata fortran
}

using namespace pamela;
using namespace pamela::trigger;

 
/**
 * Constructor.     
 */
TriggerReader::TriggerReader(void): 
  TechmodelAlgorithm(PacketType::Physics, "TechmodelTriggerReader") { 
  trigger = new TriggerEvent();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string TriggerReader::GetVersionInfo(void) const {
  return 
    "$Header: /repository/PamOffLineSW/physics/TriggerReader.cpp,v 1.4 2008-03-05 11:37:37 messineo Exp $";
}

/**
 * Initialize the algorithm with a special run. This will initialize the
 * event reader routines for all packet types.
 */
void TriggerReader::Init(PamelaRun *run) {
  run->WriteSubPacket(this, &trigger, trigger->Class());
}

/**
 * Unpack the trigger event from an input file.
 */
void TriggerReader::RunEvent(int EventNumber) {
  
}

/**
 * Unpack the Trigger data event from the physical packet.
 */
void TriggerReader::RunEvent(int EventNumber, const char subData[], long int length) {
        std::stringstream oss; 
    string msg;

    char *data = new char[length];
    memcpy(data, subData, length);
    int  ERROR = 0;
    
    /*unsigned short convdata[length];
    for (int i = 0; i<length; i++){
        convdata[i] = (unsigned short)((unsigned char)subData[i]&0xFF);
    }*/

    
    //Call to the routine that unpack trigger events
    if ( length < 500000 && length >= 0 ){
      triggerunpack_((unsigned char*)data, &length, &ERROR);
    } else {
      ERROR = 1;
    };
    

    if (ERROR != 0) {
        char *errmsg;
        switch (ERROR){
            case 1: errmsg = "GENERIC TRIGGER ERROR";
            break;
            default: errmsg = "TRIGGER ERRROR CODE UNIDENTIFIED";
        }
        oss.str("");
        oss << "Trigger: Fortran77 function triggerunpack: " <<  errmsg;
        msg=oss.str();
        PamOffLineSW::mainLogUtil->logWarning(msg);       
        
    }
    //  In case of "ERROR != 0" the calunpack will take care to set all 
    //  parameters to zero
    //} else {
    trigger->unpackError = (unsigned short)ERROR;
    trigger->evcount   = trig_.evcount;
    memcpy(trigger->pmtpl, trig_.pmtpl, sizeof(trigger->pmtpl));
    memcpy(trigger->trigrate, trig_.trigrate, sizeof(trigger->trigrate));
    memcpy(trigger->dltime, trig_.dltime, sizeof(trigger->dltime));
    memcpy(trigger->s4calcount, trig_.s4calcount, sizeof(trigger->s4calcount));
    memcpy(trigger->pmtcount1, trig_.pmtcount1, sizeof(trigger->pmtcount1));
    memcpy(trigger->pmtcount2, trig_.pmtcount2, sizeof(trigger->pmtcount2));
    memcpy(trigger->patternbusy, trig_.patternbusy, sizeof(trigger->patternbusy));
    memcpy(trigger->patterntrig, trig_.patterntrig, sizeof(trigger->patterntrig));
    trigger->trigconf = trig_.trigconf;
    delete [] data;

}
1	/** @file
2	* $Source: /repository/PamOffLineSW/physics/TriggerReader.cpp,v $
3	* $Id: TriggerReader.cpp,v 1.4 2008-03-05 11:37:37 messineo Exp $
4	* $Author: messineo $
5	*
6	* Implementation of the CalorimeterReader class.
7	*/
8
9	#include <iostream>
10	#include <string>
11	#include "TriggerReader.h"
12
13	extern "C" {
14
15
16	extern struct {
17	int evcount;
18	int pmtpl[3];
19	int trigrate[6];
20	int dltime[2];
21	int s4calcount[2];
22	int pmtcount1[24];
23	int pmtcount2[24];
24	int patternbusy[3];
25	int patterntrig[6];
26	int trigconf;
27	} trig_;
28
29	void triggerunpack_(unsigned char[], long int, int);
30
31	//Struct per il passaggio di dati da e verso la chiamata fortran
32	}
33
34	using namespace pamela;
35	using namespace pamela::trigger;
36
37
38	/**
39	* Constructor.
40	*/
41	TriggerReader::TriggerReader(void):
42	TechmodelAlgorithm(PacketType::Physics, "TechmodelTriggerReader") {
43	trigger = new TriggerEvent();
44	}
45
46	/**
47	* Get a string with the version info of the algorithm.
48	*/
49	std::string TriggerReader::GetVersionInfo(void) const {
50	return
51	"$Header: /repository/PamOffLineSW/physics/TriggerReader.cpp,v 1.4 2008-03-05 11:37:37 messineo Exp $";
52	}
53
54	/**
55	* Initialize the algorithm with a special run. This will initialize the
56	* event reader routines for all packet types.
57	*/
58	void TriggerReader::Init(PamelaRun *run) {
59	run->WriteSubPacket(this, &trigger, trigger->Class());
60	}
61
62	/**
63	* Unpack the trigger event from an input file.
64	*/
65	void TriggerReader::RunEvent(int EventNumber) {
66
67	}
68
69	/**
70	* Unpack the Trigger data event from the physical packet.
71	*/
72	void TriggerReader::RunEvent(int EventNumber, const char subData[], long int length) {
73	std::stringstream oss;
74	string msg;
75
76	char *data = new char[length];
77	memcpy(data, subData, length);
78	int ERROR = 0;
79
80	/*unsigned short convdata[length];
81	for (int i = 0; i<length; i++){
82	convdata[i] = (unsigned short)((unsigned char)subData[i]&0xFF);
83	}*/
84
85
86	//Call to the routine that unpack trigger events
87	if ( length < 500000 && length >= 0 ){
88	triggerunpack_((unsigned char*)data, &length, &ERROR);
89	} else {
90	ERROR = 1;
91	};
92
93
94	if (ERROR != 0) {
95	char *errmsg;
96	switch (ERROR){
97	case 1: errmsg = "GENERIC TRIGGER ERROR";
98	break;
99	default: errmsg = "TRIGGER ERRROR CODE UNIDENTIFIED";
100	}
101	oss.str("");
102	oss << "Trigger: Fortran77 function triggerunpack: " << errmsg;
103	msg=oss.str();
104	PamOffLineSW::mainLogUtil->logWarning(msg);
105
106	}
107	// In case of "ERROR != 0" the calunpack will take care to set all
108	// parameters to zero
109	//} else {
110	trigger->unpackError = (unsigned short)ERROR;
111	trigger->evcount = trig_.evcount;
112	memcpy(trigger->pmtpl, trig_.pmtpl, sizeof(trigger->pmtpl));
113	memcpy(trigger->trigrate, trig_.trigrate, sizeof(trigger->trigrate));
114	memcpy(trigger->dltime, trig_.dltime, sizeof(trigger->dltime));
115	memcpy(trigger->s4calcount, trig_.s4calcount, sizeof(trigger->s4calcount));
116	memcpy(trigger->pmtcount1, trig_.pmtcount1, sizeof(trigger->pmtcount1));
117	memcpy(trigger->pmtcount2, trig_.pmtcount2, sizeof(trigger->pmtcount2));
118	memcpy(trigger->patternbusy, trig_.patternbusy, sizeof(trigger->patternbusy));
119	memcpy(trigger->patterntrig, trig_.patterntrig, sizeof(trigger->patterntrig));
120	trigger->trigconf = trig_.trigconf;
121	delete [] data;
122
123	}