Histo2DActions/RigFillAction/RigFillAction.cpp

/*
 * RigFillAction.cpp
 *
 *  Created on: 14/lug/2009
 *      Author: Nicola Mori
 */

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

#include "RigFillAction.h"

RigFillAction::RigFillAction(const char *actionName, TString outFileBase, vector<float> &bins, float thresholdCoeff,
    float chargeSign, bool spilloverFlag, float spilloverLimit, float mdrMin) :
  Histo2DAction<Float_t> (actionName, "", outFileBase, "RECREATE", false, true), _thresholdCoeff(thresholdCoeff),
      _chargeSign(chargeSign), _spilloverFlag(spilloverFlag), _spilloverLimit(spilloverLimit), _mdrMin(mdrMin) {

  SetXAxis("Rc [GV]", bins);
  SetYAxis("R [GV]", bins);
}

RigFillAction::RigFillAction(const char *actionName, TString outFileBase, TString rigBinsFile, float thresholdCoeff,
    float chargeSign, bool spilloverFlag, float spilloverLimit, float mdrMin) :
  Histo2DAction<Float_t> (actionName, "", outFileBase, "RECREATE", false, true), _thresholdCoeff(thresholdCoeff),
      _chargeSign(chargeSign), _spilloverFlag(spilloverFlag), _spilloverLimit(spilloverLimit), _mdrMin(mdrMin) {

  SetXAxis("Rc [GV]", rigBinsFile);
  SetYAxis("R [GV]", rigBinsFile);
}

void RigFillAction::OnGood(PamLevel2 *event) {

  float rigThreshold = _thresholdCoeff * event->GetOrbitalInfo()->GetCutoffSVL();
  float rig = 1. / event->GetTrack(0)->GetTrkTrack()->GetDeflection();

  // Compute deflection correction due to MDR cut
  if (_mdrMin > 0.) {
    Double_t F = 0.;
    if (_mdrMin > 300.) {
      F = 2.215e-4 * (log(_mdrMin) - log(300.));
    }
    rig = 1. / (1. / rig + F);
  }

  if (rig / _chargeSign > 0) { // Check if the particle has the right sign
    Fill(rigThreshold, rig);
  }
  else if (_spilloverFlag && fabs(rig) > _spilloverLimit) { // Check if it is a spillover event
    Fill(rigThreshold, _yBins.back() + 1.); // Place it in the Y overflow
  }

}

void RigFillAction::Finalize() {

  // Save the histogram
  Histo2DAction<Float_t>::Finalize();

  // Save the zero bins
  if (_outFileBase != "") {
    ofstream outTextFile((_outFileBase + "-InfBins.txt").Data(), ios_base::out);
    // XUnderflow is the vector of events whose Rc is below the lower histogram limit
    for (unsigned int i = 0; i < GetXUnderflow().size(); i++)
      outTextFile << GetXUnderflow()[i] << "\n";
    outTextFile.close();

    // Save the (spillover + rigidity overflow)
    if (_spilloverFlag) {
      ofstream outTextFile((_outFileBase + "-Spillover.txt").Data(), ios_base::out);
      // XUnderYOverflow is the number of (spillover + rigidity overflow) events with an underflowing critical rigidity
      outTextFile << GetXUnderYOverflow() << "\n";
      for (unsigned int i = 0; i < GetXUnderflow().size(); i++)
        // YOverflow are the (spillover + rigidity overflow) bins
        outTextFile << GetYOverflow()[i] << "\n";
      outTextFile.close();
    }
  }
}
1	/*
2	* RigFillAction.cpp
3	*
4	* Created on: 14/lug/2009
5	* Author: Nicola Mori
6	*/
7
8	/! @file RigFillAction.cpp The RigFillAction class implementation file. /
9
10	#include "RigFillAction.h"
11
12	RigFillAction::RigFillAction(const char *actionName, TString outFileBase, vector<float> &bins, float thresholdCoeff,
13	float chargeSign, bool spilloverFlag, float spilloverLimit, float mdrMin) :
14	Histo2DAction<Float_t> (actionName, "", outFileBase, "RECREATE", false, true), _thresholdCoeff(thresholdCoeff),
15	_chargeSign(chargeSign), _spilloverFlag(spilloverFlag), _spilloverLimit(spilloverLimit), _mdrMin(mdrMin) {
16
17	SetXAxis("Rc [GV]", bins);
18	SetYAxis("R [GV]", bins);
19	}
20
21	RigFillAction::RigFillAction(const char *actionName, TString outFileBase, TString rigBinsFile, float thresholdCoeff,
22	float chargeSign, bool spilloverFlag, float spilloverLimit, float mdrMin) :
23	Histo2DAction<Float_t> (actionName, "", outFileBase, "RECREATE", false, true), _thresholdCoeff(thresholdCoeff),
24	_chargeSign(chargeSign), _spilloverFlag(spilloverFlag), _spilloverLimit(spilloverLimit), _mdrMin(mdrMin) {
25
26	SetXAxis("Rc [GV]", rigBinsFile);
27	SetYAxis("R [GV]", rigBinsFile);
28	}
29
30	void RigFillAction::OnGood(PamLevel2 *event) {
31
32	float rigThreshold = _thresholdCoeff * event->GetOrbitalInfo()->GetCutoffSVL();
33	float rig = 1. / event->GetTrack(0)->GetTrkTrack()->GetDeflection();
34
35	// Compute deflection correction due to MDR cut
36	if (_mdrMin > 0.) {
37	Double_t F = 0.;
38	if (_mdrMin > 300.) {
39	F = 2.215e-4 * (log(_mdrMin) - log(300.));
40	}
41	rig = 1. / (1. / rig + F);
42	}
43
44	if (rig / _chargeSign > 0) { // Check if the particle has the right sign
45	Fill(rigThreshold, rig);
46	}
47	else if (_spilloverFlag && fabs(rig) > _spilloverLimit) { // Check if it is a spillover event
48	Fill(rigThreshold, _yBins.back() + 1.); // Place it in the Y overflow
49	}
50
51	}
52
53	void RigFillAction::Finalize() {
54
55	// Save the histogram
56	Histo2DAction<Float_t>::Finalize();
57
58	// Save the zero bins
59	if (_outFileBase != "") {
60	ofstream outTextFile((_outFileBase + "-InfBins.txt").Data(), ios_base::out);
61	// XUnderflow is the vector of events whose Rc is below the lower histogram limit
62	for (unsigned int i = 0; i < GetXUnderflow().size(); i++)
63	outTextFile << GetXUnderflow()[i] << "\n";
64	outTextFile.close();
65
66	// Save the (spillover + rigidity overflow)
67	if (_spilloverFlag) {
68	ofstream outTextFile((_outFileBase + "-Spillover.txt").Data(), ios_base::out);
69	// XUnderYOverflow is the number of (spillover + rigidity overflow) events with an underflowing critical rigidity
70	outTextFile << GetXUnderYOverflow() << "\n";
71	for (unsigned int i = 0; i < GetXUnderflow().size(); i++)
72	// YOverflow are the (spillover + rigidity overflow) bins
73	outTextFile << GetYOverflow()[i] << "\n";
74	outTextFile.close();
75	}
76	}
77	}