ExtBin Struct Reference

Public Member Functions
void	addEvent (float mass, float me, float me2, float weight=1)
void	sift ()
void	adjustWeights (float refth=3.)
void	mergeBin (const ExtBin &other)
void	addEvent (SimpleEntry evt)
void	addEvent (int ev, float trueval, float recoval, float recotrackval=0, float weight=1)
void	sortByRecoMinusTrueVals (vector< pair< float, int >> &res)
void	sortByRecoVals (vector< pair< float, int >> &res)
void	sift (double threshold=0.998, bool use_recominustrue=false)
void	adjustWeights (float refth=-1)
void	mergeBin (const ExtBin &other)

Public Attributes
double	binlow
double	binhigh
vector< int >	events
vector< float >	masses
vector< float >	mevals
vector< float >	me2vals
vector< float >	weights
vector< SimpleEntry >	collection

Detailed Description

Definition at line 111 of file calcC_lintolog.c.

Member Function Documentation

addEvent() [1/3]

void ExtBin::addEvent ( float  mass, float  me, float  me2, float  weight = 1  )

inline

Definition at line 120 of file calcC_lintolog.c.

                                                                {
    masses.push_back(mass);
    mevals.push_back(me);
    me2vals.push_back(me2);
    weights.push_back(weight);
  }

addEvent() [2/3]

void ExtBin::addEvent ( int  ev, float  trueval, float  recoval, float  recotrackval = 0, float  weight = 1  )

inline

Definition at line 185 of file getResolution.c.

                                                                                           {
    SimpleEntry tmp(ev, trueval, recoval, recotrackval, weight);
    collection.push_back(tmp);
  }

addEvent() [3/3]

void ExtBin::addEvent ( SimpleEntry  evt )

inline

Definition at line 182 of file getResolution.c.

                                {
    collection.push_back(evt);
  }

adjustWeights() [1/2]

void ExtBin::adjustWeights ( float  refth = -1 )

inline

Definition at line 228 of file getResolution.c.

                                    {
    if (collection.size()==0) return;
 
    vector<pair<float, int>> weight_entry;
    for (unsigned int ev=0; ev<collection.size(); ev++) addByLowest<float, int>(weight_entry, collection.at(ev).weight, ev);
    int bin=weight_entry.size()-1;
 
    int at99p0ev = (float(weight_entry.size()))*0.99;
    if (at99p0ev==(int)weight_entry.size()) at99p0ev--;
    float threshold = weight_entry.at(at99p0ev).first;
    while (bin>at99p0ev){
      if (
        threshold*2.<weight_entry.at(bin).first
        ){
        if (weight_entry.at(bin).first==0.) cerr << "Weight=0!" << endl;
        collection.at(weight_entry.at(bin).second).weight = pow(threshold, 2)/weight_entry.at(bin).first;
        cout << "Adjusting weight for event " << bin << endl;
      }
      bin--;
    }
    float sum[2]={ 0 };
    for (unsigned int ev=0; ev<collection.size(); ev++){
      sum[0] += collection.at(ev).weight;
      sum[1] += 1;
    }
    if (sum[0]==0. || sum[1]==0.) cerr << "sum[0]=" << sum[0] << " or sum[1]=" << sum[1] << " is invalid." << endl;
    for (unsigned int ev=0; ev<collection.size(); ev++) collection.at(ev).weight *= sum[1]/sum[0];
    if (refth>0.){
      for (unsigned int ev=0; ev<collection.size(); ev++){
        if (collection.at(ev).weight<=refth) continue;
        collection.at(ev).weight = pow(refth, 2)/collection.at(ev).weight;
      }
    }
  }

adjustWeights() [2/2]

void ExtBin::adjustWeights ( float  refth = 3. )

inline

Definition at line 158 of file calcC_lintolog.c.

                                    {
    if (masses.size()==0) return;
 
    vector<pair<float, int>> weight_entry;
    for (unsigned int ev=0; ev<masses.size(); ev++) addByLowest<float, int>(weight_entry, weights.at(ev), ev);
    int bin=weight_entry.size()-1;
 
    int at99p0ev = (float(weight_entry.size()))*0.99;
    if (at99p0ev==(int)weight_entry.size()) at99p0ev--;
    float threshold = weight_entry.at(at99p0ev).first;
    while (bin>at99p0ev){
      if (
        threshold*2.<weight_entry.at(bin).first
        ){
        weights.at(weight_entry.at(bin).second) = pow(threshold, 2)/weight_entry.at(bin).first;
        cout << "Adjusting weight for event " << bin << endl;
      }
      bin--;
    }
    float sum[2]={ 0 };
    for (unsigned int ev=0; ev<weights.size(); ev++){
      sum[0] += weights.at(ev);
      sum[1] += 1;
    }
    for (unsigned int ev=0; ev<weights.size(); ev++) weights.at(ev) *= sum[1]/sum[0];
    if (refth>0.){
      for (unsigned int ev=0; ev<weights.size(); ev++){
        if (weights.at(ev)<=refth) continue;
        weights.at(ev) = pow(refth, 2)/weights.at(ev);
      }
    }
  }

mergeBin() [1/2]

void ExtBin::mergeBin ( const ExtBin &  other )

inline

Definition at line 190 of file calcC_lintolog.c.

                                    {
    if (
      other.events.size()!=other.masses.size()
      ||
      other.events.size()!=other.mevals.size()
      ||
      other.events.size()!=other.me2vals.size()
      ||
      other.events.size()!=other.weights.size()
      ){
      cerr << "Size of the following elements are not the same:" << endl;
      cerr << "\t- Events: " << other.events.size() << endl;
      cerr << "\t- Masses: " << other.masses.size() << endl;
      cerr << "\t- MEs: " << other.mevals.size() << endl;
      cerr << "\t- MEs (2): " << other.me2vals.size() << endl;
      cerr << "\t- Weights: " << other.weights.size() << endl;
    }
    for (unsigned int ev=0; ev<other.events.size(); ev++){
      events.push_back(other.events.at(ev));
      addEvent(other.masses.at(ev), other.mevals.at(ev), other.me2vals.at(ev), other.weights.at(ev));
    }
  }

mergeBin() [2/2]

void ExtBin::mergeBin ( const ExtBin &  other )

inline

Definition at line 262 of file getResolution.c.

{ for (unsigned int ev=0; ev<other.collection.size(); ev++) collection.push_back(other.collection.at(ev)); }

sift() [1/2]

void ExtBin::sift ( )

inline

Definition at line 126 of file calcC_lintolog.c.

             {
    if (masses.size()==0) return;
 
    vector<int> take_out;
    vector<pair<float, int>> me_entry[2];
    for (unsigned int ev=0; ev<masses.size(); ev++){
      addByLowest<float, int>(me_entry[0], mevals.at(ev), ev);
      addByLowest<float, int>(me_entry[1], me2vals.at(ev), ev);
    }
    for (unsigned int im=0; im<2; im++){
      int at99p8ev = (float(me_entry[im].size()))*0.998;
      if (at99p8ev==(int)me_entry[im].size()) at99p8ev--;
      int bin=me_entry[im].size()-1;
      while (bin>at99p8ev){
        if (
          me_entry[im].at(at99p8ev).first*2.<me_entry[im].at(bin).first
          ) addByLowest<int>(take_out, me_entry[im].at(bin).second, true);
        bin--;
      }
    }
    for (int bin=take_out.size()-1; bin>=0; bin--){
      int t_ev = take_out.at(bin);
      cout << "Discarding event " << events.at(t_ev) << endl;
      vector<int>::iterator itev;
      vector<float>::iterator it;
      itev=events.begin()+t_ev; events.erase(itev);
      it=masses.begin()+t_ev; masses.erase(it);
      it=mevals.begin()+t_ev; mevals.erase(it);
      it=me2vals.begin()+t_ev; me2vals.erase(it);
      it=weights.begin()+t_ev; weights.erase(it);
    }
  }

sift() [2/2]

void ExtBin::sift ( double  threshold = 0.998, bool  use_recominustrue = false  )

inline

Definition at line 200 of file getResolution.c.

                                                                 {
    if (collection.size()==0) return;
 
    vector<int> take_out;
    vector<pair<float, int>> theEntries;
    if (!use_recominustrue) sortByRecoVals(theEntries);
    else sortByRecoMinusTrueVals(theEntries);
    for (unsigned int im=0; im<2; im++){
      int at99p8ev = std::floor((float(theEntries.size()))*threshold);
      int bin=theEntries.size()-1;
      while (bin>at99p8ev){
        if (
          theEntries.at(at99p8ev).first*2.<theEntries.at(bin).first
          ) addByLowest<int>(take_out, theEntries.at(bin).second, true);
        bin--;
      }
    }
    for (int bin=take_out.size()-1; bin>=0; bin--){
      int t_ev = take_out.at(bin);
      cout << "Attempting to discard event at position " << t_ev << endl;
      if (t_ev>=(int)collection.size()) cerr << "Collection size " << collection.size() << " >= " << t_ev << endl;
      else{
        cout << "Discarding event " << collection.at(t_ev).id << " at position " << t_ev << endl;
        vector<SimpleEntry>::iterator it;
        it=collection.begin()+t_ev; collection.erase(it);
      }
    }
  }

sortByRecoMinusTrueVals()

void ExtBin::sortByRecoMinusTrueVals ( vector< pair< float, int >> &  res )

inline

Definition at line 189 of file getResolution.c.

                                                             {
    res.clear();
    for (unsigned int ev=0; ev<collection.size(); ev++){
      float val = collection.at(ev).recoval/collection.at(ev).trueval-1.;
      addByLowest<float, int>(res, val, ev);
    }
  }

sortByRecoVals()

void ExtBin::sortByRecoVals ( vector< pair< float, int >> &  res )

inline

Definition at line 196 of file getResolution.c.

                                                    {
    res.clear();
    for (unsigned int ev=0; ev<collection.size(); ev++) addByLowest<float, int>(res, collection.at(ev).recoval, ev);
  }

Member Data Documentation

binhigh

double ExtBin::binhigh

Definition at line 113 of file calcC_lintolog.c.

binlow

double ExtBin::binlow

Definition at line 112 of file calcC_lintolog.c.

collection

vector<SimpleEntry> ExtBin::collection

Definition at line 180 of file getResolution.c.

events

vector<int> ExtBin::events

Definition at line 114 of file calcC_lintolog.c.

masses

vector<float> ExtBin::masses

Definition at line 115 of file calcC_lintolog.c.

me2vals

vector<float> ExtBin::me2vals

Definition at line 117 of file calcC_lintolog.c.

mevals

vector<float> ExtBin::mevals

Definition at line 116 of file calcC_lintolog.c.

weights

vector<float> ExtBin::weights

Definition at line 118 of file calcC_lintolog.c.

---

The documentation for this struct was generated from the following files:

• MELA/test/calcC_lintolog.c
• MELA/test/getResolution.c