ChargeAssignment.h File Reference

Go to the source code of this file.

Functions
int	getCharge (int phi1, int phi2, int phi3, int phi4, int mode)

Function Documentation

int getCharge	(	int	phi1,
		int	phi2,
		int	phi3,
		int	phi4,
		int	mode
	)

Definition at line 17 of file ChargeAssignment.h.

References RecoTauCleanerPlugins::charge.

                                                               {
  
  // -1 = postive physical charge to match pdgId code (i.e. -13 is positive, anti-muon). +1 = negative physical charge.
  // Also matches DN-2015/017 format for track finder --> uGMT interface format, where 0 indicates positive, 1 negative.
  int emuCharge = 0;
  int phidiffs[6] = {phi2 - phi1, phi3 - phi1, phi4 - phi1, phi3 - phi2, phi4 - phi2, phi4 - phi3};
  
  if(mode == 15 ){
    if(phidiffs[0] > 0)
      emuCharge = 1;
    else if(phidiffs[0] == 0 && phidiffs[1] < 0)
      emuCharge = 1;
    else if(phidiffs[1] == 0 && phidiffs[2] < 0)
      emuCharge = 1;
    else
      emuCharge = -1;
  }
  
  if(mode == 14){
    if(phidiffs[0] < 0)
      emuCharge = -1;
    else if(phidiffs[0] == 0 && phidiffs[1] < 0)
      emuCharge = -1;
    else
      emuCharge = 1;
  }
  
  if(mode == 13){
    if(phidiffs[0] > 0)
      emuCharge = 1;
    else if(phidiffs[0] == 0 && phidiffs[2] < 0)
      emuCharge = 1;
    else
      emuCharge = -1;
  }
  
  if(mode == 12){
    if(phidiffs[0] > 0)
      emuCharge = 1;
    else
      emuCharge = -1;
  }
  
  if(mode == 11){
    if(phidiffs[1] > 0)
      emuCharge = 1;
    else if(phidiffs[1] == 0 && phidiffs[2] < 0)
      emuCharge = 1;
    else
      emuCharge = -1;
  }
  
  if(mode == 10){
    if(phidiffs[1] > 0)
      emuCharge = 1;
    else
      emuCharge = -1;
  }
  
  if(mode == 9){
    if(phidiffs[2] > 0)
      emuCharge = 1;
    else
      emuCharge = -1;
  }
  
  if(mode == 7){
    if(phidiffs[3] > 0)
      emuCharge = 1;
    else if(phidiffs[3] == 0 && phidiffs[4] < 0)
      emuCharge = 1;
    else
      emuCharge = -1;
  }
  
  if(mode == 6){
    if(phidiffs[3] > 0)
      emuCharge = 1;
    else
      emuCharge = -1;
  }
  
  if(mode == 5){
    if(phidiffs[4] > 0)
      emuCharge = 1;
    else
      emuCharge = -1;
  }
  
  if(mode == 3){
    if(phidiffs[5] > 0)
      emuCharge = 1;
    else
      emuCharge = -1;
  }
  
  int charge = 0;
  if(emuCharge == 1)
    charge = 1;
  
  return charge;
  
}