ProtonCandidate Class Reference

Inheritance diagram for ProtonCandidate:

Classes
enum	Cut
	This is the enum for proton cut types. More...

Public Member Functions
	ProtonCandidate (int ipart)
LorentzVector	getLorentzVector ()
boolean	cut_PID ()
boolean	cut_FORWARD ()
boolean	cut_CHI2PID ()
boolean	cut_DC_FIDUCIAL_REG1 ()
boolean	cut_DC_FIDUCIAL_REG2 ()
boolean	cut_DC_FIDUCIAL_REG3 ()
boolean	cut_DELTA_VZ ()
boolean	isproton ()
boolean	isproton (Cut ...applycuts)
Public Member Functions inherited from Candidate
	Candidate (int ipart)
void	setINBENDING ()
void	setOUTBENDING ()
void	setPxyz (Number px, Number py, Number pz)
void	setVZ (Number vz)
void	setBETA (Number beta)
void	setCHI2PID (Number chi2pid)
void	setDVZ (Number dvz)
void	setPID (Number pid)
void	setStatus (Number status)
void	setNPHE (Number nphe)
void	setPCALsector (Number sector)
void	setPCALenergy (Number energy)
void	setECINenergy (Number energy)
void	setECOUTenergy (Number energy)
void	setPCALvw (Number lv, Number lw)
void	setPCALxyz (Number x, Number y, Number z)
Integer	getPCALsector ()
Integer	getDCsector ()
Float	getDC1x ()
Float	getDC1y ()
Float	getDC1z ()
void	setDCxyz (int region, Number x, Number y, Number z)

Static Public Member Functions
static ProtonCandidate	getProtonCandidate (int ipart, Bank recbank, Bank trajbank, boolean isinbending)

Additional Inherited Members
Public Attributes inherited from Candidate
final int	ipart
Protected Attributes inherited from Candidate
MagField	field = MagField.INBENDING
	magnetic field, INBENDING by default
Float	nphe = null
	number of photoelectrons
Integer	pcal_sector = null
	pcal sector
Float	pcal_energy = null
	energy deposited in PCAL
Float	ecin_energy = 0f
	energy deposited in ECin
Float	ecout_energy = 0f
	energy deposited in ECout
Float	pcal_lv = null
	distance on V-side
Float	pcal_lw = null
	distance on W-side
Float	pcal_x = null
	X from PCAL.
Float	pcal_y = null
	Y from PCAL.
Float	pcal_z = null
	Z from PCAL.
Float	p = null
	momentum
Float	px = null
	px momentum
Float	py = null
	px momentum
Float	pz = null
	px momentum
Float	vz = null
	z vertex
Float	beta = null
	beta value
Float	chi2pid = null
	chi2pid value
Float	dvz = null
	the difference between trigger particle vz and candidate vz
Float	traj_x1 = null
	x-position of the track at the detector surface (cm) at region 1
Float	traj_y1 = null
	y-position of the track at the detector surface (cm) at region 1
Float	traj_z1 = null
	z-position of the track at the detector surface (cm) at region 1
Float	traj_x2 = null
	x-position of the track at the detector surface (cm) at region 2
Float	traj_y2 = null
	y-position of the track at the detector surface (cm) at region 2
Float	traj_z2 = null
	z-position of the track at the detector surface (cm) at region 2
Float	traj_x3 = null
	x-position of the track at the detector surface (cm) at region 3
Float	traj_y3 = null
	y-position of the track at the detector surface (cm) at region 3
Float	traj_z3 = null
	z-position of the track at the detector surface (cm) at region 3
Integer	pid = null
	particle PID
Integer	dc_sector = null
	dc sector
Short	status = null
	particle status

---

Class Documentation

uconn::utils::pid::stefan::ProtonCandidate::Cut

enum uconn::utils::pid::stefan::ProtonCandidate::Cut

This is the enum for proton cut types.

Enumerator
CHI2PID	cut on chi2pid
DC_FIDUCIAL_REG1	fiducial DC cut for region 1
DC_FIDUCIAL_REG2	fiducial DC cut for region 2
DC_FIDUCIAL_REG3	fiducial DC cut for region 3
DELTA_VZ	cut on difference between VZ of Proton candidate and trigger particle
FORWARD	only forward detector
PID	cut on PDG code

Constructor & Destructor Documentation

ProtonCandidate()

ProtonCandidate ( int  ipart )

inline

A Constructor

Parameters

ipart

particle index

                                      {
        super(ipart);
    }

Member Function Documentation

cut_CHI2PID()

boolean cut_CHI2PID ( )

inline

Returns

chi2pid cut

                                 {
        if(pid==null || chi2pid==null || p==null) return false;
        return HadronCuts.Chi2pid_cut(chi2pid, p, pid);
    }

cut_DC_FIDUCIAL_REG1()

boolean cut_DC_FIDUCIAL_REG1 ( )

inline

Returns

DC fiducial region 1

                                          {
        if(dc_sector==null || traj_x1==null || traj_y1==null || traj_z1==null || pid==null) return false;
        if(field==MagField.INBENDING)
          return HadronCuts.DC_fiducial_cut_theta_phi(dc_sector, 1, traj_x1, traj_y1, traj_z1, pid, field==MagField.INBENDING);
        return ElectronCuts.DC_fiducial_cut_XY(dc_sector, 1, traj_x1, traj_y1, pid, field==MagField.INBENDING);
    }

cut_DC_FIDUCIAL_REG2()

boolean cut_DC_FIDUCIAL_REG2 ( )

inline

Returns

DC fiducial region 2

                                          {
        if(dc_sector==null || traj_x2==null || traj_y2==null || traj_z2==null || pid==null) return false;
        if(field==MagField.INBENDING)
          return HadronCuts.DC_fiducial_cut_theta_phi(dc_sector, 2, traj_x2, traj_y2, traj_z2, pid, field==MagField.INBENDING);
        return ElectronCuts.DC_fiducial_cut_XY(dc_sector, 2, traj_x2, traj_y2, pid, field==MagField.INBENDING);
    }

cut_DC_FIDUCIAL_REG3()

boolean cut_DC_FIDUCIAL_REG3 ( )

inline

Returns

DC fiducial region 3

                                          {
        if(dc_sector==null || traj_x3==null || traj_y3==null || traj_z3==null || pid==null) return false;
        if(field==MagField.INBENDING)
          return HadronCuts.DC_fiducial_cut_theta_phi(dc_sector, 3, traj_x3, traj_y3, traj_z3, pid, field==MagField.INBENDING);
        return ElectronCuts.DC_fiducial_cut_XY(dc_sector, 3, traj_x3, traj_y3, pid, field==MagField.INBENDING);
    }

cut_DELTA_VZ()

boolean cut_DELTA_VZ ( )

inline

Returns

cut on vertex position

                                  {
        if(pid==null || dvz==null) return false;
        return HadronCuts.Delta_vz_cut(pid, dvz);
    }

cut_FORWARD()

boolean cut_FORWARD ( )

inline

Returns

if it is detected in forward

                                 {
        if(status==null) return false;
        return status>=2000 && status<4000;
    }

cut_PID()

boolean cut_PID ( )

inline

Returns

pid cut

                             {
        if(pid == null) return false;
        return pid == 2212;
    }

getLorentzVector()

LorentzVector getLorentzVector ( )

inline

Returns

LorentzVector instance

                                            {
        LorentzVector vec = null;
        if(px!=null && py!=null && pz!=null) {
            vec = new LorentzVector();
            vec.setPxPyPzM(px,py,pz,0.938272);
        }
        return vec;
    }

getProtonCandidate()

static ProtonCandidate getProtonCandidate ( int  ipart, Bank  recbank, Bank  trajbank, boolean  isinbending  )

inlinestatic

return ProtonCandidate instance

Parameters

ipart	particle index
recbank,trajbank	particle and trajectory banks
isinbending	true for inbending, false for outbending

                                                                                                                  {
        ProtonCandidate candidate = new ProtonCandidate(ipart);
        if(!isinbending) candidate.setOUTBENDING();
 
        if(recbank!=null) {
            candidate.setPID(recbank.getInt("pid",ipart));
            candidate.setStatus(recbank.getShort("status",ipart));
            candidate.setCHI2PID(recbank.getFloat("chi2pid",ipart));
            candidate.setVZ(recbank.getFloat("vz",ipart));
            candidate.setDVZ(recbank.getFloat("vz",ipart) - recbank.getFloat("vz",0));
            candidate.setPxyz(recbank.getFloat("px",ipart), recbank.getFloat("py",ipart), recbank.getFloat("pz",ipart));
        }
 
        if(trajbank!=null) IntStream.range(0,trajbank.getRows())
            .filter(i -> trajbank.getShort("pindex",i) == ipart && trajbank.getByte("detector",i) == DetectorType.DC.getDetectorId())
            .forEach(i -> {
            if(trajbank.getByte("layer",i) == 6)
                candidate.setDCxyz(1, trajbank.getFloat("x",i), trajbank.getFloat("y",i), trajbank.getFloat("z",i));
            else if(trajbank.getByte("layer",i) == 18)
                candidate.setDCxyz(2, trajbank.getFloat("x",i), trajbank.getFloat("y",i), trajbank.getFloat("z",i));
            else if(trajbank.getByte("layer",i) == 36)
                candidate.setDCxyz(3, trajbank.getFloat("x",i), trajbank.getFloat("y",i), trajbank.getFloat("z",i));
        });
 
        return candidate;
    }

isproton() [1/2]

boolean isproton ( )

inline

testing against all proton cuts

                              {
        return isproton(Cut.values());
    }

isproton() [2/2]

boolean isproton ( Cut ...  applycuts )

inline

assembly of multiple proton cuts

Parameters

applycuts

the list of cuts required to apply

                                              {
        for(Cut thiscut: applycuts) {
            if(thiscut == Cut.PID) {
                if(!cut_PID()) return false;
 
            } else if(thiscut == Cut.CHI2PID) {                                                                                                                                          
                if(!cut_CHI2PID()) return false;
 
            } else if(thiscut == Cut.FORWARD) {
                if(!cut_FORWARD()) return false;
 
            } else if(thiscut == Cut.DC_FIDUCIAL_REG1) {
                if(!cut_DC_FIDUCIAL_REG1()) return false;
 
            } else if(thiscut == Cut.DC_FIDUCIAL_REG2) {
                if(!cut_DC_FIDUCIAL_REG2()) return false;
 
            } else if(thiscut == Cut.DC_FIDUCIAL_REG3) {
                if(!cut_DC_FIDUCIAL_REG3()) return false;
 
            } else if(thiscut == Cut.DELTA_VZ) {
                if(!cut_DELTA_VZ()) return false;
 
            } else {
                return false;
            }
        }
        return true;
    }

---

The documentation for this class was generated from the following file:

• src/main/java/uconn/utils/pid/stefan/ProtonCandidate.java