TMarCutMC.C

Go to the documentation of this file.

// Class TMarCutMC
//
// Author     : F. Cassol Brunner
// Created    : 21/10/2002
// Last update: 
//          by: 
// Comment: Class to handle cut on generated particles

#include "TMarCutMC.h"

TMarCutMC::TMarCutMC()
{
  Variable="none";
  Valid=false;

};

istream& operator>>(istream& in, TMarCutMC& CutMC){

   int n_cut,out=1;
   float Min,Max;
   double Lumi;
   char nom[20];

   in >> CutMC.Variable >> nom;
   
   if(CutMC.Variable!="MELEPHO"&&
      CutMC.Variable!="MMUMU"&&
      CutMC.Variable!="Q2"&&
      CutMC.Variable!="RAPQ2"&&
      CutMC.Variable!="Y"&&
      CutMC.Variable!="RAPY"&&
      CutMC.Variable!="PHAT")
   {
      Error("TMarCutMC::istream","Wrong CutMC.Variable '%s' in card 'genCut'. Exit",CutMC.Variable.c_str());
      exit(1);
   }
    
      
   out= sscanf(nom,"%d",&n_cut);
   
   if(!out){
error:;   
     Error("TMarCutMC::istream","Wrong input format (card 'genCut'). Exit");
     exit(1);
   }
   
   for(Int_t i=0;i<n_cut;i++){
       in >> nom;
       out= sscanf(nom,"%f",&Min);
       if(!out)goto error;
       
       in >> nom;
       out= sscanf(nom,"%f",&Max);
       if(!out)goto error;      

       in >> nom;
       out= sscanf(nom,"%lf",&Lumi);
       if(!out)goto error;    
       
          
       CutMC.CutMin.push_back(Min);
       CutMC.CutMax.push_back(Max);

       if(Lumi!=0)CutMC.Weight.push_back(1/Lumi);
       else CutMC.Weight.push_back(0) ;
   }
   
  
   
   return in;
};

ostream& operator<<(ostream& out, TMarCutMC& CutMC){

   out << "variable= " << CutMC.Variable;
   for (unsigned int i=0;i<CutMC.CutMin.size();i++){
      out << " min= "   << CutMC.CutMin[i] 
          << " max= "   << CutMC.CutMax[i]    << " weight= "<< CutMC.Weight[i] ;
   }
   return out;
};


/*TMarCutMC& TMarCutMC::operator= (const TMarCutMC& other) 
{
   Variable=other.GetVariable();
   CutMin=other.GetCutMin();
   CutMax=other.GetCutMax();
   Valid=other.GetValid();
   Value=other.GetValue();
   Weight=other.GetWeight();

   return *this;

};*/

// bool TMarCutMC::Evaluate(Float_t value)
// {
// 
// //... set Valid to true if the cut is passed
// 
//    Value=value;
//    Valid=false;
//    for (unsigned int i=0;i<CutMin.size();i++){
// 
//      if (value>CutMin[i] && value<CutMax[i]){
//        Valid=true;
//        return Valid;
//      }
// 
//    }
//    return Valid;
// };

Double_t TMarCutMC::Evaluate(Double_t value)
{
  
//... set Valid to true if the cut is passed
  for (unsigned int i=0;i<CutMin.size();i++){
//      printf("evaluate cut %d min=%f max=%f w=%f \n",i,CutMin[i],CutMax[i],Weight[i]);
     if (value>CutMin[i] && value<CutMax[i]){

       return Weight[i];
     }

   }
   return 0;
};


istream& operator>>(istream& in, TMarCutMCList& List){

   int n_cut,out=1;
   char nom[20];
   in >>  nom;
   
   out= sscanf(nom,"%d",&n_cut);
   
   if(!out){
     Error("TMarCutMCList::istream","Wrong input format (card 'genCut'). Exit");
     exit(1);
   }   
   for(Int_t i=0;i<n_cut;i++){

       TMarCutMC newCut;
       in >> newCut;
       List.CutList.push_back(newCut); 

   }
   return in;
};

ostream& operator<<(ostream& out, TMarCutMCList& List){

   out << "Cuts lists: "<< endl;
   for (list<TMarCutMC>::iterator it= List.CutList.begin();it != List.CutList.end(); it++)
   {
      
      out << *it<< endl;
   }
   return out;
};

---