#include <StringBasedNTupler.h>

Public Types
typedef TreeBranch::value	value

Public Member Functions
value	operator() ()

	StringBranchHelper (const TreeBranch &B, const edm::Event &iEvent)

Private Attributes
value	value_

Detailed Description

template<typename Object, typename Collection = std::vector<Object>>
class StringBranchHelper< Object, Collection >

Definition at line 126 of file StringBasedNTupler.h.

Member Typedef Documentation

template<typename Object , typename Collection = std::vector<Object>>

typedef TreeBranch::value StringBranchHelper< Object, Collection >::value

Definition at line 128 of file StringBasedNTupler.h.

Constructor & Destructor Documentation

template<typename Object , typename Collection = std::vector<Object>>

StringBranchHelper< Object, Collection >::StringBranchHelper	(	const TreeBranch &	B,
		const edm::Event &	iEvent
	)

inline

Definition at line 131 of file StringBasedNTupler.h.

References TreeBranch::className(), TreeBranch::expr(), edm::HandleBase::failedToGet(), edm::Event::getByLabel(), mps_fire::i, edm::EventBase::isRealData(), LogDebug, TreeBranch::order(), TreeBranch::selection(), TreeBranch::src(), and parallelization::uint().

     {
       const float defaultValue = 0.;
 
       //    grab the collection
       edm::Handle<Collection> oH;
       iEvent.getByLabel(B.src(), oH);
 
 
       //empty vector if product not found
       if (oH.failedToGet()){
     if(!(iEvent.isRealData() && B.className()=="reco::GenParticle") ) {  //don't output genparticle error in data
       edm::LogError("StringBranchHelper")<<"cannot open: "<<B.src()<<"  "<<B.className();
         }
         value_.reset(new std::vector<float>());
       }
       else{
     //parser for the object expression
     StringObjectFunction<Object> expr(B.expr());
     //allocate enough memory for the data holder
         value_.reset(new std::vector<float>());
         value_->reserve(oH->size());
 
     StringCutObjectSelector<Object> * selection=nullptr;
     if (B.selection()!=""){
       //std::cout<<"trying to get to a selection"<<std::endl;
       selection = new StringCutObjectSelector<Object>(B.selection());
       //std::cout<<"got the objet"<<std::endl;
     }
     uint i_end=oH->size();
     //sort things first if requested
     if (B.order()!=""){
       StringObjectFunction<Object> order(B.order());
       // allocate a vector of pointers (we are using view) to be sorted
       std::vector<const Object*> copyToSort(oH->size());
       for (uint i=0;i!=i_end;++i)  copyToSort[i]= &(*oH)[i];
       std::sort(copyToSort.begin(), copyToSort.end(), sortByStringFunction<Object>(&order));
       //then loop and fill
       for (uint i=0;i!=i_end;++i) {
         //try and catch is necessary because ...
         try{
           if (selection && !((*selection)(*(copyToSort)[i]))) continue;
           value_->push_back((expr)(*(copyToSort)[i]));
         }catch(...){
           LogDebug("StringBranchHelper")<<"with sorting. could not evaluate expression: "<<B.expr()<<" on class: "<<B.className();
           value_->push_back(defaultValue);//push a default value to not change the indexing
         }
       }
     }
     else{
       //actually fill the vector of values
       for (uint i=0;i!=i_end;++i){
         //try and catch is necessary because ...
         try {
           if (selection && !((*selection)((*oH)[i]))) continue;
           value_->push_back((expr)((*oH)[i]));
         }catch(...){
           LogDebug("StringBranchHelper")<<"could not evaluate expression: "<<B.expr()<<" on class: "<<B.className();
           value_->push_back(defaultValue);//push a default value to not change the indexing
         }
       }
     }
     if (selection) delete selection;
       }
     }