TofCuts/TofDedxS1Cut/TofDedxS1Cut.cpp

/*
 * TofDedxS1Cut.cpp
 *
 *  Created on: 13-mar-2009
 *      Author: Sergio Ricciarini, Nicola Mori
 *  Reworked on: 21-jul-2011
 *      Author: Nicola Mori, Valeria Di Felice
 */

/*! @file TofDedxS1Cut.cpp The TofDedxS1Cut class implementation file */

#include "TofDedxS1Cut.h"

int TofDedxS1Cut::Check(PamLevel2 *event) {

  int trkSeqNo = 0;
  for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) {
    if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == event->GetTrack(0)->GetTrkTrack()->GetSeqNo())
      trkSeqNo = i;
  }

  float dEdx = 0;
  int nLayers = 0;
  float dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 0, 100);

  // Patch to exclude S115b
  ToFTrkVar *tofTrk = event->GetToFLevel2()->GetToFTrkVar(trkSeqNo);
  bool has8 = false, has9 = false;
  float dedx9 = 0.;
  for (Int_t ib = 0; ib < tofTrk->npmtadc; ib++) {
    if (tofTrk->pmtadc[ib] == 8) {
      has8 = true;
    }
    if (tofTrk->pmtadc[ib] == 9) {
      has9 = true;
      dedx9 = tofTrk->dedx[ib];
    }
  }

  if (has9){
    dEdxLayer -= (dedx9 / ((int) has8 + 1));
    dEdxLayer *= 2.;
  }
  //End of patch

  if (dEdxLayer > 0. && dEdxLayer < 4090) {
    dEdx += dEdxLayer;
    nLayers++;
  }
  dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 1, 100);
  if (dEdxLayer > 0. && dEdxLayer < 4090) {
    dEdx += dEdxLayer;
    nLayers++;
  }

  if (nLayers > 0) {
    dEdx /= nLayers;

    float R = event->GetTrack(0)->GetTrkTrack()->GetRigidity();

    if (dEdx >= (3.102 - 2.6 * TMath::Log(R)) * (R < 1.05) + 3. * (R >= 1.05))
      return 0;

    return CUTOK;
//    if (dEdx < (3.102 - 2.6 * TMath::Log(R)) * (R < 1.05) + 3. * (R >= 1.05))
//      return 0;
//
//    return CUTOK;
  }
  else {
    return 0;
  }

}
1	/*
2	* TofDedxS1Cut.cpp
3	*
4	* Created on: 13-mar-2009
5	* Author: Sergio Ricciarini, Nicola Mori
6	* Reworked on: 21-jul-2011
7	* Author: Nicola Mori, Valeria Di Felice
8	*/
9
10	/! @file TofDedxS1Cut.cpp The TofDedxS1Cut class implementation file /
11
12	#include "TofDedxS1Cut.h"
13
14	int TofDedxS1Cut::Check(PamLevel2 *event) {
15
16	int trkSeqNo = 0;
17	for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) {
18	if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == event->GetTrack(0)->GetTrkTrack()->GetSeqNo())
19	trkSeqNo = i;
20	}
21
22	float dEdx = 0;
23	int nLayers = 0;
24	float dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 0, 100);
25
26	// Patch to exclude S115b
27	ToFTrkVar *tofTrk = event->GetToFLevel2()->GetToFTrkVar(trkSeqNo);
28	bool has8 = false, has9 = false;
29	float dedx9 = 0.;
30	for (Int_t ib = 0; ib < tofTrk->npmtadc; ib++) {
31	if (tofTrk->pmtadc[ib] == 8) {
32	has8 = true;
33	}
34	if (tofTrk->pmtadc[ib] == 9) {
35	has9 = true;
36	dedx9 = tofTrk->dedx[ib];
37	}
38	}
39
40	if (has9){
41	dEdxLayer -= (dedx9 / ((int) has8 + 1));
42	dEdxLayer *= 2.;
43	}
44	//End of patch
45
46	if (dEdxLayer > 0. && dEdxLayer < 4090) {
47	dEdx += dEdxLayer;
48	nLayers++;
49	}
50	dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 1, 100);
51	if (dEdxLayer > 0. && dEdxLayer < 4090) {
52	dEdx += dEdxLayer;
53	nLayers++;
54	}
55
56	if (nLayers > 0) {
57	dEdx /= nLayers;
58
59	float R = event->GetTrack(0)->GetTrkTrack()->GetRigidity();
60
61	if (dEdx >= (3.102 - 2.6 * TMath::Log(R)) * (R < 1.05) + 3. * (R >= 1.05))
62	return 0;
63
64	return CUTOK;
65	// if (dEdx < (3.102 - 2.6 * TMath::Log(R)) * (R < 1.05) + 3. * (R >= 1.05))
66	// return 0;
67	//
68	// return CUTOK;
69	}
70	else {
71	return 0;
72	}
73
74	}