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);

#ifdef DEBUGPAMCUT
    cout << ">>> TofPatternCut CHECK <<<" << endl;
#endif

  bool TOF__OK = false;

#ifdef DEBUGPAMCUT
    cout << "Check topology:" << endl;
    for (int ip = 0; ip < 6; ip++)
      cout << "plane " << ip << " n.hit paddles " << tofl2->GetNHitPaddles(ip) << endl;
    for (int ip = 0; ip < 6; ip++)
      cout << "plane " << ip << " n.hit PMTs outside track " << GetNHitPMTsOutsideTrack(event, ip) << endl;
#endif

  // -----------------------------------------
  // at least one hit paddle on both S1 and S2
  // -----------------------------------------
  if (
  //-----------------------------------
  // not more than 1 paddle per layer
  //-----------------------------------
  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
      //-----------------------------------
      // at least 1 paddle per plane
      //-----------------------------------
      (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
      //-----------------------------------
      // n.hit pmts outside track
      //-----------------------------------
      //        GetHitPMTsOutsideTrack(event,0)<2 && //S11       //troppo (puo` succedere )
      //        GetHitPMTsOutsideTrack(event,1)<2 && //S12
      //        GetHitPMTsOutsideTrack(event,2)<2 && //S21
      //        GetHitPMTsOutsideTrack(event,3)<2 && //S22
      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) {
#ifdef DEBUGPAMCUT
      cout << "(((DISCARDED)))" << endl;
#endif
    return (0);
  };

#ifdef DEBUGPAMCUT
    cout << "--> OK" << endl;
    cout << "Check track concistency:" << endl;
#endif

  // -----------------------------------------------
  // the hit paddles MUST be associated to the track
  // (in collaborazione con wolfgang)
  // -----------------------------------------------
  //
  //  first build the map of hit PMTs along the track
  //  and discard the event if there is an 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
#ifdef DEBUGPAMCUT
      if (hitpaddle[iplane] >= 0 && hitpaddle[iplane] != ipaddle)
        cout << "ORRORE!!!!" << endl;
#endif
      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) {
#ifdef DEBUGPAMCUT
          cout << "hit plane " << iplane << " -- track PMT " << hitpaddle[iplane] << " is not hit " << endl;
          cout << "(((DISCARDED)))" << endl;
#endif
        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) {
#ifdef DEBUGPAMCUT
        cout << "hit plane " << iplane << " has no PMTs associated to the track " << endl;
        cout << "(((DISCARDED)))" << endl;
#endif
      return (0);//discard the event
    }
  }

  //
  // forse non c'e` bisogno di tutto cio`?
  //
  //     for(int iplane=0; iplane<4; iplane++){//loop over S11/S12/S21/S22
  //    //retrieve the id of the paddle traversed by the track
  //    // (...there is some tolerance...)
  //    int ipaddle = tofl2->GetPaddleIdOfTrack(tof->xtr_tof[iplane],tof->ytr_tof[iplane],iplane);
  //    //check if the traversed paddle is hit
  //    bool OK = true;
  //    if(
  //        tofl2->GetNHitPaddles(iplane)>0 &&//if there is a hit paddle in this plane...
  //        !tofl2->HitPaddle(iplane,ipaddle) && //...and the paddle traversed by the track is not hit
  //        true)OK = false;//..discard the event
  //    /// hence try to recover some events...
  //    if(
  //        tofl2->GetNHitPaddles(iplane)>0 &&//if there is a hit paddle...
  //        ipaddle==-1 &&//...and the track does not traverse any paddle...
  //        tofl2->HitPaddle(iplane,hitpaddle[iplane]) &&//... BUT there are tdc signals belonging to a hit paddle
  //        true)OK=true;//recover

  //    if(!OK)return (0);

  //     }

  //    if( !TOF__OK )return (0);

#ifdef DEBUGPAMCUT
    cout << "--> OK" << endl;
#endif

  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			#ifdef DEBUGPAMCUT
29			cout << ">>> TofPatternCut CHECK <<<" << endl;
30			#endif
31
32			bool TOF__OK = false;
33
34			#ifdef DEBUGPAMCUT
35			cout << "Check topology:" << endl;
36			for (int ip = 0; ip < 6; ip++)
37			cout << "plane " << ip << " n.hit paddles " << tofl2->GetNHitPaddles(ip) << endl;
38			for (int ip = 0; ip < 6; ip++)
39			cout << "plane " << ip << " n.hit PMTs outside track " << GetNHitPMTsOutsideTrack(event, ip) << endl;
40			#endif
41
42			// -----------------------------------------
43			// at least one hit paddle on both S1 and S2
44			// -----------------------------------------
45			if (
46			//-----------------------------------
47			// not more than 1 paddle per layer
48			//-----------------------------------
49			tofl2 && tofl2->GetNHitPaddles(0) <= 1 && //S11
50			tofl2->GetNHitPaddles(1) <= 1 && //S12
51			tofl2->GetNHitPaddles(2) <= 1 && //S21
52			tofl2->GetNHitPaddles(3) <= 1 && //S22
53			// tofl2->GetNHitPaddles(4) <= 1 && //S31
54			// tofl2->GetNHitPaddles(5) <= 1 && //S32
55			//-----------------------------------
56			// at least 1 paddle per plane
57			//-----------------------------------
58			(tofl2->GetNHitPaddles(0) == 1 \|\| tofl2->GetNHitPaddles(1) == 1) && //S1
59			(tofl2->GetNHitPaddles(2) == 1 \|\| tofl2->GetNHitPaddles(3) == 1) && //S2
60			// (tofl2->GetNHitPaddles(4)==1 \|\| tofl2->GetNHitPaddles(5)==1) && //S3
61			//-----------------------------------
62			// n.hit pmts outside track
63			//-----------------------------------
64			// GetHitPMTsOutsideTrack(event,0)<2 && //S11 //troppo (puo` succedere )
65			// GetHitPMTsOutsideTrack(event,1)<2 && //S12
66			// GetHitPMTsOutsideTrack(event,2)<2 && //S21
67			// GetHitPMTsOutsideTrack(event,3)<2 && //S22
68			GetNHitPMTsOutsideTrack(event, 0) < 3 && //S11
69			GetNHitPMTsOutsideTrack(event, 1) < 3 && //S12
70			// GetHitPMTsOutsideTrack(event,2)<3 && //S21
71			// GetHitPMTsOutsideTrack(event,3)<3 && //S22
72			true)
73			TOF__OK = true;
74
75			if (!TOF__OK) {
76			#ifdef DEBUGPAMCUT
77			cout << "(((DISCARDED)))" << endl;
78			#endif
79			return (0);
80			};
81
82			#ifdef DEBUGPAMCUT
83			cout << "--> OK" << endl;
84			cout << "Check track concistency:" << endl;
85			#endif
86
87			// -----------------------------------------------
88			// the hit paddles MUST be associated to the track
89			// (in collaborazione con wolfgang)
90			// -----------------------------------------------
91			//
92			// first build the map of hit PMTs along the track
93			// and discard the event if there is an hit paddle,
94			// but it is not the one associated to the track
95			//
96
97			int hitplane[] = { 0, 0, 0, 0, 0, 0 };
98			int hitpaddle[] = { -1, -1, -1, -1, -1, -1 };
99			for (Int_t ip = 0; ip < tof->npmttdc; ip++) {//loop over tdc hits associated to the track
100			Int_t iplane = -1;
101			Int_t ipaddle = -1;
102			Int_t ipmt = -1;
103			Int_t id = tof->pmttdc[ip];//pmt id
104			tofl2->GetPMTName(id, iplane, ipaddle, ipmt);
105			if (tof->tdcflag[ip] == 0) {
106			hitplane[iplane] = 1; //there is a true tdc signal associated to the track
107			#ifdef DEBUGPAMCUT
108			if (hitpaddle[iplane] >= 0 && hitpaddle[iplane] != ipaddle)
109			cout << "ORRORE!!!!" << endl;
110			#endif
111			hitpaddle[iplane] = ipaddle;//store the id of the paddle associated to the track
112			//
113			//
114			//
115			if (tofl2->GetNHitPaddles(iplane) > 0 && //if there is a hit paddle in this plane...
116			!tofl2->HitPaddle(iplane, hitpaddle[iplane]) && //... and the paddle associated to the track is not hit
117			true) {
118			#ifdef DEBUGPAMCUT
119			cout << "hit plane " << iplane << " -- track PMT " << hitpaddle[iplane] << " is not hit " << endl;
120			cout << "(((DISCARDED)))" << endl;
121			#endif
122			return (0);//discard the event
123			}
124			}
125			};
126			//
127			// hence check the planes with no PMTs associated to the track:
128			// the event is discarded if there is a hit paddle
129			//
130			for (int iplane = 0; iplane < 4; iplane++) {//loop over S11/S12/S21/S22
131			if (hitplane[iplane] == 0 && //if there are no PMTs associated to the track...
132			tofl2->GetNHitPaddles(iplane) > 0 && //...and the plane has hit paddles
133			true) {
134			#ifdef DEBUGPAMCUT
135			cout << "hit plane " << iplane << " has no PMTs associated to the track " << endl;
136			cout << "(((DISCARDED)))" << endl;
137			#endif
138			return (0);//discard the event
139			}
140			}
141
142			//
143			// forse non c'e` bisogno di tutto cio`?
144			//
145			// for(int iplane=0; iplane<4; iplane++){//loop over S11/S12/S21/S22
146			// //retrieve the id of the paddle traversed by the track
147			// // (...there is some tolerance...)
148			// int ipaddle = tofl2->GetPaddleIdOfTrack(tof->xtr_tof[iplane],tof->ytr_tof[iplane],iplane);
149			// //check if the traversed paddle is hit
150			// bool OK = true;
151			// if(
152			// tofl2->GetNHitPaddles(iplane)>0 &&//if there is a hit paddle in this plane...
153			// !tofl2->HitPaddle(iplane,ipaddle) && //...and the paddle traversed by the track is not hit
154			// true)OK = false;//..discard the event
155			// /// hence try to recover some events...
156			// if(
157			// tofl2->GetNHitPaddles(iplane)>0 &&//if there is a hit paddle...
158			// ipaddle==-1 &&//...and the track does not traverse any paddle...
159			// tofl2->HitPaddle(iplane,hitpaddle[iplane]) &&//... BUT there are tdc signals belonging to a hit paddle
160			// true)OK=true;//recover
161
162			// if(!OK)return (0);
163
164			// }
165
166			// if( !TOF__OK )return (0);
167
168			#ifdef DEBUGPAMCUT
169			cout << "--> OK" << endl;
170			#endif
171
172			return CUTOK; //good pattern
173
174			}
175			//===========================================================
176			//
177			//
178			//
179			//
180			//
181			//===========================================================
182			int TofPatternCut::GetNHitPMTsOutsideTrack(PamLevel2 *event, int plane) {
183
184			if (!event)
185			return (0);
186			//
187			ToFLevel2 *tofl2 = event->GetToFLevel2();
188			if (!tofl2)
189			return (0);
190			ToFTrkVar *tof = 0;
191			if (!_standalone && _notrk >= 0 && _notrk < event->GetTrkLevel2()->GetNTracks())
192			tof = event->GetTrack(_notrk)->GetToFTrack();
193			if (_standalone)
194			tof = event->GetToFStoredTrack(-1);
195			if (!tof)
196			return (0);
197
198			int nn = 0;
199			// cout << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"<<endl;
200			for (int ip = 0; ip < tofl2->npmt(); ip++) { //loop over ALL PMTs
201			Int_t iplane = -1;
202			Int_t ipaddle = -1;
203			Int_t ipmt = -1;
204			Int_t id = tofl2->GetToFPMT(ip)->pmt_id; //PMT id
205			tofl2->GetPMTName(id, iplane, ipaddle, ipmt);
206
207			if (iplane != plane)
208			continue;
209			// cout << endl<<plane << " -- "<<id;
210			float tdc = tofl2->GetToFPMT(ip)->tdc;
211			if (tdc >= 4095.)
212			continue;///to avoid adc pile-up
213			// spiegazione:
214			// quando c'e` adc pile-up, succede che un pmt abbia un segnale adc (adc<4095) ma NON tdc (tdc==4095),
215			// che e` il contrario di quanto avviene di solito. quindi, se nel livello 2 c'e` un pmt colpito
216			// e questo ha tdc==4095, implica che c'e` solo il segnale adc e che probabilmente non e` un vero segnale,
217			// ma e` il pile-up dell'evento precedente
218
219			// loop over TDC hits along the track
220			int iht = -1;
221			for (iht = 0; iht < tof->npmttdc; iht++)
222			if (id == tof->pmttdc[iht])
223			break; //match
224			if (iht >= 0 && iht < tof->npmttdc)
225			continue;
226			// loop over ADC hits along the track (credo sia inutile...)
227			// int iha = -1;
228			// for(iha=0; iha<tof->npmtadc; iha++)if(id==tof->pmtadc[iha])break; //match
229			// if( iha>=0 && iha <tof->npmtadc )continue;
230			// get here if outside track
231			nn++;
232
233			// cout << " OUTSIDE ";
234			}
235			return nn;
236
237			}
238