test_phi2170_ingo.C

Go to the documentation of this file.

{
 
 
    //acticate the plugin in advance:
    //makeDistributionManager()->Exec("exotics:init");
    //loop over std physics to generate the model objects
    //makeDistributionManager()->ActivateStdModels();
    //listParticle("phi2170");
    
    //PHadronDecayM2Mesh  *model = new PHadronDecayM2Mesh("m2@phi2170#to#phi#f0_980", "Mesh decay for phi2170", -2);
    //makeDistributionManager()->Add(model);
    
    PReaction r(29., "p", "p", "p p phi2170 [phi [K+ K-] f0_980]");
    //PReaction r(29., "p", "p", "p p phi2170 [phi f0_980]");
 
    TH1F *histo = new TH1F("histo", "phi2170 mass", 100, 1.6, 3.0);
    TH2F *histo2 = new TH2F("histo2", "phi2170 vs f0 mass", 100, 1.6, 2.6, 100, 0.6, 1.2);
    r.Do(histo, "_x = [phi2170]->M()");
    r.Do(histo2, "_x = [phi2170]->M(); _y=[f0_980]->M();");
 
    r.Print();
 
    r.Loop(100000);
 
        //get the decay indices:
    Int_t phi2170_to_phi_pipi_idx = makeStaticData()->GetDecayIdx("phi2170", "phi, pi+, pi-");
    Int_t phi2170_to_KKpipi_idx   = makeStaticData()->GetDecayIdx("phi2170", "K+, K-, pi+, pi-");
    Int_t phi2170_to_phif0_idx    = makeStaticData()->GetDecayIdx("phi2170", "phi, f0_980");
    
    //the name of the distributions as obtained with listParticle
    PBreitWigner *phi2170             = dynamic_cast<PBreitWigner*>(makeDistributionManager()->GetDistribution("phi2170_bw")->Clone());
    PBreitWigner *phi2170_to_phi_pipi = dynamic_cast<PBreitWigner*>(phi2170->Clone());
    PBreitWigner *phi2170_to_KKpipi   = dynamic_cast<PBreitWigner*>(phi2170->Clone());
    PBreitWigner *phi2170_to_phif0    = dynamic_cast<PBreitWigner*>(phi2170->Clone());
    phi2170->SetParameter(0, 0);
    phi2170->SetParameter(1, -1);
    phi2170->SetName("total");
    phi2170_to_phi_pipi->SetParameter(1, phi2170_to_phi_pipi_idx);
    phi2170_to_phi_pipi->SetName("phi_pipi");
    phi2170_to_KKpipi->SetParameter(1, phi2170_to_KKpipi_idx);
    phi2170_to_KKpipi->SetName("KKpipi");
    phi2170_to_phif0->SetParameter(1, phi2170_to_phif0_idx);
    phi2170_to_phif0->SetName("phif0");
 
    //dummy to create the meshes
    Double_t mass;
    phi2170->GetWidth(2.2,&mass);
 
    phi2170_to_phi_pipi->SetLineColor(1); //pipi black
    phi2170_to_phif0->SetLineColor(3); //phif0 green
 
    histo->Scale(1000. * 0.3 * 0.7 * 0.95 * 0.8 * 0.9);
    histo->Draw();
    phi2170_to_phif0->Draw("same");
 
    
}