TofCuts/TofPatternCut/TofPatternCut.cpp

/*
 * TofPatternCut.cpp
 *
 *  Created on: 12-mar-2009
 *      Author: Elena Vannuccini
 */

/*! @file TofPatternCut.cpp The ToFPatternCut class implementation file */

#include "TofPatternCut.h"

int TofPatternCut::Check(PamLevel2 *event) {

  if (!event)
    return (0);
  //
  ToFLevel2 *tofl2 = event->GetToFLevel2();
  if (!tofl2)
    return (0);
  ToFTrkVar *tof = 0;
  if (!_standalone && _notrk >= 0 && _notrk < event->GetTrkLevel2()->GetNTracks())
    tof = event->GetTrack(_notrk)->GetToFTrack();
  if (_standalone)
    tof = event->GetToFStoredTrack(-1);
  if (!tof)
    return (0);

  bool TOF__OK = false;

  if (
      //-----------------------------------
      // (1) no more than 1 paddle on S11, S12, S21, S22
      //-----------------------------------
      tofl2 && tofl2->GetNHitPaddles(0) <= 1 && //S11
      tofl2->GetNHitPaddles(1) <= 1 && //S12
      tofl2->GetNHitPaddles(2) <= 1 && //S21
      tofl2->GetNHitPaddles(3) <= 1 && //S22
      //        tofl2->GetNHitPaddles(4) <= 1 && //S31
      //        tofl2->GetNHitPaddles(5) <= 1 && //S32
      //-----------------------------------
      // (2) at least 1 hit paddle on S1, S2
      //-----------------------------------
      (tofl2->GetNHitPaddles(0) == 1 || tofl2->GetNHitPaddles(1) == 1) && //S1
      (tofl2->GetNHitPaddles(2) == 1 || tofl2->GetNHitPaddles(3) == 1) && //S2
      //        (tofl2->GetNHitPaddles(4)==1 || tofl2->GetNHitPaddles(5)==1) && //S3
      //-----------------------------------
      // (3) less than 3 hit pmts outside track on S11,S12
      //-----------------------------------
      GetNHitPMTsOutsideTrack(event, 0) < 3 && //S11
      GetNHitPMTsOutsideTrack(event, 1) < 3 && //S12
      //        GetHitPMTsOutsideTrack(event,2)<3 && //S21
      //        GetHitPMTsOutsideTrack(event,3)<3 && //S22
      true)
    TOF__OK = true;

  if (!TOF__OK) {
    return (0);
  };

  // -----------------------------------------------
  // (4) for S11,S12,S21,S22: if there is a hit paddle, it MUST be associated to the track (done with Wolfgang help)
  // -----------------------------------------------
  //
  //  first build the map of hit PMTs along the track
  //  and discard the event if there is a hit paddle,
  //  but it is not the one associated to the track
  //
  int hitplane[] = { 0, 0, 0, 0, 0, 0 };
  int hitpaddle[] = { -1, -1, -1, -1, -1, -1 };
  for (Int_t ip = 0; ip < tof->npmttdc; ip++) {//loop over tdc hits associated to the track
    Int_t iplane = -1;
    Int_t ipaddle = -1;
    Int_t ipmt = -1;
    Int_t id = tof->pmttdc[ip];//pmt id
    tofl2->GetPMTName(id, iplane, ipaddle, ipmt);
    if (tof->tdcflag[ip] == 0) {
      hitplane[iplane] = 1; //there is a true tdc signal associated to the track
      hitpaddle[iplane] = ipaddle;//store the id of the paddle associated to the track
      //
      //
      //
      if (tofl2->GetNHitPaddles(iplane) > 0 && //if there is a hit paddle in this plane...
          !tofl2->HitPaddle(iplane, hitpaddle[iplane]) && //... and the paddle associated to the track is not hit
          true) {
        return (0);//discard the event
      }
    }
  };
  //
  // hence check the planes with no PMTs associated to the track:
  // the event is discarded if there is a hit paddle
  //
  for (int iplane = 0; iplane < 4; iplane++) {//loop over S11/S12/S21/S22
    if (hitplane[iplane] == 0 && //if there are no PMTs associated to the track...
        tofl2->GetNHitPaddles(iplane) > 0 && //...and the plane has hit paddles
        true) {
      return (0);//discard the event
    }
  }

  return CUTOK; //good pattern

}
//===========================================================
//
//
//
//
//
//===========================================================
int TofPatternCut::GetNHitPMTsOutsideTrack(PamLevel2 *event, int plane) {

  if (!event)
    return (0);
  //
  ToFLevel2 *tofl2 = event->GetToFLevel2();
  if (!tofl2)
    return (0);
  ToFTrkVar *tof = 0;
  if (!_standalone && _notrk >= 0 && _notrk < event->GetTrkLevel2()->GetNTracks())
    tof = event->GetTrack(_notrk)->GetToFTrack();
  if (_standalone)
    tof = event->GetToFStoredTrack(-1);
  if (!tof)
    return (0);

  int nn = 0;
  //    cout << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"<<endl;
  for (int ip = 0; ip < tofl2->npmt(); ip++) { //loop over ALL PMTs
    Int_t iplane = -1;
    Int_t ipaddle = -1;
    Int_t ipmt = -1;
    Int_t id = tofl2->GetToFPMT(ip)->pmt_id; //PMT id
    tofl2->GetPMTName(id, iplane, ipaddle, ipmt);

    if (iplane != plane)
      continue;
    //  cout << endl<<plane << " -- "<<id;
    float tdc = tofl2->GetToFPMT(ip)->tdc;
    if (tdc >= 4095.)
      continue;///to avoid adc pile-up
    // spiegazione:
    // quando c'e` adc pile-up, succede che un pmt abbia un segnale adc (adc<4095) ma NON tdc (tdc==4095),
    // che e` il contrario di quanto avviene di solito. quindi, se nel livello 2 c'e` un pmt colpito
    // e questo ha tdc==4095, implica che c'e` solo il segnale adc e che probabilmente non e` un vero segnale,
    // ma e` il pile-up dell'evento precedente

    // loop over TDC hits along the track
    int iht = -1;
    for (iht = 0; iht < tof->npmttdc; iht++)
      if (id == tof->pmttdc[iht])
        break; //match
    if (iht >= 0 && iht < tof->npmttdc)
      continue;
    // loop over ADC hits along the track (credo sia inutile...)
    //  int iha = -1;
    //          for(iha=0; iha<tof->npmtadc; iha++)if(id==tof->pmtadc[iha])break; //match
    //          if( iha>=0 && iha <tof->npmtadc )continue;
    // get here if outside track
    nn++;

    //  cout << " OUTSIDE ";
  }
  return nn;

}

1	pam-fi	1.1	/*
2			* TofPatternCut.cpp
3			*
4			* Created on: 12-mar-2009
5			* Author: Elena Vannuccini
6			*/
7
8			/! @file TofPatternCut.cpp The ToFPatternCut class implementation file /
9
10			#include "TofPatternCut.h"
11
12			int TofPatternCut::Check(PamLevel2 *event) {
13
14			if (!event)
15			return (0);
16			//
17			ToFLevel2 *tofl2 = event->GetToFLevel2();
18			if (!tofl2)
19			return (0);
20			ToFTrkVar *tof = 0;
21			if (!_standalone && _notrk >= 0 && _notrk < event->GetTrkLevel2()->GetNTracks())
22			tof = event->GetTrack(_notrk)->GetToFTrack();
23			if (_standalone)
24			tof = event->GetToFStoredTrack(-1);
25			if (!tof)
26			return (0);
27
28			bool TOF__OK = false;
29
30			if (
31	pam-fi	1.2	//-----------------------------------
32			// (1) no more than 1 paddle on S11, S12, S21, S22
33			//-----------------------------------
34			tofl2 && tofl2->GetNHitPaddles(0) <= 1 && //S11
35	pam-fi	1.1	tofl2->GetNHitPaddles(1) <= 1 && //S12
36			tofl2->GetNHitPaddles(2) <= 1 && //S21
37			tofl2->GetNHitPaddles(3) <= 1 && //S22
38			// tofl2->GetNHitPaddles(4) <= 1 && //S31
39			// tofl2->GetNHitPaddles(5) <= 1 && //S32
40			//-----------------------------------
41	pam-fi	1.2	// (2) at least 1 hit paddle on S1, S2
42	pam-fi	1.1	//-----------------------------------
43			(tofl2->GetNHitPaddles(0) == 1 \|\| tofl2->GetNHitPaddles(1) == 1) && //S1
44			(tofl2->GetNHitPaddles(2) == 1 \|\| tofl2->GetNHitPaddles(3) == 1) && //S2
45			// (tofl2->GetNHitPaddles(4)==1 \|\| tofl2->GetNHitPaddles(5)==1) && //S3
46			//-----------------------------------
47	pam-fi	1.2	// (3) less than 3 hit pmts outside track on S11,S12
48	pam-fi	1.1	//-----------------------------------
49			GetNHitPMTsOutsideTrack(event, 0) < 3 && //S11
50			GetNHitPMTsOutsideTrack(event, 1) < 3 && //S12
51			// GetHitPMTsOutsideTrack(event,2)<3 && //S21
52			// GetHitPMTsOutsideTrack(event,3)<3 && //S22
53			true)
54			TOF__OK = true;
55
56			if (!TOF__OK) {
57			return (0);
58			};
59
60			// -----------------------------------------------
61	pam-fi	1.2	// (4) for S11,S12,S21,S22: if there is a hit paddle, it MUST be associated to the track (done with Wolfgang help)
62	pam-fi	1.1	// -----------------------------------------------
63			//
64			// first build the map of hit PMTs along the track
65	pam-fi	1.2	// and discard the event if there is a hit paddle,
66	pam-fi	1.1	// but it is not the one associated to the track
67			//
68			int hitplane[] = { 0, 0, 0, 0, 0, 0 };
69			int hitpaddle[] = { -1, -1, -1, -1, -1, -1 };
70			for (Int_t ip = 0; ip < tof->npmttdc; ip++) {//loop over tdc hits associated to the track
71			Int_t iplane = -1;
72			Int_t ipaddle = -1;
73			Int_t ipmt = -1;
74			Int_t id = tof->pmttdc[ip];//pmt id
75			tofl2->GetPMTName(id, iplane, ipaddle, ipmt);
76			if (tof->tdcflag[ip] == 0) {
77			hitplane[iplane] = 1; //there is a true tdc signal associated to the track
78			hitpaddle[iplane] = ipaddle;//store the id of the paddle associated to the track
79			//
80			//
81			//
82			if (tofl2->GetNHitPaddles(iplane) > 0 && //if there is a hit paddle in this plane...
83			!tofl2->HitPaddle(iplane, hitpaddle[iplane]) && //... and the paddle associated to the track is not hit
84			true) {
85			return (0);//discard the event
86			}
87			}
88			};
89			//
90			// hence check the planes with no PMTs associated to the track:
91			// the event is discarded if there is a hit paddle
92			//
93			for (int iplane = 0; iplane < 4; iplane++) {//loop over S11/S12/S21/S22
94			if (hitplane[iplane] == 0 && //if there are no PMTs associated to the track...
95			tofl2->GetNHitPaddles(iplane) > 0 && //...and the plane has hit paddles
96			true) {
97			return (0);//discard the event
98			}
99			}
100
101			return CUTOK; //good pattern
102
103			}
104			//===========================================================
105			//
106			//
107			//
108			//
109			//
110			//===========================================================
111			int TofPatternCut::GetNHitPMTsOutsideTrack(PamLevel2 *event, int plane) {
112
113			if (!event)
114			return (0);
115			//
116			ToFLevel2 *tofl2 = event->GetToFLevel2();
117			if (!tofl2)
118			return (0);
119			ToFTrkVar *tof = 0;
120			if (!_standalone && _notrk >= 0 && _notrk < event->GetTrkLevel2()->GetNTracks())
121			tof = event->GetTrack(_notrk)->GetToFTrack();
122			if (_standalone)
123			tof = event->GetToFStoredTrack(-1);
124			if (!tof)
125			return (0);
126
127			int nn = 0;
128			// cout << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"<<endl;
129			for (int ip = 0; ip < tofl2->npmt(); ip++) { //loop over ALL PMTs
130			Int_t iplane = -1;
131			Int_t ipaddle = -1;
132			Int_t ipmt = -1;
133			Int_t id = tofl2->GetToFPMT(ip)->pmt_id; //PMT id
134			tofl2->GetPMTName(id, iplane, ipaddle, ipmt);
135
136			if (iplane != plane)
137			continue;
138			// cout << endl<<plane << " -- "<<id;
139			float tdc = tofl2->GetToFPMT(ip)->tdc;
140			if (tdc >= 4095.)
141			continue;///to avoid adc pile-up
142			// spiegazione:
143			// quando c'e` adc pile-up, succede che un pmt abbia un segnale adc (adc<4095) ma NON tdc (tdc==4095),
144			// che e` il contrario di quanto avviene di solito. quindi, se nel livello 2 c'e` un pmt colpito
145			// e questo ha tdc==4095, implica che c'e` solo il segnale adc e che probabilmente non e` un vero segnale,
146			// ma e` il pile-up dell'evento precedente
147
148			// loop over TDC hits along the track
149			int iht = -1;
150			for (iht = 0; iht < tof->npmttdc; iht++)
151			if (id == tof->pmttdc[iht])
152			break; //match
153			if (iht >= 0 && iht < tof->npmttdc)
154			continue;
155			// loop over ADC hits along the track (credo sia inutile...)
156			// int iha = -1;
157			// for(iha=0; iha<tof->npmtadc; iha++)if(id==tof->pmtadc[iha])break; //match
158			// if( iha>=0 && iha <tof->npmtadc )continue;
159			// get here if outside track
160			nn++;
161
162			// cout << " OUTSIDE ";
163			}
164			return nn;
165
166			}
167