class HParticleCutRange: public TNamed

_HADES_CLASS_DESCRIPTION



 HParticleCutRange

 Helper class for cuts.
 A cut has a name (should be unique), a cut number (better unique)
 and a condition as argument of creation. The cut object owns counters
 for the number of calls and the number of failed evaluations. The
 The counters cand be used to monitor the efficiency of the cut
 If the TCutG object belonging to the cut is set  evalG(var,var,version)
 will evaluate the graph. cut. As option TF1 objects cann be used to
 cut (see the exapmple macro below)
 The cuts can be inversed by setInverse(Bool_t). In fact
 cuts will than select outside the band low-up range and
 for the graph cut isInside(var,var)==kFALSE will be used for selection.


 USAGE:
  // outside eventloop:
  HParticleCutRange cut1("BetaCut",1,0.9,1.1);
  // will create a cut object selecting beta range and

  //  inside eventloop:
  if(cut1->eval(beta,0)) { // true if condition is fullfilled, stat counted for version 0

  }
  if(cut1->eval(beta,1)) { // true if condition is fullfilled, stat counted for version 1

  }


  // after eventloop
  cut1.print(); // show cut name,number,cut statistics and range

void testCutRange()
{
    // example macro for HParticleCutRange demonstrating the
    // different features. The cut will be applied on the
    // mass of an particle in this example.

    // Cuts can be applied in 3 ways:

    // CUTMODE:
    // 1. 1-dim cuts (upper+lower limit)
    // 2. 2-dim cut using TCutG
    // 3. 2-dim cut using TF1

    // In TCutG or TF1 mode the example cuts in mass depends on the momentum.
    // The TCutG is owned by the cut and the user just has to set the
    // points. TF1 have to be contructed outside the clas due to the
    // available contructors of TF1. The user ahs to set the pointers
    // to the TF1. Remark: if the TF1 is contructed from a user function
    // the cut will not work when read back from root file.

    // TF1MODE:
    // There are 4 options :
    // 1. low+up  (1 seperate TF1 for lower and upper cut)
    // 2. low     (1 TF1 for lower cut)
    // 3. up      (1 TF1 for upper cut)
    // 3. mean+width (1 TF1 for mean of the cut, 1 TF1 for the width of the cut)


    Bool_t   invert   = kFALSE;  // invert cut : default: kFALSE
    Int_t    cutmode  = 2; //0=range,1=cutG,2=TF1
    TString  tf1mode  = "low+up"; // TF1 mode: low+up, low,up,mean+width

    Double_t slope    = .10;

    Double_t mass_pion   = HPhysicsConstants::mass(8);
    Double_t mass_proton = HPhysicsConstants::mass(14);


    HParticleCutRange cut1("massCut",1,mass_proton-50,mass_proton+50);
    cut1.setInverse(invert);


    //--------------------------------------------------
    // TF1 mode
    // Cuts in mass depend on mommentum
    // low+up , low , up mode
    TF1* flow = new TF1("flow","pol1",0,1000);
    flow->SetParameter(0,mass_pion);
    flow->SetParameter(1,-slope);

    TF1* fup = new TF1("fup","pol1"  ,0,1000);
    fup->SetParameter(0,mass_pion);
    fup->SetParameter(1,slope);
    fup->SetLineColor(4);

    // mean+width mode
    TF1* fmean = new TF1("fmean","pol1",0,1000);
    fmean->SetParameter(0,mass_pion);
    fmean->SetParameter(1,0);

    TF1* fwidth = new TF1("fwidth","pol1",0,1000);
    fwidth->SetParameter(0,0);
    fwidth->SetParameter(1,slope);
    fwidth->SetLineColor(4);

    if(tf1mode=="low+up")    cut1.setTF1(flow,fup,tf1mode);
    if(tf1mode=="low")       cut1.setTF1(flow,  0,tf1mode);
    if(tf1mode=="up")        cut1.setTF1(   0,fup,tf1mode);
    if(tf1mode=="mean+width")cut1.setTF1(flow,fup,tf1mode);
    //--------------------------------------------------


    //--------------------------------------------------
    // TCutG mode
    cut1.getTCutG()->SetPoint(0,   0,mass_proton);
    cut1.getTCutG()->SetPoint(1,1000,mass_proton - slope*1000);
    cut1.getTCutG()->SetPoint(2,1000,mass_proton + slope*1000);
    cut1.getTCutG()->SetPoint(3,   0,mass_proton);
    //--------------------------------------------------



    TH2F* hmomvsmass = new TH2F("hmassvsmom","mom vs mass",200,0,2000,200,0,1200);
    hmomvsmass->SetXTitle("mom [MeV/c]");
    hmomvsmass->SetYTitle("mass [MeV/c^{2}]");

    //--------------------------------------------------
    // create some dummy spectrum
    Double_t masses [] = {mass_pion,mass_proton};

    for(Int_t i = 0 ; i <10000; i ++){
        Double_t mom     =  gRandom->Rndm()*2000;
	Double_t massres =  mom*0.1;
        Int_t    index   = (Int_t) (gRandom->Rndm()*2);
        Double_t mass    = gRandom->Gaus(masses[index],massres);

	if(cutmode == 0 && cut1.eval(mass)) {
	    hmomvsmass->Fill(mom,mass);
	}

	if(cutmode == 1 && cut1.evalG(mom,mass)) {
	    hmomvsmass->Fill(mom,mass);
	}

	if(cutmode == 2 && cut1.evalF(mass,mom)) {
	    hmomvsmass->Fill(mom,mass);
	}
    }
    //--------------------------------------------------


    cut1.print();

    TCanvas* result = new TCanvas("result","result",1000,800);

    hmomvsmass->Draw("colz");


    if(cutmode == 1){
        cut1.getTCutG()->Draw("same");
    }

    if(cutmode == 2){
	if(cut1.getTF1Low())cut1.getTF1Low()->DrawCopy("same");
	if(cut1.getTF1Up() )cut1.getTF1Up() ->DrawCopy("same");
    }

    TFile* cutsfile = new TFile("myCuts.root","RECREATE");
    cut1.Write();
    cutsfile->Save();
    cutsfile->Close();


}