123 recHodo = pHodo.product();
134 recTDC = pTDC.product();
145 evtHeader = pEventHeader.product();
147 std::cout <<
"Error! can't get the event header " <<std::endl;
148 std::cerr <<
"Error! can't get the event header " <<std::endl;
153 if (!recTDC) {
noTdc++;
return; }
154 if (!recHodo) {
noHodo++;
return; }
155 if (!evtHeader) {
noHeader++;
return; }
164 int event = evtHeader->eventNumber();
168 int nomXtalInBeam = -999;
169 int nextXtalInBeam = -999;
172 if (xtalInBeamId==
EBDetId(0)){
return; }
173 int mySupCry = xtalInBeamId.ic();
174 int mySupEta = xtalInBeamId.ieta();
175 int mySupPhi = xtalInBeamId.iphi();
180 double x = recHodo->posX();
181 double y = recHodo->posY();
182 double sx = recHodo->slopeX();
183 double sy = recHodo->slopeY();
184 double qx = recHodo->qualX();
185 double qy = recHodo->qualY();
190 double tdcOffset = recTDC->offset();
199 for (
unsigned int icry=0; icry<49; icry++){
200 unsigned int row = icry/7;
201 unsigned int column = icry%7;
206 if ( Xtals7x7[icry].
ism() == 1){
208 crystal[icry] = Xtals7x7[icry].ic();
210 if ( Xtals7x7[icry].
ism() != 1){
212 crystal[icry] = -100;
220 crystal[icry] = -100;
228 double maxEne = -999.;
229 int maxEneCry = 9999;
230 int maxEneInMatrix = -999;
231 for (
int ii=0; ii<49; ii++){
if (
energy[ii] > maxEne){
233 maxEneCry = crystal[ii];
234 maxEneInMatrix = ii;}
244 myTree->
fillInfo(
run,
event, mySupCry, maxEneCry, nomXtalInBeam, nextXtalInBeam, mySupEta, mySupPhi, tbm,
x,
y, Xcal, Ycal, sx, sy, qx, qy, tdcOffset, allMatrix,
energy, crystal);
std::string tdcRecInfoProducer_
std::string eventHeaderProducer_
std::string hodoRecInfoProducer_
std::string tdcRecInfoCollection_
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
static const int ETAPHIMODE
void fillInfo(int run, int eve, int xnum, int maxX, int nomX, int nextX, int xeta, int xphi, int tbm, double xx, double yy, double ecalx, double ecaly, double sx, double sy, double qx, double qy, double tdcoff, int allm, double amp[], int cry[])
static const int SMCRYSTALMODE
std::string RecHitProducer_
int ism(int ieta, int iphi)
std::string EBRecHitCollection_
std::string hodoRecInfoCollection_
Public Member Functions inherited from edm::EDAnalyzer
1.8.5
